计算化学公社

标题: 优化[Tl(H2O)6]+离子时出错 [打印本页]

作者
Author: shiguang 时间: 2022-1-20 17:48
标题: 优化[Tl(H2O)6]+离子时出错
我在优化[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这种体系应该用什么泛函计算？如何解决有小虚频的问题？（沿着虚频的方向不是对称结构）

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

作者
Author: chemhou 时间: 2022-1-20 19:32
ecp没用对

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

作者
Author: 北大-陶豫 时间: 2022-1-21 22:28
你用 ub3lyp，但是初猜也要加 guess=mix 破坏 alpha beta 自旋的对称性，否则 unrestricted 也白搭。
l9999 不是 I9999

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