18碳环加电场后的结构优化为什么会出现虚频

chuan437

Greetings！大家好。
我对社长那篇对18碳环加电场的研究很感兴趣。自己也尝试着算了一下。我用的泛函和基组跟社长的文献相同，ωB97XD/Def2TZVP
这里我遇到了一个问题：如果施加x轴方向（和社长论文中相同）的0.025 a.u.的电场，进行结构优化时能够正常收敛。但是如果同时进行振动分析的话，会发现其中有一个（有且只有一个）虚频。我想问一下这其中的原因是什么。这个虚频的出现是否和施加了电场相关（即施加电场得到的优化结构都存在虚频），还是虚频说明结构优化没有收敛到极小点构型。谢谢！
附图是我的振动分析结果。

sobereva

优化时加了电场，振动分析也记得必须加电场，否则不是同一个势能面

chuan437

sobereva 发表于 2025-6-8 00:58
优化时加了电场，振动分析也记得必须加电场，否则不是同一个势能面

社长你好，我的优化和振动分析是同时进行的，输入文件如下：
%chk=./c18_x250.chk
%nproc=36
%mem=60GB
# wB97XD/Def2TZVP field=x+250 opt nosymm freq

Title Card Required

0 1
C                 -0.61025100    3.64209300    0.00000000
C                 -1.87361300    3.18226700    0.00000000
C                  0.61025100    3.64209300    0.00000000
C                  1.87361300    3.18226700    0.00000000
C                  2.80857200    2.39774300    0.00000000
C                  3.48079300    1.23342300    0.00000000
C                 -2.84902000   -2.34954000    0.00000000
C                  2.84902000   -2.34954000    0.00000000
C                 -1.81911800   -3.21373000    0.00000000
C                  1.81911800   -3.21373000    0.00000000
C                 -0.67222100   -3.63116600    0.00000000
C                  0.67222100   -3.63116600    0.00000000
C                 -2.80857200    2.39774300    0.00000000
C                 -3.45927100   -1.29255400    0.00000000
C                  3.69273100    0.03146300    0.00000000
C                  3.45927100   -1.29255400    0.00000000
C                 -3.69273100    0.03146300    0.00000000
C                 -3.48079300    1.23342300    0.00000000

sobereva

chuan437 发表于 2025-6-8 03:23
社长你好，我的优化和振动分析是同时进行的，输入文件如下：
%chk=./c18_x250.chk
%nproc=36

按下文说的消虚频
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）