We attribute the discrepancy between the two sets of computed anisotropy trends to the sensitivity of PADs to the diffuse parts of the wave functions. The importance of long-range wave function behavior in PAD calculations has been stressed before; it comes from the dipole operator 𝑟, which amplifies the contribution of the tails of Dyson orbitals into the dipole matrix elements.