Dear All,<br><br>i would like to calculate the U and J for a given material Fe based.<br>I can do LDA+U calculations without problems is such system using PW.<br><br>The paper "Phys. Rev. B
                      71,
                      035105
                      (2005)" well explains how to calculate the parameter U.<br>i have however troubles in setting the input files to do the calculations.<br>could some examples be provided?<br><br>I have also another question.<br>
<br>In the current distribution of PW in the directory PP there is a file wannier_ham.f90.<br>I believe it is referring to the paper:<br><br><a href="http://xxx.lanl.gov/pdf/0801.3500">http://xxx.lanl.gov/pdf/0801.3500</a><br>
<br>The subroutines are reading wfs and eigenvalues of a given self-consistent calculation and they calculate the Hamiltonian in Wannier basis set.<br>This gives the possibility to calculate on-site energy of given d states.<br>
If now we change the occupancy at a given site by the flag Hubbard_alpha() and we recalculate the Hamiltonian in Wannier basis set for such self-consistent calculation, the variation of eigenvalues respect to the occupations should give the U and J parameters.<br>
Is this correct?<br>Is this the procedure used in the mentioned paper?<br><br>thank you.<br><br>cheers,<br><br>Gianluca<br><br>