Strictly-correlated Electrons Research Articles

An external dual charge approach to the multimarginal optimal transport with coulomb cost

In this paper, we study the multimarginal optimal transport with Coulomb cost, also known in the physics literature as the Strictly-Correlated Electrons functional. We prove that the dual Kantorovich potential is an electrostatic potential induced by an external charge density, which we call the dual charge . We study its properties and use it to discretize the potential in one and three space dimensions.

Constrained overdamped Langevin dynamics for symmetric multimarginal optimal transportation

The Strictly Correlated Electrons (SCE) limit of the Levy–Lieb functional in Density Functional Theory (DFT) gives rise to a symmetric multi-marginal optimal transport problem with Coulomb cost, where the number of marginal laws is equal to the number of electrons in the system, which can be very large in relevant applications. In this work, we design a numerical method, built upon constrained overdamped Langevin processes to solve Moment Constrained Optimal Transport (MCOT) relaxations (introduced in [A. Alfonsi, R. Coyaud, V. Ehrlacher and D. Lombardi, Approximation of optimal transport problems with marginal moments constraints, Math. Comp. 90 (2021) 689–737; C. Villani, Optimal Transport: Old and New (Springer Science & Business Media, 2008)]) of symmetric multi-marginal optimal transport problems with Coulomb cost. Some minimizers of such relaxations can be written as discrete measures charging a low number of points belonging to a space whose dimension, in the symmetrical case, scales linearly with the number of marginal laws. We leverage the sparsity of those minimizers in the design of the numerical method and prove that there is no strict local minimizer to the resulting problem. We illustrate the performance of the proposed method by numerical examples which solves MCOT relaxations of 3D systems with up to 100 electrons.

Kohn–Sham equations with functionals from the strictly-correlated regime: investigation with a spectral renormalization method

We re-adapt a spectral renormalization method, introduced in nonlinear optics, to solve the Kohn–Sham (KS) equations of density functional theory, with a focus on functionals based on the strictly-correlated electrons (SCE) regime, which are particularly challenging to converge. Important aspects of the method are: (i) the eigenvalues and the density are computed simultaneously; (ii) it converges using randomized initial guesses; (iii) easy to implement. Using this method we could converge for the first time the Kohn–Sham equations with functionals that include the next leading term in the strong-interaction limit of density functional theory, the so called zero-point energy (ZPE) functional as well as with an interaction-strength-interpolation functional that includes both the exact SCE and ZPE terms. This work is the first building block for future studies on quantum systems confined in low dimensions with different statistics and long-range repulsions, such as localization properties of fermions and bosons with strong long-range repulsive interactions in the presence of a random external potential.

Semidefinite Relaxation of Multimarginal Optimal Transport for Strictly Correlated Electrons in Second Quantization

We consider the strictly correlated electron (SCE) limit of the fermionic quantum many-body problem in the second-quantized formalism. This limit gives rise to a multimarginal optimal transport (MM...

Smoothing of Transport Plans with Fixed Marginals and Rigorous Semiclassical Limit of the Hohenberg–Kohn Functional

We prove rigorously that the exact N-electron Hohenberg-Kohn density functional converges in the strongly interacting limit to the strictly correlated electrons (SCE) functional, and that the absolute value squared of the associated constrained-search wavefunction tends weakly in the sense of probability measures to a minimizer of the multi-marginal optimal transport problem with Coulomb cost associated to the SCE functional. This extends our previous work for N=2 [CFK11]. The correct limit problem has been derived in the physics literature by Seidl [Se99] and Seidl, Gori-Giorgi and Savin [SGS07]; in these papers the lack of a rigorous proof was pointed out. We also give a mathematical counterexample to this type of result, by replacing the constraint of given one-body density -- an infinite-dimensional quadratic expression in the wavefunction -- by an infinite-dimensional quadratic expression in the wavefunction and its gradient. Connections with the Lawrentiev phenomenon in the calculus of variations are indicated.

Exchange-Correlation Functionals via Local Interpolation along the Adiabatic Connection.

The construction of density-functional approximations is explored by modeling the adiabatic connection locally, using energy densities defined in terms of the electrostatic potential of the exchange–correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to assess the accuracy of a range of local interpolation models against accurate exchange–correlation energy densities. The importance of the strictly correlated electrons (SCE) functional describing the strong coupling limit is emphasized, enabling the corresponding interpolated functionals to treat strong correlation effects. In addition to exploring the performance of such models numerically for the helium and beryllium isoelectronic series and the dissociation of the hydrogen molecule, an approximate analytic model is presented for the initial slope of the local adiabatic connection. Comparisons are made with approaches based on global models, and prospects for future approximations based on the local adiabatic connection are discussed.

The adiabatic strictly-correlated-electrons functional: kernel and exact properties.

We investigate a number of formal properties of the adiabatic strictly-correlated electrons (SCE) functional, relevant for time-dependent potentials and for kernels in linear response time-dependent density functional theory. Among the former, we focus on the compliance to constraints of exact many-body theories, such as the generalised translational invariance and the zero-force theorem. Within the latter, we derive an analytical expression for the adiabatic SCE Hartree exchange-correlation kernel in one dimensional systems, and we compute it numerically for a variety of model densities. We analyse the non-local features of this kernel, particularly the ones that are relevant in tackling problems where kernels derived from local or semi-local functionals are known to fail.

Hydrogen Molecule Dissociation Curve with Functionals Based on the Strictly Correlated Regime.

Using the dual Kantorovich formulation, we compute the strictly correlated electrons (SCE) functional (corresponding to the exact strong-interaction limit of density functional theory) for the hydrogen molecule along the dissociation curve. We use an exact relation between the Kantorovich potential and the optimal map to compute the comotion function, exploring corrections based on it. In particular, we analyze how the SCE functional transforms in an exact way the electron-electron distance into a one-body quantity, a feature that can be exploited to build new approximate functionals. We also show that the dual Kantorovich formulation provides in a natural way the constant in the Kohn-Sham potential recently introduced by Levy and Zahariev [Phys. Rev. Lett. 2014, 113, 113002] for finite systems.