优化[Tl(H2O)6]+离子时出错

shiguang

我在优化[Tl(H2O)6]+离子时，不是出现I9999报错就是出现小虚频。
这是我的gjf文件，我用的B3LYP，加D3BJ校正就I9999报错加结构不合理。
%chk=[Tl(H2O)6]+.chk
%mem=5GB
%nprocshared=64
#p freq ub3lyp/gen pseudo=cards em=gd3bj opt=(calcall,MaxStep=5,NoTrustUpdate,Cartesian)


[Tl(H2O)6]+ opt&freq


1 1
O                 -2.52486500    0.30175300    1.55708100
H                 -2.79886300   -0.56468100    1.20155400
O                  2.53564600   -1.13715000    1.09646600
H                  2.83611900   -1.27229900    0.17764800
O                 -2.61720900   -1.47235500   -0.43154800
H                 -2.89261100   -0.72390100   -0.99366800
O                 -2.65864500    1.14126700   -0.98037000
H                 -3.05791700    1.84391300   -1.50265400
O                  2.68680500   -0.36046000   -1.44828000
H                  2.91793500    0.51973600   -1.09766900
O                  2.58328500    1.46872800    0.50487100
H                  2.81137000    0.73261600    1.10266100
H                 -3.00059200   -2.28845800   -0.76756000
H                 -2.85204800    0.40173700    2.45650900
H                  2.87402600   -1.84212600    1.65697300
H                 -2.87457700    1.26958100   -0.03837200
H                  3.11105200   -0.49041400   -2.30213900
H                  2.94577300    2.29090700    0.84911500
Tl                -0.00073800    0.00727300   -0.03861900


Tl 0
SDD
****
O H 0
6-311g(d)
****


Tl 0
SDD



请问Tl这种体系应该用什么泛函计算？如何解决有小虚频的问题？（沿着虚频的方向不是对称结构）

北大-陶豫

你用 ub3lyp，但是初猜也要加 guess=mix 破坏 alpha beta 自旋的对称性，否则 unrestricted 也白搭。
l9999 不是 I9999

sobereva

仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

chemhou

ecp没用对

paramecium86

看起来优化曲线能量一直在降低。也没振荡。 应该是的步长太小了。还没有优化到位就到最大步数了。可以用出错的位置的结构作为初始结构再继续优化，用正常的步长。 要是怕有虚频就用 opt(recalc=5) 之类的。没必要一上来就把步长调的这么小。