Classical path sudden approximations for atom—rigid rotor collisions

The classical path sudden approximation has been used to calculate cross sections for rotational excitation in atom—molecule collisions. Three different decoupling approximations for the angular momentum of the molecule are used: one effective potential approach and two different coupled states approaches, all in their semiclassical versions. Numerical results are presented and discussed for the systems Ar---N₂, Ar---Li₂, Ar---LiH, He---CO and Ar---HCl assumming intermolecular potentials of P₁ symmetry, P₂ symmetry and combinations.     

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.     

Semiclassical bound-continuum Franck-Condon factors uniformly valid at 4 coinciding critical points: 2 Crossings and 2 turning points

A uniform semiclassical approximation to bound-continuum Franck-Condon factors is derived by applying differential topological mapping techniques on a suitable three dimensional integralrepresentation. This approach even holds uniformly if two potentialcurve intersections (real or complex conjugate) and two turning points come close or coincide. The resulting Franck-Condon matrix element is expressed in terms of two generic swallowtail (A₄) integrals whose unfolding parameters are obtained from a single algebraic equation amenable to fast standard routines. Transitional approximations to this result are shown to cover all previously known approaches and to lead to a generalization of a formula of K. Sando and F. H. Mies [1]. A simple and fast trapezoidal method to evaluate the generic swallowtail integrals and other generic integrals of odd determinacy is presented, which even permits the derivation of numerical error bounds.     

Using semiclassical surface hopping for coupled partial wave calculations on systems with non-spherically symmetric potentials

It is shown that a semiclassical surface hopping (SH) approach provides a simple and efficient method for scattering calculations with non-spherically symmetric potentials. The calculations are performed by expanding the wave function in an angular momentum state basis. Since the potential is not spherically symmetric, the different angular states are coupled. The semiclassical SH method, which is typically used for problems with coupled electronic states, can, in principle, be employed for any coupled state problem. The particular SH method employed is known to provide highly accurate results for coupled electronic state problems. The method is tested on model two angular state problems using potential surfaces and couplings arising from a non-spherically symmetric scattering problem. The results for these model problems are in excellent agreement with exact quantum calculations. Full calculations, which are converged with regard to the number of angular basis states, are also performed for the non-spherically symmetric problem. It is shown that an approximation to the surface hopping amplitudes that simplifies the numerical implementation of the method provides results in excellent agreement with the full surface hopping calculation.     

Generator coordinate method in time-dependent density-functional theory: memory made simple

The generator coordinate (GC) method is a variational approach to the quantum many-body problem in which interacting many-body wave functions are constructed as superpositions of (generally nonorthogonal) eigenstates of auxiliary Hamiltonians containing a deformation parameter. This paper presents a time-dependent extension of the GC method as a new approach to improve existing approximations of the exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT). The time-dependent GC method is shown to be a conceptually and computationally simple tool to build memory effects into any existing adiabatic XC potential. As an illustration, the method is applied to driven parametric oscillations of two interacting electrons in a harmonic potential (Hooke's atom). It is demonstrated that a proper choice of time-dependent generator coordinates in conjunction with the adiabatic local-density approximation reproduces the exact linear and nonlinear two-electron dynamics quite accurately, including features associated with double excitations that cannot be captured by TDDFT in the adiabatic approximation.     

Markovian approximation in the relaxation of open quantum systems

In this paper, we examine the validity of the Markovian approximation and the slippage scheme used to incorporate short time transient memory effects in the Markovian master equations (Redfield equations). We argue that for a bath described by a spectral function, J(omega), that is dense and smoothly spread out over the range omega(d), a time scale of tau(b) approximately 1/omega(d) exists; for times of t > tau(b), the Markovian approximation is applicable. In addition, if J(omega) decays to zero reasonably fast in both the omega --> 0 and omega --> infinity limits, then the bath relaxation time, tau(b), is determined by the width of the spectral function and is weakly dependent on the temperature of the bath. On the basis of this criterion of tau(b), a scheme to incorporate transient memory effects in the Markovian master equation is suggested. Instead of using slipped initial conditions, we propose a concatenation scheme that uses the second-order perturbation theory for short time dynamics and the Markovian master equation at long times. Application of this concatenation scheme to the spin-boson model shows that it reproduces the reduced dynamics obtained from the non-Markovian master equation for all parameters studied, while the simple slippage scheme breaks down at high temperatures.     

Importance of anharmonicity, recrossing effects, and quantum mechanical tunneling in transition state theory with semiclassical tunneling. A test case: the H2 + Cl hydrogen abstraction reaction

The hydrogen abstraction reaction from H2 by the Cl atom is studied by means of the variational transition state theory with semiclassical tunneling coefficients on the BW2 potential energy surface. Vibrational anharmonicity and coupling between the bending modes are taken into account. The occurrence of trajectories that recross the transition state is estimated by means of the canonical unified statistical method and by classical trajectories calculations. Different semiclassical methods for tunneling calculations are tested. Our results show that anharmonicity has a small but nonnegligible effect on the thermal rate constants, recrossing can be neglected, and tunneling is adequately described by the least-action approximation, and less successfully by the large-curvature version 3 approximation. However, the large-curvature version 4 and small-curvature approximations lead to a severe underestimation of tunneling. Thermal rate constants calculated using transition state theory including anharmonicity and tunneling agree very well with accurate quantal thermal rate constants over a wide temperature range, although the improvement over the harmonic transition state theory with the microcanonically optimized semiclassical tunneling approximation (based on version 3 of the large-curvature tunneling method) used in a previous study of this reaction is only marginal.     

Non-equilibrium spin-boson model: counting statistics and the heat exchange fluctuation theorem

We focus on the non-equilibrium two-bath spin-boson model, a toy model for examining quantum thermal transport in many-body open systems. Describing the dynamics within the noninteracting-blip approximation equations, applicable, e.g., in the strong system-bath coupling limit and/or at high temperatures, we derive expressions for the cumulant generating function in both the Markovian and non-Markovian limits by energy-resolving the quantum master equation of the subsystem. For a Markovian bath, we readily demonstrate the validity of a steady-state heat exchange fluctuation theorem. In the non-Markovian limit a "weaker" symmetry relation generally holds, a general outcome of microreversibility. We discuss the reduction of this symmetry relation to the universal steady-state fluctuation theorem. Using the cumulant generating function, an analytic expression for the heat current is obtained. Our results establish the validity of the steady-state heat exchange fluctuation theorem in quantum systems with strong system-bath interactions. From the practical point of view, this study provides tools for exploring transport characteristics of the two-bath spin-boson model, a prototype for a nonlinear thermal conductor.     

Theory of kinetic transport coefficients for organic semi-conductors

A generalization of the small-times approximation is considered for an electron-phonon system having a fairly narrow conduction band. The result differs from the semiclassical approximation in giving an analytic expression for the generating function of the semi-invariants. Edgworth's expansion is used as the interpolation function for the electrical-conductivity tensor. More rigorous criteria for the application of the semiclassical approximation are derived for the high-temperature range. Allowance is made for temperature corrections to the generating function, and trial estimates show that this revised function may be applicable down to temperatures well below normal.     

A systematic sequence of relativistic approximations

An approach to the development of a systematic sequence of relativistic approximations is reviewed. The approach depends on the atomically localized nature of relativistic effects, and is based on the normalized elimination of the small component in the matrix modified Dirac equation. Errors in the approximations are assessed relative to four-component Dirac-Hartree-Fock calculations or other reference points. Projection onto the positive energy states of the isolated atoms provides an approximation in which the energy-dependent parts of the matrices can be evaluated in separate atomic calculations and implemented in terms of two sets of contraction coefficients. The errors in this approximation are extremely small, of the order of 0.001 pm in bond lengths and tens of microhartrees in absolute energies. From this approximation it is possible to partition the atoms into relativistic and nonrelativistic groups and to treat the latter with the standard operators of nonrelativistic quantum mechanics. This partitioning is shared with the relativistic effective core potential approximation. For atoms in the second period, errors in the approximation are of the order of a few hundredths of a picometer in bond lengths and less than 1 kJ mol(-1) in dissociation energies; for atoms in the third period, errors are a few tenths of a picometer and a few kilojoule/mole, respectively. A third approximation for scalar relativistic effects replaces the relativistic two-electron integrals with the nonrelativistic integrals evaluated with the atomic Foldy-Wouthuysen coefficients as contraction coefficients. It is similar to the Douglas-Kroll-Hess approximation, and is accurate to about 0.1 pm and a few tenths of a kilojoule/mole. The integrals in all the approximations are no more complicated than the integrals in the full relativistic methods, and their derivatives are correspondingly easy to formulate and evaluate.     

Protein sequence threading,the alignment problem,and a two-step strategy

Conventionally, protein structure prediction via “threading” relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1455–1467, 1999     

Semiclassical glory analyses in the time domain for the H + D2(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction

We make the first application of semiclassical (SC) techniques to the plane-wavepacket formulation of time-domain (T-domain) scattering. The angular scattering of the state-to-state reaction, H + D(2)(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D, is analysed, where v and j are vibrational and rotational quantum numbers, respectively. It is proved that the forward-angle scattering in the T-domain, which arises from a delayed mechanism, is an example of a glory. The SC techniques used in the T-domain are: An integral transitional approximation, a semiclassical transitional approximation, a uniform semiclassical approximation (USA), a primitive semiclassical approximation and a classical semiclassical approximation. Nearside-farside (NF) scattering theory is also employed, both partial wave and SC, since a NF analysis provides valuable insights into oscillatory structures present in the full scattering pattern. In addition, we incorporate techniques into the SC theory called "one linear fit" and "two linear fits", which allow the derivative of the quantum deflection function, Θ?(')(J), to be estimated when Θ?J exhibits undulations as a function of J, the total angular momentum variable. The input to our SC analyses is numerical scattering (S) matrix data, calculated from accurate quantum collisional calculations for the Boothroyd-Keogh-Martin-Peterson potential energy surface No. 2, in the energy domain (E-domain), from which accurate S matrix elements in the T-domain are generated. In the E-domain, we introduce a new technique, called "T-to-E domain SC analysis." It half-Fourier transforms the E-domain accurate quantum scattering amplitude to the T-domain, where we carry out a SC analysis; this is followed by an inverse half-Fourier transform of the T-domain SC scattering amplitude back to the E-domain. We demonstrate that T-to-E USA differential cross sections (DCSs) agree well with exact quantum DCSs at forward angles, for energies where a direct USA analysis in the E-domain fails.     

Reduced density matrix hybrid approach: an efficient and accurate method for adiabatic and non-adiabatic quantum dynamics

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.     

Approaches to large-scale parallel self-consistent field calculations

The availability of massively parallel computers with high computation rates but limited memory and input/output bandwidth prompts the reevaluation of appropriate solution schemes for the self-consistent field (SCF) equations. Several algorithms are considered which exhibit between linear and quadratic convergence using various approximations to the orbital Hessian. A prototype is developed to understand the computational expense of each approach. The optimal choice is found to be a conjugate–gradient method preconditioned with a level-shifted approximation to the orbital Hessian. This is a compromise between efficiency, stability, and low memory usage. Sample benchmarks on two parallel supercomputers are also reported. © 1995 John Wiley & Sons, Inc.     

First principles semiclassical calculations of vibrational eigenfunctions

Vibrational eigenfunctions are calculated on-the-fly using semiclassical methods in conjunction with ab initio density functional theory classical trajectories. Various semiclassical approximations based on the time-dependent representation of the eigenfunctions are tested on an analytical potential describing the chemisorption of CO on Cu(100). Then, first principles semiclassical vibrational eigenfunctions are calculated for the CO(2) molecule and its accuracy evaluated. The multiple coherent states initial value representations semiclassical method recently developed by us has shown with only six ab initio trajectories to evaluate eigenvalues and eigenfunctions at the accuracy level of thousands trajectory semiclassical initial value representation simulations.     

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.