Time-dependent wave packet study on trans-cis isomerization of HONO

Analytical high-dimensional operators in canonical polyadic finite basis representation (CP-FBR)

In the present work, we introduce a simple means of obtaining an analytical (i.e., grid-free) canonical polyadic (CP) representation of a multidimensional function that is expressed in terms of a set of discrete data. For this, we make use of an initial CP guess, even not fully converged, and a set of auxiliary basis functions [finite basis representation (FBR)]. The resulting CP-FBR expression constitutes the CP counterpart of our previous Tucker sum-of-products-FBR approach. However, as is well-known, CP expressions are much more compact. This has obvious advantages in high-dimensional quantum dynamics. The power of CP-FBR lies in the fact that it requires a grid much coarser than the one needed for the dynamics. In a subsequent step, the basis functions can be interpolated to any desired density of grid points. This is useful, for instance, when different initial conditions (e.g., energy content) of a system are to be considered. We show the application of the method to bound systems of increased dimensionality: H₂ (3D), HONO (6D), and CH₄ (9D).

Intramolecular Vibrational Redistribution in Formic Acid and its Deuterated Forms

The intramolecular vibrational relaxation dynamics of formic acid and its deuterated isotopologues is simulated on the full-dimensional potential energy surface of Richter and Carbonnière [F. Richter and P. Carbonnière, J. Chem. Phys. 148, 064303 (2018)] using the Heidelberg MCTDH package. We focus on couplings with the torsion vibrational modes close to the trans- cis isomerisation coordinate from the dynamics of artificially excited vibrational mode overtones. The C-O stretch bright vibrational mode is coupled to the out-of-the plane torsion mode in HCOOH, where this coupling could be exploited for laser-induced trans-to- cis isomerisation. Strong isotopic effects are observed: deuteration of the hydroxyl group, i.e., in HCOOD and DCOOD, destroys the C-O stretch to torsion mode coupling whereas in DCOOH, little to no effect is observed.

Suppression and Enhancement of Thermal Chemical Rates in a Cavity

The observed modification of thermal chemical rates in Fabry-Perot cavities remains a poorly understood effect theoretically. Recent breakthroughs explain some of the observations through the Grote-Hynes theory, where the cavity introduces friction with the reaction coordinate, thus reducing the transmission coefficient and the rate. The regime of rate enhancement, the observed sharp resonances at varying cavity frequencies, and the survival of these effects in the collective regime remain mostly unexplained. In this Letter, we consider the cis-trans isomerization of HONO atomistically using an ab initio potential energy surface. We evaluate the transmission coefficient using the reactive flux method and identify the conditions for rate acceleration. In the underdamped, low-friction regime of the reaction coordinate, the cavity coupling enhances the rate with increasing coupling strength until reaching the Kramers turnover point. Sharp resonances in this regime are related to cavity-enabled energy redistribution channels.

Analyzing Grid-Based Direct Quantum Molecular Dynamics Using Non-Linear Dimensionality Reduction

Grid-based schemes for simulating quantum dynamics, such as the multi-configuration time-dependent Hartree (MCTDH) method, provide highly accurate predictions of the coupled nuclear and electronic dynamics in molecular systems. Such approaches provide a multi-dimensional, time-dependent view of the system wavefunction represented on a coordinate grid; in the case of non-adiabatic simulations, additional information about the state populations adds a further layer of complexity. As such, wavepacket motion on potential energy surfaces which couple many nuclear and electronic degrees-of-freedom can be extremely challenging to analyse in order to extract physical insight beyond the usual expectation-value picture. Here, we show that non-linear dimensionality reduction (NLDR) methods, notably diffusion maps, can be adapted to extract information from grid-based wavefunction dynamics simulations, providing insight into key nuclear motions which explain the observed dynamics. This approach is demonstrated for 2-D and 9-D models of proton transfer in salicylaldimine, as well as 8-D and full 12-D simulations of cis-trans isomerization in ethene; these simulations demonstrate how NLDR can provide alternative views of wavefunction dynamics, and also highlight future developments.

Ab initio calculation of rovibrational states for non-degenerate double-well potentials: cis-trans isomerization of HOPO

The rovibrational spectra of metaphosphorous acid, HOPO, and its deuterated isotopologue have been studied by vibrational configuration interaction calculations, relying on the internal coordinate path Hamiltonian and the Watson Hamiltonian. Tunneling effects for the overtones of the torsional mode, which gives rise to the cis-trans isomerization, and its rovibrational transitions have been investigated in detail. Due to strong matrix effects, comparison with experimental data is hindered, and thus, the calculations provide accurate estimates for the fundamental modes of these species.

Fitting potential energy surfaces to sum-of-products form with neural networks using exponential neurons

In this paper, the use of the neural network (NN) method with exponential neurons for directly fitting ab initio data to generate potential energy surfaces (PESs) in sum-of-product form will be discussed. The utility of the approach will be highlighted using fits of CS2, HFCO, and HONO ground state PESs based upon high-level ab initio data. Using a generic interface between the neural network PES fitting, which is performed in MATLAB, and the Heidelberg multi-configuration time-dependent Hartree (MCTDH) software package, the PESs have been tested via comparison of vibrational energies to experimental measurements. The review demonstrates the potential of the PES fitting method, combined with MCTDH, to tackle high-dimensional quantum dynamics problems.

A ground state potential energy surface for HONO based on a neural network with exponential fitting functions

Using CCSD(T)-F12/cc-pVTZ-F12 and CCSD(T)/CBS ab initio energies, two different six-dimensional ground state potential energy surfaces for HONO have been fit in sum-of-products form using neural network exponential fitting functions and tested by computing vibrational energies with MCTDH.

Quantum dynamics study of fulvene double bond photoisomerization: the role of intramolecular vibrational energy redistribution and excitation energy

The double bond photoisomerization of fulvene has been studied with quantum dynamics calculations using the multi-configuration time-dependent Hartree method. Fulvene is a test case to develop optical control strategies based on the knowledge of the excited state decay mechanism. The decay takes place on a time scale of several hundred femtoseconds, and the potential energy surface is centered around a conical intersection seam between the ground and excited state. The competition between unreactive decay and photoisomerization depends on the region of the seam accessed during the decay. The dynamics are carried out on a four-dimensional model surface, parametrized from complete active space self-consistent field calculations, that captures the main features of the seam (energy and locus of the seam and associated branching space vectors). Wave packet propagations initiated by single laser pulses of 5-25 fs duration and 1.85-4 eV excitation energy show the principal characteristics of the first 150 fs of the photodynamics. Initially, the excitation energy is transferred to a bond stretching mode that leads the wave packet to the seam, inducing the regeneration of the reactant. The photoisomerization starts after the vibrational energy has flowed from the bond stretching to the torsional mode. In our propagations, intramolecular energy redistribution (IVR) is accelerated for higher excess energies along the bond stretch mode. Thus, the competition between unreactive decay and isomerization depends on the rate of IVR between the bond stretch and torsion coordinates, which in turn depends on the excitation energy. These results set the ground for the development of future optical control strategies.

Multilayer multiconfiguration time-dependent Hartree method: implementation and applications to a Henon-Heiles hamiltonian and to pyrazine

The multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method is discussed and a fully general implementation for any number of layers based on the recursive ML-MCTDH algorithm given by Manthe [J. Chem. Phys. 128, 164116 (2008)] is presented. The method is applied first to a generalized Henon-Heiles (HH) hamiltonian. For 6D HH the overhead of ML-MCTDH makes the method slower than MCTDH, but for 18D HH ML-MCTDH starts to be competitive. We report as well 1458D simulations of the HH hamiltonian using a seven-layer scheme. The photoabsorption spectrum of pyrazine computed with the 24D hamiltonian of Raab et al. [J. Chem. Phys. 110, 936 (1999)] provides a realistic molecular test case for the method. Quick and small ML-MCTDH calculations needing a fraction of the time and resources of reference MCTDH calculations provide already spectra with all the correct features. Accepting slightly larger deviations, the calculation can be accelerated to take only 7 min. When pushing the method toward convergence, results of similar quality than the best available MCTDH benchmark, which is based on a wavepacket with 4.6×10(7)time-dependent coefficients, are obtained with a much more compact wavefunction consisting of only 4.5×10(5) coefficients and requiring a shorter computation time.

Photochemical processes induced by vibrational overtone excitations: dynamics simulations for cis-HONO, trans-HONO, HNO3, and HNO3-H2O

Photochemical processes in HNO3, HNO3-H2O, and cis- and trans-HONO following overtone excitation of the OH stretching mode are studied by classical trajectory simulations. Initial conditions for the trajectories are sampled according to the initially prepared vibrational wave function. Semiempirical potential energy surfaces are used in "on-the-fly" simulations. Several tests indicate at least semiquantitative validity of the potential surfaces employed. A number of interesting new processes and intermediate species are found. The main results include the following: (1) In excitation of HNO3 to the fifth and sixth OH-stretch overtone, hopping of the H atom between the oxygen atoms is found to take place in nearly all trajectories, and can persist for many picoseconds. H-atom hopping events have a higher yield and a faster time scale than the photodissociation of HNO3 into OH and NO2. (2) A fraction of the trajectories for HNO3 show isomerization into HOONO, which in a few cases dissociates into HOO and NO. (3) For high overtone excitation of HONO, isomerization into the weakly bound species HOON is seen in all trajectories, in part of the events as an intermediate step on the way to dissociation into OH + NO. This process has not been reported previously. Well-established processes for HONO, including cis-trans isomerization and H hopping are also observed. (4) Only low overtone levels of HNO3-H2O have sufficiently long liftimes to be spectrocopically relevant. Excitation of these OH stretching overtones is found to result in the dissociation of the cluster H hopping, or dissociation of HNO3 does not take place. The results demonstrate the richness of processes induced by overtone excitation of HNO(x) species, with evidence for new phenomena. Possible relevance of the results to atmospheric processes is discussed.

Theoretical investigation of highly excited vibrational states in DFCO: Calculation of the out-of-plane bending states and simulation of the intramolecular vibrational energy redistribution

A previously developed modified Davidson scheme [C. Iung and F. Ribeiro, J. Chem. Phys. 121, 174105 (2005)] is applied to compute and analyze highly excited (nu2,nu6) eigenstates in DFCO. The present paper is also devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) initiated by an excitation of the out-of-plane bending vibration (nnu6, n=2,4,6, . . . ,18, and 20). The multiconfiguration time-dependent Hartree method is exploited to propagate the corresponding six-dimensional wave packets. A comprehensive comparison with experimental data as well as with previous simulations of IVR in HFCO [G. Pasin et al. J. Chem. Phys. 124, 194304 (2006)] is presented.

Direct quantum dynamics using variational multi-configuration Gaussian wavepackets. Implementation details and test case

We present here a direct quantum dynamics method using variational multi-configuration Gaussian wavepackets. Based on the efficient multi-configuration time-dependent Hartree wavepacket propagation algorithm, it uses on-the-fly quantum chemical calculation of the potential energy and its derivatives rather than fitted surfaces. Intermediate results are stored in a database so that expensive quantum chemical computations can be recycled. This method is intended to treat quantum effects in the photochemistry of large molecules and the use of Cartesian coordinates to perform direct dynamics is discussed with a comparison between Cartesian coordinates of Jacobi vectors and Cartesian coordinates of the nuclei, using various free and constrained approaches depending on the way the rotation is treated. As a test calculation to be compared to full quantum dynamics it is applied here to the computation of the photodissociation spectrum of nitrosyl chloride (NOCl).

The infrared-driven cis-trans isomerization of HONO. II: Vibrational relaxation and slow isomerization channel

In a recent paper [R. Schanz et al., J. Chem. Phys. 122, 044509 (2005)], we investigated the IR-driven cis-trans isomerization of HONO in a Kr matrix with the help of femtosecond IR spectroscopy. We found that isomerization occurs on a 20 ps time scale, however, with a cis-->trans quantum yield of only 10% that is significantly below the value reported in the literature (close to 100%). At the same time, we concluded that vibrational energy has not completely dissipated out of the molecule at the maximum delay time we reached in this study (500 ps). In order to verify whether additional, slower reaction channels exist, we extend the study here to delay times up to 100 ns. At a temperature of 32 K, we indeed find an additional isomerization channel on a 2 ns timescale, which increases the total cis-->trans quantum yield to approximately 30%. The trans-->cis quantum yield is approximately 7%. There is still a discrepancy between the quantum yields we observe and the literature values, however, we provide experimental evidence that this discrepancy is due to the different temperatures of our study. Vibrational cooling occurs on a 20 ns time scale, and cascades in a highly nonstatistical manner through one single normal mode (most likely the ONO bending mode nu(5)). Intermolecular energy dissipation into the rare gas matrix is more efficient than intramolecular vibrational energy redistribution and the matrix environment can certainly not be considered a weak perturbation.

Theoretical investigation of intramolecular vibrational energy redistribution in highly excited HFCO

The present paper is devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) in HFCO initiated by an excitation of the out-of-plane bending vibration [nnu(6)=2,4,6,...,18,20]. Using a full six-dimensional ab initio potential energy, the multiconfiguration time-dependent Hartree (MCTDH) method was exploited to propagate the corresponding six-dimensional wave packets. This study emphasizes the stability of highly excited states of the out-of-plane bending mode which exist even above the dissociation threshold. More strikingly, the structure of the IVR during the first step of the dynamics is very stable for initial excitations ranging from 2nu(6) to 20nu(6). This latter result is consistent with the analysis of the eigenstates obtained, up to 10nu(6), with the aid of the Davidson algorithm in a foregoing paper [Iung and Ribeiro, J. Chem. Phys. 121, 174105 (2005)]. The present study can be considered as complementary to this previous investigation. This paper also shows how MCTDH can be used to predict the dynamical behavior of a strongly excited system and to determine the energies of the corresponding highly excited states.

A femtosecond study of the infrared-driven cis-trans isomerization of nitrous acid (HONO)

We investigate the dynamics and mechanism of the IR-driven cis-trans isomerization of nitrous acid (HONO) in a low-temperature krypton matrix applying ultrafast time resolved IR spectroscopy. After excitation of the OH-stretching mode the trans HONO state decays biexponentially on a 8 and 260 ps time scale. The initially excited cis HONO state decays on a 20 ps time scale. Cis HONO isomerizes with 10% quantum yield on a 20 ps time scale to trans HONO. The quantum yield we observe is significantly smaller than the previously reported 100%, which could imply that additional, much slower reaction channels exist. We furthermore developed a four-dimensional model of the system, which includes the three proton intramolecular degrees of freedom of HONO fully quantum mechanically and one intermolecular translational degree of freedom of the molecule in the crystal cage. We find that cis-trans isomerization necessarily is accompanied by a translation of the molecule as a whole in the crystal cage. The translational degree of freedom tunes the intramolecular proton states of HONO with respect to each other. When resonances occur, the proton states might couple and transfer population. We suggest a possible reaction pathway, where the cis OH-stretch excited state first couples to a high cis torsional mode, which then may transfer almost instantaneously to the trans side. The model qualitatively explains all experimental observations.

Study of vibrational energy localization and redistribution in hydrogen peroxide H2O2 at low energy

Vibrational energy localization and/or redistribution in hydrogen peroxide H2O2 is studied at about 4000 cm(-1) above the quantum mechanical ground state using the ab initio potential energy surface of Koput, Carter, and Handy [J. Phys. Chem. A 102, 6325 (1998)]. In this work, the recently derived canonical perturbation procedure for floppy molecules serves two purposes. First, from the quantum mechanical point of view, it is shown that the energies of the lowest 130 states are reproduced with an average error smaller than 1.5 cm(-1) by a two-dimensional Hamiltonian, which is a function of the torsion and OO-stretch coordinates and momenta, while the other four degrees of freedom contribute only through powers of good quantum numbers. Moreover, the canonical perturbation procedure is also used in classical mechanics calculations, in order to define meaningful local modes, for which the ingoing and outgoing energy flows are monitored. Almost all the individual trajectories launched on the ab initio surface display the same behavior, that is, the superposition of (a) rapid (few hundreds of femtoseconds) and quasiperiodic energy exchanges between the two OH stretches and between the torsion and OO-stretch, and (b) slower (few to several picoseconds) but erratic-looking energy flows between all degrees of freedom. When averaging over large numbers of trajectories with the same local mode energies at time t=0, one observes instead a smooth and irreversible energy flow between all degrees of freedom, which usually thermalize in the range of several tens of picoseconds, that is, on time scales larger than the 5 ps period associated with the quantum density of states.

Cumulative isomerization probability studied by various transition state wave packet methods including the MCTDH algorithm. Benchmark: HCN→CNH isomerization

The 3D cumulative isomerization probability N(E) for the transfer of a light particle between two atoms is computed by one time-independent and two time-dependent versions of the transition state wave packet (TSWP) method. The time-independent method is based on the direct expansion of the microcanonical projection operator on Chebyshev polynomials. In the time-dependent TSWP methods, the propagations are carried out by the split operator scheme and the multiconfiguration time-dependent Hartree (MCTDH) algorithm. This is the very first implementation of the TSWP method in the Heidelberg MCTDH package [G. W. Worth, M. H. Beck, A. Jackle, and H.-D. Meyer, The MCDTH package, Version 8.2 (2000); H.-D Meyer, Version 8.3 (2002). See http://www.pci.uni-heidelberg.de/tc/usr/mctdh/]. The benchmark is the HCN-->CNH isomerization for zero total angular momentum. Particular insights are given into the tunneling region. In larger systems, the time-dependent version of TSWP making use of the MCTDH algorithm will permit to treat more and more modes quantum mechanically, for very accurate results. Therefore, it was important to calibrate the implementation. Besides, we also assess the efficiency of a reduced dimensionality approach by comparing the new exact 3D calculations of N(E) for the HCN-->CNH isomerization with results obtained via 1D or 2D active subspaces. This suggests that, it should be possible to take directly benefit of the present 3D approaches, adapted for triatomic Jacobi coordinates to compute N(E) for H-transfer in larger systems, via three active coordinates. The prerequisite is then the simplification of the reduced 3D kinetic energy operator with rigid constraint to take the form corresponding to a pseudo triatomic system in Jacobi coordinates with effective masses. This last step is checked in the methoxy radical and malonaldehyde. Finally, different ways to obtain reliable eigenvectors of the flux operator associated with a dividing surface are revisited.

Intramolecular vibrational energy redistribution in the highly excited fluoroform molecule: A quantum mechanical study using the multiconfiguration time-dependent Hartree algorithm

The present paper is devoted to a detailed study of the intramolecular vibrational energy redistribution in fluoroform initiated by a local mode excitation of the CH stretch [nnu(CH) (n=1,...,4)]. All nine internal degrees of freedom are explicitly taken into account and the full quantum mechanical simulation is performed by means of the multiconfiguration time-dependent Hartree algorithm. The existence of different time scales considerably complicates the dynamics. The mode-to-mode energy transfer is analyzed by calculating the evolution of the partial energies of all vibrational modes. This study emphasizes the crucial role played by the two-dimensional FCH bending modes which act as an energy reservoir. The fast energy flow into these bending modes significantly hinders an energy flow from the CH chromophore. Finally, our results are compared with those obtained previously with the wave operator sorting algorithm approach.