24小时热门版块排行榜    

查看: 1671  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

骑英奔腾

铜虫 (小有名气)

[求助] 错误ecutwfc out of range 该如何解决

输入文件
&control
    prefix='MoS2_lattice',
    calculation='scf',
    pseudo_dir='/home_soft/home/caodp/pseudo/',
    outdir='/home_soft/home/caodp/tmp/'
/
&system   
    ibrav=4,
    celldm(1) =9.55,
    celldm(3)=23.1702/9.55,
    nat=27,
    ntyp=2,
    ecutwfc=40,
    ecutrho = 400,
    occupations='smearing',
    smearing='marzari-vanderbilt',
    degauss=0.002,
/
&electrons
/
ATOMIC_SPECIES
Mo  95.94  Mo.pbe-mt_fhi.UPF
S   32.066  S.pbe-mt_fhi.UPF
ATOMIC_POSITIONS
S    0.97003    1.32501    0.01866
S    0.97003    1.32501    3.14573
S    4.15252    1.32501    0.01866
S    4.15252    1.32501    3.14573
S    7.33503    1.32501    0.01866
S    7.33503    1.32501    3.14573
S    2.56128    4.08114    0.01866
S    2.56128    4.08114    3.14573
S    5.74378    4.08114    0.01866
S    5.74378    4.08114    3.14573
S    8.92628    4.08114    0.01866
S    8.92628    4.08114    3.14573
S    4.15252    6.83726    0.01866
S    4.15252    6.83726    3.14573
S    7.33502    6.83726    0.01866
S    7.33502    6.83726    3.14573
S   10.51752    6.83726    0.01866
S   10.51752    6.83726    3.14573
Mo    2.56128    0.40630    1.58232
Mo    5.74377    0.40630    1.58232
Mo    8.92628    0.40630    1.58232
Mo    4.15252    3.16243    1.58232
Mo    7.33502    3.16243    1.58232
Mo   10.51752    3.16243    1.58232
Mo    5.74377    5.91855    1.58232
Mo    8.92627    5.91855    1.58232
Mo   12.10878    5.91855    1.58232
K_POINTS automatic
  2 2 2 1 1 1

crash内容:
task #        45
     from  system_checkin      : error #         1
      ecutwfc out of range
请问大家该怎么改
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

引用回帖:
4楼: Originally posted by 骑英奔腾 at 2013-10-11 22:18:32
你说的这些和ecutwfc out of range 有关系吗? 我改了之后还是出现同样问题。...

贴个我算MoS2的输入文件,用的赝势和你的一样,只是把名字改了。

&CONTROL
  prefix           = 'mos2'
  calculation      = 'scf'
  restart_mode     = 'from_scratch'
  pseudo_dir       = './'
  outdir           = './'
  wfcdir           = './'
  verbosity        = 'high'
  wf_collect       = .TRUE.
  nstep            =  250
/
&SYSTEM
  ibrav            = 4
  celldm(1)        = 6.025347816
  celldm(3)        = 5.645328914
  ntyp             = 2
  nat              = 3
  ecutwfc          = 45.0
  force_symmorphic = .TRUE.
/
&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
   S  32.066   S.pbe-mt.UPF
  Mo  95.940  Mo.pbe-mt.UPF
ATOMIC_POSITIONS {crystal}
   S   0.000000000   0.000000000   0.412540672
   S   0.000000000   0.000000000   0.587459328
  Mo   0.666666667   0.333333333   0.500000000
K_POINTS {automatic}
12 12 1 0 0 0
5楼2013-10-12 16:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

骑英奔腾

铜虫 (小有名气)

我已经试了ecutwfc从10到100都出错
2楼2013-10-11 16:54:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
从赝势来看,你用的是模守恒赝势。对于模守恒赝势,只需设置ecutwfc,ecutrho默认为ecutwfc4倍。只有用超软赝势时才需要改变这个默认设置。

另外,二硫化钼不用加展宽也能收敛。
3楼2013-10-11 17:15:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

骑英奔腾

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by liqizuiyang at 2013-10-11 17:15:14
从赝势来看,你用的是模守恒赝势。对于模守恒赝势,只需设置ecutwfc,ecutrho默认为ecutwfc4倍。只有用超软赝势时才需要改变这个默认设置。

另外,二硫化钼不用加展宽也能收敛。

你说的这些和ecutwfc out of range 有关系吗? 我改了之后还是出现同样问题。
4楼2013-10-11 22:18:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见