简单量子化学问题答疑专帖（旧版，已关帖）

lufenghxu · 发表于 Post on 2017-2-11 15:26:16

sobereva 发表于 2017-2-11 14:25
从头到尾仔细看一遍log文件，如果没看出原因，说明gv版本太老或者有bug，也可尝试用chemcraft打开

老师说的不错，刚把电脑重启，又能打开了。上次晶格内复制分子出现死机。看来是bug啊。谢谢老师。

erya123 · 发表于 Post on 2017-2-12 17:17:45

猜想过渡态结构时，初猜结构和真实结果越接近，越容易找到。是不是gaussian也会自己改变输入结构（自己设定路径）来确定hessian矩阵中的负值（找到虚频）。和采取scan的方式寻找过渡态相比，gaussian是不是不再自己改变结构，而是依据设定键长（键角）的改变计算，看这条路径上是否可以找到虚频？

sobereva · 发表于 Post on 2017-2-12 20:36:35

erya123 发表于 2017-2-12 17:17
猜想过渡态结构时，初猜结构和真实结果越接近，越容易找到。是不是gaussian也会自己改变输入结构（自己设定 ...

意义不明。什么叫Gaussian也会自己改变输入结构。你输入文件里是什么结构初始结构就是什么结构，程序以这个结构为起点，按照过渡态搜索算法，一点一点逼近实际过渡态。算法见
过渡态、反应路径的计算方法及相关问题
http://sobereva.com/44

除非极特殊体系，否则scan根本不能用来找过渡态。如果要帮助找合适的初猜结构，那也是柔性扫描而不是刚性扫描。

lily19830130 · 发表于 Post on 2017-2-13 00:07:21

我在文献中看到说geometry optimization were perforemed using the M06 function.但是我在gaussview中没有看到相关设置，不知道这个M06 function怎么用。还有文献中一处不懂内容如下：In basis set system I, the Los alamos relativistic effective core potentials were used for the core electrons of Ni and a (541/541/311/1)basis set was employed for its valence electrons .The 6-31G9(d) basis sets were employed for other main-group elements. A better basis set system (BS-II) was used to evaluate the potential energy changes. In BS-II, a(311111/22111/411/11)basis set was employed for Ni with the CEPs of the Stuttgart-Dresden-Bonn group.以上内容中描述的这些541,311111等，这些基组是不是可以理解为针对不同原子设定不同基组吗？代表什么意思啊？本人刚接触计算对这些不胜理解。

sobereva · 发表于 Post on 2017-2-13 03:23:53

lily19830130 发表于 2017-2-13 00:07
我在文献中看到说geometry optimization were perforemed using the M06 function.但是我在gaussview中没有 ...

不是function而是functional，如果文章确实那么写的，说明作者是个菜鸟。
直接在关键词里写M06就完了，看g09手册DFT部分http://sobereva.com/g09

那种基组表示方法表示的是基组怎么收缩的

如果还搞不懂，把这些看了http://bbs.keinsci.com/forum.php?mod=viewthread&tid=275

zc000 · 发表于 Post on 2017-2-13 09:09:16

请问老师，在HgCl 结构优化结束后（输入文件见截图1），打开文件时，出现了warning（截图2）。这个什么原因？

ssp2493187 · 发表于 Post on 2017-2-13 09:18:28

本帖最后由 ssp2493187 于 2017-2-13 09:19 编辑

请教各位老师，我在小木虫和gaussFAQ上找l801错误的解决方法没找到，只是在高斯说明书上对l801解释是双电子积分交换的初始化，希望各位给些建议，谢谢(具体信息见图)

冰释之川 · 发表于 Post on 2017-2-13 09:24:17

Gaussian常见简单报错及解决方法（Ctrl+F 关键词检索，定位至相应段落）：

         · 2070 Error (Windows)
         · segmentation violation/segmentation fault (Linux)

         · L101，End of file in ZSymb.
         · L103，Error in internal coordinate system
         · L103，Linear angle in Bend
         · L103，FormBX had a problem
         · L123，Max corrector steps exceded
         · L123，GS2 Optimization Failure.
         · L502，Convergence failure
         · L502，Inaccurate quadrature in CalDSu.
      · L801，Excessive mixing of core and valence orbitals.
         · L1002，Inaccurate quadrature in CalDSu.
         · L9999，Optimization stopped.

参看置顶帖《Gaussian FAQ；新手求助报错时的注意事项（持续更新）》http://bbs.keinsci.com/forum.php ... 4829&extra=page%3D1

lily19830130 · 发表于 Post on 2017-2-13 10:56:20

sobereva 发表于 2017-2-13 03:23
不是function而是functional，如果文章确实那么写的，说明作者是个菜鸟。
直接在关键词里写M06就完了 ...

十分感谢，是functional，我写错了。看来真是菜鸟，有好多东西得学。谢谢老师!

sobereva · 发表于 Post on 2017-2-13 16:23:55

不用管它，不影响观看结果。是gview的bug，对于直线体系经常出现

登登 · 发表于 Post on 2017-2-13 16:39:27

老师，用Mnltiwfn分析福井函数和双描述的时候，您的说明书里分析了闭壳层的情况，对于N-1，就是将电荷数和自旋多重度0,1改成1,2算单点。那对于自旋多重度为2的情况应该怎么改呢,是改成1,3还是1,1。

sobereva · 发表于 Post on 2017-2-13 17:39:48

登登发表于 2017-2-13 16:39
老师，用Mnltiwfn分析福井函数和双描述的时候，您的说明书里分析了闭壳层的情况，对于N-1，就是将电荷数和 ...

1,1

登登 · 发表于 Post on 2017-2-13 18:53:16

本帖最后由登登于 2017-2-13 19:42 编辑

sobereva 发表于 2017-2-13 17:39
1,1

谢谢老师的回复，还有一个问题，我看了您之前的回复贴说，如果要比较不同分子间的亲电亲核活性，应该要将福井函数乘上一个全局软度之后才能比较。这个全局软度在Mnltiwfn里是怎么算出来的，我想比较加载不同官能团的分子间的活性位点强弱。是不是用CFF算出来的值不需要乘全局软度就直接可以比较了？

sobereva · 发表于 Post on 2017-2-13 23:25:29

登登发表于 2017-2-13 18:53
谢谢老师的回复，还有一个问题，我看了您之前的回复贴说，如果要比较不同分子间的亲电亲核活性，应该要将 ...

全局软度用不着Multiwfn算，你在http://bbs.keinsci.com/forum.php?mod=viewthread&tid=384这里的综述中随便一搜就找到公式

lily19830130 · 发表于 Post on 2017-2-14 12:22:51

sobereva 发表于 2017-2-13 03:23
不是function而是functional，如果文章确实那么写的，说明作者是个菜鸟。
直接在关键词里写M06就完了 ...

好像链接已经失效了。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）

优化出现的问题

请教l801错误

评分 Rate