ORCA加电场对分子优化时出现诡异虚频

ljc050512

各位老师，我在对溴苯进行加电场（x方向，大小0.01 a.u.）优化时候发现，ORCA 5.0.3会给出诡异的虚频（下图）；我用高斯16算，就不会这样。请各位老师给看看，我的输入文件哪里出了问题。
所有计算文件都在附件。

ORCA的输入文件：
! B3LYP D3 def2-SVP def2/J RIJCOSX opt freq tightSCF noautostart miniprint nopop
%maxcore  10000
%pal nprocs   20 end
%scf
  efield 0.01,0,0
end
* xyz   0   1
C     -2.50980000    1.38390000    0.00020000
C     -1.74730000    0.04600000    0.00050000
C     -0.20740000    0.03740000    0.00150000
C      0.57000000    1.36670000    0.00230000
C     -0.19250000    2.70460000    0.00120000
C     -1.73240000    2.71320000    0.00030000
Br     2.48000000    1.35620000    0.00470000
H     -3.64970000    1.39020000   -0.00030000
H     -2.32300000   -0.93800000   -0.00010000
H      0.35700000   -0.95310000    0.00170000
H      0.38300000    3.68860000    0.00170000
H     -2.29680000    3.70360000   -0.00030000
*


高斯的输入文件
%chk=BenBr-Gau-3.chk
%mem=100GB
%nproc=20
# B3LYP/Def2SVP em=gd3bj opt freq nosymm field=x+100

BenBr-Gau-3

0 1
C                 -2.50980000    1.38390000    0.00020000
C                 -1.74730000    0.04600000    0.00050000
C                 -0.20740000    0.03740000    0.00150000
C                  0.57000000    1.36670000    0.00230000
C                 -0.19250000    2.70460000    0.00120000
C                 -1.73240000    2.71320000    0.00030000
Br                 2.48000000    1.35620000    0.00470000
H                 -3.64970000    1.39020000   -0.00030000
H                 -2.32300000   -0.93800000   -0.00010000
H                  0.35700000   -0.95310000    0.00170000
H                  0.38300000    3.68860000    0.00170000
H                 -2.29680000    3.70360000   -0.00030000

wzkchem5

试试用orca算数值频率？

ljc050512

wzkchem5 发表于 2023-10-27 16:58
试试用orca算数值频率？

谢谢老师指导。算数值频率就正常了。

沁溪听友

wzkchem5 发表于 2023-10-27 16:58
试试用orca算数值频率？

老师好！我也遇到了相同的情况，也是像您说的算数值频率就好了。想请教一下算数值频率和默认的freq算出来的频率有什么区别呢？由它们分别得到的热力学校正值是否都正确呢？

wzkchem5

沁溪听友 发表于 2025-2-19 02:49
老师好！我也遇到了相同的情况，也是像您说的算数值频率就好了。想请教一下算数值频率和默认的freq算出来 ...

数值频率有一定的数值误差，但如果相关的收敛限都足够严的话，数值误差可以做到对于高频振动小于1cm-1，低频振动3cm-1左右。
两者的热力学校正都可用，但不宜混用

hxmjj

请问orca里面0.01a.u.的电场相当于多少V/m呢？

wzkchem5

hxmjj 发表于 2025-4-29 01:27
请问orca里面0.01a.u.的电场相当于多少V/m呢？

ORCA手册讲加电场的章节里说了

imasen

https://en.wikipedia.org/wiki/Atomic_units
https://physics.nist.gov/cgi-bin/cuu/Value?auefld

沁溪听友

wzkchem5 发表于 2025-2-19 15:12
数值频率有一定的数值误差，但如果相关的收敛限都足够严的话，数值误差可以做到对于高频振动小于1cm-1， ...

好的，谢谢老师！