Efficient vibrationally correlated calculations using n-mode expansion-based kinetic energy operators.

Due to its efficiency and flexibility, the n-mode expansion is a frequently used tool for representing molecular potential energy surfaces in quantum chemical simulations. In this work, we investigate the performance of n-mode expansion-based models of kinetic energy operators in general polyspherical coordinate systems. In particular, we assess the operators with respect to accuracy in vibrationally correlated calculations and their effect on potential energy surface construction with the adaptive density guided approach. Our results show that the n-mode expansion-based operator variants are reliable and systematically improvable approximations of the full kinetic energy operator. Moreover, we introduce a workflow to generate the n-mode expanded kinetic energy operators on-the-fly within the adaptive density guided approach. This scheme can be applied in studies of species and coordinate systems, for which an analytical form of the kinetic energy operator is not available.

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics.

Through approximating electron-nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. A set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

Smolyak Scheme for solving the Schrödinger equation: Application to Malonaldehyde in Full Dimensionality.

In 1963 Smolyak introduced an approach to overcome the exponential scaling with respect to the number of variables of the direct product size [S. A. Smolyak Soviet Mathematics Doklady, 4, 240 (1963)]. The main idea is to replace a single large direct product by a sum of selected small direct products. It was first used in quantum dynamics in 2009 by Avila and Carrington [G. Avila and T. Carrington, J. Chem. Phys., 131, 174103 (2009)]. Since then, several calculations have been published by Avila and Carrington and by other groups. In the present study, and to push the limit to larger and more complex systems, this scheme is combined with the use of an on-the-fly calculation of the kinetic energy operator and a Block-Davidson procedure to obtain eigenstates in our home-made Fortran codes, ElVibRot and Tnum-Tana. This was applied to compute the tunneling splitting of malonaldehyde in full dimensionality (21D) using the potential of Mizukami et al. [W. Mizukami, S. Habershon, and D.P. Tew, J. Chem. Phys. 141, 1443-10 (2014)]. Our tunneling splitting calculations, 21.7±0.3 cm-1 and 2.9±0.1 cm-1 , show excellent agreement with the experimental values, 21.6 cm-1 and 2.9 cm-1 for the normal isotopologue and the mono-deuterated one, respectively.

Adiabatic and Nonadiabatic Dynamics with Interacting Quantum Trajectories.

We present a quantum dynamics method based on the propagation of interacting quantum trajectories to describe both adiabatic and nonadiabatic processes within the same formalism. The idea originates from the work of Poirier [Chem. Phys.2010,370, 4-14] and Schiff and Poirier [J. Chem. Phys.2012,136, 031102] on quantum dynamics without wavefunctions. It consists of determining the quantum force arising in the Bohmian hydrodynamic formulation of quantum dynamics using only information about quantum trajectories. The particular time-dependent propagation scheme proposed here results in very stable dynamics. Its performance is discussed by applying the method to analytical potentials in the adiabatic regime, and by combining it with the exact factorization method in the nonadiabatic regime.

The coupling of the hydrated proton to its first solvation shell.

The Zundel ([Formula: see text]) and Eigen ([Formula: see text]) cations play an important role as intermediate structures for proton transfer processes in liquid water. In the gas phase they exhibit radically different infrared (IR) spectra. The question arises: is there a least common denominator structure that explains the IR spectra of both, the Zundel and Eigen cations, and hence of the solvated proton? Full dimensional quantum simulations of these protonated cations demonstrate that two dynamical water molecules and an excess proton constitute this fundamental subunit. Embedded in the static environment of the parent Eigen cation, this subunit reproduces the positions and broadenings of its main excess-proton bands. In isolation, its spectrum reverts to the well-known Zundel ion. Hence, the dynamics of this subunit polarized by an environment suffice to explain the spectral signatures and anharmonic couplings of the solvated proton in its first solvation shell.

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.

Describing the photo-isomerization of a retinal chromophore model with coupled and quantum trajectories.

The exact factorization of the electron-nuclear wavefunction is applied to the study of photo-isomerization of a retinal chromophore model. We describe such an ultrafast nonadiabatic process by analyzing the time-dependent potentials of the theory and by mimicking nuclear dynamics with quantum and coupled trajectories. The time-dependent vector and scalar potentials are the signature of the exact factorization, as they guide nuclear dynamics by encoding the complete electronic dynamics and including excited-state effects. Analysis of the potentials is, thus, essential-when possible-to predict the time-dependent behavior of the system of interest. In this work, we employ the exact time-dependent potentials, available for the numerically exactly solvable model used here, to propagate quantum nuclear trajectories representing the isomerization reaction of the retinal chromophore. The quantum trajectories are the best possible trajectory-based description of the reaction when using the exact-factorization formalism and, thus, allow us to assess the performance of the coupled-trajectory, fully approximate schemes derived from the exact-factorization equations.

Smolyak Algorithm Adapted to a System-Bath Separation: Application to an Encapsulated Molecule with Large-Amplitude Motions.

A Smolyak algorithm adapted to system-bath separation is proposed for rigorous quantum simulations. This technique combines a sparse grid method with the system-bath concept in a specific configuration without limitations on the form of the Hamiltonian, thus achieving a highly efficient convergence of the excitation transitions for the "system" part. Our approach provides a general way to overcome the perennial convergence problem for the standard Smolyak algorithm and enables the simulation of floppy molecules with more than a hundred degrees of freedom. The efficiency of the present method is illustrated on the simulation of H2 caged in an sII clathrate hydrate including two kinds of cage modes. The transition energies are converged by increasing the number of normal modes of water molecules. Our results confirm the triplet splittings of both translational and rotational (j = 1) transitions of the H2 molecule. Furthermore, they show a slight increase of the translational transitions with respect to the ones in a rigid cage.

Quantum dynamics with curvilinear coordinates: models and kinetic energy operator.

In order to simplify the numerical solution of the time-dependent or time-independent Schrödinger equations associated with atomic and molecular motions, the use of well-adapted coordinates is essential. Usually, this set of curvilinear coordinates leads to a Hamiltonian operator that is as separable as possible. Although their corresponding kinetic energy operator (KEO) expressions can be derived analytically for small systems or special kinds of coordinates, a numerical and exact approach allows one to compute them in terms of sophisticated curvilinear coordinates. Furthermore, the numerical approach enables one to easily define reduced-dimensionality or constrained models. We present here a recent implementation of this numerical approach that allows nested coordinate transformations, therefore leading to great flexibility in the definition of the curvilinear coordinates. Furthermore, this implementation has no limitations in terms of numbers of atoms or coordinate transformations. The quantum dynamics of the cis-trans photoisomerization of part of the retinal chromophore illustrates the construction of the coordinates and KEO part of a three-dimensional model. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Electronic Structure and Excited States of the Collision Reaction O(3P) + C2H4: A Multiconfigurational Perspective.

We present a study of the O(3P) + C2H4 scattering reaction, a process that takes place in the interstellar medium and is of relevance in atmospheric chemistry as well. A comprehensive investigation of the electronic properties of the system has been carried out based on multiconfigurational ab initio CASSCF/CASPT2 calculations, using a robust and consistent active space that can deliver accurate potential energy surfaces in the key regions visited by the system. The paper discloses detailed description of the primary reaction pathways and the relevant singlet and triplet excited states at the CASSCF and CASPT2 level, including an accurate description of the critical configurations, such as minima and transition states. The chosen active space and the CASSCF/CASPT2 computational protocol are assessed against coupled-cluster calculations to further check the stability and reliability of the entire multiconfigurational procedure.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies.

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model.

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

Internal Conversion and Intersystem Crossing with the Exact Factorization.

We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.

First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel-Holstein Hamiltonian.

A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].

Vibrational Coupled Cluster Computations in Polyspherical Coordinates with the Exact Analytical Kinetic Energy Operator.

We present the first use of curvilinear vibrational coordinates, specifically polyspherical coordinates, in combination with vibrational coupled cluster theory. The polyspherical coordinates are used in the context of both the adaptive density-guided approach to potential energy surface construction and in the subsequent vibrational coupled cluster calculations of anharmonic vibrational states. Results obtained based on the polyspherical coordinate parametrization are compared to results obtained with the use of rectilinear vibrational coordinates, namely, normal coordinates and hybrid optimized and localized coordinates for the formaldehyde molecule. This comparison is carried out with the full vibrational configuration interaction model, using the respective fully coupled potential energy surfaces and untruncated kinetic energy operators. The polyspherical coordinates are shown to facilitate an acceleration of convergence for truncated methods when compared to the use of normal coordinates. We furthermore report on calculations on the hydrogen peroxide molecule in the polyspherical coordinate parametrization. The polyspherical vibrational coordinates are shown to perform very well, even for truncated methods, especially when considering the difficulty that rectilinear vibrational coordinates can exhibit in treating complicated internal molecular motion.

Spin-Orbit Interactions in Ultrafast Molecular Processes.

We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.

Intramolecular stretching vibrational states and frequency shifts of (H2)2 confined inside the large cage of clathrate hydrate from an eight-dimensional quantum treatment using small basis sets.

We report the results of calculations pertaining to the HH intramolecular stretching fundamentals of (p-H2)2 encapsulated in the large cage of structure II clathrate hydrate. The eight-dimensional (8D) quantum treatment assumes rotationless (j = 0) H2 moieties and a rigid clathrate structure but is otherwise fully coupled. The (H2)2-clathrate interaction is constructed in a pairwise-additive fashion, by combining the ab initio H2-H2O pair potential for flexible H2 and rigid H2O [D. Lauvergnat et al., J. Chem. Phys. 150, 154303 (2019)] and the six-dimensional (6D) H2-H2 potential energy surface [R. J. Hinde, J. Chem. Phys. 128, 154308 (2008)]. The calculations are performed by first solving for the eigenstates of a reduced-dimension 6D "intermolecular" Hamiltonian extracted from the full 8D Hamiltonian by taking the H2 moieties to be rigid. An 8D contracted product basis for the solution of the full problem is then constructed from a small number of the lowest-energy 6D intermolecular eigenstates and two discrete variable representations covering the H2-monomer internuclear distances. Converged results are obtained already by including just the two lowest intermolecular eigenstates in the final 8D basis of dimension 128. The two HH vibrational stretching fundamentals are computed for three hydrate domains having an increasing number of H2O molecules. For the largest domain, the two fundamentals are found to be site-split by ∼0.5 cm-1 and to be redshifted by about 24 cm-1 from the free-H2 monomer stretch frequency, in excellent agreement with the experimental value of 26 cm-1. A first-order perturbation theory treatment gives results that are nearly identical to those of the 8D quantum calculations.

A generalized vibronic-coupling Hamiltonian for molecules without symmetry: Application to the photoisomerization of benzopyran.

We present a model for the lowest two potential energy surfaces (PESs) that describe the photoinduced ring-opening reaction of benzopyran taken as a model compound to study the photochromic ring-opening reaction of indolinobenzospiropyran and its evolution toward its open-chain analog. The PESs are expressed in terms of three effective rectilinear coordinates. One corresponds to the direction between the equilibrium geometry in the electronic ground state, referred to as the Franck-Condon geometry, and the minimum of conical intersection (CI), while the other two span the two-dimensional branching space at the CI. The model correctly reproduces the topography of the PESs. The ab initio calculations are performed with the extended multiconfiguration quasidegenerate perturbation theory at second order method. We demonstrate that accounting for electron dynamic correlation drastically changes the global energy landscape since some zwitterionic states become strongly stabilized. Quantum dynamics calculations using this PES model produce an absorption spectrum that matches the experimental one to a good accuracy.

H2, HD, and D2 in the small cage of structure II clathrate hydrate: Vibrational frequency shifts from fully coupled quantum six-dimensional calculations of the vibration-translation-rotation eigenstates.

We report the first fully coupled quantum six-dimensional (6D) bound-state calculations of the vibration-translation-rotation eigenstates of a flexible H2, HD, and D2 molecule confined inside the small cage of the structure II clathrate hydrate embedded in larger hydrate domains with up to 76 H2O molecules, treated as rigid. Our calculations use a pairwise-additive 6D intermolecular potential energy surface for H2 in the hydrate domain, based on an ab initio 6D H2-H2O pair potential for flexible H2 and rigid H2O. They extend to the first excited (v = 1) vibrational state of H2, along with two isotopologues, providing a direct computation of vibrational frequency shifts. We show that obtaining a converged v = 1 vibrational state of the caged molecule does not require converging the very large number of intermolecular translation-rotation states belonging to the v = 0 manifold up to the energy of the intramolecular stretch fundamental (≈4100 cm-1 for H2). Only a relatively modest-size basis for the intermolecular degrees of freedom is needed to accurately describe the vibrational averaging over the delocalized wave function of the quantum ground state of the system. For the caged H2, our computed fundamental translational excitations, rotational j = 0 → 1 transitions, and frequency shifts of the stretch fundamental are in excellent agreement with recent quantum 5D (rigid H2) results [A. Powers et al., J. Chem. Phys. 148, 144304 (2018)]. Our computed frequency shift of -43 cm-1 for H2 is only 14% away from the experimental value at 20 K.

Does cage quantum delocalisation influence the translation-rotational bound states of molecular hydrogen in clathrate hydrate?

In this study, we examine the effect of a flexible description of the clathrate hydrate framework on the translation-rotation (TR) eigenstates of guest molecules such as molecular hydrogen. Traditionally, the water cage structure is assumed to be rigid, thus ignoring the quantum nature of hydrogen nuclei in the water framework. However, it has been shown that protons in a water molecule possess a marked delocalised character in many situations, ranging from water clusters to proton transfer in the bulk. In the case of water clathrates, all previous TR bound-state calculations of guest molecules consider that the caging water molecules are fixed at their equilibrium geometry. Only recently, a static investigation of the role of proton configurations was performed by Bacic and co-workers by sampling a very large number of different static structures of water clathrates. Here, we investigate the importance of the rotational degrees of freedom of the water cage on the TR levels of the guest molecule using an efficient adiabatic decoupling scheme. Our approach combines rigid body diffusion Monte Carlo calculations for the description of the rotational degree of freedom of water molecules surrounding the guest molecular hydrogen to an efficient Smolyak sparse-grid technique for the calculation of the TR levels. This approach allows us to take into account the highly anharmonic nature of the rotational water motions in a high-dimensional system. The clathrate-induced splittings of the j = 1 rotational levels are much more sensitive to the quantum hydrogen delocalisation than the translational transitions. This result is in good agreement with the previous static study of Bacic and co-workers.