Search

The Hydroxymethyl (•CH2OH) radical is an important intermediate species in both atmosphere and combustion reaction systems. The rate coefficients for •CH2OH + 3O2 and (CH2OH + 3O2 (+H2O) reactions were calculated using Rice−Ramsperger−Kassel−Marcus (RRKM)/master equation (ME) simulation and canonical variational transition state theory (CVT) between the temperature range of 200 to 1500 K based on the potential energy surface constructed using CCSD(T)//ωB97XD /6-311++G(3df,3pd). The results show that •CH2OH + 3O2 leads to form CH2O and HO2 in the temperature 1000K). When water molecule added to the reaction favored the formation of CH2O and HO2 in all the temperature. The calculated rate coefficient for the •CH2OH + 3O2 (2.8×10-11 cm3 molecule-1 s-1 at 298K) is in good agreement with the previous experimental values (~1×10-11 cm3 molecule-1 s-1 at 298K). The rate coefficients for the water-assisted reaction (2.4×10-16 cm3 molecule-1 s-1 at 1000K) is at least 3-4 order magnitude smaller than the water-free reaction (6.2×10-12 cm3 molecule-1 s-1 at 1000 K). This result is consistent with the similar types of reaction system. Our calculations also predict that the effect of a single water molecule favors the formation of CH2O in the combustion condition however, the water-free reaction favors the formation of CH2O in the atmospheric condition. The current study helps to understand how a single water molecule changes the reaction mechanism and chemical kinetic behaviour under atmospheric and combustion conditions.

An in-depth theoretical analysis of key chemical equilibria in Signal Amplification by Reversible Exchange (SABRE) is provided, employing density functional theory calculations to characterize the likely reaction network. For all reactions in the network, the potential energy surface is probed to identify minimum energy pathways. Energy barriers and transition states are calculated, and harmonic transition state theory is applied to calculate exchange rates that approximate experimental values. The reaction network energy surface can be modulated by chemical potentials that account for the dependence on concentration, temperature, and partial pressure of molecular constituents (hydrogen, methanol, pyridine) supplied to the experiment under equilibrium conditions. We show that, under typical experimental conditions, the Gibbs free energies of the two key states involved in pyridine-hydrogen exchange at the common Ir-IMes catalyst system in methanol are essentially the same, i. e., nearly optimal for SABRE. We also show that a methanol-containing intermediate is plausible as a transient species in the process.

8-Oxo-2'-deoxyguanosine (OG) is the most common DNA lesion. Notably, OG becomes more susceptible to oxidative damage than the undamaged nucleoside, forming mutagenic products in vivo. Herein the reactions of singlet O2 with the radical cations of 8-oxo-2'-deoxyguanosine (OG.+ ) and 9-methyl-8-oxoguanine (9MOG.+ ) were investigated using ion-molecule scattering mass spectrometry, from which barrierless, exothermic O2 -addition products were detected for both reaction systems. Corroborated by static reaction potential energy surface constructed using multi-reference CASPT2 theory and molecular dynamics simulated in the presence of the reactants' kinetic and internal energies, the C5-terminal O2 -addition was pinpointed as the most probable reaction pathway. By elucidating the reaction mechanism, kinetics and dynamics, and reaction products and energetics, this work constitutes the first report unraveling the synergetic damage of OG by ionizing radiation and singlet O2 .

A new global potential energy surface (PES) for the ground state of the SH2+(X4A′′) system is constructed using a permutation invariant polynomial neural network method. In ab initio calculations, the MRCI-F12 method with the AVTZ basis set is used. Furthermore, the dynamics calculations of the S+ + H2(v = 0–3, j = 0) → SH+ + H reaction are carried out based on the new PES. The reaction probabilities and integral cross sections are compared with available theoretical calculations. Present values are in general good agreement with the previous theoretical studies. However, some discrepancies can still be found due to different PESs used in the calculation. Furthermore, the vibrational energy of the reactant molecule can significantly enhance the reactivity compared to the translational energy. The differential cross sections indicated that the reaction mechanism is changed from the “head-on” rebound mechanism to the tripping mechanism with the increasing number of initial vibrational excitation state.

In contrast to their spontaneous deprotonation in aqueous solution, reactions of guanine and guanosine radical cations with water in the gas phase are exclusively initiated by hydration of the radical cations as reported in recent work (Y. Sun et al., Phys. Chem. Chem. Phys., 2018, 20, 27510). As gas-phase hydration reactions closely mimic the actual scenario for guanine radical cations in double-stranded DNA, exploration of subsequent reactions within their water complexes can provide an insight into the resulting oxidative damage to nucleosides. Herein guided-ion beam mass spectrometry experiment and direct dynamics trajectory simulations were carried out to examine prototype complexes of the 9-methylguanine radical cation with one and two water ligands (i.e., 9MG˙+·(H2O)1–2) in the gas phase, wherein the complexes were activated by collisional activation in the experiment and by thermal excitation at high temperatures in the simulations. Guided by mass spectroscopic measurements, trajectory results and reaction potential energy surface, three reaction pathways were identified. The first two reaction pathways start with H-atom abstraction from water by the O6 and N7 atoms in 9MG˙+ and are referred to as HAO6 and HAN7, respectively. The primary products of HAO6 and HAN7 reactions, including [9MG + HO6]+/[9MG + HN7]+ and ˙OH, react further to either form [8OH-9MG + HO6]˙+ and [8OH-9MG + HN7]˙+via C8-hydroxylation or form radical cations of 6-enol-guanine (6-enol-G˙+) and 7H-guanine (7HG˙+) via SN2-type methanol elimination. The third reaction pathway corresponds to the formation of 8OH-9MG+ by H elimination from the complex, referred to as HE. Among these product channels, [8OH-9MG + HN7]˙+ has the most favorable formation probability, especially in the presence of additional water molecules. This product may serve as a preceding structure to the 8-oxo-7,8-dihydroguanine lesion in DNA and has implications for health effects of radiation exposure and radiation therapy.

Ab initio CCSD(T)/CBS//ωB97X-D/6-311+G(d,p) calculations of the C6H6N potential energy surface were performed to investigate the reaction mechanism underlying the reaction of atomic nitrogen (2D) with benzene. Thereafter, Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of reaction rate constants and product branching ratios were performed under single-collision conditions. The results revealed that the N(2D) + C6H6 reaction in the case of statistical behavior is expected to produce hydrogen cyanide plus a cyclopentadienyl radical (91.5-88.9%), acetylene plus a pyrrole radical (5.8-7.5%), 1-cyano-2,4-cyclopentadiene + H (2.3-3.0%) and 1-ethynyl-pyrrole + H (0.4-0.6%), with the most favorable pathways being the initial adduct i1 leading to the formation of a seven-membered cyclic intermediate i12 through an exothermic ring expansion process and a multistep route i12 → i15 → i16 → C5H5 + HCN featuring an intramolecular ring-shrinking process involving a C-C bond fusion elimination channel to yield the bicyclic intermediate i15, followed by hydrogen cyanide elimination, thus forming a cyclopentadienyl radical. The calculated product branching ratios were consistent with the available experimental data; however, some quantitative deviations from the experimental results and the possible reasons are also discussed. The possible effects of the title reaction on the upper atmosphere of Titan, with critical implications for the rapid degradation of nitrogen-bearing polycyclic aromatic hydrocarbons, were compared with the mass growth processes of their polycyclic aromatic hydrocarbon counterparts produced through ring expansion.

The reaction, BrO + HO2 → HOBr + O2, is exothermic and can produce O2 in both its ground state (X3∑g-) and its first excited state (ã1Δg). As a result, this reaction can proceed on both a singlet and a triplet potential energy surface. Recently, Tsona, Tang and Du published a paper entitled "Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products (Phys. Chem. Chem. Phys. 2019, 21, 20296-203072). The results of this work showed significant differences from those published earlier on this reaction by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569). Further calculations performed in this present work, combined with higher level calculations published by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569), demonstrate that the work of Tsona et al. is flawed because the integration grid size used in their lowest singlet and triplet calculations is too small, and a closed-shell wavefunction, rather than an open-shell wavefunction, has been used for the singlet surface. The major conclusion in the work of Tsona et al. that the lowest singlet and triplet channels are barrierless is shown to be incorrect. Also, the computed rate coefficients of Tsona et al. showed a positive temperature dependence, which is inconsistent with the experimentally observed negative temperature dependence, whereas the singlet rate coefficients computed by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569) showed a negative temperature dependence consistent with experiment.

In this work, the detailed reaction mechanism of the astrochemically relevant exoergic and barrierless, H + LiH+ → H2 + Li+ , reaction is investigated by both time-dependent wave packet and quasi-classical trajectory (QCT) methods on the ab initio electronic ground state potential energy surface of Martinazzo et al. [Martinazzo et al., J. Chem. Phys., 2003, 119, 11241]. The interference terms due to the coherence between the partial waves are quantified. When plotted along the scattering angle they reveal interference of constructive or destructive nature. Significant interference found in the differential cross section (DCS) which is a symbolic of the non-statistical nature of the reaction. This is further complemented by calculating the lifetime of the collision complex by the QCT method. It is found that the reaction follows a direct stripping mechanism at higher collision energies and yields forward scattered products from collisions involving high total angular momentum. At low collision energies, the reaction follows a mixed direct/indirect mechanism but with a dominant indirect contribution. The product state-resolved DCSs reveal that two opposite mechanisms co-exist, both at low and high collision energies. The microscopic scattering mechanism of the reaction is found to be unaffected by the ro-vibrational excitation of the reagent diatom.

Ab initio G3(MP2,CC)/B3LYP/6–31G** calculations of the potential energy surface for the reaction of 3-acenaphthyl radical with acetylene were combined with Rice-Ramsperger-Kassel-Marcus Master Equation calculations of the temperature- and pressure-dependent reaction rate constants and product branching ratios with the goal to examine the transformation of an edge five-membered ring into a six-membered one. The reaction mechanism has been shown to contain two stages: (1) the formation of 3-ethynylacenaphthylene + H in the 3-acenaphthyl + C2H2 reaction via acetylene addition to the radical site and an H loss, (2) the formation of the 1-phenanthrenyl radical in the 3-ethynylacenaphthylene + H reaction via H addition to the five-membered ring followed by the ring flip. The results of kinetic modeling show this pathway being faster than the pathway proposed earlier featuring the CCH chemisorption on the six-membered ring without involvement of the five-membered ring. A kinetic test demonstrated the dominance of the former pathway for the formation of pyrene from naphthalene. It is expected that the new mechanism plays a major role in the transformation of five-membered rings on edges of PAH structures.

Reactive molecular dynamics (MD) simulation is a powerful tool to study the reaction mechanism of complex chemical systems. Central to the method is the potential energy surface (PES) that can desc...

The N2H3 + NO2 reaction plays a key role during the early stages of hypergolic ignition between N2H4 and N2O4. Here for the first time, the reaction kinetics of N2H3 in excess NO2 was studied in 2.0 Torr of N2 and in the narrow temperature range 298-348 K in a pulsed photolysis flow-tube reactor coupled to a mass spectrometer. The temporal profile of the product, HONO, was determined by direct detection of the m/z +47 amu ion signal. For each chosen [NO2], the observed [HONO] trace was fitted to a biexponential kinetics expression, which yielded a value for the pseudo-first-order rate coefficient, k', for the reaction of N2H3 with NO2. The slope of the plot of k' versus [NO2] yielded a value for the observed bimolecular rate coefficient, kobs, which could be fitted to an Arrhenius expression of (2.36 ± 0.47) × 10-12 exp((520 ± 350)/T) cm3 molecule-1 s-1. The errors are 1σ and include estimated uncertainties in the NO2 concentration. The potential energy surface of N2H3 + NO2 was investigated by advanced ab initio quantum chemistry theories. It was found that the reaction occurs via a complex reaction mechanism, and all of the reaction channels have transition state energies below that of the entrance asymptote. The radical-radical addition forms the N2H3NO2 adducts, while roaming-mediated isomerization reactions yield the N2H3ONO isomers, which undergo rapid dissociation reactions to several sets of distinct products. The RRKM multiwell master equation simulations revealed that the major product channel involves the formation of trans-HONO and trans-N2H2 below 500 K and the formation of NO + NH2NHO above 500 K, which is nearly pressure independent. The pressure-dependent rate coefficients of the product channels were computed over a wide pressure-temperature range, which encompassed the experimental data.

The potential energy surface for the CH2CCl + O2 reaction has been investigated by using the CCSD(T)/cc-pVTZ//B3LYP/6-311++G(d,p) method. Two type reaction mechanisms have been located. The H-abstraction of CH2CCl by O2 generates CHCCl + HO2 surmounting a 20.86 kcal/mol barrier. The addition between O2 and CH2CCl proceeds to an intermediate CH2CClO2 (IM1t and IM1c) without a barrier, which can further dissociate or isomerize to various products with the complicated processes. The temperature and pressure dependence rate constants for the CH2CCl + O2 reaction were computed by means of multi-channel RRKM-TST theory. Moreover, TDDFT calculations imply that IM1t, IM1c, IM2, IM4, IM5t, and IM5c will photolyze under the sunlight.

The reaction mechanism, product branching ratios, and relevant rate constants for the reaction of imidogen (NH) with sulfur monoxide (SO) over singlet and triplet potential energy surfaces are theoretically investigated. Various quantum chemical methods at the single-reference methods (PBE, M06-2X, MP2, GBS-QB3, G3MP2B3, and CCSD(T)) and the multi-reference methods of CASPT2 are carried out to examine the characteristics of the title reaction's potential energy surface. Eighteen chemically activated intermediates and more than 35 different reaction channels are predicted over the singlet surface, while less species and channels are distinguished over the triplet surface. The entrance channels for both surfaces appeared to be barrier-less association reactions to form pre-reaction energized intermediates of singlet or triplet HNSO or HNOS. OH and NS radicals are indicated as the major products for the title reaction on both surfaces in agreement with the reported experimental observations. The RRKM-steady state approximation method is used to calculate the rate constants and branching ratios of the main products. The obtained overall rate constant is in agreement with the available reported experimental data over the wide range of temperature from 300 to 3000 K. By considering single-reference calculations, the singlet and triplet total rate constants were found to be k(T) = 5.04 × 1010 and 2.47 × 1012 T-0.83 exp(-1.56 kJ mol-1/T), respectively. Also, the total rate constant for the consumption of reactants by inclusion of multi-reference calculations was found to be in the range of 3.86 × 1010 to 4.18 × 1010, depending on the level of calculations. In addition, our results revealed that the total rate constant for the NH + SO reaction is pressure-independent in the range of 0.1-2000 Torr.

Ab initio CCSD(T)/CBS//B3LYP/6-311G(d, p) calculations of the potential energy surface for possible dissociation channels of HOC $_2$ H $_3$ F, as well as Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of rate constants, were carried out, in order to predict statistical product branching ratios in dissociation of HOC $_2$ H $_3$ F at various internal energies. The most favorable reaction pathway leading to the major CH $_2$ CHO+HF products is as the following: OH+C $_2$ H $_3$ F $\rightarrow$ i2 $\rightarrow$ TS14 $\rightarrow$ i6 $\rightarrow$ TS9 $\rightarrow$ i3 $\rightarrow$ TS3 $\rightarrow$ CH $_2$ CHO+HF, where the rate-determining step is HF elimination from the CO bridging position via TS11, lying above the reactants by 3.8 kcal/mol. The CH $_2$ O+CH $_2$ F products can be formed by F atom migration from C $_\beta$ to C $_\alpha$ position via TS14, then H migration from O to C $_\alpha$ position via TS16, and C-C breaking to form the products via TS5, which is 1.8 kcal/mol lower in energy than the reactants, and 4.0 kcal/mol lower than TS11.

Ethylene polymerization mediated by methyl-bis(cyclopentadienyl)-zirconium or zirconocene catalyst, [ZrCp₂CH₃]⁺, is one of the most popular catalytic reaction for polyethylene production. Rationalizing the major effects that control the polymer growth result in a challenge for computational studies. Through quantum chemical calculations, we characterized the zirconocene ethylene polymerization reaction mechanism: chain initiation (I; first ethylene insertion) [ZrCp₂CH₂CH₂CH₃]⁺, chain propagation (P; from second (P₁) to ninth (P₉) ethylene insertion) [ZrCp₂ (CH₂)₂₀CH₃]⁺, and chain termination processes (T; β-hydrogen elimination from P₅ or P₉) [ZrHCp₂ (H₂C═CH(CH₂)₁₈CH₃]⁺ are analyzed through the potential energy surface (PES) and reaction force analysis (RFA). The RFA approach involves pulling out the portion of an activation barrier that corresponds to distorting reactants into the geometries they adopt in a transition state structure until it reaches the structural relaxation toward the equilibrium geometry of the product. Because the interactions between the zirconocene and the ethylene molecule are influenced by a combination of several kinds of steric and electronic effects, it is indispensable to understanding these interactions in order to rationalize and predict in a quantitative manner the reaction barrier heights and the concomitant polymer growth. In the present work, we employ a simple procedure within the framework of the RFA and the density functional steric energy decomposition analysis (EDA) approach to quantitatively separate the different types of interactions; steric (ΔE_s), electrostatic (ΔE_e), and quantum (ΔE_q) effects in order to predict the impact of each factor on the course of the polymerization process as well as for the polymer control and design.

The reaction mechanism of HNCO with HO2 radicals is investigated by means of the B3PW91/6-311+G(d,p) method to determine a more reasonable pathway. Four possible entrance patterns are designed; however, only three situations are finally confirmed. When HNCO and HO2 are close to each other, they first form an intermediate a2. Then, the terminal O atom in HO2 connects with the C atom of HNCO along with the H atom transferring from HO2 to the O atom of HNCO. After that the C–O bond ruptures again to form i. Finally, product P6(NCO + HOOH) is generated from i via H atom transfer. P6 is the most accessible product with the simplest steps and lowest barrier height. In addition, the energy of several of the routes is refined at the CCSD(T)/6-311+G(d,p) level based on the optimized geometries. Although there are some differences, the most favorable product is still P6. It is expected that this study would be helpful in further studies on HNCO and in understanding its reactions.

The reaction mechanism between phosphacyclopropenylidene and ethylene has been systematically investigated at the B3LYP/6-311++G(d,p) level of theory in order to better understand the reactivity of unsaturated cyclic phosphorus-bearing carbene. Geometry optimizations and vibrational analyses have been performed for the stationary points on the potential energy surface of the system. Calculations show that the spiro bicyclic intermediate could be produced through the cycloaddition process between phosphacyclopropenylidene and ethylene initially. The reaction mechanism is illustrated with the frontier molecular orbital theory. Introduction of electron-withdrawing group in phosphacyclopropenylidene will better facilitate the addition process. Through subsequent ring-expanding and hydrogen-migrating process, fuse-ring and allene compounds could be produced, respectively. Furthermore, it?s easy for spiro bicyclic intermediate and another ethylene to form a spiro tricyclic compound. This study is helps to understand the reactivity of phosphacyclopropenylidene, the evolution of phosphorus-bearing molecules in space, and to offer an alternative approach to the formation of phosphorus-bearing heterocyclic compound.

One of the fundamental goals of chemistry is to determine how molecular structure influences interactions and leads to different reaction products. Studies of isomer-selected and resolved chemical reactions can shed light directly on how form leads to function. In the following, we present the results of gas-phase reactions between acetylene cations (C2D2+) with two different isomers of C3H4: propyne (DC3D3) and allene (H2C3H2). Our highly controlled, trapped-ion environment allows for precise determination of reaction products and kinetics. From these results, we can infer details of the underlying reaction dynamics of C2H2+ + C3H4. Through the synergy of experimental results and high-level quantum chemical potential energy surface calculations, we are able to identify distinct reaction mechanisms for the two isomers. We find long-range charge exchange with no complex formation is favored for allene, whereas charge exchange leads to an intermediate reaction complex for propyne and thus, different products. Therefore, this reaction displays a pronounced isomer-selective bi-molecular reactive process.

The relative stability and kinetics of CH4 … O2 … MO and the reaction potential energy surface under the transition metal oxide catalysts (RuO/RhO/PdO) during methane oxidation reaction were studie...

The mechanism of the reaction of lithium and magnesium carbenoids with ethylene to give cyclopropane has been explained in detail in all the steps at the G4 level of theory. We explored the lithium and magnesium interaction toward πC=C and σC-C bonds in the reactants and the products. We have also investigated the reaction path by means of the force profile formalism in order to highlight the electronic and the structural rearrangements along the potential energy surface of the cyclopropanation. The results indicate that all of the reactions are stepwise, exoenergetic, with low barriers. All our findings were confirmed by dynamic simulations for chlorometal carbenoids. Furthermore, from the intrinsic reaction coordinate procedure, we were able to find out the intermediates that can take place when the reaction is descending from the transition state to the products or reactants. The reaction force analysis at B3LYP/6-311G(d,p) indicates that the energy barriers are mostly due to structural rearrangements which are produced by the approach of the carbenoid to ethylene. Quantum theory of atoms in molecules and electron localization function analyses indicate that products, reactants, and intermediates form complexes stabilized by attractive forces between Li/Mg and single/double bonds.

Recent guided ion beam experiments have revealed interesting reaction dynamics of the HBr+ + CO2 → HOCO+ + Br· reaction under different conditions. The hypothesis is that the predominant reaction mechanism depends on the collision energy between two reactants, the angular momentum of HBr+, and the spin-orbit coupling state of the system. The potential energy profile of the HBr+ + CO2 → HOCO+ + Br· reaction is studied in this research to lay the groundwork for an ab initio molecular dynamics simulation. First, a benchmark potential energy profile of this reaction is identified using coupled-cluster theory extrapolated to the complete basis set limit. A transition state connecting the previously reported intermediates is found, making the potential energy surface of the HBr+ + CO2 → HOCO+ + Br· reaction double-welled. Second, various single reference ab initio methods are compared with the benchmark potential energy profile to search for the most suitable ab initio method for the dynamics simulation. Two combinations of double-ζ basis sets (with effective core potentials) with MP2 and density functional theory have been identified to accurately represent the potential energy profile of this reaction.

A full-dimension quasi-classical trajectory study of collisions between two NH radicals is presented. Interatomic interactions are represented by a previously reported global six-dimensional potential energy surface for singlet electronic state of the N2H2 system. This study suggests that the formation of N2 from the collision of two NH radicals may occur via a one-step (NH + NH → N2 + H + H) or two-step (NH + NH → N2H + H → N2 + H + H) microscopic reaction mechanism. A fast vibrational energy redistribution is observed in the four-body complex in the latter mechanism. Excitation functions are presented and discussed. A variant of the vibrational energy quantum mechanical threshold method was used to correct the zero-point energy leakage in the classical calculations. The influence of reactant's rotational and vibrational energy on reactivity was investigated by state-specific calculations with one or both reactants excited. Reaction rate constants for the ground and some rotationally excited states are presented using an Arrhenius-Kooij-like functional form.

The mechanism and kinetics of the reaction of ortho-benzyne with vinylacetylene have been studied by ab initio and density functional CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the pertinent potential energy surface combined with Rice-Ramsperger-Kassel-Marcus - Master Equation calculations of reaction rate constants at various temperatures and pressures. Under prevailing combustion conditions, the reaction has been shown to predominantly proceed by the biradical acetylenic mechanism initiated by the addition of C4 H4 to one of the C atoms of the triple bond in ortho-benzyne by the acetylenic end, with a significant contribution of the concerted addition mechanism. Following the initial reaction steps, an extra six-membered ring is produced and the rearrangement of H atoms in this new ring leads to the formation of naphthalene, which can further dissociate to 1- or 2-naphthyl radicals. The o-C6 H4 +C4 H4 reaction is highly exothermic, by ∼143 kcal/mol to form naphthalene and by 31-32 kcal mol-1 to produce naphthyl radicals plus H, but features relatively high entrance barriers of 9-11 kcal mol-1 . Although the reaction is rather slow, much slower than the reaction of phenyl radical with vinylacetylene, it forms naphthalene and 1- and 2-naphthyl radicals directly, with their relative yields controlled by the temperature and pressure, and thus represents a viable source of the naphthalene core under conditions where ortho-benzyne and vinylacetylene are available.

We reply to the comment on our recent paper entitled “Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products” by Chow et al. In their comment, the authors raised the differences between our results and their results in an earlier paper (R. Chow, D. K. W. Mok, E. P. F. Lee and J. M. Dyke, Phys. Chem. Chem. Phys., 2016, 18, 30554–30569), in terms of kinetics and potential energy surface, and they attributed these differences to the use of a small integration grid size and closed-shell wavefunctions for geometry optimizations in our study. Indeed, in our original manuscript, we did not ensure the proper use of UHF wavefunctions for singlet states, which led the singlet states to be treated with restricted M06-2X wavefunctions during optimizations. Furthermore, the default integration grid was used. New geometry optimizations have been performed where reactant complexes on the singlet surface were treated in their open-shell singlet states (ensured by using unrestricted-spin wave-functions) and using very tight convergence criteria, and new reaction rate constants have been calculated based on new energy barriers. No barrierless hydrogen abstraction reactions were observed as reported in our previous results and, consequently, the outer rate coefficient in the two-transition state approach (given by eqn (5) in Tsona et al., 2019) was determined by the collision theory. Overall rate constants exhibit a negative temperature dependency in the 200–400 K range. Despite the changes on the reaction energies and kinetics due to wrong UHF wavefunctions, our main conclusion that water has no net effect on the BrO + HO2 → BrOH + O2 reaction is still valid.

The Diels–Alder cycloaddition of 3-bromo-1-phenylprop-2-ynone with furan and 2-methylfuran has been studied using the molecular electron density theory through DFT calculation at the B3LYB/6-311G (d,p) level of theory. The mechanism of regioselectivity for these reactions was investigated by the potential energy surface determination for the cycloaddition process. The regiostereoselectivity was studied by exploiting the local and global reactivity descriptors defined by Fukui and Parr functions; the results showed that the diene acts as a nucleophile, by two possible approaches, namely furan and 2-methylfuran, and the dienophile behaves as an electrophile. The prediction of the possible stereo-/regioisomers for this reaction is revealed by different preferable pathways of attack. The reaction mechanism has been followed by optimizing the transition state structure using the intrinsic coordinate IRC computations and the estimated activation energies. The analysis geometries of the transition states prove that the cycloaddition is established in one step according to an asynchronous mechanism for the various rotamers.