计算化学公社
标题: 轨道成份分析咨询 [打印本页]
作者Author: hzfish 时间: 2017-5-23 18:55
标题: 轨道成份分析咨询
对A(B)n（n=6）分子进行轨道成份分析，假设A为一个金属阳离子，B为水分子。
问题：
1. 进行片断分析，A和每个B为一个片断，要得到第i轨道每个片断的贡献，得分析1+n？
2. 假设第i轨道中，第2个B片断贡献为a%，而独立B轨道为1a1、2a1、1b2、3a1、1b1，那如何判断第2个B片断中哪个轨道对第i轨道产生贡献a%？即产生贡献是1a1、2a1、1b2、3a1、1b1中的哪个？
3. 假设第i轨道中，第2个B片断贡献为a%，那么第i轨道中第2个B片断贡献部分能不能单独图形化？
谢谢！
作者Author: sobereva 时间: 2017-5-24 02:41
1 是的
2 貌似你要考察片段轨道对复合物轨道的贡献，这得用CDA模块才行，可以输出复合物轨道当中各个片段轨道的成份。见
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166
3 最简单的做法，在主界面先选-3 Obtain a fragment contribution to property，然后输入要考察的那个片段的原子序号，之后你在功能0里看到的轨道图形就只保留纯粹属于那个片段的部分了。此文也有提到
使用Multiwfn观看分子轨道
http://sobereva.com/269
作者Author: hzfish 时间: 2017-5-24 14:33
谢谢社长解答。
还一点疑问，片断在第i条轨道的成份和片断轨道对复合第i轨道贡献不是一个概念？两者原理上有什么区别？
社长能不能简单用公式解释一下？
谢谢！
作者Author: sobereva 时间: 2017-5-25 02:29
不是一个概念
前者是把复合物轨道向AO基函数展开，通过加和来求指定片段的贡献
后者是把复合物轨道向片段轨道展开，直接求某个片段轨道的贡献
片段轨道是你计算单体时候得到的，其形状和复合物轨道当中直接抠出来某个片段是不同的。而且基于片段轨道讨论，才可以说片段轨道的轨道能量、不可约表示
作者Author: hzfish 时间: 2017-5-25 10:06
谢谢社长耐心回复。
以下是我们理解，社长帮着看下是否正确。
假设所有轨道均已正交归一化
复合物轨道如下：
Ψ_1=C_11*φ_1 + C_12*φ_2 + C_13*φ_3 + C_14*φ_4 （1）
Ψ_2=C_21*φ_1 + C_22*φ_2 + C_23*φ_3 + C_24*φ_4 （2）
Ψ_3=C_31*φ_1 + C_32*φ_2 + C_33*φ_3 + C_34*φ_4 （3）
Ψ_4=C_41*φ_1 + C_42*φ_2 + C_43*φ_3 + C_44*φ_4 （4）
φ_1在Ψ_1的贡献为C_11^2
分成两个片断，片断1由1、2原子组成，片断2由3、4原子组成，片断轨道如下
断片1在Ψ_1的贡献为C_11^2 + C_12^2
片断1轨道
ψ_11=C_111*φ_1 + C_112*φ_2 （5）
ψ_12=C_121*φ_1 + C_122*φ_2 （6）
片断2轨道
ψ_21=C_211*φ_3 + C_212*φ_4 （7）
ψ_22=C_221*φ_3 + C_222*φ_4 （8）
由片断轨道组合得到的复合物轨道
Ψ_1=c_11*ψ_11 + c_12*ψ_12 + c_21*ψ_21 + c_22*ψ_22 （9）
Ψ_2=c_21*ψ_11 + c_22*ψ_12 + c_23*ψ_21 + c_24*ψ_22 （10）
Ψ_3=c_31*ψ_11 + c_32*ψ_12 + c_33*ψ_21 + c_34*ψ_22 （11）
Ψ_4=c_41*ψ_11 + c_42*ψ_12 + c_43*ψ_21 + c_44*ψ_22 （12）
片断1轨道ψ_11在Ψ_1贡献为c_11^2
把（5）、（6）、（7）（8）代入（9）得，
Ψ_1= (c_11*C_111 + c_11* C_121)*φ_1 + (c_12*C_112 + c_12* C_122) *φ_2
+ (c_21* C_211 + c_21* C_221) *φ_3 + (c_22* C_212 + c_22* C_222) *φ_4 （13）
（13）式和（1）比较，
C_11= c_11*C_111 + c_11* C_121
C_12= c_12*C_112 + c_12* C_122
C_13= c_21* C_211 + c_21* C_221
C_14= c_22* C_212 + c_22* C_222
作者Author: sobereva 时间: 2017-5-25 11:08
你看一下CDA原文吧，这看着太费劲。CDA计算时候都是用矩阵形式来搞的，就是把复合物轨道从AO基变换成片段轨道基，通过片段轨道的系数矩阵很容易实现
作者Author: hzfish 时间: 2017-5-25 12:03
谢谢！
我看下CDA原文。
有问题再请教。
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |