[QE-users] How to improve the values calculated using QE?
Tarik
bitringfs at gmail.com
Thu Nov 11 19:45:09 CET 2021
Thank you for responding, Nicola.
Yes I understand, we cannot do much in QE but I was wondering if changing
ecutrho or ecutwfc might bring better results or some other parameter can
be changed.
I heard Wein2k or VASP gives more accurate results? Infact, most of the
papers (literature) say that they use VASP for computing. Why is that so? I
do not want to believe Quantum Espresso is anything less.
Also, can you please tell me how to implement different approximations in
QE? Say for example, I want to use the GGA-PBE or the PBE-sol approx., how
do I incorporate that in my input code?
Tarik
On Thu 11 Nov, 2021, 11:56 PM Nicola Marzari, <nicola.marzari at epfl.ch>
wrote:
>
> Blush :-)
>
> nicola
>
>
> On 11/11/2021 19:25, Stefano Baroni wrote:
> > Actually, at a closer look, and under any lighting, you are much more
> > handsome than Brad Pitt. SB
> >
> >> On 11 Nov 2021, at 19:21, Nicola Marzari via users
> >> <users at lists.quantum-espresso.org
> >> <mailto:users at lists.quantum-espresso.org>> wrote:
> >>
> >>
> >>
> >> Dear Tarik,
> >>
> >>
> >> DFT is not a theory that describes electronic excitations - there is
> >> only the total charge density in DFT, and it describes the ground state.
> >>
> >> If you use Kohn-Sham DFT to approximate the kinetic energy functional
> >> introducing the non-interacting Kohn-Sham particles, these lead to
> >> band dispersions that look a bit like the quasiparticle excitations
> >> you are looking for.
> >>
> >> Here, "look like" is in the same sense that I look like Brad Pitt,
> >> under proper lightning. A fleeting illusion, that can be occasionally
> >> oversold as a viable approximation. But 1.18 eV is what it is, at the
> end.
> >>
> >> There is not much that we can do, a part from switching to theories
> >> and methods that describe electronic excitations.
> >>
> >> Hope this helps,
> >>
> >> nicola
> >>
> >>
> >>
> >> On 11/11/2021 19:03, Tarik wrote:
> >>> Dear Experts,
> >>> Kindly go through my preceding mail and suggest ways to make the
> >>> obtained values better. It is really important to me.
> >>> Thanking you,
> >>> Tarik
> >>> On Thu 11 Nov, 2021, 9:43 AM Tarik, <bitringfs at gmail.com
> >>> <mailto:bitringfs at gmail.com><mailto:bitringfs at gmail.com
> >>> <mailto:bitringfs at gmail.com>>> wrote:
> >>> Dear Experts,
> >>> I am trying to find the band gap of a material (a supercell of it)
> >>> and I am getting the value 1.12 eV whereas the experimental value of
> >>> the same is 1.8 eV. Yes, DFT does underestimate the values but it is
> >>> too less! It should at least be close to 1.8 eV. So can you please
> >>> suggest ways to improve this value?
> >>> I am using PBE-sol pseudopotentials.
> >>> This is how my scf input file looks like:
> >>> &CONTROL
> >>> calculation = 'scf'
> >>> etot_conv_thr = 1.3500000000d-04
> >>> forc_conv_thr = 1.0000000000d-05
> >>> outdir = './outdir'
> >>> prefix = 'basic'
> >>> pseudo_dir = './'
> >>> tprnfor = .true.
> >>> tstress = .true.
> >>> verbosity = 'high'
> >>> /
> >>> &SYSTEM
> >>> degauss = 1.4699723600d-02
> >>> ecutrho = 3.2000000000d+02
> >>> ecutwfc = 4.0000000000d+01
> >>> ibrav = 0
> >>> nat = 135
> >>> nosym = .false.
> >>> ntyp = 3
> >>> occupations = 'smearing'
> >>> smearing = 'cold'
> >>> /
> >>> &ELECTRONS
> >>> conv_thr = 2.7000000000d-08
> >>> electron_maxstep = 80
> >>> mixing_beta = 4.0000000000d-01
> >>> /
> >>> ATOMIC_SPECIES
> >>> Cs 132.9054519 cs_pbesol_v1.uspp.F.UPF
> >>> I 126.90447 I.pbesol-n-kjpaw_psl.0.2.UPF
> >>> Pb 207.2 Pb.pbesol-dn-kjpaw_psl.0.2.2.UPF
> >>> ATOMIC_POSITIONS (crystal)
> >>> Cs 0.1666666699 0.1666666699 0.1666666699
> >>> Pb 0.0000000000 0.0000000000 0.0000000000
> >>> I 0.0000000000 0.0000000000 0.1666666699
> >>> I 0.0000000000 0.1666666699 0.0000000000
> >>> I 0.1666666699 0.0000000000 0.0000000000
> >>> Cs 0.5000000000 0.1666666699 0.1666666699
> >>> Pb 0.3333333301 0.0000000000 0.0000000000
> >>> I 0.3333333301 0.0000000000 0.1666666699
> >>> I 0.3333333301 0.1666666699 0.0000000000
> >>> I 0.5000000000 0.0000000000 0.0000000000
> >>> Cs 0.8333333301 0.1666666699 0.1666666699
> >>> Pb 0.6666666699 0.0000000000 0.0000000000
> >>> I 0.6666666699 0.0000000000 0.1666666699
> >>> I 0.6666666699 0.1666666699 0.0000000000
> >>> I 0.8333333301 0.0000000000 0.0000000000
> >>> Cs 0.1666666699 0.5000000000 0.1666666699
> >>> Pb 0.0000000000 0.3333333301 0.0000000000
> >>> I 0.0000000000 0.3333333301 0.1666666699
> >>> I 0.0000000000 0.5000000000 0.0000000000
> >>> I 0.1666666699 0.3333333301 0.0000000000
> >>> Cs 0.5000000000 0.5000000000 0.1666666699
> >>> Pb 0.3333333301 0.3333333301 0.0000000000
> >>> I 0.3333333301 0.3333333301 0.1666666699
> >>> I 0.3333333301 0.5000000000 0.0000000000
> >>> I 0.5000000000 0.3333333301 0.0000000000
> >>> Cs 0.8333333301 0.5000000000 0.1666666699
> >>> Pb 0.6666666699 0.3333333301 0.0000000000
> >>> I 0.6666666699 0.3333333301 0.1666666699
> >>> I 0.6666666699 0.5000000000 0.0000000000
> >>> I 0.8333333301 0.3333333301 0.0000000000
> >>> Cs 0.1666666699 0.8333333301 0.1666666699
> >>> Pb 0.0000000000 0.6666666699 0.0000000000
> >>> I 0.0000000000 0.6666666699 0.1666666699
> >>> I 0.0000000000 0.8333333301 0.0000000000
> >>> I 0.1666666699 0.6666666699 0.0000000000
> >>> Cs 0.5000000000 0.8333333301 0.1666666699
> >>> Pb 0.3333333301 0.6666666699 0.0000000000
> >>> I 0.3333333301 0.6666666699 0.1666666699
> >>> I 0.3333333301 0.8333333301 0.0000000000
> >>> I 0.5000000000 0.6666666699 0.0000000000
> >>> Cs 0.8333333301 0.8333333301 0.1666666699
> >>> Pb 0.6666666699 0.6666666699 0.0000000000
> >>> I 0.6666666699 0.6666666699 0.1666666699
> >>> I 0.6666666699 0.8333333301 0.0000000000
> >>> I 0.8333333301 0.6666666699 0.0000000000
> >>> Cs 0.1666666699 0.1666666699 0.5000000000
> >>> Pb 0.0000000000 0.0000000000 0.3333333301
> >>> I 0.0000000000 0.0000000000 0.5000000000
> >>> I 0.0000000000 0.1666666699 0.3333333301
> >>> I 0.1666666699 0.0000000000 0.3333333301
> >>> Cs 0.5000000000 0.1666666699 0.5000000000
> >>> Pb 0.3333333301 0.0000000000 0.3333333301
> >>> I 0.3333333301 0.0000000000 0.5000000000
> >>> I 0.3333333301 0.1666666699 0.3333333301
> >>> I 0.5000000000 0.0000000000 0.3333333301
> >>> Cs 0.8333333301 0.1666666699 0.5000000000
> >>> Pb 0.6666666699 0.0000000000 0.3333333301
> >>> I 0.6666666699 0.0000000000 0.5000000000
> >>> I 0.6666666699 0.1666666699 0.3333333301
> >>> I 0.8333333301 0.0000000000 0.3333333301
> >>> Cs 0.1666666699 0.5000000000 0.5000000000
> >>> Pb 0.0000000000 0.3333333301 0.3333333301
> >>> I 0.0000000000 0.3333333301 0.5000000000
> >>> I 0.0000000000 0.5000000000 0.3333333301
> >>> I 0.1666666699 0.3333333301 0.3333333301
> >>> Cs 0.5000000000 0.5000000000 0.5000000000
> >>> Pb 0.3333333301 0.3333333301 0.3333333301
> >>> I 0.3333333301 0.3333333301 0.5000000000
> >>> I 0.3333333301 0.5000000000 0.3333333301
> >>> I 0.5000000000 0.3333333301 0.3333333301
> >>> Cs 0.8333333301 0.5000000000 0.5000000000
> >>> Pb 0.6666666699 0.3333333301 0.3333333301
> >>> I 0.6666666699 0.3333333301 0.5000000000
> >>> I 0.6666666699 0.5000000000 0.3333333301
> >>> I 0.8333333301 0.3333333301 0.3333333301
> >>> Cs 0.1666666699 0.8333333301 0.5000000000
> >>> Pb 0.0000000000 0.6666666699 0.3333333301
> >>> I 0.0000000000 0.6666666699 0.5000000000
> >>> I 0.0000000000 0.8333333301 0.3333333301
> >>> I 0.1666666699 0.6666666699 0.3333333301
> >>> Cs 0.5000000000 0.8333333301 0.5000000000
> >>> Pb 0.3333333301 0.6666666699 0.3333333301
> >>> I 0.3333333301 0.6666666699 0.5000000000
> >>> I 0.3333333301 0.8333333301 0.3333333301
> >>> I 0.5000000000 0.6666666699 0.3333333301
> >>> Cs 0.8333333301 0.8333333301 0.5000000000
> >>> Pb 0.6666666699 0.6666666699 0.3333333301
> >>> I 0.6666666699 0.6666666699 0.5000000000
> >>> I 0.6666666699 0.8333333301 0.3333333301
> >>> I 0.8333333301 0.6666666699 0.3333333301
> >>> Cs 0.1666666699 0.1666666699 0.8333333301
> >>> Pb 0.0000000000 0.0000000000 0.6666666699
> >>> I 0.0000000000 0.0000000000 0.8333333301
> >>> I 0.0000000000 0.1666666699 0.6666666699
> >>> I 0.1666666699 0.0000000000 0.6666666699
> >>> Cs 0.5000000000 0.1666666699 0.8333333301
> >>> Pb 0.3333333301 0.0000000000 0.6666666699
> >>> I 0.3333333301 0.0000000000 0.8333333301
> >>> I 0.3333333301 0.1666666699 0.6666666699
> >>> I 0.5000000000 0.0000000000 0.6666666699
> >>> Cs 0.8333333301 0.1666666699 0.8333333301
> >>> Pb 0.6666666699 0.0000000000 0.6666666699
> >>> I 0.6666666699 0.0000000000 0.8333333301
> >>> I 0.6666666699 0.1666666699 0.6666666699
> >>> I 0.8333333301 0.0000000000 0.6666666699
> >>> Cs 0.1666666699 0.5000000000 0.8333333301
> >>> Pb 0.0000000000 0.3333333301 0.6666666699
> >>> I 0.0000000000 0.3333333301 0.8333333301
> >>> I 0.0000000000 0.5000000000 0.6666666699
> >>> I 0.1666666699 0.3333333301 0.6666666699
> >>> Cs 0.5000000000 0.5000000000 0.8333333301
> >>> Pb 0.3333333301 0.3333333301 0.6666666699
> >>> I 0.3333333301 0.3333333301 0.8333333301
> >>> I 0.3333333301 0.5000000000 0.6666666699
> >>> I 0.5000000000 0.3333333301 0.6666666699
> >>> Cs 0.8333333301 0.5000000000 0.8333333301
> >>> Pb 0.6666666699 0.3333333301 0.6666666699
> >>> I 0.6666666699 0.3333333301 0.8333333301
> >>> I 0.6666666699 0.5000000000 0.6666666699
> >>> I 0.8333333301 0.3333333301 0.6666666699
> >>> Cs 0.1666666699 0.8333333301 0.8333333301
> >>> Pb 0.0000000000 0.6666666699 0.6666666699
> >>> I 0.0000000000 0.6666666699 0.8333333301
> >>> I 0.0000000000 0.8333333301 0.6666666699
> >>> I 0.1666666699 0.6666666699 0.6666666699
> >>> Cs 0.5000000000 0.8333333301 0.8333333301
> >>> Pb 0.3333333301 0.6666666699 0.6666666699
> >>> I 0.3333333301 0.6666666699 0.8333333301
> >>> I 0.3333333301 0.8333333301 0.6666666699
> >>> I 0.5000000000 0.6666666699 0.6666666699
> >>> Cs 0.8333333301 0.8333333301 0.8333333301
> >>> Pb 0.6666666699 0.6666666699 0.6666666699
> >>> I 0.6666666699 0.6666666699 0.8333333301
> >>> I 0.6666666699 0.8333333301 0.6666666699
> >>> I 0.8333333301 0.6666666699 0.6666666699
> >>> CELL_PARAMETERS (angstrom)
> >>> 18.712038000 0.000000000 0.000000000
> >>> 0.000000000 18.712038000 0.000000000
> >>> 0.000000000 0.000000000 18.712038000
> >>> K_POINTS automatic
> >>> 5 5 5 0 0 0
> >>> Thanking you,
> >>> Tarik
> >>> IIT Indore
> >>> _______________________________________________
> >>> Quantum ESPRESSO is supported by MaX (www.max-centre.eu
> >>> <http://www.max-centre.eu>)
> >>> users mailing list users at lists.quantum-espresso.org
> >>> <mailto:users at lists.quantum-espresso.org>
> >>> https://lists.quantum-espresso.org/mailman/listinfo/users
> >>> <https://lists.quantum-espresso.org/mailman/listinfo/users>
> >>
> >>
> >> --
> >> ----------------------------------------------------------------------
> >> Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL
> >> Director, National Centre for Competence in Research NCCR MARVEL, SNSF
> >> Head, Laboratory for Materials Simulations, Paul Scherrer Institut
> >> Contact info and websites athttp://theossrv1.epfl.ch/Main/Contact
> >> <http://theossrv1.epfl.ch/Main/Contact>
> >> _______________________________________________
> >> Quantum ESPRESSO is supported by MaX (www.max-centre.eu
> >> <http://www.max-centre.eu/>)
> >> users mailing listusers at lists.quantum-espresso.org
> >> <mailto:users at lists.quantum-espresso.org>
> >> https://lists.quantum-espresso.org/mailman/listinfo/users
> >> <https://lists.quantum-espresso.org/mailman/listinfo/users>
> >
> > —
> > Stefano Baroni - Trieste — http://stefano.baroni.me
> > <http://stefano.baroni.me>
> >
> >
> >
> >
>
>
> --
> ----------------------------------------------------------------------
> Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL
> Director, National Centre for Competence in Research NCCR MARVEL, SNSF
> Head, Laboratory for Materials Simulations, Paul Scherrer Institut
> Contact info and websites at http://theossrv1.epfl.ch/Main/Contact
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20211112/b337e54b/attachment.html>
More information about the users
mailing list