计算化学公社

标题: NBO分析中donor-acceptor 的转移电子数Delta q的公式 [打印本页]

作者
Author: diaolanxinyu 时间: 2014-11-30 10:51
标题: NBO分析中donor-acceptor 的转移电子数Delta q的公式
麻烦大家看一下在NBO分析中donor-acceptor 的转移电子数的几个问题：
1. 焢萀汶式对吗? 见附件delta_q.png
2. 怎么得到的呢？我知道627.51是为了统一单位kcal/mol
3.如果公式对的话，焢萀绝对值越大，转移电子数越多的吧？

谢谢大家！

作者
Author: diaolanxinyu 时间: 2014-11-30 10:53
问题1出现了错误，不能显示特殊字符，应该是Delta q的公式对吗？见附件delta_q.png

作者
Author: sobereva 时间: 2014-11-30 11:43
不对。见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=102里面Natural Bond Orbital Methods.pdf中的公式18，那个才是对的。
不知道你的公式怎么来的，应当以那个帖子里的文档里的公式为准。

原则上数值越大转移量越大，但应当注意正确理解含义，NBO有很大局限性，而且要注意这种方式衡量电子转移数只是定性的。直接计算原子电荷求差值才是最好的作法。

作者
Author: diaolanxinyu 时间: 2014-11-30 12:57
谢谢sobereva老师！

您说直接计算原子电荷求差值才是最好的作法，我想问一下，需要计算谁和谁的电荷然后求差值呢？怎样计算呢？

我想看一下分子内部电荷转移的情况，您说的这个原子电荷求差值可以实现吗？

作者
Author: sobereva 时间: 2014-11-30 13:27
怎么求差值，这都要看实际体系，随机应变，没有实际体系没法说

作者
Author: diaolanxinyu 时间: 2014-11-30 16:01
附件中是用Gaussian算的两个不同时刻结构的.fchk文件，方法：#p b3lyp/6-31g(d) pop=nbo 电荷和自旋多重度是-1，2（贴敷一个电子）
一个是0fs，一个是60fs，我想看一下这个贴敷的电子是如何促使C-C键断裂的，或者这个贴敷电子在轨道之间是如何传递的？
我一开始的想法是分析NBO信息中的“Perturbation Theory Energy Analysis”部分，分析不同时刻结构中E(2)较强的一些轨道，所以想算出E(2)较强的一些轨道中转移电子数，来给出定性的分析。

当您提到“直接计算原子电荷求差值才是最好的作法”，我想用一个具体的例子来请教一下sobereva老师您说的这个方法该怎么做呢？

谢谢！

作者
Author: sobereva 时间: 2014-11-30 23:08
你可以考察-1和中性状态下原子电荷的差值，或者求密度差来绘图，这样多出来的一个电子怎么分布的就清楚了。你也可以计算这两个状态下的键级，以考察多一个电子是如何影响键的强弱的。
NBO的E2分析对这个问题没直接用处，不管你加不加这个电子，照样会给出CT值，何况这还是近似方式计算的，

作者
Author: diaolanxinyu 时间: 2014-12-1 10:31
非常感谢sobereva老师！谢谢！:)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)