Gaussian 常见报错及解决方法；新手求助报错时的注意事项

fly_ylx

非常好，对问题很有帮助

菜鸟一枚

请问L308的错误是什么原因？ Error termination via Lnk1e in /share/apps/gaussian/g09D01/g09/l308.exe at Sun Mar 12 04:13:09 2017.应该怎么解决？谢谢老师

liyuanhe211

菜鸟一枚 发表于 2017-3-11 14:57
请问L308的错误是什么原因？ Error termination via Lnk1e in /share/apps/gaussian/g09D01/g09/l308.exe a ...

根本没看帖子

liyuanhe211

增加对柔性扫描找过渡态的说明

冰释之川

liyuanhe211 发表于 2017-3-11 20:57
增加对柔性扫描找过渡态的说明

噫。。H原子在C原子和S原子之间徘徊的动画好鬼畜……

冰释之川

补充一个报错，有些方法无法和polar关键词联用，比如LC-wPBE在G09.D01：
InvSVY:  IOpt=1 It=  1 EMax= 6.66D-16
Solved reduced A of dimension    24 with     3 vectors.
FullF1:  Do perturbations    1 to     3.
SCF Polarizability for W=    0.000000:
                1             2             3
      1  0.125599D+02
      2  0.000000D+00  0.125599D+02
      3  0.000000D+00  0.000000D+00  0.125599D+02
Isotropic polarizability for W=    0.000000       12.56 Bohr**3.
No func 3rd derivs with HSE.
Error termination via Lnk1e in D:\Program Files (x86)\g09\G09W\l1002.exe at Sun Mar 26 11:50:16 2017.
Job cpu time:       0 days  0 hours  0 minutes  5.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

liyuanhe211

冰释之川 发表于 2017-3-26 12:02
补充一个报错，有些方法无法和polar关键词联用，比如LC-wPBE在G09.D01：
InvSVY:  IOpt=1 It=  1 EMax= 6. ...

这是因为泛函没有三阶解析导数，应当更换有三阶导数的泛函吗？

冰释之川

liyuanhe211 发表于 2017-3-26 14:45
这是因为泛函没有三阶解析导数，应当更换有三阶导数的泛函吗？

YES，不知道G16支持不支持LC-BLYP也不行

liyuanhe211

更新：增加：
            · L1002，No func 3rd derivs with XXX.

以及增加了一些极其基础的问题：
            · L101，Wanted a XXX as input. Found an XXX as input.
            · L114，ERROR IN INITNF. NUMBER OF VARIABLES (  0) INCORRECT (SHOULD BE BETWEEN 1 AND 50)
            · L301，The combination of multiplicity X and XXX electrons is impossible.
            · L301，Atomic number out of range for XXX basis set.

菜鸟一枚

老师好，我想问一下L308错误是什么原因？应该怎么解决呢？ Error termination via Lnk1e in /share/apps/gaussian/g09D01/g09/l308.exe

liyuanhe211

菜鸟一枚 发表于 2017-3-27 16:25
老师好，我想问一下L308错误是什么原因？应该怎么解决呢？ Error termination via Lnk1e in /share/apps/ga ...

先把帖子前半部分读了。

冰释之川

关于“No func 3rd derivs with HSE.”进一步研究。
对Gaussian 09 D01版里的LC-wPBE 先吐槽一下，LC-wPBE居然只有解析的一阶导数！！！
想要使用LC-wPBE算极化率，超极化率的，只能用数值差分了：

(1) 算数值的极化率（二阶导数）：polar=Numerical LC-wPBE
(2) 算数值的极化率和第一超级化率（二阶、三阶导数）：polar=DoubleNumer LC-wPBE 【此时要两次差分，可能会导致一定误差，建议测试差分步长（默认步长是0.0003 au）对数值的稳定性】
(3) 算数值的第二超极化率（四阶导数）：连解析的二阶导都没，无法通过通过情况(2)再做差分。高斯里有“Gamma”关键词但要求必须有解析的二阶导。

另外，新发现：TD LC-wPBE 不能和density联用 （别想获得该泛函下的激发态密度了，没辙）【不过G16里的新泛函LC-wHPBE，即LC-wPBE的姊妹泛函，可以保存算激发态密度】

liyuanhe211

冰释之川 发表于 2017-4-12 11:27
关于“No func 3rd derivs with HSE.”进一步研究。
对Gaussian 09 D01版里的LC-wPBE 先吐槽一下，LC-wPBE ...

感谢冰老师补充，已加入原文里
一直不知道density计算的本质是什么 (⊙_⊙)

sobereva

冰释之川 发表于 2017-4-12 11:27
关于“No func 3rd derivs with HSE.”进一步研究。
对Gaussian 09 D01版里的LC-wPBE 先吐槽一下，LC-wPBE ...

LC-wPBE有二阶解析导数，polar关键词报错的时候实际上已经输出了极化率了。而且可以正常做freq，只有用freq=raman这样需要三阶导数的量才报错。

TD计算在构建TD矩阵元时本身就需要交换-相关密度的二阶导数。density关键词产生密度又需要有一阶解析导数，所以LC-wPBE这样交换-相关密度没有三阶解析导数的泛函没法TD density一起用。

冰释之川

liyuanhe211 发表于 2017-4-12 23:44
感谢冰老师补充，已加入原文里
一直不知道density计算的本质是什么 (⊙_⊙)

李老师，麻烦根据Sob老师的补充再修正一下：

LC-wPBE 有二阶解析导数

（1）【A】算极化率 用polar=Numerical 程序正常跑完，会给出α的数值解
【B】算解析的极化率 用polar默认关键词，则会因LC-wPBE没有三阶导数而无法求得解析β值而报错，但实际上解析的极化率α已经求出，并保存在log文件中。
InvSVY:  IOpt=1 It=  1 EMax= 4.44D-16
Solved reduced A of dimension    31 with     3 vectors.
FullF1:  Do perturbations    1 to     3.
SCF Polarizability for W=    0.000000:
                1             2             3
      1  0.311786D+02
      2  0.214915D+01  0.273009D+02
      3 -0.264577D-04  0.351431D-04  0.193900D+02
Isotropic polarizability for W=    0.000000       25.96 Bohr**3.
No func 3rd derivs with HSE.
Error termination via Lnk1e in D:\Program Files (x86)\G09\G09W\l1002.exe at Thu Apr 13 08:57:29 2017.
Job cpu time:       0 days  0 hours  0 minutes 14.0 seconds.