<div dir="ltr"><div dir="ltr">On Wed, May 20, 2020 at 9:34 PM Michal Krompiec <<a href="mailto:michal.krompiec@gmail.com">michal.krompiec@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I changed io_rho_xml.f90 so that if the kinetic energy density file isn't found, the array is set to 0:<br>[...]<br>
Now I can start a TB09 calculation from a PBE density, but still my test case diverges immediately (...)<br>
I would be very grateful for any suggestions.<br></blockquote><div><br></div><div>unfortunately I haven't any. <br></div><div><br></div><div>Paolo</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Best,<br>
Michal<br>
<br>
P.S. This is my input file:<br>
&CONTROL<br>
  Â  calculation  = "scf"<br>
  Â  prefix  Â  Â  Â = "sic"<br>
/<br>
<br>
&SYSTEM<br>
  Â  a  Â  Â  Â  Â  Â =  4.34800e+00<br>
  Â  ibrav  Â  Â  Â = 1<br>
  Â ecutwfc  Â  Â  Â  Â  = 82<br>
  Â ecutrho  Â  Â  Â  Â  = 328<br>
  Â  nat  Â  Â  Â  Â = 8<br>
  Â  ntyp  Â  Â  Â  = 2<br>
  Â  occupations = "fixed"<br>
input_dft='tb09'<br>
/<br>
&ELECTRONS<br>
  Â electron_maxstep = 1000<br>
  Â conv_thr  Â  Â  Â  Â = 1e-10<br>
  Â mixing_mode  Â  Â  = 'plain'<br>
  Â mixing_beta  Â  Â  = 0.4<br>
  startingpot = 'file'<br>
  startingwfc = 'file'<br>
  diagonalization = 'cg'<br>
/<br>
<br>
K_POINTS automatic<br>
8 8 8  Â 1 1 1<br>
<br>
ATOMIC_SPECIES<br>
Si  Â  Â 28.08500  Si.upf<br>
C  Â  Â  12.01060  C.upf<br>
<br>
<br>
ATOMIC_POSITIONS {angstrom}<br>
Si  Â  Â  0.000000  Â 0.000000  Â 0.000000<br>
Si  Â  Â  0.000000  Â 2.174000  Â 2.174000<br>
Si  Â  Â  2.174000  Â 0.000000  Â 2.174000<br>
Si  Â  Â  2.174000  Â 2.174000  Â 0.000000<br>
C  Â  Â  Â 1.087000  Â 1.087000  Â 1.087000<br>
C  Â  Â  Â 1.087000  Â 3.261000  Â 3.261000<br>
C  Â  Â  Â 3.261000  Â 1.087000  Â 3.261000<br>
C  Â  Â  Â 3.261000  Â 3.261000  Â 1.087000<br>
<br>
On Fri, 8 May 2020 at 09:39, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>> wrote:<br>
><br>
> On Thu, May 7, 2020 at 10:33 PM Michal Krompiec <<a href="mailto:michal.krompiec@gmail.com" target="_blank">michal.krompiec@gmail.com</a>> wrote:<br>
><br>
>><br>
>>  it was suggested to start from a density calculated with a different functional, but when I try to read in PBE density, it complains that it cannot read the kinetic energy file (which obviously cannot be there).<br>
><br>
><br>
> this can be easily changed, I think: just disable the check and set the kinetic energy density to zero<br>
><br>
>><br>
>> Should TPSS or SCAN pseudos work better (and where do I get them from)?<br>
><br>
><br>
> I don't think the problem is in the pseudos but in the nasty numerical behavior of meta-GGAs. See for instance here: <a href="https://gitlab.com/QEF/q-e/-/issues/32" rel="noreferrer" target="_blank">https://gitlab.com/QEF/q-e/-/issues/32</a>. Note that several bugs have been fixed in the development version.<br>
><br>
> Paolo<br>
><br>
>> Thanks,<br>
>> Michal Krompiec<br>
>> Merck KGaA<br>
>><br>
>> _______________________________________________<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
>> users mailing list <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>
> --<br>
> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
> Phone +39-0432-558216, fax +39-0432-558222<br>
><br>
> _______________________________________________<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
> users mailing list <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>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <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>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>