Arxiv:1705.00238V3 [Physics.Chem-Ph] 19 May 2017
Total Page:16
File Type:pdf, Size:1020Kb
Towards a formal definition of static and dynamic electronic correlations Carlos L. Benavides-Riveros,1, ∗ Nektarios N. Lathiotakis,2 and Miguel A. L. Marques1 1Institut für Physik, Martin-Luther-Universität Halle-Wittenberg, 06120 Halle (Saale), Germany 2Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, GR-11635 Athens, Greece Some of the most spectacular failures of density-functional and Hartree-Fock theories are related to an incorrect description of the so-called static electron correlation. Motivated by recent progress on the N-representability problem of the one-body density matrix for pure states, we propose a way to quantify the static contribution to the electronic correlation. By studying several molecular systems we show that our proposal correlates well with our intuition of static and dynamic electron correlation. Our results bring out the paramount importance of the occupancy of the highest occupied natural spin-orbital in such quantification. PACS numbers: 31.15.V-, 31.15.xr, 31.70.-f I. INTRODUCTION (the orthogonally twisted ethylene C2H4 and the methy- lene CH2, for example) and bond stretching in H2O, C2 The concept of correlation and more precisely the idea and N2, are well described by such a method [6,7]. of correlation energy are central in quantum chemistry. In recent years a considerable effort has been devoted Indeed, the electron-correlation problem (or how the dy- to characterize the correlation of a quantum system in namics of each electron is affected by the others) is terms of more meaningful quantities, such as the Slater perhaps the single largest source of error in quantum- rank for two-electron systems [8,9], the entanglement chemical computations [1]. The success of Hartree-Fock classification for the three-fermion case [10], the squared theory in providing a workable upper bound for the Frobenius norm of the cumulant part of the two-particle ground-state energy is largely due to the fact that a reduced density matrix [11] or the comparison with un- single Slater determinant is usually the simplest wave correlated states [12]. Along with Christian Schilling, we function having the correct symmetry properties for a have recently stressed the importance of the energy gap system of fermions. Since the description of interacting in the understanding of the electronic correlations [13]. fermionic systems requires multi-determinantal reference Notwithstanding, these measures do not draw a distinc- wave functions, the correlation energy is commonly de- tion between qualitatively different kinds of electronic fined as the difference between the exact ground-state correlations. In quantum chemistry, for instance, it is and the Hartree-Fock energy [2,3]. Beyond Hartree-Fock customary to distinguish between static (or nondynamic) theory, numerous other methods (such as configuration and dynamic correlations. The former corresponds to interaction or coupled-cluster theory) aim at reconstruct- configurations which are nearly degenerate with respect ing the part of the energy missing from a description to the reference Slater determinant (if any), whilst the based on a single-determinantal wave function. Indeed, latter arises from the need of mixing the Hartree-Fock one common indicator of the accuracy of a model is, by state with higher-order excited states [14, 15]. Heuris- and large, the percentage of the correlation energy it is tically, one usually states that in systems with (strong) able to recover. static correlation the wavefunction differs qualitatively For small molecules, variational methods based on from the reference Slater determinant, while strong dy- configuration interaction techniques describe well elec- namic correlation implies a wavefunction including a tronic correlations. However, due to their extreme com- large number of excited determinants, all with compa- putational cost, configuration-interaction wavefunctions rable, small occupations. Some of the most spectacular are noteworthily difficult to evaluate for larger systems. failures of the Hartree-Fock theory and density functional Among other procedures at hand, the correlation can theory (with standard exchange-correlation functionals) be treated efficiently by applying a Jastrow correlation are related to an incorrect description of static correla- term to an antisymmetrized wave function (e.g. a sin- tion [16]. arXiv:1705.00238v3 [physics.chem-ph] 19 May 2017 gle Slater determinant or an antisymmetrized geminal It is commonly believed that, to a large extent, both power) within quantum Monte Carlo methods [4,5]. The static and dynamic contributions should be included in so-called Jastrow antisymmetric geminal ansatz accounts the global computation of the electronic correlation. Yet for inter-pair interactions and multiple resonance struc- there are few systems for which one can distinguish un- tures, maintaining a polynomial scaling cost, compara- ambiguously between these two types of correlations. For ble to that of the simpler Jastrow single determinant ap- instance, the ground state of helium has no excited elec- proach. Highly correlated systems, as diradical molecules tronic states nearby, leading therefore to the absence of static correlation. In the dissociation limit of H2 a state with fractional occupations arises [17] and the correla- tion is purely static. According to Hollett and Gill [18], ∗ [email protected] static correlation comes in two “flavors”: one that can be 2 captured by breaking the spin symmetry of the Hartree- vided a set of constraints on the natural occupation num- Fock wave function (like in stretched H2) and another bers, stronger than the Pauli principle [23]. Although that cannot. The measures of correlation proposed so far rudimentary schemes to construct such constraints were, purport to include both static and dynamic correlations, to some extent, routine in quantum-chemistry literature although in an uncontrolled manner [19]. [27], it was not only with the work of Klyachko that this For pure quantum states, global structural features of rich structure could be decrypted. The main goal of this the wave function can be abstracted from local infor- section is to review the physical consequences of such a mation alone. Multiparticle entanglement, for instance, generalization. N can be completely classified with the more accesible one- Given an N-fermion state jΨi 2 ^ [H1], with H1 be- particle picture [20]. Such a characterization is addressed ing the one-particle Hilbert space, the natural occupation by a finite set of linear inequalities satisfied by the eigen- numbers are the eigenvalues fnig and the natural spin- values of the single-particle states [21]. Furthermore, by orbitals are the eigenvectors j'ii of the one-body reduced using the two-particle density matrix and its deviation density matrix, from idempotency, it is possible to propose a criterion X to distinguish static from dynamic correlation, which for ρ^1 ≡ NTrN−1[jΨihΨj] = nij'iih'ij: (1) two-fermion systems only requires the occupancies of the i natural orbitals [22]. Needless to say, grasping global in- The natural occupation numbers, arranged in decreasing formation of a many-body quantum system by tackling order ni ≥ ni+1, fulfill the Pauli condition n1 ≤ 1. The only one-particle information is quite remarkable, mainly natural spin-orbitals define an orthonormal basis B1 for because in this way a linear number of degrees of freedom H1 and can also be used to generate an orthonormal basis is required. N BN for the N-fermion Hilbert space HN ≡ ^ [H1], given Recent progress on the N-representability problem of by the Slater determinants j'i :::'i i = j'i i ^ · · · ^ the one-body reduced density matrix for pure states pro- 1 N 1 j'i i: For practical purposes, the dimension of the one- vides an extension of the well-known Pauli exclusion prin- N particle Hilbert space H1 is usually finite. Henceforth, ciple [23]. This extension is important because it pro- HN;d denotes an antisymmetric N-particle Hilbert space vides stringent constraints beyond those from the Pauli with an underlying d-dimensional one-particle Hilbert principle, which can be used, among others, to improve d space. In principle, the total dimension of HN;d is N , reduced-density-matrix functional theories [24–26]. Our but symmetries usually lower it. main aim in this paper is to employ the generalized Pauli It is by now known that the antisymmetry of N- exclusion principle to establish a general criterion to dis- fermion pure quantum states not only implies the well- tinguish static and dynamic contributions to the elec- known Pauli exclusion principle, which restricts the oc- tronic correlation in fermionic systems. cupation numbers according to [28] 0 ≤ ni ≤ 1, but The paper is organized as follows. For completeness, also entails a set of so-called generalized Pauli constraints Section II summarizes the key aspects of the so-called [23, 29–31]. These take the form of independent linear generalized Pauli exclusion principle and its potential rel- inequalities evance for quantum chemistry. In SectionIII we discuss a Shull-Löwdin-type functional for three-fermion systems, d 0 X i which can be constructed by using the pertinent gen- Dj(~n) ≡ κj + κjni ≥ 0: (2) eralized Pauli constraints along with the spin symme- i=1 tries. Since this functional depends only on the occupa- i Here the coefficients κj 2 Z and j = 1; 2; : : : ; νN;d < 1. tion numbers, it is possible to distinguish the correlation Accordingly, for pure states the spectrum of a physical degree of the so-called Borland-Dennis setting (with an fermionic one-body reduced density matrix must satisfy underlining six-dimensional one-particle Hilbert space) a set of independent linear inequalities of the type (2). by using one-particle information alone. In SectionIV The total number of independent inequalities νN;d de- we discuss a formal way to distinguish static from dy- pends on the number of fermions and the dimension of namic correlation. SectionV is devoted to investigate the underlying one-particle Hilbert space.