Crossed molecular beam studies on the reaction dynamics of O(1D)+N2O

The reaction of oxygen atom in its first singlet excited state with nitrous oxide was investigated under the crossed molecular beam condition. This reaction has two major product channels, NO+NO and N2+O2. The product translational energy distributions and angular distributions of both channels were determined. Using oxygen-18 isotope labeled O(1D) reactant, the newly formed NO can be distinguished from the remaining NO that was contained in the reactant N2O. Both channels have asymmetric and forward-biased angular distributions, suggesting that there is no long-lived collision complex with lifetime longer than its rotational period. The translational energy release of the N2+O2 channel (fT = 0.57) is much higher than that of the NO+NO channel (fT = 0.31). The product energy partitioning into translational, rotational, and vibrational degrees of freedom is discussed to learn more about the reaction mechanism. The branching ratio between the two product channels was estimated. The 46N2O product of the isotope exchange channel, 18O+44N2O-->16O+46N2O, was below the detection limit and therefore, the upper limit of its yield was estimated to be 0.8%.     

Product state distribution in an idealized collinear reaction

Vibrational transition probabilities for reactive and nonreactive collinear collisions were calculated using an idealized potential surface. Numerically exact results were obtained in an energy region, where seven reactant and seven product channels are energetically accessible. The population of the final vibrational states of products and reactants is discussed.     

The infrared chemiluminescent reaction of Cl atoms with H2CO

The reaction Cl + H₂CO → HCl + HCO has been studied at 295 K. Chlorine atoms were produced via the infrared laser induced dissociation of CCl₃F, using a pulsed CO₂ TEA laser. Using HCl infrared chemiluminescence as the diagnostic, we find the rate constant to be 7.4 ± 0.7 × 10^?11 cm³/molecule sec, in good agreement with several recent studies. An evaluation of TEA laser photolysis as a technique for the generation of chlorine atoms is made, and the relationship of this experiment to recent theories of infrared laser induced chemistry is discussed.     

Quantum reaction boundary to mediate reactions in laser fields

Dynamics of passage over a saddle is investigated for a quantum system under the effect of time-dependent external field (laser pulse). We utilize the recently developed theories of nonlinear dynamics in the saddle region, and extend them to incorporate both time-dependence of the external field and quantum mechanical effects of the system. Anharmonic couplings and laser fields with any functional form of time dependence are explicitly taken into account. As the theory is based on the Weyl expression of quantum mechanics, interpretation is facilitated by the classical phase space picture, while no "classical approximation" is involved. We introduce a quantum reactivity operator to extract the reactive part of the system. In a model system with an optimally controlled laser field for the reaction, it is found that the boundary of the reaction in the phase space, extracted by the reactivity operator, is modulated with time by the effect of the laser field, to "catch" the system excited in the reactant region, and then to "release" it into the product region. This method provides new insights in understanding the origin of optimal control of chemical reactions by laser fields.     

Novel variational principles of chemical reaction

Minimum principles of chemical reaction coordinates are established. IRC (intrinsic reaction coordinate) draws the path of minimum distance from reactant to product. The distance is measured in the rigged configuration Riemannian space whose metric is determined by the distribution of the adiabatic potential energy. Moreover, minimum property of the intrinsic principle of least action is established for the intrinsic dynamism of chemical reaction. Minimum principle of the path connecting intercell boundary with cell is also discussed.     

F114双通道激光解离的动力学研究

本文用TEA CO_2激光器研究了F_(114)(氟利昂-114)激光解离的通道竞争动力学过程。实验发现, F_(114)的激光解离是C—C链和C—Cl键断裂的双通道竞争过程; 其解离过程可以用Arrhenius公式来描述, 并且激光解离的表观活化能与热解离活化能相近, 因而该解离反应具有明显的热解离行为。理论计算还表明, F_(114)的双通道解离过程基本上符合RRKM理论。     

Temperature and pressure dependent rate coefficients for the reaction of Hg with Br and the reaction of Br with Br: a pulsed laser photolysis-pulsed laser induced fluorescence study

A pulsed laser photolysis-pulsed laser induced fluorescence technique has been employed to study the recombination of mercury and bromine atoms, Hg + Br + M --> HgBr + M (1) and the self-reaction of bromine atoms, Br + Br + M --> Br2 + M (2). Rate coefficients were determined as a function of pressure (200-600 Torr) and temperature (243-293 K) in nitrogen buffer gas and as a function of pressure (200-600 Torr) in helium buffer gas at room temperature. For reaction 1, kinetic measurements were performed under conditions in which bromine atoms were the reactant in excess concentration while simultaneously monitoring the concentration of both mercury and bromine. A temperature dependent expression of (1.46 +/- 0.34) x 10(-32) x (T/298)(-(1.86+/-1.49)) cm6 molecule(-2) s(-1) was determined for the third-order recombination rate coefficient in nitrogen buffer gas. The effective second-order rate coefficient for reaction 1 under atmospheric conditions is a factor of 9 smaller than previously determined in a recently published relative rate study. For reaction 2 we obtain a temperature dependent expression of (4.31 +/- 0.21) x 10(-33) x (T/298)(-(2.77+/-0.30)) cm6 molecule(-2) s(-1) for the third-order recombination rate coefficient in nitrogen buffer gas. The rate coefficients are reported with a 2sigma error of precision only; however, due to the uncertainty in the determination of absolute bromine atom concentrations and other unidentified systematic errors we conservatively estimate an uncertainty of +/-50% in the rate coefficients. For both reactions the observed pressure, temperature and buffer gas dependencies are consistent with the expected behavior for three-body recombination.     

Temperature and pressure dependent rate coefficients for the reaction of Hg with Cl and the reaction of Cl with Cl: a pulsed laser photolysis-pulsed laser induced fluorescence study

A pulsed laser photolysis-pulsed laser induced fluorescence technique has been employed to study the recombination of mercury and chlorine atoms, Hg + Cl + M --> HgCl + M (1), and the self-reaction of chlorine atoms, Cl + Cl + M --> Cl(2) + M (2). Rate coefficients were determined as a function of pressure (200-600 Torr) and temperature (243-293 K) in N(2) buffer gas and as a function of pressure (200-600 Torr) in He buffer gas at room temperature. For reaction (1) kinetic measurements were obtained under conditions in which either mercury or chlorine atoms were the reactant in excess concentration while simultaneously monitoring the concentration of both reactants. An Arrhenius expression of (2.2 +/- 0.5) x 10(-32) exp{(680 +/- 400)((1)/(T) - (1)/(298))} cm(6) molecule(-2) s(-1) was determined for the third-order recombination rate coefficient in nitrogen buffer gas. The effective second-order rate coefficient for reaction 1 under atmospheric conditions is much smaller than prior determinations using relative rate techniques. For reaction (2) we obtain an Arrhenius expression of (8.4 +/- 2.3) x 10(-33) exp{(850 +/- 470)((1)/(T) - (1)/(298))} cm(6) molecule(-2) s(-1) for the third-order recombination rate coefficient in nitrogen buffer gas. The rate coefficients are reported with a 2sigma error of precision only; however, due to the uncertainty in the determination of absolute chlorine atom concentrations we conservatively estimate an uncertainty of +/-50% in the rate coefficients. For both reactions the observed pressure, temperature, and buffer gas dependencies are consistent with the expected behavior for three-body recombination.     

Diffusion coefficients as a monitor of reaction kinetics of biological molecules

For revealing spectrally silent dynamics in chemical reactions, a new method, the time-dependent diffusion coefficient, is presented. Principles and typical examples of this method, in particular applications to biologically related reactions, are reviewed. The pulsed laser induced transient grating signal of the photo-decomposition reaction of caged ATP showed that the diffusion coefficient increases gradually with time reflecting the molecular size decrease by the dissociation. Hence, this rate should be a direct measurement of the photo-dissociation rate of ATP from the caged state. In an application to a protein folding reaction, the time-development of the diffusion coefficient was observed during the folding reaction. This time dependence was interpreted in terms of the intermolecular interaction change; i.e., conversion from the intermolecular hydrogen bonding to intramolecular one. It was found that the change of the hydrogen bonding network occurred by the two state manner in entire refolding process of cytochrome c. The unique feature of this time-dependent diffusion coefficient method is discussed.     

QCT and QM calculations of the Cl(2P) + NH3 reaction: influence of the reactant well on the dynamics

A detailed dynamics study, using both quasi-classical trajectory (QCT) and reduced-dimensional quantum mechanical (QM) calculations, was carried out to understand the reactivity and mechanism of the Cl((2)P) + NH(3)→ HCl + NH(2) gas-phase reaction, which evolves through deep wells in the entry and exit channels. The calculations were performed on an analytical potential energy surface recently developed by our group, PES-2010 [M. Monge-Palacios, C. Rangel, J. C. Corchado and J. Espinosa-Garcia, Int. J. Quantum. Chem., 2011], together with a simplified model surface, mod-PES, in which the reactant well is removed to analyze its influence. The main finding was that the QCT and QM methods show a change of the reaction probability with collision energy, suggesting a change of the atomic-level mechanism of reaction with energy. This change disappeared when the mod-PES was used, showing that the behaviour at low energies is a direct consequence of the existence of the reactant well. Analysis of the trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) > 5 kcal mol(-1)) practically all trajectories are direct, at low energies (E(coll) < 3 kcal mol(-1)) the trajectories are indirect, i.e., with the mediation of a trapping complex in the entry and/or the exit wells. The reactant complex allows repeated encounters between the reactants, increasing the reaction probability at low energies. The differential cross section results reinforce this change of mechanism, showing also the influence of the reactant well on this reaction. Thus, the PES-2010 surface yields a forward-backward symmetry in the scattering, while when the reactant well is removed with the mod-PES the shape is more isotropic.     

Cross sections and low temperature rate coefficients for the H + CH+ reaction: a quasiclassical trajectory study

The H + CH(+) reaction is studied by quasiclassical trajectory (QCT) calculations, along with phase space theory (PST) and quantum rigid rotor calculations, employing a global single-valued potential energy surface recently derived by our group. We report QCT total cross sections for each of the three channels, for low collision energies and different reactant rotational quantum numbers. At the lowest collision energies, all cross sections exhibit a capture-like behaviour, as expected from a barrierless reaction. At higher energies, there are important dynamical effects coming from the opening of new channels in the inelastic and reactive exchange collisions. The inelastic cross sections turn out to largely increase, while the reactive abstraction cross sections are declining faster than predicted by the capture theory. A large value of the reactant rotational quantum number tends to suppress these dynamical effects. The QCT rate coefficients are reported for a temperature range from 1-700 K. Below 20 K, the abstraction and exchange QCT rate coefficients are almost constant, as predicted by the capture theory. Above this temperature, the abstraction rate coefficient declines, while the exchange and inelastic rate coefficients are increasing, due to the opening of new channels. A good agreement is observed between the experimental abstraction rate coefficient and the QCT and PST ones. The QCT inelastic results are also compared with those obtained from rigid rotor close coupling (CCRR) calculations in order to check the ability of this approach to provide a reliable estimate of the inelastic rate coefficients for a reactive system without a barrier. The laws of variation as a function of temperature are found to be very similar and the curves are parallel above 20 K. However, reaction is not allowed in the rigid rotor approximation, therefore the CCRR results are about twice as large as their QCT counterparts.     

Steepest descent reaction path integration using a first-order predictor-corrector method

The theoretical treatment of chemical reactions inevitably includes the integration of reaction pathways. After reactant, transition structure, and product stationary points on the potential energy surface are located, steepest descent reaction path following provides a means for verifying reaction mechanisms. Accurately integrated paths are also needed when evaluating reaction rates using variational transition state theory or reaction path Hamiltonian models. In this work an Euler-based predictor-corrector integrator is presented and tested using one analytic model surface and five chemical reactions. The use of Hessian updating, as a means for reducing the overall computational cost of the reaction path calculation, is also discussed.     

Theoretical studies of the reaction of hydroxyl radical with methyl acetate

The mechanisms and the kinetics of the OH (OD) radicals with methyl acetate CH3C(O)OCH3 are investigated theoretically. The dual-level direct dynamics method is employed in the calculation of the rate constants. The optimized geometries and frequencies and the gradients of the stationary points are calculated at the MP2/6-311G(d,p) level. The energetic information of potential energy surfaces is further refined by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MP2/6-311G(d,p) geometries. Four channels are found for the title reaction. The calculated results reveal that there exists an attractive well (reactant complex) in each entrance H-abstraction channel, that is, the H-abstraction reaction makes a stepwise mechanism. The rate constants are calculated by the canonical variational transition-state theory (CVT) with the interpolated single-point energies (ISPE) approach in the temperature range of 200-1200 K. The small-curvature tunneling effect (SCT) approximation is used to evaluate the transmission coefficient. The calculated rate constants are in good agreement with the experimental ones in the measured temperature range. It is shown that the "out-of-plane hydrogen abstraction" from the methoxy end is the dominant channel at the lower temperatures, and the other two H-abstraction channels should be taken into account with the temperatures increasing. The kinetic isotope effects (KIEs) for the three H-abstraction channels and the total reaction are "inverse", and these theoretically calculated KIEs as a function of temperature are expected to be useful for the future laboratory investigation.     

Development and photo-responsive peptide bond cleavage reaction of two-photon near-infrared excitation-responsive peptide

Two-photon near-infrared excitation-responsive amino acid was developed. It was incorporated into a peptide, and focused near-infrared pulsed laser irradiation-induced peptide bond cleavage reaction at the C-terminal position of the photo-responsive amino acid was observed.     

Reactant subspaces and kinetics of chemical reaction networks

This paper studies a chemical reaction network’s (CRN) reactant subspace, i.e. the linear subspace generated by its reactant complexes, to elucidate its role in the system’s kinetic behaviour. We introduce concepts such as reactant rank and reactant deficiency and compare them with their analogues currently used in chemical reaction network theory. We construct a classification of CRNs based on the type of intersection between the reactant subspace R and the stoichiometric subspace S and identify the subnetwork of S-complexes, i.e. complexes which, when viewed as vectors, are contained in S, as a tool to study the network classes, which play a key role in the kinetic behaviour. Our main results on new connections between reactant subspaces and kinetic properties are (1) determination of kinetic characteristics of CRNs with zero reactant deficiency by considering the difference between (network) deficiency and reactant deficiency, (2) resolution of the coincidence problem between the reactant and kinetic subspaces for complex factorizable kinetics via an analogue of the generalized Feinberg–Horn theorem, and (3) construction of an appropriate subspace for the parametrization and uniqueness of positive equilibria for complex factorizable power law kinetics, extending the work of Müller and Regensburger.     

Competitive reaction pathways of C2Cl3 + NO via four-membered ring and bicyclic ring intermediates

The products and mechanisms of the atmospherically and environmentally important reaction, C(2)Cl(3) + NO, are investigated comprehensively by step-scan time-resolved Fourier transform infrared emission spectroscopy and the CCSD(T)/6-311+G(d)//B3LYP/6-311G(d) level of electronic structure calculations. Vibrationally excited products of Cl(2)CO, ClNCO, CCl(3)NCO and NCO have been observed in the IR emission spectra. Cyclic intermediates are found to play important roles leading to the rich variety of the chemical transformations of the reaction. Mainly two competitive reaction pathways are revealed: the four-membered ring intermediate pathway leading to the products Cl(2)CO + ClCN which is essentially barrierless and the bicyclic ring intermediate pathway leading to the product channels of ClNCO + CCl(2,) CCl(3)NCO and CCl(3) + NCO which is rate-limited by a barrier of 42.9 kJ mol(-1) higher than the reactants. By photolyzing the precursor at 248 and 193 nm, respectively, C(2)Cl(3) radicals with different internal energy are produced to observe the product branching ratios as a function of reactant energy. The Cl(2)CO channel via the four-membered ring intermediate pathway is shown to be overwhelmingly dominant at low energy (temperature) but become less important at high energy while the ClNCO and CCl(3)NCO channels via the bicyclic ring intermediate pathway are greatly enhanced and compete effectively. The experimental observation of the products and their branching ratios varying with reactant energy is well consistent with the calculated potential energy profiles.     

An efficient molecular dynamics simulation method for calculating the diffusion-influenced reaction rates

We present a molecular dynamics (MD) simulation method for calculating the diffusion-influenced reaction rates in the limit of low reactant concentrations. To calculate the reaction rate coefficient, we use MD trajectories of a nonreactive equilibrium system that are initiated with a pair of reactant molecules in reactive configuration. Hence reaction systems involving complicated reactant molecules with geometrically restricted reactivities can be treated with comparable efficiency as the simple hard-sphere reaction system. Compared to the similar MD method proposed by Van Beijeren, Dong, and Bocquet [J. Chem. Phys. 114, 6265 (2001)], the present method has a couple of advantages. First, reactions involving more general sink functions can be treated. Second, more accurate results can be obtained when the reaction probability upon collision is less than unity. As an application, we investigate the effects of nondiffusive dynamics and hydrodynamic interaction of reactants on the reaction rate.     

激光等离子体诱导聚甲基硅氧烷的结构重组

将脉冲激光束在惰性气氛中溅射石墨产生的等离子体,与一维链状聚甲基硅氧烷的蒸气束流反应,在产物中通过真空升华和重结晶分离出了两种单晶,经X射线衍射测定为具有(CH3SiO1.5)n(n=8,10)组成的硅氧烷,它们均具有多面体的三维笼状构型。通过色-质联用还检测到其它硅氧烷化合物,其中某些产物已表征为二维多环的构型,反应结果表明,经由激光等离子体的碰撞和能量传递,反应物的链状结构发生解离和进一步反应,经结构重组形成了多种不同组成和构型的产物。     

Vibrational mode and collision energy effects on reaction of H2CO+ with CO2

The effects of collision energy (Ecol) and five different modes of H2CO+ vibration on the title reaction have been studied over the center-of-mass Ecol range from 0.1 to 3.2 eV, including measurements of product ion recoil velocity distributions. Electronic structure and Rice-Ramsperger-Kassel-Marcus calculations were used to examine properties of various complexes and transition states that might be important along the reaction coordinate. Two product channels are observed, corresponding to Hydrogen Transfer (HT) and Proton Transfer (PT). Both channels are endothermic with similar onset energies of approximately 0.9 eV; however, HT dominates over the entire Ecol range and accounts for 70-85% of the total reaction cross section. Both HT and PT occur by direct mechanisms over the entire Ecol range, and have similar dependence on reactant vibrational and collision energy. Despite these similarities, and the fact that the two channels are nearly isoenergetic and differ only in which product moiety carries the charge, their dynamics appear quite different. PT occurs primarily in large impact parameter stripping collisions, where most of the available energy is partitioned to product recoil. HT, in contrast, results in internally hot products with little recoil energy and a more forward-backward symmetric product velocity distribution. Vibration is found to affect the reaction differently in different collision energy regimes. The appearance thresholds are found to depend only on total energy, i.e., all modes of vibration are equivalent to Ecol. With increasing Ecol, vibrational energy becomes increasingly effective, relative to Ecol, at driving reaction. For HT, this transition occurs just above threshold, while for PT it begins at roughly twice the threshold energy.     

Communication: Experimental and theoretical investigations of the effects of the reactant bending excitations in the F + CHD3 reaction

The effects of the reactant bending excitations in the F+CHD(3) reaction are investigated by crossed molecular beam experiments and quasiclassical trajectory (QCT) calculations using a high-quality ab initio potential energy surface. The collision energy (E(c)) dependence of the cross sections of the F+CHD(3)(v(b)=0,1) reactions for the correlated product pairs HF(v('))+CD(3)(v(2)=0,1) and DF(v('))+CHD(2)(v(4)=0,1) is obtained. Both experiment and theory show that the bending excitation activates the reaction at low E(c) and begins to inactivate at higher E(c). The experimental F+CHD(3)(v(b)=1) excitation functions display surprising peak features, especially for the HF(v(')=3)+CD(3)(v(2)=0,1) channels, indicating reactive resonances (quantum effects), which cannot be captured by quasiclassical calculations. The reactant state-specific QCT calculations predict that the v(5)(e) bending mode excitation is the most efficient to drive the reaction and the v(6)(e) and v(5)(e) modes enhance the DF and HF channels, respectively.