On the measure of electron correlation and entanglement in quantum chemistry based on the cumulant of the second-order reduced density matrix

On the notion of strong correlation in electronic structure theory

Strong correlation has been said to have many faces, and appears to have many synonyms of questionable suitability. In this work we aim not to define the term once and for all, but to highlight one possibility that is both rigorously defined and physically transparent, and remains so in reference to molecules and quantum lattice models. We survey both molecular examples - hydrogen systems (Hn, n = 2, 4, 6), Be2, H-He-H, and benzene - and the half-filled Hubbard model over a range of correlation regimes. Various quantities are examined including the extent of spin symmetry breaking in correlated single-reference wave functions, energetic ratios inspired by the Hubbard model and the Virial theorem, and metrics derived from the one- and two-electron reduced density matrices (RDMs). The trace and the square norm of the cumulant of the two-electron reduced density matrix capture what may well be defined as strong correlation. Accordingly, strong correlation is understood as a statistical dependence between two electrons, and is distinct from the concepts of "correlation energy" and more general than entanglement quantities that require a partitioning of a quantum system into distinguishable subspaces. This work enables us to build a bridge between a rigorous and quantifiable regime of strong electron correlation and more familiar chemical concepts such as anti-aromaticity in the context of Baird's rule.

Geometrical picture of the electron-electron correlation at the large-D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron-electron correlations and might give a bound on the correlation energy. The results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron-electron correlations.

Challenges for variational reduced-density-matrix theory with three-particle N-representability conditions

The direct variational optimization of the two-electron reduced density matrix (2RDM) can provide a reference-independent description of the electronic structure of many-electron systems that naturally capture strong or nondynamic correlation effects. Such variational 2RDM approaches can often provide a highly accurate description of strong electron correlation, provided that the 2RDMs satisfy at least partial three-particle N-representability conditions (e.g., the T2 condition). However, recent benchmark calculations on hydrogen clusters [N. H. Stair and F. A. Evangelista, J. Chem. Phys. 153, 104108 (2020)] suggest that even the T2 condition leads to unacceptably inaccurate results in the case of two- and three-dimensional clusters. We demonstrate that these failures persist under the application of full three-particle N-representability conditions (3POS). A variety of correlation metrics are explored in order to identify regimes under which 3POS calculations become unreliable, and we find that the relative squared magnitudes of the cumulant three- and two-particle reduced density matrices correlate reasonably well with the energy error in these systems. However, calculations on other molecular systems reveal that this metric is not a universal indicator for the reliability of the reduced-density-matrix theory with 3POS conditions.

Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach

This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.

Signature of van der Waals interactions in the cumulant density matrix

Here we propose and implement a universal signature of the van der Waals interactions based on the cumulant part of the two-electron reduced density matrix (2-RDM). Due to the connected property of the cumulant, we can use it to detect the van der Waals interactions between two molecular moieties. In particular, we use the squared Frobenius norm of the cumulant of the 2-RDM, which has been previously shown to provide a size-extensive measure of the electron correlation. As two moieties are separated to infinity, the cumulant Frobenius norm exhibits an r-6 decay to its asymptotic limit, providing a density-based measure of the van der Waals interaction. We study this signature of van der Waals forces in a collection of small molecules of varying geometries. These computations agree with experimental trends of known literature values.

Separation of dynamic and nondynamic correlation

The account of electron correlation and its efficient separation into dynamic and nondynamic parts plays a key role in the development of computational methods. In this paper we suggest a physically-sound matrix formulation to split electron correlation into dynamic and nondynamic parts using the two-particle cumulant matrix and a measure of the deviation from idempotency of the first-order density matrix. These matrices are applied to a two-electron model, giving rise to a simplified electron correlation index that (i) depends only on natural orbitals and their occupancies, (ii) can be straightforwardly decomposed into orbital contributions and (iii) splits into dynamic and nondynamic correlation parts that (iv) admit a local version. These expressions are shown to account for dynamic and nondynamic correlation in a variety of systems containing different electron correlation regimes, thus providing the first separation of dynamic and nondynamic correlation using solely natural orbital occupancies.

Entanglement Measures for Single- and Multireference Correlation Effects

Electron correlation effects are essential for an accurate ab initio description of molecules. A quantitative a priori knowledge of the single- or multireference nature of electronic structures as well as of the dominant contributions to the correlation energy can facilitate the decision regarding the optimum quantum chemical method of choice. We propose concepts from quantum information theory as orbital entanglement measures that allow us to evaluate the single- and multireference character of any molecular structure in a given orbital basis set. By studying these measures we can detect possible artifacts of small active spaces.

Fingerprints of spin-orbital entanglement in transition metal oxides

The concept of spin-orbital entanglement on superexchange bonds in transition metal oxides is introduced and explained on several examples. It is shown that spin-orbital entanglement in superexchange models destabilizes the long-range (spin and orbital) order and may lead either to a disordered spin-liquid state or to novel phases at low temperature which arise from strongly frustrated interactions. Such novel ground states cannot be described within the conventionally used mean field theory which separates spin and orbital degrees of freedom. Even in cases where the ground states are disentangled, spin-orbital entanglement occurs in excited states and may become crucial for a correct description of physical properties at finite temperature. As an important example of this behaviour we present spin-orbital entanglement in the RV O(3) perovskites, with R = La,Pr,…,Y b,Lu, where the finite temperature properties of these compounds can be understood only using entangled states: (i) the thermal evolution of the optical spectral weights, (ii) the dependence of the transition temperatures for the onset of orbital and magnetic order on the ionic radius in the phase diagram of the RV O(3) perovskites, and (iii) the dimerization observed in the magnon spectra for the C-type antiferromagnetic phase of Y V O(3). Finally, it is shown that joint spin-orbital excitations in an ordered phase with coexisting antiferromagnetic and alternating orbital order introduce topological constraints for the hole propagation and will thus radically modify the transport properties in doped Mott insulators where hole motion implies simultaneous spin and orbital excitations.