arXiv:2012.06165 · The Journal of Chemical Physics

Local embedding of coupled cluster theory into the random phase approximation using plane waves

Tobias Schäfer, Florian Libisch, Georg Kresse, Andreas Grüneis

Published Jan 6, 2021• 53 citations

Abstract

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.

DOI: 10.1063/5.0036363

Related papers

Random phase approximation-based local natural orbital coupled cluster theory

Ruiheng Song, Xiliang Gong, Aamy Bakry +1 more

Practical applications of fragment embedding and closely related local correlation methods critically depend on a judicious choice of a low-level theory to define the local embedding subspace and to capture long-range electrostatic and correlation effects outside the embedding region. Second-order Møller-Plesset perturbation theory (MP2) is by far the most widely used correlated low-level theory; however, its applicability becomes questionable in systems where MP2 is known to fail either quantitatively or qualitatively. In this work, we present the random phase approximation (RPA) as a promising alternative low-level theory to MP2 within the local natural orbital-based coupled-cluster (LNO-CC) framework. We demonstrate that RPA-based LNO-CC closely matches the performance of its MP2-based counterpart for systems with sizable energy gaps, while delivering significantly faster convergence toward the canonical coupled-cluster limit for metallic systems, particularly as the thermodynamic limit is approached. These results highlight the critical role of the low-level theory in fragment embedding and local correlation methods and identify RPA as a compelling alternative to the commonly used MP2.

arXiv:2601.00131Dec 31, 20250 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.06424Frontiers in MaterialsJun 5, 2019138 citations

A static quantum embedding scheme based on coupled cluster theory

Avijit Shee, Fabian M. Faulstich, K. Birgitta Whaley +2 more

We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

arXiv:2404.09078The Journal of Chemical PhysicsOct 24, 202414 citations

Towards efficient and accurate ab initio solutions to periodic systems via transcorrelation and coupled cluster theory

Ke Liao, Thomas Schraivogel, Hongjun Luo +2 more

We propose a streamlined combination of the transcorrelation (TC) and coupled cluster (CC) theory, which not only increases the convergence rate with respect to the basis set but also extends the applicability of the lowest-order CC approximations to strongly correlated systems. Using the three-dimensional uniform electron gas as a model for real periodic solids and highly accurate ground-state energies from the state-of-the-art quantum Monte Carlo methods as benchmark, we show that, with the correct physical insights incorporated, our TC-CC methods gain drastic improvements in accuracy compared with canonical CC methods, and the errors of our methods remain $\ifmmode\pm\else\textpm\fi{}0.001$ Hartree per electron (Ha/electron) across a wide range of densities. With the greatly improved efficiencies and the balanced and accurate performances at both low and high densities, our methods present great promise for applications on real solids with weak and strong correlations.

arXiv:2103.03176Physical Review ResearchJul 20, 202141 citations

Correlation-Converged Virtual Orbitals for Accurate and Efficient Quantum Molecular Simulations

Qian Wang, Calvin Ku, Jyh-Pin Chou +3 more

Density functional theory with plane-wave basis sets is widely employed in computational materials science, including applications to isolated molecular systems. However, the inadequate description of electron correlation remains a fundamental limitation. Accurate correlation treatments based on many-body Hamiltonians require reliable representations of both occupied and virtual orbitals, yet virtual orbitals are often poorly described in conventional computational schemes, resulting in reduced accuracy. In this work, we introduce localized correlation-converged virtual orbitals (LCCVOs) as an efficient basis for constructing accurate many-body Hamiltonians in molecular systems. Using a substantially reduced number of orbitals, the LCCVO framework yields dissociation energies for singlet, doublet, and triplet molecules that are comparable to, and in many cases exceed, those obtained with high-level correlation-consistent basis sets such as cc-pVXZ (X = D, T, Q, 5). These results demonstrate the efficiency, scalability, and robustness of the LCCVO approach for high-accuracy quantum chemical calculations.

physics.chem-phCross-listedarXiv:2604.16928Apr 21, 2026

Wavelength‐decomposition‐based embedded cluster density approximation for systems with nonlocal electron correlation

Chen Huang

Abstract Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one‐dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short‐wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long‐wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn‐Sham orbitals. This new method is termed λ‐ECDA, where λ indicates the wavelength decomposition. The performance of λ‐ECDA is examined on a one‐dimensional model system: a H24 chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron‐electron and electron‐ion interactions, and slightly stronger nuclear charges (1.2e) are assigned to the pseudo‐H atoms. Bond stretching energies, RPA correlation potentials, and Kohn‐Sham eigenvalues predicted by λ‐ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy.

International Journal of Quantum ChemistryJul 4, 20204 citations