24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1983  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[交流] 【求助完结】Au原子的磁性计算为0,为什么不是1呢?求大侠解答

想用spin-polarized 去计算金原子,用的模守恒赝式,configuration 5d106s1
但是为什么计算得到的“   
     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =     0.00 Bohr mag/cell

总觉得是1才对。求大侠解答
输入文件如下
“&CONTROL
  calculation    = "relax",
  prefix         = "Au-30"
  nstep          = 500,
  pseudo_dir   = "/home/Au",
  outdir       = "/scr",
/
&SYSTEM
  ibrav                  = 4,
  celldm(1)              = 33.123112,
  celldm(3)              = 1.42628937, !a=25Ang
  nat                    = 1,
  ntyp                   = 1,
  ecutwfc                = 30.0D0,
  nspin                  = 2,
  starting_magnetization(1) = 0.2,
  occupations      = "smearing",
  degauss          = 0.05,
/
&ELECTRONS
  conv_thr    = 1.0D-8,
  mixing_beta = 0.3D0,
  electron_maxstep = 400,
/
&IONS
  ds                = 1.D-8,
/
ATOMIC_SPECIES
Au 196.96655    Au-wc.UPF
ATOMIC_POSITIONS { angstrom }
Au          3.755923000      7.951333000      7.475759000
K_POINTS {gamma}



[ Last edited by minmin_0082003 on 2011-1-20 at 11:26 ]
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by goldenfisher at 2011-01-13 15:53:01:
ecutwfc太小了。你需要测试,norm-conserving需要的截断往往很大

非常感谢,我测试过的,确实比较大,我以为截断能不会影响到spin的计算结果,至少siesta的计算中是这样的,大截断能太耗时,所以想用小的截断,希望能做的快点.....
您意思是因为截断能的原因才使得spin的结果不正确吗?
3楼2011-01-13 16:27:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
★ ★
zzy870720z(金币+2):谢谢指教 2011-01-13 16:32:01
ecutwfc太小了。你需要测试,norm-conserving需要的截断往往很大
2楼2011-01-13 15:53:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
zzy870720z(金币+2):谢谢指点 2011-01-13 21:18:18
引用回帖:
Originally posted by minmin_0082003 at 2011-01-13 16:27:04:


非常感谢,我测试过的,确实比较大,我以为截断能不会影响到spin的计算结果,至少siesta的计算中是这样的,大截断能太耗时,所以想用小的截断,希望能做的快点.....
您意思是因为截断能的原因才使得spin的结 ...

很有可能的,阶段能大小直接影响到所有的结果,包括各种构型、电子结构等等,所以至少要测试到能量收敛时得到的才是可靠的。
你可以逐步加大试试,不行就用diagnization = 'cg' ,这个虽然慢,但是内存占得少。此外,你算一个单独的Au原子,个人认为最多15A的supercell足够了。
4楼2011-01-13 16:36:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)


★ ★
zzy870720z(金币+2):谢谢指教 2011-01-13 21:18:38
Au的磁矩应该很小,我用lapw算了一个10Bohr和15Bohr边长盒子里的单独Au原子,结果是 -0.2377196275E-02玻尔磁子左右。
如果这样的话,感觉赝势法算准应该比较困难。
5楼2011-01-13 18:21:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +4 江淮北月 2026-04-05 4/200 2026-04-06 22:44 by chenzhimin
[考研] 0703化学调剂 348分 +7 唉我超真没招了 2026-04-06 7/350 2026-04-06 19:50 by lijunpoly
[考研] 085410人工智能 初试316分 求调剂 +7 残星拂曙 2026-03-31 7/350 2026-04-06 10:09 by 蓝云思雨
[考研] 调剂 一志愿吉林大学357分 +5 .Starry. 2026-04-04 5/250 2026-04-06 09:28 by cql1109
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 一志愿双非085502,267分,过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊!
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 08工科275求调剂,可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
信息提示
请填处理意见