计算化学公社

标题: Gaussian 常见报错及解决方法；新手求助报错时的注意事项 [打印本页]

作者
Author: liyuanhe211 时间: 2017-1-1 14:22
标题: Gaussian 常见报错及解决方法；新手求助报错时的注意事项
本帖最后由 liyuanhe211 于 2022-9-9 15:02 编辑

本文在：(1) 完整转载 (2) 在文章内容开始前注明出处、作者 (3) 非商业用途（商业用途包含本身有财务收入、或依附于盈利机构的公众号）的前提下可以任意转载。

应一些人的要求，正在于本人GitHub撰写本文后半部分的英文版，但刚刚开始。如有人愿意贡献，可以去发Pull request（机翻就算了）。

有感于群里和论坛上常有日经、周经问题浪费答疑者大量时间，降低信息密度；以及有些提问方式令人甚为捉急，令提问者、答疑者均费时费力；撰写本帖。
本帖面向新手/新人，望能做一常见回复索引，仅求常用且正确。再出现类似问题时，可直接令其看“论坛 Gaussian FAQ”，不必重复回答。

各常见报错的解决方法见本帖的后半部分。本帖将持续更新，欢迎大家补充或纠正错误。

初来者请先阅读：

《在网上求助计算化学问题时的注意事项》

——————————————————————————————————————————————————————————————

遇到报错：

         · 提问前应做的工作
         · 提问时 Do / Do NOT
         · 提问后

Gaussian常见简单报错及解决方法 （行文为通览而使用了不同顺序，仅为解决某具体问题应用 Ctrl+F 关键词检索，定位至相应段落）：

         · Severe Error Message # 2070 (Windows)
         · segmentation violation/segmentation fault (Linux)

   · “量子化学问题”报错：

         · L1 ，ntrex1
         · L1 ，Illegal IType or MSType generated by parse.
         · L1 ，QPErr --- A syntax error was detected in the input line.
         · L101，End of file in ZSymb.
         · L101，WANTED A XXX AS INPUT. FOUND AN XXX AS INPUT.
         · L101，Input Error Input Error Input Error Input Error Input Error Input Error
         · L103，Error in internal coordinate system
         · L103，Linear angle in Bend; Linear angle in Tors
         · L103，FormBX had a problem
         · L103，NRF ne Abs NRFX
         · L103，New curvilinear step not converged. Error imposing constraints.
         · L114，ERROR IN INITNF. NUMBER OF VARIABLES (  0) INCORRECT (SHOULD BE BETWEEN 1 AND 50)
         · L202，Problem with the distance matrix.
         · L202，Atoms too close.
         · L123，Max corrector steps exceded
         · L123，GS2 Optimization Failure.
         · L123，GetHes: LRWHes > LHess!
         · L202，Atoms too close.
         · L301，The combination of multiplicity X and XXX electrons is impossible.
         · L301，End of file reading basis center.
         · L301，EOF while reading ECP pointer card.
         · L301，No solvent atoms in DisRep.
         · L301，Atomic number out of range for XXX basis set.
         · L301，R6DS8: Unable to choose the S8 parameter
         · L301，R6DRCv: No RCov radius available for IA=XX
         · L301，R6DC6: No C6 coefficient available for IA=XX
         · L502，Convergence failure
         · L502，Inaccurate quadrature in CalDSu.
         · L502，Inv3 failed in PCMMkU
         · L508，Convergence failure
         · L508，Inv3 failed in PCMMkU
         · L602，GetVDW:  no radius for atom XX atomic number XX.
         · L801，Excessive mixing of core and valence orbitals.
         · L801, Fatal Problem: The smallest alpha delta epsilon is XXXXXX
         · L913，*MAX. CYCLES*
         · L914，XXXXXX words are not enough for AIAXAO.
         · L1002，Inaccurate quadrature in CalDSu.
         · L1002，No func 3rd derivs with XXX.
         · L1002，NIJ > Max2 in MMCore.
         · L1002，XXXXXX words are not enough for AIAXAO.
         · L9999，Optimization stopped.

   · “计算机问题”报错

         · （CPU）Error: illegal instruction, illegal opcode

         · （文件）Internal input file was deleted!
         · （文件）open-new-file
         · （权限）Error Message # 2066. Can't create file Temp input file 'gxx.inp'

         · （硬盘）Erroneous write
         · （硬盘）g_write
         · （内存）needs more words of memory
         · （内存）Not enough memory to run at all
         · （内存）Out-of-memory error in routine XXX
         · （内存）galloc:  could not allocate memory.

   · 疑难杂症

         · NtrErr Called from XXXXXX.
         · Gaussian Bug 的“样式”及发现Bug后的可能做法
         · 已知的、可能为Bug的报错
                           L103，Bad arguments to LAPack or BLAS routine.
                           L401，Diagonalization in DiagDN via DSPEV failed.
                           L906，Internal consistency failure #1 in GetIJB.
                           L914，Unable to match L and R vectors in BiOrth
                           L1002，NIJ > Max2 in MMCore.
                           L1014，Tx not orthogonal to T

无报错的异常情况

         · IRC提前停止或折返
         · 优化后的结构有一个或几个小虚频
         · G09 Windows 在 L1.exe 左右无报错退出
         · GaussView观看结构有异常的成键、断键
         · GaussView 中显示的 IRC 在零点附近出现异常的低能点

不是报错、但常被误以为报错的语句

         · Warning!!: The largest alpha MO coefficient is...
         · Warning!!: The smallest alpha delta epsilon is...
         · Error on total polarization charges
         · This type of calculation cannot be archived.
         · No special actions if energy rises.
         · DIIS: error          · End of XXXXXX F.D. properties file XXX does not exist.

其他常见问题

         · GaussView 打开 chk 或 fchk 文件报错 Missing or bad data: Alpha Orbital Energies
         · “这种计算大约需多长时间？”
         · 如何通过柔性扫描产生较好的过渡态初猜
         · 如何在量化计算中设定原子/片段的电荷
         · 如何计算 Activation Strain (Distortion-Interaction) Energy

——————————————————————————————————————————————————————————————

遇到报错、提问前请做到：

充分试错：

总是有很多人的问题以 “XXXX这样行不行” 为模板提问程序问题，例如“程序安装在XXX路径可以吗？”、“G09 A.01能用XXX命令吗？”、“这样写会扫描10个点还是11个点？”、“不加上XXX会报错吗？”之类。这种问题显然一试便知，没有任何询问别人的必要。提问这种问题无非是一个字：懒！
做理论的人几分钟即可试错一次，都懒得试错；做实验的人每次试错简单的也要几个小时，岂不是没人答疑就彻底没法干了？试一试电脑又不会爆炸。

化简问题：

常有人提问时包含大量不可控因素，提问模式如“我做了一个多步计算，在第二步读取第一步的chk文件作为初猜，通过队列系统放在远程集群上运行，没能成功生成wfn文件，是什么原因？”这种问题没人能回答。谁知道是提交脚本问题、集群跪了、队列还没排到、压根就看错用户/目录了、Linux权限问题、软件安装问题、关键词不对、输入文件有错、第二步读chk出错、产生wfn文件的输入格式不对等等哪个原因？
作为科研工作者，遇到问题首先进行简化，排除所有不必要变量，得到一个能重复出错误的最简单案例（MRE，minimal reproducible example）应该是一项基本技能，也应是提问者的一项义务，很多问题在制作MRE的过程中也就知道怎么解决了。（https://en.wikipedia.org/wiki/Minimal_Working_Example）例如上面的问题可以逐步做简化，例如不用队列系统提交（直接用G09命令跑）、分成两步单独计算、设法尽量不读取chk文件等等，并比较每次化简后的区别。如果有之前可正常运行的结果，则逐一比较排除二者间的差别，找出问题所在。

应做的尝试包括：

若体系较大，尝试时间较长，首先用较小的简单体系测试，例如水分子等。
确认环境正确。可首先在自己熟悉的平台上尝试，如 Windows 上的 G09W。如自己未用过相关平台，首先确保可以成功运行某个保证正确的样例输入。
去除各种读取chk文件（guess=read, geom=allcheck，etc）、Restart、多步计算、除opt freq之外的复合运算（IRC Freq之类）等等自己未掌握清楚的复杂因素后再试，如果试成，再逐一加上无妨
删除多余关键词：阅读《常见的多余的和被滥用的Gaussian关键词》（http://sobereva.com/331），并按其中所说，去除所有不必要关键词后重新运行
删除或去尝试理解自己不懂什么意思的关键词再尝试

检索：

网络上各处“解决报错方法”常有错误，新手常从百度检索到上的回答，除部分回复者的答疑很靠谱外（如yjcmwgk，beefly等），绝大多数回答坑人尤甚，甚至可做有罪推定。
较靠谱的、保证正确的资源包括本论坛及Sob的博客（有Sobereva把关），可通过如下语法在Google中检索，效率很高

site:keinsci.com 关键词/报错
site:sobereva.com 关键词/报错

复制代码

另外可检索之前整理的群聊聊天记录，可使用这个小工具（http://bbs.keinsci.com/forum.php?mod=viewthread&tid=5207）进行检索，有用信息极多，且与网上乱搜不同，几乎可保证正确。
检索报错时注意删除各种特殊字符，如Gaussian输出中常带有完整路径，类似“Error termination via /home/gaussian/g09/l1002.exe”等等，偶尔会导致问题，搜索时检索“Error termination via l1002.exe”可能效果好些

提问时：

提问时请：

上传完整输入文件。

      常有人提问半天还不告知用了什么计算路径，让答疑着瞎猜半天后才知道是用了某愚蠢的关键词导致的。
      希望上传完整的输出文件、而不是截图关键词，一是某些错误不只是在关键字里，与坐标、坐标之后补充信息的格式有关，如果关键词用了genecp，只截取关键词更是连用了什么基组都不知道；二是，与截图不同，上传文件可以使答疑者立刻在自己的机器上做测试，远比语言描述和截图高效的多。
      如果要保密，请先尝试用H2O等简单的分子代替你的真实结构，看能否重现问题，若能，照此提问；若不能，将原子坐标段落保留第一个和最后一个，进行提问。

截图输出文件的至少最后一整屏、并最好附上输出文件。有些初学者截图十分吝啬，只截取它自己感觉相关的一点输出，然而真正的报错常在其他位置、自己选取的一点有时信息不够、有时毫无意义，还要答疑者反问，不必要的浪费大量时间。

(, 下载次数 Times of downloads: 459)

希望上传完整的输出文件有几点原因，一是有些问题需要观察整个输出文件，如优化超出步数限制时，是否震荡需要观察整个优化过程中的能量变化趋势。二是有些问题的报错离输出文件最后较远，如Gaussian报错Ntrerr，会在最后输出大量非人看的数字dump，人能看懂的部分在其之上；再如 l9999.exe 错误的真正原因输出于在末尾前100行左右的 optimization stop 附近，包括过渡态虚频数量不对、超出步数限制等等。
截图的同时上传文件可以帮助答疑者一眼看出简单的问题、复杂的问题也可以进一步找到输出，而不必再要。上传较大的输出文件、或多个文件注意压缩、打包。方便下载节省空间。

恰当、简明的描述问题。

            常遇到一些提问长篇大论、从课题背景说到导师初衷、从入门历程说到在哪儿买的电脑，根本不可能关系；或者一段话六七个问题纠缠在一起，比GRE文章还绕，根本找不到到底想问什么。
            有的提问者提问则极为简短、不知所云，“请问计算激发态log和fchk哪个更合适？”“我计算时给出了五个能量，和平时的五个能量该看哪个？”，谁知道提问者到底想干什么、用的什么程序、指的什么量？
            有的提问者则极具抽象与哲学化能力，“计算精度过高，对电子轨道描述过于详细，是否会引入很多本可以忽略的因素？而这些因素对结果有很大影响，是否也说明了计算一个体系的能量并不一定是精度越高越好？” 进行抽象没有任何问题，但常常这类问题背后其实是一个十分具体的实际计算问题，提问者自己经过了好几步（并不正确的）推理得出这个推论来提问，浪费双方的时间。

            故希望提问者描述问题时做到恰当、简明、直接，使提问：（1）self-contained（2）最简。
            提问时尽量使语言描述清晰、逻辑明确；有多个纠缠在一起的问题时恰当分段、标号列出。

告知自己做过哪些尝试、结果分别如何。有些初学者不知是为了突出自己的努力还是什么原因，常说如“我都试了”“试了很久，还是报错”“就是算不对”之类的话。这种表达毫无意义，也不会博得同情促使人回答问题。相反，如何具体给出“我用了XXX关键词，报错为XXX；又尝试了XXX，报错变为XXX”则含有大量的有效信息，答疑者也可以明确看出提问者确实做出了自己应有的努力，愿意提供帮助。再如说“换了XXX还是不行”就远不如“换了XXX报错相同/报错变为XXX”有效

说出运行环境（Gaussian的大、小版本及系统，如Gaussian 09 D.01 Linux）。旧版本的Gaussian有些已知或未知的Bug，告知此信息将有帮助

提问时请不要：

不说出具体报错信息。常有初学者贴出一个输入之后即问“这个输入为什么不对”“为什么掉下来”“为什么跑不了”“到这步就死了”，答疑的人又不是程序，不能人肉编译检查。除非是特别鲜见的typo之类，没人能回答。更有甚者直接说“这个关键词为什么不行”乃至“我优化了这个结构（附结构式）为什么不行”，输入输出都不给，连人肉编译器都没得当。让人挑错犹如去医院看病，问“为什么掉下来”，就如同只给医生说“我难受”就哑巴了、还不让做检查；除了少了俩胳膊断了条腿之类鲜见的毛病，谁也不知道怎么回事。显然应当给出“具体如何跑不了”的信息，详见上一段落。

仅说/仅截图 lXXX.exe 报错；仅按照 lXXX.exe 报错解决问题。问“L1234.exe报错该怎么解决”这种话没人能回答。Gaussian报错时会给出出错link名称，但哪怕是最容易确定的L502是SCF不收敛都有例外情况，其他的link报错几乎都需要此行之前的具体报错信息才能提供解决方案。同时也切勿仅按照Link数寻找解决方案，例如本文后面说“L801，Excessive mixing of core and valence orbitals 错误”应当使用Full选项解决，但从未说“L801错误”都应该用Full解决，一定要看具体报错信息。

用手机对着屏幕拍照代替截图。请参见 Sob 对此事的讨论：http://bbs.keinsci.com/forum.php?mod=viewthread&tid=5914。简而言之，拍照常出现对焦不准、图片体积大、无关区域多、有摩尔纹、亮度对比度极不利于观察、甚至转了90°等问题。体现出提问者在态度上不觉得麻烦别人有问题。各个环境下截图的方法一搜都有。如果实在出现没有联网、且不能将文件传回，也应尽量使图片质量较好，做到旋转恰当、对焦准确、尽量摆正、适当裁剪、适当调整亮度对比度等。这在拍照、以及手机图像应用的后处理中都很容易实现，举手之劳。

在多处同时提问。有人常在论坛和群里同时提问。答疑者在某处回答一遍后常在另一处又要重新看一遍，看到最后才发现早已在别处回答过了，令人心里很不舒服。如果在某处提问，至少也要24小时没人回答之后再去别处发同样的问题。

信口胡说。这是令我个人十分不解、十分反感的情况。提问者幻想自己做了并未做过的努力，尤其是在别人批评其未做检索、未按之前的回答尝试、未按提供的博文尝试时。常见模式含：

            “说的我都试过了”    （这种情况多数是分明只试了一两个就不耐烦等不及了，根本没做到“都试过”；甚至发现有的根本就没有试、不爱看长文，这么说只是为了让人直接给出一个关键词。像诸如SCF不收敛、几何优化不收敛之类的常见错误，Sob的博文里几乎罗列了所有Gaussian内部能做到的事情。如果真的遇到了罕见的难以想象的情况，都试了什么、错误的情况分别如何，请一个一个罗列出来以便讨论。）
            “我搜过了没有找到” （并不是说这句话本身有问题，而是有很多人的问题极端无脑，诸如“GGA全称是什么”之类，Google一下就是前三条结果之内的东西，却说搜不到）
            “我看了没有看懂” （这种回复着实令人反感。博文又不是文言文写的，也不像（你本应去看的）原始文献里全是公式和方程。
                     真正认真读过博文、认真思考的人绝不可能整篇文章“没有看懂”，他们问此问题都会十分具体、并一般经过自己思考有了自己的猜想，如:

                                    “文中XXX这句话是如何得出的？”
                                    “这一步是这样操作、还是那样操作”
                                    或者已经做了尝试，“为什么这样操作与原文不同”
                                    或对文中的句子进行转述（而不是原话粘贴），询问“是否应该这样理解”
                     而很少有说“都没有看懂”、“这段话没看明白”。
                     经验表明这些 “我看了博文没有看懂”，甚至 “我看了很多遍博文没有看懂”人常常是点开博文、看到有进度条就想说 tl;dr 的人；
                     常常是连几千个汉字都懒得顺序阅读，Ctrl+F 一下自己想要的关键词、发现没有、就来说“看了没看懂”的人。）

掺入主观判断，滥用“对照实验”。常有提问者用这种句式提问“如何找到带甲基化合物的过渡态”“加了苯环就算不出来怎么办”“PBE0为什么不能算过渡态”，令人莫名其妙，这都是加入了自己的主观判断、滥用“对照实验”。这类提问者也是不提供上一段所说的具体信息的常客，实际导致报错的原因常与其“对照试验”毫无关系，这样提问误导答疑者。

试想找到一个“不带甲基的化合物就能算完、带了甲基就报错了”的例子就说明是要解决甲基的优化问题吗？B3LYP能算过渡态，PBE0不知道什么原因报错了，就说明PBE0不能算过渡态吗？很多其他情况诸如初猜不好、恰好有原子呈直线、读取了错误的文件、电脑问题（硬盘、内存）、关键词不对、乃至仅仅数值巧合都可能有相应问题。

更有甚者算了两三个例子，就开始做归纳，“带甲基的‘都’算不出来，不带的‘都’可以”，对这种归纳深信不疑的请去Google一下“P-Hacking”的概念，然后考虑一下一个分子有多少种维度可以描述它，就明白了这种归纳的错误概率是多少了。遇到报错，应多从原理上推演、而少从结果上归纳；看具体报错解决，而不是自己瞎想。

抱怨。虽有上纲上线之嫌，但仍觉的值得一说。常见模式如：

提问者：XXX如何解决
答疑者：尝试一下XXX关键词。
（2分钟后）
提问者：【“我都加上了XXX，怎么还是不行啊！”】【“可是怎么还是不行？！！！”】【“可是算了半天又掉下来了！！！”】

复制代码

               试想如果上面的追问变成“我加上XXX关键词后，报错变为了XXX，又该如何解决？”是不是感觉提问者更加虚心、耐心、认真呢？
               这实际体现出科研环境越发浮躁、，多人巴不得能一键计算、一键导出、一键发文章；期待左右问题都有神人指点、保证成功，丝毫没有试错的意识。Sob也在《有感于如今许多超懒的科研工作者》（http://sobereva.com/342）一文中说过这些问题。
               人家程序积累了几十年、你写几个字节就能调用，不该感激涕零？
               有人表示有充足理由如“导师催着要”,“我不是这个方向的”等等。没有哪个导师不催着要，不是计算专业、搞计算的人也多了去了。着急的还大多要拿来发文章、并非算着玩玩：一知半解瞎算一气、计算机吐个数不管对不对就往上写。这往大了说，发表的垃圾结果早晚被后续实验/计算否认，这些凑（一个）实验数据的计算简直是败坏了计算化学的声誉，这些人是为“计算化学不靠谱”做出‘突出贡献’的一批人；往小了说，这也是发文章人的学术不端，对自己将要发表的结果、自己和导师的学术声誉毫无责任感。

伸手党。这就不用多说了，尤其是问诸如“能不能帮我写个输入文件”、“某某手册从哪儿下载”之类的问题。

“急急急”。经验表明这种人几乎100%是不急的，其中相当大一部分又都是极端基础、搜索一下就能知道的问题，如果真的“急急急”早就自己去搜了。这些人常妄想说了“急急急”大家就会重视、快速回复，其实正相反，它让人觉得提问者既轻浮又懒惰，根本不愿理。

提问后：

对热心回答的人表示适当的谢意。
认真理解回答者提供的信息，充分调研，举一反三。例如要求去看博文、看手册，请看过了再继续提问。不要没直接得到想要的答案就不管不顾继续按自己想要的错误的套路提问。
不论是否解决了问题，给予反馈，比如“之前问的Inaccurate quadrature in CalDSu的问题，用了scf=novaracc 确实解决了问题。”这为后人和旁观者的学习十分有利，也是一种对帮助的回馈。如果你的问题使用了自己查到、想出的方法解决了，则与人分享就增加了整个社区的知识总量，对自己也常在撰写解决方法的过程中学到更多东西，实为好事。

《向“好吧”作为句子开头宣战！》（http://bbs.keinsci.com/thread-334-1-1.html）
提问者的请求被“驳回”（如被要求去查Google，看手册）、甚至受到批评，在大多数情况下都是“讲道理的”。若不同意，多反思一下再予以对事、理性的反驳也无妨。

————————————————————————————————————————————————————————————————

Gaussian常见简单报错及解决方法：

（有哪些常见提问也可回复，我补充；如果下文有错误请务必批评指正）

————————————————————————————————

Windows Severe Error Message # 2070; Linux segmentation violation/segmentation fault

Linux

Error: segmentation violation
rax 0000000000000000, rbx 000000000060c640, rcx ffffffffffffffff
rdx 0000000000007f57, rsp 0000007fbffed948, rbp 0000007fbffed970
rsi 000000000000000b, rdi 0000000000007f57, r8 0000002a9558af40
r9 0000000000000000, r10 0000007fbffed801, r11 0000000000000206
r12 0000000000615890, r13 0000000000640398, r14 0000000000640398
r15 0000000000640398

复制代码

Windows

(, 下载次数 Times of downloads: 531)

成因：提供这种报错信息毫无意义，任何原因Gaussian报错退出都会有类似的输出。网上任何“解决segmentation violation”/“解决2070错误”的方案、解释都可无视并鄙视之。应当看输出文件中的具体错误解决。、

解决：用文本编辑器打开输出文件，如上图中为“a.out”文件，拖到最后看最终报错。Linux下也是如此，可用 nano, tail, cat, vi 等命令阅读文件。看到真正报错后，可按照相应信息（如下文中涉及的这些报错）进行解决。
在Windows下，因Gaussian 09W的界面只在Link更换或等待一定秒数之后才刷新输出，默认是15秒，读取输出还需要一定时间。故刚出错的时候不会在G09W天蓝色的那个图形界面中显示真正的报错信息，需等几十秒钟或/和点确定后才有真正的报错显示出来，有时一直不会更新。

————————————————————————————————

L502，Convergence failure

Convergence failure -- run terminated.
Error termination via Lnk1e in l502.exe

复制代码

成因：SCF不收敛

解决方法：https://github.com/liyuanhe211/S ... ges#scf-convergence

————————————————————————————————

L508，Convergence failure

Density matrix breaks symmetry, PCut= 6.91D-03
Density matrix has no symmetry -- integrals replicated.
Iteration 80 EE= -1377.03506721338 Delta-E= -0.001137227868 Grad=4.268D-02
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
Convergence failure.
Error termination via Lnk1e in C:\G09W\l508.exe at Mon Jun 26 19:36:04 2017.

复制代码

成因：用 QC 方法时SCF不收敛

解决方法：去掉SCF=qc、SCF=xqc、SCF=yqc 关键词进行计算。如果继而出现 L502, Convergence Failure问题，参照上文解决，只不过绕开涉及 scf=qc 的方案即可。

————————————————————————————————

L502，L1002；Inaccurate quadrature in CalDSu.

Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in l502.exe

复制代码

Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in l1002.exe

复制代码

首先应检查：

（1）基组/赝势是否合理，是否有例如赝势基组没写赝势之类的低级错误

（2）结构（或是轨迹）是否合理

如果均没有问题，则按下面的PPT解决：

(, 下载次数 Times of downloads: 457)

————————————————————————————————

L9999

Error termination via Lnk1e in l9999.exe

复制代码

此报错需查看输出文件中的额外信息，自结尾处向上检索“Optimization stop”可见如下几种情况

Optimization stopped.
-- Wrong number of Negative eigenvalues: Desired= 1 Actual= 4
-- Flag reset to prevent archiving.

复制代码

成因：Gaussian会在过渡态优化时自动检查虚频的数量，若不为1即终止。这一检查没有必要
解决：将NoEigenTest关键词加入到opt的选项中，如 opt=(TS, CalcFC, NoEigenTest)

Optimization stopped.
-- Number of steps exceeded, NStep= 100
-- Flag reset to prevent archiving.

复制代码

成因：几何优化未能在指定步数内收敛
解决：打开GaussView，File-Open，按照如下设置：

(, 下载次数 Times of downloads: 472)

然后点击Result-Optimization，可以看到优化时的能量变化曲线。观察优化曲线应纵向放大，方法为用鼠标反复选择纵轴方向仅包含最末几个点、但横轴方向包含所有点的“宽而矮”的矩形范围，直到能看清最末几个点的变化情况为止。根据图形判断是否震荡。如下图，需纵向放大至右边才能看清确实震荡了。（注意图中偶尔有突然的、无“周期性的”某1~2个点的spike是正常现象，不是震荡，下方小幅度但有“周期性”的才是震荡。）

(, 下载次数 Times of downloads: 432)

左侧看不出是否震荡。需纵向放大至右侧才能看出。

如果没有震荡：opt=maxcycle=当前步数的2~3倍。可以读取当前最末/最优的一个构象做为新的初猜。注意检查构象、轨迹是否合理，如果结构已经不合理了，则应适当调整初猜后再做（尤其是找过渡态的时候）。

如果震荡了：http://sobereva.com/164。

文中的众多解决方法中，我个人倾向于找到能量较低、4个收敛标准离收敛限较近，且结构合理的某点为新的初猜，并首先尝试在 opt 选项中加上(MaxStep=5, NoTrustUpdate, GDIIS) 选项；如果有计算频率的资源，可再加上opt=calcfc选项。特别仅对于对称性较高的结构，如T、C3 等群的复杂分子，若反复调整关键词后仍然震荡，可以尝试构建或破坏具有对称性的初猜结构，分别结合opt=Cartesian 优化，如 opt=(MaxStep=5,NoTrustUpdate,Cartesian), 有时有效果。

注意上一段仅是我个人首先尝试的建议，不要仅仅试了这一个不行就又把问题拿出来问。http://sobereva.com/164 有的是其他办法。

————————————————————————————————

L123，Max corrector steps exceded；GS2 Optimization Failure.

Delta-x Convergence NOT Met
Maximum number of corrector steps exceded.
Error termination via Lnk1e in l123.exe at Sun Sep 6 07:30:18 2015.

复制代码

GS2 Optimization Failure.
Error termination via Lnk1e in l123.exe

复制代码

成因：Gaussian的默认IRC算法HPC为了使曲线平滑需做重校正，但常有校正步不收敛的问题。换用GS2算法时，其限制性优化有时也不收敛。

解决：增加IRC=LQA选项，例如IRC=(calcfc,LQA)。

————————————————————————————————

L103，Error in internal coordinate system；Linear angle in Bend；Linear angle in Tors；FormBX had a problem.

Bend failed for angle 1 - 11 - 3
Tors failed for dihedral 9 - 1 - 11 - 3
Tors failed for dihedral 10 - 1 - 11 - 3
Tors failed for dihedral 12 - 1 - 11 - 3
Tors failed for dihedral 14 - 3 - 11 - 1
Tors failed for dihedral 17 - 3 - 11 - 1
FormBX had a problem.
Error termination via Lnk1e in l103.exe

复制代码

(Enter /home/gauuser/g09/l103.exe)
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Linear angle in Bend.
Error termination via Lnk1e in /home/gauuser/g09/l103.exe at Sun Feb 5 05:25:49 2017.

复制代码

NTrRot= -1 NTRed= 798 NAtoms= 66 NSkip= 606 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in l103.exe

复制代码

Berny optimization.
Using GEDIIS/GDIIS optimizer.
Linear angle in Tors.
Error termination via Lnk1e in /home/program/g09/l103.exe at Thu Jun 25 22:29:51 2015.

复制代码

成因：内坐标有其内在的限制，遇到优化过程中遇到有几个原子正好排成直线等情况可能出现此问题
解决： (0) 检查轨迹、当前结构、冻结（柔性扫描）变量设置是否合理
(1) opt=cartesian
   这一方法可从原理上彻底解决此问题，但 opt=cartesian 在多数情况下会增加优化至相应极小点所需的步数，如果体系不非常耗时可以直接用其至优化收敛。
   如果体系接近自己计算能力的极限，可以用 opt=cartesian 计算两三步之后终止计算，保存坐标后重新写一个输入文件，换回默认的opt方法。
   另外注意在使用modredundant冻结、扫描了坐标的时候，不能与opt=cartesian合用
(2) 有时直接保存最末结构后，重新opt也能解决此问题，Gaussian 实际上会自动为接近直线的原子添加一些 Linear Bend，但并不总是有效。
(3) 部分情况下，人为增加一 Linear bend 冗余内坐标，输入为在冗余内坐标段落加入如 L 1 2 3 -1 B 的表述，并启用Modredundant，其中 1 2 3 为 “Bend failed for angle” 的三个原子，如以上图报错来说，若原输入文件为：

#p B3LYP/genecp opt freq
Title
0 1
[原子坐标们]
C O H
TZVP
****
Co
SDD
****
Co
SDD

复制代码

则新的输入文件应为（注意共3处不同）:

#p B3LYP/genecp opt=modredundant freq
Title
0 1
[自出错的输出文件中提取的、出错前一步结构的原子坐标们]
L 1 11 3 -1 B
C O H
TZVP
****
Co
SDD
****
Co
SDD

复制代码

其中 “[自出错的输出文件中提取的、出错前一步结构的原子坐标们]” 可以用GaussView打开之前出错的输出文件，随便另存为一个gjf输入文件，用文本编辑器打开，把坐标部分拷出来。
此处用 genecp 并添加相应自定义基组、赝势的段落是为了说明输入文件中各段落的顺序。
如果你不知道我在说什么，可阅读：http://sobereva.com/60
完整的输入文件段落顺序可见于：http://sobereva.com/g09/m_input.htm

————————————————————————————————

L101，End of file in ZSymb.
L301，End of file reading basis center.
L301，EOF while reading ECP pointer card.

C -1.21995 2.13345 0.
End of file in ZSymb.
Error termination via Lnk1e in l101.exe

复制代码

End of file reading basis center.
Error termination via Lnk1e in /home/gauuser/g09/l301.exe at Wed Mar 1 20:50:57 2017.

复制代码

======================================================================================================
Pseudopotential Parameters
======================================================================================================
Center Atomic Valence Angular Power
Number Number Electrons Momentum of R Exponent Coefficient SO-Coeffient
======================================================================================================
EOF while reading ECP pointer card.
Error termination via Lnk1e in /home/gauuser/g09/l301.exe at Wed Mar 1 20:43:41 2017.

复制代码

成因：Gaussian需要读到换行表示的段落结尾才认为相应输入段落结束，若当前段落最后一个字符后，没有空行，Gaussian会读入一个EOF（End of File）特殊字符，不明白啥意思就跪了。这两个问题的提示不同是因为输入文件中分别以原子坐标段落（没有gen）、混合基组定义（gen）、和赝势定义（genecp）的段落为末尾，从而在读取不同信息时报错了。

解决：输入文件最后多打几个空行（应总是这么做，没坏处）。

————————————————————————————————

L101，Wanted a XXX as input. Found an XXX as input.

Wanted an integer as input.
Found a string as input.
H 0. 0. 0.
?
Error termination via Lnk1e in C:\G09W\l101.exe at Sun Mar 26 17:48:52 2017.
Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

成因：输入文件格式错误。高斯会检查特定段落的变量类型是否正确，诸如电荷、多重度应为整数（如 “1”），坐标应为浮点数（如 “1.0”），还有文本类型诸如Title段落等等，若不正确即报此错误。

解决：根据提示具体分析错误在哪儿。如故意制造上面报错使用的输入文件为：

#p B3LYP/6-31G(d)
Title
H 0. 0. 0.

复制代码

由报错提示“H 0. 0. 0. ”可知Gaussian在这一行想读到其他东西，而不是坐标，仔细思考就知道是缺少了电荷和多重度，都是整型变量，读到了“H”是文本型变量，故而报错，修改正确即可：

#p B3LYP/6-31G(d)
Title
0 2
H 0. 0. 0.

复制代码

————————————————————————————————

L101，Input Error Input Error Input Error Input Error Input Error Input Error

-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge = 0 Multiplicity = 1
Input Error Input Error Input Error Input Error Input Error Input Error
There are no atoms in this input structure !
Please fix the molecule specification section of your input and try again.
Input Error Input Error Input Error Input Error Input Error Input Error
Error termination via Lnk1e in C:\G09W\l101.exe at Wed Jan 16 20:26:00 2019.
Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

成因：Gaussian没有读到结构声明

解决：检查：（1）是不是真没写结构；（2）是不是忘了写geom=check, geom=allcheck；（3）是不是多了/少了空行，使高斯读入了空段落。

————————————————————————————————

L202，Problem with the distance matrix.

Small interatomic distances encountered: 6 1 7 2 8 3 9 4 10 5
Problem with the distance matrix.
Error termination via Lnk1e in C:\G09W\l202.exe at Wed Jan 16 20:33:11 2019.

复制代码

成因：有至少两个原子间距离太近，很有可能距离是0。例如上述的报错中的一串数字写明了是：6号原子和1号原子之间太近，7号和2号太近，8号和3号太近，etc
解决：找到相应原子调整分子结构，或者删去多余的原子。一种出现这种情况的原因是在GaussView中复制粘贴了一遍自己的结构而不自知，所以在GaussView里看着特别正常，但打开原子列表就可以发现问题了。

————————————————————————————————

L202，Atoms too close.

Small interatomic distances encountered:
2 1 5.00D-02
Atoms too close.
Error termination via Lnk1e in C:\G09W\l202.exe at Wed Jan 16 20:38:12 2019.

复制代码

成因：原子间距离太近，小于高斯认为正常结构应有的阈值（但不很接近于0）。例如上述例子中报错是2号原子和1号原子之间距离太近，2号和1号的距离是 0.05A。
解决：一般情况下，应该检查结构，确认其正常，拉远相应的原子；如果确实因为某种原因需要如此近的距离，可以用geom=nocrowd关闭这个检查。

————————————————————————————————

L1，Illegal IType or MSType generated by parse.

----------
#p sp freq
----------
Illegal IType or MSType generated by parse.
Error termination via Lnk1e in C:\G09W\l1.exe at Thu Dec 07 13:58:19 2017.
Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.
File lengths (MBytes): RWF= 1 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

成因：撰写输入文件时缺乏基本常识，将不搭界的、相互矛盾的计算任务写在同一个计算路径里。例如上面的输入告诉高斯要算单点、同时又说要算频率，高斯懵逼了。
解决：别犯傻。（不要自己乱组合，不要让 sp 和 freq，force 和 freq 等矛盾的任务设定在一条计算路径里出现）

————————————————————————————————

L1，QPErr --- A syntax error was detected in the input line.

------------------
#p M06-2X/6-31G(d)
------------------
QPErr --- A syntax error was detected in the input line.
#p M06-2X/6-31G(d)
'
Last state= "GCL"
TCursr= 3656 LCursr= 7
Error termination via Lnk1e in C:\G09W\l1.exe at Wed Mar 06 14:19:36 2019.
Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.

复制代码

成因：关键词或语法错误，错误的位置由下一行的 “ ’ ” 标记给出了
解决：按着手册写关键词，纠正关键词或语法错误。例如上面的例子中，程序标记出M06-2X关键字有问题，看手册得知高斯中正确的泛函写法应为M062X，修改后即可。

————————————————————————————————

L202，Atoms too close.

Small interatomic distances encountered:
2 1 4.87D-03
Atoms too close.
Error termination via Lnk1e in /gauuser003/Program/g09/l202.exe

复制代码

成因：输出已说明，Atoms too close.
解决：“2 1 4.87D-03”表示2号原子和1号原子之间的距离只有4.87D-03，这显然是不合理的。检查自己给的结构是否有问题。如果是在例如优化等过程中走了几步变成这样了，一定是级别（方法、基组）或其他设置有低级错误，诸如某个原子未给定基组之类，应仔细检查。
如果是在进行扫描键长的任务、扫描曲线某一端刻意要求某两个原子离得很近，可以使用Geom=NoCrowd关闭这个检查，但注意即使是在扫描里，这么近的地方一般也没什么意义；在其他任务里用 Geom=NoCrowd 就更掩耳盗铃了。

————————————————————————————————

L301，Atomic number out of range for XXX basis set.

Rotational constants (GHZ): 0.0817250 0.0474806 0.0408748
Standard basis: 6-31G(d) (6D, 7F)
Atomic number out of range for 6-31G basis set.
Error termination via Lnk1e in /data2/G09/g09/l301.exe at Fri Mar 24 20:58:27 2017.

复制代码

成因：使用内置基组进行计算时，分子中含有基组不支持（超出基组定义范围）的元素。如上面的案例中，计算使用了6-31G(d)基组，其只支持H~Kr；但结构中含有Pt原子，故报错为Pt超出基组定义范围。
解决：选择合适的基组（及赝势），按照规定格式输入，可参考：

谈谈量子化学中基组的选择
http://sobereva.com/336
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60爬虫并处理后的离线EMSL基组库
http://bbs.keinsci.com/forum.php ... hlight=%BB%F9%D7%E9
在线基组和赝势数据库一览
http://sobereva.com/309

————————————————————————————————

L301，The combination of multiplicity X and XXX electrons is impossible.

(Enter C:\G09W\l301.exe)
Standard basis: 6-31G(d) (6D, 7F)
Ernie: Thresh= 0.10000D-02 Tol= 0.10000D-05 Strict=F.
The combination of multiplicity 1 and 1 electrons is impossible.
Error termination via Lnk1e in C:\G09W\l301.exe at Sun Mar 26 18:00:32 2017.
Job cpu time: 0 days 0 hours 0 minutes 1.0 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

成因：完全没正确设置/理解自旋多重度或电荷。自旋多重度定义如下，即 Alpha电子数 - Beta电子数 + 1，一个电子要么是 alpha 电子，要么是 beta 电子。故可知偶数电子的体系，自旋多重度只能是奇数（1,3,5,7...），反之亦然。

(, 下载次数 Times of downloads: 412)

解决：正确设置自旋多重度、电荷、结构。
应检查诸如：自由基的自旋多重度习惯性的写了1；正负离子体系忘了改分子净电荷数；建模时结构中多了、少了氢原子之类的问题。
另外需注意，GaussView虽然看似会“自动判断”自旋多重度，但只是提供最低可能的自旋多重度而已（即偶数电子给1，奇数电子给2），应该用结构化学的知识自己判断体系可能的自旋多重度可能是多少（如闭壳层分子、单重态双自由基为1，一般有机自由基为2，金属配合物应该用配体场理论判断等等），如果实在不确定（如对某些过渡金属配合物），可以在不同的自旋多重度下优化分子，看哪种自旋多重度的能量最低。

————————————————————————————————

L301，R6DS8: Unable to choose the S8 parameter

IExCor= 408 DFT=T Ex=B Corr=PW91 ExCW=0 ScaHFX= 0.000000
ScaDFX= 1.000000 1.000000 1.000000 1.000000 ScalE2= 1.000000 1.000000
IRadAn= 0 IRanWt= -1 IRanGd= 0 ICorTp=0 IEmpDi=141
NAtoms= 1 NActive= 1 NUniq= 1 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=F Big=F
Integral buffers will be 262144 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
R6DS8: Unable to choose the S8 parameter, IExCor= 408 IXCFnc= 0 ScaHFX= 0.000000 IDFTD=4
Error termination via Lnk1e in C:\G09W\l301.exe at Wed Jul 12 20:20:24 2017.

复制代码

成因：使用 empiricaldispersion=GD3 或 empiricaldispersion=GD3BJ 时程序未存有相应泛函的D3色散校正参数。
解决：(1) 确认在D3或D3BJ中选择的正确，例如米尼苏达系列泛函没有D3-BJ参数，不能用empiricaldispersion=GD3BJ，只能用empiricaldispersion=GD3 (2) 换一个有色散校正的泛函 (3) 自定义色散校正参数（往往来源于较新的文献），可以参考《Gaussian中非内置的理论方法和泛函的用法》（http://sobereva.com/344）一文中2.2节的做法，以及手册中的相应部分 (4) 不用色散校正
————————————————————————————————

L301，R6DRCv: No RCov radius available for IA=XX
L301，R6DC6: No C6 coefficient available for IA=XX

Integral buffers will be 262144 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
R6DC6: No C6 coefficient available for IA= 96
Error termination via Lnk1e in C:\G09W\l301.exe at Sat Mar 17 13:56:16 2018.

复制代码

Integral buffers will be 262144 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
R6DRCv: No RCov radius available for IA= 96
Error termination via Lnk1e in C:\G09W\l301.exe at Sat Mar 17 13:57:31 2018.

复制代码

成因：使用了色散校正（empiricaldispersion=GD3，empiricaldispersion=GD3BJ，empiricaldispersion=GD2，以及使用了自带DFT-D的泛函比如wB97xD），且程序中没有相应元素的参数。上文的两个例子分别是在使用DFT-D2和DFT-D3时，分子中有锔原子（IA=XX 即为相应的元素序数）时的报错。例如对 DFT-D3 来说，原文只为 1-94 号元素（H-Pu）做了参数化，对 96 号元素使用 DFT-D3 当然会缺参数。

解决：不用色散校正。

————————————————————————————————

L114，ERROR IN INITNF. NUMBER OF VARIABLES ( 0) INCORRECT (SHOULD BE BETWEEN 1 AND 50)

NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-NEF-
NUMERICAL EIGENVECTOR FOLLOWING MINIMUM SEARCH
INITIALIZATION PASS
************************************************
** ERROR IN INITNF. NUMBER OF VARIABLES ( 0) **
** INCORRECT (SHOULD BE BETWEEN 1 AND 50) **
************************************************
Error termination via Lnk1e in C:\G09W\l114.exe at Sun Mar 26 16:57:35 2017.
Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

成因：使用无解析梯度的方法做优化而为明确声明变量，如CCSD(T)方法
解决：

更换方法，除非极特殊的目的，否则用CCSD(T)之类的方法优化总是没有必要的，如果体系很小、计算量很充裕，用CCSD优化已经足够了。
按下面 Sobereva 的 PPT 处理

(, 下载次数 Times of downloads: 358)

附：Sobereva制作的关于解析导数的PPT：

(, 下载次数 Times of downloads: 357)

(, 下载次数 Times of downloads: 328)

(, 下载次数 Times of downloads: 347)

(, 下载次数 Times of downloads: 381)

(, 下载次数 Times of downloads: 346)

————————————————————————————————

L801，Excessive mixing of core and valence orbitals.

Largest valence mixing into a core orbital is 5.06D-01
Largest core mixing into a valence orbital is 4.89D-01
Excessive mixing of frozen core and valence orbitals.
Error termination via Lnk1e in l801.exe

复制代码

成因：Gaussian在post-HF、双杂化、TD等运算时会做默认的冻芯，不将其考虑到激发行列式中，这有时会导致问题
解决：
（1）最简单的解决方法是在理论方法的关键词中加入(Full)选项，如做MP2计算出问题时对应用MP2(Full)代替、做CCSD(T)时用CCSD(T,full)代替，使程序不使用冻芯近似。
（2）做其他的冻芯设定、可见Gaussian手册对Frozen Core的说明（http://sobereva.com/g09/k_fc.htm）
（3）IOp(8/11=1)可防止在检查得到冻芯不合理时即终止任务，而仅将这个警告打印出来、接着进行后续运算。但切记自己检查结果的合理性。

————————————————————————————————

L801, Fatal Problem: The smallest alpha delta epsilon is XXXXXX

(Enter /home/bjwang/g09/l801.exe)
Range of M.O.s used for correlation: 1 264
NBasis= 264 NAE= 114 NBE= 114 NFC= 0 NFV= 0
NROrb= 264 NOA= 114 NOB= 114 NVA= 150 NVB= 150
**** Fatal Problem: The smallest alpha delta epsilon is -0.78537370D-01
Error termination via Lnk1e in /home/bjwang/g09/l801.exe at Mon Jan 31 18:17:13 2011.

复制代码

成因：在TDDFT、post-HF 中做积分变换时涉及分子轨道的能量差为分母的项，若分子轨道Gap太小（“alpha delta epsilon” 即两个α轨道的△ε），分母接近0，得到一个大数会导致数值不稳定性问题。具体参看本文后面关于“Warning!!: The smallest alpha delta epsilon is...”的章节。

解决：由于未能找到可重复的算例，以下解决方案仅为推测、未经完全证实，仅供尝试
（1）IOp(8/11=1)可防止在出现此警告时即终止任务，仅将这个警告打印出来、接着进行后续运算。但切记自己检查结果的合理性。
（2）换泛函，更换HF成分较高的泛函可能能解决这个问题（增大Gap）

————————————————————————————————

L502 & L508，Inv3 failed in PCMMkU

(Enter /home/g16/l502.exe)
Integral symmetry usage will be decided dynamically.
Closed shell SCF:
Using DIIS extrapolation, IDIIS= 1040.
NGot= 18983288832 LenX= 18982597916 LenY= 18982263254
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on energy=1.00D-06.
No special actions if energy rises.
Fock matrices will be formed incrementally for 20 cycles.
Cycle 1 Pass 1 IDiag 1:
FoFJK: IHMeth= 1 ICntrl= 0 DoSepK=F KAlg= 0 I1Cent= 0 FoldK=F
IRaf= 990000000 NMat= 1 IRICut= 1 DoRegI=T DoRafI=F ISym2E= 0 IDoP0=0 IntGTp=1.
FoFCou: FMM=T IPFlag= 0 FMFlag= 100000 FMFlg1= 2001
NFxFlg= 0 DoJE=F BraDBF=F KetDBF=F FulRan=T
wScrn= 0.000000 ICntrl= 0 IOpCl= 0 I1Cent= 0 NGrid= 0
NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0
Symmetry not used in FoFCou.
FMM levels: 10 Number of levels for PrismC: 9
Inv3: Mode=1 IEnd= 55186563.
Iteration 1 A*A^-1 deviation from unit magnitude is 2.18D-11 for 1454.
Iteration 1 A*A^-1 deviation from orthogonality is 1.36D-11 for 4245 1454.
Iteration 1 A^-1*A deviation from unit magnitude is 2.55D-11 for 1454.
Iteration 1 A^-1*A deviation from orthogonality is 2.25D-10 for 4227 3067.
Iteration 2 A*A^-1 deviation from unit magnitude is 7.03D-07 for 4227.
Iteration 2 A*A^-1 deviation from orthogonality is 8.89D-07 for 2879 1454.
Iteration 2 A^-1*A deviation from unit magnitude is 1.00D-06 for 1454.
Iteration 2 A^-1*A deviation from orthogonality is 9.53D-07 for 4227 1454.
...
...
Iteration 10 A*A^-1 deviation from unit magnitude is 7.97D-07 for 2943.
Iteration 10 A*A^-1 deviation from orthogonality is 1.08D-06 for 2943 1454.
Iteration 10 A^-1*A deviation from unit magnitude is 3.30D-07 for 4227.
Iteration 10 A^-1*A deviation from orthogonality is 6.08D-07 for 4227 2061.
Inv3 failed in PCMMkU.
Error termination via Lnk1e in /home/g16/l502.exe at Mon Sep 4 11:33:20 2017.
Job cpu time: 0 days 0 hours 18 minutes 30.5 seconds.
Elapsed time: 0 days 0 hours 0 minutes 42.9 seconds.

复制代码

成因及解决方法见：http://bbs.keinsci.com/thread-6682-1-1.html

————————————————————————————————

L1002，NIJ > Max2 in MMCore.

Estimated number of processors is: 10
Inverted reduced A of dimension 2103 with in-core refinement.
NIJ > Max2 in MMCore.
Error termination via Lnk1e in l1002.exe at Thu Jun 2 11:12:43 2016.

复制代码

成因：已知、且已修复的Gaussian程序Bug，自G09 C.01版本之后已修复
（G09 C.01 Release Note：“A memory allocation bug for very large systems, which could cause a failure with the message “NIJ > Max2 in MMCore,” was fixed”）

解决：更换G09 C.01或之后版本。

————————————————————————————————

L913，*MAX. CYCLES*

CCSD、CCSD(T)：

Iteration Nr. 50
**********************
DD1Dir will call FoFMem 1 times, MxPair= 2
NAB= 1 NAA= 0 NBB= 0.
Norm of the A-vectors is 4.0812478D-03 conv= 1.00D-05.
RLE energy= -0.0380878830
DE(Corr)= -0.37790236E-01 E(CORR)= -1.1531999202 Delta=-1.55D-03
NORM(A)= 0.10071894D+01
*************
*MAX. CYCLES*
*************
Largest amplitude= 4.58D-02
Error termination via Lnk1e in C:\G09W\l913.exe at Sun Mar 05 10:12:08 2017.
Job cpu time: 0 days 0 hours 0 minutes 1.0 seconds.
File lengths (MBytes): RWF= 22 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

CISD

Iteration Nr. 50
**********************
DD1Dir will call FoFMem 1 times, MxPair= 2
NAB= 1 NAA= 0 NBB= 0.
Norm of the A-vectors is 4.3269814D+01 conv= 1.00D-05.
RLE energy= -9.6293060152
DE(CI)= 5.1143514 E(CI)= 4.3978787385
NORM(A)= 0.36740722D+02
SIZE-CONSISTENCY CORRECTION:
S.C.C.= 0.00000000D+00 E(CI,SIZE)= -0.71597715082D+00
*************
*MAX. CYCLES*
*************
***************************************************************
Dominant configurations:
***********************
Spin Case I J A B Value
AA 1 3 -0.354624D+01
......
ABAB 1 1 23 23 0.147097D+00
Largest amplitude= 1.84D+01
Error termination via Lnk1e in C:\G09W\l913.exe at Sat Jun 11 03:43:10 2016.
Job cpu time: 0 days 0 hours 0 minutes 2.0 seconds.
File lengths (MBytes): RWF= 23 Int= 0 D2E= 0 Chk= 1 Scr= 1

复制代码

成因：高斯用某种迭代方法求解 CISD 和 CCSD（CCSD(T)），而该迭代过程也会不收敛，在默认循环（50步）或设定的循环内未收敛即会报此错误。（注意CISD报出 Max Cycles 错误之后还会吐一堆组态，真正报错在它上面）

解决：
如下">>>"标出的一行所示，计算过程中会输出每步迭代后的CCSD / CISD相关能（DE(corr)），CCSD还会输出本轮迭代的能量变化（Delta，即本轮的 DE(corr) 减去上一轮的 DE(corr)），需观察前者是否有稳定的收敛至某一值、后者是否逐渐趋近于0，以判断是否有收敛趋势。

**********************
DD1Dir will call FoFMem 1 times, MxPair= 380
NAB= 190 NAA= 0 NBB= 0.
Norm of the A-vectors is 1.2424967D-05 conv= 1.00D-05.
RLE energy= -1.1041652032
>>> DE(Corr)= -1.1041652 E(CORR)= -236.55647322 Delta=-7.56D-08
NORM(A)= 0.11814199D+01
Iteration Nr. 12
**********************

复制代码

如果上述能量有收敛趋势，增加CCSD、CCSD(T)、CISD迭代循环的最大值，默认为50，应大于此值，语法如 CCSD(T) 改为 CCSD(T, maxcyc=100)
若能量在震荡，应首先检查结构、参考态是否合理，如超过稳定点的单重态双自由基是否用了对称性破缺的初猜；另外可尝试微调结构、换基组
如果最终震荡的幅度很小，例如能量的默认收敛限为1E-7，最终Delta一直在 nE-7的水平，到不了E-8，可以用诸如 CCSD(T,conver=6) 将收敛阈值提高至 1E-6。（注意如果是在热力学组合方法中遇到这个问题，就不应该写CCSD(T,conver=6)，而应该只写CCSD(conver=6)。

————————————————————————————————

L103，New curvilinear step not converged. Error imposing constraints.

Iteration 96 RMS(Cart)= 0.00000206 RMS(Int)= 0.00542712
Iteration 97 RMS(Cart)= 0.00000193 RMS(Int)= 0.00542766
Iteration 98 RMS(Cart)= 0.00000180 RMS(Int)= 0.00542817
Iteration 99 RMS(Cart)= 0.00000169 RMS(Int)= 0.00542865
Iteration100 RMS(Cart)= 0.00000158 RMS(Int)= 0.00542909
New curvilinear step not converged.
Error imposing constraints
Error termination via Lnk1e in /home/gauuser/g09/l103.exe at Mon Jul 17 10:21:34 2017.

复制代码

成因：进行需要做限制性优化的任务（如opt=modredundant 的 F 和 S、QST2等），Optimizer 不知道在当前限制条件下该把结构初猜摆成什么样子。

解决：

如果用了QST2：

换用TS（Berny）方法找过渡态，QST2不是什么好方法

如果在柔性扫描中遇到问题：

（1）如果是在从上一个扫描点（如键长为1.7A）、至下一个扫描点的过程中产生的问题（如键长为1.8A），则可以手动将初猜结构中的相应键长调整到下一点的键长（1.8A），从这一点开始后续扫描。最后把结果拼接起来就可以了
（2）扫描过程中频繁出现此问题有可能是“步子扯得太大”，可以把步长调小
（3）如果有很快速、合适的方法进行预扫描（如PM6-D3）可以先在相应级别下做扫描，产生每一点的结构，然后用多个限制性优化（F）代替扫描（S），此时由于并行化是100%的，只要预优化合理、效率还比柔性扫描高。

如果在限制性优化中（或柔性扫描的某一点中间）遇到问题：

（1）稍微调整结构重新来
（2）偶尔将结构重新保存为输入文件、直接接着运行问题也可以解决。原因不明。可以尝试

————————————————————————————————

L301，No solvent atoms in DisRep.

Solvent : n,n-DiMethylFormamide, Eps= 37.219000 Eps(inf)= 2.046330
RSolv= 0.000000 Ang.
------------------------------------------------------------------------------
Warning! Inconsistent VMol and RSolv for this solvent, using VMol= 0.00 Ang**3.
No solvent atoms in DisRep.
Error termination via Lnk1e in /home/igors/g09/l301.exe at Fri May 19 02:16:27 2017

复制代码

成因：

示例出错的输入文件如下（下面的乙酸溶剂没有缺乏相应非极性参数）：

# opt() freq m062x/6-311G(2df,p) SCRF(CPCM, solvent=AceticAcid, read)
Title Card Required
0 1
O 0.00000000 0.00000000 0.11813800
H 0.00000000 0.75681000 -0.47255200
H 0.00000000 -0.75681000 -0.47255200
dis
rep
cav

复制代码

即使用PCM溶剂模型时使用 “Dis Rep Cav” 考虑了非极性部分，而所使用的相应溶剂在高斯中并没有定义非极性部分的参数。高斯列表支持的溶剂当中（G09 或 G16 有约180种）只有 G03 年代就支持的约20种有相应参数支持以这样的方式考虑非极性部分（可参考 G03 的说明书），其他均会报此错误。

解决：高斯或许可能有方法自定义 Dis Rep Cav 的参数，但未在说明书中给出。可以用SMD模型考虑非极性部分（对自定义溶剂可见：http://bbs.keinsci.com/forum.php?mod=viewthread&tid=6683）；或是不考虑非极性部分。

————————————————————————————————

L103，NRF ne Abs NRFX

Cartesian Forces: Max 0.101934913 RMS 0. 010471640
Leave Link 716 at Sun Jul 3 00:08:55 2016, MaxMem= 1073741824 cpu:0.7
(Enter /home/export/base/fanyu/apps/g09/1103.exe)
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Internal Forces: Max 0.077268082 RMS 0.018267038
Search for a local minimum.
Step number 1 out of a maximum of 202
All quantities printed in internal units (Hartrees-Bohrs-Radians)
NRF ne Abs NRFX
Error termination via Lnkle in /home/export/base/fanyu/apps/g09/1103.exe at Sun Jul 3 00:08:57 2016.

复制代码

成因：在ONIOM计算中冻结不合理。此例中其冻结了所有 Low layer 而放开了所有 High layer。
对 Link atom 进行处理时，由于Link atom的种类在两层中不同（如 Si-O 键打断，用 Si-H 封端），两层中的键长应当不同；所以 Link atom 和它在 High layer 里连接的原子要么都冻结、要么都放开。

解决：
1. link atom和其连接的原子用相同方式处理
2. 对于中等大小的体系（<300个原子），可以使用Opt=NoMicro来取消MicroIteration

详细讨论见：http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3671

————————————————————————————————

L123，GetHes: LRWHes > LHess!

******** Start new reaction path calculation ********
RCFC Option Requested - Data Read From Chk File:
"/home/gauuser/Gaussian/PavB_Rad_Cyc/TS3/To_IM3__To_IM3_3[Complete_PBE1PBE]_Step
6[IRC_Forward].chk"
Energy From Chk = -2089.1954828
Current Structure is TS -> form Hessian eigenvectors.
Diagonalizing Hessian.
WARNING: NO IMAGINARY FREQUENCIES AT TS!
Supplied step size of 0.1000 bohr.
Integration on MW PES will use step size of 0.2173 sqrt(amu)*bohr.
Point Number: 0 Path Number: 1
...
Leave Link 716 at Wed Jun 28 18:41:00 2017, MaxMem= 18666618880 cpu: 0.5
(Enter /home/gauuser/g09/l123.exe)
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
******** Start new reaction path calculation ********
LRWHes = 0 LHess = ***
GetHes: LRWHes > LHess!
Error termination via Lnk1e in /home/gauuser/g09/l123.exe at Wed Jun 28 18:41:00 2017.
Job cpu time: 0 days 0 hours 56 minutes 39.2 seconds.

复制代码

成因：此错误较为罕见，未能获得更多的错误样本。仅有的一个例子中，经分析是由于 1.读取了前一步 opt=TS freq 的Hessian，2. 该Hessian无虚频（注意前面有“WARNING: NO IMAGINARY FREQUENCIES AT TS!”提示）。
解决：检查过渡态的频率是否正常、只有一个虚频。若虚频数不对、首先解决过渡态优化的问题再说；如果没有该问题可上传文件讨论。

————————————————————————————————

L1002，No func 3rd derivs with XXX.

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

复制代码

成因：此任务要求Gaussian用LC-wPBE泛函计算超极化率，polar任务计算超极化率需要三阶解析导数且默认认为DFT方法都有三阶解析导数，但可见于Gaussian手册对DFT方法的说明，少数泛函没有:

AVAILABILITY
Energies, analytic gradients, and analytic frequencies; ADMP calculations.
Third order properties such as hyperpolarizabilities and Raman intensities are not available for functionals for which third derivatives are not implemented: the exchange functionals Gill96, P (Perdew86), BRx, PKZB, TPSS, wPBEh and PBEh; the correlation functionals PKZB and TPSS; the hybrid functionals OHSE1PBE and OHSE2PBE.

复制代码

解决：使用具有三阶导数的方法；或令程序用数值方法计算三阶解析导数
具体至使用LC-wPBE计算极化率的问题：
（1）用 polar=Numerical 计算数值方法的极化率α
（2）可以使用默认的Polar关键词，虽然其在计算超极化率β时会因为LC-wPBE没有三阶导数而报错，但报错前极化率α实际上已经输出来了（如下面段落“>>>”所示的行）。

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.

复制代码

————————————————————————————————

L602，GetVDW: no radius for atom XX atomic number XX.

(Enter /home/software/g09D01/g09/l602.exe)
FitSet: NAtFit= 76 NAtPot= 76 NAtFrz= 0 MDM= 80 TotChg= 0.00000
Merz-Kollman atomic radii used.
GetVDW: no radius for atom 2 atomic number 29.
Error termination via Lnk1e in /home/software/g09D01/g09/l602.exe at Mon Dec 26 13:57:24 2016.
Job cpu time: 0 days 4 hours 6 minutes 15.9 seconds.
File lengths (MBytes): RWF= 856 Int= 0 D2E= 0 Chk= 31 Scr= 1

复制代码

原因：使用pop=CHELPG 拟合静电势时没有内置相应元素的半径（上面例子中是第29号元素Cu）。
解决：pop里用readradii，输入文件末尾写上元素名和指定的半径（一般用范德华半径，可以查得），例如此例对铜来说应写如下的关键词：

#p (其他关键词) pop=(CHELPG,ReadRadii)
title
坐标们
C H N P Cl 0
6-31G*
****
Cu 0
lanl2dz
****
Cu 0
lanl2dz
Cu 1.4

复制代码

这里用赝势是为了体现各段落的顺序关系

————————————————————————————————

Internal input file was deleted!

Leave Link 703 at Wed Mar 1 20:23:53 2017, MaxMem= 18983419904 cpu: 101.8
Internal input file was deleted!
Error termination via Lnk1e at Wed Mar 1 20:23:53 2017.

复制代码

该报错总发生于某一Link刚刚结束、新Link尚未开始之时，但具体发生在哪一Link并不固定。

成因：Gaussian似乎需要一个完整的（没有“@”）、当前步的输入文件来完成Link之间的交接，这个临时文件存储于 $GAUSS_SCRDIR 环境变量下，文件名和路径都会在输出文件中写明，如下面的任务中，该临时文件就是“/home/gauuser/g09/scratch/Gau-16535.inp”。这个报错的原因是这个文件somehow在运行之中被删掉了。

Initial command:
/home/gauuser/g09/l1.exe "/home/gauuser/g09/scratch/Gau-16535.inp" -scrdir="/home/gauuser/g09/scratch/"
Entering Link 1 = /home/gauuser/g09/l1.exe PID= 16536.

复制代码

解决：重算（或视任务类型，从当前步骤Restart），并尝试指认这个临时文件被删除的原因并解决它。遇到过的情况包括：使用了一些奇怪的脚本、如在某任务 A 结束后清除临时文件夹，但却误将任务 B 的临时文件也清除了；队列控制、文件系统等问题；当然还包括自己、或使用同一机器的其他人手贱→_→。

————————————————————————————————

L1，open-new-file

Gaussian 09: IA32W-G09RevD.01 24-Apr-2013
05-Nov-2017
******************************************
%chk=C:\gtest\C6H5NO2.chk
fname=C:\gtest\C6H5NO2.chk fd = -1
open-new-file

复制代码

原因：高斯无法创建、访问 %chk（或其他Link0指令）指向的文件（本例中是 C:\gtest\C6H5NO2.chk）
解决：

确认相应目录确实存在（此例中是 C:\gtest 目录）。
确认高斯对其有读写权限（Windows下赋予高斯（G09W.exe）管理员权限，Linux 下对调用高斯的用户调整相应文件夹所有权与权限）。
确认相应文件没有正被占用。

————————————————————————————————

Windows G09W，Error Message # 2066. Can't create file Temp input file 'gxx.inp'

(, 下载次数 Times of downloads: 368)

成因：Gaussian无法创建临时文件

解决：确认Gaussian 09W File-Preference 中的 Scratch Path 目录存在；确认Gaussian对该目录有读写权限（可以尝试修改至系统盘之外的其他磁盘；或给予 Gaussian09w.exe 管理员权限）；确认磁盘未满。

(, 下载次数 Times of downloads: 396)

————————————————————————————————

L1，ntrex1

%mem=130000MB
%chk=D:\Gaussian\ORCaussian_test\Parallel\ORCA_derv\calcfc_opt_TS2\TS.chk
ntrex1

复制代码

成因、解决：Gaussian无法访问指定的路径，例如在Linux上写了Windows下的路径，改正确了即可。
————————————————————————————————

Erroneous write; g_write

Erroneous write. Write -1 instead of 4096.
fd = 4
orig len = 4096 left = 4096
g_write

复制代码

成因：硬盘不够
解决：

(, 下载次数 Times of downloads: 380)

————————————————————————————————

L914 or L1002，XXXXXX words are not enough for AIAXAO.
Out-of-memory；needs more words of memory；Not enough memory；need an additional words of memory .

Out-of-memory error in routine FoFDir-SEAll (IEnd= 1 MxCore= -1974194264)
Use %mem=1915MW to provide the minimum amount of memory required to complete this step.
Error termination via Lnk1e in l914.exe at Tue Mar 08 15:32:05 2016.

复制代码

Iteration 1 Dimension 60 NMult 60
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX needs 2533811 more words of memory.
Error termination via Lnk1e in l914.exe

复制代码

Not enough memory to run at all: LenERI= 2254414
Increase memory by 1247333 words.
Error termination via Lnk1e in d:\gaussian03\l804.exe

复制代码

GetIJB would need an additional 555279372 words of memory to use all 12 processors.
JobTyp=2 Pass 1: I= 41 to 74 NPSUse= 3 ParTrn=F ParDer=T DoDerP=T.

复制代码

Generate precomputed XC quadrature information.
Solving linear equations simultaneously, MaxMat= 72.
AlAXAO: NMat= 690 NPMax= 1 NPMax1= 0 MaxMat= 72.
33484188 words are not enough for AlAXAO.
Error termination via Lnk1e in /cluster/gaussian-09.C.01-amd64-sse4a//g09/l1002.exe at Tue Jan 6 15:17:39 2015.

复制代码

Convergence on wavefunction: 0.001000000000000
Iteration 1 Dimension 128 NMult 0 NNew 128
AlAXAO: NMat= 128 NPMax= 1 NPMax1= 0 MaxMat= 0 Max3X=F.
33548910 words are not enough for AlAXAO.
Error termination via Lnk1e in /public/home/qsh/g09/l914.exe at Tue Jul 11 13:23:34 2017.

复制代码

成因：提示都非常直接，内存不够
解决：用%mem设置增加内存即可。
对某些情况，如“need an additional words of memory to use all XX processors.”，若没有足够的内存提供，可以减少运行核数（%nprocshared）。若仍然不够，加内存条 or 降低计算机别

————————————————————————————————

galloc: could not allocate memory.

Leave Link 1 at Fri Mar 8 01:49:58 2013, MaxMem=17716740096 cpu: 0.8
galloc: could not allocate memory.

复制代码

成因：输入文件中准许Gaussian程序使用的内存超过了机器当前剩余的内存量，故Gaussian无法想系统申请到想要数量的内存而报错。
解决：将输入文件中 %mem 命令后的内存量设置为当前系统剩余的内存量，并留下足够的余量（防止运行途中剩余内存减少而又出现这个问题）。在Windows下查看剩余内存可用任务管理器；在 Linux 下用 top 命令查看；若同时运行多个任务，也需额外注意内存分配问题。

————————————————————————————————

NtrErr Called from FileIO. NtrErr Called from XXXXXX.

Structure from the checkpoint file: "sideroidide_CS_1_opt.chk"
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -584 Len= 36 IPos= 0 Q= 47012767336112
dumping /fiocom/, unit = 1 NFiles = 33 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 1206272 FType=2 FMxFil=10000
Number 0 501 502 503 507 511
Base 216576 43520 72192 180736 181760 196608
End 1206272 44520 76297 180821 181925 197835
End1 1206272 44544 76800 181248 182272 198144
Wr Pntr 216576 43520 72192 180736 181760 196608
Rd Pntr 216576 44520 72192 180736 181760 196608
Length 989696 1000 4105 85 165 1227
Number 551 552 561 562 575 579
Base 213504 212480 214016 205824 183808 212992
End 213542 212501 214017 212016 196499 213016
End1 214016 212992 214528 212480 196608 213504
Wr Pntr 213504 212480 214016 205824 183808 212992
Rd Pntr 213504 212480 214016 205824 183808 212992
Length 38 21 1 6192 12691 24
Number 598 603 665 672 674 698
Base 76800 216064 181248 183296 182784 215040
End 76802 216065 181684 183506 182924 215112
End1 77312 216576 181760 183808 183296 215552
Wr Pntr 76800 216064 181248 183296 182784 215040
Rd Pntr 76800 216064 181248 183296 182784 215040
Length 2 1 436 210 140 72
Number 700 701 730 761 801 989
Base 215552 198144 182272 214528 180224 44544
End 215842 205553 182514 214529 180230 64544
End1 216064 205824 182784 215040 180736 65024
Wr Pntr 215842 198144 182272 214528 180224 44544
Rd Pntr 215552 198144 182502 214528 180224 44544
Length 290 7409 242 1 6 20000
Number 991 992 993 994 995 996
Base 65536 65024 43008 40448 42496 41472
End 72098 65033 43208 40478 42516 41672
End1 72192 65536 43520 40960 43008 41984
Wr Pntr 65536 65024 43008 40448 42496 41472
Rd Pntr 72098 65033 43208 40478 42516 41672
Length 6562 9 200 30 20 200
Number 997 998 999
Base 41984 40960 77312
End 42274 41160 179816
End1 42496 41472 180224
Wr Pntr 41984 40960 77312
Rd Pntr 42274 41160 79816
Length 290 200 102504
dumping /fiocom/, unit = 2 NFiles = 1 SizExt = 0 WInBlk = 512
defal = F LstWrd = 65536 FType=2 FMxFil=10000
Number 0
Base 40448
End 65536
End1 65536
Wr Pntr 40448
Rd Pntr 40448
Length 25088
dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 65536 FType=2 FMxFil=10000
Number 0
Base 40448
End 65536
End1 65536
Wr Pntr 40448
Rd Pntr 40448
Length 25088
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -584 Len= 36 IPos= 0 Q= 47012767336112
Error termination in NtrErr:
NtrErr Called from FileIO.

复制代码

%chk=temp.chk
Bad file opened by FileIO: Unit=2 I= 2 FPrev=40960 FCur= 0.
FileIO: IOper= 9 IFilNo(1)= 2 Len= 0 IPos= 0 Q= 135419884
dumping /fiocom/, unit = 1 NFiles = 1 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 65536 FType=2 FMxFil=10000
Number 0
Base 20480
End 65536
End1 65536
Wr Pntr 20480
Rd Pntr 20480
Length 45056
Error termination in NtrErr:
NtrErr Called from FileIO.

复制代码

NtrErr Called from XXXXXX的报错成因非常复杂，且真实报错离末尾较远，初学者常只贴 Error termination in NtrErr: NtrErr Called from XXXXXX. 乃至最后一屏幕等，实际上在其上的一堆数字之上，例如最开始的案例的实际报错是，后面一堆乱码一样的文字对用户没啥用：

Structure from the checkpoint file: "sideroidide_CS_1_opt.chk"
FileIO operation on non-existent file.

复制代码

此时就很有可能是chk文件出了什么问题。应首先检查chk文件是否存在、可读，例如可以检查：
（1）路径是否正确、完整（如路径存在中文、存在空格、存在括号、没写扩展名等低级错误）
（2）使用相对路径时，“相对”的不对，如两次运行程序时CWD（PWD）不同，导致不能找到同一个文件，可以将完整路径（Windows从盘符开始，Linux从/开始）
（3）文件权限问题。在Linux上检查路径存在、相应文件（文件夹）有正确的读写权限。在Windows上，不要将临时文件存储在需要特殊权限的目录下，或者给Gaussian程序管理员权限后再运行。
（4）文件是否存在，例如我见过某些脚本在计算完成后会自动将chk文件删除、或将chk文件转化为fchk文件，此时当然会出错。
（5）如果实在找不到本质、应当先去除所有需要用到chk文件做输入的因素后重试，例如放弃使用geom=allcheck, guess=read之类的关键词从头开始算试试，可解决一大部分报错。这在本文一开始的部分也提到过。

注意类似因chk文件问题产生的报错有时并不直接，如下面这个案例中：

******************************************
Gaussian 09: ES64L-G09RevE.01 30-Nov-2015
28-Sep-2016
******************************************
%nprocshared=14
Will use up to 14 processors via shared memory.
%mem=43522MB
----------------------------------------------------------------------
#p m062x/genecp opt=(nofreeze,noeigentest,readfc,gdiis,maxstep=10,notr
ust,ts,restart) int=ultrafine freq empiricaldispersion=gd3 5d 7f scrf=
(smd,solvent=acetonitrile)
----------------------------------------------------------------------

复制代码

......

Berny optimization.
Restoring state from the checkpoint file "/home/gauuser/Gaussian/PhS_Rad_Cyc/4_5
_TS/4_5_TS[Complete_M062X]_Step2[TS_opt].chk".
ONIOM data not found on unit 2.
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -997 Len= 20 IPos= 0 Q= 140735655680592
dumping /fiocom/, unit = 1 NFiles = 12 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 596480 FType=2 FMxFil=10000
......
dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 65536 FType=2 FMxFil=10000
Number 0
Base 40448
End 65536
End1 65536
Wr Pntr 40448
Rd Pntr 40448
Length 25088
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -997 Len= 20 IPos= 0 Q= 140735655680592
Error termination in NtrErr:
NtrErr Called from FileIO.

复制代码

按上面方法看似实际有效的报错为“ONIOM data not found on unit 2.”，但其实这根本不是个ONIOM任务，输入文件与ONIOM毫无关系。
其看似报出ONIOM错误的原因是，在正常的opt=restart任务中，程序会在尝试读取 ONIOM data 之后进行标题和Route的读取（见下面的样例），故上述错误很可能是在此处读取这些信息时出错。故也可通过不读取chk文件解决。分析出这个问题需要了解类似的正常任务中“将会有什么输出”，难度较大。

Restoring state from the checkpoint file "/home/gauuser/Gaussian/PhS_Rad_Cyc/2_Z
E_transform/TS[Complete_M062X]_Step2[TS_opt]cont.chk".
ONIOM data not found on unit 2.
Title: TS [EXTRACT_GEOM]:2,7,12
Route: #p m062x/genecp opt=(noeigentest,ts,calcfc,nofreeze) int=ultrafine freq
empiricaldispersion=gd3 geom=allcheck 5d 7f scrf=(smd,solvent=acetonitrile) gues
s=tcheck
FncErr=1.00D-07 GrdErr=1.00D-06

复制代码

还有些杂七杂八的情况，诸如我在做片段初猜的计算时，误将多重度 0 1 0 2 0 -2 输入为 0 1 0 3 0 -3，结果即有如下的错误

(Enter /home/gauuser/g09/l122.exe)
Structure from the checkpoint file: "sideroidide_CS_1_opt.chk"
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -584 Len= 36 IPos= 0 Q= 47012767336112
dumping /fiocom/, unit = 1 NFiles = 33 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 1206272 FType=2 FMxFil=10000
Number 0 501 502 503 507 511
Base 216576 43520 72192 180736 181760 196608
End 1206272 44520 76297 180821 181925 197835
End1 1206272 44544 76800 181248 182272 198144
Wr Pntr 216576 43520 72192 180736 181760 196608
Rd Pntr 216576 44520 72192 180736 181760 196608
Length 989696 1000 4105 85 165 1227
......
dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512
defal = T LstWrd = 65536 FType=2 FMxFil=10000
Number 0
Base 40448
End 65536
End1 65536
Wr Pntr 40448
Rd Pntr 40448
Length 25088
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -584 Len= 36 IPos= 0 Q= 47012767336112
Error termination in NtrErr:
NtrErr Called from FileIO.

复制代码

此时只能根据出错的Link数来推测、测试，得知真正为多重度设定问题。

总结一下我解决这个报错大致思路：

（1）去除所有需要读chk、restart的操作尝试
（2）在大堆数字的上面找到真正的报错，找有效信息
（3）输入文件中用#p代替#，观察出错的Link，思考相关原因
（4）其他杂七杂八问题，硬盘满了，权限问题，路径中有奇怪的字符等等
（5）若在优化过程中出现这个问题，可以将当前结构保存重新开始一个任务，偶尔能解决。
（6）Gaussian Bug，如果排查了上述问题、做了化简，换了新版本、机器、跑多次都能重复这个错误，参阅下文关于Gaussian Bug的说明。

如果做了上述尝试还未能解决这个问题，要向别人提问的话请提供完整的输入输出。

————————————————————————————————

Error: illegal instruction, illegal opcode

这一错误不在输出文件中报错，而是在 Linux 终端中显示，常见样式如下：

Error: illegal instruction, illegal opcode
rax 000000001b23ef84, rbx 0000000000730438, rcx 000000001b23ef84
rdx 00007fffa89b0b10, rsp 00007fffa89b0b40, rbp 00007fffa89b0b50
rsi 00000000007c29f0, rdi 00007fffa89b0b10, r8 0000003e2d78fee8
r9 0000000000000001, r10 00007fffa89b0880, r11 0000000000000202
r12 00007fffa89c7878, r13 00007fffa89c7848, r14 0000000000000000
r15 0000000000000000
/lib64/libpthread.so.0() [0x3e2e00f710]
/home/gauss/g09/l1.exe() [0x4a72e1]
/home/gauss/g09/l1.exe() [0x41d295]
/home/gauss/g09/l1.exe() [0x403745]

复制代码

输出文件中没有有效信息，仅输出初始命令

Entering Gaussian System, Link 0=g09
Input=CO.com
Output=CO.log
Initial command:
/home/gauss/g09/l1.exe "/home/gauss/g09/tmp/Gau-4412.inp" -scrdir="/home/gauss/g09/tmp/"

复制代码

成因：
Gaussian的可执行文件有不同的指令集版本，不同指令集的CPU需购买不同版本的Gaussian程序（具体可见：http://gaussian.com/g16/g16_plat.pdf）

自己的CPU支持什么指令集可以在Linux上使用命令“cat /proc/cpuinfo”查看，若Gaussian的版本是“AVX2-enabled”，则仅当 “/proc/cpuinfo” 文件的 flags 字段中支持AVX2时才能使用，否则应该更换其他指令集版本的Gaussian程序，如 Legacy 即 pre-SSE4.2 版本

解决：

更换Gaussian的可执行文件的指令集版本，如将Gaussian 09 E.01 X86_64(AVX2-enabled) （兼容性差、但性能高）更换为 Gaussian 09 E.01 X86_64(“Legacy”-pre-SSE4.2)（下兼容性好，但性能稍差）
更换机器、使用最新的CPU

————————————————————————————————

Gaussian Bug 的“样式”

触发Gaussian的Bug时常有类似这样“模式”的报错：

Raff turned off since only 57.62% of shell-pairs survive.
Unable to match L and R vectors in BiOrth.

复制代码

Coeff: 0.000D+00 0.000D+00 0.000D+00 0.100D+01
Logic error in SftOpn.

复制代码

即以某六字母变量为中心语的某些负面描述
这时就有较大的Bug嫌疑了

另外部分Ntrerr也是Bug

Error termination in NtrErr:
NtrErr Called from FileIO.

复制代码

遇到“Bug”时

反思自己有没有瞎搞高斯的临时文件。诸如：跑着跑着临时文件被删了；两个文件用同一个CHK、RWF文件名及路径在同时跑（将某个输入文件复制改写为其他结构时常常发生）；用脚本比较深层次的控制高斯可能引起的一些奇怪问题（比如用脚本调用 Gaussian external 读梯度时若两个临时文件冲突，会诱发 L601 的 RdWrB1 read garbage pointers 报错。）
首先应更换当前最新Gaussian版本（发稿时为Gaussian09 E.01）
如果仍然报同样问题，可询问高斯客服。
如果不愿询问（如没有版权等），只能更换方法/基组、twitch一下结构、更换积分格点、如果是激发态计算将TD换为TDA等等杂七杂八的方式尝试绕开它了、

附：已知的、可能为Bug的报错：
            L103，Bad arguments to LAPack or BLAS routine.
            L401，Diagonalization in DiagDN via DSPEV failed.
            L906，Internal consistency failure #1 in GetIJB. (G09 E.01(不含)之前版本)
            L914，Unable to match L and R vectors in BiOrth
            L1002，NIJ > Max2 in MMCore. (G09 C.01(不含)之前版本)
            L1014，Tx not orthogonal to T

            NtrErr Called from FileIO 的一部分情况

————————————————————————————————

无报错的异常情况

IRC提前停止或折返

(, 下载次数 Times of downloads: 380)

————————————————————————————————

优化后的结构有一个或几个小虚频

有时优化后会出现一小虚频，虚频数值在几至几十cm^-1。

(, 下载次数 Times of downloads: 398)

解决：

若满足：

(1) 使用的是Gaussian 09 E.01之前的版本，如Gaussian 09 D.01
(2) 使用了D3BJ色散校正（关键词为empiricaldispersion=GD3BJ）

则可选择在优化和频率计算中同时：
(1) 换用Gaussian 09 E.01或之后版本
(2) 换用empiricaldispersion=GD3
(3) 不用色散校正

这个Bug的具体解释见：http://bbs.keinsci.com/thread-2410-1-1.html

若不满足上述两个情况：

看《Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法》（http://sobereva.com/278）

————————————————————————————————

G09W 在 L1.exe 左右无报错退出

(, 下载次数 Times of downloads: 340)

(, 下载次数 Times of downloads: 378)

成因：目前网上能找到的G09W都是32位程序，故最多可设置使用约1300MB内存（%mem=1300MB，实际极限在1380~1600之间），4 核并行，16GB 硬盘（切割Scratch file无效）。如果超过此值可能触发不报错退出。如上图即为故意将内存设为 %mem=1500MB 时的“无报错错误输出”。

解决：

将内存、核数、硬盘设置到相应极限以下
有钱的去买64位的Windoes G09 程序（ ╮(╯▽╰)╭ 然后打个包给大家就更好了）
换用Linux下的Gaussian程序

————————————————————————————————

GaussView观看结构有异常的成键、断键

成因：首先“键”并不是原生的量化概念，量化的输入、输出也（几乎）不依赖键连信息，显示成键情况只是为了辅助理解而已。很多非经典的情况下经典意义上的“键”并不是一个好的概念。

Gaussview 中的显示的成键状态（在未人为规定时）是用原子间距离判断的，即两原子间键长在这对儿元素（比如 C 和 C ）对应的某个范围内（如145 pm-165 pm），即判定为相应键级的成键（如C-C单键），这个成键情况的判断显然是十分粗糙的，不必在意。

常见的成键、断键异常有这几种情况

下图中H2O2是优化后，O-O 键键长变得较长，显示 O-O 之间断键了，但测量 O-O 原子间距离可知仅 145pm，做进一步的其他键级分析也可以得出其中显然有一个共价键，故 “O-O断键” 只是显示问题而已，并非结构优化等过程有错。
中间的例子是一个芳环，芳环键的判定常有错误（不显示为共轭的“1.5级键”），经常判断为全是双键、或完全混乱等等，同样不必在意。
第三个例子较为特殊，其打开的是一个fchk（或chk）文件。在优化过程中，Gaussian会在fchk文件中保存第一步优化/IRC 时判断出的成键信息 / 或 geom=connectivity 定义的连接信息，故不论后续优化过程如何进行，打开时GaussView都会显示优化/IRC 刚刚开始时的成键信息。第三个例子就是一个IRC过程，在初始结构中，1-2两个氮原子很近，根据原子间距离判断图中1-2两个氮原子间有单键，IRC/优化后虽然 1-2 两个氮原子远离，但第一步的键连关系被保留，显示为一根异常长的单键。观察这个结果中的键长、键角，根据常识判断是合理的，说明其也是正确的，不说明计算有问题。如果打开的是Log文件，且没有记录geom=connectivity则没有这个显示问题

(, 下载次数 Times of downloads: 320)

解决：GaussView （及其他很多软件，如Spartan）显示的成键情况很粗糙，不论是否出现未预期的结果，都应自己通过化学常识判断几何优化过程、显示的成键情况等是否合理。另外可以用Multiwfn计算Mayer键级来辅助判断。

对于上图第三种情况，可以使用下图所示的Rebond工具，让GaussView根据当前键长重新粗略判断成键情况，可以去除一些明显错误的异常成键

(, 下载次数 Times of downloads: 304)

————————————————————————————————

GaussView 中显示的 IRC 在零点附近出现异常的低能点，即 “第一个点就掉下去了”

症状为在 GaussView 中打开 Result - IRC/Path 时，在反应坐标零点处附近出现一低能点，使图形看起来是“第一个点就掉下去了”：

(, 下载次数 Times of downloads: 338)

成因：GaussView 显示 Bug。

仔细观察可发现，如下图所示，只要将“掉下去的”点移动到曲线的一端（红叉所示，左右不定），曲线就正常了（绿线所示）。

(, 下载次数 Times of downloads: 357)

其原因在于 GaussView 在读取未正常结束的（尚未跑完、跑出错）的 IRC 时，当前点的反应坐标尚未计算出（如下面段落所示的信息缺失），但当前结构及其能量却已经得到并被GaussView读取，故GaussView 默认其横坐标为0，导致点的顺序出错，显示为这一问题。

CHANGE IN THE REACTION COORDINATE = 0.30843
NET REACTION COORDINATE UP TO THIS POINT = 2.77474
# OF POINTS ALONG THE PATH = 19
# OF STEPS = 1

复制代码

解决：如果任务正在运行，使其跑完
如果任务已经跑完，应是出了什么问题没有正常结束，看输出文件具体报错解决。

————————————————————————————————

不是报错、但常被误以为报错的语句

Warning!!: The largest alpha MO coefficient is... ；Warning!!: The smallest alpha delta epsilon is...

(Enter /home/ki/g09/l801.exe)
Windowed orbitals will be sorted by symmetry type.
ExpMin= 2.53D-02 ExpMax= 3.74D+05 ExpMxC= 1.18D+03 IAcc=3 IRadAn= 5 AccDes= 0.00D+00
HarFok: IExCor= 205 AccDes= 0.00D+00 IRadAn= 5 IDoV=-2 UseB2=F ITyADJ=14
ICtDFT= 12500011 ScaDFX= 1.000000 1.000000 1.000000 1.000000
Largest valence mixing into a core orbital is 4.82D-04
Largest core mixing into a valence orbital is 1.44D-04
Range of M.O.s used for correlation: 6 101
NBasis= 101 NAE= 9 NBE= 9 NFC= 5 NFV= 0
NROrb= 96 NOA= 4 NOB= 4 NVA= 92 NVB= 92
**** Warning!!: The largest alpha MO coefficient is 0.21049552D+02
**** Warning!!: The smallest alpha delta epsilon is 0.84363472D-01
Leave Link 801 at Mon Sep 8 14:14:12 2014, MaxMem= 805306368 cpu: 0.7
(Enter /home/ki/g09/l804.exe)

复制代码

这两组语句是在SCF结束后发现最大分子轨道的系数太大、以及分子轨道间的能量差太小。这一检查的意义是在post-HF（如MP2）、TDDFT等计算中存在以分子轨道系数为被除数、分子轨道间能量差为除数的项，所以分子轨道系数太大、或分子轨道Gap接近0时，这一除法的结果将得到一个较大数值、从而可能引发数值不稳定性，导致结果有误。故而设置这一警告。在基组存在近似的线性相关问题时，如使用了弥散函数，出现这一警告的概率较大。

但这个警告的阈值被设置的较为严格，在正常的计算中也会报出这一警告。如果确认结果基本合理，则可以忽略这两个Warning。
若结果不合理，可：

（1）自己更换或修改基组，尤其是在不需要弥散函数时去掉弥散函数，防止出现接近线性相关的问题
（2）IOp(3/59)可能可以控制检查基组线性相关问题的阈值，可将其放宽尝试
(3) 参见上文中关于“L801，Fatal problem: The smallest alpha delta epsilon is XXXXX”的解决方案。
某些错误常会在该提示后面发生，而Warning前后有空行，颇引人注意，因而常有人误以为这是错误的根源。例如如下段落的报错是“EpsInf not defined for this solvent.”，而不是几个硕大的Warning。

**** Warning!!: The largest alpha MO coefficient is 0.30204407D+03
**** Warning!!: The smallest alpha delta epsilon is 0.88447585D-01
**** Warning!!: The largest beta MO coefficient is 0.30099143D+03
**** Warning!!: The smallest beta delta epsilon is 0.91819897D-01
Leave Link 801 at Fri Nov 6 19:09:46 2015, MaxMem= 1048576000 cpu: 0.3
(Enter /software/gs09/g09/l1002.exe)
Minotr: UHF open shell wavefunction.
NEqPCM: Using equilibrium solvation (IEInf=0, Eps= 13.8200, EpsInf= 0.0000)
EpsInf not defined for this solvent.
Error termination via Lnk1e in /public/software/gauss/g09/l1002.exe at Mon May 11 07:58:05 2015.

复制代码

附：高斯官方对此问题的回复：

Gaussian tests the values of the MO coefficients and orbital energies after each SCF calculation. When you get a warning about the smallest delta epsilon, the warning is referring to the difference between orbital energies. Since some of the equations in post-HF methods include a term in the numerator which has the MO coefficients to some positive integer and a difference in orbital energies in the denominator, the calculations can become numerically unstable if the numerator is too large (e.g. MO coefficients too large) or the denominator is too small (orbital energy difference too small). So, this is a warning that the post-HF results may be affected. However, the criteria for these warnings are fairly strict, so the vast majority of the calculations that have these warning messages do not actually suffer from numerical instabilities. If you have results that don't seem to make sense, and these warnings are present, it would be a good idea to repeat the calculation with a different basis set (these problems tend to be worse for basis sets with many diffuse functions) to verify the results.

复制代码

————————————————————————————————

Error on total polarization charges；DIIS: error

Cycle 11 Pass 1 IDiag 1:
RMSU= 6.48D-09 CP: 1.00D+00 1.01D+00 1.00D+00 1.05D+00 1.06D+00
CP: 9.80D-01 9.86D-01 1.06D+00 1.10D+00 9.43D-01
E= -12097.4433058840 Delta-E= 0.000000000164 Rises=F Damp=F
DIIS: error= 4.95D-08 at cycle 11 NSaved= 11.
NSaved=11 IEnMin=10 EnMin= -12097.4433058841 IErMin=11 ErrMin= 4.95D-08
ErrMax= 4.95D-08 0.00D+00 EMaxC= 1.00D-01 BMatC= 9.84D-12 BMatP= 5.98D-11
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.204D-04 0.138D-05 0.105D-03 0.650D-03 0.539D-03-0.301D-02
Coeff-Com: -0.154D-01-0.324D-01 0.431D-01 0.322D+00 0.684D+00
Coeff: -0.204D-04 0.138D-05 0.105D-03 0.650D-03 0.539D-03-0.301D-02
Coeff: -0.154D-01-0.324D-01 0.431D-01 0.322D+00 0.684D+00
Gap= 0.206 Goal= None Shift= 0.000
RMSDP=3.88D-09 MaxDP=5.72D-07 DE= 1.64D-10 OVMax= 1.85D-06
Error on total polarization charges = 0.03783
SCF Done: E(RM052X) = -12097.4433059 A.U. after 11 cycles
NFock= 11 Conv=0.39D-08 -V/T= 2.0073
KE= 1.200978298414D+04 PE=-6.824880392971D+04 EE= 2.436187738767D+04
SMD-CDS (non-electrostatic) energy (kcal/mol) = -8.14
(included in total energy above)
Leave Link 502 at Sat Feb 11 02:31:20 2017, MaxMem= 18352963584 cpu: 8769.1

复制代码

如上面完全正常的任务输出段落所示，这两个语句都不是报错。
认为它是报错其实是英语问题。。。此处 Error 分别指总极化电荷和DIIS的误差，而不是错误。

————————————————————————————————

No special actions if energy rises.

(Enter /gpfs/share/home/g16/l502.exe)
Keep R1 ints in memory in canonical form, NReq=136241796.
FoFCou: FMM=F IPFlag= 0 FMFlag= 0 FMFlg1= 0
NFxFlg= 0 DoJE=F BraDBF=F KetDBF=F FulRan=T
wScrn= 0.000000 ICntrl= 600 IOpCl= 0 I1Cent= 0 NGrid= 0
NMat0= 1 NMatS0= 14878 NMatT0= 0 NMatD0= 1 NMtDS0= 0 NMtDT0= 0
Symmetry not used in FoFCou.
Two-electron integral symmetry not used.
Closed shell SCF:
Using DIIS extrapolation, IDIIS= 1040.
NGot= 9895542784 LenX= 9784774425 LenY= 9784735568
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on energy=1.00D-06.
【No special actions if energy rises.】
Cycle 1 Pass 1 IDiag 1:
Inv3: Mode=1 IEnd= 1030188.
Iteration 1 A*A^-1 deviation from unit magnitude is 3.00D-15 for 375.
Iteration 1 A*A^-1 deviation from orthogonality is 2.66D-15 for 435 266.
Iteration 1 A^-1*A deviation from unit magnitude is 2.66D-15 for 375.
Iteration 1 A^-1*A deviation from orthogonality is 2.66D-15 for 379 156.

复制代码

此语句不是报错，它只是在描述SCF过程中遇到该步骤的能量比上一步骤高的时候，程序该做什么。
这一语句的选项可以由 IOp(5/86) 控制，例如设置 IOp(5/86=101202) 时会输出下列内容（默认为 IOp(5/86=101100)）：

Reduce DIIS space if energy rises from previous iteration.
Dynamic level shift is off after energy rises.

复制代码

初学者常以为这句话是报错的原因是程序在输出这句话之后会进行SCF迭代，耗时较长，如果初学者不加#p输出SCF迭代的详细信息，那输出文件末尾会保留在“No special actions if energy rises”很久，令初学者误以为出错了。平时计算时都应该加上#p来输出具体SCF信息。

————————————————————————————————

This type of calculation cannot be archived.

(Enter /home/gauuser/g09/l9999.exe)
This type of calculation cannot be archived.
Job cpu time: 1 days 8 hours 55 minutes 51.4 seconds.
File lengths (MBytes): RWF= 1896 Int= 0 D2E= 0 Chk= 138 Scr= 1
Normal termination of Gaussian 09 at Sat Feb 11 03:43:47 2017.

复制代码

此“报错”出现于任务正常结束前。“Archive” 表示的是Gaussian在一般任务的最后常会输出一大段乱码一样的文本、用于存储任务的各种关键信息。这里“cannot be archived”仅表示 Gaussian 对于 IRC 等类型任务未使其输出存档段信息。

典型的Archieve段落：

1|1|UNPC-LIYUANHE-PC|FOpt|RB3LYP|6-311+G(d,p)|C1H4|LIYUANHE|10-May-201
1|0||# opt freq b3lyp/6-311+g(d,p) geom=connectivity||CH4 OPT FREQ||0,
1|C,6.705202501,1.156069305,0.|H,7.0687325183,0.127811533,-0.000000934
4|H,7.0687519538,1.6701910498,0.8904935408|H,7.0687504279,1.6701921288
,-0.8904935408|H,5.614575104,1.1560825085,0.0000009344||Version=IA32W-
G09RevA.02|State=1-A1|HF=-40.5339325|RMSD=6.153e-009|RMSF=9.920e-005|D
ipole=0.,0.,0.|Quadrupole=0.,0.,0.,0.,0.,0.|PG=TD [O(C1),4C3(H1)]||@

复制代码

————————————————————————————————

End of XXXXXX F.D. properties file XXX does not exist.

FullF1: Do perturbations 1 to 3.
Isotropic polarizability for W= 0.000000 466.35 Bohr**3.
End of Minotr F.D. properties file 721 does not exist.
End of Minotr F.D. properties file 722 does not exist.
End of Minotr F.D. properties file 788 does not exist.
Leave Link 1002 at Fri May 3 18:59:51 2019, MaxMem= 10995105792 cpu: 153536.8 elap: 4809.4
(Enter /gpfs/share/home/1501110295/g16/l601.exe)

复制代码

G2PCM: DoFxE=T DoFxN=T DoGrad=T DoDP/DQ/DG/TGxP=FFFF NFrqRd= 0 IEInf=0 SqF1=F DoCFld=F IF1Alg=4.
GePol: Maximum number of non-zero 1st derivatives = 301
End of G2Drv F.D. properties file 721 does not exist.
End of G2Drv F.D. properties file 722 does not exist.
End of G2Drv F.D. properties file 788 does not exist.
Leave Link 1110 at Fri May 3 06:07:58 2019, MaxMem= 10995105792 cpu: 17529.5 elap: 549.1
(Enter /gpfs/share/home/1501110295/g16/l1002.exe)

复制代码

这几个语句的含义不明，但不是报错。其常被用户误以为是报错的原因是其出现在计算频率时较为耗时的进程前面，对较耗时的体系没有新的输出，故而使用户误以为“跑到这里就卡死了”，但其实经常是高斯正在进行计算，确认高斯还在运行中之后等待算完即可。

————————————————————————————————

其他常见问题

GaussView 打开 chk 或 fchk 文件报错 Missing or bad data: Alpha Orbital Energies

解决：两种方法 (1) 先将chk转换为fch文件，用Multiwfn软件载入fch文件，用主功能100-选项2转换为fch文件给GaussView读入。或者：（2）先将chk转换为fch文件，然后用文本编辑器打开fch文件，看看里面的Number of basis functions和Number of independent functions是否一致。如果不同，将independent改成independant

————————————————————————————————

估计耗时，“这种计算大约需多长时间？”

常有人给出一个体系询问这个计算大概需要多长时间，这种问题往往难以回答。

计算所需的时间取决于很多因素，包括而不限于下列因素（请注意，只是举例，即使依次回答下面的问题也常不足以给出靠谱的猜测）：

体系的大小（原子数与元素构成）
方法（除post-HF和DFT/HF这样的比较外，还例如开RI时杂化泛函比纯泛函慢）
基组数量、种类（相同数量的弥散函数远比极化函数耗时多）；轨道数
电子态及电子结构的复杂程度（例如一般分子需要10~20步SCF，稀土有时需要几百步SCF；同样的基组数，一般有机物也远比金属簇容易）
任务类型
初猜的好坏
任务的顺利程度（如优化步数、SCF步数）
加速或“减速”设定（如RI, int=ultrafine等）；
电脑的CPU（核数、频率、架构）
Gaussian版本（例如avx2就比avx版本快一些，更极端的还可能涉及数值导数和解析导数的差别）
内存、硬盘情况（大小、速度）及为任务分配了多少，是否有其他任务在共享，使用N个核的超线程时是用了N个线程还是2N个线程等等等等琐碎的问题。

以上任何一个因素改变都可能显著的、在数量级层面上改变计算所需的时间，同时也不会有人了解所有的机器环境（用过各种CPU之类）、所有的任务和体系。故往往此类问题根本没法回答。

比较合理的估计耗时的方法是在自己机器上测试计算中某个步骤的耗时，用它估算。有限的可供参考的“感觉”大致如下：
（我主要用（非双杂化的）DFT，如B3LYP、M06-2X之类，下面仅对这类方法、仅对一般有机体系适用。阅读下文时请注意区分“一步SCF迭代”，“（整个）SCF过程”，“一步几何优化”，“整个几何优化过程”的区别）

一次单点的耗时，可以通过一步SCF迭代进行估计。初次SCF一般在十几步至几十步之间收敛，超过50步应开始检查SCF震荡；优化时后续优化步骤的SCF会读取之前的初猜，故优化的第二步开始，SCF所需的步数常会显著减少，一般在几次至十几次（常可见优化中第一步L502耗时最多，后续L502耗时较少）。可根据这一经验、通过一步SCF迭代大致乘出整个SCF过程的耗时。
解析梯度计算的耗时很少，故几何优化的耗时（等于一次单点+一次梯度，不含calcfc）一般只比单点高一点，可略。
几何优化任务，对小体系一般在二三十步以内优化收敛，大体系（100+原子）在50~150步都有可能，再多需谨慎考虑优化震荡的可能性。另外上面说过优化中第一次SCF过程常常比后续SCF所需的时间长很多，所以估计几何优化的耗时使用第二步几何优化的时间，乘以可能需要的几何优化步数来估计比较合适。
频率计算一般需要单点几倍至几十倍的时间，少则3倍，多则近100倍也遇到过，这一倍数似乎与几何的变量数目有关系。可用这个规律估计各处计算频率的耗时（含优化完成后、使用calcfc时优化第一步、以及calcall的每一步涉及的频率计算）

其他可能（严重）影响耗时的情况：

SCF、优化震荡：做大任务时，应设法实时监测SCF和优化的步骤，可以自己写一些小工具，不要提交个任务，跑了两个星期跑到500步报错了，才发现从20步开始就在震荡了。
其他理论方法：一些有解析梯度、频率的 post-HF方法（含双杂化泛函如B2PLYP），其解析梯度和频率耗时与单点的比例一般高于HF和DFT方法，解析梯度可至单点的几倍，解析频率
复杂电子结构：有的体系的电子结构复杂，诸如含镧系原子的分子，其SCF收敛可能需要几百步，远多于一般体系，应特殊处理（此时优化的时候用calcall甚至可能比不用的时间还短）
使用特殊选项：诸如使用 scf=xqc时，有时可以用常规方法收敛，有时不行；如果使用了qc，会使得一次SCF收敛的时间大大增长，以往的“SCF需要多少步”的经验也不再适用
数值梯度和数值频率：对于没有数值梯度的方法，诸如CCSD(T)，每步几何优化的耗时等于6N步单点能（N为原子数），此时梯度计算所需的时间显然不可忽略；对没有解析Hessian的方法，计算一次频率的耗时为6N次[单点+梯度]

对这个问题有不同或其他意见欢迎补充

————————————————————————————————

如何通过柔性扫描产生较好的过渡态初猜

该方法是寻找过渡态的备用方法，并非一定成功、且耗时甚高，这绝不是找过渡态的“标准方法”，使用者稍有提供初猜的经验时、绝大部分过渡态都可以用opt=TS寻找到。

柔性扫描是有效的获得过渡态初猜的傻瓜方法，在初学者刚接触某类反应、不了解大致的过渡态结构（而且还懒得去查文献的时候）时、用此法不必动很多脑筋、比较robust；在研究有些较难找到的过渡态，比如某些自由基过程时也比较有效。柔性扫描很简单、但却常有人问，故在此说明。请务必注意本段最后的注意事项。

我们以下列自由基氢转移反应为例：

(, 下载次数 Times of downloads: 358)

在GaussView中建模，首先创建一个甲烷，使用Add-Valence按钮点击其中一个H，将新出现的氢原子取代为SH基团。基本建模就完成了。

(, 下载次数 Times of downloads: 349)

随后需要调整键长。这个反应中涉及两个键、分别是 S-H 和 H-C 键，但我们只需扫描其一，因为在优化过程中另一个键长会随之弛豫变化。

如果自己对过渡态中的键长没有较好把握的话，建议做两次扫描，均以自己猜想的过渡态为起点，这有两点好处，一是可以在扫描过程中监测扫描过程，发现已有最高点时可立即停止计算，节约计算量；二是这样在优化过程中另一边不会“飞掉”，如本例中，若 C-H 键键长较小，这个体系就变成了分子间弱相互作用体系，只用一次扫描时容易在开始阶段优化不收敛、或因为初猜原因跑到不希望要的构型上。

这个反应中我们以扫描反应中涉及的 C-H 键为例，正常的C-H键键长在110 pm左右；过渡态中势必会较长，应当根据经验适当猜测（也可查阅文献中同类反应过渡态的结构），此处假设所猜数值为170 pm。

(, 下载次数 Times of downloads: 371)

随后在 Edit-Redundant Coordinate 中设置柔性扫描 C-H 键，扫 5 步，每步 15 pm （即 170 ~ 245 pm 的 6 个结构），用 B3LYP/6-31G(d) 为例提交运算；再重新设置每步为 -15pm，扫描 3 步（即 125 ~ 170 pm 的 4 个结构），提交运算。

(, 下载次数 Times of downloads: 376)

GaussView产生的第二个输出文件如下：

%nprocshared=2
%mem=1200MB
%chk=C:\Users\LiYuanhe\Desktop\1.chk
#p opt=modredundant freq b3lyp/6-31g(d)
Title Card Required
0 2
C -3.01654961 -1.31556721 -0.25719647
H -2.65989518 -2.32437721 -0.25719647
H -2.65987677 -0.81116902 -1.13084798
H -4.08654961 -1.31555402 -0.25719647
H -2.44987313 -0.51418690 1.13084798
S -2.01319891 0.10334733 2.20045870
H -2.44953630 1.33854341 2.20102655
B 5 1 S 5 0.150000

复制代码

其中 opt=modredundant 是进行柔性扫描的关键词，分子段落后面的 “B 5 1 S 5 -0.150000” 表示增加一个冗余内坐标，内容是一根化学键（B），在原子 5 和原子 1 之间，模式为扫描（S），共扫描 5 步，步长为向减小键长方向的 0.15 A。

计算完后可在 GaussView 中 Result-Scan 分别查看两次任务的曲线。为方便考察，这里将两个结果做到一条曲线上：

(, 下载次数 Times of downloads: 323)

可见其最高点在 1.7 A 处。故我们在 GaussView 中取 1.7 A 的扫描结果作为过渡态初猜，用opt=TS，在同一水平下找过渡态，近4步优化就找到了正确的过渡态，过渡态 C-H 键键长 1.675 A。

%nprocshared=2
%mem=1200MB
#p opt=(calcfc,ts,noeigentest) freq b3lyp/6-31g(d)
Title Card Required
0 2
柔性扫描最高点的分子坐标

复制代码

(, 下载次数 Times of downloads: 354)

（如果你完全按上述设置，会发现扫描时最左侧一点时会遇到了L9999 的 opt 未收敛错误，可以按照本帖上面提到的方法解决，但也可无视，因为此时打开输出文件可见路径上的最高点已经找到了，就不用算更多点了）

几点说明：

该方法应作为寻找过渡态的备用方法（比如上面反应的自己画初猜用 opt=TS 其实也容易找到），而不是每次找过渡态都用它。该方法较为耗时，但还可以接受。例如扫描时取10个点时，因为后续的优化是用前面的结构做初猜的，其并不等于做十个优化任务，而只有第一步是完整的优化、后续优化步大概只需要一般优化任务的一半甚至更少的步数和时间，故柔性扫描10个点一般也只是2~5倍优化任务所需的时间，对不十分大的任务还是容易接受的。
需强调使用的是柔性扫描（opt=ModRedundant），不是刚性扫描（Scan），后者对找过渡态几乎没用。
做柔性扫描的级别可以比实际寻找过渡态的级别稍低，但不能太差。例如过渡态用PBE0/def-TZVP计算，则柔性扫描最低可以在 PBE0/def2-SV(P)、或B3LYP/6-31G(d)级别下做；而如果用HF/6-31G(d)这样的级别耗时不低、得到的结果却帮助有限。
柔性扫描只是帮助产生过渡态初猜，而不是“寻找过渡态”。柔性扫描的最高点（即使步长无限小）也几乎都不是过渡态。柔性扫描的曲线有时与IRC相差甚远。
步长 10 pm 左右基本是够用的，如果势能面比较复杂，对初猜十分敏感，也可以先用较宽的步长扫描，然后对感兴趣的区间用更小的步长扫描（如上述例子中、可以最后再用0.03的步长在1.5~1.8之间扫描）
本文的扫描和过渡态寻找都用 B3LYP/6-31G(d) 只是举个例子，实际计算中应当选取合适的计算级别。扫描和寻找过渡态可以用不同的级别，但如果两种方法的势能面差别太大，可能效果不好，若体系、级别很耗时，可如先用 B3LYP/6-31G(d) 做上述扫描，最后再用高级别在最高点附近确认三个点，如在本例中用 M06-2X/def-TZVP 重新确认 1.6 A， 1.7 A， 1.8 A 三个点，确认大小关系正确（方法是在 Edit-Redundant Coordinate 中选择 Freeze 而不是 Scan）。
注意柔性扫描过程中会出现因结构跳变造成的Artifact，体现为曲线不平滑、有“悬崖”，如下图所示。注意此时曲线上的极大值点与过渡态毫无关系，需为连续、平滑的曲线（如上面的例子中所示），才体现过程中确实存在过渡态、且可作为初猜。遇到这种扫描曲线不平滑的情况应视扫描轨迹的具体情况解决问题，尝试避免跳变的发生。

(, 下载次数 Times of downloads: 335)

————————————————————————————————

如何在量化计算中设定原子/片段的电荷

看 http://sobereva.com/82 一文第六节

————————————————————————————————

如何计算 Activation Strain (Distortion-Interaction) Energy

Activation Strain 模型和 Distortion-Interaction 完全是一个东西。前者由 Bickelhaupt (1999) 基于前人的一些想法建立起来，提出 Activation Strain 模型的名称；后来 K. N. Houk 在 2007 年 JACS 上发表了一篇文章，未引用 Bickelhaupt，只引用之前更早的一些文章，称之为 Distortion-Interaction，内容与 Bickelhaupt 是完全相同的，没有任何新东西。但或许只是因为他有名，故而大家多称之为 Distortion-Interaction Model。在此鄙视这种行为。

Distortion-Interaction 模型是一种对活化能进行能量分解的方案。对于一个双/多分子反应（单分子反应不适用），distortion能量定义为过渡态各组分在过渡态构型下、无相互作用时的能量与稳定底物之间的电子能差，此值一般为正值；interaction 定义为前述无相互作用体系与真实的过渡态电子能差，此值一般为负值。二者之和应精确等于活化能。
以下述Diels-Alder反应为例。

(, 下载次数 Times of downloads: 322)

定义如下能量。E1、E2为相应物质优化后（opt）的电子能。E5为过渡态优化后的电子能（opt=TS）。随后将过渡态优化后的结构提取出来，删去其中的蓝色部分的几个原子，不优化，在同样级别下计算单点能（sp），得到能量E3；类似的删除红色部分，计算单点能得到E4。

(, 下载次数 Times of downloads: 373)

则有：

活化能 Ea = E5-E1-E2；

Distortion Energy = E3 + E4 - E1 - E2；

Interaction Energy = E5 - E3 - E4。

可见 Ea = Interaction Energy + Distortion Energy，是一个精确地能量分解。

这个能量分解方案简单，不需要任何特殊技术，耗时相对优化也很低。有时可以为分析提供一定角度。

除了对过渡态做这件事之外，也可以对IRC上的每个点做完全一样的能量分解。观察 distortion 和 interaction 在“反应过程中”的变化。

————————————————————————————————
▣

作者
Author: 978142355 时间: 2017-1-1 18:35
Li老师，L906的问题也总能看到，好像一直都是建议用G09E01版本。
Internal consistency failure #1 in GetIJB.
Error termination via Lnk1e in /usr/local/share/g09/l906.exe。
下面的这个错误：
End of file reading connectivity.
Error termination via Lnk1e in /usr/local/share/g09/l101.exe
内存给的小了：
Out-of-memory error in routine RdGeom-1 (IEnd= 6750021 MxCore= 0)
Use %mem=7MW to provide the minimum amount of memory required to complete this step.
Error termination via Lnk1e in /usr/local/share/g09/l101.exe

作者
Author: liyuanhe211 时间: 2017-1-1 20:02

978142355 发表于 2017-1-1 18:35
Li老师，L906的问题也总能看到，好像一直都是建议用G09E01版本。
Internal consistency failure #1 in Ge ...

都补充了，谢谢！

作者
Author: liyuanhe211 时间: 2017-1-1 20:02

978142355 发表于 2017-1-1 18:35
Li老师，L906的问题也总能看到，好像一直都是建议用G09E01版本。
Internal consistency failure #1 in Ge ...

都补充了，谢谢！

作者
Author: sobereva 时间: 2017-1-1 20:48
关于L801，Excessive mixing of core and valence orbitals.这个，建议说成把full写在理论方法关键词里面。要不然有人做TDDFT或CCSD的时候碰见这个错误，会傻乎乎地直接额外写个MP2(full)

作者
Author: liyuanhe211 时间: 2017-1-1 21:57

sobereva 发表于 2017-1-1 20:48
关于L801，Excessive mixing of core and valence orbitals.这个，建议说成把full写在理论方法关键词里面。 ...

已更正。谢谢指出！

作者
Author: 978142355 时间: 2017-1-2 08:31

liyuanhe211 发表于 2017-1-1 20:02
都补充了，谢谢！

客气了，Li老师~

作者
Author: wudazhuang509 时间: 2017-1-10 12:00
老师好，请问下一般计算EPR需要使用什么基组比较合适呢？我计算的是NaCl晶体掺Rh的g因子。看过您的帖子谈谈量化中基组的选择，是需要和NMR计算一样的基组么？

作者
Author: sobereva 时间: 2017-1-10 12:14

wudazhuang509 发表于 2017-1-10 12:00
老师好，请问下一般计算EPR需要使用什么基组比较合适呢？我计算的是NaCl晶体掺Rh的g因子。看过您的帖子谈谈 ...

你发的地方不对，要发到简单量化答疑专贴里

作者
Author: wudazhuang509 时间: 2017-1-10 16:20

sobereva 发表于 2017-1-10 12:14
你发的地方不对，要发到简单量化答疑专贴里

不好意思哦，谢谢提醒

作者
Author: yjcmwgk 时间: 2017-1-13 13:13
http://muchong.com/bbs/viewthread.php?tid=1566146
老李，这是我以前在emuch上弄得。要不想办法合并进来？

作者
Author: liyuanhe211 时间: 2017-1-14 04:12

yjcmwgk 发表于 2017-1-13 13:13
http://muchong.com/bbs/viewthread.php?tid=1566146
老李，这是我以前在emuch上弄得。要不想办法合并进来 ...

我初学的时候着实从那个贴子里学了不少东西，刚开始写这个帖子的时候本来想照着那贴修改错误，但最后觉得一是太庞大、二是确实有些不太合适的地方，就直接重写了。

(, 下载次数 Times of downloads: 315)

我有个还想写内容的列表，那个帖子提到、这里没有的很多内容（我了解的部分）在逐渐写，过一阵稍闲些再更新。

作者
Author: yjcmwgk 时间: 2017-1-14 09:02
本帖最后由 yjcmwgk 于 2017-1-14 09:19 编辑

liyuanhe211 发表于 2017-1-14 04:12
我初学的时候着实从那个贴子里学了不少东西，刚开始写这个帖子的时候本来想照着那贴修改错误，但最后觉得 ...

嗯嗯！太有意义了。我也想加入。因为以前在emuch上弄得太早了（那时候我是硕士生涯即将结束，博士生涯尚未开始的时候），有些错误，很多不妥，如鲠在喉。而且当时弄这东西的人，以我为主，却也有十几个参与者一起完成的，有些乱，有些回答质量参差不齐。
所以，我和你一起弄新的FAQ行么？
又及，突发奇想：老李，要不然，你也跟广州墨灵格公司和上海绎模公司打个招呼，转个正，挂到他们那儿去？首发肯定是科音，以“原创非首发”的方式挂到墨灵格那儿去？他们那边儿搞Gaussian代理的人叫肖高铿。因为买版权和买书（Exploring那本书），我跟他联系过几次，感觉挺nice的（很可能是错觉，毕竟卖货的人对客户肯定都nice）。
这样做的目的是：早年在emuch上发的GaussianFAQ，被上海源资未经同意就直接挂网了。其实这样做没啥，我本来就说了，这东西完全可以不经过我同意而直接转载，随意散布。然鹅，当我联系源资，想更正或新增内容的时候，人家根本不理我。中国公司不守规矩，让我很尴尬呃……
所以，预感到你的FAQ也有可能被如此待遇，还不如直接拿到主动权。免费给他们挂网，但是要给你开个账号，给予修改内容的权限。

作者
Author: sobereva 时间: 2017-1-14 16:41

yjcmwgk 发表于 2017-1-14 09:02
嗯嗯！太有意义了。我也想加入。因为以前在emuch上弄得太早了（那时候我是硕士生涯即将结束，博士生涯尚 ...

我倒觉得这多虑了，肖主要感兴趣的是生物制药领域，而且曾经也有科研背景，不至于干这样的事。
上海那家我就不评论了

最不爽的是早年发emuch上的文，后来想改、想扩充，却提示说什么发帖后超过一定时间，不让改了。

作者
Author: liyuanhe211 时间: 2017-1-14 18:31

yjcmwgk 发表于 2017-1-14 09:02
嗯嗯！太有意义了。我也想加入。因为以前在emuch上弄得太早了（那时候我是硕士生涯即将结束，博士生涯尚 ...

噫，这么一说突然觉得很宏大 (⊙o⊙)，得想想怎么组织这件事，也搞个话题认领？就是怕又会乱了→_→以及怕后面就没空管了

作者
Author: sobereva 时间: 2017-1-15 08:58

liyuanhe211 发表于 2017-1-14 18:31
噫，这么一说突然觉得很宏大 (⊙o⊙)，得想想怎么组织这件事，也搞个话题认领？就是怕又会乱了→_→以及 ...

我个人觉得还是liyuanhe弄，他人提出意见和补充的形式比较好，人一多，风格和管理容易乱套，反倒还不如一个人弄效率高和效果好。如果等未来什么时候liyuanhe没精力继续负责FAQ维护，我或者其他人再接手。

作者
Author: liyuanhe211 时间: 2017-1-15 19:57
本帖最后由 liyuanhe211 于 2017-1-15 21:00 编辑

更新：

修改标题名，使新手更愿意看→_→
增加“提问后”
增加目录
解释2070错误后到哪里看报错的问题
windows G09W 32bit 资源限制
增加GaussView显示键级异常的问题
小错误若干

作者
Author: yjcmwgk 时间: 2017-1-15 21:23

sobereva 发表于 2017-1-14 16:41
我倒觉得这多虑了，肖主要感兴趣的是生物制药领域，而且曾经也有科研背景，不至于干这样的事。
上海那 ...

借老李这儿说个题外话：今天跟肖高铿联系了下。他也还么有Gaussian 16卖呢，再等几周才成。

作者
Author: sobereva 时间: 2017-1-16 07:05

yjcmwgk 发表于 2017-1-15 21:23
借老李这儿说个题外话：今天跟肖高铿联系了下。他也还么有Gaussian 16卖呢，再等几周才成。

正打算春节后就找他买一套G16

作者
Author: liyuanhe211 时间: 2017-1-16 18:19

yjcmwgk 发表于 2017-1-15 21:23
借老李这儿说个题外话：今天跟肖高铿联系了下。他也还么有Gaussian 16卖呢，再等几周才成。

噫，不知道北大得多久才能有

作者
Author: xujian 时间: 2017-1-19 16:37
本帖最后由 xujian 于 2017-1-22 07:38 编辑

李老师，我用CBS-QB3优化过渡态，关键词为opt=(ts,calcfc) freq CBS-QB3 出现了下图情况，L913，请问应该如何解决
(, 下载次数 Times of downloads: 142)

作者
Author: sobereva 时间: 2017-1-21 13:55

xujian 发表于 2017-1-19 16:37
李老师，我用CBS-QB3优化过渡态，关键词为opt=(to,calcfc) freq CBS-QB3 出现了下图情况，L913，请问应该如 ...

这根本就不是优化激发态的方法
自己用合适的级别优化激发态
当前是激发态定位完了算单点才出的问题

作者
Author: xujian 时间: 2017-1-22 07:40

sobereva 发表于 2017-1-21 13:55
这根本就不是优化激发态的方法
自己用合适的级别优化激发态
当前是激发态定位完了算单点才出的问题

老师，关键词打错了，我就是单纯算过渡态，用的是b3lyp优化好的结构，再用CBS-QB3算自由能

作者
Author: sobereva 时间: 2017-1-22 08:45

xujian 发表于 2017-1-22 07:40
老师，关键词打错了，我就是单纯算过渡态，用的是b3lyp优化好的结构，再用CBS-QB3算自由能

你这就是CBS-QB3当中CCSD(T)那一步迭代求解CCSD方程不收敛。用#P计算，如果发现CCSD迭代求解有收敛趋势，是因为默认的步数上限不够才导致报错，关键词里额外写上CCSD(maxcyc=100)把步数上限设大

PS: 问零碎问题请到简单量化答疑专贴里问，此帖只是展示FAQ目的

作者
Author: xujian 时间: 2017-1-22 21:29

sobereva 发表于 2017-1-22 08:45
你这就是CBS-QB3当中CCSD(T)那一步迭代求解CCSD方程不收敛。用#P计算，如果发现CCSD迭代求解有收敛趋势 ...

谢谢sob，提前祝您新年快乐！

作者
Author: liyuanhe211 时间: 2017-2-6 07:48
Update：
内坐标问题增加自定义 Linear bend 的解决方法。
增加 Linear angle in Bend --> opt=cartesian。
解释 GaussView 中显示的 IRC 在零点附近出现异常的低能点，

作者
Author: liyuanhe211 时间: 2017-2-17 16:56
更新：
增加Error: illegal instruction, illegal opcode
增加常被误以为是报错的语句
         · Warning!!: The largest alpha MO coefficient is...
         · Warning!!: The smallest alpha delta epsilon is...
         · Error on total polarization charges
         · This type of calculation cannot be archived.
         · DIIS: error

作者
Author: sobereva 时间: 2017-2-28 16:21
Internal input file was deleted! 可以加上

作者
Author: liyuanhe211 时间: 2017-3-1 20:59
增加：
Internal input file was deleted!
L301，End of file reading basis center.
L301，EOF while reading ECP pointer card.

作者
Author: liyuanhe211 时间: 2017-3-2 14:33
增加关于耗时估计的问题。

作者
Author: fly_ylx 时间: 2017-3-10 20:26
非常好，对问题很有帮助

作者
Author: 菜鸟一枚 时间: 2017-3-11 14:57
请问L308的错误是什么原因？ Error termination via Lnk1e in /share/apps/gaussian/g09D01/g09/l308.exe at Sun Mar 12 04:13:09 2017.应该怎么解决？谢谢老师

作者
Author: liyuanhe211 时间: 2017-3-11 17:02

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

根本没看帖子

(, 下载次数 Times of downloads: 147)

作者
Author: liyuanhe211 时间: 2017-3-11 20:57
增加对柔性扫描找过渡态的说明

作者
Author: 冰释之川 时间: 2017-3-12 09:00

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

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

作者
Author: 冰释之川 时间: 2017-3-26 12:02
补充一个报错，有些方法无法和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

作者
Author: liyuanhe211 时间: 2017-3-26 14:45

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

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

作者
Author: 冰释之川 时间: 2017-3-26 15:07

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

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

作者
Author: liyuanhe211 时间: 2017-3-26 18:13
更新：增加：
         · 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.

作者
Author: 菜鸟一枚 时间: 2017-3-27 16:25
老师好，我想问一下L308错误是什么原因？应该怎么解决呢？ Error termination via Lnk1e in /share/apps/gaussian/g09D01/g09/l308.exe

作者
Author: liyuanhe211 时间: 2017-3-27 16:52

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

先把帖子前半部分读了。

作者
Author: 冰释之川 时间: 2017-4-12 11:27
本帖最后由冰释之川于 2017-4-12 11:32 编辑

关于“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的姊妹泛函，可以保存算激发态密度】

作者
Author: liyuanhe211 时间: 2017-4-12 23:44

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

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

作者
Author: sobereva 时间: 2017-4-13 01:00

冰释之川发表于 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一起用。

作者
Author: 冰释之川 时间: 2017-4-13 09:04
本帖最后由冰释之川于 2017-4-13 09:08 编辑

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.

作者
Author: clvn 时间: 2017-4-16 11:15
老师，我看了您最后部分的通过做柔性扫描找过渡态，我有几个问题，还想请老师解决一下：
1.我用了您文中的opt=(calcfc,ts,noeigentest) freq这条命令，如您文中所述，最后的结构报错L9999,不收敛，但我不明白的是您在文中写到：“描时最左侧一点时会遇到了L9999 的 opt 未收敛错误“，而我做完opt=(calcfc,ts,noeigentest) freq只有一个结构，还是不收敛的，所以，我更加不明白老师文章里写的”此时打开输出文件可见路径上的最高点已经找到了“
2.您在文中说，近4步优化就找到了正确的过渡态，过渡态 C-H 键键长 1.675 A。我不是很懂您的4步是怎么看出来的
以上是我的问题，还请老师百忙之中，解答一下我的困惑，谢谢！

作者
Author: liyuanhe211 时间: 2017-4-16 11:54

clvn 发表于 2017-4-16 11:15
老师，我看了您最后部分的通过做柔性扫描找过渡态，我有几个问题，还想请老师解决一下：
1.我用了您文中的 ...

柔性扫描根本不是opt=TS，看教程
另外最好看前半部分该如何提问。

作者
Author: clvn 时间: 2017-4-16 14:37

liyuanhe211 发表于 2017-4-16 11:54
柔性扫描根本不是opt=TS，看教程
另外最好看前半部分该如何提问。

老师，您好，通过您的文章，我知道柔性扫描，不是opt=ts,我现在就是通过柔性扫描找到了势垒点，然后通过opt=(calcfc,ts,noeigentest) freq去找过渡态，
但是我仔细地阅读了您的文章，还是不太懂您的“4步优化就找到了正确的过渡态，过渡态 C-H 键键长 1.675 A”以及不用太在意L9999报错说结构不收敛，此时打开输出文件可见路径上的最高点已经找到了”
还请老师指点一下，谢谢！

作者
Author: liyuanhe211 时间: 2017-4-16 17:16

clvn 发表于 2017-4-16 14:37
老师，您好，通过您的文章，我知道柔性扫描，不是opt=ts,我现在就是通过柔性扫描找到了势垒点，然后通过o ...

文中说不必理会L9999不是opt=ts时不必理会，而是在opt=modredundant阶段时、若已找到最高点，则不必考虑。

如果还不明白，压缩后上传你的柔性扫描及opt=ts两个输出文件。

作者
Author: clvn 时间: 2017-4-18 16:57

liyuanhe211 发表于 2017-4-16 17:16
文中说不必理会L9999不是opt=ts时不必理会，而是在opt=modredundant阶段时、若已找到最高点，则不必考虑 ...

老师，我计算的是C6H14=C2H5+C4H9,通过做柔性扫描，发现，在2A左右是势垒，但是势垒附近变化平缓，然后我按照您文章后的操作，进行过渡态搜索，结果就报错了，L9999的错误，之前试过Sob老师博文里的消多于虚频的方法，结果不是很理想，还请老师百忙之中指点一下，谢谢！

作者
Author: liyuanhe211 时间: 2017-4-18 17:27

clvn 发表于 2017-4-18 16:57
老师，我计算的是C6H14=C2H5+C4H9,通过做柔性扫描，发现，在2A左右是势垒，但是势垒附近变化平缓，然后我 ...

你没有提供柔性扫描的文件，不知做的是否恰当，但2A这个距离经验上是比较合适的。这种过程的势垒确实一般不大，也常没有过渡态
另外既然在2A左右为扫描曲线的最高点，你为什么要将过渡态初猜中的键长定为3.5A？
L9999错误和“消虚频”没关系
你的优化振荡了，用本帖解决L9999振荡的方法解决。

作者
Author: clvn 时间: 2017-4-18 17:38

liyuanhe211 发表于 2017-4-18 17:27
你没有提供柔性扫描的文件，不知做的是否恰当，但2A这个距离经验上是比较合适的。这种过程的势垒确实一般 ...

抱歉，老师，是我自己表述不当，我的意思是键长增加了2A的位置，到达3.5A埃才出现势垒高度
老师觉得这个是合理的吗
附上的表格是自己整理的数据，以及柔性扫描文件
震荡和不收敛我都出现了问题，不过我试过了那些方法，似乎对我的这个不太管用……请老师指点一下，谢谢！

作者
Author: liyuanhe211 时间: 2017-4-18 18:14

clvn 发表于 2017-4-18 17:38
抱歉，老师，是我自己表述不当，我的意思是键长增加了2A的位置，到达3.5A埃才出现势垒高度
老师觉得这个 ...

你的CCSD(T) “单点扫描曲线”上确实有最大值，但在你找过渡态的 wb97xd/cc-pvtz 水平上却没有最大值。一般的TS都寻找的是电子能势能面上的鞍点，故而你在右侧的、wb97xd/cc-pvtz 势能面上做 opt=TS 是显然找不到过渡态的。右侧图应该下的结论是 wb97xd/cc-pvtz 水平下这个过程很可能没有过渡态。

(, 下载次数 Times of downloads: 125)

这个例子倒是提醒了DFT和CCSD(T)的区别。看扫描曲线，或许 CCSD(T) 水平下应该能找到过渡态，但CCSD(T)没有解析梯度，过于耗时，如果算的动的话或许可以尝试用CCSD扫描一下看有没有极大值点，如果有的话，用CCSD找过渡态还算可行、但也不简单。双杂化泛函诸如DSD-PBEP86也可以尝试。

这个过程的过渡态即使找到，也肯定是势垒很低的过渡态，其实意义不大

作者
Author: clvn 时间: 2017-4-18 18:35

liyuanhe211 发表于 2017-4-18 18:14
你的CCSD(T) “单点扫描曲线”上确实有最大值，但在你找过渡态的 wb97xd/cc-pvtz 水平上却没有最大值。一 ...

听了老师的讲解，我还有两个疑惑：
1.我的这类断裂C-C键或者C-H键，用我自己的方法找出的势垒，势垒点的结构能够代表过渡态的结构吗？和搜索过渡态的相比，会有很大差别吗?
2.之前尝试过CCSD(T),但是代价昂贵，所以放弃了，您说的双杂化泛函诸如DSD-PBEP86，我们没有试过，现在就是按照我们之前的扫描，得到的每一帧做freq，虚频多，而且我发现我扫描其他反应物断裂C-C键，也有同样的情况，但我做C-H键断裂，就没怎么出现这样的情况，这是C-C断裂的共性吗？而且我试过Sob老师关于消虚频的文章里方法，但我的虚频还是有，请问我应该怎么办呢？

作者
Author: liyuanhe211 时间: 2017-4-18 19:48

clvn 发表于 2017-4-18 18:35
听了老师的讲解，我还有两个疑惑：
1.我的这类断裂C-C键或者C-H键，用我自己的方法找出的势垒，势垒点的 ...

1. “你自己的方法”如果指柔性扫描，然后计算曲线上各几何结构的高精度单点能，根本就不能代表反应的势垒；可能会出现很大区别。
2. 我没动懂你到底想问什么。这和“消虚频”没关系。

作者
Author: tacoiong 时间: 2017-4-20 07:45
对于L123，Max corrector steps exceded
可以用recorrect=never, 跳过校正步吗?
因为我用了LQA后曲线不平滑...

话说我最近算出来的过渡态做IRC时都出现L123...

作者
Author: clvn 时间: 2017-4-20 17:35

liyuanhe211 发表于 2017-4-18 19:48
1. “你自己的方法”如果指柔性扫描，然后计算曲线上各几何结构的高精度单点能，根本就不能代表反应的势 ...

谢谢老师指点，我的疑惑就是断键反应很多都没有明显的势垒，即按我自己的拙见就是没有过渡态，那请问做搜索过渡态有何意义呢？
还请老师解答一下，谢谢！

作者
Author: laomalao 时间: 2017-4-21 17:44
学习一下

作者
Author: EdRhodes 时间: 2017-4-22 18:26
各位老师们你们好，在这里我想请教一个关于轨迹运算的问题。之前使用Gaussian默认的ADMP条件admp=(maxpoints=50,stepsize=1000)在不改变基组等条件的基础下可以顺利结束运算。如今想改变为admp=(maxpoints=120,stepsize=10000)，计算总是出错。还请各位老师指点一下。

输入条件：# b3lyp/6-31g(d,p) geom=(connectivity,crowd) admp=(maxpoints=120,stepsize=10000)，以下上传结束前的最后一页输出信息。

作者
Author: liyuanhe211 时间: 2017-4-22 18:31

EdRhodes 发表于 2017-4-22 18:26
各位老师们你们好，在这里我想请教一个关于轨迹运算的问题。之前使用Gaussian默认的ADMP条件admp=(maxpoint ...

你看了帖子了么

作者
Author: EdRhodes 时间: 2017-4-22 18:42

liyuanhe211 发表于 2017-4-22 18:31
你看了帖子了么

看了关于L502的一些，尝试修改，还是会在相同的用时后报错。

作者
Author: liyuanhe211 时间: 2017-4-22 19:51
本帖最后由 liyuanhe211 于 2017-4-26 18:48 编辑

EdRhodes 发表于 2017-4-22 18:42
看了关于L502的一些，尝试修改，还是会在相同的用时后报错。

真看了帖子就不会这样提问

作者
Author: clvn 时间: 2017-4-24 15:14
老师，您讲的针对L9999的错误的改正，似乎不太好用啊，请问您还有其他方法吗？

作者
Author: liyuanhe211 时间: 2017-4-26 18:49

clvn 发表于 2017-4-24 15:14
老师，您讲的针对L9999的错误的改正，似乎不太好用啊，请问您还有其他方法吗？

你这次上传的文件并无错误。之前已经说过了在当前级别下很可能并无过渡态，看之前的回复。

作者
Author: LL5312 时间: 2017-5-31 16:09
这个错误是什么原因呢在网上查了是加上%rwf=
这个是应该如何修改？

作者
Author: liyuanhe211 时间: 2017-5-31 20:45

LL5312 发表于 2017-5-31 16:09
这个错误是什么原因呢在网上查了是加上%rwf=
这个是应该如何修改？

更新了

作者
Author: gchen001 时间: 2017-6-2 13:29
求一个高湿达不留滴灵异的ubuntu版和高湿微哦的ubuntu版。谢谢各位大神！

作者
Author: liyuanhe211 时间: 2017-6-15 19:21
更新Distortion-Interaction模型的计算方法
避免误解、修正《如何通过柔性扫描找过渡态》-->《如何借助柔性扫描产生较好的过渡态初猜》，并在开头做出提醒
增加Atom too close
引用 http://sobereva.com/82 解释“如何在量化计算中设定原子/片段的电荷”

作者
Author: chenboeli 时间: 2017-6-19 20:27
谢谢谢

作者
Author: 不懂计算 时间: 2017-7-13 10:00
老师好，原文说“自出错的输出文件中提取的、出错前一步结构的原子坐标”，如何找到log文件中优化最后一步的原子坐标呢？

作者
Author: liyuanhe211 时间: 2017-7-13 10:43

不懂计算发表于 2017-7-13 10:00
老师好，原文说“自出错的输出文件中提取的、出错前一步结构的原子坐标”，如何找到log文件中优化最后一步 ...

GaussView打开输出文件，随便另存为一个gjf输入文件，用文本编辑器打开，把坐标部分拷出来。

作者
Author: liyuanhe211 时间: 2017-7-17 03:02
更新：

增加报错：
         · L1 ，ntrex1
         · L914，XXXXXX words are not enough for AIAXAO.
         · L508，Convergence failure
         · L123，GetHes: LRWHes > LHess!
         · L301，R6DS8: Unable to choose the S8 parameter
         · L1002，XXXXXX words are not enough for AIAXAO.
         · L801，Excessive mixing of core and valence orbitals.
         · Missing or bad data: Alpha Orbital Energies

增加：

关于Activation-Strain（Distortion-Interaction）的计算问题
关于充分试错的建议

作者
Author: steven 时间: 2017-7-19 18:54
很受用.感谢您的辛苦整理. 菜鸟不懂的还很多特别看了Activation Strain (Distortion-Interaction) Energy后心里一直在想: 还有这种操作...

作者
Author: Faye 时间: 2017-7-24 17:27
赞！都是干货，谢谢老师！

作者
Author: 夏洛特 时间: 2017-10-24 22:26
楼主，你好
我计算优化了一个多酸结构，计算后出现 NtrErr Called from FileIO.向上寻找出错原因发现L301中各中心没有基函数。在输入文件中没有读取chk文件的命令，所以也不存在读取chk文件的问题。我将输入和输出文件以及我认为的出错位置的截图都上传在帖子上，希望楼主指点，谢谢！

作者
Author: liyuanhe211 时间: 2017-10-25 01:46

夏洛特发表于 2017-10-24 22:26
楼主，你好
我计算优化了一个多酸结构，计算后出现 NtrErr Called from FileIO.向上寻找出错原因发现L301 ...

段落间空行太多，看: 详解Gaussian中混合基组、自定义基组和赝势基组的输入 http://sobereva.com/60

(, 下载次数 Times of downloads: 176)

作者
Author: 夏洛特 时间: 2017-10-25 11:11
谢谢楼主，解决了！

作者
Author: 夏洛特 时间: 2017-10-30 14:58
老师您好：
我最近计算一个多酸结构，最后输出文件在L9999出错，我在您前边的帖子中没有找到类似的解决方法，所以来求助您一下，我截了输出文件最后的一部分，因为输出文件太大无法上传所以就没有上传，请您不吝赐教，谢谢。

作者
Author: liyuanhe211 时间: 2017-10-30 16:43

夏洛特发表于 2017-10-30 14:58
老师您好：
我最近计算一个多酸结构，最后输出文件在L9999出错，我在您前边的帖子中没有找到类似的解决方 ...

说“我在您前边的帖子中没有找到类似的解决方法”是想找的什么？怎么找的？

作者
Author: 夏洛特 时间: 2017-10-30 21:37

liyuanhe211 发表于 2017-10-30 16:43
说“我在您前边的帖子中没有找到类似的解决方法”是想找的什么？怎么找的？

我复制图片中选中的一行在帖子中搜索，没有找到类似的表述，是这样的。

作者
Author: meatball1982 时间: 2017-11-8 15:34
本帖最后由 meatball1982 于 2017-11-9 08:45 编辑

不好意思，发错地儿了，已经修改。请删。

作者
Author: liyuanhe211 时间: 2017-11-8 17:45

meatball1982 发表于 2017-11-8 15:34
关于compliance对角线的输出。
根据：
http://sobereva.com/364

本贴与之无关，我个人也不熟，请另外发帖

作者
Author: 锁塘柳烟池 时间: 2017-12-15 18:03
谢谢。原来g09 d01算gd3bj的频率有bug，学习了

作者
Author: mmjjmmjj 时间: 2017-12-18 10:15
请问：高斯计算中*.cub和*.cube有何区别。想用Molekel打开*.cube，但是发现用Multiwfn生成的是*.cub，而这个Molekel打不开。谢谢！

作者
Author: sobereva 时间: 2017-12-18 15:56

mmjjmmjj 发表于 2017-12-18 10:15
请问：高斯计算中*.cub和*.cube有何区别。想用Molekel打开*.cube，但是发现用Multiwfn生成的是*.cub，而这 ...

一码事。改后缀就完了

作者
Author: flour_wheat 时间: 2018-1-8 12:09
这里面的distortion energy其实挺重要的。很多情况下优化的反应物的结构和初始结构是不同的。

作者
Author: flour_wheat 时间: 2018-1-8 12:15

sobereva 发表于 2017-1-14 16:41
我倒觉得这多虑了，肖主要感兴趣的是生物制药领域，而且曾经也有科研背景，不至于干这样的事。
上海那 ...

我的D01就是肖代理的。当时以为要付全款。后来肖帮着看了一下，别的学院已经买过B01，所以最后是升级到D01。
一高兴用省下的前直接买了Linux的GV509。买了以后发现没啥用。因为工作站一直是INIT 3。

作者
Author: doublelbegins 时间: 2018-1-13 17:32
老师您好，请问一下在使用混合基组进行Cu(NH3)4I2的优化和频率计算时出现了一下错误需要怎么解决呢？
输入文件如下
%chk=D:\g09w\Scratch\input\mechanism study\2.chk

#p opt freq b3lyp/gen scrf=(solvent=methanol,smd) geom=connectivity temperature=373.15

title

0,2

分子坐标信息

N H
6-31G(d)
****
Cu I
LANL2DZ
****
已经尝试过opt=cartesian 出现了同样的问题
麻烦老师了！

作者
Author: 赵云跳槽 时间: 2018-3-17 09:45
本帖最后由赵云跳槽于 2018-3-17 09:47 编辑

群上的一个问题
Q：
Two-electron integral symmetry is turned off.
R6DC6: No C6 coefficient available for IA= 96
Error termination via Lnk1e in /home-yw/Soft/Gaussian.D01.Linda/g09/l301.exe at Fri Mar 16 19:34:10 2018
C6系数是什么？

Sob：
C6是DFT-D校正里用的参数
DFT-D色散校正的使用
http://sobereva.com/210

要么当前体系有的元素不在DFT-D支持范围内，要么关键词写得明显不对

作者
Author: liyuanhe211 时间: 2018-3-17 14:06

赵云跳槽发表于 2018-3-17 09:45
群上的一个问题
Q：
Two-electron integral symmetry is turned off.

已测试并增加相关报错。谢谢！

作者
Author: 见见 时间: 2018-3-28 07:40
错误显示该操作超出范围。这是怎么回事啊？看了新手报错的博文，没有看到和我一样的错误。我尝试用水分子代替进行计算，输出正常。

作者
Author: anqingpingle 时间: 2018-4-8 21:58
输入文件这样：
chk=C:\2-ch-ga gaussian spe\2ch-ga-h.chk
# M062X/def2TZVP geom=connectivity output=wfn

显示图片中的信息后就很久很久不动了，百度后没明白什么问题，请赐教！

作者
Author: sobereva 时间: 2018-4-9 02:38

anqingpingle 发表于 2018-4-8 21:58
输入文件这样：
chk=C:\2-ch-ga gaussian spe\2ch-ga-h.chk
# M062X/def2TZVP geom=connectivity output ...

靠垃圾百度不可能解决问题，搞科研的一定要用Google

如果正在跑着，继续等
对较大任务强烈建议用#P监控状态

作者
Author: sobereva 时间: 2018-4-9 02:40

见见发表于 2018-3-28 07:40
错误显示该操作超出范围。这是怎么回事啊？看了新手报错的博文，没有看到和我一样的错误。我尝试用水分子代 ...

没完整输出文件没法判断。
优化和振动分析用6-311G**足矣，没特殊情况加弥散是瞎耽误时间
仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387

作者
Author: 见见 时间: 2018-4-9 14:52
1、请问除了改基组还有什么方法吗？我的是有些结构会遇到这种问题，有些是在同一基组下能够正常输出。2、我记得之前有上传过完整的输出文件，为啥今天过来看又没了

，现重新上传一份。

作者
Author: 哈哈哈 时间: 2018-5-8 10:42

sobereva 发表于 2018-4-9 02:38
靠垃圾百度不可能解决问题，搞科研的一定要用Google

如果正在跑着，继续等

#p 监控状态，看SCF收敛信息----这个能看到什么？我通常就直接看最后结果看是scf不收敛还是什么。如果scf不收敛再根据scf不收敛的解决方法一个一个地试。有几次我都看到用#p来看具体的信息，我想问那个能看到什么？是可以根据这个信息有依据地来修改吗？

作者
Author: Shine剪水 时间: 2018-5-8 14:19

哈哈哈发表于 2018-5-8 10:42
#p 监控状态，看SCF收敛信息----这个能看到什么？我通常就直接看最后结果看是scf不收敛还是什么。如果scf ...

#p 会打印更详细的信息

作者
Author: 哈哈哈 时间: 2018-5-8 14:29

Shine剪水发表于 2018-5-8 14:19
#p 会打印更详细的信息

我看的手册上是有每一链接在开始和结束与计算机有关的信息，和SCF计算收敛的信息。但具体的作用是什么？麻烦举一个例，或者麻烦给我说说哪能看到相关介绍好吗？

作者
Author: Shine剪水 时间: 2018-5-8 14:40

哈哈哈发表于 2018-5-8 14:29
我看的手册上是有每一链接在开始和结束与计算机有关的信息，和SCF计算收敛的信息。但具体的作用是什么？ ...

自己交一次对比一下不就好了

作者
Author: sobereva 时间: 2018-5-9 00:37

哈哈哈发表于 2018-5-8 10:42
#p 监控状态，看SCF收敛信息----这个能看到什么？我通常就直接看最后结果看是scf不收敛还是什么。如果scf ...

用#P时每一次迭代的能量变化、密度矩阵变化等信息都会输出，帮助你了解当前收敛趋势是什么样的

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