Gaussian 从l1002 restart 频率计算未输出热力学量

liyuanhe211

在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之后也没有写入新的东西。

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

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

复制代码

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

sobereva

用freqchk看看能否从chk得到热力学数据，如果不行的话就重算吧。

liyuanhe211

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

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

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

复制代码

978142355

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

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

liyuanhe211

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

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

978142355

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

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