24小时热门版块排行榜    

查看: 1083  |  回复: 3

gemucai

木虫 (正式写手)

[交流] 交流一下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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

pwscf相关

» 猜你喜欢

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

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

meng303

金虫 (小有名气)


小木虫(金币+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
2楼2011-11-17 14:30:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
2楼: Originally posted by meng303 at 2011-11-17 14:30:10:
&CELL
/
这一项中不需要说明吗?例如cell_dynamics 或 cell_dofree

这是因为我发现都用默认值就行了
3楼2011-11-17 14:39:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meng303

金虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
youzhizhe(金币+2): 谢谢交流。 2011-11-17 20:27:21
测试Ecut,能量相差0.1Ry是不是太大了,一般是10mev 5mev 换算一下也要到小数点后第三位第四位吧,不清楚会不会影响你要计算的结果
4楼2011-11-17 15:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +12 晴天加油 2026-08-26 13/650 2026-08-28 11:32 by skyjudy
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见