[QE-users] How to improve the values calculated using QE?
Nicola Marzari
nicola.marzari at epfl.ch
Thu Nov 11 19:26:12 CET 2021
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
More information about the users
mailing list