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

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Origin Independent Current Density Vector Fields Induced by Time-Dependent Magnetic Field. I. The LiH molecule

The electronic current density, induced in a molecule by the optical magnetic field associated to a frequency dependent monochromatic plane wave, assumed to be spatially uniform within the electric quadrupole approximation, has been studied by a theoretical method based on continuous translation of its origin. The induced electronic current density vector field designated by this procedure, invariant of the origin for any point of the molecular domain, is obtained via a computational scheme formally annihilating the diamagnetic contribution of the conventional common-origin approach based on perturbation theory. In a preliminary application of the theoretical methods outlined in the present work, the simple molecule of lithium hydride has been investigated. Particular attention has been paid to the structure of induced electronic current density for several values of the magnetic field frequency, by investigating equilibrium points of four different types, organized in stagnation lines which constitute its stagnation graph, i.e., a topological fingerprint of the vector field conveying complete information in condensed form, to verify the fulfillment of fundamental requirements, e.g., the continuity equation and the Poincaré-Hopf theorem on spherical and toroidal surfaces.

Electronic Currents Induced by Optical Fields and Rotatory Power Density in Chiral Molecules

The electric dipole–magnetic dipole polarizability tensor κ′, introduced to interpret the optical activity of chiral molecules, has been expressed in terms of a series of density functions kαβ′, which can be integrated all over the three-dimensional space to evaluate components καβ′ and trace καα′. A computational approach to kαβ′, based on frequency-dependent electronic current densities induced by monochromatic light shining on a probe molecule, has been developed. The dependence of kαβ′ on the origin of the coordinate system has been investigated in connection with the corresponding change of καβ′. It is shown that only the trace kαα′ of the density function defined via dynamic current density evaluated using the continuous translation of the origin of the coordinate system is invariant of the origin. Accordingly, this function is recommended as a tool that is quite useful for determining the molecular domains that determine optical activity to a major extent. A series of computations on the hydrogen peroxide molecule, for a number of different HO–OH dihedral angles, is shown to provide a pictorial documentation of the proposed method.

Frequency-dependent current density tensors as density functions of dynamic polarizabilities

Relationships accounting for contributions to the first-order charge density ρ⁽¹⁾ and current density J⁽¹⁾, induced in the electrons of a molecule by a monochromatic plane wave, have been obtained via time-dependent quantum mechanical perturbation theory. Their gauge invariance and invariance in passive translations of the coordinate system have been demonstrated (i) within the (long wavelengths) electric dipole approximation, in which only the electric field contributions to these densities are required, and (ii) within the electric quadrupole approximation, in which nonseparable terms provided by the magnetic field and by the electric field gradient, assumed uniform over the molecular domain, are needed. It is shown that the physical meaning of current density tensors depending on the frequency ω of the monochromatic wave shone on the molecule, and corresponding to derivatives of the current density J⁽¹⁾(r, ω) with respect to components of the perturbing fields, is that of property density. Therefore, frequency-dependent current density tensors can be interpreted as the integrand function in 3d-space integrals defining dynamic molecular response tensors, e.g., electric dipole, mixed electric dipole-magnetic dipole, and electric dipole-quadrupole polarizabilites. Plots of current density tensors are expected to provide important information on the molecular domains giving dominant contributions to these properties.

Gauge invariance and origin independence of electronic charge density and current density induced by optical fields

Expressions for the first-order polarization charge density ρ ⁽¹⁾ and current density J ⁽¹⁾ induced in a molecule by a monochromatic plane wave, obtained by time-dependent quantum mechanical perturbation theory, have been investigated to assess their gauge invariance and independence of the coordinate system in passive and active translations. The conditions arrived at show that, within the (long wavelengths) dipole approximation, only the electric contributions to these densities are needed to rationalize the phenomenology. To the next higher quadrupole approximation, assuming that the magnetic field and the electric field gradient are uniform over the molecular dimensions, corresponding contributions to ρ ⁽¹⁾ and J ⁽¹⁾ are considered. It has been found that total densities are independent of the origin, whereas the contributions from electric and magnetic fields are not separately invariant. A magnetic contribution to J ⁽¹⁾, which is by itself origin independent, can be defined by means of an approach based on continuous translation of the origin of the coordinate system.

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.