计算化学公社

标题: 请教一下关于优化激发态一直不收敛的问题，谢谢！ [打印本页]

作者
Author: kongxp 时间: 2018-5-10 11:22
标题: 请教一下关于优化激发态一直不收敛的问题，谢谢！
我在气相下优化了一个分子T1态，加了色散矫正，#p opt BMK/6-31g(d) TD(triplet,root=1) empiricaldispersion=gd3bj，算了14天，一直不收敛。然后减小了最大允许步长，#p opt(maxstep=10,notrust) BMK/6-31g(d) TD(triplet,root=1) empiricaldispersion=gd3bj，之后又算了又一个周，还是没有收敛。所以想请教一下，还有什么解决办法，可以令其收敛，谢谢老师！
以下是优化后的基态结构：
0 1
N                -1.49664300 1.07367900 0.71767000
C                -1.37863500 0.18195500 1.78417200
C                -0.36218600 0.09613100 2.73775600
H                0.46934800 0.79239900 2.72024000
C                -0.46624200 -0.91345600 3.69870100
H                0.31475800 -1.01192900 4.45019600
C                -1.55828100 -1.80486400 3.71085000
H                -1.60475800 -2.58585300 4.46677600
C                -2.59552700 -1.67506500 2.78347300
H                -3.45651200 -2.34052000 2.81757600
C                -2.51903200 -0.65977100 1.82069600
C                -3.41036800 -0.18056800 0.77483000
C                -2.74769400 0.90939300 0.14424400
C                -3.34882600 1.63081600 -0.89695800
H                -2.83844300 2.47092000 -1.36259900
C                -4.61253100 1.22128800 -1.32077800
H                -5.09910700 1.76373400 -2.12915600
C                -5.27918400 0.13296000 -0.71533600
H                -6.26889000 -0.15291000 -1.06490400
C                -4.68117600 -0.56946200 0.33338300
H                -5.19501100 -1.40312600 0.80901400
C                -1.00788700 3.41676300 0.60826100
H                -2.00612000 3.58063200 1.00982600
C                -0.16337700 4.48310500 0.30803000
H                -0.49560300 5.50777100 0.46097400
C                1.11507400 4.21566900 -0.19387200
H                1.76965500 5.04630700 -0.45544000
C                1.57270900 2.90632200 -0.39786800
C                2.95273300 2.73737100 -1.00922900
H                2.95006800 1.97875800 -1.80007700
H                3.29260100 3.68877800 -1.43457600
H                3.69264600 2.41345100 -0.26532500
C                0.73363400 1.78600600 -0.10528000
C                -0.56237900 2.10473900 0.40417600
C                2.60783100 -0.25869600 -0.21939300
C                3.47125300 0.05119000 0.86677800
C                4.79537900 -0.41149600 0.88905500
H                5.42904500 -0.15182500 1.73819300
C                5.32496800 -1.19872600 -0.13811400
C                4.48061500 -1.51490900 -1.20984100
H                4.87064300 -2.11146300 -2.03560300
C                3.15485000 -1.06525000 -1.26358700
C                3.00050300 0.83923000 2.07839800
H                2.46104300 0.17507500 2.76818400
H                3.85566600 1.25218800 2.62704600
H                2.32812600 1.66141300 1.81666300
C                6.75619400 -1.69802100 -0.09636900
H                6.79426600 -2.76994600 0.14647600
H                7.25187600 -1.56481000 -1.06722900
H                7.34057400 -1.16337700 0.66249700
C                2.34493800 -1.42831300 -2.49986700
H                1.88479300 -0.54132100 -2.95329800
H                2.98959400 -1.89769900 -3.25366900
H                1.53354700 -2.13757300 -2.27665900
C                0.04083800 -0.65511000 -0.88404100
C                -0.85751400 -0.26760100 -2.01338400
C                -1.92631400 -1.04957400 -2.35119400
H                -2.58776700 -0.73945200 -3.16127200
C                -2.23194500 -2.27159800 -1.65758600
C                -1.32573000 -2.71931100 -0.65052300
H                -1.54186000 -3.65971200 -0.14262700
C                -0.20289100 -2.01071000 -0.29880700
C                -0.58979500 0.99688700 -2.78933000
H                0.45837900 1.03465400 -3.11865100
H                -1.24027500 1.05648000 -3.66913600
H                -0.74110600 1.89328200 -2.17575900
C                -3.48643400 -3.02616900 -1.98742900
H                -3.57347200 -3.19638200 -3.07057100
H                -3.52923800 -3.99390100 -1.47356800
H                -4.36366000 -2.43013300 -1.68393400
C                0.73823700 -2.55703700 0.74006900
H                1.00927200 -1.77370700 1.45911300
H                0.28022200 -3.39681500 1.27414800
H                1.68085700 -2.89119400 0.28277400
B                1.14129900 0.28746500 -0.36877300

作者
Author: sobereva 时间: 2018-5-10 11:34
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

T1尽量用UKS做，别用TDDFT优化，这点这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314

作者
Author: abin 时间: 2018-5-10 14:05
顺路贴个奇怪的事情。

一个水分子，用G16， TD Opt Freq B3LYP/6-31G(d) 你会发现，这个玩意居然没法收敛。

我这里是用G16A03测试的。

作者
Author: kongxp 时间: 2018-5-10 14:56

sobereva 发表于 2018-5-10 11:34
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

谢谢sob老师指点，博文很详细

作者
Author: sobereva 时间: 2018-5-10 17:41

abin 发表于 2018-5-10 14:05
顺路贴个奇怪的事情。

一个水分子，用G16， TD Opt Freq B3LYP/6-31G(d) 你会发现，这个玩意居然没 ...

越是小分子，处于激发态的时候特征相比于基态差异越大，TDDFT往往表现越不好
这种情况用多参考方法居多

作者
Author: abin 时间: 2018-5-11 09:43

sobereva 发表于 2018-5-10 17:41
越是小分子，处于激发态的时候特征相比于基态差异越大，TDDFT往往表现越不好
这种情况用多参考方法居多

大神说的对啊。
我当时忽略了这个问题。

就想着找个简单的做个算例看看。结果把自己埋坑里了。

换了一个大一点的分子，三五分钟搞定。

作者
Author: captain 时间: 2018-5-11 21:36

sobereva 发表于 2018-5-10 17:41
越是小分子，处于激发态的时候特征相比于基态差异越大，TDDFT往往表现越不好
这种情况用多参考方法居多

请问大神多参考方法最常用的是什么方法啊？
另外 UKS 是什么意思啊？帖子里没有找到相关的信息

作者
Author: sobereva 时间: 2018-5-12 04:08

captain 发表于 2018-5-11 21:36
请问大神多参考方法最常用的是什么方法啊？
另外 UKS 是什么意思啊？帖子里没有找到相关的信息

比如CASPT2，最常见的多参考二阶微扰理论，基于CASSCF多组态波函数做二阶微扰
UKS即非限制性开壳层方式做KS-DFT

作者
Author: captain 时间: 2018-5-12 16:20

sobereva 发表于 2018-5-12 04:08
比如CASPT2，最常见的多参考二阶微扰理论，基于CASSCF多组态波函数做二阶微扰
UKS即非限制性开壳层方式 ...

明白了谢大神！

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