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

sobereva · 发表于 Post on 2017-4-25 10:52:47

huangtf 发表于 2017-4-24 14:58
大家好！最近我在学习Molcas8.0(免费版)/doc/manual里面的手册，准备进行CASPT2/ANO-RCC-VTZP基态优化NUN的 ...

请单独发帖，分类选择molpro/molcas，这样才会有更多人看到

498746012 · 发表于 Post on 2017-4-25 14:54:16

sobereva 发表于 2017-4-25 10:51
约束成哪个子群，修改结构的时候就会按照哪个点群来修改结构，便于准确、方便地构建高对称体系结构

谢谢sob老师

clvn · 发表于 Post on 2017-4-26 10:04:59

老师，再麻烦您一下，我最近在学习rrkm理论，公式上要求计算过渡态和反应物的态密度，我在gaussian中使用了pop=full，kisthelp却无法识别？请教一下您，我是命令用错了吗？怎样才是正确的呢？

adia · 发表于 Post on 2017-4-26 11:48:07

老师您好！我用高斯算了一下两个甲醛分子在三重态下反应的过渡态，最后结果就是其中一个甲醛上的H转移到了另一个甲醛的氧上，生成了HC=O和CH2(OH)这么两个片段，我开始认为这两个片段都是近似于碳上带单电子的自由基，随后又做了一下这两个片段生成乙醇醛的过渡态，最后能垒很小，像是两个自由基之间的结合。但是当我做了一下类似的反应，甲醛和甲醇的反应，甲醇碳上的氢转移到了甲醛氧上，生成了两个CH2(OH)，这两个片段再生成乙二醇的能垒非常高，显然不是两自由的结合，所以我就觉得开始认为两个都是自由基的想法不正确。于是我就想问，对于两个正常分子之间进行了氢转移，怎么判断生成的片段那些是自由基，那些是那种不是正常分子的非自由基片段呢？个人理论基础很薄弱，期待老师指点，谢谢

xiaojingling016 · 发表于 Post on 2017-4-26 17:44:14

高斯计算出现l9999错误，能量震荡不收敛，加了关键词opt=calcall，或opt=calcfc，任务一提交，还没几秒钟，就不能算了，不知道咋回事？求助各位大侠

xujian · 发表于 Post on 2017-4-26 18:26:54

FAQ:http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4829

xiaojingling016 · 发表于 Post on 2017-4-26 18:46:53

我算的是小体系，就是看的这个帖子，才加了opt=calcall，或opt=calcfc，但是就是算不了，不知道是不是高斯版本低的问题还是怎样

zjxitcc · 发表于 Post on 2017-4-26 19:52:44

得说详细点儿，交代体系类型，大小，基函数数目，log结尾报错信息（不要只说9999错误···），加了opt=calcfc之后的报错信息；不然让别人去猜答案是比较低效的

xujian · 发表于 Post on 2017-4-26 20:09:17

xiaojingling016 发表于 2017-4-26 18:46
我算的是小体系，就是看的这个帖子，才加了opt=calcall，或opt=calcfc，但是就是算不了，不知道是不是高斯 ...

可以把输出文件上传看看

xiaojingling016 · 发表于 Post on 2017-4-26 20:54:37

Gaussian 09:  EM64L-G09RevA.02 11-Jun-2009
            26-Apr-2017
******************************************
%nprocshared=16
Will use up to 16 processors via shared memory.
%mem=16GB
%chk=s1.chk
----------------------------------------------------------------------
#p opt=calcfc freq td=(nstates=6,root=1) B3LYP/6-311++G(d,p) scrf=(ief
pcm,solvent=toluene)
----------------------------------------------------------------------
No analytic 2nd derivatives for this method.
Error termination via Lnk1e in /data2/G09/g09/l1.exe at Wed Apr 26 20:48:33 2017.
Job cpu time:  0 days  0 hours  0 minutes 14.7 seconds.
File lengths (MBytes):  RWF=    17 Int=    0 D2E=    0 Chk=    2 Scr=    2
这是输出文件最后几行

paradise · 发表于 Post on 2017-4-26 23:14:07

问老师一个问题，阅读了您的博客：在Origin中绘制能量折线图的方法，有几个问题：
1.如果我在线的旁边添加其他说明性的文字，这个应该怎么添加？
2.基线的长度怎么才能和折线重合？我试了半天，总是不能很好的重合。
谢谢老师~

498746012 · 发表于 Post on 2017-4-27 13:55:23

sob老师您好！能不能打扰您一下，我一直有这样一个想法不知道对不对，高斯计算给出的电子能以及自由能都是单个分子的能量，不涉及系综的概念，系综是针对大量分子，而分子模拟给出的结合自由能，是针对大量分子，涉及系综概念，这样想法对不对？

clvn · 发表于 Post on 2017-4-27 22:50:36

请教老师，为什么我在计算irc路径的时候，forward和reverse路径都是键长减小啊？

paitel · 发表于 Post on 2017-4-28 01:21:25

还有Counterpoise指什么？如何用Counterpoise改进结果？

冰释之川 · 发表于 Post on 2017-4-28 07:51:01

（1）参见《Gaussian中非内置的理论方法和泛函的用法》http://sobereva.com/344
（2）参见《大体系弱相互作用计算的解决之道》http://sobereva.com/214

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

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

高斯计算出现l9999错误，用opt=calcall，或opt=calcfc，一提交就出错

如何在Gaussian中使用 MP2.5，MP2.X 和 DSD-BLYP-D3