Time-dependent Hartree Research Articles

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.

Multielectron effects in CO ionization by orthogonal two-color laser fields

We study the single ionization of CO molecule in orthogonally polarized two-color laser fields using time-dependent Hartree–Fock method. Employing the coherent superposition consisting of a fundamental field and its second harmonic, we consider different wavelengths (780 and 390nm) and intensities (1×1014, 1.5×1014 and 2×1014 W/cm2) of the fundamental field for which different ionization mechanisms are involved. Our numerical results reveal that the contributions of the three highest occupied molecular orbitals are strongly influenced by the intensity ratio and time delay between the two components of the bichromatic field. At the fundamental wavelength of 780nm and intensity of 2×1014 W/cm2, the two outermost molecular orbitals display a dominant response to the orthogonal two-color laser field when the second harmonic field is stronger than the fundamental field.

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.

Various variational approximations of quantum dynamics

We investigate variational principles for the approximation of quantum dynamics that apply for approximation manifolds that do not have complex linear tangent spaces. The first one, dating back to McLachlan [Mol. Phys. 8, 39–44 (1964)], minimizes the residuum of the time-dependent Schrödinger equation, while the second one, originating from the lecture notes of Kramer and Saraceno [Geometry of the Time-Dependent Variational Principle in Quantum Mechanics, Lecture Notes in Physics Vol. 140 (Springer, Berlin, 1981)], imposes the stationarity of an action functional. We characterize both principles in terms of metric and symplectic orthogonality conditions, consider their conservation properties, and derive an elementary a posteriori error estimate. As an application, we revisit the time-dependent Hartree approximation and frozen Gaussian wave packets.

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.

The time-dependent Hartree–Fock–Bogoliubov equations for Bosons

We introduce the map of dynamics of quantum Bose gases into dynamics of quasifree states, which we call the “nonlinear quasifree approximation”. We use this map to derive the time-dependent Hartree–Fock–Bogoliubov (HFB) equations describing the dynamics of quantum fluctuations around a Bose–Einstein condensate. We prove global well-posedness of the HFB equations for pair potentials satisfying suitable regularity conditions, and we establish important conservation laws. We show that the space of solutions of the HFB equations has a symplectic structure reminiscent of a Hamiltonian system. This is then used to relate the HFB equations to the HFB eigenvalue equations discussed in the physics literature. We also construct Gibbs equilibrium states at positive temperature associated with the HFB equations, and we establish criteria for the appearance of Bose–Einstein condensation.

Numerically Precise Benchmark of Many-Body Self-Energies on Spherical Atoms.

We investigate the performance of beyond-GW approaches in many-body perturbation theory by addressing atoms described within the spherical approximation via a dedicated numerical treatment based on B-splines and spherical harmonics. We consider the GW, second Born (2B), and GW + second order screened exchange (GW+SOSEX) self-energies and use them to obtain ionization potentials from the quasi-particle equation (QPE) solved perturbatively on top of independent-particle calculations. We also solve the linearized Sham–Schlüter equation (LSSE) and compare the resulting xc potentials against exact data. We find that the LSSE provides consistent starting points for the QPE but does not present any practical advantage in the present context. Still, the features of the xc potentials obtained with it shed light on possible strategies for the inclusion of beyond-GW diagrams in the many-body self-energy. Our findings show that solving the QPE with the GW+SOSEX self-energy on top of a PBE or PBE0 solution is a viable scheme to go beyond GW in finite systems, even in the atomic limit. However, GW shows a comparable performance if one agrees to use a hybrid starting point. We also obtain promising results with the 2B self-energy on top of Hartree–Fock, suggesting that the full time-dependent Hartree–Fock vertex may be another viable beyond-GW scheme for finite systems.

Dynamical approximations for composite quantum systems: assessment of error estimates for a separable ansatz

Numerical studies are presented to assess error estimates for a separable (Hartree) approximation for dynamically evolving composite quantum systems which exhibit distinct scales defined by their mass and frequency ratios. The relevant error estimates were formally described in our previous work Burghardt et al (2021 J. Phys. A: Math. Theor. 54 414002). Specifically, we consider a representative two-dimensional tunneling system where a double well and a harmonic coordinate are cubically coupled. The time-dependent Hartree approximation is compared with a fully correlated solution, for different parameter regimes. The impact of the coupling and the resulting correlations are quantitatively assessed in terms of a time-dependent reaction probability along the tunneling coordinate. We show that the numerical error is correctly predicted on moderate time scales by a theoretically derived error estimate.

An efficient implementation of time-dependent linear-response theory for strongly orthogonal geminal wave function models.

We present an implementation of time-dependent linear-response equations for strongly orthogonal geminal wave function models: the time-dependent generalized valence bond (TD-GVB) perfect-pairing theory and the antisymmetrized product of strongly orthogonal geminals. The geminal wave functions are optimized using a restricted-step second-order algorithm suitable for handling many geminals, and the linear-response equations are solved in an efficient way using a direct iterative approach. The wave function optimization algorithm features an original scheme to create initial orbitals for the geminal functions in a black-box fashion. The implementation is employed to examine the accuracy of the geminal linear response for singlet excitation energies of small and medium-sized molecules. In systems dominated by dynamic correlation, geminal models constitute only a minor improvement with respect to time-dependent Hartree-Fock. Compared to the linear-response complete active space self-consistent field, TD-GVB either misses or gives large errors for states dominated by double excitations.

Adaptive fitting of potential energy surfaces of small to medium-sized molecules in sum-of-product form: Application to vibrational spectroscopy.

A semi-automatic sampling and fitting procedure for generating sum-of-product (Born-Oppenheimer) potential energy surfaces based on a high-dimensional model representation is presented. The adaptive sampling procedure and subsequent fitting rely on energies only and can be used for re-fitting existing analytic potential energy surfaces in the sum-of-product form or for direct fits from ab initio computations. The method is tested by fitting ground electronic state potential energy surfaces for small to medium sized semi-rigid molecules, i.e., HFCO, HONO, and HCOOH, based on ab initio computations at the coupled-cluster single double and perturbative triples-F12/cc-pVTZ-F12 or MP2/aug-cc-pVTZ levels of theory. Vibrational eigenstates are computed using block improved relaxation in the Heidelberg multi-configurational time dependent Hartree package and compared to available experimental and theoretical data. The new potential energy surfaces are compared to the best ones currently available for these molecules in terms of accuracy, including resulting vibrational states, required number of sampling points, and number of fitting parameters. The present procedure leads to compact expansions and scales well with the number of dimensions for simple potentials such as single or double wells.

Origin invariant electronic circular dichroism in the length dipole gauge without London atomic orbitals.

We present a method for obtaining origin-independent electronic circular dichroism (ECD) in the length-gauge representation LG(OI) without the usage of London atomic orbitals. This approach builds upon the work by Caricato [J. Chem. Phys. 153, 151101 (2020)] and is applied to rotatory strengths and ECD spectra from damped response theory. Numerical results are presented for time-dependent Hartree-Fock and density-functional theory, the second-order algebraic diagrammatic construction method, and linear-response coupled-cluster theory with singles and approximate doubles. We can support the finding that the common choice of placing the gauge origin in the center of mass of a molecule in conventional length-gauge calculations involving chiroptical properties might not be optimal and show that LG(OI) is a valuable alternative for the origin-independent calculation of ECD spectra. We show that, for a limited test set, the convergence of the rotatory strengths calculated with the LG(OI) approach toward the basis-set limit tends to be faster than for the established velocity gauge representation. Relationships between the sum-over-states expression of the optical rotation in the LG(OI) framework and its representation in terms of response functions are analyzed.

Breathing dynamics of the few-body Bose polaron in a species-selective harmonic trap

We perform an extensive numerical study on the breathing dynamics of a few-body Bose polaron setup in a one-dimensional species-selective harmonic trap. The dynamics is triggered by a quench of the impurity trap. The excitation of the background majority atoms is mediated via the majority-impurity interaction. The breathing spectrum is obtained for different numbers of majority particles, several values of the majority-component interaction strengths and trap ratios. It is further compared to the breathing spectrum of a particle-balanced few-body Bose-Bose mixture. In particular, for equal post-quench traps the employed protocol allows to couple states of different center-of-mass parity in contrast to species-symmetric trap quenches. Among the participating eigenstates we identify one having odd center-of-mass parity and even global parity. The breathing frequency induced by this state is a monotonically decreasing function of the coupling parameter. Importantly, in order to be observable it requires the entanglement between the species to be taken into account. We demonstrate this by comparing the numerically exact results obtained by means of the Multi-Layer Multi-Configuration Time-Dependent Hartree Method for Mixtures to the ones of a species mean-field ansatz. The entanglement-sensitive breathing frequency persists also for unequal post-quench traps where the center-of-mass cannot be decoupled. Finally, we analyze the impact of parity symmetry on the breathing dynamics by initializing a state of odd global parity. We evidence a striking resemblance to the ground state breathing dynamics.

A non-hierarchical correlation discrete variable representation.

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. The non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH3 demonstrate the accuracy of the non-hierarchical CDVR approach.

Spin-vibronic coupling in the quantum dynamics of a Fe(III) trigonal-bipyramidal complex.

The presence of a high density of excited electronic states in the immediate vicinity of the optically bright state of a molecule paves the way for numerous photo-relaxation channels. In transition-metal complexes, the presence of heavy atoms results in a stronger spin-orbit coupling, which enables spin forbidden spin-crossover processes to compete with the spin-allowed internal conversion processes. However, no matter how effectively the states cross around the Franck-Condon region, the degree of vibronic coupling, of both relativistic and non-relativistic nature, drives the population distribution among these states. One such case is demonstrated in this work for the intermediate-spin Fe(III) trigonal-bipyramidal complex. A quantum dynamical investigation of the photo-deactivation mechanism in the Fe(III) system is presented using the multi-configurational time-dependent Hartree approach based on the vibronic Hamiltonian whose coupling terms are derived from the state-averaged complete active space self-consistent field/complete active space with second-order perturbation theory (CASPT2) calculations and spin-orbit coupling of the scalar-relativistic CASPT2 states. The results of this study show that the presence of a strong (non-relativistic) vibronic coupling between the optically bright intermediate-spin state and other low-lying states of the same spin-multiplicity overpowers the spin-orbit coupling between the intermediate-spin and high-spin states, thereby lowering the chances of spin-crossover while exhibiting ultrafast relaxation among the intermediate-spin states. In a special case, where the population transfer pathway via the non-relativistic vibronic coupling is blocked, the probability of the spin-crossover is found to increase. This suggests that a careful modification of the complex by incorporation of heavier atoms with stronger relativistic effects can enhance the spin-crossover potential of Fe(III) intermediate-spin complexes.

Mixed-quantum-classical or fully-quantized dynamics? A unified code to compare methods.

Three methods for non-adiabatic dynamics are compared to highlight their capabilities. Multi-configurational time-dependent Hartree is a full grid-based solution to the time-dependent Schrödinger equation, variational multi-configurational Gaussian (vMCG) uses a less flexible but unrestricted Gaussian wavepacket basis, and trajectory surface hopping (TSH) replaces the nuclear wavepacket with a swarm of classical trajectories. Calculations with all methods using a model Hamiltonian were performed. The vMCG and TSH were also then run in a direct dynamics mode, with the potential energy surfaces calculated on-the-fly using quantum chemistry calculations. All dynamics calculations used the Quantics package, with the TSH calculations using a new interface to a surface hopping code. A novel approach to calculate adiabatic populations from grid-based quantum dynamics using a time-dependent discrete variable representation is presented, allowing a proper comparison of methods. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Multidimensional stochastic dissipative quantum dynamics using a Lindblad operator.

In this paper, multidimensional dissipative quantum dynamics is studied within a system-bath approach in the Markovian regime using a model Lindblad operator. We report on the implementation of a Monte Carlo wave packet algorithm in the Heidelberg version of the Multi-Configuration Time-Dependent Hartree (MCTDH) program package, which is henceforth extended to treat stochastic dissipative dynamics. The Lindblad operator is represented as a sum of products of one-dimensional operators. The new form of the operator is not restricted to the MCTDH formalism and could be used with other multidimensional quantum dynamical methods. As a benchmark system, a two-dimensional coupled oscillators model representing the internal stretch and the surface-molecule distance in the O2/Pt(111) system coupled to a Markovian bath of electron-hole-pairs is used. The simulations reveal the interplay between coherent intramolecular coupling due to anharmonic terms in the potential and incoherent relaxation due to coupling to an environment. It is found that thermalization of the system can be approximately achieved when the intramolecular coupling is weak.

The multi-configurational time-dependent Hartree approach in optimized second quantization: Thermal ensembles and statistical sampling

The recently introduced multi-layer (ML) multi-configurational time-dependent Hartree (MCTDH) in optimized second quantization representation (oSQR) approach combines the use of an optimized time-dependent orbital basis and the multi-layer MCTDH tensor contraction scheme to efficiently describe the quantum dynamics of system consisting of identical particles. Extending the original work on the subject, which studied propagation in real and imaginary time for Hamiltonians with flexible and fixed particle numbers, here the representation of thermal ensembles and the calculation of thermally averaged observables within the ML-MCTDH-oSQR is investigated. Different statistal sampling schemes for fixed and flexible particle numbers are introduced and analyzed. Numerical calculations studying a Hubbard model are used to illustrate the properties of the different schemes.