Search

Non-adiabatic coupling matrix elements (NACMEs) are important in quantum chemistry, particularly for molecular dynamics methods such as surface hopping. However, NACMEs are gauge dependent. This presents a difficulty for their calculation in general, where there are no restrictions on the gauge function except that it be differentiable. These cases are relevant for complex-valued electronic wave functions, such as those that arise in the presence of a magnetic field or spin-orbit coupling. Additionally, the Berry curvature and Berry force play an important role in molecular dynamics in a magnetic field, and are also relevant in the context of spin-orbit coupling. For methods such as surface hopping, excited-state Berry curvatures will also be of interest. With this in mind, we have developed a scheme for the calculation of continuous, differentiable NACMEs as a function of the molecular geometry for complex-valued wave functions. We demonstrate the efficacy of the method by using the H$_2$ molecule at the full configuration-interaction (FCI) level of theory. Additionally, ground- and excited- state Berry curvatures are computed for the first time using FCI theory. Finally, Berry phases are computed directly in terms of diagonal NACMEs.

The exchange-only virial relation due to Levy and Perdew is revisited. Invoking the adiabatic connection, we introduce the exchange energy in terms of the right-derivative of the universal density functional w.r.t. the coupling strength $\lambda$ at $\lambda=0$. This agrees with the Levy-Perdew definition of the exchange energy as a high-density limit of the full exchange-correlation energy. By relying on $v$-representability for a fixed density at varying coupling strength, we prove an exchange-only virial relation without an explicit local-exchange potential. Instead, the relation is in terms of a limit ($\lambda \searrow 0$) involving the exchange-correlation potential $v_\mathrm{xc}^\lambda$, which exists by assumption of $v$-representability. On the other hand, a local-exchange potential $v_\mathrm{x}$ is not warranted to exist as such a limit.

Ab initio molecular dynamics in a magnetic field requires solving equations of motion with velocity-dependent forces -- namely, the Lorentz force arising from the nuclear charges moving in a magnetic field and the Berry force arising from the shielding of these charges from the magnetic field by the surrounding electrons. In this work, we revisit two existing propagators for these equations of motion, the auxiliary-coordinates-and-momenta (ACM) propagator and the Tajima propagator (TAJ), and compare them with a new exponential (EXP) propagator based on the Magnus expansion. Additionally, we explore limits (for example, the zero-shielding limit), the implementation of higher-order integration schemes, and series truncation to reduce computational cost by carrying out simulations of a HeH$^+$ model system for a wide range of field strengths. While being as efficient as the TAJ propagator, the EXP propagator is the only propagator that converges to both the schemes of Spreiter and Walter (derived for systems without shielding of the Lorentz force) and to the exact cyclotronic motion of a charged particle. Since it also performs best in our model simulations, we conclude that the EXP propagator is the recommended propagator for molecules in magnetic fields.

The Born-Oppenheimer (BO) approximation is less accurate in the presence of a strong magnetic field than in the absence of a field. This is due to the complicated and unpredictable response of electronic structure to the field, especially in the mixed regime $B \approx B_0 = 2.35 \times 10^5\,$T. Therefore, it is desirable to explore non-BO methods in magnetic fields. In this work, the nuclear-electronic orbital (NEO) method is employed to study protonic vibrational excitation energies in the presence of a strong magnetic field. NEO Generalized Hartree-Fock theory and time-dependent Hartree-Fock theory are derived and implemented, accounting for all terms that result as a consequence of the nonperturbative treatment of molecular systems in a magnetic field. The NEO results for HCN and FHF$^-$ with clamped heavy nuclei are compared against the quadratic eigenvalue problem (QEP). Each molecule has three semi-classical modes owing to the hydrogen - two precession modes that are degenerate in the absence of a field and one stretching mode. The NEO-TDHF model is found to perform well - in particular it automatically captures the screening effects of the electrons on the nuclei, which are quantified through the difference in energy of the precession modes.

A semiclassical theory of small oscillations is developed for nuclei that are subject to velocity-dependent forces in addition to the usual interatomic forces. When the velocity-dependent forces are due to a strong magnetic field, novel effects arise -- for example, the coupling of vibrational, rotational, and translational modes. The theory is first developed using Newtonian mechanics and we provide a simple quantification of the coupling between these types of modes. We also discuss the mathematical structure of the problem, which turns out to be a quadratic eigenvalue problem rather than a standard eigenvalue problem. The theory is then re-derived using the Hamiltonian formalism, which brings additional insight, including a close analogy to the quantum-mechanical treatment of the problem. Finally, we provide numerical examples for the H$_2$, HT, and HCN molecules in a strong magnetic field.

The Berry curvature is essential in Born$-$Oppenheimer molecular dynamics, describing the screening of the nuclei by the electrons in a magnetic field. Parts of the Berry curvature can be understood as the external magnetic field multiplied by an effective charge so that the resulting Berry force behaves like a Lorentz force during the simulations. Here we investigate whether these effective charges can provide insight into the electronic structure of a given molecule or, in other words, whether we can perform a population analysis based on the Berry curvature. To develop our approach, we first rewrite the Berry curvature in terms of charges that partially capture the effective charges and their dependence on the nuclear velocities. With these Berry charges and charge fluctuations, we then construct our population analysis yielding atomic charges and overlap populations. Calculations at the Hartree$-$Fock level reveal that the atomic charges are similar to those obtained from atomic polar tensors. However, since we additionally obtain an estimate for the fluctuations of the charges and a partitioning of the atomic charges into contributions from all atoms, we conclude that the Berry population analysis is a useful alternative tool to analyze the electronic structure of molecules.

The Berry connection and curvature are key components of electronic structure calculations for atoms and molecules in magnetic fields. They ensure the correct translational behavior of the effective nuclear Hamiltonian and the correct center-of-mass motion during molecular dynamics in these environments. In this work, we demonstrate how these properties of the Berry connection and curvature arise from the translational symmetry of the electronic wave function and how they are fully captured by a finite basis set of London orbitals but not by standard Gaussian basis sets. This is illustrated by a series of Hartree-Fock calculations on small molecules in different basis sets. Based on the resulting physical interpretation of the Berry curvature as the shielding of the nuclei by the electrons, we introduce and test a series of approximations using the Mulliken fragmentation scheme of the electron density. These approximations will be particularly useful in ab initio molecular dynamics calculations in a magnetic field, since they reduce the computational cost, while recovering the correct physics and up to 95% of the exact Berry curvature.

Lieb's convex formulation of density-functional theory is presented in a pedagogical manner, emphasizing its connection to Hohenberg-Kohn theory and to Levy's constrained-search theory. The Hohenberg-Kohn and Lieb variation principles are discussed, highlighting the dual relationship between the ground-state energy and the universal density functional. Applications of the Lieb variation principle are reviewed, demonstrating how it may be utilized to calculate the Kohn-Sham potential of atoms and molecules, to study the exchange-correlation functional and the adiabatic connection by high-precision many-body methods, and to calculate the exchange-correlation hole and energy densities of atoms and molecules., Comment: Submitted to the book in honour of Elliott Lieb's 90th birthday, edited by Rupert L. Frank, Ari Laptev, Mathieu Lewin and Robert Seiringer, published by EMS Press

The diagonal nonadiabatic term arising from the Born--Oppenheimer wave-function ansatz contains contributions from a vector and scalar potential. The former is provably zero when the wave function can be taken to be real valued, and the latter, known as the diagonal Born--Oppenheimer correction (DBOC), is typically small in magnitude. Therefore, unless high accuracy is sought, the diagonal nonadiabatic term is usually neglected when calculating molecular properties. In the presence of a magnetic field, the wave function is generally complex, and the geometric vector potential gives rise to a screening force that is qualitatively important for molecular dynamics. This screening force is written in terms of the Berry curvature and is added to the bare Lorentz force acting on the nuclei in the presence of the field. In this work, we derive analytic expressions for the Berry curvature and DBOC using both first and second quantization formalisms for the case of generalized and restricted Hartree--Fock theories in a uniform magnetic field. The Berry curvature and DBOC are calculated as a function of the magnetic field strength and the bond distance for the ground-state singlets of H$_2$, LiH, BH, and CH$^+$. We also examine the stability and time-reversal symmetry of the underlying self-consistent field solutions. The character of the DBOC and Berry curvature is found to depend upon the magnetic field and varies between molecules. We also identify instances of broken time-reversal symmetry for the dissociation curves of BH and CH$^+$., Comment: 29 Pages, 14 Figures, 34 individual figures (figure panels)

Strong magnetic fields have a large impact on the dynamics of molecules. In addition to the changes of the electronic structure, the nuclei are exposed to the Lorentz force with the magnetic field being screened by the electrons. In this work, we explore these effects using ab-initio molecular dynamics simulations based on an effective Hamiltonian calculated at the Hartree-Fock level of theory. To correctly include these non-conservative forces in the dynamics, we have designed a series of novel propagators that show both good efficiency and stability in test cases. As a first application, we analyze simulations of He and H$_2$ at two field strengths characteristic of magnetic white dwarfs (0.1 $B_0 = 2.35 \times 10^4$ T and $B_0 = 2.35 \times 10^5$ T). While the He simulations clearly demonstrate the importance of electron screening of the Lorentz force in the dynamics, the extracted rovibrational spectra of H$_2$ reveal a number of fascinating features not observed in the field-free case: couplings of rotations/vibrations with the cyclotron rotation, overtones with unusual selection rules, and hindered rotations that transmute into librations with increasing field strength. We conclude that our presented framework is a powerful tool to investigate molecules in these extreme environments.

The dynamics of a molecule in a magnetic field is significantly different form its zero-field counterpart. One important difference in the presence of a field is the Lorentz force acting on the nuclei, which can be decomposed as the sum of the bare nuclear Lorentz force and a screening force due to the electrons. This screening force is calculated from the Berry curvature and can change the dynamics qualitatively. It is therefore important to include the contributions from the Berry curvature in molecular dynamics simulations in a magnetic field. In this work, we present a scheme for calculating the Berry curvature numerically, by a finite-difference technique, addressing challenges related to the arbitrary global phase of the wave function. The Berry curvature is calculated as a function of bond distance for H$_2$ at the restricted and unrestricted Hartree--Fock levels of theory and for CH$^{+}$ as a function of the magnetic field strength at the restricted Hartree--Fock level of theory. The calculations are carried out using basis sets of contracted Gaussian functions equipped with London phase factors (London orbitals) to ensure gauge-origin invariance. In the paper, we also interpret the Berry curvature in terms of atomic charges and discuss its convergence in basis sets with and without London phase factors. Calculation of the Berry curvature allows for its inclusion in \textit{ab initio} molecular dynamics simulations in a magnetic field.

We present a Gaussian-basis implementation of orbital-free density-functional theory (OF-DFT) in which the trust-region image method (TRIM) is used for optimization. This second-order optimization scheme has been constructed to provide benchmark all-electron results with very tight convergence of the particle number constraint, associated chemical potential and electron density. It is demonstrated that, by preserving the saddle-point nature of the optimization and simultaneously optimizing the density and chemical potential, an order of magnitude reduction in the number of iterations required for convergence is obtained. The approach is compared and contrasted with a new implementation of the nested optimization scheme put forward by Chan, Cohen and Handy. Our implementation allows for semi-local kinetic-energy (and exchange-correlation) functionals to be handled self-consistently in all-electron calculations. The all-electron Gaussian-basis setting for these calculations will enable direct comparison with a wide range of standard high-accuracy quantum-chemical methods as well as with Kohn-Sham density-functional theory. We expect that the present implementation will provide a useful tool for analysing the performance of approximate kinetic-energy functionals in finite systems.

A cornerstone of current-density functional theory (CDFT) in its paramagnetic formulation is proven. After a brief outline of the mathematical structure of CDFT, the lower semi-continuity and expectation valuedness of the CDFT constrained-search functional is proven, meaning that there is always a minimizing density matrix in the CDFT constrained-search universal density functional. These results place the mathematical framework of CDFT on the same footing as that of standard DFT.

We investigate the helium dimer in strong magnetic fields, focusing on the spectrum of low-lying electronic states and their dissociation curves, at the full configuration-interaction level of theory. To address the loss of cylindrical symmetry and angular momentum as a good quantum number for nontrivial angles between the bond axis and magnetic field, we introduce the almost quantized angular momentum (AQAM) and show that it provides useful information about states in arbitrary orientations. In general, strong magnetic fields dramatically rearrange the spectrum, with the orbital Zeeman effect bringing down states of higher angular momentum below the states with pure $\sigma$ character as the field strength increases. In addition, the spin Zeeman effect pushes triplet states below the lowest singlet; in particular, a field of one atomic unit is strong enough to push a quintet state below the triplets. In general, the angle between the bond axis and the magnetic field also continuously modulates the degree of $\sigma$, $\pi$, and $\delta$ character of bonds and the previously identified perpendicular paramagnetic bonding mechanism is found to be common among excited states. Electronic states with preferred skew field orientations are identified and rationalized in terms of permanent and induced electronic currents.

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.

A detailed account of the Kohn-Sham algorithm from quantum chemistry, formulated rigorously in the very general setting of convex analysis on Banach spaces, is given here. Starting from a Levy-Lieb-type functional, its convex and lower semi-continuous extension is regularized to obtain differentiability. This extra layer allows to rigorously introduce, in contrast to the common unregularized approach, a well-defined Kohn-Sham iteration scheme. Convergence in a weak sense is then proven. This generalized formulation is applicable to a wide range of different density-functional theories and possibly even to models outside of quantum mechanics.

A G{\"o}rling-Levy (GL)-based perturbation theory along the range-separated adiabatic connection is assessed for the calculation of electronic excitation energies. In comparison with the Rayleigh-Schr{\"o}dinger (RS)-based perturbation theory introduced in a previous work [E. Rebolini, J. Toulouse, A. M. Teale, T. Helgaker, A. Savin, Mol. Phys. 113, 1740 (2015)], this GL-based perturbation theory keeps the ground-state density constant at each order and thus gives the correct ionization energy at each order. Excitation energies up to first order in the perturbation have been calculated numerically for the helium and beryllium atoms and the hydrogen molecule without introducing any density-functional approximations. In comparison with the RS-based perturbation theory, the present GL-based perturbation theory gives much more accurate excitation energies for Rydberg states but similar excitation energies for valence states.

We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre--Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogues in LDFT. We prove results concerning $N$-representability, Hohenberg--Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analogue to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed.

The relationship between the densities of ground-state wave functions (i.e., the minimizers of the Rayleigh--Ritz (RR) variation principle) and the ground-state densities in density-functional theory (i.e., the minimizers of the Hohenberg--Kohn (HK) variation principle) is studied within the framework of convex conjugation, in a generic setting covering molecular systems, solid-state systems, and more. Having introduced admissible density functionals as functionals that produce the exact ground ground-state energy for a given external potential by minimizing over densities in the HK variation principle, necessary sufficient conditions on such functionals are established to ensure that the RR ground-state densities and the HK ground-state densities are identical. We apply the results to molecular systems in the BO-approximation. For any given potential $v \in L^{3/2}(\mathbb{R}^3) + L^{\infty}(\mathbb{R}^3)$, we establish a one-to-one correspondence between the mixed ground-state densities of the RR variation principle and the mixed ground-state densities of the HK variation principle when the Lieb density-matrix constrained-search universal density functional is taken as the admissible functional. A similar one-to-one correspondence is established between the pure ground-state densities of the RR variation principle and the pure ground-state densities obtained using the HK variation principle with the Levy--Lieb pure-state constrained-search functional. In other words, all physical ground-state densities (pure or mixed) are recovered with these functionals and no false densities (i.e., minimizing densities that are not physical) exist. The importance of topology (i.e., choice of Banach space of densities and potentials) is emphasized and illustrated. The relevance of these results for current-density-functional theory is examined.

We present the self-consistent implementation of current-dependent (hybrid) meta generalized gradient approximation (mGGA) density functionals using London atomic orbitals. A previously proposed generalized kinetic energy density is utilized to implement mGGAs in the framework of Kohn--Sham current density-functional theory (KS-CDFT). A unique feature of the non-perturbative implementation of these functionals is the ability to seamlessly explore a wide range of magnetic fields up to 1 a.u. ($\sim 235000$T) in strength. CDFT functionals based on the TPSS and B98 forms are investigated and their performance is assessed by comparison with accurate CCSD(T) data. In the weak field regime magnetic properties such as magnetizabilities and NMR shielding constants show modest but systematic improvements over GGA functionals. However, in strong field regime the mGGA based forms lead to a significantly improved description of the recently proposed perpendicular paramagnetic bonding mechanism, comparing well with CCSD(T) data. In contrast to functionals based on the vorticity these forms are found to be numerically stable and their accuracy at high field suggests the extension of mGGAs to CDFT via the generalized kinetic energy density should provide a useful starting point for further development of CDFT approximations.

An implementation of coupled-cluster (CC) theory to treat atoms and molecules in finite magnetic fields is presented. The main challenges stem from the magnetic-field dependence in the Hamiltonian, or, more precisely, the appearance of the angular momentum operator, due to which the wave function becomes complex and which introduces a gauge-origin dependence. For this reason, an implementation of a complex CC code is required together with the use of gauge-including atomic orbitals to ensure gauge-origin independence. Results of coupled-cluster singles--doubles--perturbative-triples (CCSD(T)) calculations are presented for atoms and molecules with a focus on the dependence of correlation and binding energies on the magnetic field.

In this paper, an alternative method to range-separated linear-response time-dependent density-functional theory and perturbation theory is proposed to improve the estimation of the energies of a physical system from the energies of a partially interacting system. Starting from the analysis of the Taylor expansion of the energies of the partially interacting system around the physical system, we use an extrapolation scheme to improve the estimation of the energies of the physical system at an intermediate point of the range-separated or linear adiabatic connection where either the electron--electron interaction is scaled or only the long-range part of the Coulomb interaction is included. The extrapolation scheme is first applied to the range-separated energies of the helium and beryllium atoms and of the hydrogen molecule at its equilibrium and stretched geometries. It improves significantly the convergence rate of the energies toward their exact limit with respect to the range-separation parameter. The range-separated extrapolation scheme is compared with a similar approach for the linear adiabatic connection, highlighting the relative strengths and weaknesses of each approach., Comment: to appear in Phys. Rev. A

We present a study of the variation of total energies and excitationenergies along a range-separated adiabatic connection. This connectionlinks the non-interacting Kohn-Sham electronic system to the physicalinteracting system by progressively switching on theelectron-electron interactions whilst simultaneously adjusting aone-electron effective potential so as to keep the ground-statedensity constant. The interactions are introduced in arange-dependent manner, first introducing predominantly long-range,and then all-range, interactions as the physical system is approached,as opposed to the conventional adiabatic connection where theinteractions are introduced by globally scaling the standard Coulomb interaction.Reference data are reported for the He and Be atoms and the H2molecule, obtained by calculating the short-range effective potentialat the full configuration-interaction level using Lieb'sLegendre-transform approach. As the strength of the electron-electroninteractions increases, the excitation energies, calculated for thepartially interacting systems along the adiabatic connection, offerincreasingly accurate approximations to the exact excitation energies.Importantly, the excitation energies calculated at an intermediatepoint of the adiabatic connection are much better approximations tothe exact excitation energies than are the corresponding Kohn-Shamexcitation energies. This is particularly evident in situationsinvolving strong static correlation effects and states with multipleexcitation character, such as the dissociating H2 molecule. Theseresults highlight the utility of long-range interacting referencesystems as a starting point for the calculation of excitation energiesand are of interest for developing and analyzing practical approximaterange-separated density-functional methodologies.

We explore the possibility of calculating electronic excited states by using perturbation theory along a range-separated adiabatic connection. Starting from the energies of a partially interacting Hamiltonian, a first-order correction is defined with two variants of perturbation theory: a straight-forward perturbation theory, and an extension of the G{\"o}rling--Levy one that has the advantage of keeping the ground-state density constant at each order in the perturbation. Only the first, simpler, variant is tested here on the helium and beryllium atoms and on the dihydrogene molecule. The first-order correction within this perturbation theory improves significantly the total ground-and excited-state energies of the different systems. However, the excitation energies are mostly deterio-rated with respect to the zeroth-order ones, which may be explained by the fact that the ionization energy is no longer correct for all interaction strengths. The second variant of the perturbation theory should improve these results but has not been tested yet along the range-separated adiabatic connection.

The universal density functional $F$ of density-functional theory is a complicated and ill-behaved function of the density-in particular, $F$ is not differentiable, making many formal manipulations more complicated. Whilst $F$ has been well characterized in terms of convex analysis as forming a conjugate pair $(E,F)$ with the ground-state energy $E$ via the Hohenberg-Kohn and Lieb variation principles, $F$ is nondifferentiable and subdifferentiable only on a small (but dense) set of its domain. In this article, we apply a tool from convex analysis, Moreau-Yosida regularization, to construct, for any $\epsilon>0$, pairs of conjugate functionals $({}^\epsilon\!E,{}^\epsilon\!F)$ that converge to $(E,F)$ pointwise everywhere as $\epsilon\rightarrow 0^+$, and such that ${}^\epsilon\!F$ is (Fr\'echet) differentiable. For technical reasons, we limit our attention to molecular electronic systems in a finite but large box. It is noteworthy that no information is lost in the Moreau-Yosida regularization: the physical ground-state energy $E(v)$ is exactly recoverable from the regularized ground-state energy ${}^\epsilon\!E(v)$ in a simple way. All concepts and results pertaining to the original $(E,F)$ pair have direct counterparts in results for $({}^\epsilon\! E, {}^\epsilon\!F)$. The Moreau-Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of ${}^\epsilon\!F$, a rigorous formulation of Kohn-Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn-Sham theory.

The $N$-representability problem is the problem of determining whether or not there exists $N$-particle states with some prescribed property. Here we report an affirmative solution to the fermion $N$-representability problem when both the density and paramagnetic current density are prescribed. This problem arises in current-density functional theory and is a generalization of the well-studied corresponding problem (only the density prescribed) in density functional theory. Given any density and paramagnetic current density satisfying a minimal regularity condition (essentially that a von Weiz\"acker-like the canonical kinetic energy density is locally integrable), we prove that there exist a corresponding $N$-particle state. We prove this by constructing an explicit one-particle reduced density matrix in the form of a position-space kernel, i.e.\ a function of two continuous position variables. In order to make minimal assumptions, we also address mathematical subtleties regarding the diagonal of, and how to rigorously extract paramagnetic current densities from, one-particle reduced density matrices in kernel form.

The alternative separation of exchange and correlation energies proposed by Toulouse et al. [Theor. Chem. Acc. 114, 305 (2005)] is explored in the context of multi-configuration range-separated density-functional theory. The new decomposition of the short-range exchange-correlation energy relies on the auxiliary long-range interacting wavefunction rather than the Kohn-Sham (KS) determinant. The advantage, relative to the traditional KS decomposition, is that the wavefunction part of the energy is now computed with the regular (fully-interacting) Hamiltonian. One potential drawback is that, because of double counting, the wavefunction used to compute the energy cannot be obtained by minimizing the energy expression with respect to the wavefunction parameters. The problem is overcome by using short-range optimized effective potentials (OEPs). The resulting combination of OEP techniques with wavefunction theory has been investigated in this work, at the Hartree-Fock (HF) and multi-configuration self-consistent-field (MCSCF) levels. In the HF case, an analytical expression for the energy gradient has been derived and implemented. Calculations have been performed within the short-range local density approximation on H2, N2, Li2 and H2O. Significant improvements in binding energies are obtained with the new decomposition of the short-range energy. The importance of optimizing the short-range OEP at the MCSCF level when static correlation becomes significant has also been demonstrated for H2, using a finite-difference gradient. The implementation of the analytical gradient for MCSCF wavefunctions is currently in progress., Comment: 5 figures

The selection of basic variables in current-density functional theory and formal properties of the resulting formulations are critically examined. Focus is placed on the extent to which the Hohenberg--Kohn theorem, constrained-search approach and Lieb's formulation (in terms of convex and concave conjugation) of standard density-functional theory can be generalized to provide foundations for current-density functional theory. For the well-known case with the gauge-dependent paramagnetic current density as a basic variable, we find that the resulting total energy functional is not concave. It is shown that a simple redefinition of the scalar potential restores concavity and enables the application of convex analysis and convex/concave conjugation. As a result, the solution sets arising in potential-optimization problems can be given a simple characterization. We also review attempts to establish theories with the physical current density as a basic variable. Despite the appealing physical motivation behind this choice of basic variables, we find that the mathematical foundations of the theories proposed to date are unsatisfactory. Moreover, the analogy to standard density-functional theory is substantially weaker as neither the constrained-search approach nor the convex analysis framework carry over to a theory making use of the physical current density.