计算化学公社

标题: Gaussian 从l1002 restart 频率计算未输出热力学量 [打印本页]

作者
Author: liyuanhe211 时间: 2015-10-10 13:35
标题: Gaussian 从l1002 restart 频率计算未输出热力学量
在Linux上跑一个4天的freq计算 [M06-2X-D3/def2-TZVP//B3LYP-D3/6-31+G(d,p)]的第二步
中间停电断掉之后又从rwf文件restart了，之后正常结束。
之后restart从l1002.exe开始，正常结束（.out文件见附件），但其中未包含想要的热力学信息，其他有些输出也不正常。
而且由于写了%nosave，RWF文件在正常结束后被删除了，.chk文件在restart之后也没有写入新的东西。

另外有些地方的输出似乎也不正常

Variable Old X -DE/DX Delta X Delta X Delta X New X
(Linear) (Quad) (Total)
R1 2.91527 -0.00745 0.00000********** -0.01177 2.90350
R2 2.92820 -0.00640 0.00000********** -0.01010 2.91811
R3 2.91304 -0.00732 0.00000********** -0.01155 2.90149
R4 2.69313 -0.01257 0.00000********** -0.01986 2.67328
R5 2.91116 -0.00618 0.00000********** -0.00976 2.90140
R6 2.93307 -0.00759 0.00000********** -0.01198 2.92110

复制代码

现在是否可能在不重算的情况下把热力学信息提出来？或者以后如何避免频率计算Restart的时候出问题。

作者
Author: sobereva 时间: 2015-10-10 17:09
用freqchk看看能否从chk得到热力学数据，如果不行的话就重算吧。

作者
Author: liyuanhe211 时间: 2015-10-11 18:08

sobereva 发表于 2015-10-10 17:09
用freqchk看看能否从chk得到热力学数据，如果不行的话就重算吧。

用FreqChk确实输出了热力学量的修正量。只不过没有加上E_ele之后的总能量，用l502的输出算一下就行了。十分感谢。

Zero-point correction= 0.318762 (Hartree/Particle)
Thermal correction to Energy= 0.336863
Thermal correction to Enthalpy= 0.337808
Thermal correction to Gibbs Free Energy= 0.272535

复制代码

作者
Author: 978142355 时间: 2015-10-11 19:34

liyuanhe211 发表于 2015-10-11 18:08
用FreqChk确实输出了热力学量的修正量。只不过没有加上E_ele之后的总能量，用l502的输出算一下就行了。十 ...

问一下，freqchk命令怎么输入啊？在终端里直接输入：
freqchk *.chk，是这样嘛？

作者
Author: liyuanhe211 时间: 2015-10-11 19:58

978142355 发表于 2015-10-11 19:34
问一下，freqchk命令怎么输入啊？在终端里直接输入：
freqchk *.chk，是这样嘛？

http://sobereva.com/221 《谈谈谐振频率校正因子》

作者
Author: 978142355 时间: 2015-10-11 20:32

liyuanhe211 发表于 2015-10-11 19:58
http://sobereva.com/221 《谈谈谐振频率校正因子》

懂了，谢谢，另外也把Gaussian09的手册中的也粘贴了出来，方便后面的人对比学习。
http://www.gaussian.com/g_tech/g_ur/u_freqchk.htm

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