Quantum dynamics simulations using Gaussian wavepackets: the vMCG method

Highly efficient surface hopping dynamics using a linear vibronic coupling model

We report an implementation of the linear vibronic coupling (LVC) model within the surface hopping dynamics approach and present utilities for parameterizing this model in a blackbox fashion. This results in an extremely efficient method to obtain qualitative and even semi-quantitative information about the photodynamical behavior of a molecule, and provides a new route toward benchmarking the results of surface hopping computations. The merits and applicability of the method are demonstrated in a number of applications. First, the method is applied to the SO2 molecule showing that it is possible to compute its absorption spectrum beyond the Condon approximation, and that all the main features and timescales of previous on-the-fly dynamics simulations of intersystem crossing are reproduced while reducing the computational effort by three orders of magnitude. The dynamics results are benchmarked against exact wavepacket propagations on the same LVC potentials and against a variation of the electronic structure level. Four additional test cases are presented to exemplify the broader applicability of the model. The photodynamics of the isomeric adenine and 2-aminopurine molecules are studied and it is shown that the LVC model correctly predicts ultrafast decay in the former and an extended excited-state lifetime in the latter. Futhermore, the method correctly predicts ultrafast intersystem crossing in the modified nucleobase 2-thiocytosine and its absence in 5-azacytosine while it fails to describe the ultrafast internal conversion to the ground state in the latter.

Modified Gaussian Wave Packet Method for Calculating Initial State Wave Functions in Photodissociation

A modified Gaussian wave packet relaxation method is proposed to calculate the ground state wave function using an expansion of frozen Gaussian wave packets. This new procedure consists of two steps. In the first step, a multidimensional Gaussian product placed at the ground state equilibrium geometry is propagated in imaginary time. The relaxation optimizes the widths of the one-dimensional Gaussians. In the second step, additional Gaussian wave packets with the same widths are placed near the equilibrium geometry, and the corresponding expansion coefficients are optimized using the same relaxation method. This new algorithm is tested in photodissociation of NOCl and NH₃, and the results show good agreement with the exact results in the energy, wave function, and absorption spectrum. In particular, the highly structured absorption spectrum of NH₃ is reproduced, underscoring the accuracy of both the initial wave packet and the excited state propagation.

Nonadiabatic dynamics: The SHARC approach

We review the Surface Hopping including ARbitrary Couplings (SHARC) approach for excited-state nonadiabatic dynamics simulations. As a generalization of the popular surface hopping method, SHARC allows simulating the full-dimensional dynamics of molecules including any type of coupling terms beyond nonadiabatic couplings. Examples of these arbitrary couplings include spin-orbit couplings or dipole moment-laser field couplings, such that SHARC can describe ultrafast internal conversion, intersystem crossing, and radiative processes. The key step of the SHARC approach consists of a diagonalization of the Hamiltonian including these couplings, such that the nuclear dynamics is carried out on potential energy surfaces including the effects of the couplings-this is critical in any applications considering, for example, transition metal complexes or strong laser fields. We also give an overview over the new SHARC2.0 dynamics software package, released under the GNU General Public License, which implements the SHARC approach and several analysis tools. The review closes with a brief survey of applications where SHARC was employed to study the nonadiabatic dynamics of a wide range of molecular systems. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and KineticsSoftware > Simulation MethodsSoftware > Quantum Chemistry.

Computation of the S1 ← S0 Vibronic Absorption Spectrum of Formaldehyde by Variational Gaussian Wavepacket and Semiclassical IVR Methods

The vibronic absorption spectrum of the electric dipole forbidden and vibronically allowed S₁(¹ A₂) ← S₀(¹ A₁) transition of formaldehyde is calculated by Gaussian wavepacket and semiclassical methods, along with numerically exact reference calculations, using the potential energy surface of Fu, Shepler, and Bowman ( J. Am. Chem. Soc. 2011, 133, 7957). Specifically, the variational multiconfigurational Gaussian (vMCG) approach and the Herman-Kluk semiclassical initial value representation (HK-SCIVR) are compared to assess the accuracy and convergence of these methods, benchmarked against numerically exact time-dependent wavepacket propagation (TDWP) on the reference potential energy surface. The vMCG calculation is shown to converge quite well with about 100 variationally evolving Gaussian functions and using a local cubic expansion instead of the conventional local harmonic approximation. By contrast, the HK-SCIVR approach with ∼10⁵ trajectories reproduces the vibrationally structured spectral envelope correctly but yields a strongly broadened spectrum. The comparison of the computed absorption spectrum with experiment shows that the relevant vibronic progressions are reasonably reproduced by all computations, but deviations of the order of 10-100 cm^-1 occur, underscoring that both electronic structure calculations and dynamical approaches remain challenging in the calculation of typical small-molecule excited-state spectra by trajectory-based methods.

The Ehrenfest method with fully quantum nuclear motion (Qu-Eh): Application to charge migration in radical cations

An algorithm is described for quantum dynamics where an Ehrenfest potential is combined with fully quantum nuclear motion (Quantum-Ehrenfest, Qu-Eh). The method is related to the single-set variational multi-configuration Gaussian approach (vMCG) but has the advantage that only a single quantum chemistry computation is required at each time step since there is only a single time-dependent potential surface. Also shown is the close relationship to the "exact factorization method." The quantum Ehrenfest method is compared with vMCG for study of electron dynamics in a modified bismethylene-adamantane cation system. Illustrative examples of electron-nuclear dynamics are presented for a distorted allene system and for HCCI⁺ where one has a degenerate Π system.

A walk through the approximations of ab initio multiple spawning

Full multiple spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the ab initio multiple spawning method which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the ab initio multiple spawning approximations on different observables and (ii) the convergence of the ab initio multiple spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying ab initio multiple spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.

On-the-Fly Ab Initio Semiclassical Evaluation of Absorption Spectra of Polyatomic Molecules beyond the Condon Approximation

To evaluate vibronic spectra beyond the Condon approximation, we extend the on-the-fly ab initio thawed Gaussian approximation by considering the Herzberg-Teller contribution due to the dependence of the electronic transition dipole moment on nuclear coordinates. The extended thawed Gaussian approximation is tested on electronic absorption spectra of the phenyl radical and benzene; calculated spectra reproduce experimental data and are much more accurate than standard global harmonic approaches, confirming the significance of anharmonicity. Moreover, the extended method provides a tool to quantify the Herzberg-Teller contribution; we show that in the phenyl radical, anharmonicity outweighs the Herzberg-Teller contribution, whereas in benzene, the Herzberg-Teller contribution is essential because the transition is electronically forbidden and the Condon approximation yields a zero spectrum. Surprisingly, both adiabatic harmonic spectra outperform those of the vertical harmonic model, which describes the Franck-Condon region better. Finally, we provide a simple recipe for orientationally averaging spectra, valid beyond the Condon approximation, and a relation among the transition dipole, its gradient, and nonadiabatic coupling vectors.

Variational nonadiabatic dynamics in the moving crude adiabatic representation: Further merging of nuclear dynamics and electronic structure

A new methodology of simulating nonadiabatic dynamics using frozen-width Gaussian wavepackets within the moving crude adiabatic representation with the on-the-fly evaluation of electronic structure is presented. The main feature of the new approach is the elimination of any global or local model representation of electronic potential energy surfaces; instead, the electron-nuclear interaction is treated explicitly using the Gaussian integration. As a result, the new scheme does not introduce any uncontrolled approximations. The employed variational principle ensures the energy conservation and leaves the number of electronic and nuclear basis functions as the only parameter determining the accuracy. To assess performance of the approach, a model with two electronic and two nuclear spacial degrees of freedom containing conical intersections between potential energy surfaces has been considered. Dynamical features associated with nonadiabatic transitions and nontrivial geometric (or Berry) phases were successfully reproduced within a limited basis expansion.

Using time-resolved photoelectron spectroscopy to unravel the electronic relaxation dynamics of photoexcited molecules

Time-resolved photoelectron spectroscopy measurements combined with quantum chemistry and dynamics calculations allow unprecedented insight into the electronic relaxation mechanisms of photoexcited molecules in the gas-phase. In this Tutorial Review, we explain the essential concepts linking photoelectron spectroscopy measurements with electronic structure and how key features on the potential energy landscape are identified using quantum chemistry and quantum dynamics calculations. We illustrate how time-resolved photoelectron spectroscopy and theory work together using examples ranging in complexity from the prototypical organic molecule benzene to a pyrrole dimer bound by a weak N-Hπ interaction and the green fluorescent protein chromophore.

Quantum modeling of ultrafast photoinduced charge separation

Phenomena involving electron transfer are ubiquitous in nature, photosynthesis and enzymes or protein activity being prominent examples. Their deep understanding thus represents a mandatory scientific goal. Moreover, controlling the separation of photogenerated charges is a crucial prerequisite in many applicative contexts, including quantum electronics, photo-electrochemical water splitting, photocatalytic dye degradation, and energy conversion. In particular, photoinduced charge separation is the pivotal step driving the storage of sun light into electrical or chemical energy. If properly mastered, these processes may also allow us to achieve a better command of information storage at the nanoscale, as required for the development of molecular electronics, optical switching, or quantum technologies, amongst others. In this Topical Review we survey recent progress in the understanding of ultrafast charge separation from photoexcited states. We report the state-of-the-art of the observation and theoretical description of charge separation phenomena in the ultrafast regime mainly focusing on molecular- and nano-sized solar energy conversion systems. In particular, we examine different proposed mechanisms driving ultrafast charge dynamics, with particular regard to the role of quantum coherence and electron-nuclear coupling, and link experimental observations to theoretical approaches based either on model Hamiltonians or on first principles simulations.

Curve crossing in a manifold of coupled electronic states: direct quantum dynamics simulations of formamide

Quantum dynamics simulations are an important tool to evaluate molecular behaviour including the, often key, quantum nature of the system. In this paper we present an algorithm that is able to simulate the time evolution of a molecule after photo-excitation into a manifold of states. The direct dynamics variational multi-configurational Gaussian (DD-vMCG) method circumvents the computational bottleneck problems of traditional grid-based methods by computing the potential energy functions on-the-fly, i.e. only where required. Unlike other commonly used direct dynamics methods, DD-vMCG is fully quantum mechanical. Here, the method is combined with a novel on-the-fly diabatisation scheme to simulate the short-time dynamics of the key molecule formamide and its acid analogue formimidic acid. This is a challenging test system due to the nature and large number of excited states, and eight coupled states are included in the calculations. It is shown that the method is able to provide unbiased information on the product channels open after excitation at different energies and demonstrates the potential to be a practical scheme, limited mainly by the quality of the quantum chemistry used to describe the excited states.

An ab initio multiple cloning approach for the simulation of photoinduced dynamics in conjugated molecules

Multidimensional wave function: a superposition of Gaussian coherent states guided by Ehrenfest trajectories suited to clone and swap their electronic amplitudes.

Kinetic isotope effects and how to describe them

We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.

Quantum Nonadiabatic Cloning of Entangled Coherent States

We propose a systematic approach to the basis set extension for nonadiabatic dynamics of entangled combination of nuclear coherent states (CSs) evolving according to the time-dependent variational principle (TDVP). The TDVP provides a rigorous framework for fully quantum nonadiabatic dynamics of closed systems; however, the quality of results strongly depends on available basis functions. Starting with a single nuclear CS replicated vertically on all electronic states, our approach clones this function when replicas of the CS on different electronic states experience increasingly different forces. Created clones move away from each other (decohere), extending the basis set. To determine a moment for cloning, we introduce generalized forces based on derivatives that maximally contribute to a variation of the total quantum action and thus account for entanglement of all basis functions.

Electron Dynamics upon Ionization of Polyatomic Molecules: Coupling to Quantum Nuclear Motion and Decoherence

Knowledge about the electronic motion in molecules is essential for our understanding of chemical reactions and biological processes. The advent of attosecond techniques opens up the possibility to induce electronic motion, observe it in real time, and potentially steer it. A fundamental question remains the factors influencing electronic decoherence and the role played by nuclear motion in this process. Here, we simulate the dynamics upon ionization of the polyatomic molecules paraxylene and modified bismethylene-adamantane, with a quantum mechanical treatment of both electron and nuclear dynamics using the direct dynamics variational multiconfigurational Gaussian method. Our simulations give new important physical insights about the expected decoherence process. We have shown that the decoherence of electron dynamics happens on the time scale of a few femtoseconds, with the interplay of different mechanisms: the dephasing is responsible for the fast decoherence while the nuclear overlap decay may actually help maintain it and is responsible for small revivals.

Accelerating direct quantum dynamics using graphical processing units

The direct dynamics variational multi-configurational Gaussian (DD-vMCG) method is combined with electronic structure calculations accelerated by Graphical Processing Units (GPUs). This is used to identify GPU acceleration will have a significant effect for both ground and excited state simulations.

Ab initio calculation of inelastic scattering

We calculate nonresonant inelastic electron and X-ray scattering cross sections for bound-to-bound transitions in atoms and molecules from ab initio electronic wavefunctions.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.