Pseudoatom databanks, collecting parameters of multipole model of electron densities for various atom types, are used to replace Independent Atom Model by the more accurate Transferable Aspherical Atom Model (TAAM) in crystal structure refinement. They may also be used to reconstruct electron density of a molecule, crystal or biomacromolecular complex in fast yet quite accurate way and compute from it various properties, like energy of electrostatic interactions, for example. Even faster, but similarly accurate in electrostatic energy estimations model exists, the aug-PROmol. Model analogous to aug-PROmol cannot be built from the current pseudoatom databanks, as they perform badly when truncated to the monopole level. Here we searched for new strategies of multipole model refinements, leading to better parametrization already at the monopole level. This would allow to create in a single route of model parametrization a pseudoatom databank, which would be suitable for both crystal structure refinement and rapid electrostatic energy calculations. Such a route does not exist yet, because as we show here, the aug-PROmol model, alternative to the current pseudoatom databanks, is not suitable for crystal X-ray structure refinement. Here we show that cumulative approach to multipole model refinements, as oppose to simultaneous or iterative refinements of all multipole model parameters (Pv, κ, Plm, κ') leads to substantially different models of electron density. Cumulative refinement of Plm first and then κ' parameters is much worse than simultaneous. It results in electron density model giving wrong estimates of electrostatic interaction energies and atomic displacement parameters. Cumulative refinement of two blocks of parameters, Pv and κ first and then Plm and κ', on the other hand, leads to the Pv𝜅|Plm𝜅’ model having promising properties. It is similarly good as University at Buffalo DataBank (UBDB) of pseudoatoms in X-ray structure TAAM refinement and electrostatic energy estimations, especially for less polar molecules. When truncated to monopole level, the Pv𝜅 model has a chance to replace the aug-PROmol in fast yet accurate electrostatics energy calculations, although some improvements in κ parametrization for polar functional groups are still needed. The Pv𝜅|Plm𝜅’ model is also a source of point charges which behave similarly to the RESP charges in electrostatic interaction energy estimations.
