“…The total quadrupole moment tensor Q νµ = 1 2 ρ(r)r ν r µ dr (ρ(r) is the total charge density of the system) and the terms χ Q1 , χ Q2 and χ Qs are well defined and origin independent only for a system with zero dipole and zero χ D , respectively. Hence to properly compute quadruple contributions from isotropic bulk regions, we used the total charge density of the supercell; we avoided summations over quadrupole contributions associated to single polar molecular units, at variance from several formulations presented in the literature [12,14]. We emphasize that it is key to use properties of the entire system, such as the total charge density, to evaluate χ Q1 , χ Q2 and χ Qs , thus insuring that, as long as the entire sample under consideration is centrosymmetric (χ D = 0), the results are origin independent.…”