24小时热门版块排行榜    

查看: 1000  |  回复: 0

zhl2967

铜虫 (小有名气)

[交流] 【求助】单斜结构的对称性

我计算一个单斜结构,原胞28个原子,symmetry operation应该是4个,现在总是2个,不知怎么回事,请高手们指教。
下面是输入卡:

&control
    calculation='vc-relax',
    restart_mode='from_scratch',
    pseudo_dir='/home/abin/work/pwscf/pseudo',
    outdir='/home/abin/work/pwscf/tmp/1',
    prefix='Mg2NiH4',
    tprnfor = .true.
    tstress = .true

   

/
&system   
    ibrav = 13,
    celldm(1) = 27.104339854187,
    celldm(2) =  0.446475632713,
    celldm(3) =  0.451997490065,
    celldm(4) = -0.399069158980,
     nat=28, ntyp=3,
     ecutwfc = 66.0,
     
   
  
/
&electrons

    diagonalization='cg'
    conv_thr = 1.0e-8
    mixing_beta = 0.7
   
/
&IONS
/
&CELL
   cell_dynamics = 'bfgs' ,
   press = 0.00 ,
/
ATOMIC_SPECIES
Mg  24.305  Mg.pz-bhs.UPF
Ni  58.6934 Ni.pz-nd-rrkjus.UPF
H   1.0079  H.pz-van_ak.UPF

ATOMIC_POSITIONS {crystal}

Mg  0.217514  0.252066  0.924578
Mg  0.014433  0.985567  0.750000
Mg  0.513031  0.486969  0.750000
Mg  0.252066  0.217514  0.424578
Mg  0.747934  0.782486  0.575422
Mg  0.782486  0.747934  0.075422
Mg  0.486969  0.513031  0.250000
Mg  0.985567  0.014433  0.250000
Ni  0.111351  0.649729  0.916788
Ni  0.350271  0.888649  0.583212
Ni  0.888649  0.350271  0.083212
Ni  0.649729  0.111351  0.416788
H   0.125635  0.135589  0.581523
H   0.088212  0.489160  0.696272
H   0.547670  0.180306  0.618877
H   0.489160  0.088212  0.196272
H   0.135589  0.125635  0.081523
H   0.276308  0.702670  0.946278
H   0.819694  0.452330  0.881123
H   0.702670  0.276308  0.446278
H   0.297330  0.723692  0.553722
H   0.180306  0.547670  0.118877
H   0.723692  0.297330  0.053722
H   0.864411  0.874365  0.918477
H   0.510840  0.911788  0.803728
H   0.874365  0.864411  0.418477
H   0.911788  0.510840  0.303728
H   0.452330  0.819694  0.381123

K_POINTS {automatic}

  4 4 4 0 0 0
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 zhl2967 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见