G16优化Li2S小分子构型错误

nemolee

用的初始构型比较接近实际，但是高斯一直努力优化为直线型构型，然后因为坐标问题报错结束。

加上guess=mix和nosymm后虽然可以正常结束，但是UM062X/6-311g(d,p)得到的键角为178.8度，使用UM062X/def2QZVP得到的键角为136度，而实际116度左右。

请各位大神帮忙看看哪里出了问题，谢谢！

#P UM062X/6-311g(d,p) opt(calcall,maxcyc=100) guess=mix nosymm

Title Card Required

0 1

Li                -1.23420250    2.67967550    0.00000000

Li                 1.43341121    0.79439252    0.00000000

S                 -0.56658879    0.79439252    0.00000000

--------------------------------------------------------------------------------------------------------

2022-08-16 更新

用CCSD(T)/cc-PVTZ优化了之后的键角是132度，暂且相信这个值吧，前人一堆文献里不到120度的不知道是咋搞出来的。。。

chands

你的初始构型是晶体还是气相分子？不明白你为甚用guess=mix。

nemolee

chands 发表于 2022-8-15 21:48
你的初始构型是晶体还是气相分子？不明白你为甚用guess=mix。

气相分子。

我知道不是自由基，各种尝试看看什么条件下可以收敛以及对键角的影响（无意义）。

mizu-bai

测了一下 wB97M-V/def2-QZVPPD 优化，两个 Li-S 均为 2.0794 Angstrom，Li-S-Li 键角 136.84 degree

mizu-bai

mizu-bai 发表于 2022-8-17 06:40
测了一下 wB97M-V/def2-QZVPPD 优化，两个 Li-S 均为 2.0794 Angstrom，Li-S-Li 键角 136.84 degree

CCSD(T)-F12a/aug-cc-pVTZ，Li-S = 2.0965 Angstrom，Li-S-Li = 137.33 degree。

nemolee

mizu-bai 发表于 2022-8-17 06:54
CCSD(T)-F12a/aug-cc-pVTZ，Li-S = 2.0965 Angstrom，Li-S-Li = 137.33 degree。

谢谢
一开始看见多篇文献写120度左右，以为自己搞错了，找了半天原因才开始怀疑，没想到真是这样。
CCSD(T)/cc-PVQZ下的结果是138.28，看来差不多就是这样了