Semiclassical calculation of the vibrational echo

The infrared echo measurement probes the time scales of the molecular motions that couple to a vibrational transition. Computation of the echo observable within rigorous quantum mechanics is problematic for systems with many degrees of freedom, motivating the development of semiclassical approximations to the nonlinear optical response. We present a semiclassical approximation to the echo observable, based on the Herman-Kluk propagator. This calculation requires averaging over a quantity generated by two pairs of classical trajectories and associated stability matrices, connected by a pair of phase-space jumps. Quantum, classical, and semiclassical echo calculations are compared for a thermal ensemble of noninteracting anharmonic oscillators. The semiclassical approach uses input from classical mechanics to reproduce the significant features of a complete, quantum mechanical calculation of the nonlinear response.     

Geometric constraints in semiclassical initial value representation calculations in Cartesian coordinates: accurate reduction in zero-point energy

An approach for the inclusion of geometric constraints in semiclassical initial value representation calculations is introduced. An important aspect of the approach is that Cartesian coordinates are used throughout. We devised an algorithm for the constrained sampling of initial conditions through the use of multivariate Gaussian distribution based on a projected Hessian. We also propose an approach for the constrained evaluation of the so-called Herman-Kluk prefactor in its exact log-derivative form. Sample calculations are performed for free and constrained rare-gas trimers. The results show that the proposed approach provides an accurate evaluation of the reduction in zero-point energy. Exact basis set calculations are used to assess the accuracy of the semiclassical results. Since Cartesian coordinates are used, the approach is general and applicable to a variety of molecular and atomic systems.     

Dynamical hydrogen atom tunneling in dichlorotropolone: a combined quantum, semiclassical, and classical study

Based on the Cartesian Reaction Surface framework we construct a four-dimensional potential for the tropolone derivative 3,7-dichlorotropolone, a molecule with an intramolecular O-H...O hydrogen bond. The reduced configuration space involves the in-plane hydrogen atom coordinates, a symmetric O-O vibrational mode, and an antisymmetric mode related to deformations of the seven-membered ring. The system is characterized in terms of quantum mechanical computations of the low-lying eigenstates as well as a classical and semiclassical analysis of spectra obtained via Fourier transforming autocorrelation functions. For the semiclassical analysis we utilize the amplitude-free correlation function method [K. Hotta and K. Takatsuka, J. Phys. A 36, 4785 (2003)]. Our results demonstrate substantial anharmonic couplings leading to highly correlated wave functions even at moderate energies. Furthermore, the importance of dynamical tunneling in tropolone is suggested since many low-lying states--including the ground state--lie above the classical saddle point but nevertheless appear as split pairs.     

Exact quantum statistics for electronically nonadiabatic systems using continuous path variables

We derive an exact, continuous-variable path integral (PI) representation of the canonical partition function for electronically nonadiabatic systems. Utilizing the Stock-Thoss (ST) mapping for an N-level system, matrix elements of the Boltzmann operator are expressed in Cartesian coordinates for both the nuclear and electronic degrees of freedom. The PI discretization presented here properly constrains the electronic Cartesian coordinates to the physical subspace of the mapping. We numerically demonstrate that the resulting PI-ST representation is exact for the calculation of equilibrium properties of systems with coupled electronic and nuclear degrees of freedom. We further show that the PI-ST formulation provides a natural means to initialize semiclassical trajectories for the calculation of real-time thermal correlation functions, which is numerically demonstrated in applications to a series of nonadiabatic model systems.     

A semiclassical density functional evaluation of the smoothly varying part of the Hartree–Fock binding energy in atoms

A semiclassical density functional approach is used to evaluate the smooth part of the variation of the Hartree–Fock (HF ) binding energy in atoms, from helium through xenon. The energy density functional is chosen with an improved form for the kinetic energy functional (KEF ). The variation of the HF binding energy is split into smooth and oscillating (shell-structure) parts, in accordance with Strutinskyøs self-consistent averaging procedure, which is equivalent to a semiclassical ?-expansion of the KEF . This enables a well-grounded evaluation of the oscillating part of the HF binding energy, which displays a quasi-periodic behavior expressing clearly the shell structure of atoms. © 1994 John Wiley & Sons, Inc.     

A comparison of methods to compute the potential of mean force.

Most processes occurring in a system are determined by the relative free energy between two or more states because the free energy is a measure of the probability of finding the system in a given state. When the two states of interest are connected by a pathway, usually called reaction coordinate, along which the free-energy profile is determined, this profile or potential of mean force (PMF) will also yield the relative free energy of the two states. Twelve different methods to compute a PMF are reviewed and compared, with regard to their precision, for a system consisting of a pair of methane molecules in aqueous solution. We analyze all combinations of the type of sampling (unbiased, umbrella-biased or constraint-biased), how to compute free energies (from density of states or force averaging) and the type of coordinate system (internal or Cartesian) used for the PMF degree of freedom. The method of choice is constraint-bias simulation combined with force averaging for either an internal or a Cartesian PMF degree of freedom.     

Independent trajectory implementation of the semiclassical Liouville method: application to multidimensional reaction dynamics

We describe an independent trajectory implementation of semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, a single ensemble of trajectories describes all semiclassical density matrix elements of a coupled electronic state problem, with the ensemble evolving classically under a single reference Hamiltonian chosen on the basis of physical grounds. In this paper, we introduce an additional uncoupled trajectory approximation, allowing the members of the ensemble to evolve independently of one another and eliminating the major computational costs of our previous coupled trajectory implementation. The accuracy of the method is demonstrated for model one-dimensional problems. In addition, the approach is applied to the chemical reaction dynamics of a collinear triatomic system, yielding excellent agreement with exact calculations. This method allows molecular dynamics involving coupled electronic surfaces to be modeled with essentially the same effort as classical molecular dynamics and ensemble averaging.     

Cartesian coupled coherent states simulations: Ne(n)Br2 dissociation as a test case

In this article, we describe coupled coherent states (CCS) simulations of vibrational predissociation of weakly bounded complexes. The CCS method is implemented in the Cartesian frame in a manner that is similar to classical molecular dynamics. The calculated lifetimes of the vibrationally excited Ne-Br(2)(ν) complexes agree with experiment and previous calculations. Although the CCS method is, in principle, a fully quantum approach, in practice it typically becomes a semiclassical technique at long times. This is especially true following dissociation events. Consequently, it is very difficult to converge the quantum calculations of the final Br(2) vibrational distributions after predissociation and of the autocorrelation functions. However, the main advantage of the method is that it can be applied with relative ease to determine the lifetimes of larger complexes and, in order to demonstrate this, preliminary results for tetra- and penta-atomic clusters are reported.     

On the efficient path integral evaluation of thermal rate constants within the quantum instanton approximation

We present an efficient path integral approach for evaluating thermal rate constants within the quantum instanton (QI) approximation that was recently introduced to overcome the quantitative deficiencies of the earlier semiclassical instanton approach [Miller, Zhao, Ceotto, and Yang, J. Chem. Phys. 119, 1329 (2003)]. Since the QI rate constant is determined solely by properties of the (quantum) Boltzmann operator (specifically, by the zero time properties of the flux-flux and delta-delta correlation functions), it can be evaluated by well-established techniques of imaginary time path integrals even for quite complex chemical reactions. Here we present a series of statistical estimators for relevant quantities which can be evaluated straightforwardly with any nonlinear reaction coordinates and general Hamiltonians in Cartesian space. To facilitate the search for the optimal dividing surfaces required by the QI approximation, we introduce a two-dimensional quantum free energy surface associated with the delta-delta correlation function and describe how an adaptive umbrella sampling can be used effectively to construct such a free energy surface. The overall computational procedure is illustrated by the application to a hydrogen exchange reaction in gas phase, which shows excellent agreement of the QI rates with those obtained from quantum scattering calculations.     

Adiabatic states derived from a spin-coupled diabatic transformation: semiclassical trajectory study of photodissociation of HBr and the construction of potential curves for LiBr+

The development of spin-coupled diabatic representations for theoretical semiclassical treatments of photodissociation dynamics is an important practical goal, and some of the assumptions required to carry this out may be validated by applications to simple systems. With this objective, we report here a study of the photodissociation dynamics of the prototypical HBr system using semiclassical trajectory methods. The valence (spin-free) potential energy curves and the permanent and transition dipole moments were computed using high-level ab initio methods and were transformed to a spin-coupled diabatic representation. The spin-orbit coupling used in the transformation was taken as that of atomic bromine at all internuclear distances. Adiabatic potential energy curves, nonadiabatic couplings and transition dipole moments were then obtained from the diabatic ones and were used in all the dynamics calculations. Nonadiabatic photodissociation probabilities were computed using three semiclassical trajectory methods, namely, coherent switching with decay of mixing (CSDM), fewest switches with time uncertainty (FSTU), and its recently developed variant with stochastic decoherence (FTSU/SD), each combined with semiclassical sampling of the initial vibrational state. The calculated branching fraction to the higher fine-structure level of the bromine atom is in good agreement with experiment and with more complete theoretical treatments. The present study, by comparing our new calculations to wave packet calculations with distance-dependent ab initio spin-orbit coupling, validates the semiclassical trajectory methods, the semiclassical initial state sample scheme, and the use of a distance-independent spin-orbit coupling for future applications to polyatomic photodissociation. Finally, using LiBr(+) as a model system, it is shown that accurate spin-coupled potential curves can also be constructed for odd-electron systems using the same strategy as for HBr.     

Abstractive dissociation of oxygen over Al(111): a nonadiabatic quantum model

The dissociation of oxygen on a clean aluminum surface is studied theoretically. A nonadiabatic quantum dynamical model is used, based on four electronically distinct potential energy surfaces characterized by the extent of charge transfer from the metal to the adsorbate. A flat surface approximation is used to reduce the computation complexity. The conservation of the helicopter angular momentum allows Boltzmann averaging of the outcome of the propagation of a three degrees of freedom wave function. The dissociation event is simulated by solving the time-dependent Schr?dinger equation for a period of 30 femtoseconds. As a function of incident kinetic energy, the dissociation yield follows the experimental trend. An attempt at simulation employing only the lowest adiabatic surface failed, qualitatively disagreeing with both experiment and nonadiabatic calculations. The final products, adsorptive dissociation and abstractive dissociation, are obtained by carrying out a semiclassical molecular dynamics simulation with surface hopping which describes the back charge transfer from an oxygen atom negative ion to the surface. The final adsorbed oxygen pair distribution compares well with experiment. By running the dynamical events backward in time, a correlation is established between the products and the initial conditions which lead to their production. Qualitative agreement is thus obtained with recent experiments that show suppression of abstraction by rotational excitation.     

Phase quantization of chaos in the semiclassical regime

Since the early stage of the study of Hamilton chaos, semiclassical quantization based on the low-order Wentzel-Kramers-Brillouin theory, the primitive semiclassical approximation to the Feynman path integrals (or the so-called Van Vleck propagator), and their variants have been suffering from difficulties such as divergence in the correlation function, nonconvergence in the trace formula, and so on. These difficulties have been hampering the progress of quantum chaos, and it is widely recognized that the essential drawback of these semiclassical theories commonly originates from the erroneous feature of the amplitude factors in their applications to classically chaotic systems. This forms a clear contrast to the success of the Einstein-Brillouin-Keller quantization condition for regular (integrable) systems. We show here that energy quantization of chaos in semiclassical regime is, in principle, possible in terms of constructive and destructive interference of phases alone, and the role of the semiclassical amplitude factor is indeed negligibly small, as long as it is not highly oscillatory. To do so, we first sketch the mechanism of semiclassical quantization of energy spectrum with the Fourier analysis of phase interference in a time correlation function, from which the amplitude factor is practically factored out due to its slowly varying nature. In this argument there is no distinction between integrability and nonintegrability of classical dynamics. Then we present numerical evidence that chaos can be indeed quantized by means of amplitude-free quasicorrelation functions and Heller's frozen Gaussian method. This is called phase quantization. Finally, we revisit the work of Yamashita and Takatsuka [Prog. Theor. Phys. Suppl. 161, 56 (2007)] who have shown explicitly that the semiclassical spectrum is quite insensitive to smooth modification (rescaling) of the amplitude factor. At the same time, we note that the phase quantization naturally breaks down when the oscillatory nature of the amplitude factor is comparable to that of the phases. Such a case generally appears when the Planck constant of a large magnitude pushes the dynamics out of the semiclassical regime.     

On nonadiabatic transition state theory

A nonseparable semiclassical transition state approximation for reactions involving more than one electronic surface is suggested. The single surface formulation in terms of quasiprobability distributions used by Miller is discussed along with a separable semiclassical approximation for the nonadiabatic rate suggested in the Soviet literature. A thermally averaged nonadiabatic rate is defined, and a semiclassical approximation is presented, wherein the surface through which flux is calculated in the transition state approach is determined by the intersection of adiabatic electronic surfaces viewed as functions of imaginary (or complex) time.     

The optical rotatory dispersion and circular dichroism of molecules containing interacting residues

We use a semiclassical formalism to derive expressions for the optical rotatory dispersion and circular dichroism associated with a dilute system of molecules represented as freely-rotating, rigid arrays of interacting residues. The relationship of the present approach to DeVoe's classical theory is discussed. Far from resonance, our results become equivalent to the equations of Moffitt, Fitts, and Kirkwood when the separation between adjacent residues is small compared to the wavelength of the light. In the resonance case, a calculation of the frequency-dependent circular dichroism per residue is carried out for a particular geometrical arrangement of the residues. These curves, for arrays both finite and infinite in extent, are compared with those obtained when a certain commonly used approximation to the orientational averaging is introduced. It is shown that this approximation is not always valid.     

Semiclassical initial value series solution of the spin boson problem

A numerical solution for the quantum dynamics of the spin boson problem is obtained using the semiclassical initial value series representation approach to the quantum dynamics. The zeroth order term of the series is computed using the new forward-backward representation for correlation functions presented in the preceding adjacent paper. This leads to a rapid convergence of the Monte Carlo sampling, as compared to previous attempts. The zeroth order results are already quite accurate. The first order term of the series is small, demonstrating the rapid convergence of the semiclassical initial value representation series. This is the first time that the first order term in the semiclassical initial value representation series has been converged for systems with the order of 50 degrees of freedom.     

Nonlinear optical properties of confined excitons in clusters

Size effects in femtosecond photon echo spectroscopy of neat clusters are calculated using a quasiparticle representation of the nonlinear response. We extend our previous study of cooperative effects on the nonlinear response of assemblies of two level molecules [J. A. Leegwater and S. Mukamel, Phys. Rev. A46, 452 (1992)] to allow for nuclear motion and to have an s-p model of polarizable atoms. Photon echos in Benzene/Argon clusters are calculated using a semiclassical phase averaging procedure [L. E. Fried and S. Mukamel, Adv. Chem. Phys. (in Press)].     

Motional narrowing in the time-averaging approximation for simulating two-dimensional nonlinear infrared spectra

The diagonal linewidth in two-dimensional infrared spectra is often narrower than the distribution of transition frequencies. The width along the antidiagonal is broader than predicted by the lifetime broadening. These effects arise from time-dependent fluctuations of the transition frequencies. They can be accounted for with a semiclassical approach. For systems with many coupled vibrational modes, this approach, however, becomes computationally too demanding to be practically applicable. A time-averaging approximation was suggested for linear infrared absorption spectra. In this paper, we demonstrate that the averaging can be optimized to fit a broader scale of frequency fluctuations by using a Gaussian weight function instead of the originally proposed box function. We further generalize the time-averaging method to allow the simulation of two-dimensional infrared spectra and demonstrate the method on a simple system. The approximation delivers a large speed-up of the calculation without losing significant accuracy.     

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.     

Semiclassical extension of the Landau-Teller theory of collisional energy transfer

A semiclassical version of the quantum coupled-states approximation for the vibrational relaxation of diatomic molecules in collisions with monatomic bath gases is presented. It is based on the effective mass approximation and a recovery of the semiclassical Landau exponent from the classical Landau-Teller collision time. For an interaction with small anisotropy, the Landau exponent includes first order corrections with respect to the orientational dependence of the collision time and the effective mass. The relaxation N(2)(v=1)-->N(2)(v=0) in He is discussed as an example. Employing the available vibrationally elastic potential, the semiclassical approach describes the temperature dependence of the rate constant k(10)(T) over seven orders of magnitude across the temperature range of 70-3000 K in agreement with experimental data and quantum coupled-states calculations. For this system, the hierarchy of corrections to the Landau-Teller conventional treatment in the order of importance is the following: quantum effects in the energy release, dynamical contributions of the rotation of N(2) to the vibrational transition, and deviations of the interaction potential from a purely repulsive form. The described treatment provides significant simplifications over complete coupled-states calculations such that applications to more complex situations appear promising.