There is a strict relation between the energy of hydrolysis of phosphate esters and the extent of interactions between the p ester oxygen lone pair and the antibonding orbitals of the rest of the molecule. Its impact on experimental charge density distribution is analyzed.