24小时热门版块排行榜

返回列表

gemucai

木虫 (正式写手)

应助: 16 (小学生)
金币: 1774.6
散金: 157
红花: 4
帖子: 693
在线: 201.7小时
虫号: 816478
注册: 2009-07-27
专业: 凝聚态物性 II ：电子结构

[交流] 交流一下vc-relax过程中出现的电荷数错误的问题已有1人参与

最近做vc-relax出错，报错如下：
from electrons : error #       1
   charge is wrong

WARNING: integrated charge= 25.58637298, expected= 52.00000000

电荷量只有预期的一半。我查了PW论坛以前出现的相似问题，他们的电荷跟预期相差无几，我这个就差的太多了。

给出输入文件：
&control
calculation='vc-relax'
restart_mode='from_scratch',
prefix='ret',
wf_collect=.true.
tstress=.true.
tprnfor=.true.
pseudo_dir = '...',
outdir='...'
/
&system
ibrav=  0, a=1, nat= 6, ntyp= 2,
ecutwfc =50.0, nspin = 2, starting_magnetization=0.5,
occupations='smearing', smearing='gaussian', degauss=0.02
/
&electrons
mixing_beta = 0.3
/
&IONS
  bfgs_ndim       = 3,
  pot_extrapolation = "second_order",
  wfc_extrapolation = "second_order",
/
&CELL
/
ATOMIC_SPECIES
Mo 95.94 Mo.NCPP_LDA.UPF
S  32.0655  S.pz-bhs.UPF
ATOMIC_POSITIONS crystal
Mo  0.6666666666666666  0.6666666666666666  0.0000000000000000
Mo  0.3333333333333333  0.3333333333333333  0.5000000000000000
S 0.3333333333333333  0.3333333333333333  0.1279036059999967
S 0.3333333333333333  0.3333333333333333  0.8720963940000033
S 0.6666666666666666  0.6666666666666666  0.6279036059999967
S 0.6666666666666666  0.6666666666666666  0.3720963940000033
K_POINTS {automatic}
9 9 5  0 0 0
CELL_PARAMETERS { cubic }
2.7366402800000000 -1.5800000000000001 0.0000000000000000
2.7366402800000000 1.5800000000000001 0.0000000000000000
0.0000000000000000 0.0000000000000000 12.2899999999999991

我赝势的阶段半径是不是太小了，不过我做了测试，在50的时候就与45时的能量相差不到0.1Ry

[ Last edited by gemucai on 2011-11-17 at 14:24 ]

回复此楼

1楼 2011-11-17 11:03:56

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

meng303

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1532.7
帖子: 84
在线: 214.2小时
虫号: 1092635
注册: 2010-09-08

★
小木虫(金币+0.5):给个红包，谢谢回帖

引用回帖:

1楼: Originally posted by gemucai at 2011-11-17 11:03:56:
最近做vc-relax出错，报错如下：
from electrons : error # 1
charge is wrong

WARNING: integrated charge= 25.58637298, expected= 52.00000000

电荷量只有预期的一半。我查了 ...

&CELL
/
这一项中不需要说明吗？例如cell_dynamics 或 cell_dofree