The issues raised in the comment by T.A. Manz are addressed through the presentation of calculated atomic charges for NaF, NaCl, MgO, SrTiO 3 and La 2 Ce 2 O 7 , using our previously presented method for calculating Hirshfeld-I charges in Solids [J. Comput. Chem.. doi: 10.1002/jcc.23088]. It is shown that the use of pseudo-valence charges is sufficient to retrieve the full all-electron Hirshfeld-I charges to good accuracy. Furthermore, we present timing results of different systems, containing up to over 200 atoms, underlining the relatively low cost for large systems. A number of theoretical issues is formulated, pointing out mainly that care must be taken when deriving new atoms in molecules methods based on "expectations" for atomic charges.
PACS numbers:Our recent extension of the Hirshfeld-I method to solids and periodic systems 1 , allows for the calculation of atoms in molecules (AIM) in solid state codes using common techniques as pseudo-potentials and plane wave basis sets. As a Hirshfeld-type method, 2 it relies on atomic reference densities that are used to construct AIM weight functions w(r) that allow to extract an AIM density function from a computed system electron density distribution (EDD). As was shown by us, when using pseudo-potentials and plane waves, some issues appear that require attention when generating atomic reference densities.1 However, we showed that the delocalization problem can be handled, and that it is actually part of the larger conceptual problem of defining reference densities for anions. In addition, we showed that the use of pseudo-potentials, and their resulting pseudo-valence charges, can provide the all-electron values for the atomic charges. Finally, we showed that the method scales favorably for large systems.In a comment to our paper, 3 T. A. Manz raises some questions which we address here in detail. His criticisms and suggestions can be summarized as:1. Hirshfeld-I atomic charges do not give "chemically feasible" atomic charges and improved performance of Atoms in Molecules (AIM) methods may be obtained by combining spherical averaging and Hirshfeld-I methods. 3. The presented computational scaling was insufficiently justified.In the following paragraphs we address each of these comments, from a theoretical as well as a computational perspective.