<div dir="ltr"><div><div>Greetings QE developers and users,<br></div>I know this question is vague, but does anyone expect reasonable U values for iron? following the linear response method of Cococcioni and Gironcoli, I got U values for Fe in Na3Fe(II)PO4CO3 and Na2Fe(III)PO4CO3 (single Fe site in the asymmetric unit) to  be 6.39  and 7.82 eV, which seems much higher than other Fe containing compounds in phosphates (vaguely speaking). This leads to positive deviations from experimental observations (e.g. XRD cell volume and redox potential). can anyone please comment on possible causes? the input file for unperturbed SCF calculations is included (diago_thr_init = last value in zero alpha along with charge-density.dat used for subsequent non zero perturbations).<br></div>Q.E. version: 5.2.1<br><div><strong></strong><br>&CONTROL<br>                       title = 'Na2FePO4CO3 U det' ,<br>                 calculation = 'scf' ,<br>                      outdir = '/etc./common/' ,<br>                      wfcdir = '/etc./common/' ,<br>                  pseudo_dir = '/etc./pseudo/' ,<br>                      prefix = 'Na2FePO4CO3U' ,<br>                   verbosity = 'high' ,<br> /<br> &SYSTEM<br>                       ibrav = 12,<br>                           A = 8.997 ,<br>                           B = 5.163 ,<br>                           C = 6.741 ,<br>                       cosAB = -0.0027925231 ,<br>                         nat = 10,<br>                        ntyp = 5,<br>                     ecutwfc = 30 ,<br>                     ecutrho = 270 ,<br>                  lda_plus_u = .true. ,<br>             lda_plus_u_kind = 0 ,<br>                Hubbard_U(1) = 1d-20,<br>            Hubbard_alpha(1) = 1d-20,<br>                 space_group = 11 ,<br>                     uniqueb = .false. ,<br> /<br> &ELECTRONS<br>                    conv_thr = 1d-9 ,<br>                 startingpot = 'atomic' ,<br>                 startingwfc = 'atomic' ,<br>             diagonalization = 'david' <br> /<br>ATOMIC_SPECIES<br>   Fe   55.93300  Fe.pbe-sp-van_mit.UPF <br>    P   30.97400  P.pbe-van_ak.UPF <br>    C   12.01100  C.pbe-van_ak.UPF <br>    O   15.99900  O.pbe-van_ak.UPF <br>   Na   22.99000  Na.pbe-sp-van_ak.UPF <br>ATOMIC_POSITIONS crystal_sg <br>   Fe      0.362400000    0.279900000    0.250000000    <br>    P      0.587500000    0.798900000    0.250000000    <br>    O      0.569400000    0.096500000    0.250000000    <br>   Na      0.740900000    0.251300000    0.002200000    <br>    C      0.060500000    0.234500000    0.250000000    <br>    O      0.123300000    0.461600000    0.250000000    <br>    O      0.143000000    0.031700000    0.250000000    <br>    O      0.919500000    0.215700000    0.250000000    <br>    O      0.321000000    0.281700000    0.565500000    <br>    O      0.431600000    0.674000000    0.250000000    <br>K_POINTS automatic <br>  3 6 5   1 1 1 <br clear="all"><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div><div><div><font size="1"><span style="color:rgb(153,153,153)"><b>Hooman Yaghoobnejad<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>Department of Chemistry<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>Missouri University of Science and Technology<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>Rolla, MO 65401<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>USA</b></span></font><br></div></div>
</div></div></div></div></div>