The density of states and its projection onto the atomic orbitals have been calculated through the periodic LCAO approach, using the CRYSTAL-06 code ^{153} and the same DFT functional ^{43} as used in all calculations described here. Oxygen atoms were described by an 8-411+*d* all electron basis set ^{154} and Xe atoms by a cc-pVTZ basis set ^{155}, specifically devised for the fully-relativistic ECP28MDF pseudopotential, with 26 electrons per Xe kept active. In order to avoid problems of numerical catastrophes ^{153} , the original cc-pVTZ(12*s*11*p*9*d*1*f*)/[5*s*4*p*3*d*1*f*] basis set ^{155} was slightly modified by removing the outermost *s*, *p* and *d* Gaussians and the *f* polarization function to yield a final (11*s*10*p*8*d*)/[5*s*4*p*3*d*] contracted basis set. The Kohn-Sham matrix was diagonalized on an isotropic 8 8 8 k-mesh and the same mesh was used in the Fermi energy calculation and density matrix reconstruction. Baderâ€™s charges and atomic quadrupole moment tensors were evaluated from the periodic LCAO electron density and using the TOPOND package ^{156} interfaced to the CRYSTAL code. TOPOND determines the boundaries of the atomic basins and the integrated properties within these basins using fully analytical and on-the-fly evaluations of the electron densities and its derivatives (up to the 4th order), i.e. no use is made of electron densities on a grid or of numerical approximations of electron density gradients ^{157}. Basin boundaries are determined through the PROMEGA algorithm ^{158}, while the basin integration is performed in spherical coordinates, using Gaussian quadrature formulas. The resulting atomic charges were very close to those obtained from the PAW wavefunction and the numerical algorithm described in ^{137}. Band gaps were recalculated within *GW* approximation as implemented in the VASP code ^{159}.