The Acetylene Bending Spectrum at ∼10000 cm-1: Quantum Assignments in the Midst of Classical Chaos

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.

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.

Topology of the distribution of zeros of the Husimi function in the LiNC/LiCN molecular system

Phase space representations of quantum mechanics constitute useful tools to study vibrations in molecular systems. Among all possibilities, the Husimi function or coherent state representation is very widely used, its maxima indicating which regions of phase space are relevant in the dynamics of the system. The corresponding zeros are also a good indicator to investigate the characteristics of the eigenstates, and it has been shown how the corresponding distributions can discriminate between regular, irregular, and scarred wave functions. In this paper, we discuss how this result can be understood in terms of the overlap between coherent states and system eigenfunctions.

Spectroscopic interpretation: The high vibrations of CDBrClF

We extract the dynamics implicit in an algebraic fitted model Hamiltonian for the deuterium chromophore's vibrational motion in the molecule CDBrClF. The original model has four degrees of freedom, three positions and one representing interbond couplings. A conserved polyad allows in a semiclassical approach the reduction to three degrees of freedom. For most quantum states we can identify the underlying motion that when quantized gives the said state. Most of the classifications, identifications, and assignments are done by visual inspection of the already available wave function semiclassically transformed from the number representation to a representation on the reduced dimension toroidal configuration space corresponding to the classical action and angle variables. The concentration of the wave function density to lower dimensional subsets centered on idealized simple lower dimensional organizing structures and the behavior of the phase along such organizing centers already reveals the atomic motion. Extremely little computational work is needed.