请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 292|回复: 6

[Multiwfn使用咨询] 轨道成份分析咨询

[复制链接]

81

帖子

0

威望

930

eV
积分
1011

Level 4 (黑子)

发表于 2017-5-23 18:55:29 | 显示全部楼层 |阅读模式
A(B)nn=6)分子进行轨道成份分析,假设A为一个金属阳离子,B为水分子。
问题:
1. 进行片断分析,A和每个B为一个片断,要得到第i轨道每个片断的贡献,得分析1+n
2. 假设第i轨道中,第2B片断贡献为a%,而独立B轨道为1a12a11b23a11b1,那如何判断第2B片断中哪个轨道对第i轨道产生贡献a%?即产生贡献是1a12a11b23a11b1中的哪个?
3. 假设第i轨道中,第2B片断贡献为a%,那么第i轨道中第2B片断贡献部分能不能单独图形化?

谢谢!


1万

帖子

16

威望

1万

eV
积分
24760

管理员

公社社长

发表于 2017-5-24 02:41:50 | 显示全部楼层
1 是的
2 貌似你要考察片段轨道对复合物轨道的贡献,这得用CDA模块才行,可以输出复合物轨道当中各个片段轨道的成份。见
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166
3 最简单的做法,在主界面先选-3 Obtain a fragment contribution to property,然后输入要考察的那个片段的原子序号,之后你在功能0里看到的轨道图形就只保留纯粹属于那个片段的部分了。此文也有提到
使用Multiwfn观看分子轨道
http://sobereva.com/269
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、Multiwfn程序培训
思想家公社的门口Blog:http://sobereva.com
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

81

帖子

0

威望

930

eV
积分
1011

Level 4 (黑子)

 楼主| 发表于 2017-5-24 14:33:01 | 显示全部楼层
sobereva 发表于 2017-5-24 02:41
1 是的
2 貌似你要考察片段轨道对复合物轨道的贡献,这得用CDA模块才行,可以输出复合物轨道当中各个片段 ...

谢谢社长解答。

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

1万

帖子

16

威望

1万

eV
积分
24760

管理员

公社社长

发表于 2017-5-25 02:29:11 | 显示全部楼层
hzfish 发表于 2017-5-24 14:33
谢谢社长解答。

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


不是一个概念
前者是把复合物轨道向AO基函数展开,通过加和来求指定片段的贡献
后者是把复合物轨道向片段轨道展开,直接求某个片段轨道的贡献
片段轨道是你计算单体时候得到的,其形状和复合物轨道当中直接抠出来某个片段是不同的。而且基于片段轨道讨论,才可以说片段轨道的轨道能量、不可约表示
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、Multiwfn程序培训
思想家公社的门口Blog:http://sobereva.com
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

81

帖子

0

威望

930

eV
积分
1011

Level 4 (黑子)

 楼主| 发表于 2017-5-25 10:06:17 | 显示全部楼层
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

1万

帖子

16

威望

1万

eV
积分
24760

管理员

公社社长

发表于 2017-5-25 11:08:39 | 显示全部楼层
hzfish 发表于 2017-5-25 10:06
谢谢社长耐心回复。
以下是我们理解,社长帮着看下是否正确。


你看一下CDA原文吧,这看着太费劲。CDA计算时候都是用矩阵形式来搞的,就是把复合物轨道从AO基变换成片段轨道基,通过片段轨道的系数矩阵很容易实现
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、Multiwfn程序培训
思想家公社的门口Blog:http://sobereva.com
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

81

帖子

0

威望

930

eV
积分
1011

Level 4 (黑子)

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

谢谢!
我看下CDA原文。
有问题再请教。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

Archiver|手机版|小黑屋|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2017-10-24 11:50 , Processed in 0.090180 second(s), 21 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表