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.
