The multi-configuration electron-nuclear dynamics method applied to LiH

Implementation of the Time-Dependent Complete-Active-Space Self-Consistent-Field Method for Diatomic Molecules

We present a computational approach that implements the time-dependent complete-active-space self-consistent-field method, as introduced in [Phys. Rev. A 88, 023402 (2013)]. Our implementation addresses the challenge of diatomic molecules subjected to an intense laser pulse by considering the full dimensionality of the problem using prolate spheroidal coordinates. The method incorporates the gauge-invariant frozen-core approximation, boosts the evaluation of the electron-electron interaction term using finite-element discrete-variable representation with Neumann expansion, and utilizes an exponential time differencing scheme tailored for the stable propagation of the stiff nonlinear orbital functions. We have successfully applied this methodology to study high-harmonic generation in diatomic molecules such as H2, LiH, and N2, shedding light on the impact of electron correlations in these systems.

Gaussian-Type Orbital Calculations for High Harmonic Generation in Vibrating Molecules: Benchmarks for H2. J Chem Theory Comput 2021;17:7353-7365. [PMID: 34747605 DOI: 10.1021/acs.jctc.1c00837]

The response of the hydrogen molecular ion, H₂⁺, to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrödinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.

Conditional Wave Function Theory: A Unified Treatment of Molecular Structure and Nonadiabatic Dynamics

We demonstrate that a conditional wave function theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wave function formally recasts the full interacting wave function of a closed system as a set of lower-dimensional (conditional) coupled "slices". We formulate a variational wave function ansatz based on a set of conditional wave function slices and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wave functions and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser-driven proton transfer, and nuclear quantum effects induced by a conical intersection. This work paves the road for the application of conditional wave function theory in equilibrium and out-of-equilibrium ab initio molecular simulations of finite and extended systems.

Multi-configuration electron-nuclear dynamics: An open-shell approach

The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.

Non-adiabatic quantum dynamics without potential energy surfaces based on second-quantized electrons: Application within the framework of the MCTDH method

A first principles quantum formalism to describe the non-adiabatic dynamics of electrons and nuclei based on a second quantization representation (SQR) of the electronic motion combined with the usual representation of the nuclear coordinates is introduced. This procedure circumvents the introduction of potential energy surfaces and non-adiabatic couplings, providing an alternative to the Born-Oppenheimer approximation. An important feature of the molecular Hamiltonian in the mixed first quantized representation for the nuclei and the SQR representation for the electrons is that all degrees of freedom, nuclear positions and electronic occupations, are distinguishable. This makes the approach compatible with various tensor decomposition Ansätze for the propagation of the nuclear-electronic wavefunction. Here, we describe the application of this formalism within the multi-configuration time-dependent Hartree framework and its multilayer generalization, corresponding to Tucker and hierarchical Tucker tensor decompositions of the wavefunction, respectively. The approach is applied to the calculation of the photodissociation cross section of the HeH⁺ molecule under extreme ultraviolet irradiation, which features non-adiabatic effects and quantum interferences between the two possible fragmentation channels, He + H⁺ and He⁺ + H. These calculations are compared with the usual description based on ab initio potential energy surfaces and non-adiabatic coupling matrix elements, which fully agree. The proof-of-principle calculations serve to illustrate the advantages and drawbacks of this formalism, which are discussed in detail, as well as possible ways to overcome them. We close with an outlook of possible application domains where the formalism might outperform the usual approach, for example, in situations that combine a strong static correlation of the electrons with non-adiabatic electronic-nuclear effects.

Zombie states for description of structure and dynamics of multi-electron systems

Canonical Coherent States (CSs) of Harmonic Oscillator have been extensively used as a basis in a number of computational methods of quantum dynamics. However, generalising such techniques for fermionic systems is difficult because Fermionic Coherent States (FCSs) require complicated algebra of Grassmann numbers not well suited for numerical calculations. This paper introduces a coherent antisymmetrised superposition of "dead" and "alive" electronic states called here Zombie State (ZS), which can be used in a manner of FCSs but without Grassmann algebra. Instead, for Zombie States, a very simple sign-changing rule is used in the definition of creation and annihilation operators. Then, calculation of electronic structure Hamiltonian matrix elements between two ZSs becomes very simple and a straightforward technique for time propagation of fermionic wave functions can be developed. By analogy with the existing methods based on Canonical Coherent States of Harmonic Oscillator, fermionic wave functions can be propagated using a set of randomly selected Zombie States as a basis. As a proof of principles, the proposed Coupled Zombie States approach is tested on a simple example showing that the technique is exact.

Pumping and probing vibrational modulated coupled electronic coherence in HCN using short UV fs laser pulses: a 2D quantum nuclear dynamical study

The coupled electronic-nuclear coherent dynamics induced by a short strong VUV fs pulse in the low excited electronic states of HCN is probed by transient absorption spectroscopy with a second weaker fs UV pulse. The nuclear time-dependent Schrodinger equation is solved on a 2D nuclear grid with several electronic states with a Hamiltonian including the dipole coupling to the pump and the probe electric fields. The two internal nuclear coordinates describe the motion of the light H atom. There is a band of several excited electronic states at about 8 eV above the ground state (GS) that is transiently accessed by the pump pulse. We tailored the pump so as to selectively populate the lowest 1A'' electronic state thereby the pulse creates an electronic coherence with the GS. Our simulations show that this electronic coherence is modulated by the nuclear motion and persists all the way to dissociation on the 1A'' state. Transient absorption spectra computed as a function of the delay time between the pump and the probe pulses provide a detailed probe of the electronic amplitude and its phase, as well as of the modulation of the electronic coherence by the nuclear motion, both bound and dissociative.

Wavepacket dynamics and the multi-configurational time-dependent Hartree approach

Multi-configurational time-dependent Hartree (MCTDH) based approaches are efficient, accurate, and versatile methods for high-dimensional quantum dynamics simulations. Applications range from detailed investigations of polyatomic reaction processes in the gas phase to high-dimensional simulations studying the dynamics of condensed phase systems described by typical solid state physics model Hamiltonians. The present article presents an overview of the different areas of application and provides a comprehensive review of the underlying theory. The concepts and guiding ideas underlying the MCTDH approach and its multi-mode and multi-layer extensions are discussed in detail. The general structure of the equations of motion is highlighted. The representation of the Hamiltonian and the correlated discrete variable representation (CDVR), which provides an efficient multi-dimensional quadrature in MCTDH calculations, are discussed. Methods which facilitate the calculation of eigenstates, the evaluation of correlation functions, and the efficient representation of thermal ensembles in MCTDH calculations are described. Different schemes for the treatment of indistinguishable particles in MCTDH calculations and recent developments towards a unified multi-layer MCTDH theory for systems including bosons and fermions are discussed.

A fully general time-dependent multiconfiguration self-consistent-field method for the electron–nuclear dynamics

A time-dependent multiconfiguration self-consistent-field method for a system consisting of arbitrarily different kinds and numbers of interacting fermions and bosons.

Quantum Nuclear Dynamics Pumped and Probed by Ultrafast Polarization Controlled Steering of a Coherent Electronic State in LiH

The quantum wave packet dynamics following a coherent electronic excitation of LiH by an ultrashort, polarized, strong one-cycle infrared optical pulse is computed on several electronic states using a grid method. The coupling to the strong field of the pump and the probe pulses is included in the Hamiltonian used to solve the time-dependent Schrodinger equation. The polarization of the pump pulse allows us to control the localization in time and in space of the nonequilibrium coherent electronic motion and the subsequent nuclear dynamics. We show that transient absorption, resulting from the interaction of the total molecular dipole with the electric fields of the pump and the probe, is a very versatile probe of the different time scales of the vibronic dynamics. It allows probing both the ultrashort, femtosecond time scale of the electronic coherences as well as the longer dozens of femtoseconds time scales of the nuclear motion on the excited electronic states. The ultrafast beatings of the electronic coherences in space and in time are shown to be modulated by the different periods of the nuclear motion.