计算有机分子的LUMO HOMO的疑问

exity · 发表于 Post on 2017-4-17 10:23:13

有请大神解惑，我试着尽量把问题描述的清楚一点。我想对一个有机分子进行计算后得到各个原子对总的LUMO HOMO的贡献，然后通过比较这些值来判断氧化剂进攻改分子时可能发生断键的位点，主要参考"Rapid Removal of Tetrabromobisphenol A by Ozonation in Water: Oxidation Products, Reaction Pathways and Toxicity Assessment"这篇文献。
文献里说可以通过计算LUMO HOMO的平方和之类的值判断断键位点并给出了下图：

于是我按照文献里方法用如下的关键词对一个三嗪类化合物进行了计算

# opt freq b3lyp/6-311g(d,p) scf=qc pop=all scrf=(solvent=water)

计算过程很顺利，也没有出错。

结果算完了以后，用“Molecular Orbital Coefficients"为关键词在输出文件里搜索，什么也找不到，我自己也不知道应该怎么去找LUMO HOMO的值，烦请大神解惑。是我的关键词不对，还是关键词参数不对？

此外，我有一系列的化合物，如下图

我在用同一种氧化剂对这些物质进行降解研究，得到的结论是都在氨基的右侧发生断键后生成苯胺和对应的烷烃，但是母体的降解速率各不相同。

我分别对这些分子进行计算，得到每个分子中的氮原子和相邻的碳链上的碳原子对各个分子的LUMO HOMO的贡献，能否通过比较这些值来对降解速率进行说明？

我是计算新手，上面肯定有一些蠢问题，还请大神不要笑我！谢谢！

PS 输入输出文件见附件

sobereva · 发表于 Post on 2017-4-17 18:29:36

Multiwfn是计算轨道成份最方便最强大的工具，这部分功能的操作和介绍看：
谈谈轨道成份的计算方法
http://sobereva.com/131
Multiwfn入门贴：
Multiwfn入门tips
http://sobereva.com/167

另外关于反应位点预测的文章看http://bbs.keinsci.com/forum.php?mod=viewthread&tid=384，特别是《亲电取代反应中活性位点预测方法的比较》（http://www.whxb.pku.edu.cn/CN/abstract/abstract28694.shtml）

exity · 发表于 Post on 2017-4-17 19:41:37

sobereva 发表于 2017-4-17 18:29
Multiwfn是计算轨道成份最方便最强大的工具，这部分功能的操作和介绍看：
谈谈轨道成份的计算方法
http:/ ...

谢谢Sob,我也正在研究你编写的这个软件，命令行的输入方式很燃，像黑客帝国。不过，关于我这个问题，您能不能就事论事的点评一下呢？我真的很希望能搞清楚这个事情。

冰释之川 · 发表于 Post on 2017-4-17 20:26:38

exity 发表于 2017-4-17 19:41
谢谢Sob,我也正在研究你编写的这个软件，命令行的输入方式很燃，像黑客帝国。不过，关于我这个问题，您能 ...

Sob老师回答得很清楚了：
（1）算各个原子对总的LUMO HOMO的贡献参考：
谈谈轨道成份的计算方法
http://sobereva.com/131

（2）不会用Multiwfn的参见：
Multiwfn入门tips
http://sobereva.com/167

（3）研究反应位点还有其他更好的方法，你可以参考使用其它方法：
反应位点预测的文章看http://bbs.keinsci.com/forum.php?mod=viewthread&tid=384，特别是《亲电取代反应中活性位点预测方法的比较》（http://www.whxb.pku.edu.cn/CN/abstract/abstract28694.shtml）

exity · 发表于 Post on 2017-4-17 20:34:21

冰释之川发表于 2017-4-17 20:26
Sob老师回答得很清楚了：
（1）算各个原子对总的LUMO HOMO的贡献参考：
谈谈轨道成份的计算方法

好吧，也就是说我用opt freq b3lyp/6-311g(d,p) scf=qc pop=all scrf=(solvent=water)这样的关键词算出来的fch文件也能通过那个软件算出各自的LUMO HOMO恭贡献？

sobereva · 发表于 Post on 2017-4-17 23:02:10

exity 发表于 2017-4-17 20:34
好吧，也就是说我用opt freq b3lyp/6-311g(d,p) scf=qc pop=all scrf=(solvent=water)这样的关键词算出来 ...

是

根本没有pop=all的关键词，pop=full写了也毫无意义。没事不要用scf=qc。

另见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=334

exity · 发表于 Post on 2017-4-17 23:30:27

sobereva 发表于 2017-4-17 23:02
是

根本没有pop=all的关键词，pop=full写了也毫无意义。没事不要用scf=qc。

谢谢sob！以后不会用好吧做句首发语词了，一定注意。

exity · 发表于 Post on 2017-4-18 12:17:58

sobereva 发表于 2017-4-17 23:02
是

根本没有pop=all的关键词，pop=full写了也毫无意义。没事不要用scf=qc。

SOB老师，再次打扰您一下！按照您给的指示，我进行了计算。得到我的有机分子FCHK文件后用Multiwfn打开了，用Hirshfeld方法进行轨道分析，那么我又不清楚那个是LUMO 那个是 HOMO 轨道，我就用GVIEW打开FCHK文件，在右键 edit 里 MOs中，找到这样的图

是不是通过这里就知道第38是LUMO 而 37是HOMO，然后回到Multiwfn中，选择对应的原子编号，然后轨道输入 37,38 输出的就是对应的LUMO HOMO贡献？
比如，我想知道氮原子(8号)对LUMO HOMO的贡献，我就输入，8 enter, 37,38 enter 得到这样的输出

Input atom index, e.g. 4
8
Input orbital range e.g.  20,25
37,38
Orb# Type    Ene(a.u.)    Occ Composition Population
37 Alpha&Beta    -0.2625 2.000    16.089%    0.321773
38 Alpha&Beta    -0.0329 0.000    3.206%    0.000000
Population of this atom in these orbitals: 0.321773

是不是Ene就是我要的能量，而Composition就是对轨道贡献的百分比？

同时Ene对应的就是这个文献插图的中的最上面结构中的蓝色和黑色小数字，不知道我这样理解对不对？

麻烦SOB老师点拨一二！

sobereva · 发表于 Post on 2017-4-18 20:21:42

exity 发表于 2017-4-18 12:17
SOB老师，再次打扰您一下！按照您给的指示，我进行了计算。得到我的有机分子FCHK文件后用Multiwfn打开了 ...

轨道成份算得对

ene是轨道能量，跟你图上黑字和蓝字完全没关系。具体什么含义在图注中应当会写

哪个是LUMO哪个是HOMO进入主功能0的时候直接在屏幕上提示了，根本用不着gview。看
使用Multiwfn观看分子轨道
http://sobereva.com/269

exity · 发表于 Post on 2017-4-19 09:35:14

sobereva 发表于 2017-4-18 20:21
轨道成份算得对

ene是轨道能量，跟你图上黑字和蓝字完全没关系。具体什么含义在图注中应当会写

SOB老师，谢谢您，我这里有一个建议，不知道对不对！
对于我这样的超级没入门者来说，有些基本概念还是欠缺不少，但是迫于进度压力，必须进行一些量化计算。
我发现有一个地方可能对初学者（超级菜鸟）有误导：例如，你的“http://sobereva.com/269”例子中，对LUMO. HOMO轨道的判断如下

图形窗口上方中，点击Orbital info.--show all，在文本窗口中就会一目了然地输出所有轨道的信息。例如：
Orb: 19 Ene(a.u./eV): -0.340271    -9.2592 Occ:  2.000000 Type: A+B
Orb: 20 Ene(a.u./eV): -0.246291    -6.7019 Occ:  2.000000 Type: A+B
Orb: 21 Ene(a.u./eV): -0.246291    -6.7019 Occ:  2.000000 Type: A+B
Orb: 22 Ene(a.u./eV):    0.003620    0.0985 Occ:  0.000000 Type: A+B
Orb: 23 Ene(a.u./eV):    0.003620    0.0985 Occ:  0.000000 Type: A+B
显示了a.u.和eV为单位的所有轨道能量、占据数和自旋类型。此例A+B就是Alpha+Beta，即闭壳层轨道。从这里也可以立刻看出21和22号分别是HOMO和LUMO，20是HOMO-1，23是LUMO+1。
刚好21,22的能量从负值过渡到正值，我开始以为就是通过这样的正负变化来判断LUMO HOMO的（大误！），刚好在我自己的例子中

Orb: 55 Ene(a.u./eV): -0.280158    -7.6235 Occ:  2.000000 Type: A+B
Orb: 56 Ene(a.u./eV): -0.251841    -6.8529 Occ:  2.000000 Type: A+B
Orb: 57 Ene(a.u./eV): -0.250066    -6.8046 Occ:  2.000000 Type: A+B
Orb: 58 Ene(a.u./eV): -0.031105    -0.8464 Occ:  0.000000 Type: A+B
Orb: 59 Ene(a.u./eV): -0.010169    -0.2767 Occ:  0.000000 Type: A+B
Orb: 60 Ene(a.u./eV):    0.016059    0.4370 Occ:  0.000000 Type: A+B
Orb: 61 Ene(a.u./eV):    0.043075    1.1721 Occ:  0.000000 Type: A+B
Orb: 62 Ene(a.u./eV):    0.049899    1.3578 Occ:  0.000000 Type: A+B

正确的LUMO HOMO应该是58,57，可是我最开始的时候以为是59,60（通过正负变化判断的，其实这是错误的判断方法），和GVIEW里的结果是不一样的。然后突然我似乎顿悟了一样，是不是通过2.00000-0.00000的变化来判断LUMO HOMO，这里的2和0就是指占据轨道的电子数，2就是有电子占据的所以57是HOMO, 58就是LUMO。这样就比较没有误会了，SOB老师可以添加一句话，解释一下2,0的意义，就会对我这样的超级菜鸟帮助特别大，不知道我说的对不对，请批评。

另外，现在基本上我已经知道了怎么看各个原子对LUMO HOMO的贡献了，那么如果我的分子结构比较大，例如有60多个原子，我想知道其中的40个原子对LUMO HOMO的贡献，岂不是要输入40*2次（原子序数和轨道数），然后输出的结果也要分别的黏贴复制到一起才能再继续下面的统计分析的工作，有没有简便一点的方法，烦请不吝赐教！

sobereva · 发表于 Post on 2017-4-19 10:00:24

exity 发表于 2017-4-19 09:35
SOB老师，谢谢您，我这里有一个建议，不知道对不对！
对于我这样的超级没入门者来说，有些基本概念还是 ...

没有任何误导，懂得HOMO、LUMO含义的人都不会犯这种错误

HOMO、LUMO序号是多少，屏幕上都直接输出了，都根本不用看占据数，没法更贴心了

Hirshfeld轨道成份分析界面里很清楚有一个定义片段的功能，定义完片段，看轨道成份的时候一下子就输出片段的贡献了

exity · 发表于 Post on 2017-4-19 11:17:15

sobereva 发表于 2017-4-19 10:00
没有任何误导，懂得HOMO、LUMO含义的人都不会犯这种错误

HOMO、LUMO序号是多少，屏幕上都直接输出了， ...

明白了，谢谢您！你真是一本真人版的量化百科全书！

exity · 发表于 Post on 2017-4-19 11:21:47

sobereva 发表于 2017-4-19 10:00
没有任何误导，懂得HOMO、LUMO含义的人都不会犯这种错误

HOMO、LUMO序号是多少，屏幕上都直接输出了， ...

Orbital 57 is HOMO, energy: -0.250066 a.u. -6.804646 eV
Orbital 58 is LUMO, energy: -0.031105 a.u. -0.846420 eV
LUMO/HOMO gap: 0.218961 a.u. 5.958227 eV 574.881600 kJ/mol

原来在这个地方，非常感激，真的是不能再贴心了。\笑中带泪

然后我用VMD输出LUMO HOMO图，发现可以实现原子编号（也是参考了你的教程），不过不能加原子名称，比如我要把C, N, S这样的原子符号标上去，请问应该怎么实现？

exity · 发表于 Post on 2017-4-19 13:52:57

sobereva 发表于 2017-4-19 10:00
没有任何误导，懂得HOMO、LUMO含义的人都不会犯这种错误

HOMO、LUMO序号是多少，屏幕上都直接输出了， ...

SOB老师，我还有一个问题想请教。

应该怎么计算有机分子中各个原子在LUMO HOMO下的原子轨道？

比如下图这样

这个计算结果就把各个分子轨道下原子轨道的贡献给出了，我想要得到这样的细节数据应该怎么输入关键词。

我现在用的关键词如下，得到OUT文件里没有上图这样的结果。

%nprocshared=8
%chk=DIA.chk
# opt freq b3lyp/6-311g(d,p) scrf=(solvent=water) scf=qc

lao7 · 发表于 Post on 2017-4-19 14:43:10

exity 发表于 2017-4-17 20:34
好吧，也就是说我用opt freq b3lyp/6-311g(d,p) scf=qc pop=all scrf=(solvent=water)这样的关键词算出来 ...

咱们一样弱，跟不上拍，呵呵！

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

[新手求助] 计算有机分子的LUMO HOMO的疑问

评分 Rate

浏览过的版块