频率分析遇到的问题-----NaN

nkallwar

前面的优化是正常的四个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


不知道是何种原因导致？

sobereva

贴出频率计算用的关键词。最好直接上传压缩后的输出文件

nkallwar

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
......

输出文件见附件：

sobereva

尝试E.01版，还不行就把输出文件直接发Gaussian中文客服的邮箱里问，是明显的bug。

nkallwar

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

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

hlmkh

Hi
I've got the same problem solution finding friends?

MrD

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

ananbaby2309

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

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

sobereva

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

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

ananbaby2309

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

好的，谢谢您

978142355

我师兄也遇到了这个问题，咨询了高斯客服，他们承认是个bug，让用Gaussian16重新算。

qwe1832518773

我用的G16也出现了这个BUG。。。几率事件吗？修改精度会有用吗

后琪

请问大家这个问题解决了吗？
是此类计算常出现这个bug吗？

guyizhang

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

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

北大-陶豫

我上次遇到这种问题，是因为输入文件里chk文件名忘记改了，导致两个任务同时在运行并且在读写同一个chk文件。你可以检查一下自己有没有类似的问题。