计算化学公社

标题: 能量分解之困惑（EDA） [打印本页]

作者
Author: jctyp 时间: 2019-7-29 17:56
标题: 能量分解之困惑（EDA）
大家好：计算一个Cr的配位化合物，优化结构自旋多重度为1，单重态。NBO所得中心原子Cr电子组态为dz2(1.4225),dx2-y2(0.5822）,dxy (1.7294),dxz(1.0303),dyz(0.9187)。ADNDP分析：单中心搜索Cr的dz2占据1.89个电子，dx2-y2和dxy 分别是1.30和1.29.想请教以下几个问题：1 NBO和ADNDP结果为什么不一样？
2做能量分解，分成两个碎片：中心原子和配体，请问怎么指定中心原子的状态：
(1) 碎片中心原子的自旋多重度应该1还是或者3,5等更高的多重度？
(2) 如果是1，五个d轨道分别的电子占据是多少？如果与上述NBO结果一致，是不是应该是dz2(1),dx2-y2(1),dxy(2)dxz(1)dyz(1)？如果要与ADNDP结果一致，dz2(2),dx2-y2(2),dxy(2)这样指定碎片中心原子的占据合适吗？

作者
Author: granvia 时间: 2019-7-29 20:04
笼统说dz2占据数没有意义，应该指出是哪个定域轨道中的dz2分量上的占据数。NBO找的定域轨道与adndp找到不一样。关于EDA，这里金属与配体形成的是强化学键，我不太推荐做EDA

作者
Author: sobereva 时间: 2019-7-29 22:16
对于你当前的数据，NBO程序做的是NPA分析，算法和AdNDP极为不同，没有可比性。对于考察单个原子组态问题，用NPA就行了，涉及到讨论成键轨道才需要用AdNDP。

看具体是什么体系，金属的自旋情况和配体有关。能量分解的时候一般取尽可能接近片段在当前体系内实际状态的状态

作者
Author: granvia 时间: 2019-7-30 12:02

sobereva 发表于 2019-7-29 22:16
对于你当前的数据，NBO程序做的是NPA分析，算法和AdNDP极为不同，没有可比性。对于考察单个原子组态问题， ...

NPA基于的是NAO, 而adNDP基于的也是NAO吧？

作者
Author: granvia 时间: 2019-7-30 12:05

sobereva 发表于 2019-7-29 22:16
对于你当前的数据，NBO程序做的是NPA分析，算法和AdNDP极为不同，没有可比性。对于考察单个原子组态问题， ...

“能量分解的时候一般取尽可能接近片段在当前体系内实际状态的状态”。我觉得这个不太好操作。金属的价态（氧化数）往往与实际原子电荷差别非常大

作者
Author: jctyp 时间: 2019-7-30 19:16

granvia 发表于 2019-7-29 20:04
笼统说dz2占据数没有意义，应该指出是哪个定域轨道中的dz2分量上的占据数。NBO找的定域轨道与adndp找到不一 ...

谢谢！

作者
Author: jctyp 时间: 2019-7-30 19:18

granvia 发表于 2019-7-30 12:05
“能量分解的时候一般取尽可能接近片段在当前体系内实际状态的状态”。我觉得这个不太好操作。金属的价 ...

确实不好操作，看有的文献直接给定了碎片的占据情况，没有说出所以然，所以较为困惑。

作者
Author: jctyp 时间: 2019-7-30 19:25

sobereva 发表于 2019-7-29 22:16
对于你当前的数据，NBO程序做的是NPA分析，算法和AdNDP极为不同，没有可比性。对于考察单个原子组态问题， ...

谢谢卢老师!是不是可以这样理解，不能根据NPA分析来决定中心原子碎片的实际状态。而是需要看配体是强场还是弱场，再看整个分子处在什么点群下，去判断中心原子的五个d轨道各自的不可约表示？

作者
Author: sobereva 时间: 2019-7-30 19:25

granvia 发表于 2019-7-30 12:05
“能量分解的时候一般取尽可能接近片段在当前体系内实际状态的状态”。我觉得这个不太好操作。金属的价 ...

可以计算氧化态
使用Multiwfn通过LOBA方法计算氧化态
http://sobereva.com/362（http://bbs.keinsci.com/thread-5118-1-1.html）

作者
Author: sobereva 时间: 2019-7-30 19:26

granvia 发表于 2019-7-30 12:02
NPA基于的是NAO, 而adNDP基于的也是NAO吧？

是
AdNDP即便是对于单中心搜索，也相当于重新又把NAO做了组合，对于考察单个原子的组态，这多此一举，结果也无意义。

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