Multiconfiguration Time-dependent Hartree Research Articles

Temperature Dependence of the Electronic Absorption Spectrum of NO2.

The nitrogen dioxide (NO2) radical is composed of the two most abundant elements in the atmosphere, where it can be formed in a variety of ways including combustion, detonation of energetic materials, and lightning. Relevant also to smog and ozone cycles, together these processes span a wide range of temperatures. Remarkably, high-resolution NO2 electronic absorption spectra have only been reported in a narrow range below about 300 K. Previously, we reported [ J. Phys. Chem. A 2021, 125, 5519-5533] the construction of quasi-diabatic potential energy surfaces (PESs) for the lowest four electronic states (X̃, Ã, B̃, and C̃) of NO2. In addition to three-dimensional PESs based on explicitly correlated MRCI(Q)-F12/VTZ-F12 ab initio data, the geometry dependence of each component of the dipoles and transition dipoles was also mapped into fitted surfaces. The multiconfigurational time-dependent Hartree (MCTDH) method was then used to compute the 0 K electronic absorption spectrum (from the ground rovibrational initial state) employing those energy and transition dipole surfaces. Here, in an extension of that work, we report an investigation into the effects of elevated temperature on the spectrum, considering the effects of the population of rotationally and vibrationally excited initial states. The calculations are complemented by new experimental measurements. Spectral contributions from hundreds of rotational states up to N = 20 and from 200 individually-characterized vibrational states were computed. A spectral simulation tool was developed that enables modeling the spectrum at various temperatures─by weighting individual spectral contributions via the partition function, or for pure excited initial states, which can be probed via transient absorption spectroscopy. We validate these results against experimental absorption spectroscopy data at high temperatures, as well as via a new measurement from the (1,0,1) initial vibrational state.

A numerically exact correlation discrete variable representation for multi-configurational time-dependent Hartree calculations.

The correlation discrete variable representation (CDVR) enables (multilayer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. The CDVR employs a set of grids corresponding to single-particle functions to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. In standard CDVR approaches, the number of grid points employed is tied to the number of corresponding single-particle functions. This limits the accuracy of the quadrature, which can be achieved for a given single-particle function basis. In this work, an extended CDVR approach that facilitates a numerically exact quadrature of all potential matrix elements is introduced. The number of grid points employed can be increased independent of the number of corresponding single-particle function to achieve any desired quadrature accuracy. The properties of the new scheme are illustrated by numerical calculations studying the photodissociation of NOCl and the vibrational states of CH3. Fast convergence with respect to the number of additional quadrature points is observed: Employing a single additional point in each physical or logical coordinate already ensures negligible quadrature errors.

Exploring disordered quantum spin models with a multilayer multiconfigurational approach

Numerical simulations of quantum spin models are crucial for a profound understanding of many-body phenomena in a variety of research areas in physics. An outstanding problem is the availability of methods to tackle systems that violate area-laws of entanglement entropy. Such scenarios cover a wide range of compelling physical situations including disordered quantum spin systems among others. In this work, we employ a numerical technique referred to as multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) to evaluate the ground state of several disordered spin models. ML-MCTDH has previously been used to study problems of high-dimensional quantum dynamics in molecular and ultracold physics but is here applied to study spin systems for the first time. We exploit the inherent flexibility of the method to present results in one and two spatial dimensions and treat challenging setups that incorporate long-range interactions as well as disorder. Our results suggest that the hierarchical multi-layering inherent to ML-MCTDH allows to tackle a wide range of quantum many-body problems such as spin dynamics of varying dimensionality.

Efficient simulation of multielectron dynamics in molecules under intense laser pulses: implementation of the multiconfiguration time-dependent Hartree–Fock method based on the adaptive finite element method

We present an implementation of the multiconfiguration time-dependent Hartree–Fock method based on the adaptive finite element method for molecules under intense laser pulses. For efficient simulations, orbital functions are propagated by a stable propagator using the short iterative Arnoldi scheme and our implementation is parallelized for distributed memory computing. This is demonstrated by simulating high-harmonic generation from a water molecule and achieves a simulation of multielectron dynamics with overwhelmingly less computational time, compared to our previous work [Sawada, R. et al. Phys. Rev. A, 2016, 93, 023434].

Confining CO2 inside sI clathrate‐hydrates: The impact of the CO2–water interaction on quantized dynamics

We report new results on the translational-rotational (T-R) states of the CO2 molecule inside the sI clathrate-hydrate cages. We adopted the multiconfiguration time-dependent Hartree methodology to solve the nuclear molecular Hamiltonian, and to address issues on the T-R couplings. Motivated by experimental X-ray observations on the CO2 orientation in the D and T sI cages, we aim to evaluate the effect of the CO2 -water interaction on quantum dynamics. Thus, we first compared semiempirical and ab initio-based pair interaction model potentials against first-principles DFT-D calculations for ascertaining the importance of nonadditive many-body effects on such guest-host interactions. Our results reveal that the rotational and translational excited states quantum dynamics is remarkably different, with the pattern and density of states clearly affected by the underlying potential model. By analyzing the corresponding the probability density distributions of the calculated T-R eigenstates on both semiempirical and ab initio pair CO2 -water nanocage potentials, we have extracted information on the altered CO2 guest local structure, and we discussed it in connection with experimental data on the orientation of the CO2 molecule in the D and T sI clathrate cages available from neutron diffraction and solid-state NMR studies, as well as in comparison with previous molecular dynamics simulations. Our calculations provide a very sensitive test of the potential quality by predicting the low-lying T-R states and corresponding transitions for the encapsulated CO2 molecule. As such spectroscopic observables have not been measured so far, our results could trigger further detailed experimental and theoretical investigations leading to a quantitative description of the present guest-host interactions.

Nonlinear Molecular Electronic Spectroscopy via MCTDH Quantum Dynamics: From Exact to Approximate Expressions.

We present an accurate and efficient approach to computing the linear and nonlinear optical spectroscopy of a closed quantum system subject to impulsive interactions with an incident electromagnetic field. It incorporates the effect of ultrafast nonadiabatic dynamics by means of explicit numerical propagation of the nuclear wave packet. The fundamental expressions for the evaluation of first- and higher-order response functions are recast in a general form that can be used with any quantum dynamics code capable of computing the overlap of nuclear wave packets evolving in different states. Here we present the evaluation of these expressions with the multiconfiguration time-dependent Hartree (MCTDH) method. Application is made to pyrene, excited to its lowest bright excited state S2 which exhibits a sub-100-fs nonadiabatic decay to a dark state S1. The system is described by a linear vibronic coupling Hamiltonian, parametrized with multiconfiguration electronic structure methods. We show that the ultrafast nonadiabatic dynamics can have a remarkable effect on the spectral line shapes that goes beyond simple lifetime broadening. Furthermore, a widely employed approximate expression based on the time scale separation of dephasing and population relaxation is recast in the same theoretical framework. Application to pyrene shows the range of validity of such approximations.

Charged ultralong-range Rydberg trimers

We show that the recently observed class of long-range ion-Rydberg molecules can be divided into two families of states, which are characterised by their unique electronic structures resulting from the ion-induced admixture of quantum defect-split Rydberg $n$P states with different low-field seeking high-$l$ states. We predict that in both cases these diatomic molecular states can bind additional ground state atoms lying within the orbit of the Rydberg electron, thereby forming charged ultralong-range Rydberg molecules (ULRM) with binding energies similar to that of conventional non-polar ULRM. To demonstrate this, we consider a Rydberg atom interacting with a single ground state atom and an ion. The additional atom breaks the system's cylindrical symmetry, which leads to mixing between states that would otherwise be decoupled. The electronic structure is obtained using exact diagonalisation over a finite basis and the vibrational structure is determined using the Multi-Configuration Time-Dependent Hartree method. Due to the lobe-like structure of the electronic density, bound trimers with both linear and nonlinear geometrical configurations of the three nuclei are possible. The predicted trimer binding energies and excitation series are distinct enough from those of the ion-Rydberg dimer to be observed using current experimental techniques.

Optimal Mode Combination in the Multiconfiguration Time-Dependent Hartree Method through Multivariate Statistics: Factor Analysis and Hierarchical Clustering.

The multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extension (ML-MCTDH) are powerful algorithms for the efficient computation of nuclear quantum dynamics in high-dimensional systems. By providing time-dependent variational orbitals and an optimal choice of layered effective degrees of freedom, one is able to reduce the computational cost to an amenable number of configurations. However, choices related to selecting properly the mode grouping and tensor tree are strongly system dependent and, thus far, subjectively based on intuition and/or experience. Therefore, herein we detail a new protocol based on multivariate statistics─more specifically, factor analysis and hierarchical clustering─for a reliable and convenient guiding in the optimal design of such complex "system-of-systems" tensor-network decompositions. The advantages of employing the new algorithm and its applicability are tested on water and two floppy protonated water clusters with large amplitude motions.

Rovibrational states calculations of the H2O-HCN heterodimer with the multiconfiguration time dependent Hartree method.

Water and hydrogen cyanide are two of the most common species in space and the atmosphere with the ability of binding to form dimers such as H2O-HCN. In the literature, while calculations characterizing various properties of the H2O-HCN cluster (equilibrium distance, vibrational frequencies and rotational constants) have been done in the past, extensive calculations of the rovibrational states of this system using a reliable quantum dynamical approach have yet to be reported. In this work, we intend to mend that by performing the first calculation of the rovibrational states of the H2O-HCN van der Waals complex on a recently developed potential energy surface. We use the block improved relaxation procedure implemented in the Heidelberg MultiConfiguration Time-Dependent Hartree (MCTDH) package to compute the states of the H2O-HCN isomer, from which we extract the transition frequencies and rotational constants of the complex. We further adapt an approach first suggested by Wang and Carrington-and supported here by analysis routines of the Heidelberg MCTDH package-to properly characterize the computed rovibrational states. The subsequent assignment of rovibrational states was done by theoretical analysis and visual inspection of the wavefunctions. Our simulations provide a Zero Point Energy (ZPE) and intermolecular vibrational frequencies in good agreement with past ab initio calculations. The transition frequencies and rotational constants obtained from our simulations match well with the available experimental data. This work has the broad aim to propose the MCTDH approach as a reliable option to compute and characterize rovibrational states of van der Waals complexes such as the current one.

Construction of Beyond Born-Oppenheimer Based Diabatic Surfaces and Generation of Photoabsorption Spectra: The Touchstone Pyrazine (C4 N2 H4 ).

We construct theoretically "exact" and numerically "accurate" Beyond Born-Oppenheimer (BBO) based diabatic potential energy surfaces (PESs) of pyrazine (C4 N2 H4 ) molecule involving lowest four excited adiabatic PESs (S1 to S4 ) and nonadiabatic coupling terms (NACTs) among those surfaces as functions of nonadiabatically active normal modes (Q1 , Q6a , Q9a and Q10a ) to compute its photoabsorption (PA) spectra. Those adiabatic PESs are calculated using CASSCF as well as MRCI based methodologies, where NACTs are obtained from CP-MCSCF approach. Employing ab initio quantities (adiabatic PESs and NACTs), it is possible to depict the conical intersections (CIs) and develop matrices of diabatic PESs over six normal mode planes. Once single-valued, smooth, symmetric and continuous 2×2 and 4×4 diabatic surface matrices are in hand for the first time, such matrices are used to perform multi-state multi-mode nuclear dynamics with the aid of Time-Dependent Discrete Variable Representation (TDDVR) methodology initializing the product type wavefunction on (S1 ) and (S2 ) states to obtain the corresponding PA spectra. TDDVR calculated spectra for those states (S1 and S2 ) obtained from BBO based 2×2 and 4×4 diabatic surface matrices show good and better agreement with the experimental results, respectively. Both of these calculated results depict better peak progression over the existing profiles of Multi-Configuration Time-Dependent Hartree (MCTDH) dynamics over 2×2 Vibronic Coupling Model (VCM) Hamiltonian.

Stochastic multi-configuration time-dependent Hartree for dissipative quantum dynamics with strong intramolecular coupling.

In this article, we explore the dissipation dynamics of a strongly coupled multidimensional system in contact with a Markovian bath, following a system-bath approach. We use in this endeavor the recently developed stochastic multi-configuration time-dependent Hartree approach within the Monte Carlo wave packet formalism [S. Mandal et al., J. Chem. Phys. 156, 094109 (2022)]. The method proved to yield thermalized ensembles of wave packets when intramolecular coupling is weak. To treat strongly coupled systems, new Lindblad dissipative operators are constructed as linear combinations of the system coordinates and associated momenta. These are obtained by a unitary transformation to a normal mode representation, which reduces intermode coupling up to second order. Additionally, we use combinations of generalized raising/lowering operators to enforce the Boltzmann distribution in the dissipation operators, which yield perfect thermalization in the harmonic limit. The two ansatz are tested using a model two-dimensional Hamiltonian, parameterized to disentangle the effects of intramolecular potential coupling, of strong mode mixing observed in Fermi resonances, and of anharmonicity.

On the UV spectroscopy and photodynamics of octatetraene

The UV absorption spectrum of all-trans-octatetraene is reinvestigated theoretically, focussing on the strongly dipole-allowed – transition. The dynamical calculations rely on ab initio multistate CASPT2 (MS-CASPT2) computations of the underlying potential energy surfaces and coupling elements, and the multiconfiguration time-dependent Hartree (MCTDH) method for solving the time-dependent Schrödinger equation for the nuclear motion. The vibronic structure of the absorption band near 2700 Å is well reproduced, better than ever before by a purely ab initio approach. The C–C single bond and double bond stretching modes dominate the vibrational structure of the transition. While the latter resembles nearly harmonic motion on weakly coupled potential energy surfaces, the time-dependent calculations reveal a subpicosecond electronic population transfer from the state to the lower-lying state proceeding on a time scale of . The details of this internal conversion process and its relation to the conical intersection between the two excited states are discussed, paying particular attention to the impact of the – energy gap on the dynamics.

Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine.

We present a viable protocol to compute vibrational resonance Raman (vRR) spectra for systems with several close-lying and potentially coupled electronic states. It is based on the parametrization of linear vibronic coupling (LVC) models from time-dependent density functional theory (TD-DFT) calculations and quantum dynamics propagations of vibronic wavepackets with the multilayer version of the multiconfiguration time-dependent Hartree (ML-MCTDH) method. Our approach is applied to thymine considering seven coupled electronic states, comprising the three lowest bright states, and all vibrational coordinates. Computed vRR at different excitation wavelengths are in good agreement with the available experimental data. Up to 250 nm the signal is dominated by the lowest HOMO → LUMO transition, whereas at 233 nm, in the valley between the two lowest energy absorption bands, the contributions of all the three bright states, and their interferences and couplings, are important. Inclusion of solvent (water) effects improves the agreement with experiment, reproducing the coalescence of vibrational bands due to CC and C=O stretchings. With our approach we disentangle and assess the effect of interferences between the contribution of different quasi-resonant states to the transition polarizability and the effect of interstate couplings. Our findings strongly suggest that in cases of close-lying and potentially coupled states a simple inclusion of interference effects is not sufficient, and a fully nonadiabatic computation should instead be performed. We also document that for systems with strong couplings and quasi-degenerate states, the use of HT perturbative approach, not designed for these cases, may lead to large artifacts.

Multilayer Multiconfiguration Time-Dependent Hartree Study on the Mode-/Bond-Specific Quantum Dynamics of Water Dissociation on Cu(111).

In this work, full-dimensional (9D) quantum dynamics calculations on mode-/bond-specific surface scattering of a water molecule on a copper (111) rigid surface are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. To easily perform the ML-MCTDH calculations on such a triatomic molecule-surface system, we first choose specific Jacobi coordinates as a set of coordinates of water. Next, to efficiently perform the 9D ML-MCTDH wavepacket propagation, the potential energy surface is transferred to a canonical polyadic decomposition form with the aid of a Monte Carlo-based method. Excitation-specific dissociation probabilities of H2O on Cu(111) are computed, and mode-/bond-specific dynamics are demonstrated by comparison with a probability curve computed for a water molecule in the ground state. The dependence of the dissociation probability of the initial state of H2O is studied, and it is found that the excitation-specific dissociation probabilities can be divided into three groups. We find that the vibrationally excited states enhance the dissociation reactivity of H2O, while the rotationally excited states hardly influence it.

Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs H-aggregate models.

We report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature using the Multi-Layer Multi-Configuration Time-Dependent Hartree method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 × 10-3cm2/s, showing a non-linear increase with temperature.

Modelling ultrafast dynamics at a conical intersection with regularized diabatic states: An approach based on multiplicative neural networks

• Neural network (NN) potentials for quantum dynamics of cis–trans photoisomerization. • Multiplicative NNs adapted to regularized diabatic states at a conical intersection. • NN potentials match a sum-of-products form required by MCTDH wavefunction propagation. • Good agreement found for highly correlated dynamics with up to 13 degrees of freedom. Neural-network (NN) potentials are employed in conjunction with the Multiconfiguration Time-Dependent Hartree (MCTDH) method in order to simulate an ultrafast photoinduced cis–trans type isomerization process induced by a conical intersection. To this end, NN potentials are fitted to a diabatic potential of regularized diabatic states type [Köppel et al., J. Chem. Phys. 115, 2377 (2001)], which entirely relies on adiabatic potential information. Multiplicative NNs are employed which match the sum-of-products form of the multiconfigurational MCTDH wavefunction. Good agreement with the reference dynamics is obtained for simulations of the highly correlated dynamics with up to 13 degrees of freedom. This study contributes to developing NN methodologies suitable for photochemical dynamics at complex excited-state topologies.

Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree.

Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene–oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10–100 fs time scale via multiple “hot-exciton dissociation” pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.

Representation of Diabatic Potential Energy Matrices for Multiconfiguration Time-Dependent Hartree Treatments of High-Dimensional Nonadiabatic Photodissociation Dynamics.

Conventional quantum mechanical characterization of photodissociation dynamics is restricted by steep scaling laws with respect to the dimensionality of the system. In this work, we examine the applicability of the multi-configurational time-dependent Hartree (MCTDH) method in treating nonadiabatic photodissociation dynamics in two prototypical systems, taking advantage of its favorable scaling laws. To conform to the sum-of-product form, elements of the ab initio diabatic potential energy matrix (DPEM) are re-expressed using the recently proposed Monte Carlo canonical polyadic decomposition method, with enforcement of proper symmetry. The MCTDH absorption spectra and product branching ratios are shown to compare well with those calculated using conventional grid-based methods, demonstrating its promise for treating high-dimensional nonadiabatic photodissociation problems.

Non-Markovian vibrational relaxation dynamics at surfaces.

Vibrational dynamics of adsorbates near surfaces plays both an important role for applied surface science and as a model lab for studying fundamental problems of open quantum systems. We employ a previously developed model for the relaxation of a D-Si-Si bending mode at a D:Si(100)-(2 × 1) surface, induced by a "bath" of more than 2000 phonon modes [Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], to extend previous work along various directions. First, we use a Hierarchical Effective Mode (HEM) model [Fischer et al., J. Chem. Phys. 153, 064704 (2020)] to study relaxation of higher excited vibrational states than hitherto done by solving a high-dimensional system-bath time-dependent Schrödinger equation (TDSE). In the HEM approach, (many) real bath modes are replaced by (much less) effective bath modes. Accordingly, we are able to examine scaling laws for vibrational relaxation lifetimes for a realistic surface science problem. Second, we compare the performance of the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) approach with that of the recently developed coherent-state-based multi-Davydov-D2 Ansatz [Zhou et al., J. Chem. Phys. 143, 014113 (2015)]. Both approaches work well, with some computational advantages for the latter in the presented context. Third, we apply open-system density matrix theory in comparison with basically "exact" solutions of the multi-mode TDSEs. Specifically, we use an open-system Liouville-von Neumann (LvN) equation treating vibration-phonon coupling as Markovian dissipation in Lindblad form to quantify effects beyond the Born-Markov approximation.

Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials.

The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.