<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-7">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Dear Stephen,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I recommend going through the tutorial that includes YouTube videos for DFT+U:</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<a href="https://sites.google.com/view/hubbard-koopmans/program">https://sites.google.com/view/hubbard-koopmans/program</a></div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Regarding your input files, try to use the QE input generator:</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<a href="https://qeinputgenerator.materialscloud.io/">https://qeinputgenerator.materialscloud.io/</a></div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
You can upload your existing input file, and this tool will suggest how to improve it. </div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
DFT+U and band gaps: If you apply +U to states that form the edges (top of the valence or bottom of the conduction), then the band gap will be changing. Otherwise, if you apply +U to states that are far from the edges, the band gap will not change or change
little (since those deep states that you move with +U will change the hybridization with those states that form the edges). I recommend reading this paper on this topic:</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<a href="https://www.mdpi.com/2076-3417/11/5/2395">https://www.mdpi.com/2076-3417/11/5/2395</a></div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<ol start="1" data-spread="false" style="text-align: left; background-color: rgb(255, 255, 255);">
<li style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<div role="presentation">Is it reasonable to apply Hubbard U only to Cu 3d states in Cu2O using this old syntax?</div>
</li></ol>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
Check the projected density of states and see where the Cu-3d states are. Are they close to the band edges?</div>
<ol start="2" data-spread="false" style="text-align: left; background-color: rgb(255, 255, 255);">
<li style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<div role="presentation">Is the use of Hubbard_alpha(1) = -0.08 appropriate in a production DFT+U calculation, or should it only be used for linear-response/perturbative calculations?</div>
</li></ol>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
Hubbard_alpha is the strength of the perturbation that is used to compute U using linear-response theory based on supercells and finite differences, described here:</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.71.035105">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.71.035105</a></div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<br>
</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
You need to vary Hubbard_alpha around zero and then compute the response occupations. You need to find the range of Hubbard_alpha values where the response is linear. More is explained in the paper above. </div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<br>
</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
There is a more efficient way to compute U. Namely, using the HP code, which is based on density functional perturbation theory. It gives the same result as above, but it is computationally faster. Check here if interested:</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.98.085127">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.98.085127</a></div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<a href="https://www.sciencedirect.com/science/article/pii/S0010465522001746">https://www.sciencedirect.com/science/article/pii/S0010465522001746</a></div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<br>
</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
However, unfortunately, for the system that you want to study (Cu2O) there is a problem with this method. The U will explode if you compute it self-consistently. Read more about this here:</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<a href="https://pubs.aip.org/aip/jcp/article/140/12/121105/211929/Communication-Comparing-ab-initio-methods-of">https://pubs.aip.org/aip/jcp/article/140/12/121105/211929/Communication-Comparing-ab-initio-methods-of</a></div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<br>
</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
You can also consider other methods to compute U for such systems (i.e. with closed shells):</div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<a href="https://arxiv.org/abs/2512.16803">https://arxiv.org/abs/2512.16803</a></div>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<br>
</div>
<ol start="3" data-spread="false" style="text-align: left; background-color: rgb(255, 255, 255);">
<li style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<div role="presentation">Are there any concerns with using occupations = 'fixed' for Cu2O in vc-relax/scf and occupations = 'tetrahedra' for nscf/DOS?</div>
</li></ol>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
It is ok since Cu2O is diamagnetic.</div>
<ol start="4" data-spread="false" style="text-align: left; background-color: rgb(255, 255, 255);">
<li style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
<div role="presentation">Is nbnd = 40 reasonable for a 6-atom Cu2O cell with 56 electrons for band-structure and DOS analysis?</div>
</li></ol>
<div class="elementToProof" style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
You should check how many occupied states there are in this system (it also depends on the pseudopotentials, i.e. whether the semicore states are included or not). Then add some number of empty states. The higher in energy you go for the empty states, the higher
nbnd must be. There is no universal number. It is system-dependent and it depends in which energy interval for PDOS you are interested.</div>
<ol start="5" data-spread="false" style="text-align: left; background-color: rgb(255, 255, 255);">
<li style="font-family: "Segoe UI", "Segoe UI Web (West European)", -apple-system, BlinkMacSystemFont, Roboto, "Helvetica Neue", sans-serif; font-size: 15px; color: rgb(36, 36, 36);">
Are there any syntax issues or methodological concerns that I should correct before finalizing the calculations?</li></ol>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Regarding the syntax, use QE input generator. Regarding the methodological concerns, you should not use DFT+U as a black box. There are many intricacies that one should be aware of before using this method. It looks simple, but at the same time it is very tricky. </div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
HTH</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Greetings,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Iurii</div>
<div id="Signature" class="elementToProof">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(102, 102, 102);">
----------------------------------------------------------</div>
<div class="elementToProof" style="font-family: Cambria, Georgia, serif; font-size: 12pt; color: rgb(102, 102, 102);">
Dr. Iurii TIMROV</div>
<div class="elementToProof" style="font-family: Cambria, Georgia, serif; font-size: 12pt; color: rgb(102, 102, 102);">
Tenure-track scientist</div>
<div class="elementToProof" style="font-family: Cambria, Georgia, serif; font-size: 12pt; color: rgb(102, 102, 102);">
Laboratory for Materials Simulations (LMS)</div>
<div class="elementToProof" style="font-family: Cambria, Georgia, serif; font-size: 12pt; color: rgb(102, 102, 102);">
Paul Scherrer Institute (PSI)</div>
<div class="elementToProof" style="font-family: Cambria, Georgia, serif; font-size: 12pt; color: rgb(102, 102, 102);">
CH-5232 Villigen, Switzerland</div>
<div class="elementToProof" style="font-size: 12pt;"><span style="font-family: Cambria, Georgia, serif; color: rgb(102, 102, 102);">Profile:
</span><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; color: rgb(0, 0, 0);"><a href="https://www.psi.ch/en/lms/people/iurii-timrov" target="_blank" rel="noopener noreferrer">www.psi.ch/en/lms/people/iurii-timrov</a></span></div>
<div class="elementToProof" style="font-size: 12pt;"><span style="font-family: Cambria, Georgia, serif; color: rgb(102, 102, 102);">Group website:
</span><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; color: rgb(0, 0, 0);"><a href="https://www.timrovresearch.com/">www.timrovresearch.com/</a></span></div>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Stephen Mwenda via users <users@lists.quantum-espresso.org><br>
<b>Sent:</b> Friday, May 15, 2026 17:46<br>
<b>To:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br>
<b>Subject:</b> [QE-users] Advice on Cu2O DFT+U input setup in Quantum ESPRESSO</font>
<div> </div>
</div>
<div>
<p class="x_isSelectedEnd"><span>Dear Quantum ESPRESSO users,</span></p>
<p class="x_isSelectedEnd"><span>I hope you are well.</span></p>
<p class="x_isSelectedEnd"><span>I am carrying out DFT calculations on cuprous oxide, Cu2O, using Quantum ESPRESSO. My aim is to study how the electronic band gap changes with Hubbard U applied to the Cu 3d states.</span></p>
<p class="x_isSelectedEnd"><span>I would like to kindly ask for advice on whether my input setup is reasonable, especially regarding the old-style DFT+U syntax, pseudopotentials, and use of the Hubbard alpha parameter.</span></p>
<p class="x_isSelectedEnd"><span>I have attached one representative input file for the vc-relax calculation at U = 1 eV, since this file includes the Hubbard U setup. The same workflow is then repeated for U = 2 to 6 eV by changing the value of Hubbard_U(1).
For U = 0, no Hubbard correction is applied.</span></p>
<p class="x_isSelectedEnd"><span>My calculation setup is as follows:</span></p>
<ul data-spread="false">
<li><span>Material: Cu2O, cubic cuprite structure</span></li><li><span>Exchange-correlation functional: PBE</span></li><li><span>Pseudopotentials:</span>
<ul data-spread="false">
<li><span>Cu.pbe-dn-rrkjus_psl.1.0.0.UPF</span></li><li><span>O.pbe-n-rrkjus_psl.1.0.0.UPF</span></li></ul>
</li><li><span>ecutwfc = 60 Ry</span></li><li><span>ecutrho = 480 Ry</span></li><li><span>k-points for vc-relax and scf: 6 6 6 0 0 0</span></li><li><span>k-points for nscf/DOS: 12 12 12 0 0 0</span></li><li><span>Number of atoms: 6</span></li><li><span>Number of bands in nscf: 40</span></li></ul>
<p class="x_isSelectedEnd"><span>For U = 1 eV, the relevant part of my &SYSTEM block is:</span></p>
<p class="x_isSelectedEnd"><span>&SYSTEM</span><br>
<span>ibrav = 1</span><br>
<span>a = 4.25220</span><br>
<span>nat = 6</span><br>
<span>ntyp = 2</span><br>
<span>ecutwfc = 60</span><br>
<span>ecutrho = 480</span><br>
<span>occupations = 'fixed'</span><br>
<span>lda_plus_u = .true.</span><br>
<span>Hubbard_U(1) = 1.0</span><br>
<span>Hubbard_alpha(1) = -0.08</span><br>
<span>/</span></p>
<p class="x_isSelectedEnd"><span>My ATOMIC_SPECIES block is:</span></p>
<p class="x_isSelectedEnd"><span>ATOMIC_SPECIES</span><br>
<span>Cu 63.54600 Cu.pbe-dn-rrkjus_psl.1.0.0.UPF</span><br>
<span>O 15.99940 O.pbe-n-rrkjus_psl.1.0.0.UPF</span></p>
<p class="x_isSelectedEnd"><span>Therefore, Hubbard_U(1) and Hubbard_alpha(1) are intended to apply to Cu.</span></p>
<p class="x_isSelectedEnd"><span>The relaxed structure remains cubic, and the calculated direct band gap at à increases from about 0.46 eV at U = 0 to about 0.64 eV at U = 6 eV. I understand that PBE/PBE+U may still underestimate the experimental Cu2O band
gap, but I would like to confirm whether the input syntax and general setup are technically correct.</span></p>
<p class="x_isSelectedEnd"><span>My specific questions are:</span></p>
<ol data-spread="false" start="1">
<li><span>Is it reasonable to apply Hubbard U only to Cu 3d states in Cu2O using this old syntax?</span></li><li><span>Is the use of Hubbard_alpha(1) = -0.08 appropriate in a production DFT+U calculation, or should it only be used for linear-response/perturbative calculations?</span></li><li><span>Are there any concerns with using occupations = 'fixed' for Cu2O in vc-relax/scf and occupations = 'tetrahedra' for nscf/DOS?</span></li><li><span>Is nbnd = 40 reasonable for a 6-atom Cu2O cell with 56 electrons for band-structure and DOS analysis?</span></li><li><span>Are there any syntax issues or methodological concerns that I should correct before finalizing the calculations?</span></li></ol>
<p class="x_isSelectedEnd"><span>I would be grateful for any advice or corrections from experienced users.</span></p>
<p><span>Kind regards,</span><br>
<span>Stephen</span></p>
<br>
</div>
</body>
</html>