Entangled atoms, exemplified by the two singlet-coupled atoms in a dissociating nitrogen molecule, are key to the unique properties of qubits, quantum sensors, and other quantum materials. Simulating entangled atoms and strong electron correlation is a prerequisite for modeling quantum materials and remains an outstanding challenge for theory. I introduce a correlated wave function approach capable of simulating strong electron correlation in large active spaces. Projected-interacting full configuration interaction plus density functional theory PiFCI+DFT combines near-exact correlated wave functions of multiple partially interacting model systems, each corrected by a formally exact density functional. This approach can access large active spaces and visualize entanglement and strong correlation while maintaining competitive accuracy for molecular properties. I demonstrate PiFCI+DFT for general main-group chemistry, benchmark dissociation curves, and a model for a large entangled-atom system: the spin-singlet ground state of a square lattice of 100 nitrogen atoms. This system’s 300-electron active space is out of reach for existing ab initio multireference methods. PiFCI+DFT accurately recovers the strongly correlated dissociation limit at a modest computational cost, motivating further development for modeling quantum materials.
Entangled atoms, exemplified by the two singlet-coupled atoms in a dissociating nitrogen molecule, are key to the unique properties of qubits, quantum sensors, and other quantum materials. Simulating entangled atoms and strong electron correlation is a prerequisite for modeling quantum materials and remains an outstanding challenge for theory. I introduce a correlated wave function approach capable of simulating strong electron correlation in large active spaces. Projected-interacting full configuration interaction plus density functional theory PiFCI+DFT combines near-exact correlated wave functions of multiple partially interacting model systems, each corrected by a formally exact density functional. This approach can access large active spaces and visualize entanglement and strong correlation while maintaining competitive accuracy for molecular properties. I demonstrate PiFCI+DFT for general main-group chemistry, benchmark dissociation curves, and a model for a large entangled-atom system: the spin-singlet ground state of a square lattice of 100 nitrogen atoms. This system’s 300-electron active space is out of reach for existing ab initio multireference methods. PiFCI+DFT accurately recovers the strongly correlated dissociation limit at a modest computational cost, motivating further development for modeling quantum materials.
The computational design of molecular quantum devices requires methods that capture "the quantum and the chemistry," approaching chemical accuracy for large numbers of entangled and/or strongly correlated electrons. Projected-interacting full configuration interaction (PiFCI) is a candidate for such simulations, providing a formally exact and systematically improvable approximation for correlation in large active spaces. PiFCI extends Kohn-Sham density functional theory by introducing multiple reference systems, each experiencing an electron-electron interaction projected onto one or more one-electron states. Compact CI expansions yield near-exact reference system wavefunctions, and projected exchange-correlation (XC) density functionals enable formally exact combinations of reference system correlation energies. This work presents a general treatment of the projected interactions in PiFCI and introduces regularized second-order many-body perturbation theory (MP2) as an approximate projected XC functional. Numerical results show that PiFCI plus regularized MP2 can accurately treat dynamical and nondynamical correlation in relatively large active spaces, including stacks of entangled singlet-coupled tetrathiafulvalene and phenalenyl organic radicals modeling molecular quantum devices.
The projection-based wave function in density functional theory (WF-in-DFT) embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [ Beran, P. J. Phys. Chem. Lett. 2023, 14, 716–722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H 20 chain and the cleavage of a triple CN bond in propionitrile.
We present a multilevel embedding framework for quantum chemistry calculations on a quantum computer. In our framework, a quantum algorithm treats the strongly correlated active space, while a high-level wave-function method such as coupled cluster theory or multireference perturbation theory recovers the remaining correlation in the surrounding region. A sampling-based quantum algorithm, quantum-selected configuration interaction, bridges the quantum and classical treatments. The entire calculation is embedded in a low-cost density functional theory description of the surrounding environment using Manby's projection technique. We apply the framework to organic, metal-organic, and metallic systems, computing bond dissociation energies, adsorption energies, and reaction barriers using only the subset of qubits of a 144-qubit superconducting quantum computer at the University of Osaka and achieving $\sim$1 kcal/mol agreement with classical references for a Menshutkin SN2 reaction inside a carbon nanotube. Our results may open the way to quantitatively reliable quantum-classical hybrid calculations for large-scale chemical systems.
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 density matrix renormalization group (DMRG) method has already proved itself as a very efficient and accurate computational method, which can treat large active spaces and capture the major part of strong correlation. Its application on larger molecules is, however, limited by its own computational scaling as well as demands of methods for treatment of the missing dynamical electron correlation. In this work, we present the first step in the direction of combining DMRG with density functional theory (DFT), one of the most employed quantum chemical methods with favorable scaling, by means of the projection-based wave function (WF)-in-DFT embedding. On two proof-of-concept but important molecular examples, we demonstrate that the developed DMRG-in-DFT approach provides a very accurate description of molecules with a strongly correlated fragment.
The projection-based wave function (WF)-in-DFT embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [Beran, P. et al. J. Phys. Chem. Lett. 2023, 14, 3, 716-722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H20 chain and the cleavage of a triple CN bond in propionitrile.