求助优化结构算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  但是报错了 我计算的分子结构如下 请教一下各位我该如何处理呢？

mfdsrax2

不要说“但是报错了”这种话，报了什么错，能不能解决

楠楠

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

mfdsrax2

楠楠 发表于 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.

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

chands

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

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

楠楠

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

北大-陶豫

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

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

乐平

楠楠 发表于 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 分之后才报错啊

楠楠

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

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

楠楠

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

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

LittlePupil

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

乐平

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

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