Efficient Formulation of Ab Initio Quantum Embedding in Periodic Systems: Dynamical Mean-Field Theory

Accurate Excitation Energies of Point Defects from Fast Particle-Particle Random Phase Approximation Calculations

We present an efficient particle-particle random phase approximation (ppRPA) approach that predicts accurate excitation energies of point defects, including the nitrogen-vacancy (NV-) and silicon-vacancy (SiV0) centers in diamond and the divacancy center (VV0) in 4H silicon carbide, with errors of ±0.2 eV compared with experimental values. Starting from the (N + 2)-electron ground state calculated with density functional theory (DFT), the ppRPA excitation energies of the N-electron system are calculated as the differences between the two-electron removal energies of the (N + 2)-electron system. We demonstrate that the ppRPA excitation energies converge rapidly with a few hundred canonical active-space orbitals. We also show that active-space ppRPA has weak DFT starting-point dependence and is significantly cheaper than the corresponding ground-state DFT calculation. This work establishes ppRPA as an accurate and low-cost tool for investigating excited-state properties of point defects and opens up new opportunities for applications of ppRPA to periodic bulk materials.

Performance of periodic EOM-CCSD for bandgaps of inorganic semiconductors and insulators

We calculate bandgaps of 12 inorganic semiconductors and insulators composed of atoms from the first three rows of the Periodic Table using periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Our calculations are performed with atom-centered triple-zeta basis sets and up to 64 k-points in the Brillouin zone. We analyze the convergence behavior with respect to the number of orbitals and number of k-points sampled using composite corrections and extrapolations to produce our final values. When accounting for electron-phonon corrections to experimental bandgaps, we find that EOM-CCSD has a mean signed error of -0.12 eV and a mean absolute error of 0.42 eV; the largest outliers are C (error of -0.93 eV), BP (-1.00 eV), and LiH (+0.78 eV). Surprisingly, we find that the more affordable partitioned EOM-MP2 theory performs as well as EOM-CCSD.

Machine Learning Many-Body Green's Functions for Molecular Excitation Spectra

We present a machine learning (ML) framework for predicting Green's functions of molecular systems, from which photoemission spectra and quasiparticle energies at quantum many-body level can be obtained. Kernel ridge regression is adopted to predict self-energy matrix elements on compact imaginary frequency grids from static and dynamical mean-field electronic features, which gives direct access to real-frequency many-body Green's functions through analytic continuation and Dyson's equation. Feature and self-energy matrices are represented in a symmetry-adapted intrinsic atomic orbital plus projected atomic orbital basis to enforce rotational invariance. We demonstrate good transferability and high data efficiency of the proposed ML method across molecular sizes and chemical species by showing accurate predictions of density of states (DOS) and quasiparticle energies at the level of many-body perturbation theory (GW) or full configuration interaction. For the ML model trained on 48 out of 1995 molecules randomly sampled from the QM7 and QM9 data sets, we report the mean absolute errors of ML-predicted highest occupied and lowest unoccupied molecular orbital energies to be 0.13 and 0.10 eV, respectively, compared to GW@PBE0. We further showcase the capability of this method by applying the same ML model to predict DOS for significantly larger organic molecules with up to 44 heavy atoms.

Block2: A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Optical Properties of Neutral F Centers in Bulk MgO with Density Matrix Embedding

The optical spectra of neutral oxygen vacancies (F0 centers) in the bulk MgO lattice are investigated using density matrix embedding theory. The impurity Hamiltonian is solved with the complete active space self-consistent field and second-order n-electron valence state perturbation theory (NEVPT2-DMET) multireference methods. To estimate defect-localized vertical excitation energies at the nonembedding and thermodynamic limits, a double extrapolation scheme is employed. The extrapolated NEVPT2-DMET vertical excitation energy value of 5.24 eV agrees well with the experimental absorption maxima at 5.03 eV, whereas the excitation energy value of 2.89 eV at the relaxed triplet defect-localized state geometry overestimates the experimental emission at 2.4 eV by only nearly 0.5 eV, indicating the involvement of the triplet-singlet decay pathway.

Nonadiabatic Potential Energy Surfaces for a Molecule on a Surface as Found by Constrained Complete Active Space Theory

In order to study electron-transfer mediated chemical processes on a metal surface, one requires not one but two potential energy surfaces (one ground state and one excited state) as in Marcus theory. In this letter, we report that a novel, dynamically weighted, state-averaged constrained CASSCF(2,2) (DW-SA-cCASSCF(2,2)) can produce such surfaces for the Anderson impurity model. Both ground and excited state potentials are smooth, they incorporate states with a charge transfer character, and the accuracy of the ground state surface can be verified for some model problems by renormalization group theory. Future development of gradients and nonadiabatic derivative couplings should allow for the study of nonadiabatic dynamics for molecules near metal surfaces.

Periodic Bootstrap Embedding

Bootstrap embedding (BE) is a recently developed electronic structure method that has shown great success at treating electron correlation in molecules. Here, we extend BE to treat surfaces and solids where the wave function is represented in periodic boundary conditions using reciprocal space sums (i.e., k-point sampling). The major benefit of this approach is that the resulting fragment Hamiltonians carry no explicit dependence on the reciprocal space sums, allowing one to apply traditional nonperiodic electronic structure codes to the fragments even though the entire system requires careful consideration of periodic boundary conditions. Using coupled cluster singles and doubles (CCSD) as an example method to solve the fragment Hamiltonians, we present minimal basis set CCSD-in-HF results on 1D conducting polymers. We show that periodic BE-CCSD can typically recover ∼99.9% of the electron correlation energy. We further demonstrate that periodic BE-CCSD is feasible even for complex donor-acceptor polymers of interest to organic solar cells─despite the fact that the monomers are sufficiently large that even a Γ-point periodic CCSD calculation is prohibitive. We conclude that BE is a promising new tool for applying molecular electronic structure tools to solids and interfaces.

Effective Reconstruction of Expectation Values from Ab Initio Quantum Embedding

Quantum embedding is an appealing route to fragment a large interacting quantum system into several smaller auxiliary "cluster" problems to exploit the locality of the correlated physics. In this work, we critically review approaches to recombine these fragmented solutions in order to compute nonlocal expectation values, including the total energy. Starting from the democratic partitioning of expectation values used in density matrix embedding theory, we motivate and develop a number of alternative approaches, numerically demonstrating their efficiency and improved accuracy as a function of increasing cluster size for both energetics and nonlocal two-body observables in molecular and solid state systems. These approaches consider the N-representability of the resulting expectation values via an implicit global wave function across the clusters, as well as the importance of including contributions to expectation values spanning multiple fragments simultaneously, thereby alleviating the fundamental locality approximation of the embedding. We clearly demonstrate the value of these introduced functionals for reliable extraction of observables and robust and systematic convergence as the cluster size increases, allowing for significantly smaller clusters to be used for a desired accuracy compared to traditional approaches in ab initio wave function quantum embedding.

Electronic Structure and Ultrafast Electron Dynamics in CuO Photocatalysts Probed by Surface Sensitive Femtosecond X-ray Absorption Near-Edge Structure Spectroscopy

CuO is often employed as a photocathode for H₂ evolution and CO₂ reduction, but observed efficiency is still far below the theoretical limit. To bridge the gap requires understanding the CuO electronic structure; however, computational efforts lack consensus on the orbital character of the photoexcited electron. In this study, we measure the femtosecond XANES spectra of CuO at the Cu M_2,3 and O L₁ edges to track the element-specific dynamics of electrons and holes. Results show that photoexcitation represents an O 2p to Cu 4s charge transfer state indicating the conduction band electron has primarily Cu 4s character. We also observe ultrafast mixing of Cu 3d and 4s conduction band states mediated by coherent phonons, with Cu 3d character of the photoelectron reaching a maximum of 16%. This is the first observation of the photoexcited redox state in CuO, and results provide a benchmark for theory where electronic structure modeling still relies heavily on model-dependent parametrization.

Periodic Coupled-Cluster Green's Function for Photoemission Spectra of Realistic Solids

We present an efficient implementation of the coupled-cluster Green's function (CCGF) method for simulating photoemission spectra of periodic systems. We formulate the periodic CCGF approach with Brillouin zone sampling in the Gaussian basis at the coupled-cluster singles and doubles (CCSD) level. To enable CCGF calculations of realistic solids, we propose an active-space self-energy correction scheme by combining CCGF with the cheaper many-body perturbation theory (GW) and implement the model order reduction (MOR) frequency interpolation technique. We find that the active-space self-energy correction and MOR techniques significantly reduce the computational cost of CCGF while maintaining the high accuracy. We apply the developed CCGF approaches to compute spectral properties and band structure of silicon (Si) and zinc oxide (ZnO) crystals using triple-ζ Gaussian basis sets and medium-size k-point sampling and find good agreement with experimental measurements.

Systematic electronic structure in the cuprate parent state from quantum many-body simulations

The quantitative description of correlated electron materials remains a modern computational challenge. We demonstrate a numerical strategy to simulate correlated materials at the fully ab initio level beyond the solution of effective low-energy models and apply it to gain a detailed microscopic understanding across a family of cuprate superconducting materials in their parent undoped states. We uncover microscopic trends in the electron correlations and reveal the link between the material composition and magnetic energy scales through a many-body picture of excitation processes involving the buffer layers. Our work illustrates a path toward a quantitative and reliable understanding of more complex states of correlated materials at the ab initio many-body level.

Correlation-Consistent Gaussian Basis Sets for Solids Made Simple

The rapidly growing interest in simulating condensed-phase materials using quantum chemistry methods calls for a library of high-quality Gaussian basis sets suitable for periodic calculations. Unfortunately, most standard Gaussian basis sets commonly used in molecular simulation show significant linear dependencies when used in close-packed solids, leading to severe numerical issues that hamper the convergence to the complete basis set (CBS) limit, especially in correlated calculations. In this work, we revisit Dunning's strategy for construction of correlation-consistent basis sets and examine the relationship between accuracy and numerical stability in periodic settings. We find that limiting the number of primitive functions avoids the appearance of problematic small exponents while still providing smooth convergence to the CBS limit. As an example, we generate double-, triple-, and quadruple-ζ correlation-consistent Gaussian basis sets for periodic calculations with Goedecker-Teter-Hutter (GTH) pseudopotentials. Our basis sets cover the main-group elements from the first three rows of the periodic table. Especially for atoms on the left side of the periodic table, our basis sets are less diffuse than those used in molecular calculations. We verify the fast and reliable convergence to the CBS limit in both Hartree-Fock and post-Hartree-Fock (MP2) calculations, using a diverse test set of 19 semiconductors and insulators.

Exploring Coupled Cluster Green's Function as a Method for Treating System and Environment in Green's Function Embedding Methods

Within the self-energy embedding theory (SEET) framework, we study the coupled cluster Green's function (GFCC) method in two different contexts: as a method to treat either the system or the environment present in the embedding construction. Our study reveals that when GFCC is used to treat the environment we do not see improvement in total energies in comparison to the coupled cluster method itself. To rationalize this puzzling result, we analyze the performance of GFCC as an impurity solver with a series of transition metal oxides. These studies shed light on the strength and weaknesses of such a solver and demonstrate that such a solver gives very accurate results when the size of the impurity is small. We investigate if it is possible to achieve a systematic accuracy of the embedding solution when we increase the size of the impurity problem. We found that in such a case, the performance of the solver worsens, both in terms of finding the ground state solution of the impurity problem and the self-energies produced. We concluded that increasing the rank of GFCC solver is necessary to be able to enlarge impurity problems and achieve a reliable accuracy. We also have shown that natural orbitals from weakly correlated perturbative methods are better suited than symmetrized atomic orbitals (SAO) when the total energy of the system is the target quantity.

Designing and controlling the properties of transition metal oxide quantum materials

This Perspective addresses the design, creation, characterization and control of synthetic quantum materials with strong electronic correlations. We show how emerging synergies between theoretical/computational approaches and materials design/experimental probes are driving recent advances in the discovery, understanding and control of new electronic behaviour in materials systems with interesting and potentially technologically important properties. The focus here is on transition metal oxides, where electronic correlations lead to a myriad of functional properties including superconductivity, magnetism, Mott transitions, multiferroicity and emergent behaviour at picoscale-designed interfaces. Current opportunities and challenges are also addressed, including possible new discoveries of non-equilibrium phenomena and optical control of correlated quantum phases of transition metal oxides.

Accurate Electronic Excitation Energies in Full-Valence Active Space via Bootstrap Embedding

Fragment embedding has been widely used to circumvent the high computational scaling of using accurate electron correlation methods to describe the electronic ground states of molecules and materials. However, similar applications that utilize fragment embedding to treat electronic excited states are comparably less reported in the literature. The challenge here is twofold. First, most fragment embedding methods are most effective when the property of interest is local, but the change of the wave function upon excitation is nonlocal in general. Second, even for local excitations, an accurate estimate of, for example, the excitation energy can still be challenging owing to the need for a balanced treatment of both the ground and the excited states. In this work, we show that bootstrap embedding (BE), a fragment embedding method developed recently by our group, is promising toward describing general electronic excitations. Numerical simulations show that the excitation energies in full-valence active space (FVAS) can be well-estimated by BE to an error of ∼0.05 eV using relatively small fragments, for both local excitations and the excitations of some large dye molecules that exhibit strong charge-transfer characters. We hence anticipate BE to be a promising solution to accurately describing the excited states of large chemical systems.

DMRG on Top of Plane-Wave Kohn-Sham Orbitals: A Case Study of Defected Boron Nitride

In this paper, we analyze the numerical aspects of the inherent multireference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

Bootstrap Embedding For Large Molecular Systems

Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N³) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.

Active space approaches combining coupled-cluster and perturbation theory for ground states and excited states

We evaluate the performance of approaches that combine coupled-cluster and perturbation theory based on a predefined active space of orbitals. Coupled-cluster theory is used to treat excitations that are internal to the active space and perturbation theory is used for all other excitations, which are at least partially external to the active space. We consider a variety of schemes that differ in how the internal and external excitations are coupled. Such approaches are presented for ground states and excited states within the equation-of-motion formalism. Results are given for the ionization potentials and electron affinities of a test set of small molecules and for the correlation energy and band gap of a few periodic solids.