Molpro计算单点出错

那丢失的

一个分子B3LYP优化有两个构象，是都是Cs对称性，用molpro计算单点
反式输入如下正常结束
***,antiCH2CHSO
memory,7000,m
print,basis,orbitals
geometry
C                  1.32715100   -1.60645400    0.00000000
H                  1.39001200   -2.68529900    0.00000000
H                  2.26133300   -1.05910100    0.00000000
C                  0.14354200   -0.98574700    0.00000000
H                 -0.80494600   -1.51046700    0.00000000
S                  0.00000000    0.75453600    0.00000000
O                 -1.45881900    1.09193700    0.00000000
end

basis=AVTZ

{hf
wf,39,1,1
occ,16,4
save,2100.2
}

uccsd(t)-f12a

顺式按照上述输入则出错，出错信息见附件

但是按照以下输入则正常结束，但是这样的话能量顺反相差57kcal/mol，B3LYP计算只差不到1kcal/mol。反式按照以下输入也会出错。
***,synCH2CHSO
memory,7000,m
print,basis,orbitals
geometry
C                 -1.76847300    0.50693800    0.00003800
H                 -2.84342400    0.39576900    0.00018300
H                 -1.35729500    1.50718200    0.00044500
C                 -0.96494100   -0.56135900   -0.00018700
H                 -1.33467100   -1.58102200   -0.00035000
S                  0.77519200   -0.47795300    0.00011500
O                  1.19160100    0.95648000   -0.00015200
end

basis=AVTZ

{hf
wf,39,1,1
occ,20
save,2100.2
}

uccsd(t)-f12a

Freeman

输出文件第101行显示“Point group  C1”，说明molpro没有识别出Cs对称性，所以就不要写成occ,16,4了，而要写成occ,20。至于能量相差太大，b3lyp本来算得就不太准，建议加上色散矫正em=gd3bj，或换个m06-2x再优化试试看。

那丢失的

Freeman 发表于 2021-1-1 20:25
输出文件第101行显示“Point group  C1”，说明molpro没有识别出Cs对称性，所以就不要写成occ,16,4了，而要 ...

好的，谢谢！