and then combined as in Eqs. 1 and 7 to obtain the B tensor. Note that both expressions in Eqs. 34-35 are missing a term that depends on the k-derivative of the perturbed density. Based on the results in section IV, this contribution seems numerically small. Although the tensors calculated within the MVG approach are origin invariant, this is not directly the case with the LG approach. However, this issue is resolved by applying the LG(OI) transformation, 2,3 which avoids the need to employ gauge-including atomic orbitals. This approach is based on the singular value decomposition (SVD) of the mixed-gauge electric dipole-electric dipole polarizability and the application of the inverse transformation to the β and A tensors. The extra tensor can be readily computed by combining P µ with the µ V integrals and the SVD of a square matrix of dimension 3 is quick. Therefore, the cost of the LG(OI) tranformation is negligible and this approach is computationally more efficient than MVG because the evaluation of the tensors in the static limit is not required.