计算化学公社

标题: 如何根据某原子周围的积分自旋密度差计算该原子的价态？ [打印本页]

作者
Author: lisu555 时间: 2017-3-10 17:37
标题: 如何根据某原子周围的积分自旋密度差计算该原子的价态？
如何根据某原子周围的积分自旋密度差计算该原子的价态？该原子的积分半径如何确定呢？谢谢！

作者
Author: lisu555 时间: 2017-3-10 17:39
对理想的二价铁，它的电子占据为 t2g3↑eg2↑t2g1↓

作者
Author: 卡开发发 时间: 2017-3-12 14:03
本帖最后由卡开发发于 2017-3-12 14:13 编辑

我猜测文章里面是通过不断修改Fe原子的RWIGS来实现的，这个可以通过写个简单的Bash脚本实现，每一次进行计算的时候只需要非自洽就可以，要不了多少时间就能算出来。不过这样的做法的科学性就很难保证了，Nup-Ndown是局域磁矩，但是不见得局域磁矩与价态就一一对应的关系，这样乱七八糟的东西还是少看为妙。

注：VASP自带的LORBIT本来就是通过划分Wigner-Seiz球并对球内的电荷进行积分来实现的。
还有，图横坐标的Anstron是什么？

作者
Author: sobereva 时间: 2017-3-12 15:27

卡开发发发表于 2017-3-12 14:03
我猜测文章里面是通过不断修改Fe原子的RWIGS来实现的，这个可以通过写个简单的Bash脚本实现，每一次进行计 ...

应该是指angstrom

作者
Author: lisu555 时间: 2017-3-23 20:19
补充：VASP 怎么做积分自旋得到某个磁性原子的价态？另外这个积分半径选择原则是什么呢？如图，附参考文献！

作者
Author: lisu555 时间: 2017-3-23 20:21

卡开发发发表于 2017-3-12 14:03
我猜测文章里面是通过不断修改Fe原子的RWIGS来实现的，这个可以通过写个简单的Bash脚本实现，每一次进行计 ...

可是，根据自旋积分求原子价态的文献还是很多啊！如我补充的附件文献，您能帮看看吗？谢谢！

作者
Author: 卡开发发 时间: 2017-3-23 22:35

lisu555 发表于 2017-3-23 20:19
补充：VASP 怎么做积分自旋得到某个磁性原子的价态？另外这个积分半径选择原则是什么呢？如图，附参考文献 ...

原则上应该是随着WS半径增加电荷收敛且包含的总电荷总磁矩应当和实际SCF的情况尽可能一致。

如果不同的价态磁矩不同，用磁矩区分价态当然没什么不可。但是如果说直接用磁矩定量价态，恐怕还是不行的，比如Co3O4那个系统我也没见Co的二价和三价与磁矩呈现定量关系。

作者
Author: lisu555 时间: 2017-3-24 16:36

卡开发发发表于 2017-3-23 22:35
原则上应该是随着WS半径增加电荷收敛且包含的总电荷总磁矩应当和实际SCF的情况尽可能一致。

如果不同 ...

您的意思是：上图中如 V 0到3Å 的积分半径，是通过通过不断修改V原子的RWIGS来实现的吗？比如取 0 0.2 0.5 0.8 1.0 1.5 2.0 2.5 3.0Å 求出对应的磁矩再画图把这些点连成像图中的线条？可是图中好像应该是连续积分得的磁矩值呢！！

作者
Author: 卡开发发 时间: 2017-3-24 22:15

lisu555 发表于 2017-3-24 16:36
您的意思是：上图中如 V 0到3Å 的积分半径，是通过通过不断修改V原子的RWIGS来实现的吗？比如取 0 ...

是啊，取点密集写看起来就比较连续了。

作者
Author: lisu555 时间: 2017-3-24 22:27

卡开发发发表于 2017-3-24 22:15
是啊，取点密集写看起来就比较连续了。

哦，谢谢您，我先试试！

作者
Author: aais-chendc 时间: 2021-9-13 22:09
做出了类似您贴的图的数据，遇到，相关问题，不知怎么理解Fe1.8+ 去对应纵坐标的数目。可否指教

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