Examining the robustness of first-principles calculations for metal hydride reaction thermodynamics by detection of metastable reaction pathways

First principles calculations have played a useful role in screening mixtures of complex metal hydrides to find systems suitable for H(2) storage applications. Standard methods for this task efficiently identify the lowest energy reaction mechanisms among all possible reactions involving collections of materials for which DFT calculations have been performed. The resulting mechanism can potentially differ from physical reality due to inaccuracies in the DFT functionals used, or due to other approximations made in estimating reaction free energies. We introduce an efficient method to probe the robustness of DFT-based predictions that relies on identifying reactions that are metastable relative to the lowest energy reaction path predicted with DFT. An important conclusion of our calculations is that in many examples DFT cannot unambiguously predict a single reaction mechanism for a well defined metal hydride mixture because two or more mechanisms have reaction energies that differ by a small amount. Our approach is illustrated by analyzing a series of single step reactions identified in our recent work that examined reactions with a large database of solids [Kim et al., Phys. Chem. Chem. Phys. 2011, 13, 7218].     

Quantum Chemical Approach for Determining Degradation Pathways of Phenol by Electrical Discharge Plasmas

This study uses density functional theory (DFT) simulations to predict the main pathways by which hydroxyl (OH) radicals oxidize phenol into monohydroxylated products during an electrical discharge directly in or contacting water. The calculated activation energies and reaction rate constants indicate that phenol ring H abstraction is less likely to occur than OH addition, which will be the fastest in the ortho and para positions. The chain propagation with molecular oxygen of such formed ortho and para radicals will result in the production of hydroquinone and catechol, which are, concurrently, the most likely products of phenol degradation by OH radicals. Electron transfer reactions between dihydroxycyclohexadienyl radicals and plasma oxidative species are another important reaction mechanism which may be contributing significantly to the formation of products. Good agreement between computed kinetic and experimental data demonstrates the feasibility of applying DFT to investigate chemical reaction mechanisms.     

OH-initiated oxidation mechanism and kinetics of organic sunscreen benzophenone-3: A theoretical study

Benzophenone-3 (BP-3), as an important organic UV filter, is widely used in the sunscreen, cosmetic, and personal care products. The chemical reaction mechanism and kinetics of BP-3 degradation initiated by hydroxyl (OH) radical was investigated in the atmosphere based on the density functional theory (DFT). The results showed that the OH radical is more easily added to the C3 position of the aromatic ring (pathway 3), while the H atom abstraction from the OH group on the aromatic ring (pathway 23) is an energetically favorable reaction pathway. At ambient temperature, 298 K, the overall rate constant for the primary reaction is about 1.50 × 10^?10 cm³ mol^?1 s^?1 with the lifetime of 1.92 h. OH addition reactions play the key role in the OH-initiated reaction of BP-3. The study is helpful for better understanding of the removal, transformation, and fate of BP-3 in the atmosphere.     

Hydroxyl radicals in ice: insights into local structure and dynamics

The hydroxyl radical and its reactivity within ice environments are crucial to many important atmospheric reactions. The associated molecular mechanisms are largely unknown due to challenges posed by direct experimental measurements and computational studies of this transient species. Here we report insights into the local structure and behaviour of the hydroxyl radical in bulk ice through an extensive study utilizing Car-Parrinello molecular dynamics simulations. Interstitial and in-lattice hydroxyl radicals in hexagonal ice were investigated at primarily 190 K. Our findings, utilizing both HCTH/120 and BLYP functionals, show that OH* can exhibit greater mobility than other ice defects (the trapping energy estimated to be only 0.09 eV). We observe the formation of a two-center three-electron hemibond structure between the hydroxyl radical and an in-lattice water molecule; while controversial, such a structure in ice may be amenable to experimental detection due to its relative stability. Our results show that interstitial water molecules can strongly influence the mobility of the hydroxyl radical in bulk ice through the displacement of the radical to an interstitial location. We also demonstrate that the H-transfer reaction from an interstitial water to the radical is a rare event in ice. Together, these results predict that the radical can be a reactive species in bulk ice, as both interstitial and in-lattice OH* can be available for reactions with other species. These microscopic insights should contribute to our understanding of the reactivity of OH* in ice and its implications to atmospheric reactions.     

Theoretical study on the reaction mechanism of the OH-initiated oxidation of CH<Subscript>2</Subscript>=C(CH<Subscript>3</Subscript>)CH<Subscript>2</Subscript>CH<Subscript>2</Subscript>OH

In the present work, the detailed reaction mechanism and possible products of the OH-initiated oxidation of CH₂=C(CH₃)CH₂CH₂OH (MBO331) have been revealed theoretically for the first time. The potential energy surfaces of various reaction channels both in the absence and presence of O₂ and NO are evaluated at the CCSD(T)/6−31++G(d,p)//MP2(full)/6−311G(d,p)+ZPE*0.95 level. The major products of HCHO + CH₃C(O)CH₂CH₂OH predicted for the title reaction in the presence of O₂ and NO are in agreement with those of similar reactions of unsaturated alcohols with OH radical.     

Mechanism for degradation of Nafion in PEM fuel cells from quantum mechanics calculations

We report results of quantum mechanics (QM) mechanistic studies of Nafion membrane degradation in a polymer electrolyte membrane (PEM) fuel cell. Experiments suggest that Nafion degradation is caused by generation of trace radical species (such as OH(●), H(●)) only when in the presence of H(2), O(2), and Pt. We use density functional theory (DFT) to construct the potential energy surfaces for various plausible reactions involving intermediates that might be formed when Nafion is exposed to H(2) (or H(+)) and O(2) in the presence of the Pt catalyst. We find a barrier of 0.53 eV for OH radical formation from HOOH chemisorbed on Pt(111) and of 0.76 eV from chemisorbed OOH(ad), suggesting that OH might be present during the ORR, particularly when the fuel cell is turned on and off. Based on the QM, we propose two chemical mechanisms for OH radical attack on the Nafion polymer: (1) OH attack on the S-C bond to form H(2)SO(4) plus a carbon radical (barrier: 0.96 eV) followed by decomposition of the carbon radical to form an epoxide (barrier: 1.40 eV). (2) OH attack on H(2) crossover gas to form hydrogen radical (barrier: 0.04 eV), which subsequently attacks a C-F bond to form HF plus carbon radicals (barrier as low as 1.00 eV). This carbon radical can then decompose to form a ketone plus a carbon radical with a barrier of 0.86 eV. The products (HF, OCF(2), SCF(2)) of these proposed mechanisms have all been observed by F NMR in the fuel cell exit gases along with the decrease in pH expected from our mechanism.     

Theoretical explanation of nonexponential OH decay in reactions with benzene and toluene under pseudo-first-order conditions

OH radical reactions with benzene and toluene have been studied in the 200-600 K temperature range via the CBS-QB3 quantum chemistry method and conventional transition-state theory. Our study takes into account all possible hydrogen abstraction and OH-addition channels, including ipso addition. Reaction rates have been obtained under pseudo-first-order conditions, with aromatic concentrations in large excess compared to OH concentrations, which is the case in the reported experiments as well as in the atmosphere. The reported results are in excellent agreement with the experimental data and reproduce the discontinuity in the Arrhenius plots in the 300 K < T < 400 K temperature range. They support the suggestion that the observed nonexponential OH decay is caused by the existence of competing addition and abstraction channels and by the decomposition of thermalized OH-aromatic adducts back to reactants. We also find that the low-temperature onset of the nonexponential decay depends on the concentration of the aromatic compounds and that the lower the concentration, the lower the temperature onset. Under atmospheric conditions, nonexponential decay was found to occur in the 275-325 K range, which corresponds to temperatures of importance in tropospheric chemistry. Branching ratios for the different reaction channels are reported. We find that for T > or = 400 K the reaction occurs exclusively by H abstraction. At 298 K, ipso addition contributes 13.0% to the overall OH + toluene reaction, while the major products correspond to ortho addition, which represents 43% of all possible channels.     

Computation on IR of Radical Intermediates of Isoprene Reacting with OH Radical

The radical reactions of isoprene, the most abundant natural volatile organic compounds(VOC), is important to understand the atmospheric activities of isoprene and to evaluate the role of VOC in atmospheric pollution. Isoprene reaction with OH radical is such important radical reaction as its contribution to the isoprene decomposing in natural atmosphere. Meantime its radical products also make big contribution to other atmospheric reactions.     

Trapping of the OH radical by alpha-tocopherol: a theoretical study

The antioxidant activity of alpha-tocopherol against the damaging hydroxyl radical was analyzed theoretically by hybrid density functional theory, following the direct dynamics method, where the thermal rate constants were calculated using variational transition-state theory with multidimensional tunneling. We found that the OH radical is only slightly or not at all selective, attacking by different mechanisms at several positions of the alpha-tocopherol molecule, giving competitive reactions. The most favorable pathways are the hydrogen abstraction reaction from the phenolic hydrogen and the electrophilic addition onto the aromatic ring. We propose a final rate constant, the sum of the competitive hydrogen abstraction and addition reactions, > or =2.7 x 10(8) M(-1) s(-1) at 298 K, where the hydrogen abstraction reaction represents only 20% of the total OH radical reaction. This result indicates that, molecule by molecule, in an apolar environment, alpha-tocopherol is less effective than coenzyme Q (which presents a rate constant of 6.2 x 10(10) M(-1) s(-1) at 298 K) as a scavenger of OH radicals. It was also found that both mechanisms are not direct but pass through intermediates in the entry channel, with little or no influence on the dynamics of the reactions. The hydrogen abstraction reaction also presents another intermediate in the exit channel, which may have a significant role in preventing the pro-oxidant effects of alpha-tocopherol, although less important than with free radicals other than OH.     

过氧亚硝酸与苯酚的反应机理理论研究

采用量子化学密度泛函理论(DFT)研究了过氧亚硝酸分解产生的自由基(•OH和•NO2)与苯酚的反应机理. 在B3LYP/6-311++G(d, p)//B3LYP/6-311G(d, p)水平上对该反应体系的反应物、中间体、过渡态及产物进行了几何构型优化并计算了振动频率和能量. 计算结果表明, 过氧亚硝酸与苯酚的反应生成两种主要产物, 分别为邻羟基苯酚和对羟基苯酚, 这一结论与实验结果一致. 此外在同一计算水平上采用SCRF(PCM)方法计算了溶剂化效应, 结果表明, 极性溶剂可以降低各反应通道的活化能, 有利于反应的进行.     

Direct Dynamics Simulation of Reaction Between F2 and Ethylene

Direct dynamics within the framework of DFT was used to study the long-time puzzling mechanism of the reaction between F2 and ethylene. Three types of reactions are widely accepted : F atom elimination reaction, HF elimination reaction and the addition reaction. Several reaction mechanisms have been proposed, but only the radical mechanism can reasonably explain the initial reaction at low temperature. In this article, our calculations support the radical mechanism and the reaction mechanisms of the three reactions, and they are described in detail by trajectory simulation. The reactions in a cryogenic matrix with the reaction mechanism were also discussed.     

Isoguanine formation from adenine

Several possible mechanisms underlying isoguanine formation when OH radical attacks the C(2) position of adenine (A?C?2) are investigated theoretically for the first time. Two steps are involved in this process. In the first step, one of two low-lying A?C?2???OH reactant complexes is formed, leading to C(2)-H(2) bond cleavage. Between the two reactant complexes there is a small isomerization barrier, which lies well below separated adenine plus OH radical. The complex dissociates to free molecular hydrogen and an isoguanine tautomer (isoG?1 or isoG?2). The local and activation barriers for the two pathways are very similar. This evidence suggests that the two pathways are competitive. After dehydrogenation, there are two possible routes for the second step of the reaction. One is direct hydrogen transfer, via enol-keto tautomerization, which has high local barriers for both tautomers and is not favored. The other option is indirect hydrogen transfer involving microsolvation by one water molecule. The water lowers the reaction barrier by over 20?kcal mol(-1) , indicating that water-mediated hydrogen transfer is much more favorable. Both A+OH(?) →isoG+H(?) reactions are exothermic and spontaneous. Among four isoguanine tautomers, isoG?1 has the lowest energy. Our findings explain why only the N(1)H and O(2)H tautomers of isolated isoguanine and isoguanosine have been observed experimentally.     

Photochemistry and Reactivity of the Phenyl Radical–Water System: A Matrix Isolation and Computational Study

The reaction of the phenyl radical 1 with water has been investigated by using matrix isolation spectroscopy and quantum chemical calculations. The primary thermal product of the reaction between 1 and water is a weakly bound complex stabilized by an OH???π interaction. This complex is photolabile, and visible‐light irradiation (λ>420 nm) results in hydrogen atom transfer from water to radical 1 and the formation of a highly labile complex between benzene and the OH radical. This complex is stable under the conditions of matrix isolation, however, continuous irradiation with λ>420 nm light results in the complete destruction of the aromatic system and formation of an acylic unsaturated ketene. The mechanisms of all reaction steps are discussed in the light of ab initio and DFT calculations.