Stereodirectional Origin of anti-Arrhenius Kinetics for a Tetraatomic Hydrogen Exchange Reaction: Born-Oppenheimer Molecular Dynamics for OH + HBr

Formation of Chlorine in the Atmosphere by Reaction of Hypochlorous Acid with Seawater

The highly reactive dihalogens play a significant role in the oxidative chemistry of the troposphere. One of the main reservoirs of these halogens is hypohalous acids, HOX, which produce dihalogens in the presence of halides (Y-), where X, Y = Cl, Br, I. These reactions occur in and on aerosol particles and seawater surfaces and have been studied experimentally and by field observations. However, the mechanisms of these atmospheric reactions are still unknown. Here, we establish the atomistic mechanism of HOCl + Cl- → Cl2 + OH- at the surface of the water slab by performing ab initio molecular dynamics (AIMD) simulations. Main findings are (1) This reaction proceeds by halogen-bonded complexes of (HOCl)···(Cl-)aq surrounded with the neighboring water molecules. (2) The halogen bonded (HOCl)···(Cl-)aq complexes undergo charge transfer from Cl- to OH- to form transient Cl2 at neutral pH. (3) The addition of a proton to one proximal water greatly facilitates the Cl2 formation, which explains the enhanced rate at low pH.

Kinetics, mechanism, and tautomerism in ametryn acid hydrolysis: From molecular structure to environmental impacts

The excessive use of pesticides and the demand for environmentally friendly compounds have driven the focus to detailed studies of the environmental destination of these compounds. Degradation by hydrolysis of pesticides, when released into the soil, can result in the formation of metabolites with potentially adverse effects on the environment. Moving in this direction, we investigated the mechanism of acid hydrolysis of the herbicide ametryn (AMT) and predicted the toxicities of metabolites through experimental and theoretical approaches. The formation of ionized hydroxyatrazine (HA) occurs with the release of the SCH₃- group and the addition of H₃O⁺ to the triazine ring. The tautomerization reactions privileged the conversion of AMT into HA. Furthermore, the ionized HA is stabilized by an intramolecular reaction that provides the molecule in two tautomeric states. Experimentally, the hydrolysis of AMT was obtained under acidic conditions and at room temperature with HA as the main product. HA was isolated in a solid state through its crystallization as organic counterions. The mechanism of conversion of AMT to HA and the experimental investigation of the reaction kinetics allowed us to determine the dissociation of CH₃SH as the rate-controlling step in the degradation process that culminates in a half-life between 7 and 24 months under typical acid soil conditions of the Brazilian Midwest - region with strong agricultural and livestock vocation. The keto and hydroxy metabolites showed substantial thermodynamic stability and a decrease in toxicity compared to AMT. We hope that this comprehensive study will support the understanding of the degradation of s-triazine-based pesticides.

Assessing Thermal Response of Redox Conduction for Anti-Arrhenius Kinetics in a Microbial Cytochrome Nanowire

A micrometers-long helical homopolymer of the outer-membrane cytochrome type S (OmcS) from Geobacter sulfurreducens is proposed to transport electrons to extracellular acceptors in an ancient respiratory strategy of biogeochemical and technological significance. OmcS surprisingly exhibits higher conductivity upon cooling (anti-Arrhenius kinetics), an effect previously attributed to H-bond restructuring and heme redox potential shifts. Herein, the temperature sensitivity of redox conductivity is more thoroughly examined with conventional and constant-redox and -pH molecular dynamics and quantum mechanics/molecular mechanics. A 30 K drop in temperature constituted a weak perturbation to electron transfer energetics, changing electronic couplings (⟨H_mn⟩), reaction free energies (ΔG_mn), reorganization energies (λ_mn), and activation energies (E_a) by at most |0.002|, |0.050|, |0.120|, and |0.045| eV, respectively. Changes in ΔG_mn reflected -0.07 ± 0.03 V shifts in redox potentials that were caused in roughly equal measure by altered electrostatic interactions with the solvent and protein. Changes in intraprotein H-bonding reproduced the earlier observations. Single-particle diffusion and multiparticle steady-state flux models, parametrized with Marcus theory rates, showed that biologically relevant incoherent hopping cannot qualitatively or quantitatively describe electrical conductivity measured by atomic force microscopy in filamentous OmcS. The discrepancy is attributed to differences between solution-phase simulations and solid-state measurements and the need to model intra- and intermolecular vibrations explicitly.

OH + HCl Reaction on the Surface of Ice: An Ab Initio Molecular Dynamics Study

We have investigated the OH + HCl reaction on the surface of ice using Born-Oppenheimer molecular dynamics (BOMD) simulation. The present work revealed that the OH + HCl reaction becomes ∼1 order of magnitude faster on the ice surface compared to the bare reaction. The BOMD simulation also indicates that the Cl radical formed on the ice surface through the title reaction can form two hydrogen bonds at a time with the water molecules present on the ice surface; hence, the Cl radical cannot escape from the ice surface easily.

Using quantum dynamics to study the effect of energy efficiency on the reactivity of the OH + DBr reaction

We report a time-dependent, full dimensional, wave-packet calculation for the reaction of OH + DBr to examine the effect of the energy efficiency on the reactivity. This study shows that the vibrational excitations of the OH and DBr enhance the reaction. However, the rotational excitations of OH and DBr both hinder the reaction. As a result, the vibrational energies of both the OH and DBr reactants are more efficient at promoting the reactivity than the translational energy, while the rotational energies of OH and DBr are less effective than the translational energy. By analyzing the state population of the vibrational and rotational states along the reaction pathway, we also developed an approach in order to explain the enhancement of the vibrational excitation and the hindrance of the rotational excitation of the reaction. We found that the initial-state selected vibrational excited states of OH and DBr are the dominant components, respectively, for surmounting the barrier. However, the initial-state selected rotational excited states of OH and DBr are no longer the dominant states for surmounting the transition state owing to their population changes in the van der Waals well. This quantitative analysis demonstrates the potential well in the entrance valley plays an important role in the energy efficiency with regards to the reactivity.

Aqueous picloram degradation by hydroxyl radicals: Unveiling mechanism, kinetics, and ecotoxicity through experimental and theoretical approaches

Pesticides are chemical compounds widely used to combat pests in crops, and they thus play a key role in agricultural production. However, due to their persistence in aquatic environments, even at low concentrations, their use has been considered an environmental problem and caused concern regarding the adverse effects on human health. This paper reports, for the first time, the mechanisms, kinetics, and an evaluation of the toxicity of picloram degradation initiated by OH radicals in the aqueous environment using quantum chemistry and computational toxicology calculations. The rate constants are calculated using a combination of formulations derived from the Transition State Theory in a realistic temperature range (250-310 K). The results indicate that the two favorable pathways (R1 and R5) of OH -based reactions occur by addition to the pyridine ring. The calculated rate constant at 298 K is compared with the overall second-order reaction rate constant, quantified herein experimentally via the competition kinetics method and data available in the literature showing an excellent agreement. The toxicity assessment and a photolysis study provide important information: i) picloram and the majority of degradation products are estimated as harmful; however, ii) these compounds can suffer photolysis in sunlight. The results of the present study can help understand the mechanism of picloram, also providing important clues regarding risk assessment in aquatic environments as well as novel experimental information.

"Transitivity": A Code for Computing Kinetic and Related Parameters in Chemical Transformations and Transport Phenomena

The Transitivity function, defined in terms of the reciprocal of the apparent activation energy, measures the propensity for a reaction to proceed and can provide a tool for implementing phenomenological kinetic models. Applications to systems which deviate from the Arrhenius law at low temperature encouraged the development of a user-friendly graphical interface for estimating the kinetic and thermodynamic parameters of physical and chemical processes. Here, we document the Transitivity code, written in Python, a free open-source code compatible with Windows, Linux and macOS platforms. Procedures are made available to evaluate the phenomenology of the temperature dependence of rate constants for processes from the Arrhenius and Transitivity plots. Reaction rate constants can be calculated by the traditional Transition-State Theory using a set of one-dimensional tunneling corrections (Bell (1935), Bell (1958), Skodje and Truhlar and, in particular, the deformed ( d -TST) approach). To account for the solvent effect on reaction rate constant, implementation is given of the Kramers and of Collins-Kimball formulations. An input file generator is provided to run various molecular dynamics approaches in CPMD code. Examples are worked out and made available for testing. The novelty of this code is its general scope and particular exploit of d -formulations to cope with non-Arrhenius behavior at low temperatures, a topic which is the focus of recent intense investigations. We expect that this code serves as a quick and practical tool for data documentation from electronic structure calculations: It presents a very intuitive graphical interface which we believe to provide an excellent working tool for researchers and as courseware to teach statistical thermodynamics, thermochemistry, kinetics, and related areas.

Nucleophilic substitution vs elimination reaction of bisulfide ions with substituted methanes: exploration of chiral selectivity by stereodirectional first-principles dynamics and transition state theory

Control of molecular orientation is emerging as crucial for the characterization of the stereodynamics of kinetics processes beyond structural stereochemistry. The special role played in chiral discrimination phenomena has been particularly emphasized by Aquilanti and collaborators after their extensive probes of experimental control of molecular alignment and orientation. In this work, the manifestation of the Aquilanti mechanism has been demonstrated for the first time in first-principles molecular dynamics simulations: stationary points characterized on potential energy surfaces have been calculated for the study of chemical reactions occurring between the bisulfide anion HS^- and oriented prototypical chiral molecules CHFXY (where X = CH₃ or CN and Y = Cl or I). The important reaction channels are those corresponding to bimolecular nucleophilic substitution (S_N2) and to bimolecular elimination (E2): their relative role has been assessed and alternative pathways due to the mirror forms of the oriented chiral molecule are revealed by the different reactivity of the two enantiomers of CHFCNI in S_N2 reaction.

Temperature Dependence of Rate Processes Beyond Arrhenius and Eyring: Activation and Transitivity

Advances in the understanding of the dependence of reaction rates from temperature, as motivated from progress in experiments and theoretical tools (e. g., molecular dynamics), are needed for the modeling of extreme environmental conditions (e.g., in astrochemistry and in the chemistry of plasmas). While investigating statistical mechanics perspectives (Aquilanti et al., 2017b, 2018), the concept of transitivity was introduced as a measure for the propensity for a reaction to occur. The Transitivity plot is here defined as the reciprocal of the apparent activation energy vs. reciprocal absolute temperature. Since the transitivity function regulates transit in physicochemical transformations, not necessarily involving reference to transition-state hypothesis of Eyring, an extended version is here proposed to cope with general types of transformations. The transitivity plot permits a representation where deviations from Arrhenius behavior are given a geometrical meaning and make explicit a positive or negative linear dependence of transitivity for sub- and super-Arrhenius cases, respectively. To first-order in reciprocal temperature, the transitivity function models deviations from linearity in Arrhenius plots as originally proposed by Aquilanti and Mundim: when deviations are increasingly larger, other phenomenological formulas, such as Vogel-Fulcher-Tammann, Nakamura-Takayanagi-Sato, and Aquilanti-Sanches-Coutinho-Carvalho are here rediscussed from the transitivity concept perspective and with in a general context. Emphasized is the interest of introducing into this context modifications to a very successful tool of theoretical kinetics, Eyring's Transition-State Theory: considering the behavior of the transitivity function at low temperatures, in order to describe deviation from Arrhenius behavior under the quantum tunneling regime, a "d-TST" formulation was previously introduced (Carvalho-Silva et al., 2017). In this paper, a special attention is dedicated to a derivation of the temperature dependence of viscosity, making explicit reference to feature of the transitivity function, which in this case generally exhibits a super-Arrhenius behavior. This is of relevance also for advantages of using the transitivity function for diffusion-controlled phenomena.

Quantum dynamics study of kinetic isotope effects of OD with HBr and DBr

Comparison of kinetic isotope effects between quantum dynamics calculations and experiments shows that they agree well with each other both qualitatively and quantitatively.

Kinetics of the OH+HCl→H2 O+Cl reaction: Rate determining roles of stereodynamics and roaming and of quantum tunneling

The OH + HCl → H₂ O + Cl reaction is one of the most studied four-body systems, extensively investigated by both experimental and theoretical approaches. Here, as a continuation of our previous work on the OH + HBr and OH + HI reactions, which manifest an anti-Arrhenius behavior that was explained by stereodynamic and roaming effects, we extend the strategy to understand the transition to the sub-Arrhenius behavior occurring for the HCl case. As previously, we perform first-principles on-the-fly Born-Oppenheimer molecular dynamics calculations, thermalized at four temperatures (50, 200, 350, and 500 K), but this time we also apply a high-level transition-state-theory, modified to account for tunneling conditions. We find that the theoretical rate constants calculated with Bell tunneling corrections are in good agreement with extensive experimental data available for this reaction in the ample temperature range: (i) simulations show that the roles of molecular orientation in promoting this reaction and of roaming in finding the favorable path are minor than in the HBr and HI cases, and (ii) dominating is the effect of quantum mechanical penetration through the energy barrier along the reaction path on the potential energy surface. The discussion of these results provides clarification of the origin on different non-Arrhenius mechanisms observed along this series of reactions. © 2018 Wiley Periodicals, Inc.

A novel assessment of the role of the methyl radical and water formation channel in the CH3OH + H reaction

A number of experimental and theoretical papers accounted almost exclusively for two channels in the reaction of atomic hydrogen with methanol: H-abstraction from the methyl (R1) and hydroxyl (R2) functional groups. Recently, several astrochemical studies claimed the importance of another channel for this reaction, which is crucial for kinetic simulations related to the abundance of molecular constituents in planetary atmospheres: methyl radical and water formation (R3 channel). Here, motivated by the lack of and uncertainties about the experimental and theoretical kinetic rate constants for the third channel, we developed first-principles Car-Parrinello molecular dynamics thermalized at two significant temperatures - 300 and 2500 K. Furthermore, the kinetic rate constant of all three channels was calculated using a high-level deformed-transition state theory (d-TST) at a benchmark electronic structure level. d-TST is shown to be suitable for describing the overall rate constant for the CH₃OH + H reaction (an archetype of the moderate tunnelling regime) with the precision required for practical applications. Considering the experimental ratios at 1000 K, k_R1/k_R2 ≈ 0.84 and k_R1/k_R3 ≈ 15-40, we provided a better estimate when compared with previous theoretical work: 7.47 and 637, respectively. The combination of these procedures explicitly demonstrates the role of the third channel in a significant range of temperatures and indicates its importance considering the thermodynamic control to estimate methyl radical and water formation. We expect that these results can help to shed new light on the fundamental kinetic rate equations for the CH₃OH + H reaction.

Kinetics of low-temperature transitions and a reaction rate theory from non-equilibrium distributions

This article surveys the empirical information which originated both by laboratory experiments and by computational simulations, and expands previous understanding of the rates of chemical processes in the low-temperature range, where deviations from linearity of Arrhenius plots were revealed. The phenomenological two-parameter Arrhenius equation requires improvement for applications where interpolation or extrapolations are demanded in various areas of modern science. Based on Tolman's theorem, the dependence of the reciprocal of the apparent activation energy as a function of reciprocal absolute temperature permits the introduction of a deviation parameter d covering uniformly a variety of rate processes, from those where quantum mechanical tunnelling is significant and d < 0, to those where d > 0, corresponding to the Pareto-Tsallis statistical weights: these generalize the Boltzmann-Gibbs weight, which is recovered for d = 0. It is shown here how the weights arise, relaxing the thermodynamic equilibrium limit, either for a binomial distribution if d > 0 or for a negative binomial distribution if d < 0, formally corresponding to Fermion-like or Boson-like statistics, respectively. The current status of the phenomenology is illustrated emphasizing case studies; specifically (i) the super-Arrhenius kinetics, where transport phenomena accelerate processes as the temperature increases; (ii) the sub-Arrhenius kinetics, where quantum mechanical tunnelling propitiates low-temperature reactivity; (iii) the anti-Arrhenius kinetics, where processes with no energetic obstacles are rate-limited by molecular reorientation requirements. Particular attention is given for case (i) to the treatment of diffusion and viscosity, for case (ii) to formulation of a transition rate theory for chemical kinetics including quantum mechanical tunnelling, and for case (iii) to the stereodirectional specificity of the dynamics of reactions strongly hindered by the increase of temperature.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Quantum dynamics study of energy requirement on reactivity for the HBr + OH reaction with a negative-energy barrier

A time-dependent, quantum reaction dynamics approach in full dimensional, six degrees of freedom was carried out to study the energy requirement on reactivity for the HBr + OH reaction with an early, negative energy barrier. The calculation shows both the HBr and OH vibrational excitations enhance the reactivity. However, even this reaction has a negative energy barrier, the calculation shows not all forms of energy are equally effective in promoting the reactivity. On the basis of equal amount of total energy, the vibrational energies of both the HBr and OH are more effective in enhancing the reactivity than the translational energy, whereas the rotational excitations of both the HBr and OH hinder the reactivity. The rate constants were also calculated for the temperature range between 5 to 500 K. The quantal rate constants have a better slope agreement with the experimental data than quasi-classical trajectory results.

Deformed transition-state theory: Deviation from Arrhenius behavior and application to bimolecular hydrogen transfer reaction rates in the tunneling regime

A formulation is presented for the application of tools from quantum chemistry and transition-state theory to phenomenologically cover cases where reaction rates deviate from Arrhenius law at low temperatures. A parameter d is introduced to describe the deviation for the systems from reaching the thermodynamic limit and is identified as the linearizing coefficient in the dependence of the inverse activation energy with inverse temperature. Its physical meaning is given and when deviation can be ascribed to quantum mechanical tunneling its value is calculated explicitly. Here, a new derivation is given of the previously established relationship of the parameter d with features of the barrier in the potential energy surface. The proposed variant of transition state theory permits comparison with experiments and tests against alternative formulations. Prescriptions are provided and implemented to three hydrogen transfer reactions: CH₄  + OH → CH₃  + H₂ O, CH₃ Cl + OH → CH₂ Cl + H₂ O and H₂  + CN → H + HCN, widely investigated both experimentally and theoretically. © 2016 Wiley Periodicals, Inc.