基组问题

wangweihua

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

sobereva

DZP这种老掉牙非主流基组，无论什么元素我都很不建议用。

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

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

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

wangweihua

我的输入文件是
# 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
****

wangweihua

老师，是不是应该在gen后面加ecp？

liyuanhe211

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

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

sobereva

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

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