计算化学公社

标题: 求助，计算产物单点能总是终止计算。是因为要再定义epsinf吗？ [打印本页]
作者
Author: Psykk    时间: 2025-7-29 09:45
标题: 求助，计算产物单点能总是终止计算。是因为要再定义epsinf吗？
各位老师好，我在计算产物（简单的有机小分子）的单点能，反应体系是120℃，溶剂是均三甲苯，查阅高斯手册均三甲苯Mesitylene: eps=2.2650。所以我直接在输入文件最后定义了eps=2.2650。但是计算终止了。但是我用同样的输入文件计算中间体的单点能却没有问题。请问大家这是什么原因？是因为要再定义epsinf吗？附件是我产物的输入和输出文件，还有过渡态的输入和输出文件。 (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)


作者
Author: 北大-陶豫    时间: 2025-7-29 10:21
这不说得很明确了吗？你没定义epsinf
一个输入文件示例如下：
  1. %rwf=CO2.rwf
  2. %NoSave
  3. %chk=CO2.chk
  4. %nProcShared=8
  5. %Mem=30GB
  6. # opt freq M062X/6-311+G** SCRF(read) Guess=Read Geom=AllCheck

  8. eps=20.436
  9. epsinf=2.002


复制代码
作者
Author: Psykk    时间: 2025-7-29 10:26
北大-陶豫 发表于 2025-7-29 10:21
这不说得很明确了吗？你没定义epsinf
一个输入文件示例如下：

好的，谢谢老师，我计算过渡态能量也没定义epsinf，只是定义了eps, 但是计算正常结束了。 为了计算统一，是不是计算过渡态能量也得定义epsinf？
作者
Author: 北大-陶豫    时间: 2025-7-29 10: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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3