Differential cross section for the H+D(2)-->HD(v(')=1,j(')=2,6,10)+D reaction as a function of collision energy

We have measured differential cross sections (DCSs) for the HD (v(')=1,j(')=2,6,10) products of the H+D(2) exchange reaction at five different collision energies in the range 1.48< or =E(coll)< or =1.94 eV. The contribution from the less energetic H atoms formed upon spin-orbit excitation of Br in the photolysis of the HBr precursor is taken into account for two collision energies, E(coll)=1.84 and 1.94 eV, allowing us to disentangle the two different channels. The measured DCSs agree well with new time-dependent quantum-mechanical calculations. As the product rotational excitation increases, the DCSs shift from backward to sideward scattering, as expected. We also find that the shapes of the DCSs show only a small overall dependence on the collision energy, with a notable exception occurring for HD (v(')=1,j(')=2), which appears bimodal at high collision energies. We suggest that this feature results from both direct recoil and indirect scattering from the conical intersection.     

Disagreement between theory and experiment in the simplest chemical reaction: collision energy dependent rotational distributions for H + D2 --> HD(nu' = 3,j') + D

We present experimental rotational distributions for the reaction H + D2 --> HD(nu' = 3,j') + D at eight different collision energies between 1.49 and 1.85 eV. We combine a previous measurement of the state-resolved excitation function for this reaction [Ayers et al., J. Chem. Phys. 119, 4662 (2003)] with the current data to produce a map of the relative reactive cross section as a function of both collision energy and rotational quantum number (an E-j' plot). To compare with the experimental data, we also present E-j' plots resulting from both time-dependent and time-independent quantum mechanical calculations carried out on the BKMP2 surface. The two calculations agree well with each other, but they produce rotational distributions significantly colder than the experiment, with the difference being more pronounced at higher collision energies. Disagreement between theory and experiment might be regarded as surprising considering the simplicity of this system; potential causes of this discrepancy are discussed.     

Crossed-beam study of Co+(3F4)+propane: experiment and density functional theory

A pulsed beam of Co+(3F4) crosses a pulsed beam of C3H8 or C3D8 gas under single collision conditions at collision energies of 0.01 eV and 0.21 eV. After a variable time delay t(ext) = 1-8 micros a fast high voltage pulse extracts product ions into a field-free flight tube for mass analysis. Consistent with earlier work, we observe prompt CoC3H6+ +H2 elimination products in 3:1 excess over CoC2H4+ +CH4 products at 0.21 eV on a 2-10 micros time scale. Long-lived CoC3H8+ complexes fragment predominantly back to Co+ +C3H8 reactants and to H2 elimination products on a 6-24 micros time scale. Density functional theory (B3LYP) calculations provide energies, geometries, and harmonic vibrational frequencies at key stationary points for use in a statistical rate model of the reaction. By adjusting two key multicenter transition state (MCTS) energies downward by 4-7 kcal mol(-1), we obtain good agreement with our decay time results and with the cross section versus collision energy of Armentrout and co-workers from 0.1-1.0 eV. B3LYP theory succeeds in finding relative energies of the MCTSs leading to CH4 and H2 in the proper order to explain the different product branching ratio for Co+ (which favors H2 over CH4) compared with its nearest neighbors Fe+ and Ni+ (which favor CH4 over H2).     

Collision energy dependence of the HD(nu' = 2) product rotational distribution of the H + D2 reaction in the range 1.30-1.89 eV

An experimental and theoretical investigation of the collision energy dependence of the HD(nu' = 2,j') rotational product state distribution for the H + D2 reaction in the collision energy range of Ecol = 1.30-1.89 eV has been carried out. Theoretical results based on time-dependent and time-independent quantum mechanical methods agree nearly perfectly with each other, and the agreement with the experiment is good at low collision energies and very good at high collision energies. This behavior is in marked contrast to a previous report on the HD(nu' = 3,j') product state rotational distribution [Pomerantz et al., J. Chem. Phys. 120, 3244 (2004)] where a systematic difference between experiment and theory was observed, especially at the highest collision energies. The reason for this different behavior is not yet understood. In addition, this study employs Doppler-free spectroscopy to resolve an ambiguity in the E, F-X resonantly enhanced multiphoton ionization transition originating from the HD(nu' = 2,j' = 1) state, which is found to be caused by an accidental blending with the transition coming from the HD(nu' = 1,j' = 14) state.     

Differential cross sections for H + D2 → HD(v' = 2, j' = 0,3,6,9) + D at center-of-mass collision energies of 1.25, 1.61, and 1.97 eV

We have measured differential cross sections (DCSs) for the reaction H + D(2) → HD(v' = 2,j' = 0,3,6,9) + D at center-of-mass collision energies E(coll) of 1.25, 1.61, and 1.97 eV using the photoloc technique. The DCSs show a strong dependence on the product rotational quantum number. For the HD(v' = 2,j' = 0) product, the DCS is bimodal but becomes oscillatory as the collision energy is increased. For the other product states, they are dominated by a single peak, which shifts from back to sideward scattering as j' increases, and they are in general less sensitive to changes in the collision energy. The experimental results are compared to quantum mechanical calculations and show good, but not fully quantitative agreement.     

Quasi-classical trajectory calculations analyzing the dynamics of the C-H stretch mode excitation in the H+CHD3 reaction

A state-to-state dynamics study was performed at a collision energy of 1.53 eV to analyze the effect of the C-H stretch mode excitation on the dynamics of the gas-phase H+CHD3 reaction, which can evolve along two channels, H-abstraction, CD3+H2, and D-abstraction, CHD2+HD. Quasi-classical trajectory calculations were performed on an analytical potential energy surface constructed previously by our group. First, strong coupling between different vibrational modes in the entry channel was observed; i.e., the reaction is non-adiabatic. Second, we found that the C-H stretch mode excitation has little influence on the product rotational distributions for both channels, and on the vibrational distribution for the CD3+H2 channel. However, it has significant influence on the product vibrational distribution for the CHD2+HD channel, where the C-H stretch excitation is maintained in the products, i.e., the reaction shows mode selectivity, reproducing the experimental evidence. Third, the C-H stretch excitation by one quantum increases the reactivity of the vibrational ground-state, in agreement with experiment. Fourth, the state-to-state angular distributions of the CD3 and CHD2 products are reported, finding that for the reactant ground-state the products are practically sideways, whereas the C-H excitation yields a more forward scattering.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

Collision-energy dependence of HD(nu'=1,j') product rotational distributions for the H + D2 reaction

Product rotational distributions for the reaction H + D2 --> HD(nu'=1,j') + D have been measured for 16 collision energies in the range of 1.43 < or = E(coll) < or = 2.55 eV. Time-dependent quantum-mechanical calculations agree well in general with the experimental results, but they consistently yield slightly colder distributions. In terms of the average energy channeled into rotation, the differences between experiment and theory amount to approximately 10% for all collision energies sampled. No peculiarity is found for E(coll)=2.55 eV at which the system has sufficient energy to access the first HD2 electronically excited state.     

Product rotational angular momentum polarization in the H+FCl(v=0-5, j=0, 3, 6, 9)→HF+Cl reaction

The product alignment and orientation of the title reaction on the ground potential energy surface of 1 (2)A' have been studied using the quasi-classical trajectory method. The calculations were carried out for case (a) at collision energies of 0.5-20 kcal mol(-1) with the initially rovibrational state of the reagent FCl molecule being at the v = 0 and j = 0 level to especially reveal in detail the dependence of the product integral cross section on collision energy. Further calculations at the collision energy of 15 kcal mol(-1) for case (b) at v = 0-5, and j = 0, and (c) at v = 0, and j = 3, 6, 9 initial states were carried out to reveal the effect of initially vibrational and rotational excitations on stereodynamics, respectively. Possessing final relative velocity k' (defined as a vector in the xz-plane), product alignment perpendicular to the reagent relative velocity vector k (defined as z- or parallel to the z-axis), for case (a) is found to be weaker at all collision energies, for case (b) is found to be vibrationally enhanced by the reactant molecule FCl, but for case (c), rather insensitive to initially rotational excitation. The rotational vector of product molecular orientation pointing to either negative or positive direction of the y-axis in the center of mass frame, e.g. origin of the coordinate system, is enhanced by collision energies regarding to 0.5-20 kcal mol(-1), while it becomes weaker at higher vibrational (v = 0-5) or rotational (j = 0, 3, 6, 9) excitation levels. Effects of collision energies and of rotational excitation at these collision energies, with 15 kcal mol(-1) as an example on the calculated PDDCSs are also shown and discussed. Detailed plots P(φ(r)) in the range of 0 ≤φ(r)≤ 360(o), and P(θ(r), φ(r)) in the ranges of 0 ≤θ(r)≤ 180° and 0 ≤φ(r)≤ 360° at collision energies 0.5-20 kcal mol(-1) have been presented. Overall, results of PDDCSs of the product alignment and product orientation at these collision energies in the title reaction are not very strongly distinguishable.     

A pulsed-field ionization photoelectron secondary ion coincidence study of the H2+ (X,upsilon+=0-15,N+=1)+He proton transfer reaction

The endothermic proton transfer reaction, H2+(upsilon+)+He-->HeH+ + H(DeltaE=0.806 eV), is investigated over a broad range of reactant vibrational levels using high-resolution vacuum ultraviolet to prepare reactant ions either through excitation of autoionization resonances, or using the pulsed-field ionization-photoelectron-secondary ion coincidence (PFI-PESICO) approach. In the former case, the translational energy dependence of the integral reaction cross sections are measured for upsilon+=0-3 with high signal-to-noise using the guided-ion beam technique. PFI-PESICO cross sections are reported for upsilon+=1-15 and upsilon+=0-12 at center-of-mass collision energies of 0.6 and 3.1 eV, respectively. All ion reactant states selected by the PFI-PESICO scheme are in the N+=1 rotational level. The experimental cross sections are complemented with quasiclassical trajectory (QCT) calculations performed on the ab initio potential energy surface provided by Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. The QCT cross sections are significantly lower than the experimental results near threshold, consistent with important contributions due to resonances observed in quantum scattering studies. At total energies above 2 eV, the QCT calculations are in excellent agreement with the present results. PFI-PESICO time-of-flight (TOF) measurements are also reported for upsilon+=3 and 4 at a collision energy of 0.6 eV. The velocity inverted TOF spectra are consistent with the prevalence of a spectator-stripping mechanism.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

Evidence for excited spin-orbit state reaction dynamics in F+H2: theory and experiment

We describe fully quantum, time-independent scattering calculations of the F+H2-->HF+H reaction, concentrating on the HF product rotational distributions in v'=3. The calculations involved two new sets of ab initio potential energy surfaces, based on large basis set, multireference configuration-interaction calculations, which are further scaled to reproduce the experimental exoergicity of the reaction. In addition, the spin-orbit, Coriolis, and electrostatic couplings between the three quasidiabatic F+H2 electronic states are included. The calculated integral cross sections are compared with the results of molecular beam experiments. At low collision energies, a significant fraction of the reaction is due to Born-Oppenheimer forbidden, but energetically allowed reaction of F in its excited (2P 1/2) spin-orbit state. As the collision energy increases, the Born-Oppenheimer allowed reaction of F in its ground (2P 3/2) spin-orbit state rapidly dominates. Overall, the calculations agree reasonably well with the experiment, although there remains some disagreement with respect to the degree of rotational excitation of the HF(v'=3) products as well as with the energy dependence of the reactive cross sections at the lowest collision energies.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

Cl+HD反应产物的极化与态分布

利用准经典轨线理论,在BW2和G3两个势能面上,研究了Cl+HD反应的动力学.计算结果表明,产物的转动取向对势能面及反应体系的质量因子非常敏感.在BW2势能面上,计算的两个产物的转动取向强于在G3势能面上计算的结果,而无论是在BW2势能面上还是在G3势能面上,DCl产物的取向都强于HCl产物的取向.计算结果还表明,在不同的势能面上反应物的转动激发对反应的影响有着显著的不同.在BW2势能面上,反应物的初始转动激发有利于Cl+HD反应的进行;而在G3势能面上,反应物的初始转动激发消弱了反应的反应性.     

Quantum mechanical calculation of energy dependence of OCl/OH product branching ratio and product quantum state distributions for the O(1D) + HCl reaction on all three contributing electronic state potential energy surfaces

OCl/OH product branching ratios are calculated as a function of total energy for the O( (1) D) + HCl reaction using quantum wavepacket methods. The calculations take account of reaction on all the three electronic state potential energy surfaces which correlate with both reactants and products. Our results show that reaction on the excited electronic state surfaces has a large effect on the branching ratio at higher energies and that these surfaces must therefore be fully taken into account. The calculations use the potential energy surfaces of Nanbu and co-workers. Product vibrational and rotational quantum state distributions are also calculated as a function of energy for both product channels. Inclusion of the excited electronic state potential energy surfaces improves the agreement of the predicted product vibrational quantum state distributions with experiment for the OH product channel. For OCl agreement between theory and experiment is retained for the vibrational quantum state distributions when the excited electronic state potential energy surfaces are included in the analysis. For the rotational state distributions good agreement between theory and experiment is maintained for energies at which experimental results are available. At higher energies, above 0.7 eV of total energy, the OCl rotational state distributions predicted using all three electronic state potential energy surfaces shift to markedly smaller rotational quantum numbers.     

Protonation processes and electronic spectra of histidine and related ions

A full structural assignment of the neutral, protonated, and deprotonated histidine conformers in the gas phase is presented. A total of 3024 unique trial structures were generated by all combinations of internal single-bond rotamers of these species and optimized at the B3LYP/6-311G* level and further optimized at the B3LYP/6-311++G** level. A set of unique conformers is found, and their relative energies, free energies, dipole moments, rotational constants, electron affinities, ionization energies, and harmonic frequencies are determined. The population ratio of histidine and its tautomer is 1:0.16 at 298 K. Massive conformational changes are observed due to protonation and deprotonation, and the intramolecular H-bonds are characterized with the atoms in molecules theory. The calculated proton dissociation energy, gas-phase acidity, proton affinity, and gas-phase basicity are in excellent agreement with the experiments. The deprotonation and protonation of gaseous histidine both occur on the imidazole ring, explaining the versatile biofunctions of histidine in large biomolecules. The UV spectra of neutral and singly and doubly protonated histidine are investigated with the TDDFT/B3LYP/6-311+G(2df,p) calculations. The S0-S1, S0-S2, and S0-S3 excitations of histidine are mixed pipi*/npi* transitions at 5.37, 5.44, and 5.69 eV, respectively. The three excitation energies for histidine tautomer are 4.85, 5.47, and 5.52 eV, respectively. The three excitations for protonated histidine are mainly npi* transitions at 5.45, 5.67, and 5.82 eV, respectively. The S0-S1 excitation of protonated histidine produces ImH-CbetaH2-CalphaH(COOH)-NH2+, while the S0-S2 and S0-S3 transitions produce ImH-CbetaH2-CalphaH(NH2)-(COOH)+. These data may help to understand the mechanisms of the UV fragmentation of biomolecules.     

H+ versus D+) transfer from HOD+ to N2: mode- and bond-selective effects

Reactions of HOD(+) with N(2) have been studied for HOD(+) in its ground state and with one quantum of excitation in each of its vibrational modes: (001)--predominately OH stretch, 0.396 eV, (010)--bend, 0.153 eV, and (100)--predominately OD stretch, 0.293 eV. Integral cross sections and product recoil velocities were recorded for collision energies from threshold to 4 eV. The cross sections for both H(+) and D(+) transfer rise slowly from threshold with increasing collision energy; however, all three vibrational modes enhance reaction much more strongly than equivalent amounts of collision energy and the enhancements remain large even at high collision energy, where the vibration contributes less than 10% of the total energy. Excitation of the OH stretch enhances H(+) transfer by a factor of ～5, but the effect on D(+) transfer is only slightly larger than that from an equivalent increase in collision energy, and smaller than the effect from the much lower energy bend excitation. Similarly, OD stretch excitation strongly enhances D(+) transfer, but has essentially no effect beyond that of the additional energy on H(+) transfer. The effects of the two stretch vibrations are consistent with the expectation that stretching the bond that is broken in the reaction puts momentum in the correct coordinate to drive the system into the exit channel. From this perspective it is quite surprising that bend excitation also results in large (factor of 2) enhancements of both H(+) and D(+) transfer channels, such that its effect on the total cross section at collision energies below ～2 eV is comparable to those from the two stretch modes, even though the bend excitation energy is much smaller. For collision energies above ～2 eV, the vibrational effects become approximately proportional to the vibrational energy, though still much larger than the effects of equivalent addition of collision energy. Measurements of the product recoil velocity distributions show that reaction is direct at all collision energies, with roughly half the products in a sharp peak corresponding to stripping dynamics and half with a broad and approximately isotropic recoil velocity distribution. Despite the large effects of vibrational excitation on reactivity, the effects on recoil dynamics are small, indicating that vibrational excitation does not cause qualitative changes in the reaction mechanism or in the distribution of reactive impact parameters.     

Direct dynamics simulations of O(3P) + HCl at hyperthermal collision energies

The dynamics of the O(3P) + HCl reaction at hyperthermal collision energies were investigated using the quasiclassical trajectory method. Stationary points on the OClH 3A" and 3A' potential energy surfaces (PESs) were also examined. The lowest transition state leading to OCl + H on the 3A" surface is 2.26 eV above the reagents at the CCSD(T)/cc-pVTZ level of theory. This saddle point is bent and product-like. Direct dynamics calculations at the MP2/cc-pVTZ level of theory were used to investigate the excitation functions for OH + Cl, OCl + H, and O + H + Cl formation. OCl is formed mainly from small-impact-parameter collisions, and the OCl + H excitation function peaks around 5 eV, where it is similar in magnitude to the OH + Cl excitation function. The shape of the OCl + H excitation function is discussed, and features are identified that should be general to hyperthermal collision dynamics.