计算化学公社

标题: 激发态优化始终有虚频 [打印本页]

作者
Author: ter20 时间: 2016-9-22 12:59
标题: 激发态优化始终有虚频
之前在论坛上求教过激发态优化的问题，现在有几个结构的S2或者S3优化完总是有一个虚频，试了各种办法，包括沿虚频振动方向不同程度地改结构，读取前一步优化结果ReadFC继续优化，改小Maxstep, 加notrust, 加nosymm，使用GDIIS等等，都消不掉这个虚频，有时会得到能量更高的含1个或2个虚频的结构。我也试了沿着虚频振动的方向手动取点扫描单点能量，发现优化得到的机构的确是能量极小点，实在无法理解为什么总是有一个虚频。这里放上其中一个例子的输入文件，结果log文件见附件。求大家的意见，多谢!

# CAM-B3LYP/aug-cc-pVDZ opt=vtight Int=UltraFine TD=(Root=3,Nstates=6) freq=noraman

Title Card Required

1 1
N                0.00000000 0.00000000 -1.27637700
C                0.00000000 1.21334500 0.75536600
C                0.00000000 1.20821800 -0.61236700
C                0.00000000 -1.20821800 -0.61236700
C                0.00000000 -1.21334500 0.75536600
N                0.00000000 0.00000000 1.30870700
H                0.00000000 2.11435200 1.36644800
H                0.00000000 2.12681700 -1.19296200
H                0.00000000 -2.12681700 -1.19296200
H                0.00000000 -2.11435200 1.36644800
H                0.00000000 0.00000000 -2.28926800

作者
Author: sobereva 时间: 2016-9-22 21:05
把弥散去了试试

作者
Author: beefly 时间: 2016-9-22 21:55
虚频模式是平面外方向的，应该打破平面对称性重新优化

作者
Author: ter20 时间: 2016-9-23 07:11

beefly 发表于 2016-9-22 21:55
虚频模式是平面外方向的，应该打破平面对称性重新优化

之前尝试过使用 nosymm，没有成功，请问是不是在建模时要加上dummy atom那样处理，同时再加上nosymm关键字？

作者
Author: ter20 时间: 2016-9-23 09:43
本帖最后由 ter20 于 2016-9-23 09:58 编辑

beefly 发表于 2016-9-22 21:55
虚频模式是平面外方向的，应该打破平面对称性重新优化

我尝试加入2个X atom，建了一个模型略扭曲的非平面结构，但是有C2对称性，没有加nosymm关键字，优化结果收敛到之前得到的一个有1个虚频的结构（因为结构和能量都一致，所以这个结构没有做频率分析），如附件所示，不知道我这种做法是否是您建议的意思？谢谢

作者
Author: sobereva 时间: 2016-9-23 10:16

ter20 发表于 2016-9-23 09:43
我尝试加入2个X atom，建了一个模型略扭曲的非平面结构，但是有C2对称性，没有加nosymm关键字，优化结果 ...

打破对称性用不着加X
直接在gview里扭曲个二面角让一开始对称性降低即可
对于小分子，基态有高对称性的情况，激发态时往往对称性会降低

(, 下载次数 Times of downloads: 64)

作者
Author: ter20 时间: 2016-9-23 11:07

sobereva 发表于 2016-9-23 10:16
打破对称性用不着加X
直接在gview里扭曲个二面角让一开始对称性降低即可
对于小分子，基态有高对称性的 ...

好的，我继续尝试，多谢Sob!

作者
Author: bomsaude 时间: 2016-9-23 18:01

ter20 发表于 2016-9-23 03:07
好的，我继续尝试，多谢Sob!

这个小分子的阳离子激发态也太难优化了，我把你的初始分子结构和调成非平面的结构（两个N原子翘起来），用你的方法，都还是有虚频。另外，我用非平面的结构，换了不带弥散的小基组，虚频更多，而且最后结构又变成了平面结构。总之优化不出来无虚频的。

作者
Author: ter20 时间: 2016-9-24 10:36

bomsaude 发表于 2016-9-23 18:01
这个小分子的阳离子激发态也太难优化了，我把你的初始分子结构和调成非平面的结构（两个N原子翘起来）， ...

是的，我试了各种办法，总是有虚频，但是S1和S2就没有，真的很抓狂，可是又不想再改用Turbomole或者Molcas，那样又要买软件还得学习。
感谢您还亲自做了计算，真的很感激

作者
Author: ter20 时间: 2016-9-25 11:06

sobereva 发表于 2016-9-23 10:16
打破对称性用不着加X
直接在gview里扭曲个二面角让一开始对称性降低即可
对于小分子，基态有高对称性的 ...

试了改变初始构型，手动将其改的完全不对称，结果还是不行，不知道是不是高斯本身的问题？

作者
Author: sobereva 时间: 2016-9-25 14:54

ter20 发表于 2016-9-25 11:06
试了改变初始构型，手动将其改的完全不对称，结果还是不行，不知道是不是高斯本身的问题？

没具体信息不好说
你可以用CIS也试试，还有虚频再结合calcall试试。如果之后没虚频了，可以拿它当初猜做TDDFT优化。
如果calcall有效，也可以等Gaussian下个大版本出来，会支持TDDFT解析Hessian，就也能用calcall了

作者
Author: ter20 时间: 2016-9-25 15:47

sobereva 发表于 2016-9-25 14:54
没具体信息不好说
你可以用CIS也试试，还有虚频再结合calcall试试。如果之后没虚频了，可以拿它当初猜做 ...

好的，我再继续试试，多谢Sob!

作者
Author: ter20 时间: 2016-9-25 15:55

sobereva 发表于 2016-9-25 14:54
没具体信息不好说
你可以用CIS也试试，还有虚频再结合calcall试试。如果之后没虚频了，可以拿它当初猜做 ...

还有一个小问题，就是我使用不对称的结构做初始，也会优化得到对称性的结果。为了避免这种平面虚频激发态的情形，需要加上Nosymm关键字吗？

作者
Author: sobereva 时间: 2016-9-25 21:02

ter20 发表于 2016-9-25 15:55
还有一个小问题，就是我使用不对称的结构做初始，也会优化得到对称性的结果。为了避免这种平面虚频激发态 ...

nosymm此时没用。初始结构没对称性的话，和用nosymm强制不用对称性是一样的
如果优化到了对称，说明对称结构本身就是极小点应有的特征。如果对称结构有破坏对称的虚频模式，说明优化过程不合理，需要用calcall，或减小步长之类收敛到真正的非对称的极小点结构

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