Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman’s ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman’s ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman's ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
It is well known that determining the energy of molecules and other quantum many-body systems reduces in the standard approximation to optimizing a simple linear functional of a 12-variable object, the two-electron reduced density matrix (2-RDM). The difficulty is that the variation ensemble for that functional has never been satisfactorily determined. This is known as the N-representability problem of quantum chemistry (which to a large extent is a problem of quantum information theory). The situation has given rise to competing research programs, typically trading more complicated functionals for simpler representability conditions. Chief among them is density functional theory, based on a three-variable object for which the N-representability is trivial, whereas the exact functional is very strange indeed, and probably forever unknowable. An intermediate position is occupied by 1-RDM functional theory. Postulated by Pauli to explain the electronic structure of atoms and molecules, the exclusion principle establishes an upper bound of 1 for the fermionic natural occupation numbers, accordingly allowing no more than one electron in each quantum state. This is a necessary and sufficient condition for a 1-RDM to be the contraction of an ensemble N-body density matrix. The fermionic one-body quantum marginal problem asks whether given natural occupation numbers can arise from an antisymmetric (ensemble or pure) N-particle state. The configuration interaction method affords optimal descriptions of quantum states of atoms and molecules by expanding the wave function in terms of orbital-based configurations of Slater determinants. For these systems, the dimension of the Hilbert space grows binomially with m, the number of spin-orbitals in the basis of the one-particle Hilbert space, and N, the number of electrons of the system. It has been observed that for the rank-six approximation of a pure-state N=3 system, the occupation numbers satisfy some additional constraints, stronger than the Pauli exclusion principle. The recent analysis by Alexander Klyachko and coworkers of the one-body marginal problem of the pure N-fermion state establishes a systematic approach to this type of constraints. In fact, for a pure quantum system of N electrons arranged in m spin-orbitals the occupation numbers satisfy a set of linear inequalities, known as generalized Pauli constraints (GPC). These inequalities define a convex polytope of allowed states in R^m. They are conditions for a 1-RDM to be the contraction of a N-body state. When one of the GPC is completely saturated, the system is said to be pinned, and it lies on one of the facets of the polytope. The nature of those conditions has been explored till now only in a few systems: a model of three spinless fermions confined to a one-dimensional harmonic potential, the lithium isoelectronic series and ground and excited states of some three- and four-electron molecules for the rank being at least twice the number of electrons. For all these systems the inequalities are (quite often) nearly saturated. This is the so-called quasipinning phenomenon. In this PhD thesis we have taken up the challenge of using numerical and analytical methods to examine pinned and quasipinned states, for atoms and molecules, starting from scratch with configuration-interaction and multiconfiguration self-consistent methods. This procedure serves to study the occurrence of quasipinning in realistic systems. A second goal is to show how the subsets of pinned states defined by GPC give rise to the most efficient approach from a computational viewpoint, yielding the leading order of the electron-electron correlations. As a consequence, we underline in this thesis the theoretical and practical importance of Klyachko's approach to the quantum marginal problem and its impact on the competing research programs to determine feasible electronic densities, 1-RDM, 2-RDM, intracular distributions or Wigner density quasiprobabilities. In relation with the above, our research provides a new variational optimization method for few-fermion ground states. We quantitatively confirm its high accuracy for quasipinned systems and derive an upper bound on the error of the correlation energy given by the ratio of the numerical value of the Klyachko inequality and the distance to the Hartree-Fock point. Depending on the details of the algorithm, we are able to reach 98%-99% of the correlation energy for such systems.
We derive necessary conditions for any two-fermion reduced density matrix (2-RDM) to be representable by a pure N-fermion density matrix $\mathrm{\ensuremath{\Psi}}{\mathrm{\ensuremath{\Psi}}}^{*}$ where $\mathrm{\ensuremath{\Psi}}$ is the wave function. These pure N-representability conditions of the 2-RDM are important because they provide stringent constraints beyond those from the Pauli and the generalized Pauli constraints on the structure of many-fermion 2-RDMs and their wave functions. The pure 2-RDM conditions are derived as generalized Pauli constraints on effective one-fermion reduced density matrices (1-RDMs) generated by the removal or addition of a fermion from the wave function. Computationally, we show for four-electron molecules that the derived pure N-representability conditions are nontrivially active for exact ground-state 2-RDMs and that they provide significant restrictions beyond the D, Q, and G ensemble N-representability conditions. Constraints on higher-order p-RDMs where p>2 are derived in a similar fashion. The constraints have potentially significant applications to computing strongly correlated many-fermion states with enhanced accuracy and decreased computational complexity.
The one-particle reduced density-matrix (1-RDM) functional theory is a promising alternative to density-functional theory (DFT) that uses the 1-RDM rather than the electronic density as a basic variable. However, long-standing challenges such as the lack of the Kohn-Sham scheme and the complexity of the pure N-representability conditions are still impeding its wild utilization. Fortunately, ensemble N-representability conditions derived in the natural orbital basis are known and trivial such that almost every functional of the 1-RDM is actually natural orbital functional, which does not perform well for all the correlation regimes. In this work, we propose a variational minimization scheme in the ensemble N-representable domain that is not restricted to the natural orbital representation of the 1-RDM. We show that splitting the minimization into the diagonal and off-diagonal parts of the 1-RDM can open the way toward the development of functionals of the orbital occupations, which remains a challenge for the generalization of site-occupation functional theory in chemistry. Our approach is tested on the uniform Hubbard model using the Müller and the Töws-Pastor functionals, as well as on the dihydrogen molecule using the Müller functional.
The direct variational optimization of the ground-state two-electron reduced-density matrix (2-RDM) is typically performed under ensemble N-representability conditions. Accordingly, variationally obtained 2-RDMs for degenerate ground states may not represent a pure state. When considering only ground-state energetics, the ensemble nature of the 2-RDM is of little consequence. However, the use of ensemble densities within an extended random phase approximation (ERPA) yields astonishingly poor estimates of excitation energies, even for simple atomic systems [H. van Aggelen et al., Comput. Theor. Chem. 1003, 50-54 (2013)]. Here, we outline an approach for the direct variational optimization of ground-state 2-RDMs that satisfy pure-state N-representability known as generalized Pauli constraints. Within the ERPA, 2-RDMs that satisfy both ensemble conditions and the generalized Pauli constraints yield much more reliable estimates of excitation energies than those that satisfy only ensemble conditions.
A crucial theorem in Reduced Density Matrix Functional Theory (RDMFT) suggests that the universal pure and ensemble functionals coincide on their common domain of pure N-representable one-matrices. We refute this by a comprehensive analysis of the geometric picture underlying Levy’s constrained search. Moreover, we then show that the ensemble functional follows instead as the lower convex envelop of the pure functional. It is particularly remarkable that the pure functional determines the ensemble functional even outside its own domain of pure N-representable one-matrices. From a general perspective, this demonstrates that relaxing pure RDMFT to ensemble RDMFT does not necessarily circumvent the complexity of the one-body pure N-representability conditions (generalized Pauli constraints). Instead, the complexity may simply be transferred from the underlying space of pure N-representable one-matrices to the structure of the universal one-matrix functional.