“…This approximation fails, e.g., for high-energy photons that have short wavelengths. In , this was originally addressed through a complete second-order multipole expansion, which was then applied to high-energy X-ray absorption and scattering processes. − However, the multipole expansion itself does not necessarily have a smooth convergence behavior toward the exact result, and is not origin independent unless using the correct length and velocity gauges. , In contrast, the plane-waveform of the wave vector, i.e., the exact semiclassical light–matter interaction operator, shows excellent stability also for small basis sets. ,− In , the operator has been implemented using the Gauss–Hermite quadrature, which makes it easy to implement both isotropic averages and defined directions of wave and polarization vectors. , This implementation has also been extended to circularly polarized light, allowing the computation of rotatory strengths and tensors beyond the dipole approximation …”