The Pauli exclusion principle is fundamental to understanding electronic quantum systems, imposing constraints on the expected occupancies n(i) of orbitals phi i, such that 0 <= n(i) <= 2. In this work, we refine the underlying one-body N-representability problem by incorporating spin symmetries and a potential degree of mixedness omega of the N-electron quantum state. Employing basic tools from representation theory, convex analysis, and discrete geometry, we derive a comprehensive solution to this problem. Specifically, we demonstrate that the set of admissible orbital one-body reduced density matrices is fully characterized by linear spectral constraints on the natural orbital occupation numbers, defining a convex polytope Sigma(N,S)(w) C [0, 2](d). These constraints are independent of the magnetization M and the number d of orbitals, while their dependence on N and the total spin S is linear, and we thus calculate them for arbitrary system sizes and spin quantum numbers. Our results provide a crucial missing cornerstone for ensemble density (matrix) functional theory.

The N-representability problem places fundamental constraints on reduced density matrices (RDMs) that originate from physical many-fermion quantum states. Motivated by recent developments in functional theories, we introduce a hierarchy of ensemble one-body N-representability problems that incorporate partial knowledge of the one-body RDMs (1RDMs) within an ensemble of N-fermion states with fixed weights wi. Specifically, we propose a systematic relaxation that reduces the refined problem-where full 1RDMs are fixed for certain ensemble elements-to a more tractable form involving only natural occupation number vectors. Remarkably, we show that this relaxed problem is related to a generalization of Horn's problem, enabling an explicit solution by combining its constraints with those of the weighted ensemble N-representability conditions. An additional convex relaxation yields a convex polytope that provides physically meaningful restrictions on lattice site occupations in ensemble density functional theory for excited states.

We develop a systematic framework for the spin adaptation of the cumulants of p-particle reduced density matrices (RDMs), with explicit constructions for p = 1 to 3. These spin-adapted cumulants enable rigorous treatment of both S (z) and S (2) symmetries in quantum systems, providing a foundation for spin-resolved electronic structure methods. We show that complete spin adaptation-referred to as complete S-representability-can be enforced by constraining the variances of S (z) and S (2), which require the 2-RDM and 4-RDM, respectively. Importantly, the cumulants of RDMs scale linearly with system size-size-extensive-making them a natural object for incorporating spin symmetries in scalable electronic structure theories. The developed formalism is applicable to density-based methods, one-particle RDM functional theories, and two-particle RDM methods. We further extend the approach to spin-orbit-coupled systems via total angular momentum adaptation. Beyond spin, the framework enables the adaptation of RDM theories to additional symmetries through the construction of suitable irreducible tensor operators.

We establish a toolbox for studying and applying spin-adapted generalized Pauli constraints (GPCs) in few-electron quantum systems. By exploiting the spin symmetry of realistic N-electron wave functions, the underlying one-body pure N-representability problem simplifies, allowing us to calculate the GPCs for larger system sizes than previously accessible. We then uncover and rigorously prove a superselection rule that highlights the significance of GPCs: whenever a spin-adapted GPC is (approximately) saturated-referred to as (quasi)pinning-the corresponding N-electron wave function assumes a simplified structure. Specifically, in a configuration interaction expansion based on natural orbitals only very specific spin configuration state functions may contribute. To assess the nontriviality of (quasi)pinning, we introduce a geometric measure that contrasts it with the (quasi)pinning induced by simple (spin-adapted) Pauli constraints. Applications to few-electron systems suggest that previously observed quasipinning largely stems from spin symmetries.

The extension of the Rayleigh-Ritz variational principle to ensemble states rho(w)equivalent to & sum,.(k)w(k)|Psi(k)><Psi(k)| with fixed weights w(k) lies ultimately at the heart of several recent methodological developments for targeting excitation energies by variational means. Prominent examples are density and density matrix functional theory, Monte Carlo sampling, state-average complete active space self-consistent field methods and variational quantum eigensolvers. In order to provide a sound basis for all these methods and to improve their current implementations, we prove the validity of the underlying critical hypothesis: Whenever the ensemble energy is well-converged, the same holds true for the ensemble state rho(w) as well as the individual eigenstates |Psi(k)> and eigenenergies E-k. To be more specific, we derive linear bounds d(-)Delta E-w <= Delta Q <= d(+)Delta E-w on the errors Delta Q of these sought-after quantities. A subsequent analytical analysis and numerical illustration proves the tightness of our universal inequalities. Our results and particularly the explicit form of d(+/-)equivalent to d(+/-)((Q))(w,E) provide valuable insights into the optimal choice of the auxiliary weights w(k) in practical applications.

It is the ultimate goal of this work to foster synergy between quantum chemistry and the flourishing field of quantum information theory. For this, we first translate quantum information concepts, such as entanglement and correlation, into the context of quantum chemical systems. In particular, we establish two conceptually distinct perspectives on 'electron correlation', leading to a notion of orbital and particle correlation. We then demonstrate that particle correlation equals total orbital correlation minimized over all orbital bases. Accordingly, particle correlation resembles the minimal, thus intrinsic, complexity of many-electron wave functions, while orbital correlation quantifies their complexity relative to a basis. We illustrate these concepts of intrinsic and extrinsic correlation complexity in molecular systems, which also manifests the crucial link between the two correlation pictures. Our results provide theoretical justification for the long-favored natural orbitals for simplifying electronic structures, and open new pathways for developing more efficient approaches towards the electron correlation problem. QIT offers a comprehensive toolbox for electron correlation analysis, and development of new methods for solving the electronic problem. QChem in turn provides a platform to realize quantum technology, and supplies the valuable resource of quantum entanglement in molecules.

The simultaneous treatment of static and dynamic correlations in strongly correlated electron systems is a critical challenge. In particular, finding a universal scheme for identifying a single-particle orbital basis that minimizes the representational complexity of the many-body wave function is a formidable and longstanding problem. As a contribution toward its solution, we show that the total orbital correlation actually reveals and quantifies the intrinsic complexity of the wave function, once it is minimized via orbital rotations. To demonstrate the power of this concept in practice, an iterative scheme is proposed to optimize the orbitals by minimizing the total orbital correlation calculated by the tailored coupled cluster singles and doubles (TCCSD) ansatz. The optimized orbitals enable the limited TCCSD ansatz to capture more nontrivial information on the many-body wave function, indicated by the improved wave function and energy. An initial application of this scheme shows great improvement of TCCSD in predicting the singlet ground state potential energy curves of the strongly correlated C-2 and Cr-2 molecule.

The goal of the present work is to guide the development of quantum technologies in the context of fermionic systems. For this, we first elucidate the process of entanglement swapping in electron systems such as atoms, molecules or solid bodies. This demonstrates the significance of the number-parity superselection rule and highlights the relevance of localized few-orbital subsystems for quantum information processing tasks. Then, we explore and quantify the entanglement between localized orbitals in two systems, a tight-binding model of non-interacting electrons and the hydrogen ring. For this, we apply the first closed formula of a faithful entanglement measure, derived in (arXiv:2207.03377) as an extension of the von Neumann entropy to genuinely correlated many-orbital systems. For both systems, long-distance entanglement is found at low and high densities eta, whereas for medium densities, eta approximate to 12 , practically only neighboring orbitals are entangled. The Coulomb interaction does not change the entanglement pattern qualitatively except for low and high densities where the entanglement increases as function of the distance between both orbitals.

We propose an effective quantum information-assisted complete active space optimization scheme (QICAS). What sets QICAS apart from other correlation-based selection schemes is (i) the use of unique measures from quantum information that assess the correlation in electronic structures in an unambiguous and predictive manner and (ii) an orbital optimization step that minimizes the correlation discarded by the active space approximation. Equipped with these features, QICAS yields, for smaller correlated molecule, sets of optimized orbitals with respect to which the complete active space configuration interaction energy reaches the corresponding complete active space self-consistent field (CASSCF) energy within chemical accuracy. For more challenging systems such as the chromium dimer, QICAS offers an excellent starting point for CASSCF by greatly reducing the number of iterations required for numerical convergence. Accordingly, our study validates a profound empirical conjecture: the energetically optimal nonactive spaces are predominantly those that contain the least entanglement.

From a geometric point of view, Pauli's exclusion principle de -fines a hypersimplex. This convex polytope describes the compatibility of 1-fermion and N-fermion density matrices,. therefore, it coincides with the convex hull of the pure N-representable 1-fermion density matrices. Consequently, the description of ground state physics through 1-fermion density matrices may not necessitate the intricate pure state generalized Pauli constraints. In this article, we study the generalization of the 1-body N-representability problem to ensemble states with fixed spectrum w, in order to describe finite-temperature states and distinctive mixtures of ex-cited states. By employing ideas from convex analysis and combinatorics, we present a comprehensive solution to the corresponding convex relaxation, thus circumventing the complexity of generalized Pauli constraints. In particular, we adapt and further develop tools such as symmetric poly-topes, sweep polytopes, and Gale order. For both fermions and bosons, generalized exclusion principles are discovered, which we determine for any number of particles and dimension of the 1-particle Hilbert space. These exclusion principles are expressed as linear inequalities satisfying hierarchies determined by the nonzero entries of w. The two families of polytopes resulting from these inequalities are part of the new class of so-called lineup polytopes.

To advance the foundation of one-particle reduced density matrix functional theory (1RDMFT), we refine and relate some of its fundamental features and underlying concepts. We define by concise means the scope of a 1RDMFT, identify its possible natural variables, and explain how symmetries could be exploited. In particular, for systems with time-reversal symmetry, we explain why there exist six equivalent universal functionals, prove concise relations among them, and conclude that the important notion of v-representability is relative to the scope and choice of variable. All these fundamental concepts are then comprehensively discussed and illustrated for the Hubbard dimer and its generalization to arbitrary pair interactions W. For this, we derive by analytical means the pure and ensemble functionals with respect to both the real- and complex-valued Hilbert space. The comparison of various functionals allows us to solve the underlying v-representability problems analytically, and the dependence of its solution on the pair interaction is demonstrated. Intriguingly, the gradient of each universal functional is found to always diverge repulsively on the boundary of the domain. In that sense, this key finding emphasizes the universal character of the fermionic exchange force, recently discovered and proven in the context of translationally invariant one-band lattice models.

We initiate the recently proposed w-ensemble one-particle reduced density matrix func-tional theory (w-RDMFT) by deriving the first functional approximations and illustrate how excitation energies can be calculated in practice. For this endeavour, we first study the symmetric Hubbard dimer, constituting the building block of the Hubbard model, for which we execute the Levy-Lieb constrained search. Second, due to the particular suit-ability of w-RDMFT for describing Bose-Einstein condensates, we demonstrate three con-ceptually different approaches for deriving the universal functional in a homogeneous Bose gas for arbitrary pair interaction in the Bogoliubov regime. Remarkably, in both systems the gradient of the functional is found to diverge repulsively at the boundary of the functional's domain, extending the recently discovered Bose-Einstein condensation force to excited states. Our findings highlight the physical relevance of the generalized exclusion principle for fermionic and bosonic mixed states and the curse of universality in functional theories.

"The leading terms in the large-R asymptotics of the functional of the one-electron reduced density matrix for the ground-state energy of the H-2 molecule with the internuclear separation R are derived thanks to the solution of the phase dilemma at the R -> infinity limit. At this limit, the respective natural orbitals (NOs) are given by symmetric and antisymmetric combinations of ""half-space "" orbitals with the corresponding natural amplitudes having the same amplitudes but opposite signs. Minimization of the resulting explicit functional yields the large-R asymptotics for the occupation numbers of the weakly occupied NOs and the C-6 dispersion coefficient. The highly accurate approximates for the radial components of the p-type ""half-space "" orbitals and the corresponding occupation numbers (that decay like R-6), which are available for the first time thanks to the development of the present formalism, have some unexpected properties."

From a geometric point of view, Pauli's exclusion principle de -fines a hypersimplex. This convex polytope describes the compatibility of 1-fermion and N-fermion density matrices,. therefore, it coincides with the convex hull of the pure N-representable 1-fermion density matrices. Consequently, the description of ground state physics through 1-fermion density matrices may not necessitate the intricate pure state generalized Pauli constraints. In this article, we study the generalization of the 1-body N-representability problem to ensemble states with fixed spectrum w, in order to describe finite-temperature states and distinctive mixtures of ex-cited states. By employing ideas from convex analysis and combinatorics, we present a comprehensive solution to the corresponding convex relaxation, thus circumventing the complexity of generalized Pauli constraints. In particular, we adapt and further develop tools such as symmetric poly-topes, sweep polytopes, and Gale order. For both fermions and bosons, generalized exclusion principles are discovered, which we determine for any number of particles and dimension of the 1-particle Hilbert space. These exclusion principles are expressed as linear inequalities satisfying hierarchies determined by the nonzero entries of w. The two families of polytopes resulting from these inequalities are part of the new class of so-called lineup polytopes.

The second quantum revolution is all about exploiting the quantum nature of atoms and molecules to execute quantum information processing tasks. To boost this growing endeavor and by anticipating the key role of quantum chemistry therein, our work establishes a framework for systematically exploring, quantifying and dissecting correlation effects in molecules. By utilizing the geometric picture of quantum states we compare-on a unified basis and in an operationally meaningful way-total, quantum and classical correlation and entanglement in molecular ground states. To unlock and maximize the quantum informational resourcefulness of molecules an orbital optimization scheme is developed, leading to a paradigm-shifting insight: a single covalent bond equates to the entanglement 2ln(2)

Motivated by the Penrose-Onsager criterion for Bose-Einstein condensation we propose a functional theory for targeting low-lying excitation energies of bosonic quantum systems through the one-particle picture. For this, we employ an extension of the Rayleigh-Ritz variational principle to ensemble states with spectrum w and prove a corresponding generalization of the Hohenberg-Kohn theorem: the underlying one-particle reduced density matrix determines all properties of systems of N identical particles in their w -ensemble states. Then, to circumvent the v-representability problem common to functional theories, and to deal with energetic degeneracies, we resort to the Levy-Lieb constrained search formalism in combination with an exact convex relaxation. The corresponding bosonic one-body w -ensemble N-representability problem is solved comprehensively. Remarkably, this reveals a complete hierarchy of bosonic exclusion principle constraints in conceptual analogy to Pauli's exclusion principle for fermions and recently discovered generalizations thereof.

State-average calculations based on a mixture of states are increasingly being exploited across chemistry and physics as versatile procedures for addressing excitations of quantum many-body systems. If not too many states should need to be addressed, calculations performed on individual states are also a common option. Here we show how the two approaches can be merged into one method, dealing with a generalized yet single pure state. Implications in electronic structure calculations are discussed and for quantum computations are pointed out.

In Phys. Rev. Lett. 2021, 127, 023001 a reduced density matrix functional theory (RDMFT) was proposed for calculating energies of selected eigenstates of interacting many-Fermion systems. Here, we develop a solid foundation for this so-called w-RDMFT and present the details of various derivations. First, we explain how a generalization of the Ritz variational principle to ensemble states with fixed weights w in combination with the constrained search would lead to a universal functional of the one-particle reduced density matrix. To turn this into a viable functional theory, however, we also need to implement an exact convex relaxation. This general procedure includes Valone's pioneering work on ground state RDMFT as the special case w = (1,0, ...). Then, we work out in a comprehensive manner a methodology for deriving a compact description of the functional's domain. This leads to a hierarchy of generalized exclusion principle constraints which we illustrate in great detail. By anticipating their future pivotal role in functional theories and to keep our work self-contained, several required concepts from convex analysis are introduced and discussed.

The operator algebra of fermionic modes is isomorphic to that of qubits, the difference between them is twofold: the embedding of subalgebras corresponding to mode subsets and multiqubit subsystems on the one hand, and the parity superselection in the fermionic case on the other. We discuss these two fundamental differences extensively, and illustrate these through the Jordan-Wigner representation in a coherent, self-contained, pedagogical way, from the point of view of quantum information theory. Our perspective leads us to develop useful new tools for the treatment of fermionic systems, such as the fermionic (quasi-)tensor product, fermionic canonical embedding, fermionic partial trace, fermionic products of maps and fermionic embeddings of maps. We formulate these by direct, easily applicable formulas, without mode permutations, for arbitrary partitionings of the modes. It is also shown that fermionic reduced states can be calculated by the fermionic partial trace, containing the proper phase factors. We also consider variants of the notions of fermionic mode correlation and entanglement, which can be endowed with the usual, local operation based motivation, if the parity superselection rule is imposed. We also elucidate some other fundamental points, related to joint map extensions, which make the parity superselection inevitable in the description of fermionic systems.

We propose and work out a reduced density matrix functional theory (RDMFT) for calculating energies of eigenstates of interacting many-electron systems beyond the ground state. Various obstacles which historically have doomed such an approach to be unfeasible are overcome. First, we resort to a generalization of the Ritz variational principle to ensemble states with fixed weights. This in combination with the constrained search formalism allows us to establish a universal functional of the one-particle reduced density matrix. Second, we employ tools from convex analysis to circumvent the too involved N-representability constraints. Remarkably, this identifies Valone's pioneering work on RDMFT as a special case of convex relaxation and reveals that crucial information about the excitation structure is contained in the functional's domain. Third, to determine the crucial latter object, a methodology is developed which eventually leads to a generalized exclusion principle. The corresponding linear constraints are calculated for systems of arbitrary size.

One-particle reduced density matrix functional theory would potentially be the ideal approach for describing Bose-Einstein condensates. It namely replaces the macroscopically complex wave function by the simple one-particle reduced density matrix, and therefore provides direct access to the degree of condensation and still recovers quantum correlations in an exact manner. We initiate and establish this theory by deriving the respective universal functional F for homogeneous Bose-Einstein condensates with arbitrary pair interaction. Most importantly, the successful derivation necessitates a particle-number conserving modification of Bogoliubov theory and a solution of the common phase dilemma of functional theories. We then illustrate this approach in several bosonic systems such as homogeneous Bose gases and the Bose-Hubbard model. Remarkably, the general form of F reveals the existence of a universal Bose-Einstein condensation force which provides an alternative and more fundamental explanation for quantum depletion.

The realization of Bose-Einstein condensation in ultracold trapped gases has led to a revival of interest in this fascinating quantum phenomenon. This experimental achievement necessitated both extremely low temperatures and sufficiently weak interactions. Particularly in reduced spatial dimensionality even an infinitesimal interaction immediately leads to a departure to quasi-condensation. We propose a system of strongly interacting bosons, which overcomes those obstacles by exhibiting a number of intriguing related features: (i) The tuning of just a single control parameter drives a transition from quasi-condensation to complete condensation, (ii) the destructive influence of strong interactions is compensated by the respective increased mobility, (iii) topology plays a crucial role since a crossover from one- to ‘infinite’-dimensionality is simulated, (iv) a ground state gap opens, which makes the condensation robust to thermal noise. Remarkably, all these features can be derived by analytical and exact numerical means despite the non-perturbative character of the system.

A recent development in quantum chemistry has established the quantum mutual information between orbitals as a major descriptor of electronic structure. This has already facilitated remarkable improvements in numerical methods and may lead to a more comprehensive foundation for chemical bonding theory. Building on this promising development, our work provides a refined discussion of quantum information theoretical concepts by introducing the physical correlation and its separation into classical and quantum parts as distinctive quantifiers of electronic structure. In particular, we succeed in quantifying the entanglement. Intriguingly, our results for different molecules reveal that the total correlation between orbitals is mainly classical, raising questions about the general significance of entanglement in chemical bonding. Our work also shows that implementing the fundamental particle number superselection rule, so far not accounted for in quantum chemistry, removes a major part of correlation and entanglement seen previously. In that respect, realizing quantum information processing tasks with molecular systems might be more challenging than anticipated.