The mechanism of the HLi + NH3 -> LiNH2 + H-2 reaction has been thoroughly analyzed using topological analysis of electron density (rho) and Electron Localization Function (eta). The pathway associated with the intrinsic reaction coordinate has been calculated at the B3LYP/aug-cc-pVTZ level. Five steps have been separated on rho and eight steps on rho evolution pathway. Each step has been characterized by a decay and formation of the critical points. This analysis reveals complementarity for eta and rho fields analysis. (C) 2010 Elsevier B.V. All rights reserved.