Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2305.06457 · The Journal of Chemical Physics

N-representability violations in truncated equation-of-motion coupled-cluster methods

Stephen H. Yuwono, A. Eugene DePrince

Published Aug 7, 2023• 3 citations
arXivHTMLPDF

Abstract

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. These observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

DOI
10.1063/5.0157854

Related papers

Generalized Pauli constraints in reduced density matrix functional theory

Iris Theophilou, Nektarios N. Lathiotakis, Miguel A. L. Marques +1 more

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.

•Jan 1, 2015•34 citations

Generalized Pauli constraints in reduced density matrix functional theory

Iris Theophilou, Nektarios N. Lathiotakis, Miguel A. L. Marques +1 more

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.

arXiv:1503.00746•The Journal of Chemical Physics•Apr 16, 2015•52 citations

Comparing real-time coupled-cluster methods through simulation of collective Rabi oscillations

Andreas S. Skeidsvoll, Henrik Koch A

The time-dependent equation-of-motion coupled-cluster singles and doubles (TD-EOM-CCSD) and time-dependent coupled-cluster singles and doubles (TDCCSD) methods are compared by simulating Rabi oscillations for distant atoms in a classical electromagnetic field. While the TD-EOM-CCSD simulations are numerically stable, the oscillating time-dependent energy expectation value displays an incorrect scaling with the number of atoms resonant with the field. On the other hand, the TDCCSD simulations exhibit the correct scaling in the initial stages of the Rabi cycle but break down when the multiatom system approaches complete population inversion. The numerical stability of the TDCCSD simulations is shown to be very sensitive to any interactions between the atoms and to left amplitudes that do not participate in the expression for time-dependent expectation values for noninteracting atoms. Moreover, we introduce a general theoretical framework for describing the two methods, where the cluster amplitude time derivatives serve as auxiliary conditions related to a shift of the time-dependent Hamiltonian matrix. In this framework, time-dependent coupled-cluster methods exhibit a shifted Hamiltonian matrix with a block upper triangular structure in terms of the number of excited noninteracting subsystems, explaining the correct scaling properties of the TDCCSD method.

arXiv:2301.05546•Physical Review A•Sep 28, 2023•6 citations

Exploring the accuracy of the equation-of-motion coupled-cluster band gap of solids

Evgeny Moerman, Henrique Pereira Coutada Miranda, Alejandro Gallo +7 more

While the periodic equation-of-motion coupled-cluster (EOM-CC) method promises systematic improvement of electronic band gap calculations in solids, its practical application at the singles and doubles level (EOM-CCSD) is hindered by severe finite-size errors in feasible simulation cells. We present a hybrid approach combining EOM-CCSD with the computationally less demanding G W approximation to estimate thermodynamic limit band gaps for several insulators and semiconductors. Our method substantially reduces required cell sizes while maintaining accuracy. Comparisons with experimental gaps and self-consistent G W calculations reveal that deviations in EOM-CCSD predictions correlate with reduced single excitation character of the excited many-electron states. Our work not only provides a computationally tractable approach to EOM-CC calculations in solids but also reveals fundamental insights into the role of single excitations in electronic-structure theory.

arXiv:2501.18667•Physical review. B./Physical review. B•Mar 26, 2025•9 citations

On the relationship between parametric two-electron reduced-density-matrix methods and the coupled electron pair approximation

A. Eugene DePrince, David A. Mazziotti

Parametric two-electron reduced-density-matrix (p-2RDM) methods have enjoyed much success in recent years; the methods have been shown to exhibit accuracies greater than coupled cluster with single and double substitutions (CCSD) for both closed- and open-shell ground-state energies, properties, geometric parameters, and harmonic frequencies. The class of methods is herein discussed within the context of the coupled electron pair approximation (CEPA), and several CEPA-like topological factors are presented for use within the p-2RDM framework. The resulting p-2RDM/n methods can be viewed as a density-based generalization of CEPA/n family that are numerically very similar to traditional CEPA methodologies. We cite the important distinction that the obtained energies represent stationary points, facilitating the efficient evaluation of properties and geometric derivatives. The p-2RDM/n formalism is generalized for an equal treatment of exclusion-principle-violating (EPV) diagrams that occur in the occupied and virtual spaces. One of these general topological factors is shown to be identical to that proposed by Kollmar [C. Kollmar, J. Chem. Phys. 125, 084108 (2006)], derived in an effort to approximately enforce the D, Q, and G conditions for N-representability in his size-extensive density matrix functional.

arXiv:1009.0047•Aug 31, 2010•0 citations

Pure- N -representability conditions of two-fermion reduced density matrices

David A. Mazziotti

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.

Physical Review A•Sep 29, 2016•67 citations
One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. These observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.
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.
The time-dependent equation-of-motion coupled-cluster singles and doubles (TD-EOM-CCSD) and time-dependent coupled-cluster singles and doubles (TDCCSD) methods are compared by simulating Rabi oscillations for distant atoms in a classical electromagnetic field. While the TD-EOM-CCSD simulations are numerically stable, the oscillating time-dependent energy expectation value displays an incorrect scaling with the number of atoms resonant with the field. On the other hand, the TDCCSD simulations exhibit the correct scaling in the initial stages of the Rabi cycle but break down when the multiatom system approaches complete population inversion. The numerical stability of the TDCCSD simulations is shown to be very sensitive to any interactions between the atoms and to left amplitudes that do not participate in the expression for time-dependent expectation values for noninteracting atoms. Moreover, we introduce a general theoretical framework for describing the two methods, where the cluster amplitude time derivatives serve as auxiliary conditions related to a shift of the time-dependent Hamiltonian matrix. In this framework, time-dependent coupled-cluster methods exhibit a shifted Hamiltonian matrix with a block upper triangular structure in terms of the number of excited noninteracting subsystems, explaining the correct scaling properties of the TDCCSD method.
While the periodic equation-of-motion coupled-cluster (EOM-CC) method promises systematic improvement of electronic band gap calculations in solids, its practical application at the singles and doubles level (EOM-CCSD) is hindered by severe finite-size errors in feasible simulation cells. We present a hybrid approach combining EOM-CCSD with the computationally less demanding G W approximation to estimate thermodynamic limit band gaps for several insulators and semiconductors. Our method substantially reduces required cell sizes while maintaining accuracy. Comparisons with experimental gaps and self-consistent G W calculations reveal that deviations in EOM-CCSD predictions correlate with reduced single excitation character of the excited many-electron states. Our work not only provides a computationally tractable approach to EOM-CC calculations in solids but also reveals fundamental insights into the role of single excitations in electronic-structure theory.
Parametric two-electron reduced-density-matrix (p-2RDM) methods have enjoyed much success in recent years; the methods have been shown to exhibit accuracies greater than coupled cluster with single and double substitutions (CCSD) for both closed- and open-shell ground-state energies, properties, geometric parameters, and harmonic frequencies. The class of methods is herein discussed within the context of the coupled electron pair approximation (CEPA), and several CEPA-like topological factors are presented for use within the p-2RDM framework. The resulting p-2RDM/n methods can be viewed as a density-based generalization of CEPA/n family that are numerically very similar to traditional CEPA methodologies. We cite the important distinction that the obtained energies represent stationary points, facilitating the efficient evaluation of properties and geometric derivatives. The p-2RDM/n formalism is generalized for an equal treatment of exclusion-principle-violating (EPV) diagrams that occur in the occupied and virtual spaces. One of these general topological factors is shown to be identical to that proposed by Kollmar [C. Kollmar, J. Chem. Phys. 125, 084108 (2006)], derived in an effort to approximately enforce the D, Q, and G conditions for N-representability in his size-extensive density matrix functional.
We derive necessary conditions for any two-fermion reduced density matrix (2-RDM) to be representable by a pure NNN-fermion density matrix $\mathrm{\ensuremath{\Psi}}{\mathrm{\ensuremath{\Psi}}}^{*}$ where $\mathrm{\ensuremath{\Psi}}$ is the wave function. These pure NNN-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 NNN-representability conditions are nontrivially active for exact ground-state 2-RDMs and that they provide significant restrictions beyond the D, Q, and G ensemble NNN-representability conditions. Constraints on higher-order ppp-RDMs where p>2p>2p>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.