In this paper we discuss possible mechanisms of PV annealing in Si. Our approach includes a combination of density functional theory and lattice kinetic Monte-Carlo (LKMC) simulations. The density functional theory is used to find the binding energies and jump barriers for P-V pair at different separations (from one to three interatomic bonds between complex constituents) and in different charge states. The mobility of the complex is simulated by LKMC with event probabilities calculated based on the energies from ab-initio calculations. .