[QE-users] ESM does not converge

Antonio Crepaldi antonio.crepaldi.lanza at gmail.com
Wed May 13 20:30:46 CEST 2020


Dear All,

I'm computing a Cu slab with esm_bc='bc3' with total charge = -1. However,
whenever the charge is !=0 i cannot converge the scf. Here you can find my
input file

&CONTROL
>    calculation      = 'scf'
>    wf_collect       = .true.
>    nstep            = 1000
>    outdir           = './Cu_esm3/'
>    prefix           = 'Cu_esm3'
>    etot_conv_thr    = 0.0001
>    forc_conv_thr    = 0.001
>    pseudo_dir       = '/homel/acrepaldi/pseudo'
> /
> &SYSTEM
>    ibrav            = 4
>    ecutwfc          = 34
>    ecutrho          = 340
>    occupations      = 'smearing'
>    degauss          = 0.002
>    smearing         = 'gaussian'
>    ntyp             = 1
>    nat              = 7
>    celldm(3)        = 11.710490425008139
>    celldm(1)        = 4.841111
>    assume_isolated='esm',
>    esm_bc='bc3'
>    tot_charge = -1
> /
> &ELECTRONS
>    electron_maxstep = 200
>    conv_thr         = 1e-08
>    mixing_mode      = 'local-TF'
>    mixing_beta      = 0.2
> /
> ATOMIC_SPECIES
> Cu 63.546 Cu_pbe_v1.2.uspp.F.UPF
>


> ATOMIC_POSITIONS angstrom
> Cu 0.0000000000 0.0000000000 -6.2637166039
> Cu 1.2809033158 0.7395306152 -4.1821166321
> Cu 0.0000000013 1.4790602298 -2.0948941353
> Cu 0.0000000000 0.0000000000 0.0009999887
> Cu 1.2809033158 0.7395306152 2.0968941353
> Cu 0.0000000013 1.4790602298 4.1841166096
> Cu 0.0000000000 0.0000000000 6.2657166039
>


> K_POINTS automatic
> 11 11 1  1 1 0


I do not understand where I am going wrong.
Thank you very much for your help.
Sincerely

Antonio Crepaldi-Lanza

--
Departmente of Physics "Aldo Moro"
University of Bari
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20200513/c7b1746f/attachment.html>


More information about the users mailing list