Time-dependent Density Matrix Renormalization Group Research Articles

Simulating Real-Time Molecular Electron Dynamics Efficiently Using the Time-Dependent Density Matrix Renormalization Group.

Compared to ground-state electronic structure optimizations, accurate simulations of molecular real-time electron dynamics are usually much more difficult to perform. To simulate electron dynamics, the time-dependent density matrix renormalization group (TDDMRG) has been shown to offer an attractive compromise between accuracy and cost. However, many simulation parameters significantly affect the quality and efficiency of a TDDMRG simulation. So far, it is unclear whether common wisdom from ground-state DMRG carries over to the TDDMRG, and a guideline on how to choose these parameters is missing. Here, in order to establish such a guideline, we investigate the convergence behavior of the main TDDMRG simulation parameters, such as time integrator, the choice of orbitals, and the choice of matrix-product-state representation for complex-valued nonsinglet states. In addition, we propose a method to select orbitals that are tailored to optimize the dynamics. Lastly, we showcase the TDDMRG by applying it to charge migration ionization dynamics in furfural, where we reveal a rapid conversion from an ionized state with a σ character to one with a π character within less than a femtosecond.

TD-DMRG Study of Exciton Dynamics with both Thermal and Static Disorders for Fenna-Matthews-Olson Complex.

Photosynthesis is a fundamental process that converts solar energy into chemical energy. Understanding the microscopic mechanisms of energy transfer in photosynthetic systems is crucial for the development of novel optoelectronic materials. Simulating these processes poses significant challenges due to the intricate interactions between electrons and phonons, compounded by static disorder. In this work, we present a numerically nearly exact study using the time-dependent density matrix renormalization group (TD-DMRG) method to simulate the quantum dynamics of the Fenna-Matthews-Olson (FMO) complex considering an eight-site model with both thermal and static disorders. We employ the thermo-field dynamics formalism for temperature effects. We merge all electronic interactions into one large matrix product state (MPS) site, boosting accuracy efficiently without increasing complexity. Previous combined experimental and computational studies indicated that the static disorders range from 30 to 90 cm-1 for different FMO sites. We employ a Gaussian distribution and the auxiliary bosonic operator approach to consider the static disorder in our TD-DMRG algorithm. We investigate the impact of different initial excitation sites, temperatures, and degrees of static disorder on the exciton dynamics and temporal coherence. It is found that under the influence of the experimentally determined static disorder strength, the exciton population evolution shows a non-negligible difference at zero temperature, while it is hardly affected at room temperature.

Kylin-V: An open-source package calculating the dynamic and spectroscopic properties of large systems.

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples.

Theoretical investigation of distal charge separation in a perylenediimide trimer

An exciton–phonon (ex–ph) model based on our recently developed block interaction product basis framework is introduced to simulate the distal charge separation (CS) process in aggregated perylenediimide (PDI) trimer incorporating the quantum dynamic method, i.e., the time-dependent density matrix renormalization group. The electronic Hamiltonian in the ex–ph model is represented by nine constructed diabatic states, which include three local excited (LE) states and six charge transfer (CT) states from both the neighboring and distal chromophores. These diabatic states are automatically generated from the direct products of the leading localized neutral or ionic states of each chromophore’s reduced density matrix, which are obtained from ab initio quantum chemical calculation of the subsystem consisting of the targeted chromophore and its nearest neighbors, thus considering the interaction of the adjacent environment. In order to quantum-dynamically simulate the distal CS process with massive coupled vibrational modes in molecular aggregates, we used our recently proposed hierarchical mapping approach to renormalize these modes and truncate those vibrational modes that are not effectively coupled with electronic states accordingly. The simulation result demonstrates that the formation of the distal CS process undergoes an intermediate state of adjacent CT, i.e., starts from the LE states, passes through an adjacent CT state to generate the intermediates (∼200 fs), and then formalizes the targeted distal CS via further charge transference (∼1 ps). This finding agrees well with the results observed in the experiment, indicating that our scheme is capable of quantitatively investigating the CS process in a realistic aggregated PDI trimer and can also be potentially applied to exploring CS and other photoinduced processes in larger systems.

Using Kondo entanglement to induce spin correlations between disconnected quantum dots

We investigate the entanglement between the spins of two quantum dots that are not simultaneously connected to the same system. Quantum entanglement among localized spins is a crucial property for the advancement of quantum computing and quantum information. Generating and controlling an entangled state between quantum dots have garnered significant attention in recent years for this reason. In this study, we demonstrate that information about the spin orientation of a quantum dot can be preserved, utilizing Kondo entanglement, within a reservoir of electrons. Subsequently, this information can be transmitted to another dot after the initial dot has been decoupled from the reservoirs. We employ a double quantum dot system in a parallel geometry to establish the initial state, where each dot interacts with reservoirs of different symmetries. A specific phase in the couplings is chosen to induce antiferromagnetic spin correlation between the dots. The time evolution of the initial state is analyzed using the time-dependent density matrix renormalization group method. Our findings reveal that a partially entangled state between the dots can be achieved, even when they are not simultaneously connected. This entangled state arises transiently and dissipates in the stationary state. The stability of the state observed during the transient phase is demonstrated. To comprehend the details of these phenomena, we employ a canonical transformation of real space.

Keldysh crossover in one-dimensional Mott insulators

Recent advancements in pulse laser technology have facilitated the exploration of nonequilibrium spectroscopy of electronic states in the presence of strong electric fields across a broad range of photon energies. The Keldysh crossover serves as an indicator that distinguishes between excitations resulting from photon absorption triggered by near-infrared multicycle pulses and those arising from quantum tunneling induced by terahertz pulses. Using a time-dependent density-matrix renormalization group, we investigate the emergence of the Keldysh crossover in a one-dimensional (1D) Mott insulator. We find that the Drude weight is proportional to photo-doped doublon density when a pump pulse induces photon absorption. In contrast, the Drude weight is suppressed when a terahertz pulse introduces doublons and holons via quantum tunneling. The suppressed Drude weight accompanies glassy dynamics with suppressed diffusion, which is a consequence of strong correlations and exhibits finite polarization decaying slowly after pulse irradiation. In the quantum tunneling region, entanglement entropy slowly grows logarithmically. These contrasting behaviors between the photon-absorption and quantum tunneling regions are a manifestation of the Keldysh crossover in 1D Mott insulators and provide a novel methodology for designing the localization and symmetry of electronic states called subcycle-pulse engineering.

A tensor network view of multilayer multiconfiguration time-dependent Hartree methods

The multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalisation group (DMRG) are powerful workhorses applied mostly in different scientific fields. Although both methods are based on tensor network states, very different mathematical languages are used for describing them. This severely limits knowledge transfer and sometimes leads to re-inventions of ideas well known in the other field. Here, we review ML-MCTDH and DMRG theory using both MCTDH expressions and tensor network diagrams. We derive the ML-MCTDH equations of motions using diagrams and compare them with time-dependent and time-independent DMRG algorithms. We further review two selected recent advancements. The first advancement is related to optimising unoccupied single-particle functions in MCTDH, which corresponds to subspace enrichment in the DMRG. The second advancement is related to finding optimal tree structures and on highlighting similarities and differences of tensor networks used in MCTDH and DMRG theories. We hope that this contribution will foster more fruitful cross-fertilisation of ideas between ML-MCTDH and DMRG.

Divergent Nonlinear Response from Quasiparticle Interactions.

We demonstrate that nonlinear response functions in many-body systems carry a sharp signature of interactions between gapped low-energy quasiparticles. Such interactions are challenging to deduce from linear response measurements. The signature takes the form of a divergent-in-time contribution to the response-linear in time in the case when quasiparticles propagate ballistically-that is absent for free bosonic excitations. We give a physically transparent semiclassical picture of this singular behavior. While the semiclassical picture applies to a broad class of systems we benchmark it in two simple models: in the Ising chain using a form-factor expansion, and in a nonintegrable model-the spin-1 Affleck-Kennedy-Lieb-Tasaki chain-using time-dependent density matrix renormalization group simulations. We comment on extensions of these results to finite temperatures.

Auger spectroscopy beyond the ultra-short core-hole relaxation time approximation

We present a time-dependent computational approach to study Auger electron spectroscopy (AES) beyond the ultra-short core-hole relaxation time approximation and, as a test case, we apply it to the paradigmatic example of a one-dimensional Mott insulator represented by a half-filled Hubbard chain. The Auger spectrum is usually calculated by assuming that, after the creation of a core-hole, the system thermalizes almost instantaneously. This leads to a relatively simple analytical expression that uses the ground-state with a core-hole as a reference state and ignores all the transient dynamics related to the screening of the core-hole. In this picture, the response of the system can be associated to the pair spectral function. On the other hand, in our numerical calculations, the core hole is created by a light pulse, allowing one to study the transient dynamics of the system in terms of the pulse duration and in the non-perturbative regime. Time-dependent density matrix renormalization group calculations reveal that the relaxation process involves the creation of a polarization cloud of doublon excitations that have an effect similar to photo-doping. As a consequence, there is a leak of spectral weight to higher energies into what otherwise would be the Mott gap. For longer pulses, these excited states, mostly comprised of doublons, can dominate the spectrum. By changing the duration of the light-pulse, the entire screening process can be resolved in time.

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

Time-Dependent Density Matrix Renormalization Group Method for Quantum Transport with Phonon Coupling in Molecular Junction.

Quantum transport in molecular junctions has attracted great attention. The charge motion in a molecular junction can cause geometric deformation, leading to strong electron phonon coupling, which was often overlooked. We have formulated a nearly exact method to assess the time-dependent current and occupation number in the molecular junction modeled by the electron-phonon coupled bridge state using the time-dependent density matrix renormalization group (TD-DMRG) method. The oscillation period and amplitude of the current are found to be dependent on the electron phonon coupling strength and energy level alignment with the electrodes. In an attempt to better understand these phenomena, we have devised a new approximation that explains the bistability phenomenon and the behavior of steady currents in the strong electron-phonon coupling regime. Comparisons have been made with the multilayer-multiconfiguration time-dependent Hartree (ML-MCTDH) method and the analytical result in the purely electronic limit. Furthermore, we explore the entropy of different orderings, extending to the electron phonon model problems. Regarding finite temperature, the thermal Bogoliubov transformation of both fermions and bosons is used and compared with imaginary time evolution results.

Minimizing non-radiative decay in molecular aggregates through control of excitonic coupling

The widely known “Energy Gap Law” (EGL) predicts a monotonically exponential increase in the non-radiative decay rate (knr) as the energy gap narrows, which hinders the development of near-infrared (NIR) emissive molecular materials. Recently, several experiments proposed that the exciton delocalization in molecular aggregates could counteract EGL to facilitate NIR emission. In this work, the nearly exact time-dependent density matrix renormalization group (TD-DMRG) method is developed to evaluate the non-radiative decay rate for exciton-phonon coupled molecular aggregates. Systematical numerical simulations show, by increasing the excitonic coupling, knr will first decrease, then reach a minimum, and finally start to increase to follow EGL, which is an overall result of two opposite effects of a smaller energy gap and a smaller effective electron-phonon coupling. This anomalous non-monotonic behavior is found robust in a number of models, including dimer, one-dimensional chain, and two-dimensional square lattice. The optimal excitonic coupling strength that gives the minimum knr is about half of the monomer reorganization energy and is also influenced by system size, dimensionality, and temperature.

Traces of electron-phonon coupling in one-dimensional cuprates

The appearance of certain spectral features in one-dimensional (1D) cuprate materials has been attributed to a strong, extended attractive coupling between electrons. Here, using time-dependent density matrix renormalization group methods on a Hubbard-extended Holstein model, we show that extended electron-phonon (e–ph) coupling presents an obvious choice to produce such an attractive interaction that reproduces the observed spectral features and doping dependence seen in angle-resolved photoemission experiments: diminished 3kF spectral weight, prominent spectral intensity of a holon-folding branch, and the correct holon band width. While extended e–ph coupling does not qualitatively alter the ground state of the 1D system compared to the Hubbard model, it quantitatively enhances the long-range superconducting correlations and suppresses spin correlations. Such an extended e–ph interaction may be an important missing ingredient in describing the physics of the structurally similar two-dimensional high-temperature superconducting layered cuprates, which may tip the balance between intertwined orders in favor of uniform d-wave superconductivity.

Density-Matrix Renormalization Group Algorithm for Simulating Quantum Circuits with a Finite Fidelity

We develop a density-matrix renormalization group (DMRG) algorithm for the simulation of quantum circuits. This algorithm can be seen as the extension of time-dependent DMRG from the usual situation of hermitian Hamiltonian matrices to quantum circuits defined by unitary matrices. For small circuit depths, the technique is exact and equivalent to other matrix product state (MPS) based techniques. For larger depths, it becomes approximate in exchange for an exponential speed up in computational time. Like an actual quantum computer, the quality of the DMRG results is characterized by a finite fidelity. However, unlike a quantum computer, the fidelity depends strongly on the quantum circuit considered. For the most difficult possible circuit for this technique, the so-called "quantum supremacy" benchmark of Google Inc. , we find that the DMRG algorithm can generate bit strings of the same quality as the seminal Google experiment on a single computing core. For a more structured circuit used for combinatorial optimization (Quantum Approximate Optimization Algorithm or QAOA), we find a drastic improvement of the DMRG results with error rates dropping by a factor of 100 compared with random quantum circuits. Our results suggest that the current bottleneck of quantum computers is their fidelities rather than the number of qubits.

A multi-fragment real-time extension of projected density matrix embedding theory: Non-equilibrium electron dynamics in extended systems.

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.

Hierarchical Mapping for Efficient Simulation of Strong System-Environment Interactions.

Quantum dynamics (QD) simulation is a powerful tool for interpreting ultrafast spectroscopy experiments and unraveling their microscopic mechanism in out-of-equilibrium excited state behaviors in various chemical, biological, and material systems. Although state-of-the-art numerical QD approaches such as the time-dependent density matrix renormalization group (TD-DMRG) already greatly extended the solvable system size of general linearly coupled exciton-phonon models with up to a few hundred phonon modes, the accurate simulation of larger system sizes or strong system-environment interactions is still computationally highly challenging. Based on quantum information theory (QIT), in this work, we realize that only a small number of effective phonon modes couple to the excitonic system directly regardless of a large or even infinite number of modes in the condensed phase environment. On top of the identified small number of direct effective modes, we propose a hierarchical mapping (HM) approach through performing block Lanczos transformations on the remaining indirect modes, which transforms the Hamiltonian matrix to a nearly block-tridiagonal form and eliminates the long-range interactions. Numerical tests on model spin-boson systems and realistic singlet fission models in a rubrene crystal environment with up to 7000 modes and strong system-environment interactions indicate HM can reduce the system size by 1-2 orders of magnitude and accelerate the calculation by ∼80% without losing accuracy.

Tuning the Topological θ -Angle in Cold-Atom Quantum Simulators of Gauge Theories

The topological θ-angle in gauge theories engenders a series of fundamental phenomena, including violations of charge-parity (CP) symmetry, dynamical topological transitions, and confinement-deconfinement transitions. At the same time, it poses major challenges for theoretical studies, as it implies a sign problem in numerical simulations. Analog quantum simulators open the promising prospect of treating quantum many-body systems with such topological terms, but, contrary to their digital counterparts, they have not yet demonstrated the capacity to control the θ-angle. Here, we demonstrate how a tunable topological θ-term can be added to a prototype theory with U(1) gauge symmetry, a discretized version of quantum electrodynamics in one spatial dimension. As we show, the model can be realized experimentally in a single-species Bose-Hubbard model in an optical superlattice with three different spatial periods, thus requiring only standard experimental resources. Through numerical calculations obtained from the time-dependent density-matrix renormalization group method and exact diagonalization, we benchmark the model system, and illustrate how salient effects due to the θ-term can be observed. These include charge confinement, the weakening of quantum many-body scarring, as well as the disappearance of Coleman’s phase transition due to explicit breaking of CP symmetry. This work opens the door towards studying the rich physics of topological gauge-theory terms in large-scale cold-atom quantum simulators.4 MoreReceived 15 April 2022Revised 19 September 2022Accepted 26 September 2022DOI:https://doi.org/10.1103/PRXQuantum.3.040316Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasColor confinementColor deconfinementConfinementQuantum simulationPhysical Systems1-dimensional systemsBose gasesTechniquesBose-Hubbard modelDensity matrix renormalization groupExact diagonalizationGauge symmetriesLattice gauge theoryPerturbation theoryAtomic, Molecular & OpticalCondensed Matter, Materials & Applied Physics

Formation of spin and charge ordering in the extended Hubbard model during a finite-time quantum quench

We investigate the formation of charge and spin ordering by starting from a non-interacting state and studying how it evolves in time under a Hamiltonian with finite electronic interactions. We consider the one-dimensional, half-filled extended Hubbard model, which we solve within time-dependent density matrix renormalization group. By employing linear finite-time quenches in the onsite and nearest-neighbor interactions, we find the existence of impulse, intermediate, and adiabatic regimes of time evolution. For the quenches we analyze, we observe that the adiabatic regime is reached with distinct ramping time scales depending on whether the charge density wave (CDW) or the spin density wave (SDW) is formed. The former needs to be slower than the latter to prevent entangled excited states from being accessed during the quench. More interestingly, in the intermediate regime, we observe an enhancement of the entanglement entropy with respect to its initial value, which precedes the formation of the CDW ordering; a similar enhancement is not seen in the quench towards SDW. Our findings also show that the breaking of the system integrability, by turning on the nearest-neighbor interactions, does not give rise to significant changes in the non-equilibrium behavior within the adiabatic approximation.

Computational Method for Evaluating the Thermoelectric Power Factor for Organic Materials Modeled by the Holstein Model: A Time-Dependent Density Matrix Renormalization Group Formalism.

Organic/polymeric materials are of emerging importance for thermoelectric conversion. The soft nature of these materials implies strong electron-phonon coupling, often leading to carrier localization. This poses great challenges for the conventional Boltzmann transport description based on relaxation time approximation and band structure calculations. In this work, combining the Kubo formula with the finite-temperature time-dependent density matrix renormalization group (FT-TD-DMRG) in the grand canonical ensemble, we developed a nearly exact algorithm to calculate the thermoelectric power factor PF = α2 σ, where α is the Seebeck coefficient and σ is the electrical conductivity, and apply the algorithm to Holstein Hamiltonian with electron-phonon coupling to model organic materials. Our algorithm can provide a unified description covering the weak coupling limit described by the bandlike Boltzmann transport to the strong coupling hopping limit.

Glassy dynamics of the one-dimensional Mott insulator excited by a strong terahertz pulse

The elucidation of nonequilibrium states in strongly correlated systems holds the key to emergence of novel quantum phases. The nonequilibrium-induced insulator-to-metal transition is particularly interesting since it reflects the fundamental nature of competition between itinerancy and localization of the charge degrees of freedom. We investigate pulse-excited insulator-to-metal transition of the half-filled one-dimensional extended Hubbard model. Calculating the time-dependent optical conductivity with the time-dependent density-matrix renormalization group, we find that strong mono- and half-cycle pulses inducing quantum tunneling strongly suppress spectral weights contributing to the Drude weight σD, even if we introduce a large number of carriers Δnd. This is in contrast to a metallic behavior of σD∝Δnd induced by photon absorption and chemical doping. The strong suppression of σD in quantum tunneling is a result of the emergence of the Hilbert-space fragmentation, which makes pulse-excited states glassy.Received 9 February 2022Revised 15 May 2022Accepted 18 July 2022DOI:https://doi.org/10.1103/PhysRevResearch.4.L032019Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasLocalizationNonequilibrium statistical mechanicsOptical conductivityPhotoinduced effectQuantum statistical mechanicsPhysical Systems1-dimensional systemsGlassy systemsMott insulatorsCondensed Matter, Materials & Applied PhysicsStatistical PhysicsQuantum Information