Quantum generalized Langevin equation approach to gas/solid collisions

The classical generalized Langevin equation (GLE) approach to gas/solid collisions is generalized to quantum scattering. Using Feynman's method of partial path integration, the full gas/solid propagator is reduced to a form in which only the dynamics of the incident atom and the surface oscillator(s) directly struck appear explicitly. Solving this effective dynamical problem in the semiclassical limit yields a stationary phase equation of motion identical in form to the classical GLE. The noise, however, is distributed according to quantum rather than classical statistics. From the GLE a quantum phase can be constructed and an S-matrix computed. The resulting theory is capable of describing inelastic-diffractive scattering which has been seen experimentally by Williams.     

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.     

Rainbows and glories in the angular scattering of the state-to-state F + H2 reaction at E(trans)=0.04088 eV

State-of-the-art differential cross sections (DCSs) have been reported by Wang et al. [Proc. Nat. Acad. Sci. (U.S.), 2008, 105, 6227] for the state-to-state F + H(2)→ FH + H reaction using fully quantum-state-selected crossed molecular beams. We theoretically analyze the angular scattering of this reaction, in order to quantitatively understand the physical content of structure in the DCSs. Three transitions are studied, v(i)=0, j(i)=0, m(i)=0 → v(f)=3, j(f)=0, 1, 2, m(f)=0 at a translational energy of 0.04088 eV, where v, j, m are the vibrational, rotational and helicity quantum numbers respectively for the initial and final states. The input to our analyses consists of accurate quantum scattering (S) matrix elements computed for the Fu-Xu-Zhang potential energy surface, as used by Wang et al. in a computational simulation of their experimental DCSs. We prove that the pronounced peak at forward angles observed in the experimental and simulated DCSs for all three transitions is a glory. At larger angles, it is demonstrated that the 000 → 300 and 000 → 310 DCSs both possess a broad farside rainbow, which is accompanied by diffraction oscillations. We confirm the conjecture of Wang et al. that these diffraction oscillations arise from nearside-farside (NF) interference. We find that the reaction is N dominant for all three transitions. The theoretical techniques used to analyze the angular scattering include uniform semiclassical theories of glory and of rainbow scattering. We also make the first application of a semiclassical formula that is uniform for both glory + rainbow scattering. In addition, structure in the DCSs is analyzed using NF theory and local angular momentum theory, in both cases with three resummations of the partial wave series for the scattering amplitude. We make the first explicit application of the Thiele rational interpolation formula to extract the position and residue of the leading Regge pole from a set of S matrix elements, thereby making contact with complex angular momentum theories of DCSs, which interpret the angular scattering in terms of Regge resonances. Our calculations complement the exit-valley vibrationally-adiabatic analysis of Wang et al.     

Scattering of O2 from AL111

Experimental results are presented for the scattering of well-defined beams of molecular oxygen incident on clean Al(111). The data consist of scattered angular distributions measured as a function of incident angle, and for fixed incident angle, the dependence on surface temperature of the angular distributions. The measurements are interpreted in terms of a scattering theory that treats the exchange of energy between the translational and rotational motions of the molecule and the phonons of the surface using classical dynamics. The dependence of the measured angular distributions on incident beam angle and temperature is well explained by the theory. Rotational excitation and quantum excitation of the O(2) internal stretching mode are briefly discussed.     

Scattering of Xe from graphite

Recently a series of experimental measurements for the scattering of Xe atoms from graphite has been reported for both energy-resolved spectra and angular distributions. This system is of fundamental interest because the projectile Xe atoms are considerably more massive than the carbon atoms making up the graphite surface. These measurements were initially analyzed using the hard cubes model and molecular dynamics simulations, and both treatments indicated that the scattering process was a single collision in which the incoming Xe atom interacted strongly with a large number of carbon atoms in the outermost graphite layer. In this work we analyze the data using a single scattering theory that has been shown to explain a number of other experiments on molecular beam scattering from surfaces. These calculations confirm that the scattering process is a single collision with an effective surface mass that is substantially larger than that of the basic graphite ring.     

Semiclassical theory of atom-solid surface collisions: Elastic scattering

The elastic scattering of atoms from solid surfaces is examined within the semiclassical framework. Explicit expression for diffraction intensities are obtained which utilize classical trajectory information as the computational device. Effect of lattice disorder are examined.The results of Beeby and Weinberg are considered.An alternative method for extracting the surface atom displacements and the attractive well depth of the atom-surface interaction is discussed.     

Angular scattering using parameterized S matrix elements for the H + D2(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction: an example of Heisenberg's S matrix programme

A neglected topic in the theory of reactive scattering is the use of parameterized scattering (S) matrix elements to calculate differential cross sections (DCSs). We construct four simple parameterizations, whose moduli are smooth step-functions and whose phases are quadratic functions of the total angular momentum quantum number. Application is made to forward glory scattering in the DCS of the H + D(2)(v(i) = 0, j(i) = 0) → HD(v(f) = 3, j(f) = 0) + D reaction at a translational energy of 1.81 eV, where v and j are vibrational and rotational quantum numbers respectively. The parameterized S matrix elements can reproduce the forward scattering for centre-of-mass reactive scattering angles up to 30° and can identify the total angular momenta (equivalently, impact parameters) that contribute to the glory. The theoretical techniques employed to analyze structure in the DCS include: nearside-farside theory, local angular momentum theory--in both cases incorporating resummations of the partial wave series representation of the scattering amplitude--and the uniform semiclassical theory of forward glory scattering. Our approach is an example of Heisenberg's S matrix programme, in which no potential energy surface is used. Our calculations for the DCS using the four parameterized S matrix elements are counterexamples to the following universal statements often found in the chemical physics literature: "every molecular scattering investigation needs detailed information about the interaction potential," and "an accurate potential energy surface is an essential element in carrying out simulations of a chemical reaction". Both these statements are false.     

On the origin of the forward peak and backward oscillations in the F + H(2)(v = 0) --> HF(v' = 2) + H reaction

We present a semiclassical complex angular momentum (CAM) analysis of the forward scattering peak which occurs at a translational collision energy around 32 meV in the quantum mechanical calculations for the F + H(2)(v = 0, j = 0) --> HF(v' = 2, j' = 0) + H reaction on the Stark-Werner potential energy surface. The semiclassical CAM theory is modified to cover the forward and backward scattering angles. The peak is shown to result from constructive/destructive interference of the two Regge states associated with two resonances, one in the transition state region and the other in the exit channel van der Waals well. In addition, we demonstrate that the oscillations in the energy dependence of the backward differential cross section are caused by the interference between the direct backward scattering and the decay of the two resonance complexes returning to the backward direction after one full rotation.     

Semiclassical theory of predissociation induced by rotational (Coriolis) coupling

The dynamical-state representation and the multichannel quantum-defect theory of predissociation are used to formulate a semiclassical theory of predissociation of a diatomic molecule induced by Coriolis coupling. It is demonstrated that the simple perturbation theory becomes invalid when the rotational angular momentum quantum number is large, and that the non-perturbative formula developed in this paper should be employed.     

Quantum trajectories in elastic atom-surface scattering: threshold and selective adsorption resonances

The elastic resonant scattering of He atoms off the Cu(117) surface is fully described with the formalism of quantum trajectories provided by Bohmian mechanics. Within this theory of quantum motion, the concept of trapping is widely studied and discussed. Classically, atoms undergo impulsive collisions with the surface, and then the trapped motion takes place covering at least two consecutive unit cells. However, from a Bohmian viewpoint, atom trajectories can smoothly adjust to the equipotential energy surface profile in a sort of sliding motion; thus the trapping process could eventually occur within one single unit cell. In particular, both threshold and selective adsorption resonances are explained by means of this quantum trapping considering different space and time scales. Furthermore, a mapping between each region of the (initial) incoming plane wave and the different parts of the diffraction and resonance patterns can be easily established, an important issue only provided by a quantum trajectory formalism.     

Measurements and simulations of high energy O(3P) + Ar(1S) angular scattering: single and multi-collision regimes

We present differential angular cross sections for O(3P) + Ar(1S) scattering at collision energies near 90 kcal mol(-1) (approximately 8 km s(-1) relative velocity) from molecular beam measurements and high-level theoretical calculations. Beams of hyperthermal O(3P) are now being used to investigate novel gas-phase and gas-surface chemistries, and the comparison of theory and measurements on this simple system will be a stringent test of the experimental methodology. Potential energy curves were generated for O(3P) + Ar(1S) using a large cc-pVQZ basis within a valence multi-configuration plus perturbation theory treatment. These curves were then used in quantum scattering calculations to generate differential cross sections. Agreement between experiment and theory is excellent. In addition to these comparisons, the cross sections were used in direct simulation Monte Carlo calculations to investigate effects of increasing the Ar flux above the "single-collision" regime. As the Ar flux increases, the observed differential angular cross sections change in two ways. In addition to the main "single-scatter" peak along the incident O-atom beam direction, a secondary O-atom peak appears in the direction of the incident Ar beam, and the multiple-scattered O-atom translational energy starts to reflect the energy of the relatively slow moving Ar beam.     

Collision dynamics of non-integrable systems: Validity of classical scaling

The classical collision dynamics of a model atom—molecule non-integrable collision system is studied, and the energy transfer (ET) moment is examined as a function of the initial semiclassical level of the molecule. A recently derived classical scaling theory is shown to be valid in the case when the molecular motion remains regular throughout the collision, and the ET variation is then characterized by a polynomial dependence on the initial (semiclassical) quantum numbers. When chaotic motions participate, the ET no longer follows the scaling law. The utility of the scaling theory in providing the proper interpolation form for extending classical trajectory data in non-integrable collision systems is discussed.     

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.     

Harmonic analysis and discrete polynomials. From semiclassical angular momentum theory to the hyperquantization algorithm

Orthogonal polynomials of a discrete variable have been widely investigated as fundamental tools of numerical analysis. This work aims to propose the extension of their use to quantum mechanical problems. By exploiting both their connection with coupling and recoupling coefficients of angular momentum theory and their asymptotic relationships (semiclassical limit) with spherical and hyperspherical harmonics, a discretization procedure, the hyperquantization algorithm, has been developed and applied to the study of anisotropic interactions and of reactive scattering. One of the most appealing features of this method turns out to be a drastic reduction of memory requirements and computing time for extensive dynamical calculations. Examples of the application of this technique to stereodirected dynamics via an exact representation for the S matrix as well as to the characterization of molecular beam polarization are also illustrated. Received: 17 September 1999 / Accepted: 3 February 2000 / Published online: 5 June 2000     

Semiclassical self-consistent field perturbation theory for the hydrogen atom in a magnetic field

A recently developed perturbation theory for solving self-consistent field equations is applied to the hydrogen atom in a strong magnetic field. This system has been extensively studied using other methods and is therefore a good test case for the new method. The perturbation theory yields summable large-order expansions. The accuracy of the self-consistent field approximation varies according to field strength and quantum state but is often higher than the accuracy from adiabatic approximations. A new derivation is presented for the asymptotic adiabatic approximation, the most useful of the adiabatic approaches. This derivation uses semiclassical perturbation theory without invoking an adiabatic hypothesis. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 183–192, 1998     

Collisional line shapes for low frequency vibrations of adsorbates on a metal surface

The dynamics of atoms or molecules adsorbed on a metal surface, and excited by collisions with an atomic beam, are treated within a theory that includes energy dissipation into lattice vibrations by means of a frequency and temperature dependent friction function. The theory provides dynamic structure factors for energy transfer derived from collisional time correlation functions. It describes the relaxation of a vibrationally excited atom or molecule within a model of a damped quantum harmonic oscillator bilinearly coupled to a bath of lattice oscillators. The collisional time correlation function is generalized to include friction effects and is applied to the vibrational relaxation of the frustrated translation mode of Na adsorbed on a Cu(001) surface, CO on Cu(001), and CO on Pt(111), following excitation by collisions with He atoms. Results for the frequency shift and width of line shapes versus surface temperature are in very good agreement with experimental measurements of inelastic He atom scattering. Our interpretation of the experimental results provides insight on the relative role of phonon versus electron-hole relaxation.     

Semiclassical quantum unimolecular reaction rate theory revisited

We describe a semiclassical quantum unimolecular reaction rate theory derived from the corresponding classical theory developed by Davis, Gray, Rice and Zhao (DGRZ). The analysis retains the intuitively useful mechanistic distinctions between intramolecular energy transfer and reaction, with the consequence that the semiclassical quantum theory version neglects some interference effects in the reaction dynamics. In the limiting case that intramolecular energy transfer is very fast compared to the rate of reaction we show that the DGRZ representation of the rate constant can be transformed, using the Weyl correspondence between quantum operators and classical variables, to the quantum flux–flux correlation function representation of the rate constant. In the more general case that the rate of intramolecular energy transfer influences the reaction dynamics, the semiclassical representation of the Wigner function for a classical system with both quasiperiodic and chaotic motion is used to obtain the reaction rate constant. Our analysis identifies the quantum analogue of the classical bottleneck to intramolecular energy transfer with the scars of unstable periodic orbits; it leads to a flux–flux correlation function representation of the rate constant for intramolecular energy transfer.     

Steric effects and quantum interference in the inelastic scattering of NO(X) + Ar

Rotationally inelastic collisions of NO(X) with Ar are investigated in unprecedented detail using state-to-state, crossed molecular beam experiments. The NO(X) molecules are selected in the Ω = 0.5, j = 0.5, f state and then oriented such that either the ‘N’ or ‘O’ end of the molecule is directed towards the incoming Ar atom. Velocity map ion imaging is then used to probe the scattered NO molecules in well-defined quantum states. We show that the fully quantum state-resolved differential steric asymmetry, which quantifies how the relative efficiency for scattering off the ‘O’ and the ‘N’ ends of the molecule varies with scattering angle, is strongly affected by quantum interference. Significant changes in both integral and differential cross sections are found depending on whether collisions occur with the N or O ends of the molecule. The results are well accounted for by rigorous quantum mechanical calculations, in contrast to both classical trajectory calculations and more simplistic models that provide, at best, an incomplete picture of the dynamics.     

Modification of the static polarizability by a perfect metallic surface

We calculate the modifications of the static polarizability of a general system (atom or molecule) induced by a perfect metallic surface. The calculations use the Feynman propagator and take retardation effects into account. For any distance from the mirror the induced polarizability is composed of two terms: a classical one without retardation and a quantum one which is always retarded. Two limiting cases are then considered: the Casimir limit where the classical term dominates and the near-zone limit (London limit) where the quantum correction must be added to the classical term.