计算化学公社

标题: 频率分析遇到的问题-----NaN [打印本页]

作者
Author: nkallwar 时间: 2016-7-27 18:29
标题: 频率分析遇到的问题-----NaN

前面的优化是正常的四个Yes结束，做频率分析中出现：

Full mass-weighted force constant matrix:
Low frequencies ---    NaN    NaN    NaN    NaN    NaN    NaN
Low frequencies ---    NaN    NaN    NaN
Diagonal vibrational polarizability:
            NaN          NaN          NaN
Harmonic frequencies (cm**-1), IR intensities (KM/Mole), Raman scattering
activities (A**4/AMU), depolarization ratios for plane and unpolarized
incident light, reduced masses (AMU), force constants (mDyne/A),
and normal coordinates:
                  1                   2                   3
                  A                   A                   A
Frequencies --       NaN                   NaN                   NaN
Red. masses --       NaN                   NaN                   NaN
Frc consts  --       NaN                   NaN                   NaN
IR Inten --       NaN                   NaN                   NaN

中间的Matrix也都是NaN

最后得到的correction是这样的：
Zero-point correction=                               NaN (Hartree/Particle)
Thermal correction to Energy=                   0.002833
Thermal correction to Enthalpy=                0.003777
Thermal correction to Gibbs Free Energy=       -0.147324
Sum of electronic and zero-point Energies=                NaN
Sum of electronic and thermal Energies=          -1977.250562
Sum of electronic and thermal Enthalpies=          -1977.249618
Sum of electronic and thermal Free Energies=       -1977.400718

最后的报错信息：
Maximum of*** iterations exceeded in RedStp.
Error termination via Lnk1e in /home-yw/Soft/Gaussian.D01.Linda//g09/l103.exe at Wed Jul 27 17:13:44 2016.
Job cpu time:    0 days  5 hours 54 minutes 30.0 seconds.
File lengths (MBytes):  RWF= 1156 Int=    0 D2E=    0 Chk=    23 Scr=    1

不知道是何种原因导致？

作者
Author: sobereva 时间: 2016-7-27 22:07
贴出频率计算用的关键词。最好直接上传压缩后的输出文件

作者
Author: nkallwar 时间: 2016-7-27 22:23

sobereva 发表于 2016-7-27 22:07
贴出频率计算用的关键词。最好直接上传压缩后的输出文件

输入文件：
%nprocshared=12
%mem=40GB
%chk=TPPNTf.chk
#p freq b3lyp/6-31g(d) scrf=(iefpcm,solvent=toluene)

ldc

0 1

P                0.53259300 -0.02497000 -0.11540700
C                1.58864200 -1.35036000 -0.77881400
C                1.05376200 -2.26227000 -1.70006900
C                2.92394400 -1.47445400 -0.36578800
C                1.85079500 -3.29571500 -2.19429600
......

输出文件见附件：

作者
Author: sobereva 时间: 2016-7-27 23:01
尝试E.01版，还不行就把输出文件直接发Gaussian中文客服的邮箱里问，是明显的bug。

作者
Author: nkallwar 时间: 2016-7-28 17:06

sobereva 发表于 2016-7-27 23:01
尝试E.01版，还不行就把输出文件直接发Gaussian中文客服的邮箱里问，是明显的bug。

没有E.01, 换成B.01算了一下过去了，看来是个bug

作者
Author: hlmkh 时间: 2016-8-31 15:56
Hi
I've got the same problem solution finding friends?

作者
Author: MrD 时间: 2017-5-10 09:52
我也有这个问题，用d01和e01都不行，也尝试换泛函，换基组，增大积分格点，存出坐标直接freq，能尝试的方法几乎都尝试了，无药可救。刚已经给高斯技术发了信，看他作何解释。

作者
Author: ananbaby2309 时间: 2017-6-8 21:22

MrD 发表于 2017-5-10 09:52
我也有这个问题，用d01和e01都不行，也尝试换泛函，换基组，增大积分格点，存出坐标直接freq，能尝试的方法 ...

我最近也遇到了这个问题，请问高斯客服是怎么回复您的呢？或者说您最后是怎么解决这个问题的，能给我说一下吗？

作者
Author: sobereva 时间: 2017-6-28 18:20

ananbaby2309 发表于 2017-6-8 21:22
我最近也遇到了这个问题，请问高斯客服是怎么回复您的呢？或者说您最后是怎么解决这个问题的，能给我说一 ...

如果%chk引用了之前的chk文件，把那个chk先删掉再试。有时候之前的chk会捣蛋引发此现象

作者
Author: ananbaby2309 时间: 2017-11-16 15:58

sobereva 发表于 2017-6-28 18:20
如果%chk引用了之前的chk文件，把那个chk先删掉再试。有时候之前的chk会捣蛋引发此现象

好的，谢谢您

作者
Author: 978142355 时间: 2017-11-16 18:01
我师兄也遇到了这个问题，咨询了高斯客服，他们承认是个bug，让用Gaussian16重新算。

作者
Author: qwe1832518773 时间: 2018-9-22 08:41
我用的G16也出现了这个BUG。。。几率事件吗？修改精度会有用吗

作者
Author: 后琪 时间: 2019-2-21 16:56
请问大家这个问题解决了吗？
是此类计算常出现这个bug吗？

作者
Author: guyizhang 时间: 2024-1-30 11:14

nkallwar 发表于 2016-7-27 22:23
输入文件：
%nprocshared=12
%mem=40GB

问题解决了吗，会是内存不足的问题吗，我在使用vasp的时候也遇到了同样的问题

作者
Author: 北大-陶豫 时间: 2024-1-30 14:31
我上次遇到这种问题，是因为输入文件里chk文件名忘记改了，导致两个任务同时在运行并且在读写同一个chk文件。你可以检查一下自己有没有类似的问题。

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