<div dir="ltr">Dear Lorenzo,<div>Thank you. I am sure the system is not metallic, the band gap (calculated with PBE) should be around 0.59 eV (vs. experimental 1.41 eV). I tried your suggestions, but with no luck (energy oscillates instead of converging). I tried changing mixing_beta from 0.05 to 0.9, increasing and decreasing number of k-points - still the same.</div><div>In the meantime, I did a calculation on exactly the same system in CASTEP without any problems (default settings, 3x3x3 Monkhorst-Pack grid), with on-the-fly USPs and also with norm-conserving ones. So clearly I am doing something wrong or there is something not right with my QE build (6.3, compiled with Intel 2017 compilers and Intel MPI). </div><div>I would be very grateful for any further suggestions.</div><div><br></div><div>Best wishes,</div><div><br></div><div>Michal</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, 9 Oct 2018 at 10:57, Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On 09/10/2018 11:52, Michal Krompiec wrote:<br>
> Hello,<br>
> I'm struggling with SCF convergence of something that should be a simple <br>
> calculation: InP (zincblende structure taken from a CIF). I tried <br>
> ultra-soft pseudopotentials, PAW, changing the k-point grid, increasing <br>
> ecutwfc and ecutrho, decreasing mixing_beta and every time the SCF <br>
> energy is oscillating instead of converging. I would be grateful for any <br>
> advice.<br>
<br>
Hello, are you sure you system does not turn out to be metallic? That <br>
would make convergence problematic.<br>
<br>
Try this the following input (I have removed all input variables that <br>
you MUST NOT use, unless you really know what you are doing). If it <br>
works remember to check convergence with ecutwfc/wcutrho (your values <br>
look very high to me) and of the k-points w.r.t. the degauss value.<br>
<br>
cheers<br>
<br>
&CONTROL<br>
    calculation = "scf",<br>
    pseudo_dir  = ".",<br>
    outdir      = ".",<br>
    prefix      = "InP",<br>
    verbosity='high'<br>
    restart_mode = "from_scratch"<br>
/<br>
&SYSTEM<br>
    ibrav       = 2,<br>
    nat         = 2,<br>
    ntyp        = 2,<br>
    ecutwfc     = 80,<br>
    ecutrho = 800,<br>
    occupations = "smearing",<br>
      smearing="mv"<br>
      degauss=0.02<br>
    A=5.8687,<br>
<br>
<br>
/<br>
&ELECTRONS<br>
   mixing_beta=0.2,<br>
   mixing_mode="plain"<br>
/<br>
<br>
<br>
ATOMIC_SPECIES<br>
In 49 In.pbe-dn-kjpaw_psl.1.0.0.UPF<br>
P 15 P.pbe-n-kjpaw_psl.1.0.0.UPF<br>
<br>
ATOMIC_POSITIONS {alat}<br>
In   0.00000   0.00000   0.00000<br>
  P   0.25 0.25 0.25<br>
<br>
K_POINTS {automatic}<br>
  9 9 9 0 0 0<br>
> <br>
> Thanks in advance,<br>
> <br>
> Michal<br>
> <br>
> <br>
> *Dr. Michal Krompiec*<br>
> <br>
> Adjunct Professor<br>
> <br>
> School of Chemistry, University of Southampton<br>
> <br>
> Highfield, Southampton SO17 1BJ, UK<br>
> <br>
> <br>
> and<br>
> <br>
> Head of Computational Modelling | Performance Materials | Early Research <br>
> and Business Development<br>
> <br>
> Merck<br>
> <br>
> <br>
> <br>
> _______________________________________________<br>
> users mailing list<br>
> <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
<br>
-- <br>
Lorenzo Paulatto - Paris<br>
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>