|
|
本帖最后由 liyuanhe211 于 2015-11-23 10:24 编辑
由于原程序有Bug，本帖于20151106更新一次（下面内容，包括附件，已替换原来内容）：
全文附件（所有已整理聊天记录，下载时请确认本帖为我最后发布的聊天记录帖（不一定是日期范围最早/最晚的））：
（附件内容与下文一致）
2015.10.31 05:33:05
Q:
是不是计算溶液中的热力学性质时
目前主流的溶剂化模型就是这个SMD啊？
嗯，谢谢
我以为SMD模型计算的热力学性质更好呢
那一般的氨基酸在水溶液中的热力学计算，SMD更准确还是，以前的COSMORS、CPCM IEFPCM这些结果好呢？
我在想，会不会有这样的benchmark。。。
A:
@大漠孤狼 SMD一定会比CPCM好，但是可靠度<=COSMO-RS
SMD的极性部分和PCM本质一样，但是在非极性部分做了更多考量，溶解自由能算得比PCM结合普通方式的非极性部分处理更准。
----------------------------------------------------
2015.10.31 05:33:05
Q:
想问下在作NBO的轨道图时，我想显示出两个不同构象的steric能量密度差该怎么办呢？
A:
@flybird 没法显示这样的东西，NBO框架里没有能量密度
----------------------------------------------------
2015.10.31 05:33:05
Q:
[Image]
彭国伦F95教材画人头像的练习题，测试后的确可以得到人头像，但里面的updatescreen()子程序有两个语句理解不了
第一个语句：write( sSSStr, "('(',I3.3,'A1)')" ) ScreenWidth
经常见write(*,"()")这种格式，没见过把第一个星号改成其他符号的格式
以上是第一个问题，第二个问题：在textgrphlib.f90的内容中全篇搜sSSStr关键词，也只发现出现3次，均在updatescreen()子程序中。但似乎sSSStr的含义是一种输出格式，而不是变量
screenwidth本身是一个数字， write( sSSStr, "('(',I3.3,'A1)')" ) ScreenWidth这句话里面要求用I3.3的格式输出这个数字screenwidth，但后边又有个‘A1’表示要求以1个字符宽输出字符串，这个以‘A1’要求的字符串还是screenwidth这个整数么？ 程序输出结果中并没有给出screenwidth对应的数据。
第三个问题：textgrphlib.f90中的DestroyScreen这个子程序似乎也没有被调用过？ 我把它删除了，似乎“人头”还是出现了
A:
@廖荣宝 write后面第一项的sSSStr是指写入到字符串
----------------------------------------------------
2015.10.31 05:33:05
Q:
旁边组买的机器用 i7-4790K （4核8线程，4.0 GHz-->4.4 GHz）才￥2000 一颗，相比E5-2650 v2 （8核16线程，2.6 GHz-->3.4GHz）大概￥7000 一颗，除了多耗电之外，考虑到很多时候都需要同时算数个结构，是不是买好几台独立的i7-4790K更划算？
A:
@LiYuanhe 2650 v2你报得太贵了。如果用2640v3，8核2.6，才5200
而且不光是耗电，你攒两台机子，硬盘、主板、机箱电源等等都得额外多花钱。
----------------------------------------------------
2015.10.31 05:33:06
Q:
老师，我要是想把两种构型的光谱叠加在一起，因为可能溶液中会有两种构型，应该怎么做呢，需要看的能量是哪个呢？
A:
@Nicole origin里把两条曲线对应的列
按一定比例加和再作曲线图即可
----------------------------------------------------
2015.10.31 06:42:48
Q:
# Question Not Found Markdown
@Sobereva 谢谢！但是我看到文献里面有人生成啊
你的意思是说energy只是一个点吗？
A:
--------Markdown Markdown--------
先搞清楚定义
本来NBO里就没有这个东西，定义清楚了才能考虑怎么去算
有文献算不代表任何问题
NBO里基于NLMO的steric energy计算公式都是现成的，只是计算个能量而已。若在其基础上搞什么其它的量，必须有明确的计算公式或方法说明。不能随意弄出一个新词。
----------------------------------------------------
2015.10.31 09:55:06
Q:
用NBO计算出来的原子上的电荷数，有单位么？是 e 么？
请求帮助
@206 这个有文献参考么？我现在正在找，但是还没有找到，谢谢！
A:
怎么老想要文献呢
1+1=2这种事不需要文献证明
都是搞计算的人的常识知识
e都没必要写，我根本不写，除非编辑或者审稿人找碴
----------------------------------------------------
2015.10.31 10:48:00
Q:
请教一个问题，其再算双自由基的反应，利用uhf/6-31+G**可以得到想要的过渡态，然后再用这个结构再ub3lyp/6-31+G**下计算就得不到想要的过渡态了，请问这是怎么回事，dft方法再这里有什么问题呢？
A:
@R2D2 HF成份越高往往势垒越高、越容易出现对称破缺。尝试M06-2X，HF成份比B3LYP高得多
[Image]
----------------------------------------------------
2015.10.31 10:50:34
Q:
恩恩，sob老师，还有一个问题，我用uhf计算的时候用的opt,freq。然后用b3lyp在优化出来的结构继续opt freq,是不是opt这里出了问题，我可以直接freq么
A:
不能
----------------------------------------------------
2015.10.31 10:52:14
Q:
就是说opt和freq必须是同一计算水平才有意义的，是这个意思吧
A:
y
----------------------------------------------------
2015.10.31 13:59:48
Q:
[图片]这种势能面扫描需要自己找出对应态的能量？
中间有个cross
在扫描过程中怎么确定哪个激发态才是想要的激发态？
你是说这个？[图片]
如果基态是A'，是不是激发态也是算A’的电子态就行了？
涉及到基态和激发态跃迁的过程
sob，A'和分子光谱项的Σ-是不是相同意思@Sobereva
貌似文献上是只有A'的@华电ding @Sobereva
然而我想问的是一个Cs对称的构型，从A"到A'的势能面的跃迁是禁阻的吗？@Sobereva
哦，谢谢sob，就是只能是不可约表示相同的才能跃迁？@Sobereva
此处的x y z 表示坐标？
A:
一个基态一个激发态，作曲线
看看Hg的p基函数的布居数
恩
图中都这么表示了，就算A'就完了
不可约表示不同，必然禁阻
恩，跃迁偶极矩必为0，振子强度也因此为0。
如果两个电子态的不可约表示的直积不含x、y、z中的任何一个，则说明跃迁电偶极矩的模一定为0，因此振子强度为0，此时对称禁阻。
@杉演杰方 看特征标表。不理解的话可以看cotton的群论书
----------------------------------------------------
2015.10.31 14:07:14
Q:
请群里的老师们帮我看看，这个问题该怎么解决？
A:
装python
----------------------------------------------------
2015.10.31 14:21:10
Q:
@Sobereva sob老师，对于H原子，计算时对其加极化函数还是很必要而加弥散函数意义不大是吗？即对于含H体系，6-31G**还是明显优于6-31G*，而6-311++G**不如6-311+G**划算，对吗？
A:
看考察什么问题，比如氢键、氢转移之类的，氢最好还是加极化，多耗不了多少时间。其它情况氢不加一般也无所谓。
弥散函数是那样。
----------------------------------------------------
2015.10.31 14:26:12
Q:
sob老师，还有一个问题请教您，上午的时候我说用ub3lyp泛函计算双自由基是出错，您建议我换m06-2x，我的问题是我还怎么体现非限制性计算，使我算的双自由基处于单重态呢？
A:
@R2D2 就是这帖子里的做法
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
----------------------------------------------------
2015.10.31 14:26:22
Q:
@Sobereva 老师您好，请问做了pop=npa后，想把每个原子的电荷标注在原子上面，有哪个软件可以实现啊？
所以想问除了高斯还有什么软件可以做
A:
这点事自行ps就好了
----------------------------------------------------
2015.10.31 14:30:59
Q:
最近看Franck-Condon原理，跃迁偶极矩的计算中经过BO近似，轨旋分离近似后，得到三个项，其中包括一项空间重叠积分[图片]。[图片]对这一项不是特别清楚，想请问一下。
（假如基础表述有问题，请不吝赐教）
1、psi和psi'分别对应于基态轨道和电子跃迁到达的轨道，假设分别是π和π*轨道；
2、根据FC原理，mu_e为常数，可以提到积分号外面，得到mu_e*<psi'|psi>，积分部分即为π轨道和π*轨道的内积；
可是，
问题1、分子轨道不是正交的吗？上述内积岂不是等于0？那么还怎么据此讨论轨道重叠对光谱的影响呢？（比如π和π*重叠大，跃迁容易；CT激发轨道重叠小，跃迁困难）
可能是我有哪个地方理解错了，请sob老师指点~~
问题2、上述积分的数值如何导出呢？
谢谢~~
A:
mu_e提不到外面啊，实际上其形式就是坐标变量(x,y,z)，积分的时候得带着它
----------------------------------------------------
2015.10.31 14:32:37
Q:
@Sobereva 老师，IRC路径为什么要用质量权重坐标呢，为什么不用一般的笛卡尔坐标？二者差别大吗，还是用质量权重坐标是为了理论处理上的方便？请您指教
那也就是说考虑质量权重是必要的，IRC比MEP更合理吗？
哦，好的。另外通过IRC给出的是电子能量差，而通常我们报道能垒时要对电子能量加热力学矫正，这样怎么在文章中将二者统一起来呢？
A:
说得粗略一些，受同样的力，C原子走得快还是H走得快？显然是H
对于化学问题上是这样
你可以每个点都算热力学校正量加上去，当然这就比较耗时了
通常只是对于TS和极小点考虑热力学校正
----------------------------------------------------
2015.10.31 14:44:38
Q:
sob老师，我计算的结构本身就是非对称的，分割成片段的，不能保证连个片段的电子数相等，这个该怎么解决呢？
[图片]
我的意思是分成的两个片段如果按我的这个画法是一个8，一个16，，，
triplet?双自由基不是单重态么
[图片]
是这样么？sob老师？
A:
本来就不可能相等嘛
两个都设triplet。unpaired spin一个设alpha一个设beta
没关系啊
双自由基是整体啊
你先在设片段，两个片段都是3重态，自旋方向相反，整体就成单重态了。
是
----------------------------------------------------
2015.10.31 14:48:13
Q:
[图片]似乎有点明白了，这里讲FC近似下mu_e是个常量，是指关于R（核坐标）。mu_e实际上依赖于电子坐标。请问是这样吗@Sobereva 老师？这里对τ_e进行积分，这个τ_e是指电子坐标吗？那自旋部分的τ_s积分又是指什么变量？
A:
是电子坐标
另一个就是自旋坐标，就两个值，分别对应自旋函数为alpha和beta
----------------------------------------------------
2015.10.31 14:52:21
Q:
好的，谢谢老师，对于自旋坐标我还要再理解一下~
之前还提了一个问题，就是<psi'|mu_e|psi>这个积分的数值能如何计算？能用multiwfn导出吗？
不怕说错，请 @Sobereva 老师看看我想的是不是正确？
1、mu_e其实是电子跃迁偶极矩算符，这么<psi'|mu_e|psi>得到的也就是一组轨道间跃迁的偶极矩了，所有的跃迁偶极矩 Σ 之后，就是总的跃迁偶极矩，也就是gaussian输出的那个值（x,y,z三个分量）了。
2、此外，如果先将所有的psi和psi'分别组合成基态轨道的总波函数PSI和激发轨道的总波函数PSI'，再求<PSI'|mu_e|PSI>，得到的也就是总的跃迁偶极矩即gaussian输出的结果了？（不知道mu_e需不需要组合？）
谢谢老师指点~~
A:
就是两个电子态的跃迁偶极矩，高斯做激发态计算直接就输出了
@莫哈莫哈 高斯激发态计算输出各种轨道跃迁。你对每对儿轨道计算跃迁偶极矩，乘上组态系数，对所有的加和，就是两个电子态的跃迁偶极矩。
如果PSI、PSI'是两个电子态，应当为<PSI'|Σmu_e_i|PSI>，其中mu_e_i是i电子的偶极矩算符
----------------------------------------------------
2015.10.31 15:28:23
Q:
sob老师，，我都不好意思了，计算出错了，
[Image]
输入文件是%chk=SII1-I.chk
%mem=200MW
%nprocshared=2
#p opt freq m062x/6-31+g(d,p) guess(fragment=2) pop=nbo
Title Card Required
0 1 0 3 0 -3
C(Fragment=1) -1.21236500 -0.23437800 -0.03046900
H(Fragment=1) -2.12499000 0.32626400 -0.12438200
H(Fragment=1) -1.25843300 -1.28427600 0.19193000
C(Fragment=2) 0.09993400 0.48290900 0.03018400
H(Fragment=2) 0.18980500 1.21912100 -0.76695600
O(Fragment=2) 1.20733800 -0.34845600 -0.03305000
H(Fragment=2) 0.20949900 1.03535600 0.96551200
我么？去掉opt?意思是不用优化结构了么
是说去掉opt后计算，如果计算没出错，就在这基础上在用opt?
先这样#p m062x/6-31+g(d,p) guess(fragment=2)
sob老师，您的意思是先去掉Opt得到一个合适的初猜，合适的波函数，因为它对分子片段进行分割的缘故吗？
然后#p opt freq m062x/6-31+g(d,p)
A:
去掉opt
优化。你先弄到收敛的波函数，然后优化的时候读取
开发者考虑不了那么周全，这明显就是片段初猜和opt一起用的时候出bug
y
然后#p opt guess=read m062x/6-31+g(d,p)
这样更稳妥
#p opt guess=read um062x/6-31+g(d,p) nosymm
PS:此体系用弥散没意义
----------------------------------------------------
2015.10.31 15:38:29
Q:
老师，前一段有位做有机的朋友说起离子对的问题，说有机上很常见，那计算时该怎么考虑？离子对一边边带负电，一边边带正电，应该向R2D2的一样去分片吗？
@Sobereva 那怎么考虑有机上所说的离子对？
是啊，所以搞得自己有些郁闷，的确有机体系在搞有机的朋友看来是带正电的，而算出来是负的，更有霸气的朋友会说我不管，你必须给我算出个正的来，结果不知道该怎么弄！
请问老师们，对于和有机人理解有些出入的结果，有没有推荐的好文章或帖子？求推荐，尤其是电荷相关的
我只试过NPA和MA
其他有哪些更好的？
A:
分那个没用，收敛完了都一样
算出来电荷怎么分布的就是怎么分布的
有机上很多概念纯属臆测
比如有机上说某原子带+1电荷，实际上算出来往往发现根本不是那么回事，也就比如0.1电荷甚至负电荷都有。
看群共享里原子电荷计算方法的对比
也考察不同的理论方法，包括M06-2X、MP2
低HF成份泛函的离域性误差某些情况造成电荷分布不很合理
----------------------------------------------------
2015.10.31 18:45:14
Q:
[图片]老师，原先箭头指的地方色彩渲染不出来，我按照RDG-maxrho设置为100的方法调整，结果变这样了@Sobereva
那0.10怎么样？
A:
你把RDG_maxrho设小点，比默认的大一点就足够显示出那部分了。你设成100的话，成键区域都出来了
或者你在后处理界面设定sign(lambda2)rho在某个范围以外时RDG为100，这样既能让箭头部分显示出来，化学键部分又不显示出来
看散点图判断是最好的
手册4.100.1节靠后部分有图例
比如从散点图上看那个spike差不多在0.08位置，那就比如把-9999到-0.05部分RDG值设100，0.1到9999部分RDG也设100。这样就只有sign(lambda2)rho在-0.05~0.08区域是保留的
----------------------------------------------------
2015.10.31 20:40:59
Q:
@Sobereva 老师，请问二群的群号是多少啊
论坛其实很好的，sob回复挺快的。还有很多好心人
qq群消息一刷就没了，而且后加入群的也不知道
A:
二群时机未到
副群是闲聊用的
二群!=副群
二群是必然的，但尽量让讨论质量基本与本群保持一致
N个月后才有呢
PS：群的人数是有限的，但是计算化学公社的注册人数是无限的
是的，所以我也更愿意在那里回复时花更多精力
----------------------------------------------------
2015.10.31 20:42:30
Q:
各位大侠：本人做一个激发态频率计算，发现运行了两星期了，还没有收敛。大家给看下命令有啥问题没？ #p freq m052x/6-311++g(3d,2p) td=(singlets,root=1) integral=ultrafine
A:
@HERO 不是没收敛，而是在做有限差分计算，慢得很。基组完全不需要这么大，太浪费时间。
----------------------------------------------------
2015.10.31 20:55:18
Q:
不设管理员也行，招募志愿者对聊天记录进行整理，在论坛发布
A:
有人愿意的话可以私下联系
----------------------------------------------------
2015.11.01 06:22:12
Q:
老师，我又来问问题了，下午呢按照先#p m062x/6-31+g(d,p) guess(fragment=2)，再#p opt guess=read um062x/6-31+g(d,p) nosymm计算了一个基态的物质，第一步打开log文件显示[图片]
继续按第二步算也没有出错，但是一算过渡态第二步用#p opt=(calcfc,ts,noeigen) um062x/6-31+g(d,p) nosymm guess=read，优化出来的结构就不是想要的了，这个是什么原因呢？
高斯不是默认计算对称性的么
sob老师建议我用nosymm的
多试试，我的分子小，计算很快的
然而去掉了nosymm并没有什么用@一只耳！ /大哭
broken symmetry 说的是对称破缺么
unrestricted是加在泛函前面的，非限制性，如uhf,对称破缺的话，我理解是破坏分子的对称性，以得到想要的构型，理解的很浅，还不造对不对
哦哦，原来是这样，谢谢虎王拨乱反正
A:
@R2D2 这根对称性没关系。
简单来说，如够用U结合guess=mix就能得到对称破缺波函数，就别用片段初猜，因为那样更省事
----------------------------------------------------
2015.11.01 06:24:32
Q:
为庆祝《物理化学学报》出版发行30周年，应主编邀请，系统总结我本人和合作者近年工作的综述文章 “密度泛函活性理论中的信息论方法” 于昨天上线。文尾特别致谢了该群的群主！近年来与群主的合作卓有成效。欢迎有兴趣的群友免费下载批评斧正。该文链接是：http://www.whxb.pku.edu.cn/CN/abstract/abstract29268.shtml
A:
这篇综述非常有益/强
----------------------------------------------------
2015.11.01 07:07:27
Q:
B3LYP(VWN5)和B3LYP(VWN3)的计算结果会有大的差别吗？
A:
有一点差别，但不很大
此文专门有对比数据：Chemical Physics Letters 268 (1997) 345-351
----------------------------------------------------
2015.11.01 08:27:12
Q:
lapack里面dsyev最后三个参数怎么用？
A:
Multiwfn里把dsyev封装成了简单形式可以参考，调用起来方便多了
[图片]
----------------------------------------------------
2015.11.01 08:31:22
Q:
In particular, he is grateful to Dr. LU Tian of Beijing
Kein Research Center for Natural Sciences for coding all the formulas in the text in the MultiWFN package,MULTIwfn可以做这方面的计算吗？
所谓的“密度泛函活性理论中的信息论方法，就是用于量化研究预测分子的化学反应活性？
老师，是这个意思吗
基于量化计算获取波函数时有啥普遍的原则吗，比如量化计算的方法水平等该如何确定？
A:
当然
一系列文章的计算都是基于Multiwfn做的
基于量化计算得到的波函数，用Multiwfn计算出Hirshfeld电荷以及各种局部形式的信息熵及相关的量，预测反应位点、反应活性
M06-2X、B3LYP等常用的就行，基组用6-311G**就够了
对于计算级别不是很敏感
----------------------------------------------------
2015.11.01 09:21:14
Q:
谢谢sob老师的解答，6点多就开始工作了，，老师真是辛苦。我现在有这样一个疑问，m062x的hf成分有54%，b3lyp的hf成分是20%，我想问的是这个hf成分的多少是怎么来影响我的计算的，比如我算的双自由基，您建议我选hf成分高的。如果是正常的体系中hf成分的多少又有多大影响呢？
恩恩，为什么会有这些趋势呢？我是想能更深入点，（我有可能看不懂/大哭，还是想知道，）
A:
[图片]
要想了解深入，就得自己看DFT专著
不是群里几句话能讲明白的
----------------------------------------------------
2015.11.01 10:00:37
Q:
群里的老师们，麻烦帮下忙，昨天出现这种情况，有个老师说要装python，我装好了。还是这个样子，怎么办呢
还是这样呢
[Image]
好的，我试试，谢谢老师
这样看吗？这是没有那个字符吧？
A:
如果装好，直接输入python就应当能启动之
[图片]
which python看看到底装哪去了
如果comblike.gala是个脚本，你用cat查看之，看看里面的/usr/bin/python后面有没有^M字符，有的话用dos2unix转换一下comblike.gala
先dos2unix转换一下再说
----------------------------------------------------
2015.11.01 11:10:35
Q:
老师可以推荐一下比较好的DFT专著嘛
A:
最入门的A Chemist's Guide to Density Functional Theory 2ed
----------------------------------------------------
2015.11.01 11:24:44
Q:
各位老师我有个问题，想问个问题，虽然这时候问好像，，，不是时候，我想问的是我计算过渡态时，1.用低精度的优化出构型，在高精度计算却得不到想要的构型了，怎么办？2.如果低精度的优化不出来是不是就肯定优化不出来？
我用的是hf/6-31+G**和 b3lyp/6-31+G**，m062x/6-31+G**，后边两个是不是您说的中精度啊，我理解者CCSD啥的才是高精度，我刚才说的是个相对，，@SET
A:
@R2D2 低精度优化得出来，高精度优化不出来是正常情况
高精度死活无法重现就说明不存在那个过渡态
以m062x/6-311G**为准
HF无意义
对于有机分子，B3LYP不会比M062X更好
弥散通常对于优化过渡态没意义，远不如用3-zeta
----------------------------------------------------
2015.11.01 11:30:26
Q:
@Sobereva 老师，我的老板让我去重复一个稿子里的结果，但是我实在没办法优化出那个结构，我想按照那个构型进行单点计算，这个可以吗
@Sobereva 好的，谢谢老师
A:
可以要啊
很多文章都把计算中得到的结构放到补充材料里，grimme也总是这样做，很鼓励这样
----------------------------------------------------
2015.11.01 11:33:59
Q:
恩恩，我理解@SET 老师的意思了。@Sobereva 老师的意思是一般都不加弥散函数是吧，我在小木虫上读了您关于这个的帖子，里面说的计算弱相互作用、极化率、反应势垒、旋光性、高阶激发态、研究孤对电子等用加弥散函数，我想问一下，关于反应势垒，我目前就是只用freq得到的值算算热力学的焓值差，吉布斯自由能差，这是不是您说的反应势垒？
A:
我特别讨厌看到小木虫这三个字，而且还是在那上看我的文章...
谈谈弥散函数和“月份”基组
http://sobereva.com/119
优化过渡态、做振动分析不需要弥散函数。高精度算势垒的时候，单点计算加上弥散会改进结果，是值得的。
----------------------------------------------------
2015.11.01 11:36:18
Q:
M062X可不可以算第三周期大体系结构
A:
@galois 第三周期都是主族的，没问题
----------------------------------------------------
2015.11.01 11:50:51
Q:
各位好，请教个问题，你们用哪个软件计算物质的pKa啊，ACD i-lab登录不上去了啊
A:
http://www.vcclab.org/lab/alogps
http://www.chemaxon.com/marvin/sketch/index.jsp
http://apps.ideaconsult.net:8080/ToxPredict
----------------------------------------------------
2015.11.01 15:54:18
Q:
请问，我用A triple-ζ Slater-type basis set, extended with a single polarization function, was used to describe all atoms. 描述def2tzvp基组对不对？
哦，谢谢老师
A:
不对
def2TZVP是基于高斯函数的基组
极化函数描述也错误
def2TZVP的极化函数不少
----------------------------------------------------
2015.11.01 21:02:38
Q:
multiwfn计算Hirshfeld surfaces时候，如何得到百分数吗？@Sobereva
[图片]
类似这样的C-H相互作用所占的百分含量
A:
什么百分数？
没有现成的这个功能
记错了，有这功能，看手册4.12.6节，有例子
----------------------------------------------------
2015.11.02 06:50:27
Q:
关于高斯怎么计算分析自旋密度分布的问题，懂得老师克服给说一下
A:
@十三 用Multiwfn计算，主功能3、4、5可以方便地分别绘制自旋密度曲线图、平面图和等值面图，群共享Multiwfn与波函数分析简介(2015-Jun-11)里有图例，使用方法手册第4章有一堆例子。也可以用Multiwfn计算自旋布居，手册4.A.2节有说明
----------------------------------------------------
2015.11.02 06:50:27
Q:
请问有办法计算某个化合物的不同共振式结构的稳定性吗？多谢
就是通过计算比较不同共振式结构的稳定性
A:
@华仔 可以做NRT分析，看不同共振式的比例来判断
----------------------------------------------------
2015.11.02 08:30:54
Q:
请问，对于稀土体系，自旋多重度为2，优化结果有两套轨道能量 Alpha和Beta，我应该看哪套呢？还是两套都要描述？
A:
都考虑
----------------------------------------------------
2015.11.02 17:18:58
Q:
各位老师，我感觉我从最基础的地方就没搞清，就是化学动力学和热力学的本质含义，有没有相关论文和文献推荐我研读一下，我想从源头上开始纠正自己
那能否推荐一下书籍呢
土豪...
A:
不看书学好基础，越看文献越糊涂
买本Levine的quantum chemistry 6ed，市面就有
你买影印版啊
如果没有影印版、没有电子版，卖1w我也买，知识无价
不是土豪，是价值观的问题
影印版才100出头
正版的是115元，链接那个是蒙人的
淘宝一搜就有啊，有什么难的。
100出头的都没问题
[图片]
[图片]
第二张图不是，那个是习题解
要看目录看群共享里第5版，虽然第6版章节发生了变动
----------------------------------------------------
2015.11.02 17:18:58
Q:
各位老师，我感觉我从最基础的地方就没搞清，就是化学动力学和热力学的本质含义，有没有相关论文和文献推荐我研读一下，我想从源头上开始纠正自己
A:
不看书学好基础，越看文献越糊涂
买本Levine的quantum chemistry 6ed，市面就有
影印版才100出头
正版的是115元，链接那个是蒙人的
淘宝一搜就有啊，有什么难的。
100出头的都没问题
[Image]
----------------------------------------------------
2015.11.02 17:32:22
Q:
IMOMM模型得到的能量=QM+MM吗？
[图片]
一篇文献里的化合物能量，直接是两层的能量加和
那么对于ONIOM方法里的呢？
A:
你可以当成等于
IMOMM是ONIOM的特例
----------------------------------------------------
2015.11.02 18:08:55
Q:
各位老师，在Gamess中 6-311g **基组怎么输啊，我是刚学做Gamess
A:
STO-3G：GBASIS=STO NGAUSS=3
3-21G：GBASIS=N21 NGAUSS=3
6-31G：GBASIS=N31 NGAUSS=6
6-31G*：GBASIS=N31 NGAUSS=6 NDFUNC=1
6-31G**：GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1
6-31+G**：GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.T.
6-31++G**：GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.T. DIFFS=.T.
6-311系列改为GBASIS=N311，如
6-311++G(3df,3p)：GBASIS=N311 NGAUSS=6 NDFUNC=3 NFFUNC=1 NPFUNC=3 DIFFSP=.T. DIFFS=.T.
----------------------------------------------------
2015.11.02 18:26:57
Q:
请教group-additivity 如何用中文讲？
A:
@心凡 基团加和性
----------------------------------------------------
|
评分 Rate
-
查看全部评分 View all ratings
|
发表于 Post on 2015-11-4 02:08:20
收藏 Add to favorites
然后 看到我问的问题了