Gauge-Invariant Formulation of Circular Dichroism.

Standard formulations of magnetic response properties, such as circular dichroism spectra, are plagued by gauge dependencies, which can lead to unphysical results. In this work, we present a general gauge-invariant and numerically efficient approach for the calculation of circular dichroism spectra from the current density. First we show that in this formulation the optical rotation tensor, the response function from which circular dichroism spectra can be obtained, is independent of the origin of the coordinate system. We then demonstrate that its trace is independent of the gauge origin of the vector potential. We also show how gauge invariance can be retained in practical calculations with finite basis sets. As an example, we explain how our method can be applied to time-dependent current-density-functional theory. Finally, we report gauge-invariant circular dichroism spectra obtained using the adiabatic local-density approximation. The circular dichroism spectra we thus obtain are in good agreement with experiment.

Gauge-invariant calculation of static and dynamical magnetic properties from the current density.

In this work we solve two problems related to the calculation of static and dynamical magnetic properties with ab initio theories. First, we show that the dependence of the dynamical magnetic dipole moment on the reference point of the multipole expansion and on the gauge origin of the vector potential have a clear physical significance. They are due to a dynamical electric dipole moment and an electric field, respectively. Both are fully determined by the experimental setup and do not pose any fundamental problem, contrary to what is commonly assumed. Second, in the static case, any dependence on the gauge origin is an artifact of the computational method. We show that the artificial dependence on the gauge origin can be removed in an elegant way by the introduction of a sum rule that puts the diamagnetic and paramagnetic contributions on equal footing. Our approach can be applied to calculate any magnetic observable that can be derived from the current density, and can be used in combination with any ab initio theory from which it can be obtained. To illustrate our method we apply it here to time-dependent current-density-functional theory for the calculation of static and dynamical magnetizabilities of molecules.