[QE-users] How to improve the values calculated using QE?
Nicola Marzari
nicola.marzari at epfl.ch
Thu Nov 11 19:21:40 CET 2021
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>> 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)
> users mailing list users at lists.quantum-espresso.org
> 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 at http://theossrv1.epfl.ch/Main/Contact
More information about the users
mailing list