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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 【求助】请教，为什么明明是铁磁性体系，abinit的计算结果没有磁矩
【调剂】北京石油化工学院2024年16个专业接受调剂
71/1返回列表
查看: 1052  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 identation 的 4 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会

identation

金虫 (著名写手)

[交流] 【求助】请教，为什么明明是铁磁性体系，abinit的计算结果没有磁矩

输入文件如下
occopt 1
nsppol 2
spinat
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 1.5
0.0 0.0 1.5
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0

toldfe   1.0d-10      
iscf   5  
diemac 8.0           

prtdos 2         

acell 5.4858 5.4858 7.7518 Angstroms
angdeg 90 90 90
spgroup 83
natrd 9
brvltt -1

ntypat 4            
znucl  
     
natom 20               
typat 1 2 2 3 3 4 4 4 4
xred   
0.00000   0.50000   0.25000
0.00000   0.00000   0.00000
0.50000   0.50000   0.50000
0.50000   0.50000   0.00000
0.00000   0.00000   0.50000
0.00000   0.00000   0.25040
0.25700   0.24290   0.00000
0.24290   0.25700   0.50000
0.50000   0.50000   0.24960

nband  60
ecut 72
ecutsm 0.5
kptopt 1
ngkpt 6 6 4
nshiftk 1
shiftk 0.5 0.5 0.5
nstep 400

然而得到的结果
Total charge density [el/Bohr^3]
,     Maximum=    9.9364E-01  at reduced coord.    0.3125    0.8125    0.5250
,     Minimum=    8.8712E-06  at reduced coord.    0.5000    0.0000    0.2500
Spin up density      [el/Bohr^3]
,     Maximum=    5.2799E-01  at reduced coord.    0.4875    0.5625    0.9667
,     Minimum=    4.4354E-06  at reduced coord.    0.5000    0.0000    0.2500
Spin down density    [el/Bohr^3]
,     Maximum=    5.2799E-01  at reduced coord.    0.9375    0.9875    0.5333
,     Minimum=    4.4358E-06  at reduced coord.    0.5000    0.0000    0.2500
Magnetization (spin up - spin down) [el/Bohr^3]
,     Maximum=    5.1055E-01  at reduced coord.    0.4875    0.5625    0.9667
,     Minimum=   -5.1055E-01  at reduced coord.    0.0625    0.0125    0.4667
Relative magnetization (=zeta, between -1 and 1)
,     Maximum=    9.3606E-01  at reduced coord.    0.4875    0.5625    0.9667
,     Minimum=   -9.3606E-01  at reduced coord.    0.0625    0.0125    0.4667

体系的DOS自旋上和自旋下的电子密度基本一样
用的是FHI的GGA赝势，IXC=11
请问这是为什么呢？很奇怪呀，之前用LDA赝势abinit计算磁性把绝缘体计算成金属
还想问一下，ABINIT的PAW方法能否计算介电极化呢？
谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-10-25 15:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★
zzy870720z(金币+1):谢谢回复 2010-10-25 18:32:39
identation(金币+5): 2010-10-25 21:15:39
没有用过abinit计算自旋极化，不太清楚。
其实状态是铁磁的吗？如果起始状态是非极化的，那么可能弛豫之后还是非极化的
赞一下(1人) 回复此楼
2楼2010-10-25 16:17:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★
ym23(金币+1):感谢回答 2010-10-25 21:06:55
identation(金币+10): 2010-10-25 21:15:43
可能是用的赝势有问题。另外可以考虑将初始的磁矩设置的大些。abinit提供的FHI的pp的质量并不高。用之前，最好对单质的晶体性质进行检查。abinit里的PAW好像也是可以计算介电性质。在它的手册上找找看。
赞一下(1人) 回复此楼
3楼2010-10-25 19:58:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主
★ ★
identation(金币+10): 2010-10-26 13:24:28
cenwanglai(金币+2):谢谢回复！ 2010-10-28 19:38:06
(1)spinat的单位的是 hba/2，你这个地方的磁性原子可能是Mn，1.5太小了，可以设置大一些。
(2)建议使用ABINIT的PAW赝势，可以同时做自旋极化计算。
赞一下(1人) 回复此楼
4楼2010-10-26 08:13:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangkunspace

铜虫 (小有名气)
★
identation(金币+1):谢谢参与
没有一样吧  你看总的额 是不相等啊
赞一下(1人) 回复此楼
5楼2010-10-28 17:06:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bestone

木虫 (正式写手)
★
identation(金币+1):谢谢参与
自旋极化不太清楚 好复杂
赞一下(1人) 回复此楼
6楼2010-10-28 17:26:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenjlu

木虫 (正式写手)
★ ★
identation(金币+1):谢谢参与
zxzj05(金币+1):谢谢回帖交流 2010-10-29 18:29:02
查看磁矩，在out文件中，找关键词 magnetization。

此外，你的input输入的晶格参数，单位 Angstroms ，错误。 abinit不认识。
赞一下(2人) 回复此楼
7楼2010-10-29 12:42:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
71/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 兄弟们帮我看看这两个结构怎么合成 +5 xl2088131 2024-04-17 6/300 2024-04-20 12:25 by if今年的吧
[教师之家] 发现督导在听课，需要跟TA招呼示意吗？有同学看手机、课件有疏漏算大问题吗？ +6 河西夜郎 2024-04-19 6/300 2024-04-20 09:31 by xli1984
[基金申请] 颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求 +3 zju2000 2024-04-18 3/150 2024-04-20 00:48 by cococat951
[基金申请] 申请省自然科学基金，研究区能否是省外区域 100+3 喜欢兔兔的我 2024-04-15 13/650 2024-04-19 21:49 by 喜欢兔兔的我
[考博] 24/25读博求博导 +4 宝23 2024-04-16 4/200 2024-04-19 17:18 by wangzhe_bs
[教师之家] 美国教授开除中国研究生：我就不该录取你！ +7 苏东坡二世 2024-04-15 7/350 2024-04-19 16:20 by feng6531
[论文投稿] 求助 100+3 guhangyu 2024-04-15 3/150 2024-04-19 08:48 by bobvan
[考研] 267求调剂 +3 工科材料267 2024-04-15 3/150 2024-04-18 19:45 by lature00
[论文投稿] 投稿求助 5+3 我是洲洲啊 2024-04-17 5/250 2024-04-18 17:13 by topedit
[博后之家] 博后换方向可行吗？ +3 越越不暴躁 2024-04-15 3/150 2024-04-18 10:58 by ciompman
[考研] 一志愿北化085600 328求线上调剂 +10 Mrshark 2024-04-13 10/500 2024-04-18 08:30 by 刘国宁
[论文投稿] MDPI初审7天，是不是被拒 20+6 applegirl5762 2024-04-15 8/400 2024-04-17 08:49 by bobvan
[基金申请] 请问教育部人文社科难度大吗 +9 锦衣卫寒战 2024-04-15 14/700 2024-04-16 23:10 by 锦衣卫寒战
[考博] 25年申博求助，射频方向 +4 SherlockAH 2024-04-15 10/500 2024-04-16 18:06 by ljysudaee
[考研] 求调剂 +4 鹿萌月yy 2024-04-15 5/250 2024-04-16 16:16 by chenweiwade
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[考研] 求调剂 +4 桃岸雪 2024-04-15 5/250 2024-04-15 18:49 by mthwyj
[考研] 347求调剂 +4 安谷先生 2024-04-15 4/200 2024-04-15 18:27 by huao连
[考研] 300求调剂 +5 @Xqy668800 2024-04-14 6/300 2024-04-14 23:03 by lincunhui
[考研] 290求调剂 +3 杨yhr 2024-04-14 5/250 2024-04-14 21:50 by coco1981
信息提示
关闭
请填处理意见
关闭