Can a multipole analysis faithfully reproduce topological descriptors of a total charge density?

Total charge densities ρ(r) of solid NH3 have been derived using an ab initio crystalline molecular-orbital approach and also from multipole refinement of the structure factors obtained from the same charge density. Comparison of the topological features of these charge densities, as defined by the quantum theory of atoms in molecules, has been used to probe the ability of the multipole analysis to reproduce exactly known total charge-density distributions. For the most part, multipole refinement satisfactorily returns the features of the original density, although the fit to theoretical data is not as good as that to the experimental data. The one topological parameter that is poorly reproduced is the Laplacian \nabla^2\rho({\bf r}_b) at NH bond critical points.