<div dir="auto">Thank you for responding, Nicola.<div dir="auto"><br></div><div dir="auto">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.</div><div dir="auto"><br></div><div dir="auto">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.</div><div dir="auto"><br></div><div dir="auto">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?</div><div dir="auto"> </div><div dir="auto"><br></div><div dir="auto">Tarik</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu 11 Nov, 2021, 11:56 PM Nicola Marzari, <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Blush :-)<br>
<br>
nicola<br>
<br>
<br>
On 11/11/2021 19:25, Stefano Baroni wrote:<br>
> Actually, at a closer look, and under any lighting, you are much more <br>
> handsome than Brad Pitt. SB<br>
> <br>
>> On 11 Nov 2021, at 19:21, Nicola Marzari via users <br>
>> <<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a> <br>
>> <mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a>>> wrote:<br>
>><br>
>><br>
>><br>
>> Dear Tarik,<br>
>><br>
>><br>
>> DFT is not a theory that describes electronic excitations - there is <br>
>> only the total charge density in DFT, and it describes the ground state.<br>
>><br>
>> If you use Kohn-Sham DFT to approximate the kinetic energy functional <br>
>> introducing the non-interacting Kohn-Sham particles, these lead to <br>
>> band dispersions that look a bit like the quasiparticle excitations <br>
>> you are looking for.<br>
>><br>
>> Here, "look like" is in the same sense that I look like Brad Pitt, <br>
>> under proper lightning. A fleeting illusion, that can be occasionally <br>
>> oversold as a viable approximation. But 1.18 eV is what it is, at the end.<br>
>><br>
>> There is not much that we can do, a part from switching to theories <br>
>> and methods that describe electronic excitations.<br>
>><br>
>> Hope this helps,<br>
>><br>
>> nicola<br>
>><br>
>><br>
>><br>
>> On 11/11/2021 19:03, Tarik wrote:<br>
>>> Dear Experts,<br>
>>> Kindly go through my preceding mail and suggest ways to make the <br>
>>> obtained values better. It is really important to me.<br>
>>> Thanking you,<br>
>>> Tarik<br>
>>> On Thu 11 Nov, 2021, 9:43 AM Tarik, <<a href="mailto:bitringfs@gmail.com" target="_blank" rel="noreferrer">bitringfs@gmail.com</a> <br>
>>> <mailto:<a href="mailto:bitringfs@gmail.com" target="_blank" rel="noreferrer">bitringfs@gmail.com</a>><mailto:<a href="mailto:bitringfs@gmail.com" target="_blank" rel="noreferrer">bitringfs@gmail.com</a> <br>
>>> <mailto:<a href="mailto:bitringfs@gmail.com" target="_blank" rel="noreferrer">bitringfs@gmail.com</a>>>> wrote:<br>
>>> Dear Experts,<br>
>>> I am trying to find the band gap of a material (a supercell of it)<br>
>>> and I am getting the value 1.12 eV whereas the experimental value of<br>
>>> the same is 1.8 eV. Yes, DFT does underestimate the values but it is<br>
>>> too less! It should at least be close to 1.8 eV. So can you please<br>
>>> suggest ways to improve this value?<br>
>>> I am using PBE-sol pseudopotentials.<br>
>>> This is how my scf input file looks like:<br>
>>> &CONTROL<br>
>>> calculation = 'scf'<br>
>>> etot_conv_thr = 1.3500000000d-04<br>
>>> forc_conv_thr = 1.0000000000d-05<br>
>>> outdir = './outdir'<br>
>>> prefix = 'basic'<br>
>>> pseudo_dir = './'<br>
>>> tprnfor = .true.<br>
>>> tstress = .true.<br>
>>> verbosity = 'high'<br>
>>> /<br>
>>> &SYSTEM<br>
>>> degauss = 1.4699723600d-02<br>
>>> ecutrho = 3.2000000000d+02<br>
>>> ecutwfc = 4.0000000000d+01<br>
>>> ibrav = 0<br>
>>> nat = 135<br>
>>> nosym = .false.<br>
>>> ntyp = 3<br>
>>> occupations = 'smearing'<br>
>>> smearing = 'cold'<br>
>>> /<br>
>>> &ELECTRONS<br>
>>> conv_thr = 2.7000000000d-08<br>
>>> electron_maxstep = 80<br>
>>> mixing_beta = 4.0000000000d-01<br>
>>> /<br>
>>> ATOMIC_SPECIES<br>
>>> Cs 132.9054519 cs_pbesol_v1.uspp.F.UPF<br>
>>> I 126.90447 I.pbesol-n-kjpaw_psl.0.2.UPF<br>
>>> Pb 207.2 Pb.pbesol-dn-kjpaw_psl.0.2.2.UPF<br>
>>> ATOMIC_POSITIONS (crystal)<br>
>>> Cs 0.1666666699 0.1666666699 0.1666666699<br>
>>> Pb 0.0000000000 0.0000000000 0.0000000000<br>
>>> I 0.0000000000 0.0000000000 0.1666666699<br>
>>> I 0.0000000000 0.1666666699 0.0000000000<br>
>>> I 0.1666666699 0.0000000000 0.0000000000<br>
>>> Cs 0.5000000000 0.1666666699 0.1666666699<br>
>>> Pb 0.3333333301 0.0000000000 0.0000000000<br>
>>> I 0.3333333301 0.0000000000 0.1666666699<br>
>>> I 0.3333333301 0.1666666699 0.0000000000<br>
>>> I 0.5000000000 0.0000000000 0.0000000000<br>
>>> Cs 0.8333333301 0.1666666699 0.1666666699<br>
>>> Pb 0.6666666699 0.0000000000 0.0000000000<br>
>>> I 0.6666666699 0.0000000000 0.1666666699<br>
>>> I 0.6666666699 0.1666666699 0.0000000000<br>
>>> I 0.8333333301 0.0000000000 0.0000000000<br>
>>> Cs 0.1666666699 0.5000000000 0.1666666699<br>
>>> Pb 0.0000000000 0.3333333301 0.0000000000<br>
>>> I 0.0000000000 0.3333333301 0.1666666699<br>
>>> I 0.0000000000 0.5000000000 0.0000000000<br>
>>> I 0.1666666699 0.3333333301 0.0000000000<br>
>>> Cs 0.5000000000 0.5000000000 0.1666666699<br>
>>> Pb 0.3333333301 0.3333333301 0.0000000000<br>
>>> I 0.3333333301 0.3333333301 0.1666666699<br>
>>> I 0.3333333301 0.5000000000 0.0000000000<br>
>>> I 0.5000000000 0.3333333301 0.0000000000<br>
>>> Cs 0.8333333301 0.5000000000 0.1666666699<br>
>>> Pb 0.6666666699 0.3333333301 0.0000000000<br>
>>> I 0.6666666699 0.3333333301 0.1666666699<br>
>>> I 0.6666666699 0.5000000000 0.0000000000<br>
>>> I 0.8333333301 0.3333333301 0.0000000000<br>
>>> Cs 0.1666666699 0.8333333301 0.1666666699<br>
>>> Pb 0.0000000000 0.6666666699 0.0000000000<br>
>>> I 0.0000000000 0.6666666699 0.1666666699<br>
>>> I 0.0000000000 0.8333333301 0.0000000000<br>
>>> I 0.1666666699 0.6666666699 0.0000000000<br>
>>> Cs 0.5000000000 0.8333333301 0.1666666699<br>
>>> Pb 0.3333333301 0.6666666699 0.0000000000<br>
>>> I 0.3333333301 0.6666666699 0.1666666699<br>
>>> I 0.3333333301 0.8333333301 0.0000000000<br>
>>> I 0.5000000000 0.6666666699 0.0000000000<br>
>>> Cs 0.8333333301 0.8333333301 0.1666666699<br>
>>> Pb 0.6666666699 0.6666666699 0.0000000000<br>
>>> I 0.6666666699 0.6666666699 0.1666666699<br>
>>> I 0.6666666699 0.8333333301 0.0000000000<br>
>>> I 0.8333333301 0.6666666699 0.0000000000<br>
>>> Cs 0.1666666699 0.1666666699 0.5000000000<br>
>>> Pb 0.0000000000 0.0000000000 0.3333333301<br>
>>> I 0.0000000000 0.0000000000 0.5000000000<br>
>>> I 0.0000000000 0.1666666699 0.3333333301<br>
>>> I 0.1666666699 0.0000000000 0.3333333301<br>
>>> Cs 0.5000000000 0.1666666699 0.5000000000<br>
>>> Pb 0.3333333301 0.0000000000 0.3333333301<br>
>>> I 0.3333333301 0.0000000000 0.5000000000<br>
>>> I 0.3333333301 0.1666666699 0.3333333301<br>
>>> I 0.5000000000 0.0000000000 0.3333333301<br>
>>> Cs 0.8333333301 0.1666666699 0.5000000000<br>
>>> Pb 0.6666666699 0.0000000000 0.3333333301<br>
>>> I 0.6666666699 0.0000000000 0.5000000000<br>
>>> I 0.6666666699 0.1666666699 0.3333333301<br>
>>> I 0.8333333301 0.0000000000 0.3333333301<br>
>>> Cs 0.1666666699 0.5000000000 0.5000000000<br>
>>> Pb 0.0000000000 0.3333333301 0.3333333301<br>
>>> I 0.0000000000 0.3333333301 0.5000000000<br>
>>> I 0.0000000000 0.5000000000 0.3333333301<br>
>>> I 0.1666666699 0.3333333301 0.3333333301<br>
>>> Cs 0.5000000000 0.5000000000 0.5000000000<br>
>>> Pb 0.3333333301 0.3333333301 0.3333333301<br>
>>> I 0.3333333301 0.3333333301 0.5000000000<br>
>>> I 0.3333333301 0.5000000000 0.3333333301<br>
>>> I 0.5000000000 0.3333333301 0.3333333301<br>
>>> Cs 0.8333333301 0.5000000000 0.5000000000<br>
>>> Pb 0.6666666699 0.3333333301 0.3333333301<br>
>>> I 0.6666666699 0.3333333301 0.5000000000<br>
>>> I 0.6666666699 0.5000000000 0.3333333301<br>
>>> I 0.8333333301 0.3333333301 0.3333333301<br>
>>> Cs 0.1666666699 0.8333333301 0.5000000000<br>
>>> Pb 0.0000000000 0.6666666699 0.3333333301<br>
>>> I 0.0000000000 0.6666666699 0.5000000000<br>
>>> I 0.0000000000 0.8333333301 0.3333333301<br>
>>> I 0.1666666699 0.6666666699 0.3333333301<br>
>>> Cs 0.5000000000 0.8333333301 0.5000000000<br>
>>> Pb 0.3333333301 0.6666666699 0.3333333301<br>
>>> I 0.3333333301 0.6666666699 0.5000000000<br>
>>> I 0.3333333301 0.8333333301 0.3333333301<br>
>>> I 0.5000000000 0.6666666699 0.3333333301<br>
>>> Cs 0.8333333301 0.8333333301 0.5000000000<br>
>>> Pb 0.6666666699 0.6666666699 0.3333333301<br>
>>> I 0.6666666699 0.6666666699 0.5000000000<br>
>>> I 0.6666666699 0.8333333301 0.3333333301<br>
>>> I 0.8333333301 0.6666666699 0.3333333301<br>
>>> Cs 0.1666666699 0.1666666699 0.8333333301<br>
>>> Pb 0.0000000000 0.0000000000 0.6666666699<br>
>>> I 0.0000000000 0.0000000000 0.8333333301<br>
>>> I 0.0000000000 0.1666666699 0.6666666699<br>
>>> I 0.1666666699 0.0000000000 0.6666666699<br>
>>> Cs 0.5000000000 0.1666666699 0.8333333301<br>
>>> Pb 0.3333333301 0.0000000000 0.6666666699<br>
>>> I 0.3333333301 0.0000000000 0.8333333301<br>
>>> I 0.3333333301 0.1666666699 0.6666666699<br>
>>> I 0.5000000000 0.0000000000 0.6666666699<br>
>>> Cs 0.8333333301 0.1666666699 0.8333333301<br>
>>> Pb 0.6666666699 0.0000000000 0.6666666699<br>
>>> I 0.6666666699 0.0000000000 0.8333333301<br>
>>> I 0.6666666699 0.1666666699 0.6666666699<br>
>>> I 0.8333333301 0.0000000000 0.6666666699<br>
>>> Cs 0.1666666699 0.5000000000 0.8333333301<br>
>>> Pb 0.0000000000 0.3333333301 0.6666666699<br>
>>> I 0.0000000000 0.3333333301 0.8333333301<br>
>>> I 0.0000000000 0.5000000000 0.6666666699<br>
>>> I 0.1666666699 0.3333333301 0.6666666699<br>
>>> Cs 0.5000000000 0.5000000000 0.8333333301<br>
>>> Pb 0.3333333301 0.3333333301 0.6666666699<br>
>>> I 0.3333333301 0.3333333301 0.8333333301<br>
>>> I 0.3333333301 0.5000000000 0.6666666699<br>
>>> I 0.5000000000 0.3333333301 0.6666666699<br>
>>> Cs 0.8333333301 0.5000000000 0.8333333301<br>
>>> Pb 0.6666666699 0.3333333301 0.6666666699<br>
>>> I 0.6666666699 0.3333333301 0.8333333301<br>
>>> I 0.6666666699 0.5000000000 0.6666666699<br>
>>> I 0.8333333301 0.3333333301 0.6666666699<br>
>>> Cs 0.1666666699 0.8333333301 0.8333333301<br>
>>> Pb 0.0000000000 0.6666666699 0.6666666699<br>
>>> I 0.0000000000 0.6666666699 0.8333333301<br>
>>> I 0.0000000000 0.8333333301 0.6666666699<br>
>>> I 0.1666666699 0.6666666699 0.6666666699<br>
>>> Cs 0.5000000000 0.8333333301 0.8333333301<br>
>>> Pb 0.3333333301 0.6666666699 0.6666666699<br>
>>> I 0.3333333301 0.6666666699 0.8333333301<br>
>>> I 0.3333333301 0.8333333301 0.6666666699<br>
>>> I 0.5000000000 0.6666666699 0.6666666699<br>
>>> Cs 0.8333333301 0.8333333301 0.8333333301<br>
>>> Pb 0.6666666699 0.6666666699 0.6666666699<br>
>>> I 0.6666666699 0.6666666699 0.8333333301<br>
>>> I 0.6666666699 0.8333333301 0.6666666699<br>
>>> I 0.8333333301 0.6666666699 0.6666666699<br>
>>> CELL_PARAMETERS (angstrom)<br>
>>> 18.712038000 0.000000000 0.000000000<br>
>>> 0.000000000 18.712038000 0.000000000<br>
>>> 0.000000000 0.000000000 18.712038000<br>
>>> K_POINTS automatic<br>
>>> 5 5 5 0 0 0<br>
>>> Thanking you,<br>
>>> Tarik<br>
>>> IIT Indore<br>
>>> _______________________________________________<br>
>>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a> <br>
>>> <<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">http://www.max-centre.eu</a>>)<br>
>>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a> <br>
>>> <mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a>><br>
>>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a> <br>
>>> <<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a>><br>
>><br>
>><br>
>> --<br>
>> ----------------------------------------------------------------------<br>
>> Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
>> Director, National Centre for Competence in Research NCCR MARVEL, SNSF<br>
>> Head, Laboratory for Materials Simulations, Paul Scherrer Institut<br>
>> Contact info and websites athttp://<a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer noreferrer" target="_blank">theossrv1.epfl.ch/Main/Contact</a> <br>
>> <<a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a>><br>
>> _______________________________________________<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a> <br>
>> <<a href="http://www.max-centre.eu/" rel="noreferrer noreferrer" target="_blank">http://www.max-centre.eu/</a>>)<br>
>> users mailing <a href="mailto:listusers@lists.quantum-espresso.org" target="_blank" rel="noreferrer">listusers@lists.quantum-espresso.org</a> <br>
>> <mailto:<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a>><br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a> <br>
>> <<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a>><br>
> <br>
> —<br>
> Stefano Baroni - Trieste — <a href="http://stefano.baroni.me" rel="noreferrer noreferrer" target="_blank">http://stefano.baroni.me</a> <br>
> <<a href="http://stefano.baroni.me" rel="noreferrer noreferrer" target="_blank">http://stefano.baroni.me</a>><br>
> <br>
> <br>
> <br>
> <br>
<br>
<br>
-- <br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, SNSF<br>
Head, Laboratory for Materials Simulations, Paul Scherrer Institut<br>
Contact info and websites at <a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a><br>
</blockquote></div>