is also origin invariant. The LG(OI) transformation is a post-processing treatment of the β tensor that is trivial to implement and apply. In HF or DFT, one needs only to calculate the perturbed density for the µ R operator and then contract it with the appropriate multipole integrals to obtain all relevant tensors: the standard electric dipole-electric dipole polarizability in the length and mixed gauges, the electric dipolemagnetic dipole polarizability β, and even the electric dipole-electric quadrupole polarizability necessary to calculate the full OR tensor for oriented systems. In this sense, the LG(OI)-OR is a free subproduct of a standard electric polarizability calculation. For CC methods, our current implementation with the symmetric form of the LR function requires an explicit evaluation of the µ P perturbed densities. However, a non-symmetric formulation of the LR function would only require the evaluation of the µ R perturbed amplitudes, similar to the DFT case.