We performed Hartree-Fock (HF) calculations with the aug-cc-pVDZ basis set 32 with a 589.3 nm perturbation wavelength to avoid resonance conditions. It is not our intention to compute experimentally accurate OR, but rather simply to demonstrate the origin invariance of the SAvg ia and SHemi ia definitions, see Eqs. 12 and 14, respectively. Furthermore, the implementation of the SHemi ia definition is based on a simple Python code that forms the entire response matrix in Eq. 4, then applies the Cholesky decomposition in Eq. 13 to form SHemi ia . The Python code reads two-electron integrals and orbital energies from a modified version of the GAUSSIAN suite of programs to form the appropriate matrices and it is very time and memory intensive. Consequently, calculations are performed at HF (rather than density functional theory) level to avoid the need to also read in 4-center, 1electron exchange-correlation integrals when computing SHemi ia . However, the results obtained here equally apply to calculations using approximate density functionals. All OR and SAvg ia calculations were performed with the same development version of GAUSSIAN .