24小时热门版块排行榜    

查看: 1744  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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方法能否计算介电极化呢?
谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangkunspace

铜虫 (小有名气)

★
identation(金币+1):谢谢参与
没有一样吧  你看总的额 是不相等啊
5楼2010-10-28 17:06:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

bingmou

金虫 (著名写手)

★
zzy870720z(金币+1):谢谢回复 2010-10-25 18:32:39
identation(金币+5): 2010-10-25 21:15:39
没有用过abinit计算自旋极化,不太清楚。
其实状态是铁磁的吗?如果起始状态是非极化的,那么可能弛豫之后还是非极化的
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好像也是可以计算介电性质。在它的手册上找找看。
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赝势,可以同时做自旋极化计算。
4楼2010-10-26 08:13:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见