[QE-users] Why my Hubbard parameter is so much bigger than what I find in the litterature?

[NAIMI SALMA]

Dear Iurii,

Actually in the first attempt  I had run just one scf calculation, but after running hp.x I got the following message in the output file:

      WARNING: The Fermi energy shift is zero or too big!
      This may happen in two cases:
      1. The DOS at the Fermi level is too small:
         DOS(E_Fermi) =   -0.3416E-35
         This means that most likely the system has a gap,
         and hence it should NOT be treated as a metal
         (otherwise numerical instabilities will appear).
      2. Numerical instabilities due to too low cutoff
         for hard pseudopotentials.

      Stopping...

      Solution: Try to use the 2-step scf procedure as in HP/example02


So based on that I followed HP/example02 then and ran the second scf calculation then hp.x




-------------------------------------------------------
Salma NAIMI
PhD student - Mohamed 5 University-Rabat


________________________________
From: users <users-bounces at lists.quantum-espresso.org> on behalf of Iurii TIMROV via users <users at lists.quantum-espresso.org>
Sent: Tuesday, February 28, 2023 12:44 PM
To: Quantum ESPRESSO users Forum <users at lists.quantum-espresso.org>
Subject: Re: [QE-users] Why my Hubbard parameter is so much bigger than what I find in the litterature?


Dear Salma,


There is no need to run two SCF calculations in your case (it is only needed for magnetic insulators, see README in example02 of HP).


Cu is a problematic element for HP (linear response theory) because the 3d shell is full and in this case the approach is not suitable. For more details see J. Chem. Phys. 140, 121105 (2014).


HTH


Iurii


--
Dr. Iurii TIMROV
Senior Research Scientist
Theory and Simulation of Materials (THEOS)
Swiss Federal Institute of Technology Lausanne (EPFL)
CH-1015 Lausanne, Switzerland
+41 21 69 34 881
http://people.epfl.ch/265334
________________________________
From: users <users-bounces at lists.quantum-espresso.org> on behalf of NAIMI SALMA <salma_naimi at um5.ac.ma>
Sent: Tuesday, February 28, 2023 12:29:27 PM
To: Quantum ESPRESSO users Forum
Subject: Re: [QE-users] Why my Hubbard parameter is so much bigger than what I find in the litterature?

Dear Giovanni,

You`ll find attached the input files of scf1, scf2 and hp that I used to compute the hubbard parameter.

Can you help me to know the exact error please!




-------------------------------------------------------
Salma NAIMI
PhD student - Mohamed 5 University-Rabat

________________________________
From: users <users-bounces at lists.quantum-espresso.org> on behalf of Giovanni Cantele <giovanni.cantele at spin.cnr.it>
Sent: Tuesday, February 28, 2023 12:18 PM
To: Quantum ESPRESSO users Forum <users at lists.quantum-espresso.org>
Subject: Re: [QE-users] Why my Hubbard parameter is so much bigger than what I find in the litterature?

Dear Salma,

hard to say without knowing which material and atomic species you're dealing with. Generally speaking, if people expect a 5-6 eV U and you find more than 3 times more, I would suspect that an error is somwhere in the simulations.
You should certainly check whether you used appropriate calculation parameters, e.g. BZ sampling, cutoff(s), and so on.

Giovanni
--

Giovanni Cantele, PhD
CNR-SPIN
c/o Dipartimento di Fisica
Universita' di Napoli "Federico II"
Complesso Universitario M. S. Angelo - Ed. 6
Via Cintia, I-80126, Napoli, Italy
e-mail: giovanni.cantele at spin.cnr.it<mailto:giovanni.cantele at spin.cnr.it>
Phone: +39 081 676910
Skype contact: giocan74

ResearcherID: http://www.researcherid.com/rid/A-1951-2009
Web page: https://sites.google.com/view/giovanni-cantele/home


Il giorno mar 28 feb 2023 alle ore 12:04 NAIMI SALMA <salma_naimi at um5.ac.ma<mailto:salma_naimi at um5.ac.ma>> ha scritto:
Dear experts,
I'm working on inorganic p-type semiconductor using qe-7.0 . I want to use GGA+U, so I computed the U value by following the steps from the HP directory/example 02 for NiO (scf1 then scf2 then then hp) without including the magnetization because my material is not magnetic. My hubbard parameter was around 20 eV. While the used parameter for the same material in previous papers is around 5-6 eV ( Knowing that in previous papers they were extracting the U parameter from the litterature). The difference is big between their hubbard parameter and mine is that normal or did I do something wrong?





-------------------------------------------------------
Salma NAIMI
PhD student - Mohamed 5 University-Rabat

_______________________________________________
The Quantum ESPRESSO community stands by the Ukrainian
people and expresses its concerns about the devastating
effects that the Russian military offensive has on their
country and on the free and peaceful scientific, cultural,
and economic cooperation amongst peoples
_______________________________________________
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20230228/5a23d6d4/attachment.html>