简单Multiwfn使用问题和波函数分析问题答疑专贴

Jameslzy

sobereva 发表于 2019-2-26 00:31
绘制轨道相互作用图的界面里有对某个片段的能级进行整体平移的功能，设置移动参数让能级接近即可
Multiw ...

Sob老师好，我按照您说的操作了一下，确实解决了很多问题，谢谢老师的指导。现在我想使用CDA功能时想用里面的子功能2，看一下片段1的某条轨道对配合物中所有轨道的贡献程度。但是因为涉及轨道特别多，屏幕上只能显示出来最后的一小部分，里面没有想要看的数据。想问下Sob老师，这部分功能的完整输出结果可以输出出来看到吗？谢谢老师。

sobereva

Jameslzy 发表于 2019-2-26 11:20
Sob老师好，我按照您说的操作了一下，确实解决了很多问题，谢谢老师的指导。现在我想使用CDA功能时想用里 ...

Multiwfn FAQ
http://sobereva.com/452

看Q25

Jameslzy

sobereva 发表于 2019-2-26 13:51
Multiwfn FAQ
http://sobereva.com/452

谢谢Sob老师回复。还想再请教您一个问题，我用CDA的子功能2看片段1的某条轨道对配合物中所有轨道的贡献，发现片段1的轨道序号只能输入到Alpha轨道的序号上限，如果输入Beta轨道序号就会显示超出范围。想问下老师是不是因为输入Alpha轨道序号时本身就会显示这条轨道对配合物的Alpha，Beta两套轨道贡献，其中对Beta的贡献对应的就是片段1的Beta轨道，所以没有必要单独将Beta轨道的序号列出来。不知道这样理解对不对，谢谢老师

sobereva

Jameslzy 发表于 2019-2-26 15:30
谢谢Sob老师回复。还想再请教您一个问题，我用CDA的子功能2看片段1的某条轨道对配合物中所有轨道的贡献， ...

输入x，就会输出alpha x和beta x复合物轨道的成份，因此试图单独选择alpha或beta是没有意义的，例如：

Input the index of complex orbital you are interested in, e.g. 6
To exit, input 0
Note: If you want to obtain the composition of an orbital (e.g. orbital 5) of a
fragment (e.g. fragment 2) in all complex orbitals, you can input e.g. 2,5
16
Note: Only the fragment orbitals with contribution >  1.0% will be shown below
                            **** Alpha orbitals ****
Occupation number of alpha orbital    16 of the complex:  0.00000000
Alpha orbital    11 of fragment  1, Occ: 0.00000    Contribution:    1.9335%
Alpha orbital    13 of fragment  1, Occ: 0.00000    Contribution:    5.9076%
Alpha orbital    14 of fragment  1, Occ: 0.00000    Contribution:   85.8532%
Alpha orbital    15 of fragment  1, Occ: 0.00000    Contribution:    1.0990%
Alpha orbital     2 of fragment  2, Occ: 0.00000    Contribution:    4.8107%
                            **** Beta orbitals ****
Occupation number of beta orbital    16 of the complex:  0.00000000
Beta orbital    10 of fragment  1, Occ: 0.00000    Contribution:    2.6735%
Beta orbital    12 of fragment  1, Occ: 0.00000    Contribution:    8.1104%
Beta orbital    13 of fragment  1, Occ: 0.00000    Contribution:    1.3983%
Beta orbital    14 of fragment  1, Occ: 0.00000    Contribution:   81.4563%
Beta orbital     2 of fragment  2, Occ: 0.00000    Contribution:    6.0164%

Jameslzy

sobereva 发表于 2019-2-26 16:38
输入x，就会输出alpha x和beta x复合物轨道的成份，因此试图单独选择alpha或beta是没有意义的，例如：

...

好的，谢谢Sob老师回复

heheda

sob老师，下面是Multiwfn计算的跃迁偶极矩，那么i=j=1计算的是不是激发态偶极矩呢？
Ground state dipole moment in X,Y,Z:   -2.401910    0.525221    0.049533 a,u,

Transition dipole moment between excited states (a.u.):
     i     j         X             Y             Z        Diff.(eV)   Oscil.str
     1     1   -10.1646347    -0.5865941     0.0635409     0.00000     0.00000
     1     2     0.7071085     0.5383246    -0.0210329     0.52850     0.01023
     1     3     0.6315688     0.0792784     0.0061669     0.67880     0.00674
     1     4    -0.1861636    -0.0269183     0.0357128     0.75420     0.00068

在gjf文件加density计算的基态和激发态偶极矩如下：
Dipole moment (field-independent basis, Debye):
    X=            -6.1051    Y=             1.3350    Z=              0.1259  Tot=             6.2506
Dipole moment (field-independent basis, Debye):
    X=            -20.2461    Y=             -0.2695    Z=              0.0007  Tot=             20.2479

两种方法基态偶极矩是一样的，但是Multiwfn的激发态偶极矩偏大。我是新手，不太懂，sob老师不要生气我问这么简单的问题。

heheda

heheda 发表于 2019-3-1 09:34
sob老师，下面是Multiwfn计算的跃迁偶极矩，那么i=j=1计算的是不是激发态偶极矩呢？
Ground state dipole ...

我找到答案了。
对于i=j的情况来说，其数值<i|-r|i>对应的是第i激发态的偶极矩的由电子贡献的部分（注意这不等于这个态的偶极矩，因为还有另一部分，即原子核电荷对偶极矩的贡献）。

sobereva

heheda 发表于 2019-3-1 09:34
sob老师，下面是Multiwfn计算的跃迁偶极矩，那么i=j=1计算的是不是激发态偶极矩呢？
Ground state dipole ...

仔细看
使用Multiwfn计算激发态间的跃迁偶极矩和各个激发态的偶极矩
http://sobereva.com/227
里面该交代的都详细交代了，也给出了和Gaussian输出的对比

Jameslzy

Sob老师好，最近一直用您Multiwfn中轨道双正交化功能来分析分子轨道，现在文章进入收尾阶段了。想问您一下关于双正交化方法是否有比较推荐的，可以用来说明双正交化方法的原理，还有就是体现双正交化方法应用的一些文献？谢谢Sob老师

sobereva

Jameslzy 发表于 2019-3-3 17:49
Sob老师好，最近一直用您Multiwfn中轨道双正交化功能来分析分子轨道，现在文章进入收尾阶段了。想问您一下 ...

没有很好的专门说这个的文献，之后我打算自行写一篇专门介绍双正交化方法的原理和应用的文章，目前直接引用Multiwfn 3.6(dev)手册3.100.12节即可

Jameslzy

sobereva 发表于 2019-3-3 23:46
没有很好的专门说这个的文献，之后我打算自行写一篇专门介绍双正交化方法的原理和应用的文章，目前直接引 ...

好的好的，谢谢Sob老师回复

erya123

计算CO+的轨道能量，计算了两次，仅仅一个钴离子，却发现在multiwfn里面显示的轨道的能量和轨道图像都不一样啊？

lufenghxu

老师，当电子是奇数的时候，α和β电子的轨道是分开计算的，现在分析HOMO-3和LUMO的相互用，但是HOMO-3有α和β两个电子的能量，文章中要列巨额HOMO-3轨道能量如何描述？

sobereva

lufenghxu 发表于 2019-3-7 02:15
老师，当电子是奇数的时候，α和β电子的轨道是分开计算的，现在分析HOMO-3和LUMO的相互用，但是HOMO-3有α ...

我不知道你具体怎么分析HOMO-3和LUMO的相互作用。根据具体情况考虑。也可能有更好的角度去说明你要说明的问题而不基于轨道这个本来就是虚构的概念

sobereva

erya123 发表于 2019-3-7 01:06
计算CO+的轨道能量，计算了两次，仅仅一个钴离子，却发现在multiwfn里面显示的轨道的能量和轨道图像都不一 ...

对比Gaussian输出文件里的能量，如果不同，说明当前两次计算恰好收敛到了不同的波函数，低的那个更合理。