The doping of LiNbO3 with V2+, V3+, V4+ and V5+ as well as Mo3+, Mo4+, Mo5+ and Mo6+ ions is of interest in enhancing its photorefractive properties. In this paper, possible incorporation mechanisms for these ions in LiNbO3 are modelled, using a new set of interaction potentials fitted to the oxides VO, V2O3, VO2, V2O5 and to LiMoO2, Li2MoO3, LiMoO3, Li2MoO4.