计算化学公社

 找回密码 Forget password
 注册 Register

求助优化结构算freq有虚频

查看数: 2607 | 评论数: 11 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-10-25 08:48

正文摘要:

请教一下各位 我基于优化好的结构算freq 命令行是 #p freq b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj   但是如图有负值;看了帖子http://sobereva.com/278 采取了在算opt ...

回复 Reply

乐平 发表于 Post on 2022-10-26 14:50:22
本帖最后由 乐平 于 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 可以写在一条命令里

LittlePupil 发表于 Post on 2022-10-25 23:47:56
当且仅当opt加入“calcall”选项收敛以后程序会自动输出振动频率和热力学校正量,此时不需要再写"freq"(写了只会白白浪费机时)。
楠楠 发表于 Post on 2022-10-25 10:37:23
乐平 发表于 2022-10-25 10:09
首先,几何优化和频率分析可以写在一起,Gaussian 会先做几何优化,结束之后自动根据几何优化后的构型做 ...

老师您的意思是opt和freq一起算是吗?那我的命令行是写成这样吗?#p freq b3lyp/genecp scrf=(solvent=acetonitrile,pcm) empiricaldispersion=gd3bj  opt=calcall
楠楠 发表于 Post on 2022-10-25 10:10:14
北大-陶豫 发表于 2022-10-25 10:07
1)端正态度,别写错别字。
2)为什么 opt 之后又来个 opt=CalcAll,正常来说优化到极小点是不需要 Calc ...

好的老师;
那我将命令行改为#p  b3lyp/genecp scrf=(solvent=acetonitrile,pcm)  empiricaldispersion=gd3bj  opt=calcall 再计算一下;
老师我将n=1了  没有视为周期性体系算
乐平 发表于 Post on 2022-10-25 10:09:48
楠楠 发表于 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 分之后才报错啊
北大-陶豫 发表于 Post on 2022-10-25 10:07:06
楠楠 发表于 2022-10-25 10:03
明白老师 这是我的输入文件 可是我就是基于这个文件之前算都是正确的 加上这个关键词以后就这样报错 我的 ...

1)端正态度,别写错别字。
2)为什么 opt 之后又来个 opt=CalcAll,正常来说优化到极小点是不需要 CalcAll 的,直接 opt 就行。
3)高斯算周期性体系不能算频率,不能 freq 也不能 CalcAll,不知道你怎么算的。
楠楠 发表于 Post on 2022-10-25 10:03:00
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


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

这个估计是输入格式的问题,在该读取电荷数和自旋多重度的地方出现其他字符,于是报错。
mfdsrax2 发表于 Post on 2022-10-25 09:57:20
楠楠 发表于 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.

这个应该看的懂吧,电荷和自选多重度这行有问题,你看看输入文件
楠楠 发表于 Post on 2022-10-25 09:44:16
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 发表于 Post on 2022-10-25 09:25:24
不要说“但是报错了”这种话,报了什么错,能不能解决

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-3-1 05:09 , Processed in 0.219952 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list