Multicomponent systems are defined as chemical systems that require a quantum mechanical description of two or more different types of particles. Non-Born-Oppenheimer electron-nuclear interactions in molecules, electron-hole interactions in electronically excited nanoparticles, and electron-positron interactions are examples of physical systems that require a multicomponent quantum mechanical formalism. The central challenge in the theoretical treatment of multicomponent systems is capturing the many-body correlation effects that exist not only between particles of identical types (electron-electron) but also between particles of different types (electron-nuclear and electron-hole). In this work, the development and implementation of multicomponent coupled-cluster (mcCC) theory for treating particle-particle correlation in multicomponent systems is presented. This method provides a balanced treatment of many-particle correlation in a general multicomponent system while maintaining a size-consistent and size-extensive formalism. The coupled-cluster ansatz presented here is the extension of the electronic structure CCSD formulation for multicomponent systems and is defined as $\vert Ψ_\mathrm{mcCC} \rangle = e^{T_1^\mathrm{I}+T_2^\mathrm{I}+T_1^\mathrm{II}+T_2^\mathrm{II}+T_{11}^\mathrm{I,II}+T_{12}^\mathrm{I,II}+T_{21}^\mathrm{I,II}+T_{22}^\mathrm{I,II}}\vert 0^\mathrm{I}0^\mathrm{II}\rangle$. The applicability of the mcCC method was demonstrated by computing biexciton binding energies for multiexcitonic systems and benchmarking the results against full configuration interaction calculations. The results demonstrated that connected cluster operators that generate simultaneous excitation in type I and type II space are critical for capturing electron-hole correlation in multiexcitonic systems.
Multicomponent systems are defined as chemical systems that require a quantum mechanical description of two or more different types of particles. Non-Born-Oppenheimer electron-nuclear interactions in molecules, electron-hole interactions in electronically excited nanoparticles, and electron-positron interactions are examples of physical systems that require a multicomponent quantum mechanical formalism. The central challenge in the theoretical treatment of multicomponent systems is capturing the many-body correlation effects that exist not only between particles of identical types (electron-electron) but also between particles of different types (electron-nuclear and electron-hole). In this work, the development and implementation of multicomponent coupled-cluster (mcCC) theory for treating particle-particle correlation in multicomponent systems is presented. This method provides a balanced treatment of many-particle correlation in a general multicomponent system while maintaining a size-consistent and size-extensive formalism. The coupled-cluster ansatz presented here is the extension of the electronic structure CCSD formulation for multicomponent systems and is defined as ∣ΨmcCC⟩=eT1I+T2I+T1II+T2II+T11I,II+T12I,II+T21I,II+T22I,II∣0I0II⟩. The applicability of the mcCC method was demonstrated by computing biexciton binding energies for multiexcitonic systems and benchmarking the results against full configuration interaction calculations. The results demonstrated that connected cluster operators that generate simultaneous excitation in type I and type II space are critical for capturing electron-hole correlation in multiexcitonic systems.
The variational quantum eigensolver (VQE) algorithm combined with the unitary coupled cluster (UCC) ansatz has been developed for the quantum computation of molecular energies and wave functions within the Born–Oppenheimer approximation. Herein, this approach is extended to multicomponent systems to enable the quantum mechanical treatment of more than one type of particle, such as electrons and positrons or electrons and nuclei, without invoking the Born–Oppenheimer approximation. Specifically, we introduce the multicomponent unitary coupled cluster (mcUCC) method combined with the VQE algorithm for the calculation of ground-state energies and wave functions as well as the multicomponent equation-of-motion (mcEOM) method for the calculation of excitation energies. These methods are developed within the nuclear-electronic orbital (NEO) framework and are formulated in the qubit basis to enable implementations on quantum computers. Moreover, these methods are used to calculate the ground-state energy and excitation energies of positronium hydride, where both electrons and the positron are treated quantum mechanically, as well as the H 2 molecule, where both electrons and one proton are treated quantum mechanically. These applications validate the implementation and provide benchmark data for future calculations. The errors due to Trotterization of the mcUCC ansatz are also analyzed. This formalism, as well as the accompanying computer code, will serve as the basis for applications to more complex multicomponent systems, such as simulations of photoinduced nonadiabatic molecular processes, on both classical and quantum computers.
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.
We present in this contribution the first application of local correlation in the context of multicomponent methods. Multicomponent approaches allow for the targeted simulation of electrons together with other Fermions (most commonly protons) as quantum particles. These methods have become increasingly popular over the last years, particularly for the description of nuclear quantum effects (in strong hydrogen bonds, proton tunneling, and many more). However, most implementations are still based on canonical formulations of wave function theory, which we know for decades to be computationally inefficient for capturing dynamical correlation effects. Local correlation approaches, particularly with the use of pair natural orbitals (PNOs), enable asymptotically linear scaling of computational costs with very little impact on the overall accuracy. In this context, the efficient use of density fitting approximations in the integral calculation proves essential. We start by discussing our implementation of density-fitted NEO-MP2 and NEO-PNO-LMP2, upgrading the electronic correlation treatment up to PNO local coupled cluster level of theory. Several challenging examples are provided to benchmark the method in terms of accuracy as well as computational cost scaling. Following appropriate protocols, anharmonic corrections to localized X-H stretches can be applied routinely with little computational overhead.
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
We introduce a multicomponent unitary coupled cluster framework for quantum simulations of molecular systems that incorporate both electronic and nuclear quantum effects beyond the Born-Oppenheimer approximation. Using the nuclear-electronic orbital formalism, we construct mcUCC ansätze for positronium hydride and molecular hydrogen with a quantum proton, and analyze hardware requirements for different excitation truncations. To further reduce resource costs effectively, we employ the local unitary cluster Jastrow ansatz and implement it experimentally on IBM Q's Heron superconducting hardware. With the Physics-Inspired Extrapolation error mitigation protocol, the computed ground-state energies remain within chemical accuracy, consistent with the stated uncertainty level. These results provide the first demonstration of error-mitigated multicomponent correlated simulations on quantum hardware and outline a path toward scalable algorithms unifying electronic and nuclear degrees of freedom.
Scope The book focuses primarily on many-body (or better, many-electron) methods for electron correlation. These include Rayleigh–Schrödinger perturbation theory (RSPT) , particularly in its diagrammatic representation (referred to as many-body perturbation theory , or MBPT ), and coupled-cluster (CC) theory; their relationship to configuration interaction (CI) is included. Further extensions address properties other than the energy, and also excited states and multireference CC and MBPT methods. The many-body algebraic and diagrammatic methods used in electronic structure theory have their origin in quantum field theory and in the study of nuclear matter and nuclear structure. The second-quantization formalism was first introduced in a treatment of quantized fields by Dirac (1927) and was extended to fermion systems by Jordan and Klein (1927) and by Jordan and Wigner (1928). This formalism is particularly useful in field theory, in scattering problems and in the study of infinite systems because it easily handles problems involving infinite, indefinite or variable numbers of particles. The diagrammatic approach was introduced into field theory by Feynman (1949a,b) and applied to many-body systems by Hugenholtz (1957) and by Goldstone (1957). Many-body perturbation theory and its linked-diagram formalism were first introduced by Brueckner and Levinson (1955) and by Brueckner (1955), and were formalized by Goldstone (1957). Other important contributions to the methodology, first in field theory and then in the theory of nuclear structure, are due to Dyson (1949a,b), Wick (1950), Hubbard (1957, 1958a,b) and Frantz and Mills (1960). Applications to the electronic structure of atoms and molecules began with the work of Kelly (1963, 1964a,b, 1968), and molecular applications using finite analytical basis sets appeared in the work of Bartlett and Silver (1974a, b).