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

Physical review. B, Condensed matter

Coupled-cluster approach to electron correlation in one dimension: Cyclic polyene model in delocalized basis

Josef Paldus, M. Takahashi, R. W. H. Cho

Published Oct 15, 1984• 138 citations
Publisher

Abstract

The many-electron correlation problem for metalliclike systems with Born---von K\'arm\'an boundary conditions, as modeled by the cyclic polyene Pariser-Parr-Pople and Hubbard Hamiltonians, is examined over the entire range of the coupling constant using the coupled-pair (CP) many-electron theory based on the exponential-cluster ansatz for the exact wave function. It is shown that the standard CP theory breaks down not only in the highly correlated, but even in the intermediately correlated, regions, and the nature of its singular behavior in these regions is examined. This breakdown is linked with the increasingly important role played by the connected quadruply excited clusters, which invalidate the basic assumption of the CP theory, as the highly correlated limit and/or the extended character of the model are approached. The contribution from the quadruply excited clusters is then taken into account using the new version of the approximate coupled-pair theory corrected for connected quadruply excited clusters, called ACPQ. This approximation is almost identical with the standard approximate coupled-pair (ACP) -theory approach, in which only factorizable (with respect to hole pairs) nonlinear terms are retained, and differs from it only by a numerical factor associated with one of the nonlinear diagrams. It is shown that the ACPQ not only removes the singularities and associated convergency problems of the standard CP approaches, but, in fact, provides excellent quantitative results over the entire range of the coupling constant, yielding the exact correlation energy in the strongly correlated limit.

DOI
10.1103/physrevb.30.4267

Related papers

Coupled-cluster approach to electron correlation in one dimension. II. Cyclic polyene model in localized basis

M. Takahashi, Josef Paldus

The many-electron correlation problem for one-dimensional metalliclike systems with Born--von K\'arm\'an boundary conditions, represented by the Pariser-Parr-Pople and the Hubbard Hamiltonian cyclic polyene models, ${\mathrm{C}}_{\mathrm{N}}$${\mathrm{H}}_{\mathrm{N}}$, N=2n=4\ensuremath{\nu}+2, \ensuremath{\nu}=1,2,..., is studied using the coupled-cluster approach in the localized Wannier basis representation. Various truncation schemes for the pair clusters are examined. It is shown that already the intracell pair-cluster approximation, which can be handled analytically and yields the same expression for the correlation energy as the variational approach of Ukrainskii, provides a reasonable approximation in the entire range of the coupling constant. Using all doubly excited clusters composed of locally excited particle-hole pairs, one obtains the exact correlation energy in the fully correlated limit assuming that the coupled-pair equations are corrected for the connected quadruply excited cluster contributions. This is achieved by using the recently developed approximate coupled-pair approach with corrections for the quadruply excited clusters (ACPQ), which is almost identical, except for a numerical factor of certain diagrams, with the standard approximate coupled-pair approach. The ACPQ approach removes the singularities which otherwise plague the standard coupled-pair approach and yields very good correlation energies in the entire range of the coupling constant even when the ${n}^{3}$ doubly excited pair clusters are truncated to only n+10 locally and quasilocally excited pair clusters in the localized Wannier basis representation.

Physical review. B, Condensed matter•Apr 15, 1985•96 citations

Perturbation theory and electron correlation in extended systems: Cyclic polyene model

M. Takahashi, Josef Paldus, Jakub Čı́žek

Abstract The applicability of the finite‐order many‐body perturbation theory to the electron correlation problem in extended one‐dimensional systems is examined. The cyclic polyenes CNHN, N = 4ν + 2, ν = 1, 2, …, with the DNh geometry as described by both the Pariser–Parr–Pople and Hubbard Hamiltonians, are employed to model the metallic‐like one‐dimensional systems. The second‐order perturbation theory contributions to the correlation energy are obtained with three different partitionings of the Hamiltonian (Hückel, M⊘ller–Plesset, and Epstein–Nesbet). The third‐ and fourth‐order contributions are also calculated in special cases. A comparison with other methods is given wherever available. For the Hubbard Hamiltonian the asymptotic behavior of the perturbation theory expansion is examined numerically. It is shown that the finite‐order perturbation expansion can provide reliable results for the correlation energy of one‐dimensional systems even in the correlation region which corresponds to the spectroscopically determined physical value of the coupling constant.

International Journal of Quantum Chemistry•Dec 1, 1983•43 citations

Coupled-cluster approaches with an approximate account of triexcitations and the optimized-inner-projection technique. III. Lower bounds to the ground-state correlation energy of cyclic-polyene model systems

Piotr Piecuch B, Sohrab Zarrabian, Josef Paldus +1 more

The optimized-inner-projection (OIP) method is implemented and applied to determine the lower bounds to the ground-state correlation energy for the Pariser-Parr-Pople and Hubbard models of cyclic polyenes ${\mathrm{C}}_{\mathit{N}}$${\mathrm{H}}_{\mathit{N}}$, N=4\ensuremath{\nu}+2, \ensuremath{\nu}=3--5, in the whole range of the coupling constant. In actual calculations, the orthogonally spin-adapted form of the OIP equations is used, and full advantage is taken of the high symmetry of the cyclic polyene model, so that systems with N\ensuremath{\ge}14 can be explored. Comparison is made with the results of both full and limited configuration-interaction calculations and with the correlation energies obtained with certain approximate coupled-cluster schemes. Contrary to the six- and ten-membered rings, the OIP lower bounds to the ground-state correlation energy of cyclic polyenes with larger number of sites are very poor and cannot serve as a meaningful source of information. Since the eigenvalues of the intermediate Hamiltonian, in terms of which the OIP technique is defined, are highly degenerate even in the weakly correlated region, and in view of the random behavior of the OIP bracketing function, the OIP method becomes computationally very demanding for N\ensuremath{\ge}14 rings. Moreover, simple iterative schemes usually diverge and sophisticated root searching procedures have to be applied. The usefulness of the OIP technique for the determination of lower bounds to energy eigenvalues in larger systems is thus questioned.

Physical Review A•Nov 1, 1990•23 citations

Coupled-cluster approach to electron correlations in the two-dimensional Hubbard model

Yoshihiro Asai, Hideki Katagiri

We have studied electron correlations in the doped two-dimensional (2D) Hubbard model by using the coupled-cluster method (CCM) to investigate whether or not the method can be applied to correct the independent particle approximations actually used in ab initio band calculations. The double excitation version of the CCM, implemented using the approximate coupled pair (ACP) method, accounts for most of the correlation energies of the 2D Hubbard model in the weak $(U/t\ensuremath{\simeq}1)$ and the intermediate $U/t$ regions $(U/t\ensuremath{\simeq}4).$ The error in this case is always less than $1%.$ The ACP approximation gets less accurate for large $U/t$ $(U/t\ensuremath{\simeq}8)$ and near half-filling. Further incorporation of electron correlation effects is necessary in this region. The accuracy depends neither on the system size, nor the gap between the lowest unoccupied level and the highest occupied level in the finite size effect. Hence, the CCM may be favorably applied to ab initio band calculations on metals as well as semiconductors and insulators.

arXiv:cond-mat/9904400•Physical review. B, Condensed matter•Nov 15, 1999•6 citations

Coupled-cluster singles and doubles for extended systems

So Hirata, Rafał Podeszwa, Motoi Tobita +1 more

Coupled-cluster theory with connected single and double excitation operators (CCSD) and related approximations, such as linearized CCSD, quadratic configuration interaction with single and double excitation operators, coupled-cluster with connected double excitation operator (CCD), linearized CCD, approximate CCD, and second- and third-order many-body perturbation theories, are formulated and implemented for infinitely extended one-dimensional systems (polymers), on the basis of the periodic boundary conditions and distance-based screening of integrals, density matrix elements, and excitation amplitudes. The variation of correlation energies with the truncation radii of short- and long-range lattice sums and with the number of wave vector sampling points in the first Brillouin zone is examined for polyethylene, polyacetylene, and polyyne, and is shown to be a function of the degree of pi-electron conjugation or the fundamental band gaps. The t2 and t1 amplitudes in the atomic orbital (AO) basis are obtained by first computing the t amplitudes in the Bloch-orbital basis and subsequently back-transforming them into the AO basis. The plot of these AO-based t amplitudes as a function of unit cells also indicates that the t2 amplitudes of polyacetylene and polyyne exhibit appreciably slower decay than those of polyethylene, although the asymptotic decay behavior is invariably 1/r3. The AO-based t1 amplitudes appear to correlate strongly with the electronic structure, and they decay seemingly exponentially for polyethylene whereas they stay at a constant magnitude across the seventh nearest neighbors of polyacetylene and polyyne, which attests to far reaching effects of nondynamical electron correlation mediated by orbital rotation. Nonetheless, the unit cell contributions to the correlation energies taper below 10(-6) hartree after 15 A for all three polymers. The basis set dependence of the decay behavior of t2 amplitudes is also examined for linear hydrogen fluoride polymer (HF)infinity and linear beryllium polymer (Be)infinity employing the STO-3G, 6-31G, and 6-31G* basis sets, and proves to be rather small.

The Journal of Chemical Physics•Feb 4, 2004•162 citations

Coupled Cluster Theory in Materials Science

Igor Ying Zhang, Andreas Grüneis

The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.

arXiv:2004.06424•Frontiers in Materials•Jun 5, 2019•138 citations
The many-electron correlation problem for metalliclike systems with Born---von K'arm'an boundary conditions, as modeled by the cyclic polyene Pariser-Parr-Pople and Hubbard Hamiltonians, is examined over the entire range of the coupling constant using the coupled-pair (CP) many-electron theory based on the exponential-cluster ansatz for the exact wave function. It is shown that the standard CP theory breaks down not only in the highly correlated, but even in the intermediately correlated, regions, and the nature of its singular behavior in these regions is examined. This breakdown is linked with the increasingly important role played by the connected quadruply excited clusters, which invalidate the basic assumption of the CP theory, as the highly correlated limit and/or the extended character of the model are approached. The contribution from the quadruply excited clusters is then taken into account using the new version of the approximate coupled-pair theory corrected for connected quadruply excited clusters, called ACPQ. This approximation is almost identical with the standard approximate coupled-pair (ACP) -theory approach, in which only factorizable (with respect to hole pairs) nonlinear terms are retained, and differs from it only by a numerical factor associated with one of the nonlinear diagrams. It is shown that the ACPQ not only removes the singularities and associated convergency problems of the standard CP approaches, but, in fact, provides excellent quantitative results over the entire range of the coupling constant, yielding the exact correlation energy in the strongly correlated limit.
The many-electron correlation problem for one-dimensional metalliclike systems with Born--von K\'arm\'an boundary conditions, represented by the Pariser-Parr-Pople and the Hubbard Hamiltonian cyclic polyene models, CN{\mathrm{C}}_{\mathrm{N}}CN​HN{\mathrm{H}}_{\mathrm{N}}HN​, N=2n=4\ensuremath{\nu}+2, \ensuremath{\nu}=1,2,..., is studied using the coupled-cluster approach in the localized Wannier basis representation. Various truncation schemes for the pair clusters are examined. It is shown that already the intracell pair-cluster approximation, which can be handled analytically and yields the same expression for the correlation energy as the variational approach of Ukrainskii, provides a reasonable approximation in the entire range of the coupling constant. Using all doubly excited clusters composed of locally excited particle-hole pairs, one obtains the exact correlation energy in the fully correlated limit assuming that the coupled-pair equations are corrected for the connected quadruply excited cluster contributions. This is achieved by using the recently developed approximate coupled-pair approach with corrections for the quadruply excited clusters (ACPQ), which is almost identical, except for a numerical factor of certain diagrams, with the standard approximate coupled-pair approach. The ACPQ approach removes the singularities which otherwise plague the standard coupled-pair approach and yields very good correlation energies in the entire range of the coupling constant even when the n3{n}^{3}n3 doubly excited pair clusters are truncated to only n+10 locally and quasilocally excited pair clusters in the localized Wannier basis representation.
Abstract The applicability of the finite‐order many‐body perturbation theory to the electron correlation problem in extended one‐dimensional systems is examined. The cyclic polyenes CNHN, N = 4ν + 2, ν = 1, 2, …, with the DNh geometry as described by both the Pariser–Parr–Pople and Hubbard Hamiltonians, are employed to model the metallic‐like one‐dimensional systems. The second‐order perturbation theory contributions to the correlation energy are obtained with three different partitionings of the Hamiltonian (Hückel, M⊘ller–Plesset, and Epstein–Nesbet). The third‐ and fourth‐order contributions are also calculated in special cases. A comparison with other methods is given wherever available. For the Hubbard Hamiltonian the asymptotic behavior of the perturbation theory expansion is examined numerically. It is shown that the finite‐order perturbation expansion can provide reliable results for the correlation energy of one‐dimensional systems even in the correlation region which corresponds to the spectroscopically determined physical value of the coupling constant.
The optimized-inner-projection (OIP) method is implemented and applied to determine the lower bounds to the ground-state correlation energy for the Pariser-Parr-Pople and Hubbard models of cyclic polyenes CN{\mathrm{C}}_{\mathit{N}}CN​HN{\mathrm{H}}_{\mathit{N}}HN​, N=4\ensuremath{\nu}+2, \ensuremath{\nu}=3--5, in the whole range of the coupling constant. In actual calculations, the orthogonally spin-adapted form of the OIP equations is used, and full advantage is taken of the high symmetry of the cyclic polyene model, so that systems with N\ensuremath{\ge}14 can be explored. Comparison is made with the results of both full and limited configuration-interaction calculations and with the correlation energies obtained with certain approximate coupled-cluster schemes. Contrary to the six- and ten-membered rings, the OIP lower bounds to the ground-state correlation energy of cyclic polyenes with larger number of sites are very poor and cannot serve as a meaningful source of information. Since the eigenvalues of the intermediate Hamiltonian, in terms of which the OIP technique is defined, are highly degenerate even in the weakly correlated region, and in view of the random behavior of the OIP bracketing function, the OIP method becomes computationally very demanding for N\ensuremath{\ge}14 rings. Moreover, simple iterative schemes usually diverge and sophisticated root searching procedures have to be applied. The usefulness of the OIP technique for the determination of lower bounds to energy eigenvalues in larger systems is thus questioned.
We have studied electron correlations in the doped two-dimensional (2D) Hubbard model by using the coupled-cluster method (CCM) to investigate whether or not the method can be applied to correct the independent particle approximations actually used in ab initio band calculations. The double excitation version of the CCM, implemented using the approximate coupled pair (ACP) method, accounts for most of the correlation energies of the 2D Hubbard model in the weak $(U/t\ensuremath{\simeq}1)$ and the intermediate U/tU/tU/t regions $(U/t\ensuremath{\simeq}4).$ The error in this case is always less than 11%.1 The ACP approximation gets less accurate for large U/tU/tU/t $(U/t\ensuremath{\simeq}8)$ and near half-filling. Further incorporation of electron correlation effects is necessary in this region. The accuracy depends neither on the system size, nor the gap between the lowest unoccupied level and the highest occupied level in the finite size effect. Hence, the CCM may be favorably applied to ab initio band calculations on metals as well as semiconductors and insulators.
Coupled-cluster theory with connected single and double excitation operators (CCSD) and related approximations, such as linearized CCSD, quadratic configuration interaction with single and double excitation operators, coupled-cluster with connected double excitation operator (CCD), linearized CCD, approximate CCD, and second- and third-order many-body perturbation theories, are formulated and implemented for infinitely extended one-dimensional systems (polymers), on the basis of the periodic boundary conditions and distance-based screening of integrals, density matrix elements, and excitation amplitudes. The variation of correlation energies with the truncation radii of short- and long-range lattice sums and with the number of wave vector sampling points in the first Brillouin zone is examined for polyethylene, polyacetylene, and polyyne, and is shown to be a function of the degree of pi-electron conjugation or the fundamental band gaps. The t2 and t1 amplitudes in the atomic orbital (AO) basis are obtained by first computing the t amplitudes in the Bloch-orbital basis and subsequently back-transforming them into the AO basis. The plot of these AO-based t amplitudes as a function of unit cells also indicates that the t2 amplitudes of polyacetylene and polyyne exhibit appreciably slower decay than those of polyethylene, although the asymptotic decay behavior is invariably 1/r3. The AO-based t1 amplitudes appear to correlate strongly with the electronic structure, and they decay seemingly exponentially for polyethylene whereas they stay at a constant magnitude across the seventh nearest neighbors of polyacetylene and polyyne, which attests to far reaching effects of nondynamical electron correlation mediated by orbital rotation. Nonetheless, the unit cell contributions to the correlation energies taper below 10(-6) hartree after 15 A for all three polymers. The basis set dependence of the decay behavior of t2 amplitudes is also examined for linear hydrogen fluoride polymer (HF)infinity and linear beryllium polymer (Be)infinity employing the STO-3G, 6-31G, and 6-31G* basis sets, and proves to be rather small.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.