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

The Journal of Chemical Physics

Extracting dominant pair correlations from many-body wave functions

Gregory J. O. Beran, Martin Head‐Gordon

Published Jun 16, 2004• 12 citations
Publisher

Abstract

The singular value decomposition of the n-particle excitation operator as determined by coupled cluster or perturbation theory is used to extract the dominant and interesting electron-electron correlations from complex molecular wave functions. As an example of the very general formalism, the decomposition of the T(2) operator obtained from coupled cluster doubles calculations is used to analyze the strength and character of pair correlations in a variety of molecules with interesting electronic structure. The magnitude of the largest singular value(s) determines the strength of the correlation(s), and the corresponding right- and left-hand singular vectors characterize the physical and spatial nature of the correlations. The primary advantage of this tool over natural orbital analysis is that it provides direct associations between the occupied and virtual geminals involved in the correlations.

DOI
10.1063/1.1756860

Related papers

Tensor decomposition methods for correlated electron pairs

Airi Kawasaki, Osamu Sugino

We analyze wave functions constructed as a sum of product of two-electron functions, or as a polynomial of geminals, to investigate their ability to represent the ground state of a strongly correlated few-body system. The known difficulty associated with variational determination of the total energy is overcome by applying a tensor decomposition method called Waring decomposition. Convergence speed of the total energy is compared for various polynomial types. The result provides information bridging between geminal product wave functions and the full-CI in the strongly correlated regime, thereby enriching knowledge on the hierarchy of molecular orbital theories of electron pairs.

arXiv:1805.06138•May 16, 2018•1 citation

Singular value decomposition approach for the approximate coupled-cluster method

Tomoko Kinoshita, Osamu Hino, Rodney J. Bartlett

We present a method for the approximate solution of the coupled-cluster (CC) equation, based upon the singular value decomposition (SVD). The key idea of this method is that we use SVD for the cluster amplitudes to exploit the physically important states in the CC equation. This method enables us to significantly reduce the computational requirements for CC calculations without losing the essence of the method. Relationships to the density matrix renormalization group theory and the local correlation methods are mentioned. We perform pilot calculations on some atoms and molecules to investigate the applicability of the method.

The Journal of Chemical Physics•Oct 15, 2003•87 citations

Analysis of two-orbital correlations in wave functions restricted to electron-pair states

Katharina Bogusławski, Paweł Tecmer, Örs Legeza

Wave functions constructed from electron-pair states can accurately model strong electron correlation effects and are promising approaches especially for larger many-body systems. In this article, we analyze the nature and the type of electron correlation effects that can be captured by wave functions restricted to electron-pair states. We focus on the pair-coupled-cluster doubles (pCCD) ansatz also called the antisymmetric product of the 1-reference orbital geminal (AP1roG) method, combined with an orbital optimization protocol presented in Boguslawski et al. [Phys. Rev. B 89, 201106(R) (2014)], whose performance is assessed against electronic structures obtained form density-matrix renormalization-group reference data. Our numerical analysis covers model systems for strong correlation: the one-dimensional Hubbard model with a periodic boundary condition as well as metallic and molecular hydrogen rings. Specifically, the accuracy of pCCD/AP1roG is benchmarked using the single-orbital entropy, the orbital-pair mutual information, as well as the eigenvalue spectrum of the one-orbital and two-orbital reduced density matrices. Our study indicates that contributions from singly occupied states become important in the strong correlation regime which highlights the limitations of the pCCD/AP1roG method. Furthermore, we examine the effect of orbital rotations within the pCCD/AP1roG model on correlations between orbital pairs.

arXiv:1606.08503•Physical review. B./Physical review. B•Oct 14, 2016•55 citations

Coupled clusters tailored by matrix product state wave functions

Andrej Antalík

The central problem in the modern electronic structure theory is the calculation of cor- relation energy, possibly by an approach that would account for both static and dynamic correlation in an efficient, balanced and accurate way. In this thesis, I present a collection of methods that combine the effective treatment of dynamic correlation by the coupled cluster theory with density matrix renormalization group, a well-established technique for calculations of strongly correlated systems. The connection between them is achieved via the tailored coupled clusters (TCC) ansatz, which conveniently does not impose any ad- ditional computational costs. After the successful initial assessment, we developed more efficient implementations of these methods by employing the local approaches based on pair natural orbitals. This way, we extended the range of possible applications to larger systems with thousands of basis functions. To assess the accuracy of TCC as well as its local counterparts, we performed a variety of benchmark calculations ranging from small, yet challenging systems such as the nitrogen molecule or tetramethyleneethane diradical, to larger molecules like oxo-Mn(Salen) or Fe(II)-porphyrin model. 1

•Mar 30, 2021•0 citations

Coupled cluster-inspired geminal wavefunctions

Pratiksha B. Gaikwad, Taewon David Kim, Michelle Richer +5 more

Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.

arXiv:2310.01764•The Journal of Chemical Physics•Apr 10, 2024•20 citations

An Alternative Algebraic Framework for the Simplification of Coupled Cluster Type Expectation Values

Daniel Pape

The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.

•Nov 12, 2013•0 citations
The singular value decomposition of the n-particle excitation operator as determined by coupled cluster or perturbation theory is used to extract the dominant and interesting electron-electron correlations from complex molecular wave functions. As an example of the very general formalism, the decomposition of the T(2) operator obtained from coupled cluster doubles calculations is used to analyze the strength and character of pair correlations in a variety of molecules with interesting electronic structure. The magnitude of the largest singular value(s) determines the strength of the correlation(s), and the corresponding right- and left-hand singular vectors characterize the physical and spatial nature of the correlations. The primary advantage of this tool over natural orbital analysis is that it provides direct associations between the occupied and virtual geminals involved in the correlations.
We analyze wave functions constructed as a sum of product of two-electron functions, or as a polynomial of geminals, to investigate their ability to represent the ground state of a strongly correlated few-body system. The known difficulty associated with variational determination of the total energy is overcome by applying a tensor decomposition method called Waring decomposition. Convergence speed of the total energy is compared for various polynomial types. The result provides information bridging between geminal product wave functions and the full-CI in the strongly correlated regime, thereby enriching knowledge on the hierarchy of molecular orbital theories of electron pairs.
We present a method for the approximate solution of the coupled-cluster (CC) equation, based upon the singular value decomposition (SVD). The key idea of this method is that we use SVD for the cluster amplitudes to exploit the physically important states in the CC equation. This method enables us to significantly reduce the computational requirements for CC calculations without losing the essence of the method. Relationships to the density matrix renormalization group theory and the local correlation methods are mentioned. We perform pilot calculations on some atoms and molecules to investigate the applicability of the method.
Wave functions constructed from electron-pair states can accurately model strong electron correlation effects and are promising approaches especially for larger many-body systems. In this article, we analyze the nature and the type of electron correlation effects that can be captured by wave functions restricted to electron-pair states. We focus on the pair-coupled-cluster doubles (pCCD) ansatz also called the antisymmetric product of the 1-reference orbital geminal (AP1roG) method, combined with an orbital optimization protocol presented in Boguslawski et al. [Phys. Rev. B 89, 201106(R) (2014)], whose performance is assessed against electronic structures obtained form density-matrix renormalization-group reference data. Our numerical analysis covers model systems for strong correlation: the one-dimensional Hubbard model with a periodic boundary condition as well as metallic and molecular hydrogen rings. Specifically, the accuracy of pCCD/AP1roG is benchmarked using the single-orbital entropy, the orbital-pair mutual information, as well as the eigenvalue spectrum of the one-orbital and two-orbital reduced density matrices. Our study indicates that contributions from singly occupied states become important in the strong correlation regime which highlights the limitations of the pCCD/AP1roG method. Furthermore, we examine the effect of orbital rotations within the pCCD/AP1roG model on correlations between orbital pairs.
The central problem in the modern electronic structure theory is the calculation of cor- relation energy, possibly by an approach that would account for both static and dynamic correlation in an efficient, balanced and accurate way. In this thesis, I present a collection of methods that combine the effective treatment of dynamic correlation by the coupled cluster theory with density matrix renormalization group, a well-established technique for calculations of strongly correlated systems. The connection between them is achieved via the tailored coupled clusters (TCC) ansatz, which conveniently does not impose any ad- ditional computational costs. After the successful initial assessment, we developed more efficient implementations of these methods by employing the local approaches based on pair natural orbitals. This way, we extended the range of possible applications to larger systems with thousands of basis functions. To assess the accuracy of TCC as well as its local counterparts, we performed a variety of benchmark calculations ranging from small, yet challenging systems such as the nitrogen molecule or tetramethyleneethane diradical, to larger molecules like oxo-Mn(Salen) or Fe(II)-porphyrin model. 1
Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.
The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.