计算化学公社

标题: 有虚频该如何消除 [打印本页]

作者
Author: 兰兰爱学习 时间: 2022-10-29 15:21
标题: 有虚频该如何消除
想请教一下各位，我对如下分子优化结果正常。然后算freq有虚频看了帖子尝试过在做opt优化时加上关键词int=ultrafine int=superfine  opt=tight freq 然后算freq还是有虚频；然后尝试opt=calcall体系震荡了，于是又加上关键词opt=calcall=GDIIS 还是不行。又将pcm换成IEFPCM也不行。考虑到是否是对称性问题加上了关键词nosymm 也还是不行。以下是我的部分输入文件，请问老师们还有什么方法可以一试的吗？我看了一下模型主要是侧链部分带动整体震荡%chk=BR.chk
%nprocshared=24
%mem=24GB
#p OPT  b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj

opt

0 1
O                4.81864200 -4.11494100 0.23683200
O                7.90368500 -0.81375600 -0.33717300
N                6.36347800 -2.46848300 -0.08886300
C                5.06300500 -2.96404900 -0.11112100

----
N                -7.03904400 -2.78213700 -2.77854700
Br             -1.70322600 7.22060400 1.74135000
Br             -6.92202900 -4.76119900 -2.78220300

C H O N S
6-31G*
****
Br 0
LANL2DZ
****

Br 0
LANL2DZ

作者
Author: wzkchem5 时间: 2022-10-29 15:29
你这里Br是负离子吧？如果是的话必须加弥散。否则你即使虚频消掉了结果也不能用。

作者
Author: 兰兰爱学习 时间: 2022-10-29 15:39
本帖最后由兰兰爱学习于 2022-10-29 16:33 编辑

wzkchem5 发表于 2022-10-29 15:29
你这里Br是负离子吧？如果是的话必须加弥散。否则你即使虚频消掉了结果也不能用。

老师是溴离子但体系整体是中性的也要考虑弥散吗？如果考虑弥散那学生将后面针对溴的赝势基组LANL2DZ改为LANL2DZdp可以吗？

作者
Author: niobium 时间: 2022-10-29 16:26

兰兰爱学习发表于 2022-10-29 15:39
老师那学生将后面针对溴的赝势基组LANL2DZ改为LANL2DZdp可以吗？

可以用ma-SVP

作者
Author: wzkchem5 时间: 2022-10-29 16:35

兰兰爱学习发表于 2022-10-29 08:39
老师那学生将后面针对溴的赝势基组LANL2DZ改为LANL2DZdp可以吗？

也许可以，但我觉得还是用ma-def2-SVP之类的全电子基组比较好。这个体系只有两个Br，用赝势省不了什么计算时间，还引入不小的误差，恐怕得不偿失

作者
Author: sobereva 时间: 2022-10-29 16:46

兰兰爱学习发表于 2022-10-29 15:39
老师是溴离子但体系整体是中性的也要考虑弥散吗？如果考虑弥散那学生将后面针对溴的赝势基组LANL2DZ改 ...

弥散函数是对带显著负电荷的原子加的，不是看整体是不是带负电

作者
Author: 兰兰爱学习 时间: 2022-10-29 16:47

wzkchem5 发表于 2022-10-29 16:35
也许可以，但我觉得还是用ma-def2-SVP之类的全电子基组比较好。这个体系只有两个Br，用赝势省不了什么计 ...

老师再请教一下您的意思是对于我这个体系因为溴只有两个且为第四周期所以不用赝势也可以是吗？那我全部都用6-31+G*基组您看可以吗？因为目前只是做优化想着考虑机时问题。

作者
Author: sobereva 时间: 2022-10-29 16:49
此问题跟nosymm毫无关系
根本没有opt=calcall=GDIIS这种写法
int=ultrafine int=superfine opt=tight freq 莫名其妙，int怎么能写两次
把问题表述具体，不要有任何歧义。

PCM和IEFPCM关键词完全是一码事，写PCM完全多此一举

Br改成ma-SVP或者6-31+G*，关键词用b3lyp/gen scrf=solvent=acetonitrile em=gd3bj opt(recalc=3,gdiis) freq。如果是G09，再带上int=ultrafine。若还有虚频或者震荡了，把输出文件压缩后上传。

有效的解决虚频的办法下文都说了，没有可补充的，也没有按博文的做法消不掉的，没有例外。没提的别瞎试
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

作者
Author: 兰兰爱学习 时间: 2022-10-29 16:50

sobereva 发表于 2022-10-29 16:49
此问题跟nosymm毫无关系
根本没有opt=calcall=GDIIS这种写法
int=ultrafine int=superfine opt=tight fr ...

好的老师谢谢

作者
Author: 兰兰爱学习 时间: 2022-11-1 08:36

sobereva 发表于 2022-10-29 16:49
此问题跟nosymm毫无关系
根本没有opt=calcall=GDIIS这种写法
int=ultrafine int=superfine opt=tight fr ...

老师按照关键词计算了一下貌似震荡了麻烦您给看一下我看优化过程整个分子都在震动

作者
Author: sobereva 时间: 2022-11-2 08:13

兰兰爱学习发表于 2022-11-1 08:36
老师按照关键词计算了一下貌似震荡了麻烦您给看一下我看优化过程整个分子都在震动

把N131和Br133摆远点再优化，理应不该成键。另外，Br-应该摆在氨基氢的方向，因为氨基氢带明显正电

如果还这样，opt加上maxstep=3,notrust再试

作者
Author: 兰兰爱学习 时间: 2023-4-24 19:30
体系依赖性问题，请管理员删除帖子，谢谢

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