What these profiles also demonstrate is the lack of a need to include lone-pair repulsion sites when computing E Pauli . If instead we were using a model that was some increasing function of electron density overlap, then LP sites would also contribute to an increase in Pauli repulsion. As the slopes of AmberFD match very nicely to that of the DFT energies, a spherical LP-less representation is a sufficient substitute.