Simen Kvaal

Ab initio quantum dynamics using coupled-cluster.

The curse of dimensionality (COD) limits the current state-of-the-art ab initio propagation methods for non-relativistic quantum mechanics to relatively few particles. For stationary structure calculations, the coupled-cluster (CC) method overcomes the COD in the sense that the method scales polynomially with the number of particles while still being size-consistent and extensive. We generalize the CC method to the time domain while allowing the single-particle functions to vary in an adaptive fashion as well, thereby creating a highly flexible, polynomially scaling approximation to the time-dependent Schrödinger equation. The method inherits size-consistency and extensivity from the CC method. The method is dubbed orbital-adaptive time-dependent coupled-cluster, and is a hierarchy of approximations to the now standard multi-configurational time-dependent Hartree method for fermions. A numerical experiment is also given.

Choice of basic variables in current-density-functional theory

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 (or 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.

Differentiable but exact formulation of density-functional theory

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. While 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) subset of its domain. In this article, we apply a tool from convex analysis, Moreau-Yosida regularization, to construct, for any ε > 0, pairs of conjugate functionals ((ε)E, (ε)F) that converge to (E, F) pointwise everywhere as ε → 0(+), and such that (ε)F is (Fréchet) 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 (ε)E(v) in a simple way. All concepts and results pertaining to the original (E, F) pair have direct counterparts in results for ((ε)E, (ε)F). The Moreau-Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of (ε)F, a rigorous formulation of Kohn-Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn-Sham theory.

Absorbing boundary conditions for dynamical many-body quantum systems

In numerical studies of the dynamics of unbound quantum mechanical systems, absorbing boundary conditions are frequently applied. Although this certainly provides a useful tool in facilitating the description of the system, its applications to systems consisting of more than one particle are problematic. This is due to the fact that all information about the system is lost upon the absorption of one particle; a formalism based solely on the Schrodinger equation is not able to describe the remainder of the system as particles are lost. Here we demonstrate how the dynamics of a quantum system with a given number of identical fermions may be described in a manner which allows for particle loss. A consistent formalism which incorporates the evolution of sub-systems with a reduced number of particles is constructed through the Lindblad equation. Specifically, the transition from an N-particle system to an (N − 1)-particle system due to a complex absorbing potential is achieved by relating the Lindblad operators to annihilation operators. The method allows for a straight forward interpretation of how many constituent particles have left the system after interaction. We illustrate the formalism using one-dimensional two-particle model problems.

Ab initio computation of the energies of circular quantum dots

We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons.

Fermion

TheN-representability problem is the problem of determining whether there existsN-particle states with some prescribed property. Here we report an affirmative solution to the fermionN-representability problem when both the density and the 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 canonical kinetic energy density is locally integrable), we prove that there exists a correspondingN-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.

N

We give a thorough analysis of the convergence properties of the configuration-interaction method as applied to parabolic quantum dots among other systems, including a priori error estimates. The method converges slowly in general, and in order to overcome this, we propose to use an effective two-body interaction well known from nuclear physics. Through numerical experiments we demonstrate a significant increase in accuracy of the configuration-interaction method.

-representability for prescribed density and paramagnetic current density

The relationship between the densities of ground-state wave functions (i.e., the minimizers of the Rayleigh-Ritz variation principle) and the ground-state densities in density-functional theory (i.e., the minimizers of the Hohenberg-Kohn 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-state energy for a given external potential by minimizing over densities in the Hohenberg-Kohn variation principle, necessary and sufficient conditions on such functionals are established to ensure that the Rayleigh-Ritz ground-state densities and the Hohenberg-Kohn ground-state densities are identical. We apply the results to molecular systems in the Born-Oppenheimer approximation. For any given potential v ∈ L(3/2)(ℝ(3)) + L(∞)(ℝ(3)), we establish a one-to-one correspondence between the mixed ground-state densities of the Rayleigh-Ritz variation principle and the mixed ground-state densities of the Hohenberg-Kohn 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 Rayleigh-Ritz variation principle and the pure ground-state densities obtained using the Hohenberg-Kohn 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.

Harmonic oscillator eigenfunction expansions, quantum dots, and effective interactions

Double eigenvalues are not generic for matrices without any particular structure. A matrix depending linearly on a scalar parameter, A+μB, will, however, generically have double eigenvalues for some values of the parameter μ. In this paper, we consider the problem of finding those values. More precisely, we construct a method to accurately find all scalar pairs (λ,μ) such that A+μB has a double eigenvalue λ, where A and B are given arbitrary complex matrices. The general idea of the globally convergent method is that if μ is close to a solution, then A+μB has two eigenvalues which are close to each other. We fix the relative distance between these two eigenvalues and construct a method to solve and study it by observing that the resulting problem can be stated as a two-parameter eigenvalue problem, which is already studied in the literature. The method, which we call the method of fixed relative distance (MFRD), involves solving a two-parameter eigenvalue problem which returns approximations of all soluti...

Ground-state densities from the Rayleigh–Ritz variation principle and from density-functional theory

It is demonstrated how a numerical approach based on absorbing boundaries may be used to describe the process of non-sequential two-photon double ionization of helium. Contrary to any method based on solving the Schrodinger equation alone, this numerical scheme is able to reconstruct the remaining particles as one particle is absorbed. This may be used to distinguish between single and double ionization. A model of reduced dimensionality, which describes the process at a qualitative level, has been used. The results have been compared with a more conventional method in which the time-dependent Schrodinger equation is solved and the final wavefunction is analysed in terms of projection onto eigenstates of the uncorrelated Hamiltonian, i.e. with no electron–electron interaction included in the final states. It is found that the two methods indeed produce the same total cross sections for the process.