计算化学公社

标题: 基组问题 [打印本页]

作者
Author: wangweihua 时间: 2016-2-26 15:25
标题: 基组问题
老师，我的体系是带一个正电荷的离子，用高斯程序，M06L方法，对Te原子用DZP基组，现在审稿人让Te原子用SDD基组，Te data have to berecomputed at nonrelativistic with a relativistically adjustedpseudopotential (such as SDD)，但是换成SDD基组后我的构型在502出错，构型完全变了，扭曲的很厉害，请问老师怎么解决？

作者
Author: sobereva 时间: 2016-2-26 15:49
DZP这种老掉牙非主流基组，无论什么元素我都很不建议用。

Te是需要考虑相对论效应的，直接用DZP由不考虑相对论哈密顿，结果是不准确的。

审稿人让你用SDD，是赝势基组，可以等效地把标量相对论体现出来。

如果构型变得明显不合理，多半原因是你只定义了赝势基组，但是没有定义赝势。输入文件一看就清楚了。

作者
Author: wangweihua 时间: 2016-2-26 19:50
我的输入文件是
# M06L/Genecp opt freq=noraman

1 2
O                0.01978500 1.06786500 -1.16427500
Te             -0.22002500 -1.68460700 -1.07835800
C                -1.92287300 -0.81701300 -0.20438800
C                -2.05165400 0.59014600 -0.06180100
C                -1.08755200 1.54284100 -0.46994300
C                1.29500700 1.45508000 -0.79470800
C                1.07785200 -1.41625500 0.52595800
H                -2.81172700 -2.74158800 0.07730600
C                -3.40768300 2.49054200 0.72580000
C                -4.25278700 0.18861900 0.94991900
C                -4.10534700 -1.16313800 0.77986600
C                -2.92805300 -1.66719500 0.19923300
C                -3.24815000 1.09837100 0.54295500
H                -4.32163400 2.85226100 1.19052900
H                -5.15773500 0.59023000 1.40097100
H                -4.88503500 -1.85155000 1.09287000
H                -0.46056900 -0.96042300 1.97375400
C                -1.26336300 2.89258000 -0.28920700
C                -2.43452000 3.36844800 0.32181000
C                2.27331700 1.29998200 -1.77470600
C                3.59026500 1.60502900 -1.46635000
C                3.92769800 2.04934100 -0.18828400
C                2.93862100 2.19033400 0.78146900
C                1.61491800 1.88751400 0.49001400
C                2.44934300 -1.59196000 0.30778700
C                3.34132800 -1.41918900 1.35295000
C                2.86976000 -1.08867900 2.62266700
C                1.50319600 -0.92439300 2.84437400
C                0.60618500 -1.08602300 1.80054200
H                -0.49635400 3.58062800 -0.63484500
H                -2.56527400 4.43822600 0.45753600
H                1.98185200 0.95773500 -2.76454500
H                4.35709100 1.50114500 -2.22921400
H                4.95982500 2.29253000 0.04795200
H                3.19796800 2.53401300 1.77973700
H                0.83994800 1.98758000 1.24446400
H                2.81924200 -1.84735200 -0.68325100
H                4.40658800 -1.54490800 1.17853300
H                3.57056500 -0.96463200 3.44442100
H                1.13524000 -0.67716900 3.83684300

C H O
6-311g(d)
****
Te
SDD
****

作者
Author: wangweihua 时间: 2016-2-26 19:51
老师，是不是应该在gen后面加ecp？

作者
Author: liyuanhe211 时间: 2016-2-26 23:17

wangweihua 发表于 2016-2-26 19:51
老师，是不是应该在gen后面加ecp？

仔细阅读：http://sobereva.com/60

作者
Author: sobereva 时间: 2016-2-26 23:35

wangweihua 发表于 2016-2-26 19:50
我的输入文件是
# M06L/Genecp opt freq=noraman

很明显你只定义了Te的SDD赝势基组但没给它指定SDD赝势，结果明显不合理。输入文件怎么写看5L的链接。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)