版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] "starting and expected charges differ"有虫友碰到过这个错误吗？已有1人参与

如题。

最近在算一个奇数电子体系。scf计算没问题，nscf计算一直报错。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   task #       6
   from potinit : error #       1
   starting and expected charges differ
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

不知哪位高人知道如何解决？

附scf和nscf计算的输入文件。

scf.in
----------------------------------------------------------------------------------
&CONTROL
  prefix          = 'pristine'
  calculation    = 'scf'
  restart_mode    = 'from_scratch'
  pseudo_dir    = './'
  outdir          = './'
  wfcdir          = './'
  verbosity       = 'high'
  wf_collect    = .TRUE.
/
&SYSTEM
  ibrav          = 0
  celldm(1)       = 1.889726125
  ntyp          = 3
  nat             = 16
  ecutwfc       = 60
  force_symmorphic = .TRUE.

occupations    = 'smearing'
  smearing       = 'gaussian'
  degauss       = 0.02
  nbnd          = 36

  nspin          = 2
  starting_magnetization(1) = 0.0
  starting_magnetization(2) = 0.0
  starting_magnetization(3) = 0.3
/
&ELECTRONS
  electron_maxstep = 250
  conv_thr       = 1.0d-10
  mixing_mode    = 'plain'
  mixing_beta    = 0.7
  mixing_ndim    = 8
  diagonalization  = 'david'
  diago_david_ndim = 4
  diago_full_acc = .TRUE.
/
ATOMIC_SPECIES
   O 15.999  O.pbe-mt_fhi.UPF
   Mg 24.3 Mg.pbe-mt_fhi.UPF
   Sc  44.956  Sc.pbe-mt_fhi.UPF
CELL_PARAMETERS {alat}
4.305914683 4.305914683 0.000000000
0.000000000 4.305914683 4.305914683
4.305914683 0.000000000 4.305914683
ATOMIC_POSITIONS {crystal}
O       0.250000000 0.250000000 0.250000000
O       0.740711820 0.259288180 0.259288180
O       0.259288180 0.740711820 0.259288180
O       0.740711820 0.740711820 0.259288180
O       0.259288180 0.259288180 0.740711820
O       0.740711820 0.259288180 0.740711820
O       0.259288180 0.740711820 0.740711820
O       0.750000000 0.750000000 0.750000000
Mg    0.500000000 0.000000000 0.000000000
Mg    0.000000000 0.500000000 0.000000000
Mg    0.500000000 0.500000000 0.000000000
Mg    0.000000000 0.000000000 0.500000000
Mg    0.500000000 0.000000000 0.500000000
Mg    0.000000000 0.500000000 0.500000000
Mg    0.500000000 0.500000000 0.500000000
Sc    0.000000000 0.000000000 0.000000000
K_POINTS {automatic}
6 6 6 0 0 0

nscf.in
----------------------------------------------------------------------
&CONTROL
  prefix          = 'pristine'
  calculation    = 'nscf'
  restart_mode    = 'from_scratch'
  pseudo_dir    = './'
  outdir          = './'
  wfcdir          = './'
  verbosity       = 'high'
  wf_collect    = .TRUE.
/
&SYSTEM
  ibrav          = 0
  celldm(1)       = 1.889726125
  ntyp          = 3
  nat             = 16
  ecutwfc       = 60
  force_symmorphic = .TRUE.
occupations    = 'smearing'
  smearing       = 'gaussian'
  degauss       = 0.02
  nbnd          = 480

  nspin          = 2
/
&ELECTRONS
  electron_maxstep = 250
  conv_thr       = 1.0d-10
  mixing_mode    = 'plain'
  mixing_beta    = 0.7
  mixing_ndim    = 8
  diagonalization  = 'cg'
  diago_david_ndim = 4
  diago_full_acc = .TRUE.
/
ATOMIC_SPECIES
Mg  24.3 Mg.pbe-mt_fhi.UPF
O 15.999  O.pbe-mt_fhi.UPF
Sc  44.956  Sc.pbe-mt_fhi.UPF
CELL_PARAMETERS {alat}
4.305914683 4.305914683 0.000000000
0.000000000 4.305914683 4.305914683
4.305914683 0.000000000 4.305914683
ATOMIC_POSITIONS {crystal}
O       0.250000000 0.250000000 0.250000000
O       0.740711820 0.259288180 0.259288180
O       0.259288180 0.740711820 0.259288180
O       0.740711820 0.740711820 0.259288180
O       0.259288180 0.259288180 0.740711820
O       0.740711820 0.259288180 0.740711820
O       0.259288180 0.740711820 0.740711820
O       0.750000000 0.750000000 0.750000000
Mg    0.500000000 0.000000000 0.000000000
Mg    0.000000000 0.500000000 0.000000000
Mg    0.500000000 0.500000000 0.000000000
Mg    0.000000000 0.000000000 0.500000000
Mg    0.500000000 0.000000000 0.500000000
Mg    0.000000000 0.500000000 0.500000000
Mg    0.500000000 0.500000000 0.500000000
Sc    0.000000000 0.000000000 0.000000000
K_POINTS {automatic}
6 6 6 0 0 0

回复此楼

» 猜你喜欢

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

《Metal Matrix Composites_Custom-made Materials for Automotive & Aerospace》已经有73人回复
【英美经典书籍】《Dictionary of Engineering》(第2版)【已搜索，无重复】已经有81人回复
《Data Mining: Concepts, Models, Methods, and Algorithms》英文版&中文版已经有103人回复
《Chemical Engineering Dynamics: An Intro to Modelling&Computer Simulation》已经有51人回复
【英美经典书籍】《Chemistry from First Principles》【已搜索，无重复】已经有149人回复
《Fluidization, Solids Handling, and Processing: Industrial Applications》已经有34人回复
【英美经典书籍】《Neural Networks in Chemistry and Drug Design》(第2版) 已经有38人回复
A subtle Chinglish phrase, "Mind acts upon mind" 已经有5人回复
【英美经典】《Goldfrank's Manual of Toxicologic Emergencies (中毒急诊手册)》已经有30人回复
【课件】【英美小学课件083】《The Grass is Singing》【课本 + 教案 + 课堂活动】已经有26人回复
【英美高中生物教材002】《Biology: An Illustrated Guide to Science》【无重复】已经有247人回复
【2014年将要出版的书籍】《Quantum Theory: A Two-Time Success Story》已经有321人回复
【职场英语技巧19】求职英语：五种技巧"摆平"面试官已经有0人回复
【手机阅读小说】《Complete Ghost Stories (鬼故事大全)》【已搜索，无重复】已经有12人回复
给<外文书籍求助板块>版主".net", 提供的证据已经有4人回复
From FACEBOOK--"Interesting Facts Around The World" O(∩_∩)O 已经有21人回复
SCOPUS2011寻找青年科学之星"第二轮评审入围名单已经有1人回复
"My name is Khan" 已经有6人回复
求助“校样"时"AUTHOR QUERY FORM"中一句话的意思已经有6人回复

1楼 2014-11-07 14:43:25

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

qingeraaa

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 10.7
散金: 18
红花: 5
帖子: 219
在线: 195.9小时
虫号: 1503611
注册: 2011-11-21
性别: MM
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ljw4010: 金币+1, 谢谢指导解答！ 2014-11-08 10:02:34
liqizuiyang: 金币+30, ★有帮助 2014-11-10 14:52:00

看一下：http://trac-foundry.lbl.gov/trac ... potinit.f90?rev=527
你的问题出现在第184行，是由于电荷数设置问题，如果设置tot_charge应该有帮助。

赞一下

回复此楼

对于未来，我一点也不不担心，因为我相信，时光会把我变得更好！

3楼2014-11-07 21:41:33

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

查看全部 9 个回答

qingeraaa

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 10.7
散金: 18
红花: 5
帖子: 219
在线: 195.9小时
虫号: 1503611
注册: 2011-11-21
性别: MM
专业: 凝聚态物性 II ：电子结构

★
ljw4010: 金币+1, 谢谢交流！ 2014-11-08 10:02:16

是不是nbnd设置太高啊？

赞一下

回复此楼

对于未来，我一点也不不担心，因为我相信，时光会把我变得更好！

2楼2014-11-07 21:36:05

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

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

3楼: Originally posted by qingeraaa at 2014-11-07 21:41:33
看一下：http://trac-foundry.lbl.gov/trac/browser/ESPRESSO/branches/shirley_QE5.0/PW/src/potinit.f90?rev=527
你的问题出现在第184行，是由于电荷数设置问题，如果设置tot_charge应该有帮助。

tot_charge设置的是整个体系净电荷。
我是在2x2x2MgO超胞中用一个Sc原子替换一个Mg原子，净电荷还是0，不过电子数变成了奇数。

赞一下

回复此楼

4楼2014-11-08 10:11:38

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