Parson, R.P. & Heller, E.J. Classical and semiclassical aspects of dissipative molecular processes: Radiative and nonradiative vibrational relaxation in polyatomic molecules. J. Chem. Phys. 85, 2581 (1986). Publisher's VersionAbstract

The time dependence of spontaneous emission from isolated, highly vibrationally excited molecules is studied from the point of view of the classical and semiclassical mechanics of dissipatively perturbed Hamiltonian systems. A correlation function expression for the time‐dependent emission spectrum serves as a starting point for quasiclassical approximations. These in turn lead to an investigation of the classical dynamics of model molecular systems in which small, energy‐nonconserving terms have been added to Hamilton’s equations. Numerical calculations show rich dynamical behavior which can be qualitatively understood in terms of the resonance structure of the unperturbed system. For example, trajectories tend to be captured by zones of nonlinear resonance. This ‘‘mode locking’’ produces a characteristic cleanup of an emission spectrum that otherwise is rather congested at these energies. The close connection between spectra and dynamics suggests that the classical dynamics of dissipatively perturbed systems may provide a useful language for interpreting radiative and condensed‐phase vibrational relaxation, and possibly some types of intramolecular relaxation phenomena as well.

Heller, E.J. Qualitative properties of eigenfunctions of classically chaotic Hamiltonian systems. Quantum Chaos and Statistical Nuclear Physics: Proceedings of the 2nd International Conference on Quantum Chaos and the 4th International Colloquium on Statistical Nuclear Physics, Held at Cuernavaca, México, January 6–10, 1986 162–181 (1986). Publisher's VersionAbstract

We begin by discussing the properties expected of eigenfunctions of a classically chaotic Hamiltonian system, using simple Correspondence Principle arguments. The properties involve nodal surfaces, coordinate and momentum space amplitude distribution, and phase space distribution. The eigenfunctions of the stadium billiard are examined, and it is found that the periodic orbits of shortest periods and smallest stability parameter profoundly affect the eigenfunctions: “scars” of higher wavefunction density surround the special periodic orbits. Finally a theory is presented for the scars, showing that they must exist, and relating them directly to the special periodic orbits. These same periodic orbits cause level density fluctuations.

Heller, E.J. & Reimers, J.R. Recent Progress in Semiclassical Wavepacket Methods for Chemical Physics Applications. Dynamics of Wave Packets in Molecular and Nuclear Physics: Proceedings of the International Meeting Held in Priorij Corsendonck, Belgium July 2–4, 1985 79-81 (1986). Publisher's Version
Parson, R.P. & Heller, E.J. Semiclassical mechanics of dissipative molecular processes: The correspondence principle and population relaxation. J. Chem. Phys. 85, 2569 (1986). Publisher's VersionAbstract

The relaxation of a primary system coupled weakly to a bath of environmental modes is examined from the standpoint of recent developments in the semiclassical theory of molecular bound states. Emphasis is placed upon highly excited, strongly nonlinear (but quasiperiodic) primary systems and zero temperature baths. The starting point for the analysis is a master equation for the populations of the eigenstates of the primary system. The correspondence principle provides semiclassical approximations to the transition rates, allowing quantum state populations to be calculated from classical trajectories. A second semiclassical approximation leads to an equation of motion for a probability density in the classical action variables. As h→0, this density agrees with the density generated by running an ensemble of damped classical trajectories and averaging out the angle variables; retention of terms of order h provides smoothed quantum corrections. Numerical examples of both semiclassical approximations are presented.

Reimers, J.R., Wilson, K.R., Heller, E.J. & Langhoff, S.R. CASSCF-wave packet ab initio prediction of electronic and vibrational spectra: Application to the A(²Pi) reverse arrow X(². sigma. /sup +/) absorption of C2H at 3000 K. J. Chem. Phys. 82, 5064 (1985). Publisher's VersionAbstract

The combination of a b i n i t i o calculation of electronic wave functions with a wave packet calculation of the nuclear motion is used, within the Born–Oppenheimer approximation to compute the vibrational and electronic absorption of a polyatomic molecule. A particular virtue of this approach is that high as well as low temperature spectra are both calculable. This method is applied to C2H, for which the complete active space self‐consistent field (CASSCF) method is used to determine full Born–Oppenheimer potential surfaces. Using the assumption that the A(2Π) ← X(2Σ+) absorption can be written as the sum of the A(2A’) ← X and A(2A‘) ← Xabsorptions, the spectra are determined to 60 cm1 resolution at a temperature of 3000 K. As a result of the large thermal bending amplitude at 3000 K, the calculated spectra are broad and have little resolved structure. Two bands are resolvable, one is due to the A(2A‘) ← Xabsorption and is centered at 5500 cm1, while the other is due to A(2A’) ← Xabsorption and is centered at 9500 cm1. The dramatic blue shift of the A(2A’) ← X band results from the combination of the large X state thermal bending amplitude and high bending frequency of the A(2A’) state. We also determine the X state pure vibrational absorptionspectrum and show it to be of much lower intensity than the pure electronic spectrum.

Reimers, J.R. & Heller, E.J. The exact eigenfunctions and eigenvalues of a two-dimensional rigid rotor obtained using Gaussian wave packet dynamics. J. Chem. Phys. 83, 511 (1985). Publisher's VersionAbstract

Exact eigenfunctions for a two‐dimensional rigid rotor are obtained using Gaussian wave packet dynamics. The wave functions are obtained by propagating, without approximation, an infinite set of Gaussian wave packets that collectively have the correct periodicity, being coherent states appropriate to this rotational problem. This result leads to a numerical method for the semiclassical calculation of rovibrational, molecular eignestates. Also, a simple, almost classical, approximation to full wave packet dynamics is shown to give exact results: this leads to an a p o s t e r i o r i justification of the De Leon–Heller spectral quantization method.

Reimers, J.R. & Heller, E.J. The exact thermal rotational spectrum of a two-dimensional rigid rotor obtained using Gaussian wave packet dynamics. J. Chem. Phys. 83, 516, (1985). Publisher's VersionAbstract


The exact thermal rotational spectrum of a two‐dimensional rigid rotor is obtained using Gaussian wave packet dynamics. The spectrum is obtained by propagating, without approximation, infinite sets of Gaussian wave packets. These sets are constructed so that collectively they have the correct periodicity, and indeed, are coherent states appropriate to this problem. Also, simple, almost classical, approximations to full wave packet dynamics are shown to give results which are either exact or very nearly exact. Advantages of the use of Gaussian wave packet dynamics over conventional linear response theory are discussed.


Heller, E.J. & Sundberg, R.L. Quantum Ergodicity and Intensity Fluctuations. Chaotic Behavior in Quantum Systems 255-292 (1985). Publisher's VersionAbstract

Flow in phase space is intimately related to the concept of ergodicity. An ergodic system will sample all regions of phase space democratically in the time average, subject only to known a priori constants of the motion. For convenience, we shall refer to the accessible phase space as the “energy shell.” This shell may have some thickness if a distribution of energies are initially populated An energy envelope is defined by the initial distribution of energy and the envelope plays a role both in classical and quantum phase Space flow, as discussed before and as will become evident in the theory and examples to follow.

Reimers, J.R. & Heller, E.J. The exact thermal rotational spectrum of a two-dimensional rigid rotor obtained using Gaussian wave packet dynamics. Journal of Chemical Physics 83, 516 (1985). Publisher's VersionAbstract
The exact thermal rotational spectrum of a two-dimensional rigid rotor is obtained using Gaussian wave packet dynamics. The spectrum is obtained by propagating, without approximation, infinite sets of Gaussian wave packets. These sets are constructed so that collectively they have the correct periodicity, and indeed, are coherent states appropriate to this problem. Also, simple, almost classical, approximations to full wave packet dynamics are shown to give results which are either exact or very nearly exact. Advantages of the use of Gaussian wave packet dynamics over conventional linear response theory are discussed. [ABSTRACT FROM AUTHOR]
Davis, M.J. & Heller, E.J. Comparisons of classical and quantum dynamics for initially localized states. J. Chem. Phys. 80, 5036 (1984). Publisher's VersionAbstract

We compare the dynamics of quantum wave packets with the dynamics of classical trajectory ensembles. The wave packets are Gaussian with expectation values of position and momenta which centers them in phase space. The classical trajectory ensembles are generated directly from the quantum wave packets via the Wigner transform. Quantum and classical dynamics are then compared using several quantum measures and the analogous classical ones derived from the Wigner equivalent formalism. Comparisons are made for several model potentials and it is found that there is generally excellent classical–quantum correspondence except for certain specific cases of tunneling and interference. In general, this correspondence is also very good in regions of phase space where there is classical chaos.

Eaker, C.W., Schatz, G.C., Leon, N.D. & Heller, E.J. Fourier transform methods for calculating action variables and semiclassical eigenvalues for coupled oscillator systems. J. Chem. Phys. 81, 5913 (1984). Publisher's VersionAbstract

Two methods for calculating the good action variables and semiclassical eigenvalues for coupled oscillator systems are presented, both of which relate the actions to the coefficients appearing in the Fourier representation of the normal coordinates and momenta. The two methods differ in that one is based on the exact expression for the actions together with the EBK semiclassical quantization condition while the other is derived from the Sorbie–Handy (SH) approximation to the actions. However, they are also very similar in that the actions in both methods are related to the same set of Fourier coefficients and both require determining the perturbed frequencies in calculating actions. These frequencies are also determined from the Fourier representations, which means that the actions in both methods are determined from information entirely contained in the Fourier expansion of the coordinates and momenta. We show how these expansions can very conveniently be obtained from fast Fourier transform (FFT) methods and that numerical filtering methods can be used to remove spurious Fourier components associated with the finite trajectory integration duration. In the case of the SH based method, we find that the use of filtering enables us to relax the usual periodicity requirement on the calculated trajectory. Application to two standard Henon–Heiles models is considered and both are shown to give semiclassical eigenvalues in good agreement with previous calculations for nondegenerate and 1:1 resonant systems. In comparing the two methods, we find that although the exact method is quite general in its ability to be used for systems exhibiting complex resonant behavior, it converges more slowly with increasing trajectory integration duration and is more sensitive to the algorithm for choosing perturbed frequencies than the SH based method. The SH based method is less straightforward to use in studying resonant systems, but good results are obtained for 1:1 resonant systems using actions defined in terms of the complex coordinates Q 1±i Q 2. The SH based method is also shown to be remarkably accurate in determining high energy eigenvalues (about three‐quarters of the dissociation energy).

Sundberg, R.L. & Heller, E.J. Predissociation rates as a probe of intramolecular dynamics. J. Chem. Phys. 80, 3680 (1984). Publisher's VersionAbstract

We investigate the consequences of quasiperiodic and ergodic classical dynamics on predissociation rates from bound vibrational states into a continuum. We find a strong correlation between the quasiperiodic/ergodic motion and rate constants for certain locations of the predissociationsurface. Other locations for the predissociative surface are insensitive to the difference between classical motion which is confined to small regions of phase space and motion which samples a large portion of its available phase space.

Leon, N.D., Davis, M.J. & Heller, E.J. Quantum manifestations of classical resonance zones. J. Chem. Phys. 80, 794 (1984). Publisher's VersionAbstract

We examine the concept of nodal breakup of wave functions as a criterion for quantum mechanical ergodicity. We find that complex nodal structure of wave functions is n o t sufficient to determine quantum mechanical ergodicity. The influence of classical resonances [which manifest themselves as classical resonance zones (CRZ)] may also be responsible for the seeming complexity of nodal structure. We quantify this by reexamining one of the two systems studied by Stratt, Handy, and Miller [J. Chem. Phys. 7 1, 3311 (1974)] from both a quantum mechanical and classical point of view. We conclude that quasiperiodic classical motion can account for highly distorted quantum eigenstates. One should always keep this in mind when addressing questions regarding quantum mechanical ergodicity.

Leon, N.D. & Heller, E.J. Semiclassical spectral quantization: Application to two and four coupled molecular degrees of freedom. J. Chem. Phys. 81, 5957 (1984). Publisher's VersionAbstract

Semiclassical quantization of the quasiperiodic vibrational motion of molecules is usually based on Einstein–Brillouin–Keller (EBK) conditions for the quantization of the classical actions. Explicit use of the EBK conditions for molecular systems of K degrees of freedom requires K quantization conditions. Therefore, explicit use of the EBK conditions becomes increasingly difficult if not impossible for polyatomic systems of three or more degrees of freedom. In this paper we propose a semiclassical quantization method which makes explicit use of phase coherence of the de Broglie wave associated with the trajectory rather than the EBK conditions. We show that taking advantage of phase coherence reduces the K quantization conditions to a single quantum condition—regardless of the number of degrees of freedom. For reasons that will become obvious we call this method ‘‘spectral quantization.’’ Polyatomic vibrational wave functions and energy eigenvalues are generated from quasiperiodic classical trajectories. The spectral method is applied to an ABA linear triatomic molecule with two degrees of freedom and to an anharmonic model of the molecule cyanoacetylene. The usefulness of the technique is demonstrated in this latter calculation since the cyanoacetylene model will have four coupled vibrational degrees of freedom.

Reimers, J.R., Wilson, K.R. & Heller, E.J. Complex time dependent wave packet technique for thermal equilibrium systems: Electronic spectra. J. Chem. Phys. 79, 4749 (1983). Publisher's VersionAbstract

A time dependent wave packet method is presented for the rapid calculation of the properties of systems in thermal equilibrium and is applied, as an illustration, to electronic spectra. The thawed Gaussian approximation to quantum wave packet dynamics combined with evaluation of the density matrix operator by imaginary time propagation is shown to give exact electronic spectra for harmonic potentials and excellent results for both a Morse potential and for the band contours of the three transitions of the visible electronic absorptionspectrum of the iodine molecule. The method, in principle, can be extended to many atoms (e.g., condensed phases) and to other properties (e.g., infrared and Raman spectra and thermodynamic variables).

Bergsma, J.P., Berens, P.H., Wilson, K.R., Fredkin, D.R. & Heller, E.J. Electronic spectra from molecular dynamics: A simple approach. J. Phys. Chem. 88, 3, 612-619 (1983). Publisher's Version
Heller, E.J. & Brown, R.C. Radiationless transitions in a new light. J. Chem. Phys. 79, 3336 (1983). Publisher's VersionAbstract

Radiationless transitions in polyatomic molecules prove to be quite amendable to a semiclassical treatment both below and above crossings between the potential surfaces involved in the transition. Below such crossings, tunneling integrals are easily performed which give good estimates of the dependence of the nonradiative rate on the energy gap and excess energy in the electronic state. Above the surface crossing, the transitions become classically allowed and a Tully–Preston surface hopping model suffices. We find that a nonlinear dependence of ln(k n r ) vs E plots is the rule rather than the exception. The ln(k n r ) vs E plots tend to flatten out with increasing energy. This effect can occur below surface crossings, but is most dramatic when a surface crossing is reached. The recent beam results of Smalley and co‐workers on pyrazine and pyrimidine are seen to be a possible case of this simple behavior.

Leon, N.D. & Heller, E.J. Semiclassical quantization and extraction of eigenfunctions using arbitrary trajectories. J. Chem. Phys. 78, 4005 (1983). Publisher's VersionAbstract

We present in this paper a coordinate independent semiclassical quantization method. We demonstrate that in order to extract accurate eigenvalues and eigenfunctions the trajectory does not necessarily have to reside on the quantizing torus, rather, one can use information obtained on arbitrary tori. Because the method is coordinate independent, no difficulty is encountered in quantizing within classical resonance zones. Furthermore, nearby eigenstates and eigenvalues (nearby in action space) may be extracted from the same trajectory—this is especially convenient when the density of states becomes large.

Myers, A.B., Mathies, R.A., Tannor, D.J. & Heller, E.J. Excited state geometry changes from preresonance Raman intensities: Isoprene and hexatriene. J. Chem. Phys. 77, 3857 (1982). Publisher's VersionAbstract

A method is presented for using a single preresonance Raman spectrum and an absorptionspectrum to obtain changes in equilibrium geometry upon electronic excitation. The relative displacements along each of the vibrational normal coordinates are obtained from the Raman intensities, while the overall scaling of the displacements is determined by the absorption band shape. The absorption spectra, as well as Raman excitation profiles, are calculated using either a sum over vibronic states or a formally equivalent time‐dependent method [S.‐Y. Lee and E. J. Heller, J. Chem. Phys. 71, 4777 (1979)]. The time‐dependent method is computationally much faster than the vibronic sum for large multidimensional systems. Our analysis, which assumes isolated molecules and separable, harmonic surfaces, yields a good fit to the vapor phase absorptionspectrum of t r a n s‐hexatriene with a Lorentzian linewidth of 175 cm−1. However, the diffuse absorptionspectrum of isoprene cannot be adequately reproduced using Lorentzian line shapes, even when all 33 normal modes are included. Finite temperature and excited state frequency changes are also found to have little effect on the calculated band shapes. These results suggest that inhomogeneous broadening may be a major factor, but calculations using Gaussian broadening fail to accurately reproduce the experimental spectrum.

Bado, P., et al. Picosecond dynamics of I2 photodissociation. Picosecond Phenomena III: Proceedings of the Third International Conference on Picosecond Phenomena Garmisch-Partenkirchen, Fed. Rep. of Germany June 16--18, 1982 23, 260-263 (1982). Publisher's VersionAbstract

While liquid solution reactions are much more important in chemistry, gas phase reactions are much better understood. Given the central importance of solution reactions to inorganic, organic, industrial and biochemistry, it is rather surprising that, as yet, there is not a single such reaction whose molecular dynamics are understood in detail. Theoretical and experimental evidence already makes clear that much of the important molecular dynamic action in solution reactions occurs on the picosecond and subpicosecond time scales. The dihalogen photodissociation and recombination reactions, X 2 + hvX + XX 2, involving the simplest possible molecular reactants and products, diatomics, and in rare gas solution involving only two elements, seem excellent candidates for study.