Intramolecular vibrational energy redistribution in toluene: a nine-dimensional quantum mechanical study using the MCTDH algorithm

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).

Femtochemistry of bimolecular reactions from weakly bound complexes: computational study of the H + H'OD → H'OH + D or HOD + H' exchange reactions

A full-dimensional wavepacket propagation describing the bimolecular exchange reactions H + H'OD → H'OH + D or HOD + H' initiated by photolysis of HCl in the hydrogen-bound complex (HCl)⋯(HOD) is reported. The dynamics of this reaction is carried out with the MCTDH method on an ab initio potential energy surface (PES) of H₃O and the initial state is derived from the ground state wavefunction of the complex obtained by relaxation on its own electronic ground state ab initio PES. The description of the system makes use of polyspherical coordinates parametrizing a set of Radau and Jacobi vectors. The calculated energy- and time-resolved reaction probabilities show, owing to the large collision energies at play stemming from the (almost full) photolysis of HCl, that the repulsion between oxygen in the H'OD molecule and the incoming hydrogen atom is the main feature of the collision and leads to non-reactive scattering. No abstraction reaction products are observed. However, both exchange processes are still observable, with a preference in O-H' bond dissociation over that of O-D. The selectivity is reversed upon vibrational pre-excitation of the O-D stretching mode in the H'OD molecule. It is shown that, after the collision, the hydrogen atom of HCl does most likely not encounter the almost stationary chlorine atom again but we also consider the limit case where the H atom is forced to collide multiple times against H'OD as a result of being pushed back by the Cl atom.

Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective

The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.

Rotationally inelastic scattering of O3–Ar: state-to-state rates with the multiconfigurational time dependent Hartree method

The rates of state-changing collisions are compared for different isotopologues of ozone from quantum scattering calculations with the MCTDH method.

The Role of Ozone Vibrational Resonances in the Isotope Exchange Reaction 16O16O + 18O → 18O16O + 16O: The Time-Dependent Picture

We consider the time-dependent dynamics of the isotope exchange reaction in collisions between an oxygen molecule and an oxygen atom: ¹⁶O¹⁶O + ¹⁸O → ¹⁶O¹⁸O + ¹⁶O. A theoretical approach using the multiconfiguration time-dependent Hartree method was employed to model the time evolution of the reaction. Two potential surfaces available in the literature were used in the calculations, and the results obtained with the two surfaces are compared with each other as well as with results of a previous theoretical time-independent approach. A good agreement for the reaction probabilities with the previous theoretical results is found. Comparing the results obtained using two potential energy surfaces allows us to understand the role of the reef/shoulder-like feature in the minimum energy path of the reaction in the isotope exchange process. Also, it was found that the distribution of final products of the reaction is highly anisotropic, which agrees with experimental observations and, at the same time, suggests that the family of approximated statistical approaches, assuming a randomized distribution over final exit channels, is not applicable to this case.

Quantum dynamics of the intramolecular vibrational energy redistribution in OCS: From localization to quasi-thermalization

We report a fully quantum-dynamical study of the intramolecular vibrational energy redistribution (IVR) in the electronic ground state of carbonyl sulfide, which is a prototype of an isolated many-body quantum system with strong internal couplings and non-Rice-Ramsperger-Kassel-Marcus (RRKM) behavior. We pay particular attention to the role of many-body localization and the approach to thermalization, which currently are topics of considerable interest, as they pertain to the very foundations of statistical mechanics and thermodynamics. We employ local-mode (valence) coordinates and consider initial excitations localized in one local mode, with energies ranging from low to near the dissociation threshold, where the classical dynamics have been shown to be chaotic. We propagate the nuclear wavepacket on the potential energy surface by means of the numerically exact multiconfiguration time-dependent Hartree method and employ mean local energies, time-dependent and time-averaged populations in quantum number space, energy distributions, entanglement entropies, local population distributions, microcanonical averages, and dissociation probabilities, as diagnostic tools. This allows us to identify a continuous localization → delocalization transition in the energy flow, associated with the onset of quantum chaos, as the excitation energy increases up to near the dissociation threshold. Moreover, we find that at this energy and ∼1 ps the molecule nearly thermalizes. Furthermore, we observe that IVR is so slow that the molecule begins to dissociate well before such quasi-thermalization is complete, in accordance with earlier classical-mechanical predictions of non-RRKM behavior.

ON THE MULTIDIMENSIONALITY OF INTRAMOLECULAR HYDROGEN BOND DYNAMICS: HYDROGEN TRANSFER AND IVR IN 3,7-DICHLOROTROPOLONE

In this work, we investigate the intramolecular vibrational energy redistribution associated with the hydrogen transfer in a derivative of tropolone, namely 3,7-dichlorotropolone. Our quantum simulation is based on the Cartesian reaction surface Hamiltonian together with the multi-configurational time-dependent Hartree approach for the wave-packet propagation. We compare results for two model systems with 6 and 14 dimensions, respectively. The 6D model accounts for the most strongly coupled modes, whereas the 14D model includes further modes with significantly weaker couplings. The linear absorption spectrum of both models shows the development of an OH-stretching band. Furthermore the results show that despite the fact, that the additional modes in the 14D system couple significantly weaker, there are qualitative differences in the decay behavior of an OH-stretching excitation. Limitations of the present reaction surface approach are also discussed.

Rovibrational energy transfer in ortho-H2+para-H2 collisions

We present the results of a full-dimensional quantum mechanical study of the rovibrational energy transfer in the collision between ortho-H2 and para-H2 in the energy range of 0.1-1.0 eV. The multiconfiguration time-dependent Hartree algorithm has been used to propagate the wave packets on the global potential energy surface by Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] and on a modification of this surface where the short range anisotropy is reduced. State-to-state attributes such as probabilities or integral cross sections are obtained using the formalism of Tannor and Weeks [J. Chem. Phys. 98, 3884 (1993)] by Fourier transforming the correlation functions. The effect of initial rotation of the diatoms on the inelastic and de-excitation processes is investigated.

Variational Calculation of Specific, Highly Excited Vibrational States in DFCO: Comparison with Experimental Data

A stimulated emission pumping spectra of jet-cooled DFCO performed by Crane et al. (J. Mol. Spectrosc. 1997, 183, 273) has provided a great number of ro-vibrational lines up to 9000 cm(-1) of excitation energy. By combining a Jacobi-Wilson (JW) approach with a Davidson scheme, we calculate the lines provided by the experiment up to 9000 cm(-1) using an ab initio global potential energy surface (PES) developed by Kato et al. (J. Chem. Phys. 1997, 107, 6114). Comparisons between experimental and calculated data provide a critical test of the quality of the PES used. We show that the variational calculated energies can be efficiently corrected by taking into account the error observed for the A' fundamental transitions nu(i) (i = 1, ..., 5) and the first overtone 2nu(6). A detailed analysis of the eigenstates obtained by the calculation allows one to quantify the coupling between the different modes. Such an information is essential to understand and predict the energy flow through a DFCO molecule that is initially excited.

Using neural networks to represent potential surfaces as sums of products

By using exponential activation functions with a neural network (NN) method we show that it is possible to fit potentials to a sum-of-products form. The sum-of-products form is desirable because it reduces the cost of doing the quadratures required for quantum dynamics calculations. It also greatly facilitates the use of the multiconfiguration time dependent Hartree method. Unlike potfit product representation algorithm, the new NN approach does not require using a grid of points. It also produces sum-of-products potentials with fewer terms. As the number of dimensions is increased, we expect the advantages of the exponential NN idea to become more significant.

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.

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.

Time-dependent perturbation theory for vibrational energy relaxation and dephasing in peptides and proteins

Without invoking the Markov approximation, we derive formulas for vibrational energy relaxation (VER) and dephasing for an anharmonic system oscillator using a time-dependent perturbation theory. The system-bath Hamiltonian contains more than the third order coupling terms since we take a normal mode picture as a zeroth order approximation. When we invoke the Markov approximation, our theory reduces to the Maradudin-Fein formula which is used to describe the VER properties of glass and proteins. When the system anharmonicity and the renormalization effect due to the environment vanishes, our formulas reduce to those derived by and Mikami and Okazaki [J. Chem. Phys. 121, 10052 (2004)] invoking the path-integral influence functional method with the second order cumulant expansion. We apply our formulas to VER of the amide I mode of a small amino-acid like molecule, N-methylacetamide, in heavy water.

Rotational excitation cross sections of para-H2+para-H2 collisions. A full-dimensional wave-packet propagation study using an exact form of the kinetic energy

A full-dimensional quantum dynamical study of the rotational excitation in para-para H2 + H2 collisions using the potential-energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] is reported. The multiconfiguration time-dependent Hartree algorithm is utilized to propagate wave packets and the cross sections for collision energies up to 1.2 eV are determined by a flux analysis through the interaction of the wave packet with a complex absorbing potential. Calculations for a collection of total angular momenta up to J = 70 are performed; the missing channels are obtained with a J-interpolation algorithm.

Wave packet study of the UV photodissociation of the Ar2HBr complex

The photodissociation dynamics of the Ar2HBr van der Waals molecule is studied using the multiconfiguration time-dependent Hartree method. Standard Jacobian coordinates are used to describe the molecule. Two four-dimensional calculations are carried out where the rotation of the Ar2 molecule and, in addition, either the vibration of Ar2-Br or that of Ar2 are frozen. The time-evolution of the probability density in the different modes and the calculation of the dissociative flux show that the dissociating hydrogen atom preferentially moves out of the plane defined by Ar2 and Br. A comprehensive study of the cage effect in the process is presented.

Multidimensional density operator propagations in open systems: Model studies on vibrational relaxations and surface sticking processes

An efficient method for the numerical treatment of multidimensional dynamics of open systems is presented: the multiconfiguration time-dependent Hartree (MCTDH) method extended to the propagation of density operators. With this method we investigate the relaxation process of a CO molecule adsorbed on a copper surface, i.e., CO/Cu(100), after the excitation with an infrared (IR) pulse. The interaction potential was taken from the literature. Lifetime estimations and thermalization studies were performed on the IR excited CO molecule. We were able to treat this system with all six degrees of freedom (DOF) and thus 12 dynamical variables, but most of our studies used a two or four DOF model. Finally, we demonstrate the applicability of MCTDH to the analysis of scattering processes in an open environment. We calculate sticking coefficients of a scattered particle to a model surface, the latter acting as heat bath. The surface corrugation and the initial particle energy have been varied, and six different relaxation strengths have been studied. These calculations were done under the inclusion of three DOFs: the two surface coordinates and the distance between the particle and the surface.