24小时热门版块排行榜    

查看: 530  |  回复: 2

子虚乌有5388

木虫 (正式写手)

[求助] 锯齿石墨烯纳米带计算没有磁性??

大家好!
我在试着计算zigzag-6石墨烯纳米带的磁性时,在relax和scf的结果里,都显示absolute磁矩为零,这是怎么回事呢?
输入如下:
&SYSTEM
      ibrav=8,
      celldm(1)=28.355,
      celldm(2)=2.133,
      celldm(3)=0.1607,
      nat=14,
      ntyp=4,
      nspin=2,
      ecutwfc=45.0,
      occupations='smearing',
      smearing='gauss',
      degauss=0.02,
      starting_magnetization(1)=0.5,
      starting_magnetization(2)=-0.5,
/
&ELECTRONS
    electron_maxstep = 1000,
    mixing_beta = 0.3,
    conv_thr =  1.0d-8,
/
&ions
/
ATOMIC_SPECIES
C1   12.011    C.pz-vbc.UPF
C2   12.011    C.pz-vbc.UPF
C    12.011    C.pz-vbc.UPF
H    1.008     H.pz-vbc.UPF
ATOMIC_POSITIONS (crystal)
C1       0.487005440   0.666688115   0.822252166
C2       0.487298760   0.312952278   0.821548039
C        0.487113582   0.534102074   0.821980559
C        0.487075914   0.578638355   0.822116343
C        0.487059909   0.600605565   0.322153627
C        0.487024733   0.645532001   0.322229975
C        0.487225142   0.400989425   0.821609578
C        0.487188953   0.445523911   0.821727786
C        0.487171450   0.467570471   0.321772706
C        0.487133573   0.512053646   0.321933849
C        0.487283278   0.334104014   0.321548976
C        0.487244537   0.379026089   0.321582090
H        0.486959303   0.700973118   0.822388785
H        0.487315425   0.278670936   0.821645522
K_POINTS (automatic)
1 1 18  0 0 0

我把2个边缘的C原子加了初始磁矩。

我的问题出在哪里呢?
谢谢大家指点。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

房无一间,地无一亩,身无分文,胸无大志。
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

浪子必回头

新虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
identation: 金币+1, 谢谢回帖交流! 2013-09-09 14:31:48
你这个算的纳米带,我觉得应该是左右对称导致磁矩抵消了~
2楼2013-09-09 12:27:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

子虚乌有5388

木虫 (正式写手)

nobody knows?
房无一间,地无一亩,身无分文,胸无大志。
3楼2013-09-09 16:22:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 子虚乌有5388 的主题更新
信息提示
请填处理意见