dear pw users <br>i want to do a calculation of V chain&nbsp; absorb on a plane with DFT+U, but i donot know how to determine the value of U.<br>i found a PPT on the pwscf's media site, they said the value of U could be decised by a linear response approach.<br>but i really donot understand that like the input file resp_mat.in<br><br>&nbsp;&amp;input_mat<br>&nbsp;&nbsp; ntyp = 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # number of types of atoms<br>&nbsp;&nbsp; na(1) = 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # number of atoms per type<br>if we have different type atoms with different number, how we set these?<br><br>&nbsp;&nbsp; nalfa = 5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # number of perturbations (alpha)<br>&nbsp;&nbsp; filepos = 'pos_sc1'&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # file containing the atomic positions and the vectors of the unit cell<br>&nbsp;&nbsp; back = 'neutral'&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # to add a neutralizing ¡°background¡±. see PRB 71 35105 (2005)<br>what does this mean<br>&nbsp;&nbsp; filednda = 'file_sc1'&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # file containing the names of the files dn*dat <br>&nbsp;&nbsp; n1 = 4<br>&nbsp;&nbsp; n2 = 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # number of unit cells in each direction to construct the supercell for<br>&nbsp;&nbsp; n3 = 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # extrapolation (in this case a 4x4x4 supercell of the SC one)<br>&nbsp;&amp;end<br><br>and they get the U is about 2.4ev for the iron, and they use 4.3ev for the FeO calculation, <br>we always need a linear response approach to get the&nbsp; U for every system <br>or just a U for a certain pseudopotential (i mean if we get a U for V bulk, then if i could use this value for another system with V atom, like V2O3&nbsp; etc.)<br>any suggestions or comment would be OK.<br>thanks <br><br>Du<br>xiaguanying university<br><br><br><span title="neteasefooter"/></span>