Assignment and Extraction of Dynamics of a Small Molecule with a Complex Vibrational Spectrum: Thiophosgene

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.

Teaching vibrational spectra to assign themselves

A new paradigm for assigning vibrational spectra is described. Instead of proceeding from potential to Hamiltonian to eigenvalues/eigenvectors/intensities to spectrum, the new method goes "backwards" directly from spectrum to eigenvectors. The eigenvectors then "assign" the spectrum, in that they identify the basis states that contribute to each eigenvalue. To start, we demonstrate an algorithm that can obtain useful estimates of the eigenvectors connecting a real, symmetric Hamiltonian to its eigenvalues even if the only available information about the Hamiltonian is its diagonal elements. When this algorithm is augmented with information about transition intensities, it can be used to assign a complex vibrational spectrum using only information about (1) eigenvalues (the peak centers of the spectrum) and (2) a harmonic basis set (taken to be the diagonal elements of the Hamiltonian). Examples will be discussed, including application to the vibrationally complex spectral region of the formic acid dimer.

Assigning the low lying vibronic states of CH3O and CD3O

The assignment of lines in vibrational spectra in strongly mixing systems is considered. Several low lying vibrational states of the ground electronic X∼²E state of the CH₃O and CD₃O radicals are assigned. Jahn-Teller, spin-orbit, and Fermi couplings mix the normal mode states. The mixing complicates the assignment of the infrared spectra using a zero-order normal mode representation. Alternative zero-order representations, which include specific Jahn-Teller couplings, are explored. These representations allow for definitive assignments. In many instances it is possible to plot the wavefunctions on which the assignments are based. The plots, which are shown in the adiabatic representation, allow one to visualize the effects of various higher order couplings. The plots also enable one to visualize the conical seam and its effect on the wavefunctions. The first and the second order Jahn-Teller couplings in the rocking motion dominate the spectral features in CH₃O, while first order and modulated first order couplings dominate the spectral features in CD₃O. The methods described here are general and can be applied to other Jahn-Teller systems.

Large scale vibrational calculations on IVR in S0 thiophosgene

In this work, using a recently derived refined potential energy surface for [Formula: see text] thiophosgene, we perform large scale vibrational calculations to explore the IVR characteristics and vibrational mixing at very high vibrational excitation energies, ranging up to the dissociation limit (at [Formula: see text]20,000[Formula: see text]cm[Formula: see text]). The results from our calculations have been compared to the conclusions based on the available experimentally measured dataset (obtained from SEP and LIF spectra) as well as to the conclusions from the analyses by other authors using local coupling models.

A refined quartic potential energy surface and large scale vibrational calculations for S0 thiophosgene

In this work we present a full 6D quartic potential energy surface (PES) for S0 thiophosgene in curvilinear symmetrized bond-angle coordinates. The PES was refined starting from an ab initio field derived from acc-pVTZ basis set with CCSD(T) corrections for electron correlation. In the present calculations we used our variational method that was recently tested on formaldehyde and some of its isotopomers, along with additional improvements. The lower experimentally known vibrational levels for 35Cl2CS were reproduced quite well in the calculations, which can be regarded as a test for the feasibility of the obtained quartic PES.

Dynamical traps lead to the slowing down of intramolecular vibrational energy flow

The phenomenon of intramolecular vibrational energy redistribution (IVR) is at the heart of chemical reaction dynamics. Statistical rate theories, assuming instantaneous IVR, predict exponential decay of the population with the properties of the transition state essentially determining the mechanism. However, there is growing evidence that IVR competes with the reaction timescales, resulting in deviations from the exponential rate law. Dynamics cannot be ignored in such cases for understanding the reaction mechanisms. Significant insights in this context have come from the state space model of IVR, which predicts power law behavior for the rates with the power law exponent, an effective state space dimensionality, being a measure of the nature and extent of the IVR dynamics. However, whether the effective IVR dimensionality can vary with time and whether the mechanism for the variation is of purely quantum or classical origins are issues that remain unresolved. Such multiple power law scalings can lead to surprising mode specificity in the system, even above the threshold for facile IVR. In this work, choosing the well-studied thiophosgene molecule as an example, we establish the anisotropic and anomalous nature of the quantum IVR dynamics and show that multiple power law scalings do manifest in the system. More importantly, we show that the mechanism of the observed multiple power law scaling has classical origins due to a combination of trapping near resonance junctions in the network of classical nonlinear resonances at short to intermediate times and the influence of weak higher-order resonances at relatively longer times.

The Fock space method of vibrational analysis

A reformulation of a semiclassical theory that presently seems uniquely capable of interpreting generic complex multiresonant vibrational spectra is presented. Once given the spectroscopic Hamiltonian which reveals the set of possible resonant couplings and its eigenstates, the new and old formulations both yield without any further computation level by level dynamical assignments for the spectra. Computing a simple trajectory in phase space reveals the motions that when quantized yield the assigned levels. The reformulation introduces two new projected representations of the wave functions. The first is in action space and the second in angle space. The projected representations often allow the reduced angle space, where nodal searches are made, to be of lower dimension than formally occurred. In addition the action representation is a similarly lower dimension lattice representation whose discreteness and regularity allow higher reduced dimensions to be studied. The lattice representation is used to produce a significantly more complete and detailed assignment of the thiophosgene spectrum than previously published.

Assignment and Extracting Dynamics from Experimentally and Theoretically Obtained Spectroscopic Hamiltonians in the Complex Spectral and Classically Chaotic Regions

An analysis of existing algebraic multiresonance spectroscopic Hamiltonians, derived by fitting to experimental data or from classical canonical or quantum Van Vleck perturbation theory, allows without any significant further classical or quantum calculation the assignment of quantum numbers and motions to states observed in spectra that were previously thought to be irregular or just unassignable. In such cases, inspection of the amplitude and phase of eigenfunctions previously calculated in the Hamiltonians derivation process but now transformed to a reduced dimension semiclassical action-angle representation reveals extremely simple albeit unfamiliar topologies that give quantum numbers by simply counting nodes and phase advances. The topology allows these simple wave functions to be sorted into dynamically different excitation ladders or classes of states which are associated with different regions of phase space. The rungs of these ladders or the states in the classes intersperse in energy causing the spectral complexity. No experimental procedure allows such sorting. The success of the work stems from (1) the qualitative insights of nonlinear dynamics, (2) the conversion of the quantum problem in full dimension to a semiclassical one in reduced dimension by use of a canonical transform that takes advantage of the polyad and other constants of the motion, and (3) the judicious choice of the reduced angle variables to reflect rational ratio resonance frequency conditions. This leads to localization of those semiclassical wave functions, which are affected by the particular resonance. In reverse, the localized appearance of the reduced dimension wave function reveals which resonances govern it and makes sorting visually simple. The success of the work also stems from (4) the revealing use of plots of phase advances as well as the usual densities of the eigenstates for sorting and assignment purposes. Even in classically chaotic regions, organizing trajectories, which correspond to averages over regional phase space structures that need not be computed, can easily be drawn as the structure about which eigenfunction localization takes place. The organizing trajectories when transformed back to the full dimensional configuration space reveal the internal molecular motions. The complexity of the usual quantum stationary and propagated wave functions and associated classical trajectories forbids most often such assignments and sorting. This procedure brings the ability to interpret complex vibrational spectra to a degree previously thought possible only for lower excitations. The new methodology replaces and extends the computationally more difficult parts of a procedure used by the authors that was applied successfully to several molecules during the past few years. The new methodology is applied to DCO, CHBrClF, and the bending of acetylene.