计算化学公社

标题: 求助优化结构算freq有虚频 [打印本页]

作者
Author: 楠楠 时间: 2022-10-25 08:48
标题: 求助优化结构算freq有虚频
请教一下各位我基于优化好的结构算freq 命令行是 #p freq b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj 但是如图有负值；看了帖子http://sobereva.com/278 采取了在算opt时加上 int=ultrafine 命令行是 #p opt b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj int=ultrafine（g09计算）然后再算freq 命令行是#p freq b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj 仍有虚频；又采取#p b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj opt=calcall 但是报错了我计算的分子结构如下请教一下各位我该如何处理呢？

作者
Author: mfdsrax2 时间: 2022-10-25 09:25
不要说“但是报错了”这种话，报了什么错，能不能解决

作者
Author: 楠楠 时间: 2022-10-25 09:44

mfdsrax2 发表于 2022-10-25 09:25
不要说“但是报错了”这种话，报了什么错，能不能解决

报告老师命令行是这个 #p opt b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersio
n=gd3bj opt=calcall
报错信息是这个
Charge and Multiplicity card seems defective:
Wanted an integer as input.

?
Error termination via Lnk1e in /share/apps/g09.d01/g09/l101.exe at Mon Oct 24 14:06:33 2022.
Job cpu time:    0 days  2 hours 56 minutes 32.6 seconds.
File lengths (MBytes):  RWF= 2196 Int=    0 D2E=    0 Chk=    81 Scr=    1

作者
Author: mfdsrax2 时间: 2022-10-25 09:57

楠楠发表于 2022-10-25 09:44
报告老师命令行是这个 #p opt b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersio
n ...

Charge and Multiplicity card seems defective:
Wanted an integer as input.

这个应该看的懂吧，电荷和自选多重度这行有问题，你看看输入文件

作者
Author: chands 时间: 2022-10-25 09:59

楠楠发表于 2022-10-25 09:44
报告老师命令行是这个 #p opt b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersio
n ...

这个估计是输入格式的问题，在该读取电荷数和自旋多重度的地方出现其他字符，于是报错。

作者
Author: 楠楠 时间: 2022-10-25 10:03

mfdsrax2 发表于 2022-10-25 09:57
Charge and Multiplicity card seems defective:
Wanted an integer as input.

明白老师这是我的输入文件可是我就是基于这个文件之前算都是正确的加上这个关键词以后就这样报错我的体系就是中性态自选多重度为1
%chk=F.chk
%nprocshared=24
%mem=24GB
#p opt b3lyp/genecp scrf=(solvent=acetonitrile,pcm)  empiricaldispersion=gd3bj  opt=calcall

opt

0 1
O                -5.21881700 2.39553600 -2.44429400
~~~
Br             -1.41050099 4.20727403 0.23931682

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

Br 0
LANL2DZ

作者
Author: 北大-陶豫 时间: 2022-10-25 10:07

楠楠发表于 2022-10-25 10:03
明白老师这是我的输入文件可是我就是基于这个文件之前算都是正确的加上这个关键词以后就这样报错我的 ...

1）端正态度，别写错别字。
2）为什么 opt 之后又来个 opt=CalcAll，正常来说优化到极小点是不需要 CalcAll 的，直接 opt 就行。
3）高斯算周期性体系不能算频率，不能 freq 也不能 CalcAll，不知道你怎么算的。

作者
Author: 乐平 时间: 2022-10-25 10:09

楠楠发表于 2022-10-25 09:44
报告老师命令行是这个 #p opt b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersio
n ...

首先，几何优化和频率分析可以写在一起，Gaussian 会先做几何优化，结束之后自动根据几何优化后的构型做频率分析。不需要你分开算两次。

其次，即使你分开算，也最好用 geom=check guess=read 读取几何优化得到的 .chk 文件里的信息。比你自己重新保存优化后的结构，重新计算 Hessian 矩阵稍稍快一些。

再次，你目前的报错很可能是自己重新写输入文件时不仔细，把电荷，自旋多重度那一行的位置弄错了。

不过很奇怪，L101 是检查输入文件，应该很快就报错，不至于 2 小时 56 分之后才报错啊

作者
Author: 楠楠 时间: 2022-10-25 10:10

北大-陶豫发表于 2022-10-25 10:07
1）端正态度，别写错别字。
2）为什么 opt 之后又来个 opt=CalcAll，正常来说优化到极小点是不需要 Calc ...

好的老师；
那我将命令行改为#p b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj opt=calcall 再计算一下；
老师我将n=1了没有视为周期性体系算

作者
Author: 楠楠 时间: 2022-10-25 10:37

乐平发表于 2022-10-25 10:09
首先，几何优化和频率分析可以写在一起，Gaussian 会先做几何优化，结束之后自动根据几何优化后的构型做 ...

老师您的意思是opt和freq一起算是吗？那我的命令行是写成这样吗？#p freq b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj opt=calcall

作者
Author: LittlePupil 时间: 2022-10-25 23:47
当且仅当opt加入“calcall”选项收敛以后程序会自动输出振动频率和热力学校正量，此时不需要再写"freq"（写了只会白白浪费机时）。

作者
Author: 乐平 时间: 2022-10-26 14:50
本帖最后由乐平于 2022-10-26 14:54 编辑

楠楠发表于 2022-10-25 10:37
老师您的意思是opt和freq一起算是吗？那我的命令行是写成这样吗？#p freq b3lyp/genecp scrf=(solvent=ac ...

不是你最后 opt=calcall 和 freq 一起算…… opt=calcall 不用加 freq，最后会自动算频率
我的“首先”说的是你最开始的时候 opt freq 可以写在一条命令里

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