| 查看: 1399 | 回复: 10 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】如何用counterpoise计算BSSE?
|
|||
| 【求助】如何用counterpoise计算BSSE?具体怎么操作这个关键词，再问这个关键词能用于计算金属的BSSE吗？十分感谢！ |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有223人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
七嗪类物质合成求助
已经有0人回复
-
分子动力学模拟合作需求
已经有3人回复
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
2楼2009-02-22 13:11:09
loovfnd
至尊木虫 (著名写手)
- 应助: 35 (小学生)
- 贵宾: 1.111
- 金币: 14361.4
- 散金: 236
- 红花: 4
- 帖子: 2790
- 在线: 187小时
- 虫号: 219674
- 注册: 2006-03-14
- 专业: 物理有机化学
3楼2009-02-22 16:40:38
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
★ ★ ★ ★ ★
xuefei06(金币+5,VIP+0):thanks！ 2-23 00:22
xuefei06(金币+5,VIP+0):thanks！ 2-23 00:22
|
举个例子，比如计算AB体系的BSSE矫正： 需要得出的五个能量分别是： 第一个是整体的能量E(AB) 第二个是单独计算A的能量E(A) 第三个是单独计算B的能量E(B) 第四个是计算在AB基组下，A的能量E(A/AB) 第五个是计算在AB基组下，B的能量E(B/AB) 那么，BSSE energy =E(A)+E(B)-E(A/AB)-E(B/AB)。计算E(A/AB)的基函数是计算AB的基函数，而计算E(A)的基函数只是计算AB基函数的一部分，基函数越多，能量计算的越准确，能量计算越准确，就越低，所以E(A/AB) 输出文件里面5个scf第一个是整体的能量(记为E(AB))，后面4个是单体能量，分别是2个部分在整体基组和各自基组下的能量。 我当时是用-Bq原子算的BSSE矫正， 估计你是用CP方法吧。但原理都是一样的。 http://muchong.com/bbs/viewthread.php?tid=1099266 |
4楼2009-02-22 18:03:05
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
★ ★ ★
xuefei06(金币+3,VIP+0):thanks！ 2-23 00:23
xuefei06(金币+3,VIP+0):thanks！ 2-23 00:23
|
用counterpoise关键字： #uB3LYP/6-31G Counterpoise=2 SCF=(tight) C2 + NO --> C2NO 0 2 0 1 0 2 C 0.000000 1.600184 0.000000 1 C 0.000000 0.960265 0.000000 1 N 0.208401 3.031919 0.000000 2 O 1.164068 3.666273 0.000000 2 0 2 0 1 0 2分别代表总体系电荷和多重度，片段一的电荷和多重度，片段二的电荷和多重度。 分子坐标后的1和2是定义的片段 用Bq原子： # b3lyp/3-21G** nosymm scf=tight test HBr + H2O manual counterpoise calculation, H2O removed 0 1 H 0.685176 -0.004924 -0.026973 Br -0.771917 0.000050 0.001967 O-Bq 2.536864 -0.000136 -0.051401 H-Bq 3.015128 0.789231 0.184042 H-Bq 3.021888 -0.784986 0.185282 这就相当于用HBr + H2O的基组计算HBr |
5楼2009-02-22 18:11:38
loovfnd
至尊木虫 (著名写手)
- 应助: 35 (小学生)
- 贵宾: 1.111
- 金币: 14361.4
- 散金: 236
- 红花: 4
- 帖子: 2790
- 在线: 187小时
- 虫号: 219674
- 注册: 2006-03-14
- 专业: 物理有机化学
6楼2009-02-25 09:20:11
7楼2009-02-27 08:20:40
zhangji3013
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1313
- 散金: 63
- 帖子: 232
- 在线: 107.5小时
- 虫号: 532956
- 注册: 2008-03-25
- 专业: 物理有机化学
8楼2009-02-27 08:59:58
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
9楼2009-02-27 09:15:25
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
10楼2009-02-27 09:23:37
