Dynamical Tunneling in More than Two Degrees of Freedom

Recent progress towards understanding the mechanism of dynamical tunneling in Hamiltonian systems with three or more degrees of freedom (DoF) is reviewed. In contrast to systems with two degrees of freedom, the three or more degrees of freedom case presents several challenges. Specifically, in higher-dimensional phase spaces, multiple mechanisms for classical transport have significant implications for the evolution of initial quantum states. In this review, the importance of features on the Arnold web, a signature of systems with three or more DoF, to the mechanism of resonance-assisted tunneling is illustrated using select examples. These examples represent relevant models for phenomena such as intramolecular vibrational energy redistribution in isolated molecules and the dynamics of Bose-Einstein condensates trapped in optical lattices.

E-Z Isomerization in Guanidine: Second-order Saddle Dynamics, Non-statisticality, and Time-frequency Analysis

Our recent work on the E-Z isomerization reaction of guanidine using ab initio chemical dynamics simulations [Rashmi et al., Regul. Chaotic Dyn. 2021, 26, 119] emphasized the role of second-order saddle (SOS) in the isomerization reaction; however, we could not unequivocally establish the non-statistical nature of the dynamics followed in the reaction. In the present study, we performed thousands of on-the-fly trajectories using forces computed at the MNDO level to investigate the influence of second-order saddle in the E-Z isomerization reaction of guanidine and the role of intramolecular vibrational energy redistribution (IVR) on the reaction dynamics. The simulations reveal that while majority of the trajectories follow the traditional transition state pathways, 15 % of the trajectories follow the SOS path. The dynamics was found to be highly non-statistical with the survival probabilities of the reactants showing large deviations from those obtained within the RRKM assumptions. In addition, a detailed analysis of the dynamics using time-dependent frequencies and the frequency ratio spaces reveal the existence of multiple resonance junctions that indicate the existence of regular dynamics and long-lived quasi-periodic trajectories in the phase space associated with non-RRKM behavior.

Electric Fields Influence Intramolecular Vibrational Energy Relaxation and Line Widths

Intramolecular vibrational energy relaxation (IVR) is fundamentally important to chemical dynamics. We show that externally applied electric fields affect IVR and vibrational line widths by changing the anharmonic couplings and frequency detunings between modes. We demonstrate this effect in benzonitrile for which prior experimental results show a decrease in vibrational line width as a function of applied electric field. We identify three major channels for IVR that depend on electric field. In the dominant channel, the electric field affects the frequency detuning, while in the other two channels, variation of anharmonic couplings as a function of field is the underlying mechanism. Consistent with experimental results, we show that the combination of all channels gives rise to reduced line widths with increasing electric field in benzonitrile. Our results are relevant for controlling IVR with external or internal fields and for gaining a more complete interpretation of line widths of vibrational Stark probes.

Sub-picosecond to Sub-nanosecond Vibrational Energy Transfer Dynamics in Pentaerythritol Tetranitrate

The time scale associated with shock-induced detonation is a key property of energetic materials that remains poorly understood. Herein, we test aspects of one potential mechanism, the phonon up-pumping mechanism, where shock compression excites lattice phonon modes, transferring energy to intramolecular vibrations leading to chemical bond cleavage and reaction. Using ultrafast infrared pump-probe spectroscopy on pentaerythritol tetranitrate (PETN), we reveal sub-picosecond vibrational energy transfer (VET) from the photoexcited band at 1660 cm^-1 into every other infrared-active mode in the probed frequency range 800-1800 cm^-1. Energy transfer processes remain incomplete at 150 ps. Computational predictions from density functional theory are used in tandem to elucidate VET pathways in PETN.

Stable chaos and delayed onset of statisticality in unimolecular dissociation reactions

Statistical models provide a powerful and useful class of approximations for calculating reaction rates by bypassing the need for detailed, and often difficult, dynamical considerations. Such approaches invariably invoke specific assumptions about the extent of intramolecular vibrational energy flow in the system. However, the nature of the transition to the statistical regime as a function of the molecular parameters is far from being completely understood. Here, we use tools from nonlinear dynamics to study the transition to statisticality in a model unimolecular reaction by explicitly visualizing the high dimensional classical phase space. We identify generic features in the phase space involving the intersection of two or more independent anharmonic resonances and show that the presence of correlated, but chaotic, intramolecular dynamics near such junctions leads to nonstatisticality. Interestingly, akin to the stability of asteroids in the Solar System, molecules can stay protected from dissociation at the junctions for several picoseconds due to the phenomenon of stable chaos.

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.

Relevance of the Resonance Junctions on the Arnold Web to Dynamical Tunneling and Eigenstate Delocalization

We study the competition and correspondence between the classical and quantum routes to intramolecular vibrational energy redistribution (IVR) in a three degrees of freedom model effective Hamiltonian. Specifically, we focus on the classical and the quantum dynamics near the resonance junctions on the Arnold web that are formed by an intersection of independent resonances. The regime of interest models the IVR dynamics from highly excited initial states near dissociation thresholds of molecular systems wherein both classical and purely quantum, involving dynamical tunneling, routes to IVR coexist. In the vicinity of a resonance junction, classical chaos is inevitably present, and hence one expects the quantum IVR pathways to have a strong classical component as well. We show that with increasing resonant coupling strengths the classical component of IVR leads to a transition from coherent dynamical tunneling to incoherent dynamical tunneling. Furthermore, we establish that the quantum IVR dynamics can be predicted based on the structures on the classical Arnold web. In addition, we investigate the nature of the highly excited eigenstates to identify the quantum signatures of the multiplicity-2 junctions. For the parameter regimes studies herein, by projecting the eigenstates onto the Arnold web, we find that eigenstates in the vicinity of the junctions are primarily delocalized due to dynamical tunneling.

Molecules and the Eigenstate Thermalization Hypothesis

We review a theory that predicts the onset of thermalization in a quantum mechanical coupled non-linear oscillator system, which models the vibrational degrees of freedom of a molecule. A system of N non-linear oscillators perturbed by cubic anharmonic interactions exhibits a many-body localization (MBL) transition in the vibrational state space (VSS) of the molecule. This transition can occur at rather high energy in a sizable molecule because the density of states coupled by cubic anharmonic terms scales as N³, in marked contrast to the total density of states, which scales as exp(aN), where a is a constant. The emergence of a MBL transition in the VSS is seen by analysis of a random matrix ensemble that captures the locality of coupling in the VSS, referred to as local random matrix theory (LRMT). Upon introducing higher order anharmonicity, the location of the MBL transition of even a sizable molecule, such as an organic molecule with tens of atoms, still lies at an energy that may exceed the energy to surmount a barrier to reaction, such as a barrier to conformational change. Illustrative calculations are provided, and some recent work on the influence of thermalization on thermal conduction in molecular junctions is also discussed.