应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 某个原子的轨道角动量怎么算？
31/1返回列表
查看: 5244  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

swuhaizhilan

金虫 (小有名气)

[求助] 某个原子的轨道角动量怎么算？ 已有1人参与

这是复旦Hua Wu老师组一篇计算磁各向异性的文章，他计算的是过渡金属吸附在MgO的O原子上的磁各向异性，用的软件是wein2k。
比如Co吸附，不加SOC的时候，按晶场作用，理论上Co原子中电子按图中的（b）分布，spin为3玻尔磁子，加了SOC之后两个简并的轨道发生劈裂（图中没有标注自旋向上的态），两个自旋向下的电子占据d1和d2，从而spin仍为3玻尔磁子，轨道角动量理论上也为3玻尔磁子，他在文章中求出了轨道角动量的值为2.48玻尔磁子。
我想知道这里的轨道角动量2.48玻尔磁子是怎么计算出的，请求知道的人帮助一下，不胜感激！某个原子的轨道角动量怎么算？
能级劈裂.PNG


某个原子的轨道角动量怎么算？-1
捕获.PNG


某个原子的轨道角动量怎么算？-2
soc+pdos.PNG
回复此楼

» 本帖附件资源列表

» 猜你喜欢
1楼 2016-04-19 15:48:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CEA_SACLAY

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
swuhaizhilan: 金币+10, ★★★很有帮助 2016-04-20 17:18:21
L = r ^ p (L是orbital angular momentum, p是momentum operator)，求轨道角动量值就是求<L>，也就是ml = sum_i <Li>=sum_n fn<psi n | Li |psi n>, ml是系统总轨道角动量，i是原子site，n是band index，fn是Fermi-dirac，如果你用DFT，那psi就是Kohn Sham orbital。
如果想得到local orbital moment，你可以用mulliken charge。
希望能帮到你。

PS：我没用过WIEN2k，但很可能算轨道角动量已经implemented。
赞一下 回复此楼
Physicsislikesex!
2楼2016-04-19 21:03:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swuhaizhilan

金虫 (小有名气)
引用回帖:
2楼: Originally posted by CEA_SACLAY at 2016-04-19 21:03:12
L = r ^ p (L是orbital angular momentum, p是momentum operator)，求轨道角动量值就是求<L>，也就是ml = sum_i <Li>=sum_n fn<psi n | Li |psi n>, ml是系统总轨道角动量，i是原子site，n是band ...

我也不用WIEN2K，我用的是VASP。
请问一下VASP能算mulliken charge吗？通过电荷的分布怎么确定轨道角动量呢？我没算过mulliken charge，还请指点一下。谢谢。
赞一下 回复此楼
3楼2016-04-20 17:18:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 swuhaizhilan 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭