计算化学公社

标题: 18碳环加电场后的结构优化为什么会出现虚频 [打印本页]

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

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

作者
Author: chuan437 时间: 2025-6-8 03:23

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

作者
Author: sobereva 时间: 2025-6-9 16:55

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）

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