For decades, scientists have studied the intricate charge migration dynamics, where after ionization a localized charge distribution (‘hole’) migrates across the molecule on a femtosecond timescale. This has the potential for controlling electrons in molecules, yet a comprehensive understanding of the many aspects of charge migration is still missing. In this work, we analyze charge migration using an extension of localized intrinsic bond orbitals (IBOs). These orbitals lead to a compact representation of the dynamics and map the complex, correlated many-electron charge migration to chemical concepts such as curly arrows and orbital-orbital interactions. By analyzing multiple challenging scenarios, we show how IBOs enable us to identify key mechanisms in charge migration. For example, we show that different mechanisms are responsible for converting a π-shaped hole to a σ-shaped hole and vice versa. We explain these in terms of hyperconjugation interactions and configurations that couple orbitals with different symmetries. We further demonstrate how IBOs can be used to find molecules with high charge migration efficiency. We carry out all simulations using an efficient set up of the time-dependent density matrix renormalization group (TDDMRG), correlating as many as 45 electrons in 50 orbitals. We believe that our results will be useful to design future experiments. The proposed IBO analysis is applicable to other types of real-time electron dynamics and spectroscopy.

Tree tensor network states (TTNSs) combined with the density matrix renormalization group (DMRG) are emerging as powerful tools for vibrational and vibronic structure simulations in molecules with strong coupling and fluxionality. In this Perspective, we discuss how TTNS methods enable accurate, full-dimensional computations of thousands of eigenstates for molecular systems ranging from quartic-force-field benchmarks to molecules with strong vibronic coupling and protonated water clusters as large as the 33-dimensional Eigen ion, H3O+⋅(H2O)3. We emphasize the close connection and interoperability between DMRG-based TTNS methods and the multilayer multiconfiguration time-dependent Hartree method (ML-MCTDH), which share the same underlying ansatz. We also highlight practical challenges of predictive simulations, including robust error estimation, convergence of observables such as infrared intensities, and optimization of tensor network tree structures. Finally, we outline recent advances toward direct targeting of excited states and discuss opportunities for broader applications in molecular spectroscopy and quantum dynamics.

To understand the dynamics of quantum many-body systems, it is essential to study excited eigenstates. While tensor network states have become a standard tool for computing ground states in computational many-body physics, obtaining accurate excited eigenstates remains a significant challenge. In this work, we develop an approach that combines the inexact Lanczos method, which is designed for efficient computation of excited states, with tree tensor network states (TTNSs). We demonstrate our approach by computing excited vibrational states for three challenging problems: (1) 122 states in two different energy intervals of acetonitrile (12-dimensional), (2) Fermi resonance states of the fluxional Zundel ion (15-dimensional),and (3) selected excited states of the fluxional and very correlated Eigen ion (33-dimensional). The proposed TTNS inexact Lanczos method is directly applicable to other quantum many-body systems.

Accurate vibrational spectra are essential for understanding how molecules behave, yet their computation remains challenging and benchmark data to reliably compare different methods are sparse. Here, we present high-accuracy eigenstate computations for the six-atom, 12-dimensional acetonitrile molecule, a prototypical, strongly coupled, anharmonic system. Using a density matrix renormalization group (DMRG) algorithm with a tree-tensor-network-state (TTNS) ansatz, a refinement using TTNSs as basis set, and reliable procedures to estimate energy errors, we compute up to 5,000 vibrational states with error estimates below 0.0007 cm-1. Our analysis reveals that previous works underestimated the energy error by up to two orders of magnitude. Our data serve as a benchmark for future vibrational spectroscopy methods and our new method offers a path toward similarly precise computations of large, complex molecular systems.

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 non-singlet 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.

The nitrate radical NO3 plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic X~2A′2 state. Using recently developed tensor network state methods, we are able to compute more than 2500 vibronic states, thereby increasing the number of computed full-dimensional states by a factor of 50, compared to previous work. While we obtain good agreement with experiment for most of the assigned vibronic levels, for several others, we observe striking disagreement. Further, for the antisymmetric bending motion we find remarkably large symmetry-induced level splittings that are larger than the zero-order reference. We discuss non-negligible nonadiabatic effects and show that the Born-Oppenheimer approximation leads to significant errors in the spectrum.

The multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization 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 optimizing unoccupied single-particle functions in MCTDH, which corresponds to subspace enrichment in the DMRG. The second one is related to finding optimal tree structures and on highlighting similarities and differences of MCTDH and DMRG theories. We hope that this contribution will foster more fruitful cross-fertilization of ideas between ML-MCTDH and DMRG.

The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water dimer (Zundel ion) giving rise to the characteristic doublet peak is both more complex and more sensitive to subtle energetic changes than previously thought. In particular, hitherto overlooked low- intensity satellite peaks in the experimental spectrum are now unveiled and mechanistically assigned. Our findings rely on the comparison of IR spectra obtained using two highly accurate potential energy surfaces (PESs) in conjunction with highly accurate state-resolved quantum simulations. We demonstrate that these high-accuracy simulations are important for providing definite assignments of the complex IR signals of fluxional molecules.

The complex electronic structure and unusual potential energy curve of the chromium dimer have fascinated scientists for decades, with agreement between theory and experiment so far elusive. Here, we present a new ab initio simulation of the potential energy curve and vibrational spectrum that significantly improves on all earlier estimates. Our data support a shift in earlier experimental assignments of a cluster of vibrational frequencies by one quantum number. The new vibrational assignment yields an experimentally derived potential energy curve in quantitative agreement with theory across all bond lengths and across all measured frequencies. By solving this long-standing problem, our results raise the possibility of quantitative quantum chemical modeling of transition metal clusters with spectroscopic accuracy.

We explore various ways to group orbitals into clusters in a matrix product state (MPS). We explain how a generic cluster MPS can often lead to an increase in computational cost and instead propose a special cluster structure, involving only the first and last orbitals/sites, with a wider scope for computational advantage. This structure is a natural formalism to describe correlated multireference (MR) theories. We demonstrate the flexibility and usefulness of this approach by implementing various uncontracted MR configuration interaction, perturbation and linearized coupled cluster theories using an MPS with large cluster sites. Applications to the nitrogen dimer, the chromium dimer, and benzene, including up to triple excitations in the external space, demonstrate the utility of an MPS with up to two large sites. We use our results to analyze the quality of different multireference approximations.

Our article sheds light on the control of the motion of the ionized electrons during double ionization. To first order, the field only drives the electron into the same direction (front-to-back motion). In this paper, our goal was to achieve the opposite, back-to-back motion, where both electrons are ejected into opposite directions. This is only possible via an interplay between the inter-particle forces and the potential. By extensive usage of four different control procedures (local control, derivative-free optimizations of basis expansions, Krotov method, and control of the classical equations of motions), we could identify pulse characteristics that lead to the desired back-to-back motion. Additionally, we were able to identify a simple, classical two-step mechanism. Besides these new physical results, this contribution also discusses new methodologies for using local control procedures for this difficult system and employs a new method to efficiently propagate the wavefunction in time [H. R. Larsson, B. Hartke and D. J. Tannor, J. Chem. Phys., 145, 204108 (2016)].

During my postdoctoral studies in the group of Prof. G. Chan (Caltech), I explored new ways to apply the density matrix renormalization group (DMRG) to bigger systems. In this article we explore a class of simple, qualitative wavefunctions, the minimal matrix product state (MMPS), which combines projection operators and the DMRG in order to yield a low-scaling method for quick explorations of the potential energy landscape of strongly correlated systems. Without jeopardizing low computational cost, MMPSs generalize well-established methods such as the projected Hartree-Fock-Bogoliubov (HFB)/antisymmetrized geminal power (AGP) method and the generalized valence bond (GVB) method. For the systems we studied, we show that MMPSs give correct qualitative behavior across the whole PES, often significantly improving aforementioned ansätze.

During my postdoctoral studies in the group of Prof. G. Chan (Caltech), as side project, I worked on the computation of vibrational spectra using tree-tensor network states (TTNS). There, I applied methods developed in the condensed matter community to molecular systems. Compared to established approaches for computing vibrational spectra, such as the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method, the new approach is faster and much more robust and the diagrammatic language from condensed matter physics I use gains much more insights. I am currently working on using the TTNS method to compute high-lying vibrational states of the Zundel ion.