求助：opt+freq计算脲类分子出现虚频的问题

AnnieXu · 发表于 Post on 2026-4-3 17:15:32

各位老师好，我用PBE1PBE/6-311G* opt freq em=GD3 scrf=solvent=water进行了20种脲类分子的计算，高斯版本为G09 D0.1，但有15个都出现虚频，且大于-200cm-1。通过阅读Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元，PBE1PBE/6-311G* opt=(calcfc,tight) int=ultrafine freq em=GD3 scrf=solvent=water这样能成功消除虚频。有以下疑问：
1、帖子里说取最后的结构重新优化，并在opt里写上calcfc，但我没有用最后的结构重新优化，而是直接采用opt=(calcfc,tight) int=ultrafine成功消除虚频，这样做可行吗？还是说先用opt=tight int=ultrafine再用opt=calcfc？
2、在卢老师回复http://bbs.keinsci.com/forum.php ... 7&fromuid=83569帖子时说“不管你最后用什么算能量，虚频还是得消（虽然对计算结合能计算其实没有影响。如果你不消的话别让审稿人知道）” 那我到底需要消除虚频吗？我用sobtop和multiwfn生成itp文件及RESP电荷，虚频对这个有影响吗？
3、那我原本没有虚频的分子是否也要用当前的方案？

wal · 发表于 Post on 2026-4-3 17:55:14

1. 顺序随意消掉就行
2. 感觉sob老师这么讲是因为原帖那种几个波数的小虚频对能量影响很轻微，实在消不掉睁只眼闭只眼也能过去。对电荷有没有影响不清楚留给大佬回答吧
3. int=ultrafine改了格点的话所有分子理论上是要保持一致的

AnnieXu · 发表于 Post on 2026-4-3 19:58:05

本帖最后由 AnnieXu 于 2026-4-3 19:59 编辑

wal 发表于 2026-4-3 17:55
1. 顺序随意消掉就行
2. 感觉sob老师这么讲是因为原帖那种几个波数的小虚频对能量影响很轻微，实在消不掉 ...

感谢您的回复，截取最后一个结构，做opt=calcfc, int(ultrafine)
%mem=40GB
%nprocshared=40
%oldchk=a-urea.chk
%chk=new-urea.chk
# PBE1PBE/6-311G* opt=calcfc int=ultrafine freq em=GD3 scrf=solvent=water geom=allcheck guess=read

这样设置是正确的吗，用%oldchk以及geom=allcheck guess=read

wal · 发表于 Post on 2026-4-3 20:13:13

AnnieXu 发表于 2026-4-3 19:58
感谢您的回复，截取最后一个结构，做opt=calcfc, int(ultrafine)
%mem=40GB
%nprocshared=40 ...

应该没啥问题

AnnieXu · 发表于 Post on 2026-4-3 20:51:18

wal 发表于 2026-4-3 20:13
应该没啥问题

谢谢老师

sobereva · 发表于 Post on 2026-4-5 03:33:48

200i cm-1这样的虚频已属于显著虚频，能消就消，也不可能超级难消掉（和几个波数小虚频的情况完全不同）。虽然算RESP电荷对结构优化精度并不很敏感，但但凡消虚频不费劲的时候我都建议消了，除非你凭化学直觉根据振动模式判定这对于结构合理性不造成显著影响

AnnieXu · 发表于 Post on 2026-4-5 15:39:05

本帖最后由 AnnieXu 于 2026-4-5 22:57 编辑

sobereva 发表于 2026-4-5 03:33
200i cm-1这样的虚频已属于显著虚频，能消就消，也不可能超级难消掉（和几个波数小虚频的情况完全不同）。 ...

感谢老师
1、我的一部分分子加上opt=tight int=ultrafine能消除，有一些得在这个基础上继续用opt=(calcfc,tight) int=ultrafine，但calcfc 只是优化算法，不改变最终精度，所以我这一批分子应该都是在同一精度下的吧
2、此时算单点能是不是要加上int=ultrafine才行呢？
3、老师您之前在Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法http://bbs.keinsci.com/forum.php ... d=633&fromuid=83569提到：默认的int=fine和int=ultrafine对结构影响甚微，所以只需要在相同格点设定下算能量就够了。我的硫酸根离子用opt=tight int=ultrafine就跑不出来，单纯写opt freq或opt freq int=ultrafine均能成功结束，如果硫酸根不用opt=tight会不会与别的分子使用精度不一样呢，这样可行吗？

sobereva · 发表于 Post on 7 day ago

AnnieXu 发表于 2026-4-5 15:39
感谢老师
1、我的一部分分子加上opt=tight int=ultrafine能消除，有一些得在这个基础上继续用opt=(calcf ...

1 是
2 G16本来就默认ultrafine。单点能用的积分格点必须统一
3 是否用tight不明显影响可比性，因为默认的几何优化收敛限就已经比较严了。PS：即便用tight，若结合恰当的关键词，也并非无法达到收敛限

AnnieXu · 发表于 Post on 7 day ago

sobereva 发表于 2026-4-6 02:55
1 是
2 G16本来就默认ultrafine。单点能用的积分格点必须统一
3 是否用tight不明显影响可比性，因为默 ...

谢谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 求助：opt+freq计算脲类分子出现虚频的问题

浏览过的版块