Phenomenological rules play a central role in the design of chemical reactions and materials with targeted properties. Typically, these are formulated heuristically in terms of non-interacting orbitals and bands yet show remarkable accuracy in predicting the complex behavior of intrinsically interacting many-body systems. While their non-interacting formulation makes them easy to interpret, it potentially hinders the development of new rules for systems governed by strong correlations, such as transition-metal-based materials. In this work, we present a rigorous framework that allows bridging between fully interacting, even strongly correlated, systems and an effective one-body picture in terms of quasiparticles. Furthermore, we present a computational strategy to efficiently and accurately access the main components of such a description: the embedding approximation of the ghost Gutzwiller ansatz. We illustrate the capabilities of this quasiparticle formulation on the Woodward-Hoffmann rules and apply their reformulated version to toy "reactions," which exemplify the main scenarios covered by them.
Phenomenological rules play a central role in the design of chemical reactions and materials with targeted properties. Typically, these are formulated heuristically in terms of non-interacting orbitals and bands yet show remarkable accuracy in predicting the complex behavior of intrinsically interacting many-body systems. While their non-interacting formulation makes them easy to interpret, it potentially hinders the development of new rules for systems governed by strong correlations, such as transition-metal-based materials. In this work, we present a rigorous framework that allows bridging between fully interacting, even strongly correlated, systems and an effective one-body picture in terms of quasiparticles. Furthermore, we present a computational strategy to efficiently and accurately access the main components of such a description: the embedding approximation of the ghost Gutzwiller ansatz. We illustrate the capabilities of this quasiparticle formulation on the Woodward-Hoffmann rules and apply their reformulated version to toy "reactions," which exemplify the main scenarios covered by them.
Strong/static electronic correlation mediates the emergence of remarkable phases of matter, and underlies the exceptional reactivity properties in transition metal-based catalysts. Modeling strongly correlated molecules and solids calls for multi-reference Ansätze, which explicitly capture the competition of energy scales characteristic of such systems. With the efficient computational screening of correlated solids in mind, the ghost Gutzwiller (gGut) Ansatz has been recently developed. This is a variational Ansatz which can be formulated as a self-consistent embedding approach, describing the system within a non-interacting, quasiparticle model, yet providing accurate spectra in both low and high energy regimes. Crucially, small fragments of the system are identified as responsible for the strong correlation, and are therefore enhanced by adding a set of auxiliary orbitals, the ghosts. These capture many-body correlations through one-body fluctuations and subsequent out-projection when computing physical observables. gGut has been shown to accurately describe multi-orbital lattice models at modest computational cost. In this work, we extend the gGut framework to strongly correlated molecules, for which it holds special promise. Indeed, despite the asymmetric embedding treatment, the quasiparticle Hamiltonian effectively describes all major sources of correlation in the molecule: strong correlation through the ghosts in the fragment, and dynamical correlation through the quasiparticle description of its environment. To adapt the gGut Ansatz for molecules, we address the fact that, unlike in the lattice model previously considered, electronic interactions in molecules are not local. Hence, we explore a hierarchy of approximations of increasing accuracy capturing interactions between fragments and environment, and within the environment, and discuss how these affect the embedding description of correlations in the whole molecule. We will compare the accuracy of the gGut model with established methods to capture strong correlation within active space formulations, and assess the realistic use of this novel approximation to the theoretical description of correlated molecular clusters.
We present benchmark calculations of the Anderson lattice model based on the recently developed ``ghost Gutzwiller approximation''. Our analysis shows that, in some parameters regimes, the predictions of the standard Gutzwiller approximation can be incorrect by orders of magnitude for this model. We show that this is caused by the inability of this method to describe simultaneously the Mott physics and the hybridization between correlated and itinerant degrees of freedom---whose interplay often governs the metal-insulator transition in real materials. Finally, we show that the ghost Gutzwiller approximation solves this problem, providing us with results in remarkable agreement with dynamical mean-field theory throughout the entire phase diagram, while being much less computationally demanding. We provide an analytical explanation of these findings and discuss their implications within the context of ab-initio computation of strongly correlated matter.
Quantum embedding (QE) methods such as the ghost Gutzwiller approximation (gGA) offer a powerful approach to simulating strongly correlated systems, but come with the computational bottleneck of computing the ground state of an auxiliary embedding Hamiltonian (EH) iteratively. In this work, we introduce an active learning (AL) framework integrated within the gGA to address this challenge. The methodology is applied to the single-band Hubbard model and results in a significant reduction in the number of instances where the EH must be solved. Through a principal component analysis (PCA), we find that the EH parameters form a low-dimensional structure that is largely independent of the geometric specifics of the systems, especially in the strongly correlated regime. Our AL strategy enables us to discover this low-dimensionality structure on the fly, while leveraging it for reducing the computational cost of gGA, laying the groundwork for more efficient simulations of complex strongly correlated materials. Published by the American Physical Society 2024
Dynamical and variational frameworks have long been viewed as distinct paradigms. In particular, in quantum embedding (QE) frameworks, dynamical mean-field theory (DMFT) captures nonperturbative dynamical correlations through a frequency-dependent self-energy, while the Gutzwiller approximation (GA) is formulated in terms of a variationally optimized ground-state wavefunction. Here we bridge these perspectives, proving that the ghost-Gutzwiller approximation (ghost-GA), which also admits a density-matrix-matching QE formulation known as ghost density matrix embedding theory (ghost-DMET), becomes strictly equivalent to DMFT in the limit of infinitely many auxiliary bath modes. This formal unification has immediate consequences. In particular, it yields a rigorous finite-temperature extension of ghost-GA and shows that the physical Green's function can be determined from static expectation values of the embedding Hamiltonians, providing a route to computational studies of competing phases in strongly correlated matter with DMFT-level accuracy, while bypassing the need to calculate dynamical spectra with conventional impurity solvers. More broadly, it shows that the variational ghost-GA, the density-matrix-matching ghost-DMET formulation, and the dynamical DMFT description are not separate constructions, but complementary formulations of the same QE structure, thereby providing a concrete formal basis for future controlled extensions beyond DMFT.
Conspectus Ab initio modeling of matter has become a pillar of chemical research: with ever-increasing computational power, simulations can be used to accurately predict, for example, chemical reaction rates, electronic and mechanical properties of materials, and dynamical properties of liquids. Many competing quantum mechanical methods have been developed over the years that vary in computational cost, accuracy, and scalability: density functional theory (DFT), the workhorse of solid-state electronic structure calculations, features a good compromise between accuracy and speed. However, approximate exchange-correlation functionals limit DFT's ability to treat certain phenomena or states of matter, such as charge-transfer processes or strongly correlated materials. Furthermore, conventional DFT is purely a ground-state theory: electronic excitations are beyond its scope. Excitations in molecules are routinely calculated using time-dependent DFT linear response; however applications to condensed matter are still limited. By contrast, many-electron wavefunction methods aim for a very accurate treatment of electronic exchange and correlation. Unfortunately, the associated computational cost renders treatment of more than a handful of heavy atoms challenging. On the other side of the accuracy spectrum, parametrized approaches like tight-binding can treat millions of atoms. In view of the different (dis-)advantages of each method, the simulation of complex systems seems to force a compromise: one is limited to the most accurate method that can still handle the problem size. For many interesting problems, however, compromise proves insufficient. A possible solution is to break up the system into manageable subsystems that may be treated by different computational methods. The interaction between subsystems may be handled by an embedding formalism. In this Account, we review embedded correlated wavefunction (CW) approaches and some applications. We first discuss our density functional embedding theory, which is formally exact. We show how to determine the embedding potential, which replaces the interaction between subsystems, at the DFT level. CW calculations are performed using a fixed embedding potential, that is, a non-self-consistent embedding scheme. We demonstrate this embedding theory for two challenging electron transfer phenomena: (1) initial oxidation of an aluminum surface and (2) hot-electron-mediated dissociation of hydrogen molecules on a gold surface. In both cases, the interaction between gas molecules and metal surfaces were treated by sophisticated CW techniques, with the remainder of the extended metal surface being treated by DFT. Our embedding approach overcomes the limitations of conventional Kohn-Sham DFT in describing charge transfer, multiconfigurational character, and excited states. From these embedding simulations, we gained important insights into fundamental processes that are crucial aspects of fuel cell catalysis (i.e., O2 reduction at metal surfaces) and plasmon-mediated photocatalysis by metal nanoparticles. Moreover, our findings agree very well with experimental observations, while offering new views into the chemistry. We finally discuss our recently formulated potential-functional embedding theory that provides a seamless, first-principles way to include back-action onto the environment from the embedded region.
The ghost Gutzwiller ansatz (gGut) embedding technique was shown to achieve comparable accuracy to the gold standard dynamical mean-field theory method in simulating real material properties, yet at a much lower computational cost. Despite that, gGut is limited by the algorithmic bottleneck of computing the density matrix of the underlying effective embedding model, a quantity which must be converged within a self-consistent embedding loop. We develop a hybrid quantum-classical gGut technique which computes the ground state properties of embedding Hamiltonians with the help of a quantum computer, using the sample-based quantum-selected configuration interaction (QSCI) algorithm. We study the applicability of SCI-based methods to the evaluation of the density of states for single-band Anderson impurity models within gGut and find that such ground states of interest become sufficiently sparse in the CI basis as the number of ghost orbitals is increased. Further, we investigate the performance of QSCI using local unitary cluster Jastrow (LUCJ) variational quantum states in combination with a circuit cutting technique, prepared on IQM's quantum hardware for system sizes of up to 11 ghost orbitals, equivalent to 24 qubits. We report converged gGut calculations which correctly capture the metal-to-insulator phase transition in the Fermi-Hubbard model on the Bethe lattice by using quantum samples to build an SCI basis with as little as 1% of the total CI basis states.