dear pw users <br>i want to do a calculation of V chain  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> &input_mat<br>   ntyp = 1                           # number of types of atoms<br>   na(1) = 2                         # number of atoms per type<br>if we have different type atoms with different number, how we set these?<br><br>   nalfa = 5                         # number of perturbations (alpha)<br>   filepos = 'pos_sc1'          # file containing the atomic positions and the vectors of the unit cell<br>   back = 'neutral'               # to add a neutralizing ¡°background¡±. see PRB 71 35105 (2005)<br>what does this mean<br>   filednda = 'file_sc1'        # file containing the names of the files dn*dat <br>   n1 = 4<br>   n2 = 4                             # number of unit cells in each direction to construct the supercell for<br>   n3 = 4                             # extrapolation (in this case a 4x4x4 supercell of the SC one)<br> &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  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  etc.)<br>any suggestions or comment would be OK.<br>thanks <br><br>Du<br>xiaguanying university<br><br><br><span title="neteasefooter"/></span>