Hellmann-Feynman Forces within the DFT+U in Wannier functions basis
  The most general way to describe localized atomic-like electronic states in strongly correlated compounds is to utilize Wannier functions. In the present paper we continue the development of widely-spread DFT+U method onto Wannier function basis set and propose the technique to calculate the Hubbard contribution to the forces. The technique was implemented as a part of plane-waves pseudopotential code Quantum-ESPRESSO and successfully tested on a charge transfer insulator NiO. 
D. Novoselov, Dm. M. Korotin, V.I. Anisimov
cond-mat.str-el