求助，计算产物单点能总是终止计算。是因为要再定义epsinf吗？

Psykk · 发表于 Post on 2025-7-29 09:45:04

各位老师好，我在计算产物（简单的有机小分子）的单点能，反应体系是120℃，溶剂是均三甲苯，查阅高斯手册均三甲苯Mesitylene: eps=2.2650。所以我直接在输入文件最后定义了eps=2.2650。但是计算终止了。但是我用同样的输入文件计算中间体的单点能却没有问题。请问大家这是什么原因？是因为要再定义epsinf吗？附件是我产物的输入和输出文件，还有过渡态的输入和输出文件。

产物-P-Linear-sp.gjf (3.9 KB, 下载次数 Times of downloads: 0)

产物-P-Linear-sp.log (18.97 KB, 下载次数 Times of downloads: 0)

过渡态-TS-2aa-sp.gjf (12.16 KB, 下载次数 Times of downloads: 0)

过渡态-TS-2aa-sp.log (148.16 KB, 下载次数 Times of downloads: 0)

北大-陶豫 · 发表于 Post on 2025-7-29 10:21:48

这不说得很明确了吗？你没定义epsinf
一个输入文件示例如下：

%rwf=CO2.rwf
%NoSave
%chk=CO2.chk
%nProcShared=8
%Mem=30GB
# opt freq M062X/6-311+G** SCRF(read) Guess=Read Geom=AllCheck
eps=20.436
epsinf=2.002

复制代码

Psykk · 发表于 Post on 2025-7-29 10:26:14

北大-陶豫发表于 2025-7-29 10:21
这不说得很明确了吗？你没定义epsinf
一个输入文件示例如下：

好的，谢谢老师，我计算过渡态能量也没定义epsinf，只是定义了eps, 但是计算正常结束了。为了计算统一，是不是计算过渡态能量也得定义epsinf？

北大-陶豫 · 发表于 Post on 2025-7-29 10:35:35

Psykk 发表于 2025-7-29 10:26
好的，谢谢老师，我计算过渡态能量也没定义epsinf，只是定义了eps, 但是计算正常结束了。为了计算统一， ...

是的

如果只是做一般基态计算，只需要定义eps就够了。如果做TDDFT计算、含频(超)极化率计算之类，还得同时定义epsinf，这会影响结果。而振动分析、NMR等计算也得定义epsinf，但数值可以随意设，不影响结果。epsinf一般查不到，可以用折射率的平方估算，折射率比较好查。如果连折射率也没有，那么对于非极性溶剂，可以假定epsinf的数值等于eps；对于极性溶剂，epsinf可以姑且设1.9（各种溶剂的epsinf都在这附近）。

出处：http://bbs.keinsci.com/thread-3345-1-1.html

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 求助，计算产物单点能总是终止计算。是因为要再定义epsinf吗？

浏览过的版块