公社群聊Q&A整理：2016.12.08 02 ~ 2016.12.10 09 Concate - 计算化学公社

2016.12.08 01:46:54
Q:
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Eigenvalue 6 is 0.00D+00 should be greater than 0.000001 Eigenvector:
A2 A1 R4 R2 R1
1 -0.66065 0.66065 -0.30938 0.17215 0.04157
R3
1 0.00000
NTrRot= -1 NTRed= 6 NAtoms= 4 NSkip= 0 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in C:\G09W\l103.exe at Wed Dec 07 17:38:02 2016. 请问计算甲醛裂解成氢气和一氧化碳的过渡态，l103报错啥意思？
Job cpu time: 0 days 0 hours 0 minutes 3.0 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1
A:
opt=cartesian
----------------------------------------------------
2016.12.08 01:46:54
Q:
[图片][图片]老师，刚看到你们的讨论想问一个问题，乙烯与氟化氢的加成，从能量图上看已经有两个挺明显的极小点，但是氟与氢没有明显成键，而且F-C间的键长只变化了0.2，这能判断是找对了过渡态吗？
[图片]振动模式，氟原子的振动很小，所以我算了下IRC，但是仍然没看到成键
没懂老师，怎么判断是找对了过渡态呢？现在频率为1839.29，这个频率处于好几百到上千的范围，所以我可以说找对了过渡态吗？老师
[图片]老师这是振动模式，氟几乎没动，我是不是要拿最后一帧去做优化看碳氟能不能成键，成键才能判断是找对了过渡态/
这是个加成反应氢和氟都要加上去，老师
A:
明显这是合理的过渡态
Q:
[图片]老师您这句话什么意思没太懂，因为碳氟原子之间振动模式不明显，甚至没有，所以我做了最后一帧的优化，难道现在还确定不了找对了过度态吗？
A:
怀疑过渡态找没找对走个IRC不就完了
光看振动频率也足够有80%~90%的判断可信性
Q:
老师我也走了irc,没看到碳氟原子之间成键，振动模式，氟原子也几乎没动
A:
这个就是过渡态结构
[图片]
Q:
[图片]
irc跑的是这样的
A:
两边都跑啊
Q:
这就是两边都跑的结果啊
A:
老生常谈
[图片]
----------------------------------------------------
2016.12.08 01:46:54
Q:
gassianD版本之后对IRC改了什么，老是出现Maximum number of corrector steps exceded
算法改成了什么？
09A也改了？
09A感觉没这么容易出现问题啊
A:
没改
09相对于03改了
IRC里写LQA
09都是HPC算法
Q:
走了啊
A:
没说你
Q:
[图片]
哦，我只是换D版本后遇到过这个情况
----------------------------------------------------
2016.12.08 01:46:54
Q:
那老师我的那个呢
A:
没时间看，太罗嗦
----------------------------------------------------
2016.12.08 01:46:54
Q:
加了opt=calcall还是L9999错误
A:
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278
自行看收敛曲线
calcall不是万能的
----------------------------------------------------
2016.12.08 02:08:50
Q:
请问，使用oniom方法找到的过渡态，计算能垒时，一定要做零点能校正吗？
A:
算不动频率就用电子能量求差
算得动频率就用自由能求差
Q:
嗯，算得动频率，谢谢老师。
我的情况是，如果不加零点能校正，过渡态的能量比反应物和产物都高一点（7kcal/mol），很合理。但是一加上零点能校正，过渡态的能量比反应物和产物都低，反而变得不合理了。
A:
过渡态若比反应物低就不是真正意义的过渡态，只是理论上的过渡态
Q:
就是说，零点能校正是必须的？
能垒是指零点能校正后算出来的能差吗？
A:
电子能垒较低时必须考虑
能垒可以是电子能垒也可以是自由能垒
你说的那是U(0)垒，没什么用
Q:
您的意思是可以不考虑？
A:
考虑啊
Q:
请问U(0)垒具体指的是什么？反应速率常数中的能垒是过渡态与反应物Sum of electronic and zero-point Energies之差吗
A:
[图片]
Q:
请问老师TST计算反应速率中的势垒是 [图片]之差还是[图片]之差？
不是绝对能量啊
A:
用自由能，怎么算这里说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327
因为图里是基于配分函数形式的TST，热力学形式的要用自由能
仔细看表格里的注释
基于过渡态理论计算反应速率常数的Excel表格
http://sobereva.com/310
要用自由能
要用自由能
要用自由能
----------------------------------------------------
2016.12.08 03:52:56
Q:
AIM可以做弱相互作用能量分解吗
A:
跟能量分解是两码事
----------------------------------------------------
2016.12.08 05:45:46
Q:
老师用M06-2X对分子进行优化计算单点能，然后如何计算零点能校正和色散力校正
A:
做freq
DFT-D色散校正的使用
http://sobereva.com/210
----------------------------------------------------
2016.12.08 05:59:57
Q:
请问计算出了tddft后怎么画光吸收光谱？
A:
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD和VCD光谱图
http://sobereva.com/224
----------------------------------------------------
2016.12.08 06:17:42
Q:
老师我对几何优化和Frep分的不是很清楚
# opt b3lyp/6-311++g(d,p) scrf=(iefpcm,solvent=generic,read) pop=(nbo,savenbo) 这是我之前做的优化，当然这里的B3lyp之后会换成M06-2X。麻烦老师帮我解释一下我这个输入文件是 “结构优化”还是frep，亦或者是两者一起优化了？
A:
freq关键词都没写怎么可能做freq
scf done
Q:
老师怎么做frep？
可以举个简单的例子吗
A:
RTFM
----------------------------------------------------
2016.12.08 06:24:20
Q:
老师有什么免费的软件可以搜索全局最小
谢谢
A:
使用molclus程序做团簇构型搜索和分子构象搜索
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=577
----------------------------------------------------
2016.12.08 06:46:57
Q:
我想问一下，如果算的吸收光谱的振子强度是0.0007是不是就证明此跃迁不能进行呢？谢谢老师
A:
极弱
----------------------------------------------------
2016.12.08 06:47:07
Q:
请问: 四聚体分子有两种排列形式，类似于口字型排列。请问MT方法算得的分子体积可否作为四聚体聚集紧密程度的一个指标。我算得的体积二者相差0.5%，可以说前者聚集程度稍逊于后者吗？谢谢！
A:
姑且可以
----------------------------------------------------
2016.12.08 07:29:00
Q:
各位老师计算的体系是orca下用blyp-D3/Def2-SVP算的，要用做NBO分析时要用到高斯是不是还得用blyp-D3/Def2-SVP.
A:
没这个要求
NBO对级别很不敏感
----------------------------------------------------
2016.12.08 07:40:51
Q:
老师，我在优化RS泛函的参数w的时候，区间选择为0.05-0.5，发现J2值随着w的减小而一直成线性减小，体系是一个有MLCT特性的金属配合物，请问这种结果合理吗，谢谢
A:
实际是怎么样就怎么样
----------------------------------------------------
2016.12.08 17:24:41
Q:
请问各个元素的范德华半径在哪里可以找到
A:
简谈原子半径
http://sobereva.com/255
----------------------------------------------------
2016.12.08 17:41:06
Q:
我想问一下 HOMO和LUMO可以用来分析分子间的弱相互作用吗？
OA：
分子间没有发生化学反应，仅仅有物理作用
A:
HOMO和LUMO用来分析分子间的弱相互作用是邪道
Q:
老师，我想计算溶剂分子和溶质分子之间的相互作用？应该通过计算什么数据来表征？
各位老师，我是想计算某种溶剂，比如二甲亚砜和甲苯之间的相互作用情况。他们之间并没有发生化学反应，只是存在单纯的物理作用。
老师，之前在Gaussian09中算过 “一个二甲亚砜和一个甲苯做优化”然后考察了单点能
老师，通过看文献，觉得M06-2X再加上色散力校正比较适合我的研究
OA:
对的，我就是通过计算单点能，然后再计算结合能
Q:
老师，你觉得我这个体系中会同时存在静电作用和非静电作用吗？
之前算的结合能，没BSSE校正是-1.10599
基组的话，如何选择呢？
A: 简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272
----------------------------------------------------
2016.12.08 18:28:23
Q:
各位老师今年gromacs培训具体时间是啥时候啊？
A:
预计1月18日，待定
----------------------------------------------------
2016.12.08 18:40:59
Q:
请问gaussian 可以计算荧光强度和荧光发射光谱吗？用G view能够看荧光光谱吗
我用[图片]
第一激发态跃迁回到基态的发射光谱就是荧光光谱？
[图片]
这个第一激发态第二激发态是选择最上面scale by后免的那个数字吗
请问怎么选择
OA:
好像gaussian view自动画出来的就是第一激发态的
Q:
[图片]
OA:
他这个下面有激发能量和振子强度
Q:
我对比了一下是第一激发态后面的数字
[图片]
OA:
所以应该是直接就是画的第一激发态的
Q:
对的
我错了
看来还是得分开算
[图片]对了，顺便问一下，这个excited state 1: singlet 表示的就是单线态第一激发态吗
[图片]我想问一下，发生内转换和振动弛豫后到达第一激发态，在此态上还会发生其他跃迁形式到达基态吗？谢谢啦老师@米豆豆
A:
Gery Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314
发荧光啊，当然会
----------------------------------------------------
2016.12.08 18:48:26
Q:
[图片]各位大神，请教一个问题，为啥我用multiwfn找不到价basin啊
说错了，找不到两单体间吸引子
[图片]各位大神，请教一个问题，为啥我用multiwfn找不到两单体间的吸引子啊老师，帮帮忙看下吧
A:
对什么函数做的盆分析？
----------------------------------------------------
2016.12.08 18:52:46
Q:
各位老师，我计算的是30到60个原子的有机体系，适合用什么基组计算达到比较好的精度呢？
A:
看机子
一般M06-2X-D3/6-311+G**考虑BSSE校正算的单点足够用了。还想更好用ma-TZVP
----------------------------------------------------
2016.12.08 18:56:00
Q:
老师，做势能面的话该怎么设置参数呢？
[图片]
A:
做势能面二维扫描，两个自旋态分别做
----------------------------------------------------
2016.12.08 19:04:48
Q:
[图片]请问这个是用什么关键词我开始以为是pop=npa，但结果不对
A:
就是pop=NPA
最稳定结构自己找
----------------------------------------------------
2016.12.08 19:05:04
Q:
刚才我在看这篇博文http://sobereva.com/336，里面提到弱相互作用分色散主导和静电主导，这怎么区分啊？
A:
这里说了哪些是色散哪些是静电主导
乱谈DFT-D
http://sobereva.com/83
科音基础班里再会更详细讲
----------------------------------------------------
2016.12.08 19:06:26
Q:
老师好，我算的体系涉及到氢原子的迁移，最好选用什么方法和基组？
A:
字号不要超过12
方法用M062X，基组看此文
谈谈量子化学中基组的选择
http://sobereva.com/336
----------------------------------------------------
2016.12.08 19:08:14
Q:
请问老师，输出电荷是在这吧[图片]
同样稳定结构输出结果[图片]和文献中结果[图片]差的很大
A:
要么用的结构不用要么用的理论方法不同
----------------------------------------------------
2016.12.08 19:09:31
Q:
老师，我使用了M06-2X-D3/6-311+G**+LANL2DZ考虑BSSE校正算的单点，还使用了ma-TZVP计算单点。有必要再使用B2PLYP/def2-TZVP+LANL2TZ(f)再计算单点吗？
A:
意义完全不明。ma-TZVP比def2-TZVP+LANL2TZ(f)都大，而普通泛函对基组质量要求比双杂化要低
----------------------------------------------------
2016.12.08 19:09:56
Q:
老师咱们在北京什么时间有教学班？
A:
什么教学班？所有培训班的预告都在www.keinsci.com上面
----------------------------------------------------
2016.12.08 19:10:29
Q:
请问老师，我在公社里面没有找到liunx系统下的g09软件包
老师您有链接吗
A:
liyuanhe前头已经给你回复了，勿再问
----------------------------------------------------
2016.12.08 19:14:39
Q:
老师，我是对wfn波函数做的分析
A:
你对轨道波函数做的盆分析？分析这有什么用？
Q:
我是想做分子生成密度差积分值
A:
那你对波函数做什么盆分析
要对密度差做盆分析
语言描述要严谨
Q:
哦哦，不好意思啊，多谢老师指正
A:
先看密度差等值面图，从原理上弄清楚哪里该有极小点和极大点出现，本来就没有的地方自然没有吸引子
我特别烦问题描述不清给回答者造成歧义和理解困难的情况
本来就很忙，一天回复问题无数，还得猜对方到底想问什么
----------------------------------------------------
2016.12.08 19:20:19
Q:
[图片] 【潜水】低调的华丽分享文件 11:20:19
"优干公示名单.docx" 下载
这是什么东西？
A:
删了
----------------------------------------------------
2016.12.08 19:22:14
Q:
老师您好，我的阅读权限不够[图片]
A:
看置顶的新人必读直到攒够积分获得阅读权限
要么就自己google或者用网盘搜索工具从别的地方找
----------------------------------------------------
2016.12.08 19:23:01
Q:
请问老师，GView可以画出Spin density吗，请问如何画
A:
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353
----------------------------------------------------
2016.12.08 19:28:22
Q:
大家好，谁有MS软件操作方面的资料？
发给我一份吧，谢谢！
A:
先看自带的help里的教程
其余的google一搜一把
群专用邮箱里也有的是
----------------------------------------------------
2016.12.08 19:29:25
Q:
老师我刚看了你发给我的帖子，[图片]
隐式模型不会考虑氢键作用吗？我的研究体系中应该是存在氢键作用的
A:
当然考虑不了溶质和溶剂的氢键作用
考虑必须用显式模型或者杂化模型
靠molclus产生溶剂-溶质的二聚体或多聚体构型通过量化静态地研究相互作用是可以的
Q:
我刚才也想用molclus多画出几个溶剂，然后做RDG或者BCP什么的，这样合理么社长
A:
合理
Q:
老师，还是采用结合能来表征作用力大小？
A:
除能量以外的表征方式还很多
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252
Q:
是不是构建一个溶质分子和多个溶剂分子，然后考察多聚体的相互作用?
需要做你说的RDG吗
A:
仔细看完博文再说
----------------------------------------------------
2016.12.08 19:41:05
Q:
[图片]老师您好，连接Na原子越多的Al原子上的电荷是比较大的，但是在我做的Al5Na2结构中，连接两个Na原子的Al原子上的电荷反而是最小的，是-0.05，连接一个Na原子上最底部的Al上的电荷是-0.98，我做出了elf图，这个怎么解释啊？
A:
ELF图和电子密度分布没有直接关系
如果你用的原子电荷计算方法合理，则绘制变形密度等值面图比ELF平面图明显更能解释原因
----------------------------------------------------
2016.12.08 19:45:40
Q:
老师，我觉得这个图做出了好丑啊。有没有更好的建议啊？RDG vs Sign(Lamda2)
[图片]
感觉有点复杂啊
A:
用VMD敲一个命令加载.tcl文件就出来，不比chemcraft方便多了
----------------------------------------------------
2016.12.08 20:25:21
Q:
老师，使用了VMD得到RDG vs Sign(Lamda2)等值面的图了，如何可以显示原子名称啊？
还有这个标度如何得到？
[图片]
A:
在VMD中显示原子序号的方法
http://sobereva.com/197
[图片]
之后自行ps一下刻度条
----------------------------------------------------
2016.12.08 20:57:57
Q:
老师，请问分子动力学培训，什么时候开始
A:
1月中旬，看官网通知
Q:
多少天，
A:
3
----------------------------------------------------
2016.12.08 21:44:37
Q:
能和其他个培训放到一块培训就好了
外地的没有在北京的方便
A:
没有什么好处
1：会同时参加两种内容不同的培训的人不多。即便是量化初级班+基础班，放到一起也没意义，因为参加过初级班之后需要回去先充分消化
2：同时参加两个班，能请假那么长时间的人极少
3：连续讲课一个多礼拜吃不消，备课也来不及很充分，还会导致很多其它事务诸如打理论坛和QQ群无法正常进行
4：一次领不了那么多空白发票，不够给两次培训的人开的
----------------------------------------------------
2016.12.08 21:50:06
Q:
老师用molclus 做构象搜索的时候输入./molclus 时出现权限不够怎么办
A:
chmod改
Q:
是在molclus这个文件夹下修改？
A:
google chmod用法
Q:
chmod 750 -R *
[图片]
老师现在出现了这个错误
这个应该怎么解决
A:
settings.ini格式不对
改的时候仔细些
Q:
[图片]
A:
别删掉其它的
Q:
老师这是我的settings.ini 执行./molclus
A:
改哪个就改哪个，其它别动
Q:
[图片] 老师出现这种无法获取的字样有关系吗
A:
明显没设对，上面提示了你指定的是g09目录而不是g09可执行文件
Q:
[图片]这是我的g09 环境变量的设置
[图片]
A:
和这个无关，设的是settings.ini里g09可执行文件目录
和这个无关，设的是settings.ini里g09可执行文件
要么直接设"g09"就完了
Q:
[图片]
老师这样可以了吗
A:
自行查看.out文件看怎么回事
明显不正常
----------------------------------------------------
2016.12.08 22:12:26
Q:
量化基础班啥时候再有呢？没有赶上12月的[表情]
A:
12月是初级班
明年上半年有第三届基础班
Q:
初级班啥时候还有？
A:
明年上半年
[图片]http://www.keinsci.com/
----------------------------------------------------
2016.12.08 22:12:26
Q:
其实听老师课的感觉，就是觉得滔滔不觉，哪怕溜号一会，可能都错过一些知识。不过，听起来真是过瘾。就是下午吃饱有点困，以后中午不吃了。
A:
喝红牛
----------------------------------------------------
2016.12.08 22:26:05
Q:
[图片]各位老师请问我在做静电势的时候这个颜色怎么改啊？我看了sob老师的帖子http://sobereva.com/253 最后那个设置颜色的位置没有找到请各位老师指点一下谢谢啦
A:
再仔细看，每一步都写得很清楚
Q:
[图片]老师这段儿我看了但是没有找到图的上面有所以没找到
A:
再往上看啊
----------------------------------------------------
2016.12.08 22:29:26
Q:
相比于第二届的基础班，第三届会加入那些内容呢？
A:
各方面都有增加
第二届记得大概1400页，第三届2000页
----------------------------------------------------
2016.12.08 22:37:24
Q:
大家好，我想问一下，我用vmd做的图如何保持它的分辨率？[图片]
A:
VMD渲染图像问题看1.4节
用Multiwfn+VMD做RDG分析时的一些要点和常见问题
http://sobereva.com/291
----------------------------------------------------
2016.12.08 22:37:38
Q:
老师，直接用B3LYP/6-311G（d，p）优化过度态和计算能垒结果可信吧？
A:
可信
Q:
论坛上说至少计算能量的时候要加弥散
您的博客
A:
我从来没这么说过！！！
我专批乱加弥散的行为
Q:
http://sobereva.com/119
里面有一句话[表情]
A:
哪？
[图片]
Q:
如果体系是阴离子那句啊
A:
你哪说的你的体系是阴离子了？默认当然是中性
Q:
哦哦阴离子如果不加弥散算能垒可接受吗
A:
不可
----------------------------------------------------
2016.12.08 22:45:07
Q:
Another possibility is delayed fluorescence, the transition back to the first excited singlet level, leading to the emitting transition to the ground electronic state
我谷歌看看ISC，发现有一种情况下可以发出荧光的，但是荧光发射的时间被延迟了，原文说的是跃迁到三重激发态回到了第一激发单态上进而发射的荧光，[图片]那么是否此时的三重激发态要比第一激发单态高呢？谢谢老师
我没加载出来，所以发给你网站让你试试
最近网速奇慢
A:
T1比S1低，但由于这类体系反向系间窜越比较容易，故可以还跑到S1再发荧光
----------------------------------------------------
2016.12.08 22:52:15
Q:
各位老师，算频率的时候出现了一个-4.41的小虚频能不能忽略不计？
A:
一般不能忽略
Q:
老师应该怎么解决一下，我试了几种方法都没消掉~
A:
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278
----------------------------------------------------
2016.12.08 23:05:28
Q:
请问老师，我用Gamess做能量计算的结果里MP2能量是正值，而且在Gamess里通过整体减片段计算的相互作用能也是正值，这是什么原因造成的？
A:
没在相同级别下优化过结构，或者优化时用的级别不合理
或者GAMESS-US做HF时候没收敛到正确基态波函数（可以和高斯的比较）
或者关键词不对
----------------------------------------------------
2016.12.08 23:19:51
Q:
请教下您，如何获得笛卡尔坐标的力常数矩阵？谢谢！
A:
参考
基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq
http://sobereva.com/328
----------------------------------------------------
2016.12.08 23:21:45
Q:
Sob老师，我怎么看gamess里有没有收敛到正确基态波函数？
A:
把HF能量跟高斯的对比。高斯更容易收敛到基态波函数，和高斯不一样时大多数都是GAMESS-US错
Q:
谢谢Sob老师！我还有一个问题是，我把高斯轨道函数输入到Gamess中计算后，结果显示的是MP2与HF的能量差值怎么是正的呢？[图片]
A:
要比单点能，MP2单点能一定比HF负
相互作用能是另一码事
----------------------------------------------------
2016.12.08 23:22:22
Q:
老师，你好！请问这个多中心键级怎么算呢
[图片]
A:
迗使看Multiwfn手册4.9.2节的例子
----------------------------------------------------
2016.12.08 23:25:06
Q:
我这是Gaussian 09 A版？
A:
字号不要超过12
----------------------------------------------------
2016.12.08 23:31:48
Q:
请问老师，电子传输速率怎样计算？谢谢
@Sobereva
A:
看marcus理论相关知识
----------------------------------------------------
2016.12.09 00:02:37
Q:
[图片]＠Ｓｏｂ老师，我想问一下做ＭＦＤＤ，等值线图上面ＮＣＨ部分明显有电子密度差变化，但是做彩图通过设置色彩刻度上下限，但是ＮＣＨ部分就是不显示，怎么回事呢
[图片]
A:
MFDD是啥？
Q:
分子生成密度差
[图片]
A:
很烂的名词
如果是变形密度就说变形密度
这图8成又是核坐标不对应
Q:
两种类型的图为什么彩图不显示上面ＮＣＨ电子密度变化呢
A:
核坐标对上再说
色彩刻度范围太大时填色图也很难展现数值很小区域特征
Q:
ＮＣＨ在Ｚ轴，环在ｘｙ轴
可是色彩变化度太小环部分就区分不了
A:
没跟你说轴的是，而是求差时候波函数文件里头核坐标必须对应
Q:
其中的连个片段都是从二体复合物剥离出来的，也要从新调坐标吗
两个片段
A:
调坐标干嘛，调了就麻烦了
用nosymm
Q:
[图片]
每一个计算ｗｆｎ文件都有ｎｏｓｙｍｍ
A:
自行用文本编辑器打开wfn文件，仔细检查开头的坐标部分是否对应
另外，如果不是>=D.01版，必须加上density
怎么改输出格式博文里都说了
----------------------------------------------------
2016.12.09 00:39:14
Q:
但是感觉VMD的字体不好看
A:
嫌VMD的字体不美观就自行ps
----------------------------------------------------
2016.12.09 00:43:48
Q:
老师，bond order类型很多，含有氢键的选择哪种来分析呢
A:
弱相互作用不适合用键级来讨论
合适的方法都在这里
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252
----------------------------------------------------
2016.12.09 01:16:30
Q:
老师，您好。使用VMD做Tachyon渲染后，发现原子的label有发生位移，这个如何能够解决？
A:
自行ps
----------------------------------------------------
2016.12.09 01:17:07
Q:
[表情]
老师，请问promolecular RDG好还是直接使用RDG好啊？
A:
直接用不动RDG才用promolecular近似
----------------------------------------------------
2016.12.09 01:31:48
Q:
问一下， Actual scratch disk usage=564470925504 words 这个说明占了多少的硬盘啊？搞不懂怎么换算~~
A:
1word=8字节
----------------------------------------------------
2016.12.09 01:49:56
Q:
老师，有文献使用NBO analysis来研究弱相互作用，使用E2来分析作用强弱，这个分析的意义大么？
A:
毫无意义，而且误导
----------------------------------------------------
2016.12.09 01:59:05
Q:
用高斯算出来H+ + OH- =H2O的δ_rGm 和用查物化书标注生成吉布斯自由能得到的δ_rGm会差很大。高斯算出来是-330kcal/mol，查手册数据是-79.9kcal/mol
我用H3O+ + H2O---OH-=2H2O 高斯算出来的更离谱，
H3O+ E=-76.68 H2O--OH-是-152.3 2个水是2*-76.5
这样算起来都70多Hartree了。。。
问题出在哪儿了呢
A:
1 用热力学组合方法算，免得中途步骤搞错，而且结果还准
2 注意物质状态，是液态还是气态，相差蒸发焓，免得和实验没可比性
Q:
第二个坐标是1个水合H+和一个OH-和水的复合物。
第二个反应
得到2个水
都是PCM溶剂模型下算的
您说的热力学组合方法是什么？
A:
G4、CBS-QB3之类
而且，第二个左右两边带的电荷都不同
----------------------------------------------------
2016.12.09 02:05:50
Q:
听说引用过千可以评院士了吧
A:
院士哪这么好评
在我来看过万还差不多
过千只算是“好意思和人打招呼”的level而已
----------------------------------------------------
2016.12.09 02:12:21
Q:
请教个问题，我用orca的cosmo溶剂模型优化了分子，但是换成smd模型优化很难收敛，记得sob老师
A:
那你就用默认的IEFPCM呗，收敛只会比SMD容易，优化结构的话不会有什么差别
Q:
sob 老师既然是溶剂对构型影响下，我为了比较溶剂效应，我直接用真空中的构型在溶剂模型中算个单点可以不
A:
可
----------------------------------------------------
2016.12.09 02:18:20
Q:
[图片]
老师这是一个TS找过渡态的能量变化图
正常吗
A:
看最终结构
Q:
也就是说跟这里的能量是否降低或升高无关对吧
[图片]
这也是一个找过渡态的
这两个都找到了而且是对的
不应该TS优化的时候能量应该是降低的吗
A:
没这回事
看曲线没用，干嘛关心这个
能量降低也不代表就到找了TS
是升高还是降低取决于初猜结构
TS又不是极小点
Q:
我懂了
[图片]
老师出现这样的不收敛怎么办
尝试了博客上的方法
还是出错
A:
别一出错就问，再多试试，调整初猜和用博文里的做法
本来过渡态就远不像做几何优化那么简单
Q:
以上都是我试过之后才来问的
并不是一有问题就问他
A:
有的人为一个过渡态话一个礼拜，甚至一个月时间
反复这里问也白搭，再牛的人也不敢说给个反应，尝试3次以内就找到过渡态
反复发图还会造成刷屏，令管理员头疼
Q:
我第一个过渡态花了8个月
A:
找过渡态谈何易事
----------------------------------------------------
2016.12.09 02:28:39
Q:
老师，这种图可以直接使用Multiwfn做出来么？还是需要其他软件？
[图片]
A:
散点图能直接作出来。然后自行ps一下弄成这种彩色图
----------------------------------------------------
2016.12.09 02:36:42
Q:
老师，请问下最后不能正常结束，这是个计算小分子过渡态的，这是输入和输出[图片][图片]看了下l913是计算后-SCF能量和梯度，那我应该怎么修改呢？
A:
甭用G4
改用CBS-QB3之类
----------------------------------------------------
2016.12.09 02:40:47
Q:
老师，请问VMD能否使一部分原子cpk模式显示，其他的使用line显示？
A:
selected atoms里进行选择即可
----------------------------------------------------
2016.12.09 02:40:57
Q:
windows下有什么软件可以打开.agr类型的文件吗？
.agr 是二维数据文件，如下图示意
[图片]
A:
windows的直接打开那个的大抵没有。
一般都是删了前面那些行在origin里作图
----------------------------------------------------
2016.12.09 02:46:11
Q:
老师，用CBS-QB3也计算过了，显示的是错误结束不收敛，我把CalcFC改成了CalcALL,依旧是不收敛，这个过渡态结构是我用DFT B3LYP优化成功后用来计算的。
A:
那你就甭用热力学组合方法，省得没法控制细节
calcall和那个毛关系也没有，那个是CC/QCI迭代不收敛，又不是几何优化不收敛
Q:
老师，明白了，谢谢！
老师，在此情况下我应该怎么修改呢？在CBS-QB3情况下，迭代不收敛。
A:
直接用普通方式计算单点
Q:
好的，那我可以尝试下用这篇http://sobereva.com/61[图片]解决不收敛的问题吗额？
A:
CC/QCI迭代和SCF迭代不一码事
[图片]
----------------------------------------------------
2016.12.09 04:06:52
Q:
那我要看电子密度变化用第一种操作还是第二种操作好一些呢
两种图那个说明问题更好一些呢
A:
两种函数显然差老远
仔细看手册2.6节对函数的说明
比较密度的变化对密度的梯度模做差作甚
理清思路别乱算
----------------------------------------------------
2016.12.09 04:27:33
Q:
[图片]
请问sob老师写过这篇《力场瞎总结》吗？怎么在思想家公社门口找不到这篇
A:
力场瞎总结大抵是以前发在内部交流论坛中的
Q:
sob老师，这篇能发出来吗？想学习一下[表情]
----------------------------------------------------
2016.12.09 04:49:19
Q:
谢谢白菜鸟只是概念理解障碍症，在优化出现震荡问题中，怎么理解每步位移不是很大而受力很大这句话呢？
A:
甭纠结那个，直接把步长上限改成3~5再试
----------------------------------------------------
2016.12.09 05:12:34
Q:
老师：您好！在算分子体积的时候[图片]，为什么有些有这个密度，有些没有呢？[图片]
还有请问老师。下面的这个最小值是什么值呢？[图片]
A:
老版本没有对密度的估计，后来版本加上去的
就是分子表面上你研究的那个函数（往往是ESP或ALIE）的最小值
----------------------------------------------------
2016.12.09 05:13:56
Q:
[图片]请问sob老师，Hess2freq 程序中这样处理从原理上来说是不正确的吧？
A:
怎么不正确？
Q:
将3N 个频率按照绝对值从小到大排序，将最低的6 个视为整体运动模式，振动低频不一定比转动和平动频率高啊？
A:
怎么不一定？结构都经过优化了，这是一定的，除非根本就没优化
----------------------------------------------------
2016.12.09 08:16:28
Q:
高能了
[图片][图片]
这个是不是比潘姆巴现象更有意思
特殊情况没条件看原文。仅根据题目猜是把水塞到纳米管里形成窄列就倾向于变得特别整齐？感觉并没有什么特殊性。水中加个非极性分子，周围一圈的熵也都很小，水气界面上表面几层分子的光谱(拉曼？忘了)也很接近冰的光谱。
A:
如果是计算研究的话，这种研究感觉都老掉牙，做烂了
水在各种受限空间中的特征、行为的MD模拟已经不计其数了
----------------------------------------------------
2016.12.09 17:25:57
Q:
旋轨耦合只能通过cas来计算吗？如果应用其他计算软件是否需要使用cas算法？
[图片]
我来啦
A:
诸如ADF可以计算TDDFT下的旋轨耦合
----------------------------------------------------
2016.12.09 18:15:02
Q:
请教个各位老师个常识问题，分子前线轨道和静电势是分子的什么性质，是electronic characteristics
A:
静电势相关资料看
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=219
----------------------------------------------------
2016.12.09 18:26:30
Q:
老师好，请问计算弱相互作用的charge最好使用哪一种电荷呢？我看到有人使用Mulliken,有人使用NPA的。
A:
必须用能够较好重现分子表面静电势的，NPA、mulliken这种电荷烂到爆。看群共享里原子电荷计算方法的对比，写得极为详细
----------------------------------------------------
2016.12.09 18:50:35
Q:
EOMCCSD，39个原子的计算激发能耗时如何
A:
几乎算不动
----------------------------------------------------
2016.12.09 18:51:43
Q:
昨天听了一个报告，说纯水的模拟，TIPxP力场的准确度比DFT要高的多
凝聚相性质为标准
比如蒸发焓，沸点，密度之类
DFT算得水的密度偏低
MP2 模拟水的密度偏高
只有力场的恰好....恰好...
双杂的模拟未见报道
[图片]优化一个泛函，对水的密度优化
赶快做，一片JPCL诞生了，再用它算算水的蒸发焓一片，再跟别的比较benchmark又一片，这样三篇文章够一个博士毕业了
A:
关于前面讨论模拟水的，不专门给水优化参数的DFT、后HF，模拟结果一定比不上诸如用TIP4P-2005、TIP4P-FB这样专门拟合实验结果的水模型...
专门拟合套力场参数，模拟结果完爆DFT的AIMD太容易了
----------------------------------------------------
2016.12.09 23:34:27
Q:
老师，我计算好单点能，要在这个基础上用同等级别的方法和基组去计算别的，是不是就不用输入理论方法和基组了？
A:
当然得输入啊，要不然高斯哪知道用什么级别
----------------------------------------------------
2016.12.09 23:42:02
Q:
请问，用高斯怎么可以看原子的价态呢
A:
高斯没法算。
在multiwfn手册里搜LOBA，用LOBA方法可以算
----------------------------------------------------
2016.12.10 01:59:34
Q:
打开bashrc之后怎么设置，输入啥
A:
molclus根本不用设置环境变量
----------------------------------------------------
2016.12.10 02:20:10
Q:
sob老师，用td-b3lyp 进行bomd（gradientonly…）分子动力学模拟，不知道怎么回事就断了？
[图片]
A:
问高斯客服
或者换基组之类，反复尝试找对策
----------------------------------------------------
2016.12.10 02:54:47
Q:
还想请教一个问题，gaussion09和03可以混着用吗？着急算几个任务，但是03的结点不够用想用09算一部分，可以吗？结果有没有可比性？[表情]
A:
如果结果要横向对比，不要混用
----------------------------------------------------
2016.12.10 03:10:38
Q:
过几天初音的演唱会要开始了，大博士别忘了照照片
A:
洒家就不去了
一是太贵，二是没时间
Q:
[表情]本以为我们不去了，大博士肯定得去呢
A:
要是去了，怕是动力学班就得延期到春节后了
Q:
[表情]才一个晚上，至于吗？
A:
最低票价还680
Q:
[表情]原来大博士也是节俭主义者
A:
不是节俭，是合理主义者
上限480我能接受
680除非是水树或者miu's来
写repo也得花好多时间
----------------------------------------------------
2016.12.10 03:12:33
Q:
老师，有两个个问题想请教一下。（1）我用材料的垂直电离势和垂直电子亲和势表征材料与阳离子和阴离子的相互作用能力是否合理（考虑到吸附阳离子是与阳离子发生电子交换，材料供给电子；吸附阴离子是材料与阴离子发生电子交换，材料得电子），这样得到整体的性质，再对局部进行计算；（2）我对一个碳基模型掺杂杂原子后，发现其垂直电离势降低，垂直电子亲和势提高，表明改性后材料得电子和失去电子能力同时增强，不知道是否可能出现这种情况。
A:
直接计算与实际离子的相互作用能多好。IP/EA只是间接的描述
Q:
谢谢老师。确实更好，一来想补充数据，二来模型比较大，考虑计算每个位点需要太多计算量，于是选择了几个位点进行计算，想通过IP/EA间接描述材料整体的性质。但就是不知道是否会由于不能用于描述对离子的吸附强弱以及出现基础错误（对一个碳基模型掺杂杂原子后，发现其垂直电离势降低，垂直电子亲和势提高，表明改性后材料得电子和失去电子能力同时增强），导致被别人质疑。
A:
你也算算吸附氨气，和文献里比比，确认计算方式合理性。合理的话和咪唑作用的情况就是那样。而且别太关注gview是否显示成键，要看距离
----------------------------------------------------
2016.12.10 03:25:42
Q:
请问老师，Gaussian用户手册中说 ”CBS-4M 和CBS-QB3 只能用于第一和第二行原子”，这其中的第一行和第二行原子指的的是什么？
A:
一般意义的第二、三周期
----------------------------------------------------
2016.12.10 03:27:36
Q:
老师，关于分析氢键的成分（种类），有没有类似的帖子
A:
看相关资料：
书
The nature of the hydrogen bond : outline of a comprehensive hydrogen
文章
http://dx.doi.org/10.1080/00268976.2012.661478
dx.doi.org/10.1021/cr800346f
----------------------------------------------------
2016.12.10 03:38:03
Q:
动力学班大概时间确定了吗
A:
下周一发通知，1月18日
----------------------------------------------------
2016.12.10 03:59:35
Q:
[图片]请教sob老师及各位前辈，此错误如何避免？
A:
[图片]
----------------------------------------------------
2016.12.10 04:03:50
Q:
老师，我用命令pop=saveNLMO计算NBO电荷，但是NATURAL LOCALIZED MOLECULAR ORBITAL (NLMO) ANALYSIS:
Highest occupied NBOs are not at the beginning of the NBO list;
The NLMO program is not currently set up to handle this.
Illegal range for read: IUnit= 7 Address= -256 Length= -256 File length= 18830336.
Error termination in NtrErr:
NtrErr called from NtrIO.请问老师可以怎样处理
A:
就用pop=NPA就完了
那不叫NBO电荷
量子化学中的一些常见不良写法和用词
http://sobereva.com/298
----------------------------------------------------
2016.12.10 04:54:20
Q:
用gaussian做金属配位化合物的charge distribution计算，发现mulliken会显著低估金属的charge，和基于atomic polar tensors（APT）的分析虽有低估，但情况稍好一些。比如Fe（II）到Fe（III），mulliken的结果是0.7→1.4，而APT是1.3→2.6。请问Soberva社长，金属离子电荷被低估是一个普遍的状况吗？是否有相关文献报道其原因呢？在网上很难找到APT的相关资料，所以只好叨扰一下社长大人@Sobereva
但是在论文写作是，审稿人会要求计算结果与化学常识相符合，不然就需要解释其偏差的原因；若是有相关文献报道，引证一下也是必需的。可以具体解释一下或者推荐几篇文献吗？有劳啦@离庸
我明白，真实电荷应该是离域的是population或distribution，但是导师让我千万别跟审稿人干架[表情]
所以现在有些无所适从。。。
A:
GAPT电荷的原理和性质
http://sobereva.com/208
GAPT电荷很烂，甭用
想让电荷数值较大用NPA或者CHELPG
原子电荷跟价态根本没有任何直接对应关系，妄图让原子电荷与价态对应纯属胡来
这如同让真的东西非要和假的东西一致
算价态有专门的方法，去最新版multiwfn手册里搜LOBA
----------------------------------------------------
2016.12.10 05:06:44
Q:
用每个负极大值/负极大值之和的比值来比较（相当于标准化）？
----写漏了
A:
可以直接比较。
求和没有意义
----------------------------------------------------
2016.12.10 05:15:12
Q:
请问sob老师3号邮箱的密码改了吗？怎么登陆不上去呢？
A:
貌似被某个卑鄙小人改了密码，现在正在进行找回密码的申诉
----------------------------------------------------
2016.12.10 05:30:25
Q:
老师，也就是说ADCH或MK电荷貌似是一个很好的选择，对于涉及离子的弱相互作用？
A:
是
----------------------------------------------------
2016.12.10 05:34:38
Q:
[表情]好吧，祝愿sob老师早日成功
A:
本群发言无用“好吧”做句子开头，看群共享里群规
----------------------------------------------------
2016.12.10 05:34:56
Q:
谢谢社长
用不同方法计算结果：
APT：1.3→2.7
Mulliken：0.7→1.1
NPA（NBO）：0.8→1.4
ADCH：0.4→
现在看上去是APT最大，我再试一下您提到的LOBA
A:
大决不代表正确。AIM电荷甚至更大
博文里说了，GAPT电荷烂得很
----------------------------------------------------
2016.12.10 05:37:39
Q:
老师，我算的体系里有Fe,铁用lanl2dz，其他用6-31g**，铁有五重态和三重态两种，但对这两种自旋进行初猜。。优化，结果都是三重态的结果。不知道和选的基组有没有关系，还是结构或其他的原因
A:
怎么可能。设几重态就算几重态
Q:
我分碎片了，是指铁上设置是6，氧上是-1，但优化结果是铁上3，氧上零点几
不好意思，刚才没说清楚
A:
本来SCF之后就不可能正好和你的初猜一样
除非是做CDFT
Q:
老师，刚看了您的博文，CDFT在NWChem中能用，在高斯中能用吗
A:
不能
Q:
我也不是想设定几就是几，但现在铁上自旋感觉差太多
A:
当前是极为正常的
不如说你之前的错误观念根深蒂固
量化代表正义
----------------------------------------------------
2016.12.10 05:38:14
Q:
那两分子的原子电荷可以直接比较大小吗？
A:
可
Q:
我利用两个分子的表面静电势的正极值点比较是想判断哪个分子是更好的氢键给体，可以吗？这两个分子结构不同的。
A:
可以
----------------------------------------------------
2016.12.10 05:44:15
Q:
老师我想确定一下:(1)在计算某些位点与离子的相互作用能后计算垂直电离势和电子电子亲和势能表征整体的对阴阳离子的吸附能力(得失电子能力)是否会受到质疑?(2)对材料改性，发现与原材料相比垂直电离势降低(说明更易于失电子)，同时垂直电子亲和势增加(说明也更易于得电子)，这是否违背基本常识或原理。谢谢老师。
A:
1 没必要算这个。
2 不违背
----------------------------------------------------
2016.12.10 06:02:15
Q:
[图片]
这是啥？
A:
不要发乱七八糟的图
已记违规一次
已踢。之前就发过什么培训教育的链接
----------------------------------------------------
2016.12.10 06:22:03
Q:
极化终于做好了~
OA：
啊
看你的双电子积分在
全是微分递推
Q:
我更新了
OA:
更新到哪里了
GitHub?
Q:
可以算任意的角动量
没放
OA:
角动量?
Q:
而且没有依赖
极化的东西吗还是啥
通用的双电子积分程序
spdfgh都能算
OA:
咋搞
用的rys?
高科技啊
Q:
恩
OA:
rys
我靠
Q:
一天看文献，一天写
OA:
啊
Q:
除了速度慢，没别的毛病
OA:
加油加油
建议再用OS写一个，OS可以算到任意高角动量，只要内存够
Q:
恩恩
必须的
OA:
经验是，越难写的程序效率常常越高[表情]
话说，OS做了的话。。。HGP干脆就做了吧。。。prism也能做了。。。
Q:
恩恩！！
强烈同意
OA:
小白求问双电子积分一般在具体计算中起什么作用
没有双电子积分，就没有量化软件
Q:
盖房子砖头什么作用
G(0,0) = I00
if( Lab > 0 ) then
G(1,0) = C00*G(0,0)
end if
if( Lcd > 0 ) then
G(0,1)=C00p*G(0,0)
end if
if( Lab > 1 ) then
do i = 1,Lab-1
G(i+1,0) = i*B10*G(i-1,0) + C00*G(i,0)
end do
end if
if( Lcd > 1 ) then
do j = 1,Lcd-1
G(0,j+1) = j*B10p*G(0,j-1) + C00p*G(0,j)
end do
end if
if( Lab /= 0 .and. lcd /= 0) then
do i = 0,Lab-1
G(i+1,1) = (i+1)*B00*G(i,0)+C00p*G(i+1,0)
do j = 1,Lcd-1
G(i+1,j+1)=j*B10p*G(i+1,j-1)+(i+1)*B00*G(i,j)+C00p*G(i+1,j)
end do
end do
end if
Ir = 0.0d0
do i = 0,Ld
Temp2(Lc+i) = 0.0d0
do j = 0,Lb
Temp2(Lc+i) = Temp2(Lc+i) + Comb(Lb,j)*(Ra-Rb)**(Lb-j)*G(La+j,Lc+i)
end do
Ir = Ir + Comb(Ld,i)*(Rc-Rd)**(Ld-i)*Temp2(Lc+i)
end do
return
end subroutine Rys_Ixyz
!=======================================================================
薛定谔方程中，哈密顿算符有一项是双电子项
Rys的递推部分
OA:
啊
就这么一点东西
Q:
这个周末做做批处理
OA:
rys也有递推目测
Q:
就能快一些了
有
很简单
OS麻烦
如果我水一点用递推
目测很快
OA:
你可以试试OS
OS
HGP好在能提前收缩
它是OS的升级版
每次被大神们碾压
Rys和OS递推的区别是
Rys
[图片]
Rys的xyz分的很清楚
Q:
很抱歉再次打断提个问题你们说的RYs 和 OS 指的是算法吗？
OA:
OS的麻烦在不是那么分明
Q:
问题是大神用的通项嘛
OA:
算法
Q:
好像rys写起来还是好点的
推公式还是有点累的
双电子积分主要有MD,Rys,OS,HGP
OA:
都是积分方法嘛
Q:
Rys简单是因为有人做好了权重和根
话说为啥没人穷举一下
OA:
像你那样
就像boys函数
Q:
啊
OA:
有日本人做了个新算法叫ACE，公式恶心无比
你看看libint
Q:
因缺思厅
OA:
大神加油
Q:
化学里的boy 函数是不是类似于物理里面的wannier 函数
opt巨恶心
其实还有中国人做过，算法是直接积分球GTO，用到3j，6j符号，似乎也只停留在论文上
OA:
这个就没了解了
Boys函数就是个纯数学函数，erf，没啥物理意义
Q:
恩恩
所以我要是用rys递推公式
也是可以的吧
而且最后应该只要用到一个boy函数
rys用不到boys函数
啊那那个四个S轨道的积分
OA:
都不用
Q:
但是巨慢
OA:
好高端
我晕
算不动嘛会
角动量越低，Rys越慢
天
表示只是刚刚学Fortran写个程序练手
怎么感觉坑好大
Q:
不急，你才研一。
OA:
你也研一[表情]
话说今天导师说
他那里有用python写的量化程序
而且功能类似于高斯
Q:
都有什么功能
OA:
我当时问了他说和高斯差不多
CC之类的高精度也有
好像没写密度泛函
Q:
那叫差不多。。。。
OA:
老师还说Python用numpy模块
效率和c一样[表情]
又不是我说的
[表情]
Q:
怎么可能一样
OA:
比比就知道了
因为Python就是c写的[表情]
A:
真能扯
Q:
matlab还是fortran写的
我也没见比gfortran快
也不是我说的
[表情]
CPython是用C写的，python还有其他多种实现
高斯，目前没有任何一款量化软件的综合性能做得过他
[图片]
话说敢这么说也是够有自信的～
全世界没几个做量化软件的敢这么说的
可怕貌似撞枪口了
就算ORCA发展那么好，也不敢这么说～
orca好像要出4.0了？
一直造势，说4.0里有FOD分析什么的
A:
根据放出来的消息，比较有意义的也就FOD和PBEh-3c
看当前苗头，G16无望了
眼看就圣诞节了
Q:
G03 G09 按照6年间隔，应该出新版了
Gxx里应该有些能量分解才好
不然总得用别的软件
为什么能量分解对机子要求高
比如GAMESS，原子数一多。。。内存要求就很大
A:
不算高
比如Gaussian+Multiwfn那种简单能量分解和单点耗时基本一样
GAMESS-US那个在定域化过程消耗太大
如果你用PSI4做过SAPT0发现速度还是相当快的
所以主要都看程序、所用算法
像IQA能量分解巨慢
molpro做SAPT0也特慢
Q:
PSI是SAPT那帮人的“亲儿子”，自然快
不过最近那帮人发了个文章，SAPT(DFT)前端是ORCA
A:
最初搞SAPT那帮人开发的是他们自己的SAPT程序
PSI是Sherril那伙人搞的，还不支持DFT-SAPT
----------------------------------------------------
2016.12.10 06:55:57
Q:
FOD一般用来干嘛的
A:
可以看看grimme这篇最新的和FOD有关的
http://dx.doi.org/10.1002/chem.201604682
PS：引了Multiwfn
----------------------------------------------------
2016.12.10 07:16:03
Q:
[图片]老师这是硝基苯的氯代反应，报错L9999，现在结构已经是邻位取代的结构了，我是加calcall再优化好些？还是直接用loose好些呢？
这算震荡吗？感觉不太明显
A:
用最后的结构，加上calcfc再继续跑
这也不是硝基苯啊
算得动的话直接用calcall更省心
----------------------------------------------------
2016.12.10 07:20:13
Q:
casscf算出来的自旋应该非常可靠了吧..
A:
自旋具体指什么？
Q:
磁矩一类的..
磁性分子之类的..
含Fe的分子..
A:
这得有自旋密度
但一般CASSCF计算给不出自旋密度，只能给总的密度矩阵
Q:
是不给自旋，自旋密度得自己手动弄了..
不过依靠占据数来判断自旋..这个靠谱不..
比如说d轨道上单占据的电子数量
A:
够呛
还不如用UCCSD计算得到自旋密度，很准
Q:
嗯..
可以用UCCSD的密度来做自旋布局分析吗...
我一直以为只能用SCF来做
A:
可以。比如在高斯里用UCCSD的时候写上pop=NOAB就给你UCCSD级别的Mulliken自旋布居，可以进一步用Multiwfn分析
Q:
嗯.[表情]
那么多参考，多组态的情况.. 用uccsd可以吗.。
A:
虽然不算最理想，但比DFT还是可靠性更高的，而且又是黑箱方法，我认为足够了。何况自旋分布也不像能量那么对级别敏感
----------------------------------------------------
2016.12.10 07:52:33
Q:
AIM中BCP分析结果和分子间结合能有没可比性？@Sobereva
A:
看具体分析方式、体系特征
诸如体系只由氢键维持，V(BCP)/2和相互作用能有不错对应
Q:
只有范德华力呢？
这个体系是不是没有AIM分析的必要？老师，请指教。@Sobereva
A:
没必要
----------------------------------------------------
2016.12.10 07:58:38
Q:
[图片][图片]老师这个用了calcall,结果仍显示L9999报错，能量和平均受力图差别很小很小，最后一帧的构型有些小小的差别，这是不是不合理呀？还有请问老师，“算不动”是啥表现呀？
A:
不要写cartesian，笛卡尔坐标下优化通常没有冗余内坐标下容易
opt后头别有空格
继续跑下去有希望收敛，继续跑，可以一次性把maxcyc设大点比如150
----------------------------------------------------
2016.12.10 08:01:06
Q:
我没看到有这个cell文件吗，是我哪里出错了吗
A:
算不动就是机子慢，耗时太多
----------------------------------------------------

复制代码