In this work, we derive a multi-fragment real-time extension of the projected density matrix embedding theory (pDMET) designed to treat non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static pDMET, the real time pDMET partitions the total system into many fragments; the coupling between each fragment and the rest of the system is treated through a compact representation of the environment in terms of a quantum bath. The real-time pDMET involves simultaneously propagating the wavefunctions for each separate fragment-bath embedding system along with an auxiliary mean-field wavefunction of the total system. The equations of motion are derived by (i) projecting the time-dependent Schrödinger equation in the fragment and bath space associated with each separate fragment and by (ii) enforcing the pDMET matching conditions between the global 1-particle reduced density matrix (1-RDM) obtained from the fragment calculations and the mean-field 1-RDM at all points in time. The accuracy of the method is benchmarked through comparisons to time-dependent density-matrix renormalization group and time-dependent Hartree-Fock (TDHF) theory; the methods were applied to a one- and two-dimensional single-impurity Anderson model and multi-impurity Anderson models with ordered and disordered distributions of the impurities. The results demonstrate a large improvement over TDHF and rapid convergence to the exact dynamics with an increase in fragment size. Our results demonstrate that the real-time pDMET is a promising and flexible method that balances accuracy and efficiency to simulate the non-equilibrium electron dynamics in heterogeneous systems of large size.
In this work, we derive a multi-fragment real-time extension of the projected density matrix embedding theory (pDMET) designed to treat non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static pDMET, the real time pDMET partitions the total system into many fragments; the coupling between each fragment and the rest of the system is treated through a compact representation of the environment in terms of a quantum bath. The real-time pDMET involves simultaneously propagating the wavefunctions for each separate fragment-bath embedding system along with an auxiliary mean-field wavefunction of the total system. The equations of motion are derived by (i) projecting the time-dependent Schrödinger equation in the fragment and bath space associated with each separate fragment and by (ii) enforcing the pDMET matching conditions between the global 1-particle reduced density matrix (1-RDM) obtained from the fragment calculations and the mean-field 1-RDM at all points in time. The accuracy of the method is benchmarked through comparisons to time-dependent density-matrix renormalization group and time-dependent Hartree-Fock (TDHF) theory; the methods were applied to a one- and two-dimensional single-impurity Anderson model and multi-impurity Anderson models with ordered and disordered distributions of the impurities. The results demonstrate a large improvement over TDHF and rapid convergence to the exact dynamics with an increase in fragment size. Our results demonstrate that the real-time pDMET is a promising and flexible method that balances accuracy and efficiency to simulate the non-equilibrium electron dynamics in heterogeneous systems of large size.
We introduce real-time density matrix embedding theory (DMET), a dynamical quantum embedding theory for computing non-equilibrium electron dynamics in strongly correlated systems. As in the previously developed static DMET, real-time DMET partitions the system into an impurity corresponding to the region of interest coupled to the surrounding environment, which is efficiently represented by a quantum bath of the same size as the impurity. In this work, we focus on a simplified single-impurity time-dependent formulation as a first step toward a multi-impurity theory. The equations of motion of the coupled impurity and bath embedding problem are derived using the time-dependent variational principle. The accuracy of real-time DMET is compared to that of time-dependent complete active space self-consistent field (TD-CASSCF) theory and time-dependent Hartree-Fock (TDHF) theory for a variety of quantum quenches in the single impurity Anderson model (SIAM), in which the Hamiltonian is suddenly changed (quenched) to induce a non-equilibrium state. Real-time DMET shows a marked improvement over the mean-field TDHF, converging to the exact answer even in the non-trivial Kondo regime of the SIAM. However, as expected from analogous behavior in static DMET, the constrained structure of the real-time DMET wavefunction leads to a slower convergence with respect to active space size, in the single-impurity formulation, relative to TD-CASSCF. Our initial results suggest that real-time DMET provides a promising framework to simulate non-equilibrium electron dynamics in which strong electron correlation plays an important role, and lays the groundwork for future multi-impurity formulations.
The simulation of non-equil. electron dynamics in real systems provides a challenge for theor. methods due to the need to treat both large system sizes and electron correlation. Towards this goal, we present an extension of the d. matrix embedding theory (DMET) for the simulation of real-time electron dynamics in strongly correlated systems. As in the previously developed static DMET, real-time DMET partitions the system into an impurity corresponding the region of interest coupled to the surrounding environment, which is efficiently represented by a quantum bath of the same size as the impurity. The dynamics of the coupled impurity and bath embedding problem are obtained through use of the time-dependent variational principle. The methodol. allows for the efficient and accurate simulation of non-equil. electron dynamics in the presence of strong correlation, reaching total system sizes unobtainable by conventional methodol.
In recent studies by Yalouz et al. [J. Chem. Phys. 157, 214112 (2022)] and Sekaran et al. [Phys. Rev. B 104, 035121 (2021) and Computation 10, 45 (2022)], density matrix embedding theory (DMET) has been reformulated through the use of the Householder transformation as a novel tool to embed a fragment within extended systems. The transformation was applied to a reference non-interacting one-electron reduced density matrix to construct fragments' bath orbitals, which are crucial for subsequent ground state calculations. In the present work, we expand upon these previous developments and extend the utilization of the Householder transformation to the description of multiple electronic states, including ground and excited states. Based on an ensemble noninteracting density matrix, we demonstrate the feasibility of achieving exact fragment embedding through successive Householder transformations, resulting in a larger set of bath orbitals. We analytically prove that the number of additional bath orbitals scales directly with the number of fractionally occupied natural orbitals in the reference ensemble density matrix. A connection with the regular DMET bath construction is also made. Then, we illustrate the use of this ensemble embedding tool in single-shot DMET calculations to describe both ground and first excited states in a Hubbard lattice model and an ab initio hydrogen system. Finally, we discuss avenues for enhancing ensemble embedding through self-consistency and explore potential future directions.
Quantum embedding methods enable the study of large, strongly correlated quantum systems by (usually self-consistent) decomposition into computationally manageable subproblems, in the spirit of divide-and-conquer methods. Among these, Density Matrix Embedding Theory (DMET) is an efficient approach that enforces self-consistency at the level of one-particle reduced density matrices (1-RDMs), facilitating applications across diverse quantum systems. However, conventional DMET is constrained by the requirement that the global 1-RDM (low-level descriptor) be an orthogonal projector, limiting flexibility in bath construction and potentially impeding accuracy in strongly correlated regimes. In this work, we introduce a generalized DMET framework in which the low-level descriptor can be an arbitrary 1-RDM and the bath construction is based on optimizing a quantitative criterion related to the maximal disentanglement between different fragments. This yields an alternative yet controllable bath space construction for generic 1-RDMs, lifting a key limitation of conventional DMET. We demonstrate its consistency with conventional DMET in appropriate limits and exploring its implications for bath construction, downfolding (impurity Hamiltonian construction), low-level solvers, and adaptive fragmentation. We expect that this more flexible framework, which leads to several new variants of DMET, can improve the robustness and accuracy of DMET.
We extend our density matrix embedding theory (DMET) [Phys. Rev. Lett.2012, 109, 186404] from lattice models to the full chemical Hamiltonian. DMET allows the many-body embedding of arbitrary fragments of a quantum system, even when such fragments are open systems and strongly coupled to their environment (e.g., by covalent bonds). In DMET, empirical approaches to strong coupling, such as link atoms or boundary regions, are replaced by a small, rigorous quantum bath designed to reproduce the entanglement between a fragment and its environment. We describe the theory and demonstrate its feasibility in strongly correlated hydrogen ring and grid models; these are not only beyond the scope of traditional embeddings but even challenge conventional quantum chemistry methods themselves. We find that DMET correctly describes the notoriously difficult symmetric dissociation of a 4 × 3 hydrogen atom grid, even when the treated fragments are as small as single hydrogen atoms. We expect that DMET will open up new ways of treating complex strongly coupled, strongly correlated systems in terms of their individual fragments.
Density matrix embedding theory (DMET) [Phys. Rev. Lett. 109, 186404 (2012)], introduced an approach to quantum cluster embedding methods whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost but was inherently limited by the construction of a bath designed to reproduce ground-state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency-dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the one- (1D) and two-dimensional (2D) Hubbard model, where we obtain zero temperature and thermodynamic limit spectral functions, and show the trivial extension to two-particle Green's functions. This advance therefore extends the scope and applicability of DMET in condensed-matter problems as a computationally tractable route to correlated spectral functions of extended systems and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.