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

Ultralow-Friction at Cryogenic Temperature Induced by Hydrogen Correlated Quantum Effect

Temperature is one of the governing factors affecting friction of solids. Undesired high friction state has been generally reported at cryogenic temperatures due to the prohibition of thermally activated processes, following conventional Arrhenius equation. This has brought huge difficulties to lubrication at extremely low temperatures in industry. Here, the study uncovers a hydrogen-correlated sub-Arrhenius friction behavior in hydrogenated amorphous carbon (a-C:H) film at cryogenic temperatures, and a stable ultralow-friction over a wide temperature range (103-348 K) is achieved. This is attributed to hydrogen-transfer-induced mild structural ordering transformation, confirmed by machine-learning-based molecular dynamics simulations. The anomalous sub-Arrhenius temperature dependence of structural ordering transformation rate is well-described by a quantum mechanical tunneling (QMT) modified Arrhenius model, which is correlated with quantum delocalization of hydrogen in tribochemical reactions. This work reveals a hydrogen-correlated friction mechanism overcoming the Arrhenius temperature dependence and provides a new pathway for achieving ultralow friction under cryogenic conditions.

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.

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.

Effect of microsolvation on the mode specificity of the OH˙(H2O) + HCl reaction

The present study investigates the mode specificity in the microsolvated OH˙(H₂O) + HCl reaction using on-the-fly direct dynamics simulation. To the best of our knowledge, this is the first study which aims to gain insights into the effect of microsolvation on the mode selectivity. Our investigation reveals that, similar to the gas phase OH˙ + HCl reaction, the microsolvated reaction is also predominantly affected by the vibrational excitation of the HCl mode, whereas the OH vibrational mode behaves as a spectator. Interestingly, in contrast to the behavior of the bare reaction, the integral cross section at the ground state of the microsolvated reaction decreases with an increase in translational energy. However, for the vibrational excited states, the reactivity of the microsolvated reaction is found to be higher than that of the bare reaction within the selected range of translational energies.

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.

Reply to the 'Comment on "Topography of the Free Energy Landscape on the Claisen-Schmidt Condensation: Solvent and Temperature Effect in the Rate-Controlling Step"' by N. D. Coutinho, H. G. Machado, V. H. Carvalho-Silva and W. A. da Silva, Phys. Chem. Chem. Phys., 2021, , 6738

In the Comment on our paper on the description of the Gibbs Free energy profile of Claisen-Schmidt condensation, it is claimed that our calculations are flawed due to inconsistencies with experimental isotope effects in aqueous acetonitrile. In this Reply, we presented rigorous arguments, ambiguities in the Comment and new calculations confirming the consistency of our results: (i) small differences in the relative energetic barriers in both experimental and theoretical curves make the assignment of the rate-limiting step debatable, making the concept of RCS questionable; (ii) it is shown how the misinterpretation of the elementary steps and of the overall processes rate constants led the Comment to incorrect conclusions about the behavior of the inverse isotopic effect; (iii) neglect in the Comment of the inverse kinetic isotope effect in step R2 due to the hybridization conversion, and of the inverse equilibrium isotopic effect for step R1 to describe an overall iKIE > 1, (iv) an erroneous suggestion in the Comment that the disagreement between experimental kinetic parameters is due to the fact that acetonitrile is not used in previous experimental works, when contradictorily the literature recommends it as being indispensable to allow kinetic accuracy; and (v) new calculations improved by explicit-implicit hybrid treatment again ensure that step R4, and not step R5, can assume the role of RCS in protic solvents. Recognizing that questioning is an excellent path for promoting understanding, we hope that the answers provided here will help to clarify and expand the pertinent topics under discussion.

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.

Kinetics Study of the Hydrodeoxygenation of Xylitol over a ReOx-Pd/CeO2 Catalyst

In this study, we elucidate the reaction kinetics for the simultaneous hydrodeoxygenation of xylitol to 1,2-dideoxypentitol and 1,2,5-pentanetriol over a ReOx-Pd/CeO2 (2.0 weight% Re, 0.30 weight% Pd) catalyst. The reaction was determined to be a zero-order reaction with respect to xylitol. The activation energy was elucidated through an Arrhenius relationship as well as non-Arrhenius kinetics. The Arrhenius relationship was investigated at 150–170 °C and a constant H2 pressure of 10 bar resulting in an activation energy of 48.7 ± 10.5 kJ/mol. The investigation of non-Arrhenius kinetics was conducted at 120–170 °C and a sub-Arrhenius relation was elucidated with activation energy being dependent on temperature, and ranging from 10.2–51.8 kJ/mol in the temperature range investigated. Internal and external mass transfer were investigated through evaluating the Weisz–Prater criterion and the effect of varying stirring rate on the reaction rate, respectively. There were no internal or external mass transfer limitations present in the reaction.

From the Kinetic Theory of Gases to the Kinetics of Rate Processes: On the Verge of the Thermodynamic and Kinetic Limits

A variety of current experiments and molecular dynamics computations are expanding our understanding of rate processes occurring in extreme environments, especially at low temperatures, where deviations from linearity of Arrhenius plots are revealed. The thermodynamic behavior of molecular systems is determined at a specific temperature within conditions on large volume and number of particles at a given density (the thermodynamic limit): on the other side, kinetic features are intuitively perceived as defined in a range between the extreme temperatures, which limit the existence of each specific phase. In this paper, extending the statistical mechanics approach due to Fowler and collaborators, ensembles and partition functions are defined to evaluate initial state averages and activation energies involved in the kinetics of rate processes. A key step is delayed access to the thermodynamic limit when conditions on a large volume and number of particles are not fulfilled: the involved mathematical analysis requires consideration of the role of the succession for the exponential function due to Euler, precursor to the Poisson and Boltzmann classical distributions, recently discussed. Arguments are presented to demonstrate that a universal feature emerges: Convex Arrhenius plots (super-Arrhenius behavior) as temperature decreases are amply documented in progressively wider contexts, such as viscosity and glass transitions, biological processes, enzymatic catalysis, plasma catalysis, geochemical fluidity, and chemical reactions involving collective phenomena. The treatment expands the classical Tolman’s theorem formulated quantally by Fowler and Guggenheim: the activation energy of processes is related to the averages of microscopic energies. We previously introduced the concept of “transitivity”, a function that compactly accounts for the development of heuristic formulas and suggests the search for universal behavior. The velocity distribution function far from the thermodynamic limit is illustrated; the fraction of molecules with energy in excess of a certain threshold for the description of the kinetics of low-temperature transitions and of non-equilibrium reaction rates is derived. Uniform extension beyond the classical case to include quantum tunneling (leading to the concavity of plots, sub-Arrhenius behavior) and to Fermi and Bose statistics has been considered elsewhere. A companion paper presents a computational code permitting applications to a variety of phenomena and provides further examples.

"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.