计算化学公社

标题: 优化激发态总出现虚频，怎么用虚频对应原子的振动坐标与原坐标相加的方法消虚频 [打印本页]

作者
Author: 珊珊来迟 时间: 2016-10-29 22:21
标题: 优化激发态总出现虚频，怎么用虚频对应原子的振动坐标与原坐标相加的方法消虚频
在溶剂效应下优化激发态时分子结构的四个甲基上总是有四个小虚频（附件中有详细信息），尝试了很多方法，一直消不下去这四个小虚频，试过用更严格的收敛限opt=tight int=ultrafine，也尝试过调四个甲基的二面角，可是都不好使，还是会有四个小虚频。想试一下将虚频所对应的原子的振动坐标与原坐标相加的方法，感觉这个方法消虚频更有针对性，可是不知道怎么找虚频所对应的振动坐标。。 (, 下载次数 Times of downloads: 39) (, 下载次数 Times of downloads: 34)

频率下边对应三组XYZ的值，请问怎么根据这个判断振动模式的？虚频对应的振动坐标是哪个呢？对于我这个文件来说，我应该把虚频对应的振动坐标加到原输入文件的哪个坐标上呢？希望老师能指点一二~先行谢过啦

作者
Author: liyuanhe211 时间: 2016-10-30 02:13
本帖最后由 liyuanhe211 于 2016-10-30 02:14 编辑

GaussView观看振动的界面下面直接有按照位移移动后相加并保存新结构的功能。不过除非你的结构被对称性限制了，否则这种做法没什么用，不建议尝试。
另外你的输出文件里没有opt，也没有int=ultrafine，可以考虑直接读取现在计算的freq继续opt，看是否可以；另外注意如果opt里带了int=ultrafine，freq任务也必须带

作者
Author: bomsaude 时间: 2016-10-30 05:50
Are you sure you are calculating the frequency for the excited state when you use the keywords:
# freq cam-b3lyp/6-31g(d,p) scrf=(solvent=n-hexane) ?

作者
Author: sobereva 时间: 2016-10-30 08:30
当前体系呈镜面对称性，从虚频上看明显边缘的SiH3倾向于破坏这种对称性。
所以，你当前的情况正是我博文里说的由于初始结构对称性太高导致的虚频，这种情况用ultrafine、tight根本没用。按照虚频模式略微调节破坏掉镜面对称再优化。
哪种解决虚频的方式好使，一定要看具体情况。

作者
Author: 冰释之川 时间: 2016-10-30 09:20

sobereva 发表于 2016-10-30 08:30
当前体系呈镜面对称性，从虚频上看明显边缘的SiH3倾向于破坏这种对称性。
所以，你当前的情况正是我博文里 ...

搭车咨询一下：如果优化一个含有57个原子(含有S、F、B、O、N、C)的体系的激发态, opt td freq 联用，在4核3.5G内存配额下，一般能跑的动吗？

作者
Author: sobereva 时间: 2016-10-30 10:06

冰释之川发表于 2016-10-30 09:20
搭车咨询一下：如果优化一个含有57个原子(含有S、F、B、O、N、C)的体系的激发态, opt td freq 联用，在4 ...

TD下opt凑合，freq就没戏了。我是指6-31G*基组下。

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