轨道成份分析咨询

hzfish

对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片断贡献部分能不能单独图形化？

谢谢！

sobereva

1 是的
2 貌似你要考察片段轨道对复合物轨道的贡献，这得用CDA模块才行，可以输出复合物轨道当中各个片段轨道的成份。见
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166
3 最简单的做法，在主界面先选-3 Obtain a fragment contribution to property，然后输入要考察的那个片段的原子序号，之后你在功能0里看到的轨道图形就只保留纯粹属于那个片段的部分了。此文也有提到
使用Multiwfn观看分子轨道
http://sobereva.com/269

hzfish

sobereva 发表于 2017-5-24 02:41
1 是的
2 貌似你要考察片段轨道对复合物轨道的贡献，这得用CDA模块才行，可以输出复合物轨道当中各个片段 ...

谢谢社长解答。

还一点疑问，片断在第i条轨道的成份和片断轨道对复合第i轨道贡献不是一个概念？两者原理上有什么区别？
社长能不能简单用公式解释一下？
谢谢！

sobereva

hzfish 发表于 2017-5-24 14:33
谢谢社长解答。

还一点疑问，片断在第i条轨道的成份和片断轨道对复合第i轨道贡献不是一个概念？两者原 ...

不是一个概念
前者是把复合物轨道向AO基函数展开，通过加和来求指定片段的贡献
后者是把复合物轨道向片段轨道展开，直接求某个片段轨道的贡献
片段轨道是你计算单体时候得到的，其形状和复合物轨道当中直接抠出来某个片段是不同的。而且基于片段轨道讨论，才可以说片段轨道的轨道能量、不可约表示

hzfish

sobereva 发表于 2017-5-25 02:29
不是一个概念
前者是把复合物轨道向AO基函数展开，通过加和来求指定片段的贡献
后者是把复合物轨道向 ...

谢谢社长耐心回复。
以下是我们理解，社长帮着看下是否正确。

假设所有轨道均已正交归一化
复合物轨道如下：
Ψ_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

sobereva

hzfish 发表于 2017-5-25 10:06
谢谢社长耐心回复。
以下是我们理解，社长帮着看下是否正确。

你看一下CDA原文吧，这看着太费劲。CDA计算时候都是用矩阵形式来搞的，就是把复合物轨道从AO基变换成片段轨道基，通过片段轨道的系数矩阵很容易实现

hzfish

sobereva 发表于 2017-5-25 11:08
你看一下CDA原文吧，这看着太费劲。CDA计算时候都是用矩阵形式来搞的，就是把复合物轨道从AO基变换成片 ...

谢谢！
我看下CDA原文。
有问题再请教。