[QE-users] TB09, TPSS
Paolo Giannozzi
p.giannozzi at gmail.com
Thu May 21 08:46:22 CEST 2020
On Wed, May 20, 2020 at 9:34 PM Michal Krompiec <michal.krompiec at gmail.com>
wrote:
> I changed io_rho_xml.f90 so that if the kinetic energy density file isn't
> found, the array is set to 0:
> [...]
> Now I can start a TB09 calculation from a PBE density, but still my test
> case diverges immediately (...)
> I would be very grateful for any suggestions.
>
unfortunately I haven't any.
Paolo
Best,
> Michal
>
> P.S. This is my input file:
> &CONTROL
> calculation = "scf"
> prefix = "sic"
> /
>
> &SYSTEM
> a = 4.34800e+00
> ibrav = 1
> ecutwfc = 82
> ecutrho = 328
> nat = 8
> ntyp = 2
> occupations = "fixed"
> input_dft='tb09'
> /
> &ELECTRONS
> electron_maxstep = 1000
> conv_thr = 1e-10
> mixing_mode = 'plain'
> mixing_beta = 0.4
> startingpot = 'file'
> startingwfc = 'file'
> diagonalization = 'cg'
> /
>
> K_POINTS automatic
> 8 8 8 1 1 1
>
> ATOMIC_SPECIES
> Si 28.08500 Si.upf
> C 12.01060 C.upf
>
>
> ATOMIC_POSITIONS {angstrom}
> Si 0.000000 0.000000 0.000000
> Si 0.000000 2.174000 2.174000
> Si 2.174000 0.000000 2.174000
> Si 2.174000 2.174000 0.000000
> C 1.087000 1.087000 1.087000
> C 1.087000 3.261000 3.261000
> C 3.261000 1.087000 3.261000
> C 3.261000 3.261000 1.087000
>
> On Fri, 8 May 2020 at 09:39, Paolo Giannozzi <p.giannozzi at gmail.com>
> wrote:
> >
> > On Thu, May 7, 2020 at 10:33 PM Michal Krompiec <
> michal.krompiec at gmail.com> wrote:
> >
> >>
> >> 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).
> >
> >
> > this can be easily changed, I think: just disable the check and set the
> kinetic energy density to zero
> >
> >>
> >> Should TPSS or SCAN pseudos work better (and where do I get them from)?
> >
> >
> > I don't think the problem is in the pseudos but in the nasty numerical
> behavior of meta-GGAs. See for instance here:
> https://gitlab.com/QEF/q-e/-/issues/32. Note that several bugs have been
> fixed in the development version.
> >
> > Paolo
> >
> >> Thanks,
> >> Michal Krompiec
> >> Merck KGaA
> >>
> >> _______________________________________________
> >> Quantum ESPRESSO is supported by MaX (
> www.max-centre.eu/quantum-espresso)
> >> users mailing list users at lists.quantum-espresso.org
> >> https://lists.quantum-espresso.org/mailman/listinfo/users
> >
> >
> >
> > --
> > Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,
> > Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
> > Phone +39-0432-558216, fax +39-0432-558222
> >
> > _______________________________________________
> > Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso
> )
> > users mailing list users at lists.quantum-espresso.org
> > https://lists.quantum-espresso.org/mailman/listinfo/users
> _______________________________________________
> Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso)
> users mailing list users at lists.quantum-espresso.org
> https://lists.quantum-espresso.org/mailman/listinfo/users
>
--
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
Phone +39-0432-558216, fax +39-0432-558222
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20200521/f41b416f/attachment.html>
More information about the users
mailing list