24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yitahutu

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 483.2
帖子: 143
在线: 2.8小时
虫号: 434785
注册: 2007-08-19
性别: MM
专业: 表面、形态与形貌分析

[交流] 【求助】关于BSSE校正问题已有2人参与

根据以前论坛里的方法，用笛卡尔坐标分成了两个部分，给大家看看输入文件

0,1
O,-3.3139046751,3.274642247,6.6961469946,1
Cl,-1.9921210153,2.8284254822,6.2108903119,1
N,-2.0560431796,1.8169501812,3.1196713198,1
O,-0.9469754243,2.9246414934,7.243712536,1
O,-2.0968552773,1.3742998002,5.7368930013,1
O,-1.5889649023,3.6643874296,4.9988283616,1
H,-2.9276057887,1.8117616447,2.5893536071,1
H,-2.1621562143,1.3158494939,4.0472681451,1
H,-1.2973775872,1.4226590574,2.5460911575,1
H,-1.8309890433,2.7726167925,3.4901079337,1
C,0.1791707932,-1.5342724205,-3.789827674,2
N,1.5954603038,-1.3089088062,-3.9970643158,2
N,-0.1569909123,-2.0755839399,-2.4724737832,2
H,-0.364806944,-0.5905686841,-3.8712871814,2
H,-0.1419847962,-2.2058946396,-4.5855031932,2
C,2.4856006783,-0.7290727289,-2.9709574988,2
N,2.2523893761,-2.094542983,-4.9475905792,2
N,0.1623602984,-3.4152510878,-2.2078664299,2
C,-0.9221957836,-1.3619080904,-1.4898910626,2
O,1.5728782221,-2.6589789345,-5.7922516241,2
O,3.4780468363,-2.0973823387,-4.8608512932,2
O,-0.1199242469,-3.8343142737,-1.0918382862,2
O,0.7206030991,-4.018004671,-3.113832146,2
N,1.7326019313,-0.0050998003,-1.974846693,2
N,-0.0385486629,-0.7335952371,-0.4750940218,2
H,3.1895025574,-0.0572073498,-3.4587005746,2
H,3.0483188023,-1.5197806496,-2.4671457528,2
H,-1.6013802153,-2.04093408,-0.9770263095,2
H,-1.5056755919,-0.5872981496,-1.9936729838,2

命令用的，# MP2/6-31G Counterpoise=2 Opt

但是不能进行计算
不知道是不是哪里有问题，请教高手帮忙看看。第一次做BSSE校正，很多不明白。

[ Last edited by yjcmwgk on 2009-12-14 at 10:11 ]

回复此楼

» 猜你喜欢

华南师范大学先进光电子研究院电子纸团队诚聘2026级博士研究生已经有0人回复
QE利用声子计算的一系列lambda和Tc如何选已经有8人回复
物理化学论文润色/翻译怎么收费? 已经有169人回复
理论计算合作已经有2人回复
小木虫的论文辅导靠谱吗？有没有用过的同学？已经有1人回复
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有已经有0人回复
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后已经有81人回复
书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复
东方理-中科大联合博士生招聘已经有0人回复
半透性容器进行高湿或者低湿研究已经有0人回复

1楼2009-08-06 09:54:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yongleli

木虫 (正式写手)

应助: 93 (初中生)
金币: 3769.1
散金: 369
红花: 24
帖子: 764
在线: 101.3小时
虫号: 303595
注册: 2006-12-03
专业: 原子和分子物理

★ ★ ★ ★ ★ ★
mingdong(金币+5,VIP+0):感谢交流！ 9-1 22:28
yitahutu(金币+1,VIP+0):谢谢 9-2 09:24

转述别人的观点：

用counterpoise关键词（默认取counterpoise=2）自动计算的好处在于方便，不用自己写虚原子，而且还能比手动计算节约磁盘空间。但不利之处在于它对复合体系的能量又重新进行了单点计算（一般我们在优化构型时已经得到了此能量），耽误时间。所以，如果你很在意计算时间，手动计算BSSE倒是不错的选择。在手动计算counterpoise的BSSE时，最好用“元素符号-Bq”的方式（如C-Bq, H-Bq, ...）来定义虚原子，这样做的好处是程序能自动指派基组，而不必费力去使用Gen关键词。

      在使用G03计算counterpoise的BSSE时，最常见的程序出错是“Inaccurate quadrature in CalDSu”。其原因往往是在link401中无法确定猜测波函数。在G03中默认的猜测波函数是Harris范函，与G98中的不一样。

      另外，如果使用Bq虚原子或Gen基组手动计算BSSE时，也可能出现上述出错。其原因是因为程序的积分格点没有落在虚原子上，因此虚原子上的基函数无法正确进行积分。

   其解决办法是：

   1.  添加Guess=Indo选项。

2. 去掉对称性限制。如果你的复合体系有对称性或者近似对称且你使用了SYMM=loose选项，那么，有时也会导致上述出错。此时，如果去掉对称性（使用NOSYMM关键词），往往能解决问题。