Franck-Condon theory of reactive scattering

A Franck-Condon theory of reactive scattering is introduced which generalizes previous works in a number of respects. We consider the collinear AB + C → A + BC reaction where A, B, and C may be polyatomic species and show how the multi-dimensional continuum-continuum Franck-Condon integral can exactly be reduced to a two-dimensional one involving the nonorthogonal reaction coordinates on the reactant and product diabatic surfaces. These integrals are then written in a rapidly convergent series of products of one-dimensional bound-continuum integrals of a form closely related to those studied in recent theories of photodissociation where accurate analytical and numerical methods are available for treating them. The theory is applied to the case of the D + Hl isotope exchange for a model surface to enable a comparison of exact quantum calculations with those of the Franck-Condon theory for the identical surface in both cases. The calculations are all performed at energies where the reaction is classically forbidden. The relative Dl vibrational distributions (as a function of initial Hl state) are accurately reproduced in a fashion that is fairly insensitive to the choice of the reactant and product diabatic surfaces, but the absolute probabilities are shown to be sensitively dependent on this choice.     

Quantum reactive scattering with a transmission-free absorbing potential

A recently derived transmission-free absorbing potential is applied to the study of atom-diatom chemical reactions. This absorbing potential only depends on a single parameter--the width of the absorbing region--and its reflection properties are guaranteed to improve as this parameter is increased. Converged results can therefore be obtained very easily, as we illustrate with time-dependent wave packet calculations on the H + H2,F + H2, and H + O2 reactions.     

Quantum hydrodynamics: capturing a reactive scattering resonance

The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.     

Quantum reactive scattering of H + hydrocarbon reactions

A practical quantum-dynamical method is described for predicting accurate rate constants for general chemical reactions. The ab initio potential energy surfaces for these reactions can be built from a minimal number of grid points (average of 50 points) and expressed in terms of analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimised using the MP2 method with a cc-pVTZ basis set. Single point energies are calculated on the optimised geometries at the CCSD(T) level of theory with the same basis set. The dynamics of these reactions occur on effective reduced dimensionality hyper-surfaces accounting for the zero-point energy of the optimised degrees of freedom. Bonds being broken and formed are treated with explicit hyperspherical time independent quantum dynamics. Application of the method to the H + CH(4)--> H(2)+ CH(3), H + C(2)H(6)--> H(2)+ C(2)H(5), H + C(3)H(8)--> H(2)+n-C(3)H(7)/H(2)+i-C(3)H(7) and H + CH(3)OH --> H(2)+ CH(3)O/H(2)+ CH(2)OH reactions illustrate the potential of the approach in predicting rate constants, kinetic isotope effects and branching ratios. All studied reactions exhibit large quantum tunneling in the rate constants at lower temperatures. These quantum calculations compare well with the experimental results.     

Quantum hydrodynamics: application to N-dimensional reactive scattering

The quantum hydrodynamic equations associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a new methodology which gives an accurate, unitary, and stable propagation of a time dependent quantum wave packet [B. K. Kendrick, J. Chem. Phys. 119, 5805 (2003)]. The methodology is applied to an N-dimensional model chemical reaction with an activation barrier. A parallel version of the methodology is presented which is designed to run on massively parallel supercomputers. The computational scaling properties of the parallel code are investigated both as a function of the number of processors and the dimension N. A decoupling scheme is introduced which decouples the multidimensional quantum hydrodynamic equations into a set of uncoupled one-dimensional problems. The decoupling scheme dramatically reduces the computation time and is highly parallelizable. Furthermore, the computation time is shown to scale linearly with respect to the dimension N=2,...,100.     

Hyperspherical coordinates in quantum mechanical collinear reactive scattering

A new hyperspherical coordinate method for performing atom—diatom quantum mechanical collinear reactive scattering calculations is described. The method is applicable at energies for which breakup channels are open. Comparison with previous results and new results at high energies for H H₂ are given. The usefulness of this approach is discussed.     

Quantum theory of (femtosecond) time-resolved stimulated Raman scattering

We present a complete perturbation theory of stimulated Raman scattering (SRS), which includes the new experimental technique of femtosecond stimulated Raman scattering (FSRS), where a picosecond Raman pump pulse and a femtosecond probe pulse simultaneously act on a stationary or nonstationary vibrational state. It is shown that eight terms in perturbation theory are required to account for SRS, with observation along the probe pulse direction, and they can be grouped into four nonlinear processes which are labeled as stimulated Raman scattering or inverse Raman scattering (IRS): SRS(I), SRS(II), IRS(I), and IRS(II). Previous FSRS theories have used only the SRS(I) process or only the "resonance Raman scattering" term in SRS(I). Each process can be represented by an overlap between a wave packet in the initial electronic state and a wave packet in the excited Raman electronic state. Calculations were performed with Gaussian Raman pump and probe pulses on displaced harmonic potentials to illustrate various features of FSRS, such as high time and frequency resolution; Raman gain for the Stokes line, Raman loss for the anti-Stokes line, and absence of the Rayleigh line in off-resonance FSRS from a stationary or decaying v=0 state; dispersive line shapes in resonance FSRS; and the possibility of observing vibrational wave packet motion with off-resonance FSRS.     

Quantum mechanical theory of dynamic nuclear polarization in solid dielectrics

Microwave driven dynamic nuclear polarization (DNP) is a process in which the large polarization present in an electron spin reservoir is transferred to nuclei, thereby enhancing NMR signal intensities. In solid dielectrics there are three mechanisms that mediate this transfer--the solid effect (SE), the cross effect (CE), and thermal mixing (TM). Historically these mechanisms have been discussed theoretically using thermodynamic parameters and average spin interactions. However, the SE and the CE can also be modeled quantum mechanically with a system consisting of a small number of spins and the results provide a foundation for the calculations involving TM. In the case of the SE, a single electron-nuclear spin pair is sufficient to explain the polarization mechanism, while the CE requires participation of two electrons and a nuclear spin, and can be used to understand the improved DNP enhancements observed using biradical polarizing agents. Calculations establish the relations among the electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) frequencies and the microwave irradiation frequency that must be satisfied for polarization transfer via the SE or the CE. In particular, if δ, Δ < ω(0I), where δ and Δ are the homogeneous linewidth and inhomogeneous breadth of the EPR spectrum, respectively, we verify that the SE occurs when ω(M) = ω(0S) ± ω(0I), where ω(M), ω(0S) and ω(0I) are, respectively, the microwave, and the EPR and NMR frequencies. Alternatively, when Δ > ω(0I) > δ, the CE dominates the polarization transfer. This two-electron process is optimized when ω(0S(1))-ω(0S(2)) = ω(0I) and ω(M)~ω(0S(1)) or ω(0S(2)), where ω(0S(1)) and ω(0S(2)) are the EPR Larmor frequencies of the two electrons. Using these matching conditions, we calculate the evolution of the density operator from electron Zeeman order to nuclear Zeeman order for both the SE and the CE. The results provide insights into the influence of the microwave irradiation field, the external magnetic field, and the electron-electron and electron-nuclear interactions on DNP enhancements.     

Quantum mechanical embedding theory based on a unique embedding potential

We remove the nonuniqueness of the embedding potential that exists in most previous quantum mechanical embedding schemes by letting the environment and embedded region share a common embedding (interaction) potential. To efficiently solve for the embedding potential, an optimized effective potential method is derived. This embedding potential, which eschews use of approximate kinetic energy density functionals, is then used to describe the environment while a correlated wavefunction (CW) treatment of the embedded region is employed. We first demonstrate the accuracy of this new embedded CW (ECW) method by calculating the van der Waals binding energy curve between a hydrogen molecule and a hydrogen chain. We then examine the prototypical adsorption of CO on a metal surface, here the Cu(111) surface. In addition to obtaining proper site ordering (top site most stable) and binding energies within this theory, the ECW exhibits dramatic changes in the p-character of the CO 4σ and 5σ orbitals upon adsorption that agree very well with x-ray emission spectra, providing further validation of the theory. Finally, we generalize our embedding theory to spin-polarized quantum systems and discuss the connection between our theory and partition density functional theory.     

Quantum wave packet study of N(2D) + H2 reactive scattering

N(²D) + H₂ → NH + H reaction at zero total angular momentum is studied by using a time dependent quantum wave packet method. State‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energy are calculated. The probabilities show many sharp peaks that ascribed to reactive scattering resonances. The probabilities for J > 0 are estimated by using the J‐shifting method. The integral cross sections and thermal rate constants are then calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theories of reactive scattering

This paper is an overview of the theory of reactive scattering, with emphasis on fully quantum mechanical theories that have been developed to describe simple chemical reactions, especially atom-diatom reactions. We also describe related quasiclassical trajectory applications, and in all of this review the emphasis is on methods and applications concerned with state-resolved reaction dynamics. The review first provides an overview of the development of the theory, including a discussion of computational methods based on coupled channel calculations, variational methods, and wave packet methods. Choices of coordinates, including the use of hyperspherical coordinates are discussed, as are basis set and discrete variational representations. The review also summarizes a number of applications that have been performed, especially the two most comprehensively studied systems, H+H2 and F+H2, along with brief discussions of a large number of other systems, including other hydrogen atom transfer reactions, insertion reactions, electronically nonadiabatic reactions, and reactions involving four or more atoms. For each reaction we describe the method used and important new physical insight extracted from the results.     

A quantum mechanical theory for single molecule-single nanoparticle surface enhanced Raman scattering

Single molecule-single nanoparticle surface enhanced Raman scattering event is analyzed using a quantum mechanical approach, resulting in an analytical expression for the electromagnetic enhancement factor that succinctly elucidates the fundamental aspects of SERS. The nanoparticle is treated as a dielectric spherical cavity, and the resulting increase in the spontaneous emission rate of a molecule adsorbed onto the surface of the nanoparticle is examined. The overall enhancement in Raman scattering is due to both the increased local electromagnetic field and the Purcell effect. The predictions of the present model are in agreement with the simulation results of the classical model.     

Decoherence effects in reactive scattering

Decoherence effects on quantum and classical dynamics in reactive scattering are examined using a Caldeira-Leggett type model. Through a study of the dynamics of the collinear H + H2 reaction and the transmission over simple one-dimensional barrier potentials, we show that decoherence leads to improved agreement between quantum and classical reactions and transmission probabilities, primarily by increasing the energy dispersion in a well-defined way. Increased potential nonlinearity is seen to require larger decoherence in order to attain comparable quantum-classical agreement.     

Quantum trajectories for resonant scattering

Time‐dependent wave packet resonant scattering for the double square barrier has been studied in terms of Bohm quantum trajectories. The high transmission probability for the wave packet with a resonant energy can be explained by the behavior of the quantum trajectories under the influence of the relatively slow formation of a node within the first barrier. This node splits the trajectories into reflected and transmitted components. During this stage, many particle trajectories pass through the double‐barrier region and contribute to the transmitted part of the wave packet. Due to the transient nature of the nodes, trajectories in the reflected wave packet bunch together between the nodes for a finite period of time so that temporary structure (localization of particles and accompanying increase in the probability density) develops on small length scales. These calculations also show that the particles gain high momentum near the nodal points, and they reach a uniform momentum distribution after transmitting the barrier region. We have found that the presence of a node between the two barriers influences the different lifetimes of the quasi‐bound states. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 206–213, 2001     

Quantum transition state theory

We derive and discuss upper bound, of transition state form, to the quantum rate constant for bimolecular reactions.     

Quantum mechanical models in catalysis

A wave function of electrons of a catalytic complex taken as a series of products of wave functions of the reactants and catalyst was suggested for use in modeling potential energy surfaces of catalytic reactions and for analysis of catalytic activity. Quantum mechanical criteria at which catalytic transformations become possible were formulated on the basis of this concept. The character of correlations between the activity and physical properties of catalysts was explained, and a general procedure for theoretical analysis of such correlations was described.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya. No. 3, pp. 545–550, March, 1996.     

Quantum mechanical polar surface area

A correlation has been established between the absorbed fraction of training-set molecules after oral administration in humans and the Quantum Mechanical Polar Surface Area (QMPSA). This correlation holds for the QMPSA calculated with structures where carboxyl groups are deprotonated. The correlation of the absorbed fraction and the QMPSA calculated on the neutral gas phase optimized structures is much less pronounced. This suggests that the absorption process is mainly determined by polar interactions of the drug molecules in water solution. Rules are given to derive the optimal polar/apolar ranges of the electrostatic potential.