Infrared chemiluminescence from water-forming reactions: Characterization of dynamics and mechanisms

Vibrational relaxation of HOD by collisions with Ar atoms

Vibrational relaxation of HOD(v12, v3) molecules by collisions with Ar was studied at 298 K (v12 denotes coupled bending, v2, and OD stretching, v1, vibrational modes and v3 denotes OH stretching mode). The vibrationally excited HOD molecules were generated by exothermic abstraction reactions of OD radicals with 13 different RH reactants and observed by infrared emission from a fast-flow reactor as a function of Ar pressure and reaction time. State-specific relaxation rate constants were obtained by comparison of the time evolution of the experimental vibrational distributions with numerical kinetic calculations for vibrational populations. The relaxation mechanism was based on the relaxation scheme of H2O studied earlier with the addition of specific channels for HOD(v12, v3). Unlike H2O, energy in stretching and bending vibrations of HOD cannot be separated due to close ν1 and 2ν2 energies, which leads to fast collisional equilibration between these Fermi-resonant levels. For relaxation of the only pure bending state (10), a rate constant of (1.5 ± 0.3) × 10-13 cm3 molecule-1 s-1 was obtained. The relaxation rate of higher v12 states linearly increases with quantum number and very likely includes transfer of population from OD stretch levels, v1, to a lower energy bend level. The average rate constants for the loss of population from (01), (02), and (03) stretching states are (1.1 ± 0.3) × 10-14, (3.2 ± 1.0) × 10-14, and (5.6 ± 1.2) × 10-14 cm3 molecule-1 s-1, respectively.

First-principles mode-specific reaction dynamics

Controlling the outcome of chemical reactions by exciting specific vibrational and/or rotational modes of the reactants is one of the major goals of modern reaction dynamics studies. In the present Perspective, we focus on first-principles vibrational and rotational mode-specific dynamics computations on reactions of neutral and anionic systems beyond six atoms such as X + C2H6 [X = F, Cl, OH], HX + C2H5 [X = Br, I], OH- + CH3I, and F- + CH3CH2Cl. The dynamics simulations utilize high-level ab initio analytical potential energy surfaces and the quasi-classical trajectory method. Besides initial state specificity and the validity of the Polanyi rules, mode-specific vibrational-state assignment for polyatomic product species using normal-mode analysis and Gaussian binning is also discussed and compared with experiment.

Vibrational Mode-Specific Dynamics of the OH + C2H6 Reaction

We investigate the effects of the initial vibrational excitations on the dynamics of the OH + C2H6 → H2O + C2H5 reaction using the quasi-classical trajectory method and a full-dimensional analytical ab initio potential energy surface. Excitation of the initial CH, CC, and OH stretching modes enhances, slightly inhibits, and does not affect the reactivity, respectively. Translational energy activates the early-barrier title reaction more efficiently than OH and CC stretching excitations, in accord with the Polanyi rules whereas CH stretching modes have similar or higher efficacy than translation, showing that these rules are not always valid in polyatomic processes. Scattering angle, initial attack angle, and product translational energy distributions show the dominance of direct stripping with increasing collision energy, side-on OH and isotropic C2H6 attack preferences, and substantial reactant-product translational energy transfer without any significant mode specificity. The reactant vibrational excitation energy of OH and C2H6 flows into the H2O and C2H5 product vibrations, respectively, whereas product rotations are not affected. The computed mode-specific H2O vibrational distributions show that initial OH excitation appears in the asymmetric stretching vibration of the H2O product and allow comparison with experiments.

High-level ab initio mapping of the multiple H-abstraction pathways of the OH + glycine reaction

We perform a systematic search in the transition-state (TS) and product-channel complex (MIN) regions of the multi-channel OH + glycine → H₂O + H₂N-CH-COOH (CH)/HN-CH₂-COOH (NH)/H₂N-CH₂-COO (COOH) reactions. Geometry optimizations reveal {7, 3, 3} CH-TS, {2, 2, 2} CH-MIN, {17, 10, 5} NH-TS, {35, 19, 19} NH-MIN, and {6, 5, 5} COOH-TS conformers at the {MP2/3-21G, MP2/aug-cc-pVDZ, CCSD(T)-F12b/aug-cc-pVDZ} levels of theory as well as 2 additional CH-TSs based on chemical intuition. The benchmark relative energies of the TS, MIN, and product conformers are obtained by considering basis set effects up to aug-cc-pVQZ using the explicitly-correlated CCSD(T)-F12b method, post-(T) correlation up to CCSDT(Q), core correlation, scalar relativistic effects, spin-orbit coupling, and zero-point energy corrections. All the CH [ΔE_e(ΔH₀) = -38.54(-38.61) kcal mol^-1], NH [ΔE_e(ΔH₀) = -16.72(-17.98) kcal mol^-1], and COOH [ΔE_e = -4.98 kcal mol^-1] reactions are exothermic and proceed via shallow, usually negative, classical(adiabatic) barriers of -0.37(-0.95), -1.91(-2.48), and 1.02(-0.57) kcal mol^-1, respectively. In the entrance channel MRCI/aug-cc-pVTZ computations reveal several complexes with reactive(non-reactive) arrangements and binding energies of 1.0, 1.6, 3.3, (5.2 and 5.9) kcal mol^-1, stabilized by CH⋯OH, NH⋯OH, COOH⋯OH, (OH⋯OC and OH⋯N) hydrogen bonds, respectively.

Analytical potential energy surface and dynamics for the OH + CH3OH reaction

Using as functional form a combination of valence bond and mechanic molecular terms a new full-dimensional potential energy surface was developed for the title reaction, named PES-2022, which was fitted to high-level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples-F12 explicitly correlated level on a representative number of points describing the reactive system. This surface simultaneously describes the two reaction channels, hydrogen abstraction from the methyl group [(R1) path] and from the alcohol group [(R2) path] of methanol to form water. PES-2022 is a smooth and continuous surface, which reasonably describes the topology of this reactive system from reactants to products, including the intermediate complexes present in the system. Based on PES-2022 an exhaustive dynamics study was performed using quasi-classical trajectory calculations under two different initial conditions: at a fixed room temperature, for direct comparison with the experimental evidence and at different collision energies, to analyze possible mechanisms of reaction. In the first case, the available energy was mostly deposited as water vibrational energy, with the vibrational population inverted in the stretching modes and not inverted in the bending modes, reproducing the experimental evidence. In the second case, the analysis of different dynamics magnitudes (excitation functions, product energy partitioning, and product scattering distributions), allows us to suggest different mechanisms for both (R1) and (R2) paths: a direct mechanism for the (R2) path vs an indirect one, related with "nearly trapped" trajectories in the intermediate complexes, for the (R1) path.

Full-dimensional automated potential energy surface development and dynamics for the OH + C2H6 reaction

We develop a full-dimensional analytical potential energy surface (PES) for the OH + C2H6 reaction using the Robosurfer program system, which automatically (1) selects geometries from quasi-classical trajectories, (2) performs ab initio computations using a CCSD(T)-F12/triple-zeta-quality composite method, (3) fits the energies utilizing the permutationally-invariant monomial symmetrization approach, and iteratively improves the PES via steps (1)−(3). Quasi-classical trajectory simulations on the new PES reveal that hydrogen abstraction leading to H2O + C2H5 dominates in the collision energy range of 10−50 kcal/mol. The abstraction cross sections increase and the dominant mechanism shifts from rebound (small impact parameters and backward scattering) to stripping (larger impact parameters and forward scattering) with increasing collision energy as opacity functions and scattering angle distributions indicate. The abstraction reaction clearly favors side-on OH attack over O-side and the least-preferred H-side approach, whereas C2H6 behaves like a spherical object with only slight C−C-perpendicular side-on preference. Collision energy efficiently flows into the relative translation of the products, whereas product internal energy distributions show only little collision energy dependence. H2O/C2H5 vibrational distributions slightly/significantly violate zero-point energy and are nearly independent of collision energy, whereas the rotational distributions clearly blue-shift as collision energy increases.

Current Status of the X + C2H6 [X ≡ H, F(2P), Cl(2P), O(3P), OH] Hydrogen Abstraction Reactions: A Theoretical Review

This paper is a detailed review of the chemistry of medium-size reactive systems using the following hydrogen abstraction reactions with ethane, X + C₂H₆ → HX + C₂H₅; X ≡ H, F(²P), Cl(²P), O(³P) and OH, and focusing attention mainly on the theoretical developments. These bimolecular reactions range from exothermic to endothermic systems and from barrierless to high classical barriers of activation. Thus, the topography of the reactive systems changes from reaction to reaction with the presence or not of stabilized intermediate complexes in the entrance and exit channels. The review begins with some reflections on the inherent problems in the theory/experiment comparison. When one compares kinetics or dynamics theoretical results with experimental measures, one is testing both the potential energy surface describing the nuclei motion and the kinetics or dynamics method used. Discrepancies in the comparison may be due to inaccuracies of the surface, limitations of the kinetics or dynamics methods, and experimental uncertainties that also cannot be ruled out. The paper continues with a detailed review of some bimolecular reactions with ethane, beginning with the reactions with hydrogen atoms. The reactions with halogens present a challenge owing to the presence of stabilized intermediate complexes in the entrance and exit channels and the influence of the spin-orbit states on reactivity. Reactions with O(³P) atoms lead to three surfaces, which is an additional difficulty in the theoretical study. Finally, the reactions with the hydroxyl radical correspond to a reactive system with ten atoms and twenty-four degrees of freedom. Throughout this review, different strategies in the development of analytical potential energy surfaces describing these bimolecular reactions have been critically analyzed, showing their advantages and limitations. These surfaces are fitted to a large number of ab initio calculations, and we found that a huge number of calculations leads to accurate surfaces, but this information does not guarantee that the kinetics and dynamics results match the experimental measurements.

Multichannel dynamics in the OH+ n-butane reaction revealed by crossed-beam slice imaging and quasiclassical trajectory calculations

Multidimensional reactions present various channels that can exhibit very different dynamics and give products of varying subsequent reactivity. Here, we present a combination of experiment and theory to reveal the dynamics of hydrogen abstraction by OH radical at primary and secondary sites in n-butane at a collision energy of 8 kcal/mol. Crossed molecular beam slice imaging experiments unequivocally probe the secondary abstraction channel showing backward angular distributions with mild energy release to product translation, which are accurately captured by trajectory calculations using a specific-reaction-parameter Hamiltonian. Experiments containing both reaction channels indicate a less marked backward character in the angular distribution, whose origin is shown by trajectory calculations to appear as an evolution toward more sideways scattering from the secondary to primary channel. While the two channels have markedly different angular distributions, their energy release is largely comparable, showing ample energy release into the water product. The synergistic combination of crossed-beam imaging and trajectories opens the door to detailed reaction-dynamics studies of chemical reactions with ever-increasing complexity.

Benchmark ab initio characterization of the abstraction and substitution pathways of the OH + CH4/C2H6 reactions

We report benchmark ab initio stationary-point properties for the hydrogen-abstraction, hydrogen-substitution, and methyl-substitution pathways of the OH + CH₄/C₂H₆ reactions.

Kinetics and dynamics study of the OH + C2H6 → H2O + C2H5 reaction based on an analytical global potential energy surface

To describe the gas-phase hydrogen abstraction reaction between the hydroxyl radical and the ethane molecule, an analytical full-dimensional potential energy surface was developed within the Born-Oppenheimer approximation. This reactive process is a ten-body system with 24 degrees of freedom, which represents a theoretical challenge. The new surface, named PES-2020, presents low barrier, 3.76 kcal mol-1, high exothermicity, -16.20 kcal mol-1, and intermediate complexes in the entrance and exit channels. To test the quality and accuracy of the analytical surface several stringent tests were performed and, in general, PES-2020 reasonably simulates the theoretical information used as input, it is a continuous and smooth potential, without spurious minima, it presents great versatility and a reasonable description of this ten-body reaction. Based on this surface, an exhaustive kinetics and dynamics study was performed with a double objective: to analyze the capacity of the new surface to simulate the experimental evidence, and to help understand the mechanism of reaction and the role of the ethyl group in the reaction. In the kinetics study, three approaches were used: variational transition-state theory with multidimensional tunnelling (VTST/MT), ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) results, in the temperature range 200-2000 K. There is general agreement between the three approaches and they reasonably simulate the experimental behaviour, which gives confidence to the fitness of the new surface to describe the system. In the dynamics study, QCT calculations were performed at 298 K for a direct comparison with the only experimental result reported. We found that ethyl fragment presents a noticeable internal energy (∼20%) and so cannot be considered as a spectator. The water product vibrational energy is reasonably reproduced, though when a level-by-level distribution is analyzed the agreement is only qualitative.

Infrared Chemiluminescence Study of the Reaction of Hydroxyl Radical with Formamide and the Secondary Unimolecular Reaction of Chemically Activated Carbamic Acid

Reactions of OH and OD radicals with NH₂CHO and ND₂CHO were studied by Fourier transform infrared emission spectroscopy of the product molecules from a fast-flow reactor at 298 K. Vibrational distributions of the HOD and H₂O molecules from the primary reactions with the C-H bond were obtained by computer simulation of the emission spectra. The vibrational distributions resemble those for other direct H atom abstraction reactions, such as with acetaldehyde. The highest observed level gives an estimate of the C-H bond dissociation energy in formamide of 90.5 ± 1.3 kcal mol^-1. Observation of CO₂, ammonia, and secondary water chemiluminescence gave evidence that recombination of OH and NH₂CO forms carbamic acid (NH₂COOH) with excitation energy of 103 kcal mol^-1, which decomposes through two pathways forming either NH₃ + CO₂ or H₂O + HNCO. The branching fraction for ammonia formation was estimated to be 2-3 times larger than formation of water. This observation was confirmed by RRKM calculation of the decomposition rate constants. A new simulation method was developed to analyze infrared emission from NH₃, NH₂D, ND₂H, and ND₃. Dynamical aspects of the primary and secondary reactions are discussed based on the vibrational distributions of CO₂ and those of H/D isotopes of water and ammonia.

QCT dynamics study of OH/OD + GeH4 reactions. The problem of water bending excitation

The difficulties in the description of the water bending population could be related to the zero-point violation problem in QCT calculations.

Branching Ratios and Vibrational Distributions in Water-Forming Reactions of OH and OD Radicals with Methylamines

Reactions of OH and OD radicals with (CH3)3N, (CH3)2NH, and CH3NH2 were studied by Fourier transform infrared emission spectroscopy (FTIR) of the water product molecules from a fast-flow reactor at 298 K. The rate constants (4.4 ± 0.5) × 10(-11), (5.2 ± 0.8) × 10(-11), and (2.0 ± 0.4) × 10(-11) cm(3) molecule(-1) s(-1) were determined for OD + (CH3)3N, (CH3)2NH, and CH3NH2, respectively, by comparing the HOD emission intensities to the HOD intensity from the OD reaction with H2S. Abstraction from the nitrogen site competes with abstraction from the methyl group, as obtained from an analysis of the HOD and D2O emission intensities from the OD reactions with the deuterated reactants, (CD3)2NH and CD3NH2. After adjustment for the hydrogen-deuterium kinetic isotope effect, the product branching fractions of the hydrogen abstraction from the nitrogen for di- and monomethylamine were found to be 0.34 ± 0.04 and 0.26 ± 0.05, respectively. Vibrational distributions of the H2O, HOD, and D2O molecules are typical for direct hydrogen atom abstraction from polar molecules, even though activation energies are negative because of the formation of pre-transition-state complexes. Comparison is made to the reactions of hydroxyl radicals with ammonia and with other compounds with primary C-H bonds to discuss specific features of disposal of energy to water product.

Rate constant calculations of the GeH4 + OH/OD → GeH3 + H2O/HOD reactions using an ab initio based full-dimensional potential energy surface

2D representation of the analytical potential energy surface. The saddle point and the complexes in the entry and exit channels are included.

Dynamics study of the OH + NH3 hydrogen abstraction reaction using QCT calculations based on an analytical potential energy surface

To understand the reactivity and mechanism of the OH + NH3 → H2O + NH2 gas-phase reaction, which evolves through wells in the entrance and exit channels, a detailed dynamics study was carried out using quasi-classical trajectory calculations. The calculations were performed on an analytical potential energy surface (PES) recently developed by our group, PES-2012 [Monge-Palacios et al. J. Chem. Phys. 138, 084305 (2013)]. Most of the available energy appeared as H2O product vibrational energy (54%), reproducing the only experimental evidence, while only the 21% of this energy appeared as NH2 co-product vibrational energy. Both products appeared with cold and broad rotational distributions. The excitation function (constant collision energy in the range 1.0-14.0 kcal mol(-1)) increases smoothly with energy, contrasting with the only theoretical information (reduced-dimensional quantum scattering calculations based on a simplified PES), which presented a peak at low collision energies, related to quantized states. Analysis of the individual reactive trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) ≥ 6 kcal mol(-1)) all trajectories are direct, at low energies about 20%-30% of trajectories are indirect, i.e., with the mediation of a trapping complex, mainly in the product well. Finally, the effect of the zero-point energy constraint on the dynamics properties was analyzed.

Ab initio based potential energy surface and kinetics study of the OH + NH3 hydrogen abstraction reaction

A full-dimensional analytical potential energy surface (PES) for the OH + NH3 → H2O + NH2 gas-phase reaction was developed based exclusively on high-level ab initio calculations. This reaction presents a very complicated shape with wells along the reaction path. Using a wide spectrum of properties of the reactive system (equilibrium geometries, vibrational frequencies, and relative energies of the stationary points, topology of the reaction path, and points on the reaction swath) as reference, the resulting analytical PES reproduces reasonably well the input ab initio information obtained at the coupled-cluster single double triple (CCSD(T)) = FULL/aug-cc-pVTZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical PES we perform an extensive kinetics study using variational transition-state theory with semiclassical transmission coefficients over a wide temperature range, 200-2000 K. The forward rate constants reproduce the experimental measurements, while the reverse ones are slightly underestimated. However, the detailed analysis of the experimental equilibrium constants (from which the reverse rate constants are obtained) permits us to conclude that the experimental reverse rate constants must be re-evaluated. Another severe test of the new surface is the analysis of the kinetic isotope effects (KIEs), which were not included in the fitting procedure. The KIEs reproduce the values obtained from ab initio calculations in the common temperature range, although unfortunately no experimental information is available for comparison.

Near-infrared chemiluminescence tunable from 900 nm to 1700 nm from narrow-bandgap compounds and polymers

Near-infrared chemiluminescence tunable from 900 nm to 1700 nm was realized from the reaction of the narrow-bandgap compounds and polymers with oxalyl chloride and hydrogen peroxide. A prolonged chemiluminescent process can be achieved in solutions at low temperature or in polymer gels at room temperature.

Characterization of the unimolecular water elimination reaction from 1-propanol, 3,3,3-propan-1-ol-d3, 3,3,3-trifluoropropan-1-ol, and 3-chloropropan-1-ol

The unimolecular reactions of 1-propanol, 3,3,3-propan-1-ol-d3, 3,3,3-trifluoropropan-1-ol, and 3-chloropropan-1-ol have been studied by the chemical activation technique. The recombination of CH3, CD3, CF3, and CH2Cl radicals with CH2CH2OH radicals at room temperature was used to generate vibrationally excited CH3CH2CH2OH, CD3CH2CH2OH, CF3CH2CH2OH, and CH2ClCH2CH2OH molecules. The principal unimolecular reaction for propanol and propanol-d3 with 90 kcal mol(-1) of vibrational energy is 1,2-H2O elimination with rate constants of 3.4 x 10(5) and 1.4 x 10(5) s(-1), respectively. For CH2ClCH2CH2OH also with 90 kcal mol(-1) of energy, 2,3-HCl elimination with a rate constant of 3.0 x 10(7) s(-1) is more important than 1,2-H2O elimination; the branching fractions are 0.95 and 0.05. For CF3CH2CH2OH with an energy of 102 kcal mol(-1), 1,2-H2O elimination has a rate constant of 7.9 x 10(5) and 2,3-HF elimination has a rate constant of 2.6 x 10(5) s(-1). Density functional theory was used to obtain models for the molecules and their transition states. The frequencies and moments of inertia from these models were used to calculate RRKM rate constants, which were used to assign threshold energies by comparing calculated and experimental rate constants. This comparison gives the threshold energy for H2O elimination from 1-propanol as 64 kcal mol(-1). The threshold energies for 1,2-H2O and 2,3-HCl elimination from CH2ClCH2CH2OH were 59 and 54 kcal mol(-1), respectively. The threshold energies for H2O and HF elimination from CF3CH2CH2OH are 62 and 70 kcal mol(-1), respectively. The structures of the transition states for elimination of HF, HCl, and H2O are compared.

Methanesulfinic Acid Reaction with OH: Mechanism, Rate Constants, and Atmospheric Implications

The mechanism for the atmospheric oxidation of methanesulfinic acid (MSIA) has been studied. This is the first theoretical study of the reaction between MSIA and the OH radical. All the possible channels in this reaction have been studied theoretically, and their corresponding rate constants have been evaluated under the variational transition-state theory (VTST) formalism. Two different products can be formed: the CH3S(O)2 radical (which had been experimentally proposed as the only one), and sulfurous acid (H2SO3). The CH3S(O)2 radical can be formed directly or can form via an intermediate adduct, which yields to the radical through the elimination of a water molecule. For the first time, it is theoretically demonstrated that SO2 is formed in the addition channel of the DMS + OH reaction. The consequences of this result in the interpretation of the T-dependence of the SO4(2-)/MSA (methanesulfonic acid) quocient are analyzed. The competition between the unimolecular dissociation of the CH3S(O)2 radical and OH-addition to yield MSA is proposed as one of the possible multiple branching points (along the DMS + OH degradation scheme) influencing the T-dependence of the SO4(2-)/MSA relation.

Laser Probes of Conformational Isomerization in Flexible Molecules and Complexes

Molecules with several flexible coordinates have potential energy surfaces with a large number of minima and many transition states separating them. A general experimental protocol is described that is capable of studying conformational isomerization in such circumstances, measuring the product quantum yields following conformation-specific infrared excitation, and measuring energy thresholds for isomerization of specific X --> Y reactant-product isomer pairs following excitation via stimulated emission pumping (SEP). These methods have been applied to a series of molecules of varying size and conformational complexity, including 3-indolepropionic acid (IPA), meta-ethynylstyrene, N-acetyltryptophan methyl amide (NATMA), N-acetyltryptophan amide (NATA), and melatonin. Studies of isomerization in solute-solvent complexes are also described, including a measurement of the barrier to isomerization in the IPA-H2O complex, and a unique isomerization reaction in which a single water molecule is shuttled between H-bonding sites on the trans-formanilide (TFA) molecule.

The OH* + CH3SH reaction: support for an addition-elimination mechanism from ab initio calculations

Several intermediates for the CH(3)SH + OH(*) --> CH(3)S(*) + H(2)O reaction were identified using MP2(full) 6-311+g(2df,p) ab initio calculations. An adduct, CH(3)S(H)OH(*), I, with electronic energy 13.63 kJ mol(-1) lower than the reactants, and a transition state, II(double dagger), located 5.14 kJ mol(-1) above I, are identified as the entrance channel for an addition-elimination reaction mechanism. After adding zero-point and thermal energies, DeltaH(r,298) ( degrees )(reactants --> I) = -4.85 kJ mol(-1) and DeltaH(298) (double dagger)(I --> II(double dagger)) = +0.10 kJ mol(-1), which indicates that the potential energy surface is broad and flat near the transition state. The calculated imaginary vibrational frequency of the transition state, 62i cm(-1), is also consistent with an addition-elimination mechanism. These calculations are consistent with experimental observations of the OH(*) + CH(3)SH reaction that favored an addition-elimination mechanism rather than direct hydrogen atom abstraction. An alternative reaction, CH(3)SH + OH(*) --> CH(3)SOH + H(*), with DeltaH(r,298) ( degrees ) = +56.94 kJ mol(-1) was also studied, leading to a determination of DeltaH(f,298) ( degrees )(CH(3)SOH) = -149.8 kJ mol(-1).

Near-infrared chemiluminescence from the oxidation of ammonia in aqueous alkaline solution

The chemiluminescent oxidation of ammonia with hypobromite in aqueous alkaline solution evokes a broadly distributed emission in the near-infrared region, with intensity maxima at 1055 nm and 1270 nm.