24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

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

mengfc

金虫 (正式写手)

应助: 20 (小学生)
金币: 1301.2
散金: 1807
红花: 5
帖子: 927
在线: 524.4小时
虫号: 818484
注册: 2009-07-30
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】BSSE请教已有3人参与

第一次接触BSSE，请问如果我用MP2/6-311++G(d,p)//b3lyp/6-31+g**进行计算是否需要进行校正？计算内容是能量及一些几何参数。看论坛帖子上说优化完再用counterpoise校正，此时输入文件怎么写呢？

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2010-06-07 07:23:36

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

qzhaosdu

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 1401.6
散金: 2052
红花: 7
帖子: 1969
在线: 263.5小时
虫号: 493351
注册: 2008-01-10
性别: GG
专业: 理论和计算化学
管辖: 量子化学

mengfc(金币+3):谢谢，我现在计算的是团簇，只含一种原子，不知道如何用counterpoise分层，请指教~ 2010-06-07 21:36:10

引用回帖:

Originally posted by mengfc at 2010-06-07 07:23:36:
第一次接触BSSE，请问如果我用MP2/6-311++G(d,p)//b3lyp/6-31+g**进行计算是否需要进行校正？计算内容是能量及一些几何参数。看论坛帖子上说优化完再用counterpoise校正，此时输入文件怎么写呢？

只要是2个分子的复合物都要矫正，要是你用CCSD(T),QCISD之类的方法加大机组的话不矫正也成

赞一下

回复此楼

Anewday,anewhour,anewminute,anewpeople.

4楼2010-06-07 16:22:20

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

查看全部 5 个回答

LuPeng5366

新虫 (著名写手)

牛牛

应助: 1 (幼儿园)
金币: 0.5
散金: 1512
红花: 3
帖子: 1562
在线: 234.2小时
虫号: 738149
注册: 2009-04-02
性别: GG
专业: 原子和分子物理

★ ★
heyo_123(金币+2):鼓励交流！！ 2010-06-07 11:18:50
mengfc(金币+3):谢谢~但为什么是用用#MP2/6-31G 这个泛函和基组呢？该如何选取？ 2010-06-07 16:03:23

平衡计算的输入。这个例子使用Z-矩阵(左边)和笛卡尔坐标(右边):

# MP2/6-31G Counterpoise=2 Opt

Counterpoise with Z-matrix

0,1,0,3,1,2
O,0.0,0.0,0.0,1         结构从这里开始
O,1,ROO,2
X,1,1.,2,X3O
H,1,RO1H,3,HOX3,2,90.,0,1
H,1,RO1H,3,HOX3,2,-90.,0,1
X,2,1.,1,52.5,3,180.,0
H,2,RO2H1,6,H7OX,1,180.,0,2
H,2,RO2H2,6,H8OX,1,0.,0,2

Z-矩阵变量...

# MP2/6-31G Counterpoise=2 Opt

Counterpoise with Cartesian

0,1
1  0.00  0.00  0.92  1
9  0.17  0.00  2.73  2
1  0.77  0.00  3.43  2
9  0.00  0.00  0.00  1

注意，Z-矩阵输入数据需要在二面角数值/变量之后要有一个0（表示最后的角是二面
角），然后才是片段编号。另外，Z-矩阵的第一个原子必须用笛卡尔坐标。对于这样的任务，
显然用笛卡尔坐标定义片段编号更简单。
前面的Z-矩阵也描述指定片段使用的电荷与自旋多重度。在这个例子中，相应输入行的
格式是：

- 66 -

总电荷, 总自旋, 片段1的电荷, 片段1的自旋多重度, 片段2的电荷, 片段2的自旋多重度

使用ECP的平衡优化例子：

# hf/lanl2dz counterpoise=2 nosymm opt test

HBr + HF, optimization with counterpoise correction using ECP basis

0 1
H -0.046866  0. 0.586860 1
Br  -0.331864  0.  -0.801000 1
F 0.396755  0. 2.739275 2
H 0.584835  0. 3.641534 2

平衡计算的输出。这里是Counterpoise计算输出的例子：

Counterpoise: corrected energy = -2660.083831739527
Counterpoise: BSSE energy =    0.003902746890

这几行分别给出修正能量和基组重迭误差。

赞一下(1人)

回复此楼

自信，坚强，勇敢。

2楼2010-06-07 10:50:17

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