计算化学公社

标题: ORCA加电场对分子优化时出现诡异虚频 [打印本页]

作者
Author:
ljc050512    时间: 2023-10-27 16:12
标题: ORCA加电场对分子优化时出现诡异虚频
各位老师,我在对溴苯进行加电场(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




作者
Author:
wzkchem5    时间: 2023-10-27 16:58
试试用orca算数值频率?
作者
Author:
ljc050512    时间: 2023-10-27 19:50
wzkchem5 发表于 2023-10-27 16:58
试试用orca算数值频率?

谢谢老师指导。算数值频率就正常了。
作者
Author:
沁溪听友    时间: 2025-2-19 09:49
wzkchem5 发表于 2023-10-27 16:58
试试用orca算数值频率?

老师好!我也遇到了相同的情况,也是像您说的算数值频率就好了。想请教一下算数值频率和默认的freq算出来的频率有什么区别呢?由它们分别得到的热力学校正值是否都正确呢?
作者
Author:
wzkchem5    时间: 2025-2-19 15:12
沁溪听友 发表于 2025-2-19 02:49
老师好!我也遇到了相同的情况,也是像您说的算数值频率就好了。想请教一下算数值频率和默认的freq算出来 ...

数值频率有一定的数值误差,但如果相关的收敛限都足够严的话,数值误差可以做到对于高频振动小于1cm-1,低频振动3cm-1左右。
两者的热力学校正都可用,但不宜混用
作者
Author:
hxmjj    时间: 2025-4-29 01:27
请问orca里面0.01a.u.的电场相当于多少V/m呢?
作者
Author:
wzkchem5    时间: 2025-4-29 08:55
hxmjj 发表于 2025-4-29 01:27
请问orca里面0.01a.u.的电场相当于多少V/m呢?

ORCA手册讲加电场的章节里说了
作者
Author:
imasen    时间: 2025-4-29 09:52
https://en.wikipedia.org/wiki/Atomic_units
https://physics.nist.gov/cgi-bin/cuu/Value?auefld
作者
Author:
沁溪听友    时间: 2025-5-14 11:44
wzkchem5 发表于 2025-2-19 15:12
数值频率有一定的数值误差,但如果相关的收敛限都足够严的话,数值误差可以做到对于高频振动小于1cm-1, ...

好的,谢谢老师!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3