24小时热门版块排行榜    

查看: 1447  |  回复: 4

ddang100

银虫 (小有名气)

[求助] siesta计算磁性的问题

现在计算一下关于磁性的问题,如何设置才能看到每个原子的磁矩?我设置WriteMullikenPop  1 也只能得到total spin?另外如何设置反铁磁性的情况,是不是可以同过block设置两个磁性原子+,-就可以?谢谢!在网上也没有搜到相关的回答,请有经验的给指点一下。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhsh

木虫 (小有名气)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教! 2011-10-12 08:47:01
要得到每个原子的磁矩,把每个原子的上自旋和下自旋的电子数减一下就可以了。当然一个一个减会比较麻烦,特别是当原子数目很多的时候。这时候,可以自己写一个shell脚本,直接得到每个原子的磁矩。铁磁或者反反铁都可以通过block设置。

» 本帖已获得的红花(最新10朵)

2楼2011-10-11 16:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
贺仪(金币+5): 多谢指教! 2011-10-12 08:47:16
比如对于你的体系,第4,5,6,7号原子是磁性原子,磁矩分别是2,2,3,3;进行铁磁性设置,考虑共线磁性
%block DM.InitSpin
4 +2.0 0.0 0.0
5 +2.0 0.0 0.0
6 +3.0 0.0 0.0
7 +3.0 0.0 0.0
%endblock DM.InitSpin
如果进行亚铁磁设置
%block DM.InitSpin
4 -2.0 0.0 0.0
5 -2.0 0.0 0.0
6 +3.0 0.0 0.0
7 +3.0 0.0 0.0
%endblock DM.InitSpin

» 本帖已获得的红花(最新10朵)

3楼2011-10-11 16:38:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)

送鲜花一朵
引用回帖:
3楼: Originally posted by identation at 2011-10-11 16:38:53:
比如对于你的体系,第4,5,6,7号原子是磁性原子,磁矩分别是2,2,3,3;进行铁磁性设置,考虑共线磁性
%block DM.InitSpin
4 +2.0 0.0 0.0
5 +2.0 0.0 0.0
6 +3.0 0.0 0.0
7 +3.0 0.0 0.0
%endblock DM.InitSp ...

非常感谢
4楼2011-10-12 17:01:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by xhsh at 2011-10-11 16:26:59:
要得到每个原子的磁矩,把每个原子的上自旋和下自旋的电子数减一下就可以了。当然一个一个减会比较麻烦,特别是当原子数目很多的时候。这时候,可以自己写一个shell脚本,直接得到每个原子的磁矩。铁磁或者反反铁 ...

非常感谢
5楼2011-10-12 17:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ddang100 的主题更新
信息提示
请填处理意见