<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:large">Dear Dr. Lurii,</div><div class="gmail_default" style="font-size:large">Thank you for your important suggestions.</div><div class="gmail_default" style="font-size:large">I have started HP code and got some errors which I have posted in a separate email.</div><div class="gmail_default" style="font-size:large">I am looking into more details and will get back to you with some updates.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Regards</div><div class="gmail_default" style="font-size:large">Bhamu</div><div class="gmail_default" style="font-size:large"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Nov 15, 2021 at 9:55 PM Iurii TIMROV via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div id="gmail-m_-3196822762173715662divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p><span style="font-size:12pt;font-family:Arial,Helvetica,sans-serif;color:rgb(0,0,0)">Dear
</span><font size="4"><font size="4"><span style="font-size:12pt;font-family:Arial,Helvetica,sans-serif;color:rgb(0,0,0)">K C Bhamu</span></font></font><span style="font-size:12pt;font-family:Arial,Helvetica,sans-serif;color:rgb(0,0,0)">,</span></p>
<p><br>
</p>
<p><font size="4"><font size="4">> Do you think my input file (see below) has any issues? I have defined all tetra-Fe sites dn-spin and all octahedra Fe sites up-spin.
<br>
</font></font></p>
<p><br>
</p>
<p>I recommend to check with the QE input generator: <a href="https://www.materialscloud.org/work/tools/qeinputgenerator" id="gmail-m_-3196822762173715662LPlnk63756" target="_blank">
https://www.materialscloud.org/work/tools/qeinputgenerator</a></p>
<p><br>
</p>
<p><font size="4"><font size="4">> What other things I should take care of to reproduce the band gap (0.14 eV to 0.3 eV)?</font></font><br>
</p>
<p><br>
</p>
<p>I never studied this system but I can give general comments:</p>
<p>- you can try to determine Hubbard parameters using the HP code of QE<br>
</p>
<p>- you can try to use "ortho-atomic" Hubbard projectors instead of "atomic" (<span>U_projection_type</span>):
<a href="https://www.mdpi.com/2076-3417/11/5/2395" id="gmail-m_-3196822762173715662LPlnk854714" target="_blank">
https://www.mdpi.com/2076-3417/11/5/2395</a></p>
<p>- maybe try DFT+U+V? <a href="https://journals.aps.org/prmaterials/abstract/10.1103/PhysRevMaterials.5.104402" id="gmail-m_-3196822762173715662LPlnk263762" target="_blank">
https://journals.aps.org/prmaterials/abstract/10.1103/PhysRevMaterials.5.104402</a> (see Fig. 6) Tutorial:
<a href="https://www.youtube.com/watch?v=WSABAqPWNH0&t=5s" id="gmail-m_-3196822762173715662LPlnk318874" target="_blank">
https://www.youtube.com/watch?v=WSABAqPWNH0&t=5s</a></p>
<p><br>
</p>
<p><font size="4"><font size="4">> In the </font><a href="https://iopscience.iop.org/article/10.1088/1468-6996/15/4/044202" rel="noopener noreferrer" target="_blank"><font size="4">research paper [3]</font></a><font size="4">, the authors used VASP and the
valence configuration for the Fe atom was taken as 3d64s1. I could not find any PP from the list I mention below which has this configuration. Does this may be an issue? </font>If anyone has PP of Fe with 3d64s1 configuration, please provide me.</font><br>
</p>
<p><br>
</p>
<p>Check SSSP: <a href="https://www.materialscloud.org/discover/sssp/table/efficiency" id="gmail-m_-3196822762173715662LPlnk506580" target="_blank">
https://www.materialscloud.org/discover/sssp/table/efficiency</a></p>
<p><br>
</p>
<p>HTH</p>
<p><br>
</p>
<p>Iurii<br>
</p>
<p><br>
</p>
<div id="gmail-m_-3196822762173715662Signature">
<div id="gmail-m_-3196822762173715662divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div name="divtagdefaultwrapper">
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
<span style="font-family:Cambria,Georgia,serif">Dr. Iurii TIMROV</span><br>
<span style="font-family:Cambria,Georgia,serif">Senior Research Scientist</span></font></div>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">Theory and Simulation of Materials (THEOS)</span></font></div>
<div name="divtagdefaultwrapper">
<span style="font-size:12pt;color:rgb(128,128,128);font-family:Cambria,Georgia,serif">Swiss Federal Institute of Technology Lausanne (EPFL</span><font style="font-size:12pt;color:rgb(128,128,128);font-family:"Times New Roman",Times,serif" color="808080"><font face="'Times New Roman', Times, serif"><span style="font-family:Cambria,Georgia,serif">)</span></font></font><br>
</div>
<span style="font-family:Cambria,Georgia,serif"></span><font color="808080"></font><span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper">
<span style="font-family:Cambria,Georgia,serif"></span><font size="3" face="'Times New Roman', Times, serif" color="808080"><span style="font-family:Cambria,Georgia,serif">CH-1015 Lausanne, Switzerland</span><br>
<span style="font-family:Cambria,Georgia,serif">+41 21 69 34 881</span></font></div>
<span style="font-family:Cambria,Georgia,serif"></span>
<div name="divtagdefaultwrapper">
<span style="font-family:Cambria,Georgia,serif"></span><a href="http://people.epfl.ch/265334" id="gmail-m_-3196822762173715662LPNoLP" target="_blank"><span style="font-family:Cambria,Georgia,serif">http://people.epfl.ch/265334</span></a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-3196822762173715662divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-espresso.org</a>> on behalf of Dr. K. C. Bhamu <<a href="mailto:kcbhamu85@gmail.com" target="_blank">kcbhamu85@gmail.com</a>><br>
<b>Sent:</b> Monday, November 15, 2021 11:47:16 AM<br>
<b>To:</b> Quantum Espresso users Forum<br>
<b>Subject:</b> [QE-users] Need some suggestion regarding the band gap issue of Fe3O4</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div style="font-size:large">Dear QE Users,<br>
</div>
<div class="gmail_quote">
<div dir="ltr">
<div style="font-size:large"><br>
</div>
<div style="font-size:large">This email may be large for you<span class="gmail_default" style="font-size:large"> all </span>as I am trying to put all the information so that I do not miss any information in your response. </div>
<div style="font-size:large">
<div style="font-size:large">I sincerely apologize for this long email.
</div>
<br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">I am dealing with <span class="gmail_default" style="font-size:large">
ferrimagnetic </span>Fe3O4<span class="gmail_default" style="font-size:large"> (with QE6.6)</span> and trying to reproduce its band gap.<span class="gmail_default" style="font-size:large"> In the literature, its experimental band gap was reported to be in the
range of 0.14-0.3eV, <a href="https://aip.scitation.org/doi/10.1063/1.5138941" target="_blank">Table SI-1 [1]</a> for the up-spin state. For Fe3O4 nanoparticles, the band gap is reported ~2eV.</span></div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">I have tried with all available PPs (mentioned below in different sets) and U values from 1eV to 7eV<span class="gmail_default" style="font-size:large">. </span>At lower U (<span class="gmail_default">=</span>1-2eV), I am getting
this system a direct band-gap half-metal (a gap is opened in the spin-up channel) and with increasing the U value, the direct band-gap (<span class="gmail_default">=</span>~1eV) nature changes to the indirect (<span class="gmail_default">=</span>~2eV) band
gap. At U=7eV with the rrkjus PPs, it showed a metallic character. </div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">I have also carried out the band structure, by keeping the lattice parameters fixed at experimental values and only relaxing ionic position, with <span class="gmail_default" style="font-size:large">varying</span> U-values. But still,
the nature of the bang-gap is similar to what I just mentioned above(<span class="gmail_default" style="font-size:large">band gap varies from
</span>~1eV to ~2eV with U=1eV to 7eV, respectively).</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">The lattice parameters, magnetic moments, and the charge on Fe_tet and Fe_oct sites I am getting reasonably in accordance with the
<a href="https://aip.scitation.org/doi/10.1063/1.5138941" target="_blank">Table SI-1</a><span class="gmail_default" style="font-size:large">.</span></div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">A sample of my input file is also shown below (the only PPs are different with respective cutoffs).</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">Information collected from other research papers:</div>
<div style="font-size:large">1. <span class="gmail_default" style="font-size:large">
</span><a href="https://www.sciencedirect.com/science/article/abs/pii/S0039602812000544" target="_blank">See [2]</a> (In section 3.1. Bulk Fe3O4: 0.33 eV band gap with VASP with U=3.8eV ). </div>
<div style="font-size:large">2. <span class="gmail_default" style="font-size:large">
</span><a href="https://iopscience.iop.org/article/10.1088/1468-6996/15/4/044202" target="_blank">See [3]</a> <font color="#000000">(Page-4, bottom left para: 2.1 eV using VASP with U=4 eV). </font><font color="#ff0000" style="font-weight:bold">This paper
only shows the </font><font style="font-weight:bold" color="#0000ff">direct band-gap of</font><font style="font-weight:bold" color="#0000ff"> 2.1eV</font><font color="#ff0000" style="font-weight:bold">.</font></div>
<div style="font-size:large">
<div style="font-size:large"></div>
<div style="font-size:large">Based on my already finished calculations and the details that I have covered above, I have a few queries:</div>
</div>
<div style="font-size:large">
<ol>
<li><span class="gmail_default" style="font-size:large">Do you think my input file (see below) has any issues? I have defined all tetra-Fe sites dn-spin and all octahedra Fe sites up-spin. </span></li><li><span class="gmail_default" style="font-size:large">What other things I should take care of to reproduce the band gap (0.14 eV to 0.3 eV)?</span></li><li><span class="gmail_default" style="font-size:large">In the <a href="https://iopscience.iop.org/article/10.1088/1468-6996/15/4/044202" target="_blank">
research paper [3]</a>, the authors used VASP and the valence configuration for the Fe atom was taken as 3d64s1. I could not find any PP from the list I mention below which has this configuration. Does this may be an issue? </span>If anyone has PP of Fe with 3d64s1<span class="gmail_default">
configuration, please provide me.</span></li><li><span class="gmail_default" style="font-size:large">As mentioned in the<b> </b><a href="https://pubs.acs.org/doi/full/10.1021/acs.jpcc.7b09387" style="font-weight:400" target="_blank">research paper [4]</a><span style="font-weight:400">, the authors have mentioned that
"</span><span class="gmail_default" style="font-weight:400;font-family:georgia,serif"></span><span style="font-weight:400;font-family:georgia,serif"><font color="#000000">Our results show that, upon releasing the symmetry constraint on the density but not
on the geometry, charge disproportionation (Fe2+/Fe3+) is observed, resulting in a band gap of around 0.2 eV at the Fermi level.</font>" How can I release symmetry constrain on the density? I used nonsym=True also but the calculation was not converging in
this case. </span></span><b></b></li></ol>
</div>
</div>
</div>
<blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px">
<blockquote style="margin:0px 0px 0px 40px;border:none;padding:0px">
<div class="gmail_quote">
<div dir="ltr">
<div style="font-size:large">
<div>
<div style="font-size:large">Supporting information for this email:</div>
</div>
</div>
</div>
</div>
</blockquote>
</blockquote>
<div class="gmail_quote">
<div dir="ltr">
<div style="font-size:large">
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><b><span class="gmail_default"></span>Combinations of my PPs:</b></div>
<br>
</div>
<div style="font-size:large">Set-1 (ecutwfc/ecutrho = 40/320Ry)<br>
</div>
<div style="font-size:large">Fe.pbe-nd-rrkjus.UPF<br>
O.pbe-rrkjus.UPF<br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">Set-2 (ecutwfc/ecutrho = 45/450Ry)</div>
<div style="font-size:large">Fe.pbe-sp-van_mit.UPF<br>
O.pbe-rrkjus.UPF<br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><b>Set-3 (ecutwfc/ecutrho = 65/785Ry) >> similar to yours</b></div>
<div style="font-size:large"><b>Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
O.pbe-n-rrkjus_psl.1.0.0.UPF</b><br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">set-4 (ecutwfc/ecutrho = 64/782Ry)</div>
<div style="font-size:large">Fe.pbe-spn-kjpaw_psl.0.2.1.UPF<br>
O.pbe-n-kjpaw_psl.0.1.UPF<br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><b>QE input file:</b></div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">&CONTROL<br>
calculation = '<span class="gmail_default" style="font-size:large">scf</span>'<br>
etot_conv_thr = 1.0000000000d-04<br>
forc_conv_thr = 1.0000000000d-03<br>
! disk_io = 'none'<br>
outdir = './tmp/'<br>
prefix = 'pwscf'<br>
pseudo_dir = '~/PPs/'<br>
<span class="gmail_default" style="font-size:large">!</span> tprnfor = .true.<br>
<span class="gmail_default" style="font-size:large">!</span> tstress = .true.<br>
! verbosity = 'high'<br>
/<br>
&SYSTEM<br>
degauss = 0.01<br>
ecutrho = 785<br>
ecutwfc = 65<br>
ibrav = 0<br>
nat = 56<br>
nosym = .false.<span class="gmail_default" style="font-size:large"> !! I have tried with .True. also but my calculations were not converging.
</span><br>
ntyp = 3<br>
occupations = 'smearing' , smearing = 'mp'<br>
starting_magnetization(1) = 1<br>
starting_magnetization(2) = -1<br>
nspin = 2<br>
lda_plus_u=.true.<br>
Hubbard_U(1)=X (tried 1eV to 7eV)<br>
Hubbard_U(2)=X (tried 1eV to 7eV)<br>
<br>
<br>
/<br>
&ELECTRONS<br>
conv_thr = 1.000000000d-6<br>
electron_maxstep = 200<br>
mixing_beta = 3.0000000000d-01<br>
/<br>
&IONS<br>
ion_dynamics = 'bfgs'<br>
/<br>
&CELL<br>
/<br>
<br>
ATOMIC_SPECIES<br>
Fe1 55.845 Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
Fe2 55.845 Fe.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
O 15.9994 O.pbe-n-rrkjus_psl.1.0.0.UPF<br>
ATOMIC_POSITIONS crystal<br>
Fe1 0.6250000000 0.6250000000 0.6250000000 <br>
Fe1 0.6250000000 0.3750000000 0.3750000000 <br>
Fe1 0.3750000000 0.6250000000 0.3750000000 <br>
Fe1 0.3750000000 0.3750000000 0.6250000000 <br>
Fe1 0.6250000000 0.8750000000 0.8750000000 <br>
Fe1 0.8750000000 0.6250000000 0.8750000000 <br>
Fe1 0.8750000000 0.8750000000 0.6250000000 <br>
Fe1 0.6250000000 0.1250000000 0.1250000000 <br>
Fe1 0.1250000000 0.6250000000 0.1250000000 <br>
Fe1 0.1250000000 0.1250000000 0.6250000000 <br>
Fe1 0.1250000000 0.3750000000 0.8750000000 <br>
Fe1 0.1250000000 0.8750000000 0.3750000000 <br>
Fe1 0.3750000000 0.1250000000 0.8750000000 <br>
Fe1 0.8750000000 0.1250000000 0.3750000000 <br>
Fe1 0.3750000000 0.8750000000 0.1250000000 <br>
Fe1 0.8750000000 0.3750000000 0.1250000000 <br>
Fe2 0.0000000000 0.0000000000 0.0000000000 <br>
Fe2 0.2500000000 0.2500000000 0.2500000000 <br>
Fe2 0.0000000000 0.5000000000 0.5000000000 <br>
Fe2 0.5000000000 0.0000000000 0.5000000000 <br>
Fe2 0.5000000000 0.5000000000 0.0000000000 <br>
Fe2 0.2500000000 0.7500000000 0.7500000000 <br>
Fe2 0.7500000000 0.2500000000 0.7500000000 <br>
Fe2 0.7500000000 0.7500000000 0.2500000000 <br>
O 0.3750000000 0.3750000000 0.3750000000 <br>
O 0.3750000000 0.6250000000 0.6250000000 <br>
O 0.6250000000 0.3750000000 0.6250000000 <br>
O 0.6250000000 0.6250000000 0.3750000000 <br>
O 0.8750000000 0.8750000000 0.8750000000 <br>
O 0.8750000000 0.6250000000 0.6250000000 <br>
O 0.6250000000 0.8750000000 0.6250000000 <br>
O 0.6250000000 0.6250000000 0.8750000000 <br>
O 0.3750000000 0.8750000000 0.8750000000 <br>
O 0.8750000000 0.3750000000 0.8750000000 <br>
O 0.8750000000 0.8750000000 0.3750000000 <br>
O 0.3750000000 0.1250000000 0.1250000000 <br>
O 0.8750000000 0.6250000000 0.1250000000 <br>
O 0.8750000000 0.1250000000 0.6250000000 <br>
O 0.6250000000 0.8750000000 0.1250000000 <br>
O 0.1250000000 0.3750000000 0.1250000000 <br>
O 0.1250000000 0.8750000000 0.6250000000 <br>
O 0.6250000000 0.1250000000 0.8750000000 <br>
O 0.1250000000 0.6250000000 0.8750000000 <br>
O 0.1250000000 0.1250000000 0.3750000000 <br>
O 0.8750000000 0.3750000000 0.3750000000 <br>
O 0.3750000000 0.8750000000 0.3750000000 <br>
O 0.3750000000 0.3750000000 0.8750000000 <br>
O 0.8750000000 0.1250000000 0.1250000000 <br>
O 0.3750000000 0.6250000000 0.1250000000 <br>
O 0.3750000000 0.1250000000 0.6250000000 <br>
O 0.6250000000 0.3750000000 0.1250000000 <br>
O 0.1250000000 0.8750000000 0.1250000000 <br>
O 0.1250000000 0.3750000000 0.6250000000 <br>
O 0.6250000000 0.1250000000 0.3750000000 <br>
O 0.1250000000 0.6250000000 0.3750000000 <br>
O 0.1250000000 0.1250000000 0.8750000000<br>
K_POINTS automatic<br>
3 3 3 0 0 0<br>
CELL_PARAMETERS angstrom<br>
8.3840000000 0.0000000000 0.0000000000<br>
0.0000000000 8.3840000000 0.0000000000<br>
0.0000000000 0.0000000000 8.3840000000<br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">Thank you very much for having patience while going through this detailed email.</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">I look forward to hearing from you.<br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">Regards</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large">
<div style="font-size:large">K C Bhamu</div>
<div style="font-size:large">Postdoctoral Fellow</div>
<div style="font-size:large">School of Chemical Engineering</div>
<div style="font-size:large">University of Ulsan</div>
<div style="font-size:large">South Korea</div>
</div>
<div style="font-size:large"><br>
</div>
<div style="font-size:large"><br>
</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div></div>