Non-perturbative calculation of orbital- and spin effects in molecules subject to non-uniform magnetic fields

External non-uniform magnetic fields acting on molecules induce non-collinear spin-densities and spin-symmetry breaking. This necessitates a general two-component Pauli spinor representation. In this paper, we report the implementation of a General Hartree-Fock method, without any spin constraints, for non-perturbative calculations with finite non-uniform fields. London atomic orbitals are used to ensure faster basis convergence as well as invariance under constant gauge shifts of the magnetic vector potential. The implementation has been applied to an investigate the joint orbital and spin response to a field gradient---quantified through the anapole moments---of a set of small molecules placed in a linearly varying magnetic field. The relative contributions of orbital and spin-Zeeman interaction terms have been studied both theoretically and computationally. Spin effects are stronger and show a general paramagnetic behaviour for closed shell molecules while orbital effects can have either direction. Basis set convergence and size effects of anapole susceptibility tensors have been reported. The relation of the mixed anapole susceptibility tensor to chirality is also demonstrated

Guaranteed Convergence of a Regularized Kohn-Sham Iteration in Finite Dimensions

The exact Kohn-Sham iteration of generalized density-functional theory in finite dimensions witha Moreau-Yosida regularized universal Lieb functional and an adaptive damping step is shown toconverge to the correct ground-state density.Comment: 3 figures, contains erratum with additional author Paul E. Lammer

Density-Wavefunction Mapping in Degenerate Current-Density-Functional Theory

We show that the particle density, $\rho(\mathbf{r})$, and the paramagnetic current density, $\mathbf{j}^{p}(\mathbf{r})$, are not sufficient to determine the set of degenerate ground-state wave functions. This is a general feature of degenerate systems where the degenerate states have different angular momenta. We provide a general strategy for constructing Hamiltonians that share the same ground state density, yet differ in degree of degeneracy. We then provide a fully analytical example for a noninteracting system subject to electrostatic potentials and uniform magnetic fields. Moreover, we prove that when $(\rho,\mathbf{j}^p)$ is ensemble $(v,\mathbf{A})$-representable by a mixed state formed from $r$ degenerate ground states, then any Hamiltonian $H(v',\mathbf{A}')$ that shares this ground state density pair must have at least $r$ degenerate ground states in common with $H(v,\mathbf{A})$. Thus, any set of Hamiltonians that shares a ground-state density pair $(\rho,\mathbf{j}^p)$ by necessity has at least have one joint ground state

A local tensor that unifies kinetic energy density and vorticity dependent exchange-correlation functionals

We present a kinetic energy tensor that unifies a scalar kinetic energy density commonly used in meta-Generalized Gradient Approximation functionals and the vorticity density that appears in paramagnetic current-density-functional theory. Both types of functionals can thus be subsumed as special cases of a novel functional form that is naturally placed on the third rung of Jacob's ladder. Moreover, the kinetic energy tensor is related to the exchange hole curvature, is gauge invariant, and has very clearcut $N$-representability conditions. The latter conditions enable the definition of effective number of non-negligible orbitals. Whereas quantities such as the Electron Localization Function can discriminate effective one-orbital regions from other regions, the present kinetic energy tensor can discriminate between one-, two-, three-, and four-or-more orbital regions

Revisiting density-functional theory of the total current density

Density-functional theory requires an extra variable besides the electron density in order to properly incorporate magnetic-field effects. In a time-dependent setting, the gauge-invariant, total current density takes that role. A peculiar feature of the static ground-state setting is, however, that the gauge-dependent paramagnetic current density appears as the additional variable instead. An alternative, exact reformulation in terms of the total current density has long been sought but to date a work by Diener is the only available candidate. In that work, an unorthodox variational principle was used to establish a ground-state density-functional theory of the total current density as well as an accompanying Hohenberg-Kohn-like result. We here reinterpret and clarify Diener's formulation based on a maximin variational principle. Using simple facts about convexity implied by the resulting variational expressions, we prove that Diener's formulation is unfortunately not capable of reproducing the correct ground-state energy and, furthermore, that the suggested construction of a Hohenberg-Kohn map contains an irreparable mistake

Kohn-Sham theory with paramagnetic currents: compatibility and functional differentiability

Recent work has established Moreau-Yosida regularization as a mathematical tool to achieve rigorous functional differentiability in density-functional theory. In this article, we extend this tool to paramagnetic current-density-functional theory, the most common density-functional framework for magnetic field effects. The extension includes a well-defined Kohn-Sham iteration scheme with a partial convergence result. To this end, we rely on a formulation of Moreau-Yosida regularization for reflexive and strictly convex function spaces. The optimal $L^p$-characterization of the paramagnetic current density $L^1\cap L^{3/2}$ is derived from the $N$-representability conditions. A crucial prerequisite for the convex formulation of paramagnetic current-density-functional theory, termed compatibility between function spaces for the particle density and the current density, is pointed out and analyzed. Several results about compatible function spaces are given, including their recursive construction. The regularized, exact functionals are calculated numerically for a Kohn-Sham iteration on a quantum ring, illustrating their performance for different regularization parameters