24小时热门版块排行榜    

查看: 268  |  回复: 1
当前主题已经存档。

armadillo

铜虫 (初入文坛)

[交流] 【求助】G03计算氢键体系BSSE校正的问题

计算一个硫脲-甲醛的双氢键体系(分别用A,B,AB表示),A,B,AB均用B3LYP/6-31G(d)进行优化并得到TCE焓校正项,然后用MP2/6-31+G(d,p)计算A,B单点,AB单点;用MP2/6-31+G(d,p)和Counterpoise参数计算校正BSSE后的体系能量Counterpoise: corrected energy = -661.438917540144,Counterpoise: BSSE energy = 0.001963618146。但是从未经BSSE校正得到AB单点能为:-660.521304361

如果根据经过BSSE校正的单点能计算氢键能量则明显偏大(达到上千kJ/mol),但根据未经过BSSE校正的单点能计算氢键能量则相对合理(十几kJ/mol)。而且通常不是认为BSSE会使体系能量过分降低么?为什么修正后体系能量不升反降?还是由于这个“Counterpoise: corrected energy”包含其它能量?

请不吝赐教,谢谢!!

[ Last edited by zzgyb on 2008-8-5 at 10:09 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tutu2000

金虫 (正式写手)

建议楼主重新检查计算Counterpoise的输入文件,从结果来看,可能是计算的输入有问题,因为:
未经BSSE校正的AB单点能=corrected energy-BSSE energy
2楼2008-08-06 11:48:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 armadillo 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见