Next Article in Journal
Inuloxin E, a New Seco-Eudesmanolide Isolated from Dittrichia viscosa, Stimulating Orobanche cumana Seed Germination
Next Article in Special Issue
From the Kinetic Theory of Gases to the Kinetics of Rate Processes: On the Verge of the Thermodynamic and Kinetic Limits
Previous Article in Journal
Chemoenzymatic Total Synthesis of (+)-10-Keto-Oxycodone from Phenethyl Acetate
Previous Article in Special Issue
Supramolecular Fuzziness of Intracellular Liquid Droplets: Liquid–Liquid Phase Transitions, Membrane-Less Organelles, and Intrinsic Disorder
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Hugo G. Machado
1,2,*,
Flávio O. Sanches-Neto
1,2,*,
Nayara D. Coutinho
3,*,
Kleber C. Mundim
2,
Federico Palazzetti
3,* and
Valter H. Carvalho-Silva
1,2,*
1
Grupo de Química Teórica e Estrutural de Anápolis, Centro de Pesquisa e Pós-Graduação. Universidade Estadual de Goiás, 75132-400 Anápolis, GO, Brazil
2
Instituto de Química, Universidade de Brasília, Caixa Postal 4478, 70904-970 Brasília, Brazil
3
Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, 06123 Perugia, Italy
*
Authors to whom correspondence should be addressed.
Molecules 2019, 24(19), 3478; https://doi.org/10.3390/molecules24193478
Submission received: 16 August 2019 / Revised: 4 September 2019 / Accepted: 8 September 2019 / Published: 25 September 2019
(This article belongs to the Special Issue The Fuzziness in Molecular, Supramolecular, and Systems Chemistry)

Abstract

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

Graphical Abstract

1. Introduction

Recent applications of chemical kinetics to a variety of complex systems involves the accurate dealing of properties to be described by techniques, which treat a series of processes beyond elementary chemical quantum dynamics or even approximate classical and semiclassical approaches. We analyze in a companion paper [1] the state of our approaches to these problems from a general viewpoint. Here, we deal with explicit computational calculations that allow moving directly in a simple way to global applications.
Recourse needs to be made at one stage or another to statistical treatments [1,2,3] among which there is modern insurgence with respect to more traditional ones—exemplary is, in particular, this topical collection essentially dedicated to “Fuzzy Logic” [4,5].
Information on the kinetic and related parameters in chemical transformations and transport phenomena and their role in complex mechanisms is needed: Particularly the temperature dependence of rate processes, k ( T ) and often in the low temperature range, where deviations from linearity of Arrhenius plots are revealed. The phenomenology of curvatures in Arrhenius plot span all of chemistry: From the long list that is continuously updated, we refer here to some selected cases, such as combustion chemistry [6], condensed-phase [7], atmospheric and astrochemical reactions [8,9], processes involved in the preservation and aging of food and drugs [10,11], as well as in basic geochemical [12] and biochemical environments [13,14]. The current status of the phenomenology is classified emphasizing case studies, specifically (i) super-Arrhenius kinetics, convex curvature in the Arrhenius plot, where transport phenomena brakes the processes as temperature decreases; (ii) sub-Arrhenius kinetics, concave curvature in the Arrhenius plot, where quantum mechanical tunneling propitiates low temperature reactivity; (iii) anti-Arrhenius kinetics, negative apparent activation energy, where processes are limited by stereodynamic requirements.
The curvature in the Arrhenius plot promotes a temperature dependence on the apparent activation energy E a , the definition of which has been recommended by the International Union for Pure and Applied Chemistry (IUPAC) [13,15] as follows,
E a = k B d ln k ( T ) d   ( 1 T ) = d ln k ( β ) d   β .
In Equation (1), k B is the Boltzmann constant, T is the absolute temperature and β = 1 / k B T . The so-called Tolman’s Theorem [16] proposed an interpretation for E a as the difference between the average over the energy of all reacting systems and energy of all (reactive and nonreactive) systems— propitiating a connection between canonical quantities with microcanonical features of the potential energy surface. Independently, the same spirit lead to transition-state formulations [15,16,17,18]. Recently, it was found expedient to introduce the reciprocal of the apparent activation energy in order to define the Transitivity function [2,3], Equation (2),
γ ( β ) 1 E a ( β )
to construct an appropriate scaling plane, γ ( β )   vs .   β , where the regular curvatures on the Arrhenius plane are approximately linearized: The Transitivity plot. The function γ ( β ) can be interpreted as a measure of the propensity for the reaction to proceed and permit uniformly to account for experimental and theoretical rate processes, such as quantum tunneling, transport properties, and diffusion in the neighborhoods of phase transitions.
For cases where experimental and theoretical rate processes are difficult to interpret—by molecular complexity, extreme conditions or impossibility of exact solution to the Schrodinger’s equation—phenomenological and semiclassical theoretical approaches are of increasing utility. Over the years, the paramount useful phenomenological tool for studying the kinetics of physicochemical processes has been the Arrhenius law. It fails when the temperature range of the rate process becomes large: Empirical laws have been proposed, involving empirical parameters often lacking of physical interpretation, e.g., Kooij [19], power law [20], Vogel–Fulcher–Tamman (VFT) [21,22,23], Nakamura–Takayanagi–Sato (NTS) [24], deformed ( d ) Aquilanti–Mundim (AM) [25] and Aquilanti–Sanchez–Coutinho–Carvalho (ASCC) [26]. In Reference [27], it was presented how the transitivity concept can sort and interpret all these empirical laws, permitting a microscopic interpretation of the phenomenological parameters.
For understanding and predicting a wide variety of kinetic processes with large molecular complexity and presenting deviation from Arrhenius law, Transition-State Theory (TST) remains an excellent formulation. The TST triggered the development of a variety of improved approaches: Important variants include the variational TST [28], the Marcus theory of electron transfer [29], and quantum [30] and path integral [31,32] versions. However, the variants of TST require additional information of the potential energy surface rendering formidable the task of estimating kinetic parameters as the size of molecular systems increases: Examples are the degradation kinetics of organic pollutants in aquatic, soil, and atmospheric environments [33,34,35,36,37,38].
The late twentieth century saw the advent, beyond the Transition State formulation, of molecular dynamics simulations, which is nowadays a powerful theoretical tool in understanding mechanisms of physical and chemical processes. Recently, there is a vast activity evaluating whether advances in molecular dynamics simulations can provide quantitatively rate constants [39,40,41]. However, the extraction of quantitative information on rate constants from molecular dynamics simulations is an important issue but a very difficult one to tackle [42,43,44].
Software tools to calculate and interpret rate constants are enormously useful in material, biochemical, and geochemical research to permit the exploitation of intense progress in computational hardware. Widely diverse codes have been developed to calculate rate constants in a gas-, liquid-, and solid-phase using Transition-State Theory and its variants: Polyrate [45], TheRate [46], MultiWell [47], TAMkin [48], Mesmer [49], RMG [50], APUAMA [51], KiSThelP [52], FRIGUS [53] and Eyringpy [54] codes are excellent options available to work out the kinetics of chemical reactions. Herein, we describe a new code with a user-friendly graphical interface able to perform various procedures for computing kinetic and related parameters in chemical transformations and transport phenomena: The code is referred to as “Transitivity” (see Figure 1), a concept related to the function defined in Equation (2); written in Python, it is a free open-source code compatible with Windows, Linux, and MacOS platforms. It offers the possibility of estimating phenomenological parameters from Arrhenius and Transitivity plots using a stochastic optimization method, Generalized Simulated Annealing (GSA) [55], with several options: Arrhenius, Aquilanti–Mundim (AM) [25], Vogel–Fulcher–Tammann (VFT) [22], Nakamura–Takayanagi–Sato (NTS) [24,56] and Aquilanti–Sanchez–Coutinho–Carvalho (ASCC) [26] formulas. This code also allows the estimation of unimolecular and bimolecular reaction rate constants with traditional TST using Bell35 [57], Bell58 [58], Skodje–Truhlar (ST) [59] and deformed ( d -TST) corrections [60]. Solvent effects can be accounted for by the Collins–Kimball [61] and Kramers [62] models. Finally, input files for different first-principles molecular dynamics—Born–Oppenheimer Molecular Dynamics (BOMD), Car–Parrinello Molecular Dynamics (CPMD), Metadynamics, Path-Integral Molecular Dynamics (PIMD), and Trajectory Surface Hopping (TSH)—can be generated to run by the CPMD code [63].
We propose this code, as we expect it serves as a quick and practical tool for documentation data from electronic structure calculations. Additionally, it presents a very intuitive graphical interface which we believe to provide a useful working tool for the general public and researchers and also as courseware to teach statistical thermodynamics, thermochemistry, kinetics, and related areas. The body of this article is structured in two parts, one that deals with a brief theoretical and operational description of the code, and a second one that shows a variety of examples. A final section is devoted to additional and concluding remarks and an Appendix A collects used symbols and their meaning.

2. Theoretical Background

2.1. Phenomenology of Temperature Dependence of the Reaction Rate Constant

The theoretical apparatus to connect the Transitivity function,   γ ( T ) , and phenomenological reaction rate constant formulas and vice-versa is built and extensively discussed in Reference [27] and references therein. Classical and recent phenomenological reaction rate constant formulas to account for sub-, super-, and anti-Arrhenius behavior are applied in Section 4.1 by use of: (i) Deformed ( d ) Aquilanti–Mundim [25], (ii) VFT (Vogel–Fulcher–Tammann) [22], (iii) NTS (Nakamura–Takayanagi–Sato) [24,56], and (iv) ASCC (Aquilanti–Sanchez–Coutinho–Carvalho) [26]. More details of the formulas are presented in Table 1.

2.2. Calculation of Reaction Rate Constant

The Transition-State Theory, TST, is the most popular tool used to study the kinetics of chemical reactions with a well-defined activated complex—as customary, a double dagger ( ) denotes the properties pertaining to the transition-state complex. For a general bimolecular reaction, such as Reactants   T S Products ,   it is necessary to compute the Q 1 ,       Q 2   and Q partition functions of reactants R 1 and R 2 and of the transition state, respectively. At absolute temperature T , the rate constant is given by:
k T S T ( T ) = k B T h Q Q 1 Q 2   exp ( ε k B T ) ,
where h is Planck’s constant; and ε is the effective height of the energy barrier, eventually with the addition of the zero-point energy correction in Equation (3). The tunneling correction introduced in TST can be calculated by Skodje–Truhlar [59], Bell35 [57], and Bell58 [58] corrections and by the deformed ( d ) tunneling formulation presented in the next section [60].

2.2.1. Deformed Transition-State Theory ( d -TST)

The deformed Transition-State Theory ( d -TST) [60] formulation:
k d ( T ) = k B T h Q Q 1 Q 2 ( 1 d ε k B T ) 1 / d ,     d = 1 3 ( h ν 2 ε ) 2 ,
where ν is frequency for crossing the barrier, uniformly covers the range from classical to moderate tunneling regimes but needs an amendment for deep tunneling in exothermic reactions, a relatively rare case. The proposed variant of the transition-state theory was obtained from the transitivity concept and deformed ( d ) Aquilanti–Mundim law permit comparison with experiments and tests against alternative formulations. The nomenclature in Equation (4) is the same as used in Equation (3). A popular formulation in the literature is to use the Wigner tunneling correction: However, in our previous study [60,64], we have shown that d -TST is a more satisfactory approximation also in view of uniform behavior across the height of the energy barrier.

2.2.2. Bell35 and Bell58

To cover the moderate-to-deep tunneling transition in exothermic reactions, we applied both Bell35 correction [57,65], Equation (5):
κ B e l l 35 =   [ 1 ν   1 k B T exp ( ε k B T   ε ν ) ] 1 ν   1 k B T ,
and Bell58 [58,65] correction truncated at the second term (2T)
κ B e l l 58 2 T =   ( h ν 2 k B T ) sin ( h ν 2 k B T )   exp ( ε k B T   ε ν ) ( k B T ν   1   ) ,
that, although non-uniform across the transition between negligible and moderate tunneling regimes, were found to behave smoothly enough to adequately perform practically astride the whole range.
Truncating Equation (6) to the first term, we recover the usual formula used in the literature to describe the tunneling under deep regime Equation (7):
κ B e l l 58 = ( h ν 2 k B T ) sin ( h ν 2 k B T )
However, this formulation presents divergence at the crossover temperature T c = h ν / π k B . The tunneling regimes can be delimited within four temperatures ranges—negligible ( T > 2 T c ), small ( 2 T c > T > T c ), moderate ( T c > T > T c / 2 ), and deep ( T > T c / 2 ) [60,66,67].

2.2.3. Skodje and Truhlar, ST

To avoid spurious divergence at T c in the Bell58 formulation, in 1981, Skodje and Truhlar [59] gave a generalization extending the parabolic barrier treatment. In their approximation for tunneling correction Equations (8a) and (8b)
κ S T =   ( h ν 2 k B T ) sin ( h ν 2 k B T )   exp [ ( 1 k B T   1 ν ) ( ε Δ H ) ] ( k B T ν   1   ) ,           β 𝕔 ,
and
κ S T =   1 ( k B T ν   1   ) { exp [ ( 1 k B T   1 ν ) ( ε Δ H ) ] 1 } ,           𝕔 β ,
where 𝕔 = 1 / k B T c = π / h ν and Δ H is the enthalpy of reaction.

2.3. Solvent Effect on Reaction Rate Constant

2.3.1. Collins–Kimball Formulation

Treatment of chemical reactions in liquid-phase requires accounting for the solvent effect, considering the ability of the reagents to diffuse and to lead them to effective reactive collisions. According to the Onsager solvent reaction field model [68,69,70], the solvent creates a solvation layer, as a cage around the molecular entities that participate in the reaction. The reactive process between molecular entities A and B is represented by the sequential equation,
A s o l v +   B s o l v k D k D { A B } s o l v k P r o d u c t .
The “solv” symbol indicates that the molecular entity is surrounded by a roughly spherical solvation layer, the cage, having a specific radius. k D is the diffusion kinetic rate constant which the reagent A travels in the solvent to find B ; vice versa, k D is the kinetic rate constant for reverse diffusion; k denotes the reaction rate constant due to effective reactive collisions: It can be estimated from Equations (3) or (4). Assuming the Steady-State Approximation [70] for Equation (9), the Smoluchowski expression for the diffusion kinetic rate constant [71], takes the form Equation (10):
k D = 4 π r A B D A B ,
where r A B is a reaction radial distance and D A B is the sum of the diffusion constants for each reagent in the solvent. The generalization of Collins and Kimball [61,72] for the irreversible bimolecular diffusion-controlled reactions at infinite reaction rate, k D k D   [73], yields an overall reaction rate constant, k O b s Equation (11)
1 k O b s = 1 k T S T + 1 k D .
In the code, the diffusion is accounted for by the Stokes–Einstein formulation and the temperature dependence of viscosity   η of the solvent is estimated through the Aquilanti–Mundim formula (see details in Section 3).

2.3.2. Kramers’ Formulation

To account for dynamical effects of the solvent in a reactive process and to generalize to unimolecular and pseudo-unimolecular processes, Kramers’ model considers a stochastic motion of the system, where the solvent effect is added considering Brownian movements along the reaction path [70].
Assuming that the friction constant, μ (see below), is independent of time, the overall reaction rate constant k O b s can be calculated as Equation (12)
k O b s = κ K r k T S T ,
where κ K r is the transmission factor obtained by Kramers [62] as:
κ K r = 1 ω ( μ 2 4 + ω 2 μ 2 ) .
For the transition-state theory rate constant k T S T and variants see Section 2.2. In Equation (13), ω is the imaginary frequency of the transition state and the friction constant is given by μ = ( 6 π r A B / M ) η , where r A B and M are the radius of the cage and the molecular mass of the transition-state, respectively. Again, the viscosity is calculated in the code through the Aquilanti–Mundim formula.

3. Handling the Transitivity Code

Several aspects regarding nomenclature and the theoretical background used in this article are presented in already cited references and in Reference [27]. Furthermore, the code with the manual, examples (see next section), and installation video can be freely downloaded in the www.vhcsgroup.com/transitivity web page.
In the main window of the code, the user can choose between three options: (i) Fitting reaction rate constant data as a function of temperature; (ii) predicting the reaction rate constant in the gas- and liquid-phase; and (iii) creating the input for the calculation of first-principles molecular dynamics by the CPMD code. If the option chosen is “Kinetic and Related Parameters”, the program needs the electronic structure output files provided by the Gaussian program for the structures of the reagents, of the transition state and of the products. If the system under study requires the calculation of the electronic energy using a specifically high level method, the energy values must be included separately in the box. The code provides the reaction rate constants for both unimolecular and bimolecular reactions and the users should indicate if the molecular entity corresponds to an atom or whether the molecule is linear or not. In a new window, the user can choose between TST or d -TST [60] and Bell35 [57], Bell58 [58], Skodje–Truhlar (ST) [59] tunneling corrections. In addition, the program provides a visualization of the Arrhenius plot with the possibility of including any experimental and/or theoretical data available for validation. In the same window, reaction properties (such as internal energy, enthalpy, Gibbs-free energy, barrier height, d parameter, imaginary frequency ( ν ), crossover temperature ( T c ), and the parameter of the Skodje–Truhlar model) are exhibited.
If the VOLUME keyword is used in the calculation input file of G09, the option including solvent effects is available through the Kramers and Collins–Kimball formulations. It is necessary to make the choice of solvent parameters for the estimation of the viscosity or use the default (water). The parameters for the viscosity estimation ( η / P o i s e ) of a solvent other than water should be those fitted by the Aquilanti–Mundim formula [25] with η o in Poise and ε in J·mol−1 and inserted in the “Solvent Type” option. The temperature dependence of viscosity of water using the experimental data [74] is expressed as η ( T ) = 2.7024.10 4 P o i s e ( 1 213.0543 / T ) 2.75634 . When the Kramers formulation is chosen, the friction coefficient of the solvent ( μ / s 1 ) , the Kramers transmission coefficient, and the overall reaction rate constant ( k O b s / cm 3 · mol 1 · s 1   or s 1 ) are provided. If the Collins–Kimball formulation is selected, the overall reaction rate constant ( k O b s / cm 3 · mol 1 . s 1 ), separate diffusion coefficients ( cm 2 · s 1 ) for the reactants and the Smoluchowski diffusion rate constant ( k D / cm 3 · mol 1 · s 1 ) are provided.
With the "GSA Fitting" option, the values of the rate constants and temperatures will be needed. In the software, the fitting in the Arrhenius plot is implemented for the rate constant data using the Arrhenius, Aquilanti–Mundim, VFT, NTS, and ASCC formulas by the stochastic optimization algorithm GSA. It is possible to insert the guess parameters d , ε , E a , E 0 , E ν , B and T 0   (see formulas in Table 1). Furthermore, information of the fit are available in “Fit [FormulaName].dat” output file. Internal parameters of GSA can be also controlled [55]. The fitting in the Transitivity plot is implemented of the reaction rate constant data using only Arrhenius, Aquilanti–Mundim, and VFT. The Transitivity plot is calculated using numerical differentiation with the option “Preview” and if necessary the smoothing of the data can be applied with the option “Apply SG” enabling the Savitzky–Golay filter [75]. In addition, the program provides a visualization of the Arrhenius and Transitivity plots with the fitting model chosen step-by-step.
The present code also offers the possibility of creating input to first-principles molecular dynamics simulation by the CPMD computational code. After selecting the input file of the system to be studied, which contains the molecular geometry of the system in a specific extension (*.xyz, *.gjf, *.out and *.log), the user should choose the molecular dynamics method. In addition, the user has the option to choose the Density Functional Theory (DFT) functional, pseudopotential, the temperature, the charge, the simulation time, and the integration time step. The size of the simulation box can be changed in the “Lattices” section. Additionally, the Transitivity code generates another output file with a *.gjf extension, where the user can check if the geometry of the system to be simulated is correct.

4. Examples

4.1. Fitting Mode—Arrhenius and Transitivity Plots

To illustrate the use of most of functionalities and to validate the accuracy of the running of the Transitivity code, a fitting of the reaction rate constants in the Arrhenius plot is performed as a function of temperature for four different systems in different regimes: sub-Arrhenius, corresponding to both deep and moderate tunneling, super-Arrhenius, and anti-Arrhenius behavior.
The first example concerns the keto–enol tautomerization of 2-(2′-hydroxy-4′-methylphenyl) benzoxazole (MeBO), a well-known case of a deep-tunneling regime [7]: For this case we performed a fitting using Arrhenius, Aquilanti–Mundim, Nakamura–Takayanagi–Sato and ASCC formulas. In order to present the results for the sub-Arrhenius behavior under the moderate-tunneling regime, we evaluated the OH + H2 H + H2O reaction [76] using the Arrhenius and Aquilanti–Mundim formulas. The third example is centered on the investigations that have revealed the super-Arrhenius behavior for the rates of the processes promoted by enzymatic catalysis [77,78,79]. Here, we regard the reaction rate constant of the hydride transfer between the substrate and NAD+, catalyzed by F147L, which exhibits a strong convex curvature in the temperature range of 5 to 65 °C [78]: A fitting is performed with the Arrhenius, Aquilanti–Mundim and VTF formulas. The final case is aimed at showing a case of anti-Arrhenius behavior, which is characterized by the decrease of the reaction rate constant with the increase in the temperature. The OH + HBr Br + H2O reaction is prototypical in studies of an example, both from a theoretical and an experimental point of view: It exhibits negative temperature dependence of the reaction rate constant [80,81,82]. Fitting is performed with the Arrhenius and the Aquilanti–Mundim formulas.
The fitting parameters and statistical analysis of the quality of the fits for all the systems are given in Table 1. The statistical measure χ 2   is used in tests to compare the quality of the formula with the reference. It is noteworthy that in a fitting process, in order to obtain parameters with a physical meaning, it is necessary to avoid compensation effects [83] exploiting any prior knowledge of the system, since the formulas with more than two parameters can lead to multiple solutions.
Figure 2 shows the experimental data and fitted formulas for each system in the Arrhenius plot. Both NTS and ASCC formulas were satisfactory for fitting sub-Arrhenius behavior under deep-tunneling regime, while Aquilanti–Mundim is satisfactory just within a specific range of temperature. Under moderate-tunneling regime, the Aquilanti–Mundim formula was an excellent option to describe experimental data for OH + H2 reaction. The VTF and Aquilanti–Mundim formulas seem adequate to fit super-Arrhenius behavior. The Aquilanti–Mundim parameters that were obtained for the OH + HBr reaction indicate that they provide an excellent option for the anti-Arrhenius behavior. As expected, the Arrhenius formula is clearly inadequate to account for deviations at low temperature for all the reactions presented.
Curvature is easily identified through the Arrhenius plot in non-Arrhenius processes; however, the application of the phenomenological formulas to fit this behavior can lead to multiple solutions that make the physical interpretation of the obtained parameters difficult. In Reference [27], we show that curvature in the Arrhenius plot can be linearized using the transitivity plot. Analogous definition to transitivity function were proposed at 1980s to evaluate glass transition of supercooled materials [85,86,87]. However, the definitions failed to propose a connection with the Tolman theorem, Equation (1), being important only as scaling tools for relaxation processes. Figure 3 shows the Transitivity plot for the temperature dependence of relaxation time of the propylene carbonate [86], where a linear behavior is observed for a certain temperature range as expected by the deformed Aquilanti–Mundim law. However, at 198 K a transition to another linear regime is observed—a break not perceived in the Arrhenius plot. With the linearization of the data in the transitivity plot, the fitting process becomes much simpler: In the first high temperature range the Aquilanti–Mundim parameters are 1 ε = 1.71   K 1 and d = 0.32; for the second lower temperature range, 1 ε = 0.53 K 1 and d = 0.08. The interpretation of these results will be published in a future paper considering a large variety of examples.

4.2. Reaction Rate Constants’ Mode

The estimation of the rate constants for the OH + HCl H2O + Cl reaction in the gas-phase was performed to validate TST and tunneling corrections, implemented in the Transitivity code. Furthermore, using the suggestion of the Eyringpy code [54], the NH3 + OH NH2 + H2O reaction was selected to demonstrate the accuracy of the Collins–Kimball and Kramers models to estimate reaction rate constants in an aqueous solution.

4.2.1. The OH + HCl H2O + Cl Reaction

The reaction rate constants for the reaction between hydroxyl radical and hydrogen chloride are only slightly dependent on temperature in the range 138–300 K, although, as the range increases further, a temperature dependence is observed. This strong concave curvature detected in the Arrhenius plot (sub-Arrhenius behavior) is very convincing evidence of the role of deep quantum tunneling for this reaction [88]. We employed the d -TST, Bell35, Bell58, and ST tunneling corrections to calculate the rate constants for the reaction of OH + HCl over a wide range of temperatures (200–2000 K). The electronic structure properties of the reactants, of the products, and of the transition state were calculated employing the MP2/aug-cc-pVDZ calculation level using Gaussian 09 [89]. A complete study of this reaction using our methodology can be found in Reference [44].
Figure 4 obtained by the Transitivity code shows the comparison between the calculated reaction rate constant and the experimental data [90]. No major differences were found using either the Bel58-2T and ST tunneling corrections. Conversely, as expected, a divergence is observed of the Bell58 formula at T c , and d -TST does not describe the range of experimental data at low temperature, where the deep-tunneling regime becomes dominant: This confirms that its validity is limited to weak tunneling. Traditional TST is presented for comparative purposes and visualizes the need for corrections.

4.2.2. The NH3 + OH NH2 + H2O Reaction

The NH3 + OH NH2 + H2O reaction permits to illustrate the accuracy of the methodology in the liquid-phase. Energies, geometries, and frequencies of stationary points were extracted at the same level of calculation used in the Eyringpy code [54].
The upper panels in Figure 5 show the temperature dependence of the overall reaction rate constant, k O b s , for NH3 + OH NH2 + H2O reaction from 273 to 4000 K using Kramers and Collins–Kimball models. The Smoluchowski diffusion rate constant k D , which evaluates the diffusion limit for a bimolecular reactive process including the solvent effect, is shown in the lower right panel. The Kramers transmission correction, which evaluates the interference of the friction effect of the solvent in the reactive process as a function of temperature, is shown in the lower left panel. At 298.15 K, Kramers’ formulation gives for the reaction rate constant the value 6.73 × 1011 cm3 mol−1 s−1, while the Collins-Kimball formulation yields 6.77 × 1011 cm3 mol−1 s−1, (the experimental value indicates ~1011 cm3 mol−1 s−1) [91,92,93]. The value of   k D from Smoluchowski (Collins–Kimball) is 3.73 × 1012 cm3 mol−1 s−1, in accordance with that calculated in Reference [54], k D = 3.60 × 1012 cm3 mol−1 s−1.

4.3. CPMD Input Files Generator

Coordinates of specific molecules can be selected to test the function of creating a first-principles molecular dynamics input to run by the CPMD code (Figure 6). When selecting this option, the user must choose between the CPMD, PIMD, TSH, MTD and BOMD approaches and must provide the molecular geometry of the system in *.xyz format. Additional information (DFT functional, simulation temperature, system charge, the maximum number of steps, integration time step, and the size of the box) can be inserted through the indicated boxes. Furthermore, it is possible to generate several aleatory configurations by selecting specific initial conditions for the geometric parameters (For more information, see the input files in the Example directory on the www.vhcsgroup.com/transitivity web page).

5. Final Remarks

The Transitivity code presented in this article tends to systematize the tools developed over the last ten years [25] to handle the kinetics of processes beyond Arrhenius, providing options for: (i) The phenomenological fitting and estimation of reaction rate constants in gas- and liquid-phases and (ii) the preparation of first-principles molecular dynamics in order to evaluate the parameters relative to a variety of reactive processes. In addition, the code provides an easy user-friendly interface and may be relevant for didactic purposes.
A characterizing feature of the code is the consistent use of the d -formulation, which recently culminated in a series of successful applications from phenomenological to first-principles descriptions of pure and applied chemical kinetics and material science. Examples are available:
  • From the phenomenology of elementary processes (such as the H2 + F [94], OH + HBr [80], F + HD [65] and C + CH+ [95] reactions) to complex processes (such as food systems [96], plant respiration [25], plasma chemistry [97], and solid-state diffusive reaction [98]);
  • Calculation of the kinetic rate constants for chemical reactions from the potential energy surface features profile, such as the CH4 + OH [60], CH3OH + H [99], OH + HCl [44], OH + HI [43], to proton rearrangement of enol forms of curcumin [100], OH + H2 [101], and chiral nucleophilic substitution reaction [102].
Regarding the evolution and consistency of the code, efforts are being made to introduce d predictive formulations to transport properties, as well as post-processing of trajectories obtained by first-principles molecular dynamics simulations. Furthermore, the Eckart tunneling correction and the variational Transition-State Theory represent focuses for the future implementation of the Transitivity code.

Author Contributions

V.H.C.-S. and F.P. conceived and supervised the study. V.H.C.-S. and N.D.C. wrote the paper. H.G.M. and F. O.S.-N. developed the python code, performed stationary electronic structure calculations and fitting procedures. K.C.M. developed and implemented the GSA Fortran code. The manuscript was discussed and reviewed by all authors.

Funding

This research was funded by Brazilian agency CNPq, grant number Universal 01/2016 - Faixa A - 406063/2016-8], by Organizzazione Internazionale Italo-latino Americana (IILA) for a Biotechnology Sector-2019 scholarship and by the Italian Ministry for Education, University and Research, MIUR, for financial support: SIR 2014 ‘‘Scientific Independence for Young Researchers’’ (RBSI14U3VF).

Acknowledgments

The authors are grateful for the support given by Brazilian agency CAPES. This research is also supported by the High-Performance Computing Center at the Universidade Estadual de Goiás, Brazil. Valter H. Carvalho-Silva thanks Brazilian agency CNPq for the research funding programs [Universal 01/2016 - Faixa A - 406063/2016-8] and Organizzazione Internazionale Italo-latino Americana (IILA) for a Biotechnology Sector-2019 scholarship. Federico Palazzetti and Nayara D. Coutinho acknowledge the Italian Ministry for Education, University and Research, MIUR, for financial support: SIR 2014 ‘‘Scientific Independence for Young Researchers’’ (RBSI14U3VF). We thank professor Vincenzo Aquilanti for fruitful discussions.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Table A1. List of Symbols and Nomenclatures.
Table A1. List of Symbols and Nomenclatures.
SymbolsNomenclature
k Rate constant
TTemperature
k B Boltzmann constant
β Lagrange multiplier
γ Transitivity function
d Deformed parameter
h Planck’s constant
Q Partition functions
AMAquilanti-Mundim
d H Enthalpy of reaction
ASCCAquilanti–Sanchez–Coutinho–Carvalho
NTSNakamura–Takayanagi–Sato
TSTTransition-State Theory
GSAGeneralized Simulated Annealing
d -TSTDeformed Transition-State Theory
STSkodje and Truhlar tunneling correction
Bell35Bell’s tunneling correction of 1935
Bell58Bell’s tunneling correction of 1958
ε barrier height (Eyring’s parameter)
E a Apparent Activation Energy
E 0 Energy parameter from NTS formula
E υ   Energy parameter from ASCC formula
BTemperature parameter from VFT formula
T 0 Temperature parameter from NTS and VFT formulas.
T c Crossover temperature
k D Diffusion rate constant
ν Imaginary frequency
k O b s Overall reaction rate constant
κ K r Transmission factor from Kramers’ model
μ Friction constant
η Viscosity
DFTDensity functional theory
BOMDBorn-Oppenheimer molecular dynamics
CPMDCar-Parrinello molecular dynamics
PIMDPath-Integral molecular dynamics
MTDMetadynamics
TSHTrajectory Surface Hopping

References

  1. Valter Henrique, C.-S.; Coutinho, N.D.; Aquilanti, V. From the Kinetic Theory of Gases to the Kinetics of Chemical Reactions: On the Verge of the Thermodynamical and the Kinetic Limits. Molecules 2019. to be submitted. [Google Scholar]
  2. Aquilanti, V.; Coutinho, N.D.; Carvalho-Silva, V.H. Kinetics of Low-Temperature Transitions and Reaction Rate Theory from Non-Equilibrium Distributions. Philos. Trans. R. Soc. London A 2017, 375, 20160204. [Google Scholar] [CrossRef]
  3. Aquilanti, V.; Borges, E.P.; Coutinho, N.D.; Mundim, K.C.; Carvalho-Silva, V.H. From statistical thermodynamics to molecular kinetics: the change, the chance and the choice. Rend. Lincei. Sci. Fis. e Nat. 2018, 28, 787–802. [Google Scholar] [CrossRef]
  4. Gentili, P.L. The fuzziness of the molecular world and its perspectives. Molecules 2018, 23, 2074. [Google Scholar] [CrossRef]
  5. Gentili, P.L. Untangling Complex Systems: A Grand Challenge for Science; CRC Press: Boca Raton, FL, USA, 2018; ISBN 9780429847547. [Google Scholar]
  6. Atkinson, R. Kinetics and mechanisms of the gas-phase reactions of the hydroxyl radical with organic compounds under atmospheric conditions. Chem. Rev. 1986, 86, 69–201. [Google Scholar] [CrossRef]
  7. Limbach, H.-H.; Miguel Lopez, J.; Kohen, A. Arrhenius curves of hydrogen transfers: tunnel effects, isotope effects and effects of pre-equilibria. Philos. Trans. R. Soc. B Biol. Sci. 2006, 361, 1399–1415. [Google Scholar] [CrossRef] [Green Version]
  8. Smith, I.W.M. The temperature-dependence of elementary reaction rates: Beyond Arrhenius. Chem. Soc. Rev. 2008, 37, 812–826. [Google Scholar] [CrossRef]
  9. Sims, I.R. Low-temperature reactions: Tunnelling in space. Nat. Chem. 2013, 5, 734–736. [Google Scholar] [CrossRef] [PubMed]
  10. Peleg, M.; Normand, M.D.; Corradini, M.G. The Arrhenius Equation Revisited. Crit. Rev. Food Sci. Nutr. 2012, 52, 830–851. [Google Scholar] [CrossRef] [PubMed]
  11. Darrington, R.T.; Jiao, J. Rapid and Accurate Prediction of Degradant Formation Rates in Pharmaceutical Formulations Using High-Performance Liquid Chromatography-Mass Spectrometry. J. Pharm. Sci. 2004, 93, 838–846. [Google Scholar] [CrossRef]
  12. Giordano, D.; Russell, J.K. Towards a structural model for the viscosity of geological melts. Earth Planet. Sci. Lett. 2018, 501, 202–212. [Google Scholar] [CrossRef]
  13. Klinman, J.P.; Kohen, A. Hydrogen Tunneling Links Protein Dynamics to Enzyme Catalysis. Annu. Rev. Biochem. 2013, 82, 471–496. [Google Scholar] [CrossRef] [Green Version]
  14. Warshel, A.; Bora, R.P. Perspective: Defining and quantifying the role of dynamics in enzyme catalysis. J. Chem. Phys. 2016, 144, 180901. [Google Scholar] [CrossRef]
  15. Laidler, K.J. A Glossary of Terms Used in Chemical Kinetics, Including Reaction Dynamics. Pure Appl. Chem. 1996, 68, 149–192. [Google Scholar] [CrossRef]
  16. Tolman, R.C. Statistical Mechanics Applied to Chemical Kinetics. J. Amer. Chem. Soc. 1920, 42, 2506–2528. [Google Scholar] [CrossRef]
  17. Glasstone, S.; Laidler, K.J.; Eyring, H. The Theory of Rate Processes: The Kinetics of Chemical Reactions, Viscosity, Diffusion and Electrochemical Phenomena; International chemical series; McGraw-Hill: New York, NY, USA, 1941. [Google Scholar]
  18. Truhlar, D.G. Current Status of Transition-State Theory. J. Phys. Chem. 1983, 2664–2682. [Google Scholar] [CrossRef]
  19. Kooij, D.M. Über die Zersetzung des gasförmigen Phosphorwasserstoffs. Zeitschrift für Phys. Chemie 1893, 12, 155–161. [Google Scholar]
  20. Bělehrádek, J. A unified theory of cellular rate processes based upon an analysis of temperature action. Protoplasma 1957, 48, 53–71. [Google Scholar] [CrossRef]
  21. Vogel, H. Das temperature-abhangigketsgesetz der viskositat von flussigkeiten. Phys. Z 1921, 22, 645–646. [Google Scholar]
  22. Fulcher, G.S. Analysis of Recent Measurements of the Viscosity of Glasses. J. Am. Ceram. Soc. 1925, 8, 339–355. [Google Scholar] [CrossRef]
  23. Tammann, G.; Hesse, W. Die Abhängigkeit der Viscosität von der Temperatur bie unterkühlten Flüssigkeiten. Zeitschrift für Anorg. und Allg. Chemie 1926, 156, 245–257. [Google Scholar] [CrossRef]
  24. Nakamura, K.; Takayanagi, T.; Sato, S. A modified arrhenius equation. Chem. Phys. Lett. 1989, 160, 295–298. [Google Scholar] [CrossRef]
  25. Aquilanti, V.; Mundim, K.C.; Elango, M.; Kleijn, S.; Kasai, T. Temperature dependence of chemical and biophysical rate processes: Phenomenological approach to deviations from Arrhenius law. Chem. Phys. Lett. 2010, 498, 209–213. [Google Scholar] [CrossRef]
  26. Coutinho, N.D.; Silva, Y.S.; de Fazio, D.; Cavalli, S.; Carvalho-Silva, V.H.; Aquilanti, V. Chemical Kinetics under Extreme Conditions: Exact, Phenomenological and First-Principles Computational Approaches. In The Astrochemical Observatory: Focus on Chiral Molecules; Accademia Nazionale delle Scienze detta dei XL: Rome, Italy, 2018; pp. 1–15. [Google Scholar]
  27. Carvalho-Silva, V.H.; Coutinho, N.D.; Aquilanti, V. Temperature dependence of rate processes beyond Arrhenius and Eyring: Activation and Transitivity. Front. Chem. 2019, 7, 380. [Google Scholar] [CrossRef] [PubMed]
  28. Fernandez-Ramos, A.; Ellingson, B.A.; Garrett, B.C.; Truhlar, D.G. Variational Transition State Theory with Multidimensional Tunneling. In Reviews in Computational Chemistry; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2007; pp. 125–232. ISBN 9780470116449. [Google Scholar]
  29. Marcus, R.A. Electron transfer reactions in chemistry. Theory and experiment. Rev. Mod. Phys. 1993, 65, 599–610. [Google Scholar] [CrossRef] [Green Version]
  30. Miller, W. Quantum mechanical transition state theory and a new semiclassical model for reaction rate constants. J. Chem. Phys. 1974, 61, 1823. [Google Scholar] [CrossRef]
  31. Richardson, J.O. Ring-Polymer Approaches to Instanton Theory. Ph.D. Thesis, University of Cambridge, Cambridge, UK, 2012. [Google Scholar]
  32. Zhang, Y.; Stecher, T.; Cvitaš, M.T.; Althorpe, S.C. Which Is Better at Predicting Quantum-Tunneling Rates: Quantum Transition-State Theory or Free-Energy Instanton Theory? J. Phys. Chem. Lett. 2014, 5, 3976–3980. [Google Scholar] [CrossRef]
  33. Xiao, R.; Gao, L.; Wei, Z.; Spinney, R.; Luo, S.; Wang, D.; Dionysiou, D.D.; Tang, C.J.; Yang, W. Mechanistic insight into degradation of endocrine disrupting chemical by hydroxyl radical: An experimental and theoretical approach. Environ. Pollut. 2017, 231, 1446–1452. [Google Scholar] [CrossRef]
  34. Luo, S.; Gao, L.; Wei, Z.; Spinney, R.; Dionysiou, D.D.; Hu, W.P.; Chai, L.; Xiao, R. Kinetic and mechanistic aspects of hydroxyl radical‒mediated degradation of naproxen and reaction intermediates. Water Res. 2018, 137, 233–241. [Google Scholar] [CrossRef]
  35. Gao, Y.; Ji, Y.; Li, G.; An, T. Mechanism, kinetics and toxicity assessment of OH-initiated transformation of triclosan in aquatic environments. Water Res. 2014, 49, 360–370. [Google Scholar] [CrossRef]
  36. De Sainte Claire, P. Degradation of PEO in the solid state: A theoretical kinetic model. Macromolecules 2009, 42, 3469–3482. [Google Scholar] [CrossRef]
  37. Ahubelem, N.; Shah, K.; Moghtaderi, B.; Page, A.J. Formation of benzofuran and chlorobenzofuran from 1,3-dichloropropene: A quantum chemical investigation. Int. J. Quantum Chem. 2015, 115, 1739–1745. [Google Scholar] [CrossRef]
  38. Zavala-Oseguera, C.; Galano, A.; Merino, G. Computational study on the kinetics and mechanism of the carbaryl + OH reaction. J. Phys. Chem. A 2014, 118, 7776–7781. [Google Scholar] [CrossRef] [PubMed]
  39. Döntgen, M.; Przybylski-Freund, M.-D.; Kröger, L.C.; Kopp, W.A.; Ismail, A.E.; Leonhard, K. Automated discovery of reaction pathways, rate constants, and transition states using reactive molecular dynamics simulations. J. Chem. Theory Comput. 2015, 11, 2517–2524. [Google Scholar] [CrossRef] [PubMed]
  40. Piccini, G.; McCarty, J.; Valsson, O.; Parrinell, M. Variational Flooding Study of a S2 Reaction. J. Phys. Chem. A 2017, 8, 580–583. [Google Scholar]
  41. Fleming, K.L.; Tiwary, P.; Pfaendtner, J. New Approach for Investigating Reaction Dynamics and Rates with Ab Initio Calculations. J. Phys. Chem. A 2016, 120, 299–305. [Google Scholar] [CrossRef] [PubMed]
  42. Lancar, I.T.; Mellouki, A.; Poulet, G. Kinetics of the reactions of hydrogen iodide with hydroxyl and nitrate radicals. Chem. Phys. Lett. 1991, 177, 554–558. [Google Scholar] [CrossRef]
  43. Coutinho, N.D.; Carvalho-Silva, V.H.; de Oliveira, H.C.B.; Aquilanti, V. The HI + OH → H2O + I Reaction by First-Principles Molecular Dynamics: Stereodirectional and Anti-Arrhenius Kinetics, 2017; Volume 10408 LNCS, ISBN 9783319624037.
  44. Coutinho, N.D.; Sanches-Neto, F.O.; Carvalho-Silva, V.H.; de Oliveira, H.C.B.; Ribeiro, L.A.; Aquilanti, V. Kinetics of the OH+HCl→H2O+Cl reaction: Rate determining roles of stereodynamics and roaming and of quantum tunneling. J. Comput. Chem. 2018, 39, 2508–2516. [Google Scholar] [CrossRef] [PubMed]
  45. Isaacson, A.D.; Truhlar, D.G.; Rai, S.N.; Steckler, R.; Hancock, G.C.; Garrett, B.C.; Redmon, M.J. POLYRATE: A general computer program for variational transition state theory and semiclassical tunneling calculations of chemical reaction rates. Comput. Phys. Commun. 1987, 47, 91–102. [Google Scholar] [CrossRef]
  46. Duncan, W.T.; Bell, R.L.; Truong, T.N. TheRate: Program forab initio direct dynamics calculations of thermal and vibrational-state-selected rate constants. J. Comput. Chem. 1998, 19, 1039–1052. [Google Scholar] [CrossRef]
  47. Barker, J.R. Multiple-Well, multiple-path unimolecular reaction systems. I. MultiWell computer program suite. Int. J. Chem. Kinet. 2001, 33, 232–245. [Google Scholar] [CrossRef] [Green Version]
  48. Ghysels, A.; Verstraelen, T.; Hemelsoet, K.; Waroquier, M.; Van Speybroeck, V. TAMkin: A versatile package for vibrational analysis and chemical kinetics. J. Chem. Inf. Model. 2010, 50, 1736–1750. [Google Scholar] [CrossRef] [PubMed]
  49. Glowacki, D.R.; Liang, C.-H.; Morley, C.; Pilling, M.J.; Robertson, S.H. MESMER: An Open-Source Master Equation Solver for Multi-Energy Well Reactions. J. Phys. Chem. A 2012, 116, 9545–9560. [Google Scholar] [CrossRef] [PubMed]
  50. Gao, C.W.; Allen, J.W.; Green, W.H.; West, R.H. Reaction Mechanism Generator: Automatic construction of chemical kinetic mechanisms. Comput. Phys. Commun. 2016, 203, 212–225. [Google Scholar] [CrossRef] [Green Version]
  51. Euclides, H.O.; Barreto, P.R. APUAMA: a software tool for reaction rate calculations. J. Mol. Model. 2017, 23, 176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Canneaux, S.; Bohr, F.; Henon, E. KiSThelP: A program to predict thermodynamic properties and rate constants from quantum chemistry results. J. Comput. Chem. 2014, 35, 82–93. [Google Scholar] [CrossRef] [PubMed]
  53. Coppola, C.M. Mher V Kazandjian; Matrix formulation of the energy exchange problem of multi-level systems and the code FRIGUS. Rend. Lincei Sci. Fis. e Nat. 2019, in press. [Google Scholar]
  54. Dzib, E.; Cabellos, J.L.; Ortíz-Chi, F.; Pan, S.; Galano, A.; Merino, G. Eyringpy: A program for computing rate constants in the gas phase and in solution. Int. J. Quantum Chem. 2018, 119, 11–13. [Google Scholar] [CrossRef]
  55. Mundim, K.C.; Tsallis, C. Geometry optimization and conformational analysis through generalized simulated annealing. Int. J. Quantum Chem. 1998, 58, 373–381. [Google Scholar] [CrossRef]
  56. Sato, S. Tunneling in bimolecular reactions. Chem. Phys. 2005, 315, 65–75. [Google Scholar] [CrossRef]
  57. Bell, R.P. Quantum Mechanical Effects in Reactions Involving Hydrogen. Proc. R. Soc. London. Ser. A, Math. Phys. Sci. 1935, CXLVIII.A, 241–250. [Google Scholar]
  58. Bell, R.P. The Tunnel Effect Correction For Parabolic Potential Barriers. Faraday Soc. Contrib. 1958, 1–4. [Google Scholar] [CrossRef]
  59. Skodje, R.T.; Truhlar, D.G. Parabolic tunneling calculations. J. Phys. Chem. 1981, 85, 624–628. [Google Scholar] [CrossRef]
  60. Carvalho-Silva, V.H.; Aquilanti, V.; de Oliveira, H.C.B.; Mundim, K.C. Deformed transition-state theory: Deviation from Arrhenius behavior and application to bimolecular hydrogen transfer reaction rates in the tunneling regime. J. Comput. Chem. 2017, 38, 178–188. [Google Scholar] [CrossRef] [PubMed]
  61. Collins, F.C.; Kimball, G.E. Diffusion-controlled reaction rates. J. Colloid Sci. 1949, 4, 425–437. [Google Scholar] [CrossRef]
  62. Kramers, H.A. Brownian motion in a field of force and the diffusion model of chemical reactions. Phys. 1940, 7, 284–304. [Google Scholar] [CrossRef]
  63. CPMDversion 3.17.1; CPMD, version 4.; CPMDversion 4.1; CPMDversion 3.17.1 Copyright IBM 2012.
  64. Claudino, D.; Gargano, R.; Carvalho-Silva, V.H.; E Silva, G.M.; Da Cunha, W.F. Investigation of the Abstraction and Dissociation Mechanism in the Nitrogen Trifluoride Channels: Combined Post-Hartree-Fock and Transition State Theory Approaches. J. Phys. Chem. A 2016, 120, 5464–5473. [Google Scholar] [CrossRef]
  65. Cavalli, S.; Aquilanti, V.; Mundim, K.C.; De Fazio, D. Theoretical reaction kinetics astride the transition between moderate and deep tunneling regimes: The F + HD case. J. Phys. Chem. A 2014, 118, 6632–6641. [Google Scholar] [CrossRef]
  66. Bell, R.P. The Tunnel Effect in Chemistry; Champman and Hall: London, UK, 1980. [Google Scholar]
  67. Christov, S.G. The Characteristic (Crossover) Temperature in the Theory of Thermally Activated Tunneling Processes. Mol. Eng. 1997, 7, 109–147. [Google Scholar] [CrossRef]
  68. Onsager, L. Electric Moments of Molecules in Liquids. J. Am. Chem. Soc. 1936, 58, 1486–1493. [Google Scholar] [CrossRef]
  69. Wong, M.W.; Wiberg, K.B.; Frisch, M.J. Solvent effects. 3. Tautomeric equilibria of formamide and 2-pyridone in the gas phase and solution: an ab initio SCRF study. J. Am. Chem. Soc. 1992, 114, 1645–1652. [Google Scholar] [CrossRef]
  70. Henriksen, N.E.; Hansen, F.Y. Theories of Molecular Reaction Dynamics: The Microscopic Foundation of Chemical Kinetics; Oxford University Press: New York, NY, USA, 2008; ISBN 9780191708251. [Google Scholar]
  71. Smoluchowski, M.V. Drei Vortrage uber Diffusion, Brownsche Bewegung und Koagulation von Kolloidteilchen. Phys. Zeit. 1916, 17, 557–585. [Google Scholar]
  72. Collins, F.C.; Kimball, G.E. Diffusion-Controlled Reactions in Liquid Solutions. Ind. Eng. Chem. 1949, 41, 2551–2553. [Google Scholar] [CrossRef]
  73. Eigen, M. Acid-Base Catalysis, and Enzymatic Hydrolysis. Part I: Elementary Processes. Angew. Chemie Int. Ed. English 1964, 3, 1–19. [Google Scholar] [CrossRef]
  74. Hallett, J. The Temperature Dependence of the Viscosity of Supercooled Water. Proc. Phys. Soc. 1963, 82, 1046–1050. [Google Scholar] [CrossRef]
  75. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  76. Ravishankara, A.R.; Nicovich, J.M.; Thompson, R.L.; Tully, F.P. Kinetic study of the reaction of hydroxyl with hydrogen and deuterium from 250 to 1050 K. J. Phys. Chem. 1981, 85, 2498–2503. [Google Scholar] [CrossRef]
  77. Kohen, A.; Cannio, R.; Bartolucci, S.; Klinman, J.P. Enzyme dynamics and hydrogen tunnelling in a thermophilic alcohol dehydrogenase. Nature 1999, 399, 496–499. [Google Scholar] [CrossRef]
  78. Liang, Z.X.; Tsigos, I.; Bouriotis, V.; Klinman, J.P. Impact of protein flexibility on hydride-transfer parameters in thermophilic and psychrophilic alcohol dehydrogenases. J. Am. Chem. Soc. 2004, 126, 9500–9501. [Google Scholar] [CrossRef]
  79. Truhlar, D.; Kohen, A. Convex Arrhenius plots and their interpretation. Proc. Nat. Acad. Sci. USA 2001, 98, 848–851. [Google Scholar] [CrossRef] [Green Version]
  80. Coutinho, N.D.; Aquilanti, V.; Silva, V.H.C.; Camargo, A.J.; Mundim, K.C.; De Oliveira, H.C.B. Stereodirectional Origin of anti-Arrhenius Kinetics for a Tetraatomic Hydrogen Exchange Reaction: Born-Oppenheimer Molecular Dynamics for OH + HBr. J. Phys. Chem. A 2016, 120, 5408–5417. [Google Scholar] [CrossRef] [PubMed]
  81. de Oliveira-Filho, A.G.S.; Ornellas, F.R.; Bowman, J.M. Quasiclassical Trajectory Calculations of the Rate Constant of the OH + HBr → Br + H2O Reaction Using a Full-Dimensional Ab Initio Potential Energy Surface Over the Temperature Range 5 to 500 K. J. Phys. Chem. Lett. 2014, 5, 706–712. [Google Scholar] [CrossRef] [PubMed]
  82. Coutinho, N.D.; Silva, V.H.C.; de Oliveira, H.C.B.; Camargo, A.J.; Mundim, K.C.; Aquilanti, V. Stereodynamical Origin of Anti-Arrhenius Kinetics: Negative Activation Energy and Roaming for a Four-Atom Reaction. J. Phys. Chem. Lett. 2015, 6, 1553–1558. [Google Scholar] [CrossRef] [PubMed]
  83. Zuniga-Hansen, N.; Silbert, L.E.; Calbi, M.M. Breakdown of kinetic compensation effect in physical desorption. Phys. Rev. E 2018, 98, 032128. [Google Scholar] [CrossRef] [Green Version]
  84. Sims, I.R.; Smith, I.W.M.; Clary, D.C.; Bocherel, P.; Rowe, B.R. Ultra-low Temperature Kinetics of Neutral-neutral Reactions - New Experimental and Theoretical Results For OH + HBr Between 295 K and 23 K. J. Chem. Phys. 1994, 101, 1748–1751. [Google Scholar] [CrossRef]
  85. Souletie, J.; Tholence, J.L. Critical slowing down in spin glasses and other glasses: Fulcher versus power law. Phys. Rev. B 1985, 32, 516. [Google Scholar] [CrossRef]
  86. Stickel, F.; Fischer, E.W.; Richert, R. Dynamics of glass-forming liquids. II. Detailed comparison of dielectric relaxation, de-conductivity, and viscosity data. J. Chem. Phys. 1996, 104, 2043. [Google Scholar] [CrossRef]
  87. Drozd-Rzoska, A. Universal behavior of the apparent fragility in ultraslow glass forming systems. Sci. Rep. 2019, 9, 6816. [Google Scholar] [CrossRef]
  88. Silva, V.H.C.; Aquilanti, V.; de Oliveira, H.C.B.; Mundim, K.C. Uniform description of non-Arrhenius temperature dependence of reaction rates, and a heuristic criterion for quantum tunneling vs classical non-extensive distribution. Chem. Phys. Lett. 2013, 590, 201–207. [Google Scholar] [CrossRef] [Green Version]
  89. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian09 Revision D.01 2009; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  90. Ravishankara, A.R.; Wine, P.H.; Wells, J.R.; Thompson, R.L. Kinetic study of the reaction of OH with HCl from 240 to 1055 K. Chem. Phys. Lett. 1985, 17, 1281–1297. [Google Scholar] [CrossRef]
  91. Hickel, B.; Sehested, K. Reaction of hydroxyl radicals with ammonia in liquid water at elevated temperatures. Int. J. Radiat. Appl. Instrum. Part C. Radiat. Phys. Chem. 1992, 39, 355–357. [Google Scholar] [CrossRef]
  92. Men’kin, V.B.; Makarov, I.E.; Pikaev, A.K. Pulse radiolysis study of reaction rates of OH and O-radicals with ammonia in aqueous solutions. High Energy Chem. (Engl. Transl.) 1989, 22, 333–336. [Google Scholar]
  93. Neta, P.; Maruthamuthu, P.; Carton, P.M.; Fessenden, R.W. Formation and reactivity of the amino radical. J. Phys. Chem. 1978, 82, 1875–1878. [Google Scholar] [CrossRef]
  94. Aquilanti, V.; Mundim, K.C.; Cavalli, S.; De Fazio, D.; Aguilar, A.; Lucas, J.M. Exact activation energies and phenomenological description of quantum tunneling for model potential energy surfaces. The F+H2 reaction at low temperature. Chem. Phys. 2012, 398, 186–191. [Google Scholar] [CrossRef]
  95. Rampino, S.; Pastore, M.; Garcia, E.; Pacifici, L.; Laganà, A. On the temperature dependence of the rate coefficient of formation of C2+ from C + CH+. Mon. Not. R. Astron. Soc. 2016, 460, 2368–2375. [Google Scholar] [CrossRef]
  96. Coutinho, N.D.; Silva, V.H.C.; Mundim, K.C.; de Oliveira, H.C.B. Description of the effect of temperature on food systems using the deformed Arrhenius rate law: deviations from linearity in logarithmic plots vs. inverse temperature. Rend. Lincei 2015, 26, 141–149. [Google Scholar] [CrossRef]
  97. Capitelli, M.; Pietanza, L.D. Past and present aspects of Italian plasma chemistry. Rend. Lincei. Sci. Fis. e Nat. 2019. [Google Scholar] [CrossRef]
  98. Agreda, N.J.L. Aquilanti–Mundim deformed Arrhenius model in solid-state reactions: Theoretical evaluation using DSC experimental data. J. Therm. Anal. Calorim. 2016, 126, 1175–1184. [Google Scholar] [CrossRef]
  99. Sanches-Neto, F.O.; Coutinho, N.D.; Silva, V. A novel assessment of the role of the methyl radical and water formation channel in the CH3OH + H reaction. Phys. Chem. Chem. Phys. 2017, 19, 24467–24477. [Google Scholar] [CrossRef] [PubMed]
  100. Santin, L.G.; Toledo, E.M.; Carvalho-Silva, V.H.; Camargo, A.J.; Gargano, R.; Oliveira, S.S. Methanol Solvation Effect on the Proton Rearrangement of Curcumin’s Enol Forms: An Ab Initio Molecular Dynamics and Electronic Structure Viewpoint. J. Phys. Chem. C 2016, 120, 19923–19931. [Google Scholar] [CrossRef]
  101. Carvalho-Silva, V.H.; Vaz, E.C.; Coutinho, N.D.; Kobayashi, H.; Kobayashi, Y.; Kasai, T.; Palazzetti, F.; Lombardi, A.; Aquilanti, V. The Increase of the Reactivity of Molecular Hydrogen with Hydroxyl Radical from the Gas Phase versus an Aqueous Environment: Quantum Chemistry and Transition State-Theory Calculations. In Lecture Notes in Computer Science; Springer: Cham, Switzerland, 2019; pp. 450–459. [Google Scholar]
  102. Rezende, M.V.C.S.; Coutinho, N.D.; Palazzetti, F.; Lombardi, A.; Carvalho-Silva, V.H. Nucleophilic substitution vs elimination reaction of bisulfide ions with substituted methanes: exploration of chiral selectivity by stereodirectional first-principles dynamics and transition state theory. J. Mol. Model. 2019, 25, 227. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Not available.
Figure 1. Logo and main windows of the Transitivity code.
Figure 1. Logo and main windows of the Transitivity code.
Molecules 24 03478 g001
Figure 2. Arrhenius plots comparing the experimental reaction rate constant and fitted formulas for keto–enol tautomerization reaction (sub-Arrhenius behavior under deep tunneling), OH + H2 ⟶ H2O + H reaction (sub-Arrhenius behavior under moderate tunneling), hydride transfer with enzymatic catalysis (super-Arrhenius behavior) and OH + HBr ⟶ H2O + Br reaction (anti-Arrhenius behavior). NTS and ASCC formulas were of use for sub-Arrhenius behavior under deep-tunneling regime. The Aquilanti–Mundim formula was of use for sub-Arrhenius cases under moderate-tunneling regime, for super-Arrhenius and for anti-Arrhenius behaviors. VFT also was of use for super-Arrhenius situations. The references of experimental data can be found in Table 1.
Figure 2. Arrhenius plots comparing the experimental reaction rate constant and fitted formulas for keto–enol tautomerization reaction (sub-Arrhenius behavior under deep tunneling), OH + H2 ⟶ H2O + H reaction (sub-Arrhenius behavior under moderate tunneling), hydride transfer with enzymatic catalysis (super-Arrhenius behavior) and OH + HBr ⟶ H2O + Br reaction (anti-Arrhenius behavior). NTS and ASCC formulas were of use for sub-Arrhenius behavior under deep-tunneling regime. The Aquilanti–Mundim formula was of use for sub-Arrhenius cases under moderate-tunneling regime, for super-Arrhenius and for anti-Arrhenius behaviors. VFT also was of use for super-Arrhenius situations. The references of experimental data can be found in Table 1.
Molecules 24 03478 g002
Figure 3. The Arrhenius (upper panel) and Transitivity (lower panel) planes of the temperature dependence of relaxation time of the propylene carbonate. The diamond symbols represent the transitivity values obtained numerically and smoothing with the Savitzky–Golay filter. Red lines emphasize two regions where the temperature dependence of the transitivity is linearized, as expected by the Aquilanti–Mundim law.
Figure 3. The Arrhenius (upper panel) and Transitivity (lower panel) planes of the temperature dependence of relaxation time of the propylene carbonate. The diamond symbols represent the transitivity values obtained numerically and smoothing with the Savitzky–Golay filter. Red lines emphasize two regions where the temperature dependence of the transitivity is linearized, as expected by the Aquilanti–Mundim law.
Molecules 24 03478 g003
Figure 4. Arrhenius plot obtained from the Transitivity code for the OH + HCl Cl + H2O reaction using TST with Bell35, Bell58, ST tunneling correction, and d -TST. Experimental data in the literature [90] are available for comparison and shown as full dots.
Figure 4. Arrhenius plot obtained from the Transitivity code for the OH + HCl Cl + H2O reaction using TST with Bell35, Bell58, ST tunneling correction, and d -TST. Experimental data in the literature [90] are available for comparison and shown as full dots.
Molecules 24 03478 g004
Figure 5. Upper panels present the Arrhenius plots as given by the program for the NH3 + OH NH2 + H2O reaction using Kramers’ and Collins–Kimball formulations. The lower panels show the Kramers transmission and Smoluchowski diffusion limit constant as a function of inverse temperature.
Figure 5. Upper panels present the Arrhenius plots as given by the program for the NH3 + OH NH2 + H2O reaction using Kramers’ and Collins–Kimball formulations. The lower panels show the Kramers transmission and Smoluchowski diffusion limit constant as a function of inverse temperature.
Molecules 24 03478 g005
Figure 6. An exemplary view of the input generation function of the transitivity program. Details of input files can be found in the www.vhcsgroup.com/transitivity web page.
Figure 6. An exemplary view of the input generation function of the transitivity program. Details of input files can be found in the www.vhcsgroup.com/transitivity web page.
Molecules 24 03478 g006
Table 1. Fitted parameters for the Arrhenius, AM, ASCC, NTS and VFT formulas, using the Transitivity code for keto–enol tautomerization [7], OH + H2 [76], enzymatic catalysis [78] and OH + Br [84] reactions. Energy ( E a , ε , E υ and E 0 ) is in cal/mol and temperature ( T 0   and B ) in K. Pre-factor units can be identified in the references.
Table 1. Fitted parameters for the Arrhenius, AM, ASCC, NTS and VFT formulas, using the Transitivity code for keto–enol tautomerization [7], OH + H2 [76], enzymatic catalysis [78] and OH + Br [84] reactions. Energy ( E a , ε , E υ and E 0 ) is in cal/mol and temperature ( T 0   and B ) in K. Pre-factor units can be identified in the references.
Formula Chemical Processes
Fitted
Parameters
Keto-enol Tautomerization [7] Sub -Arrhenius
(Deep-Tunneling)
OH + H2 → H + H2
[76] Sub-Arrhenius
(Moderate Tunneling)
Enzymatic Catalysis
[78] Super-Arrhenius
OH + HBr → Br + H2O
[84] Anti-Arrhenius
Arrhenius
k ( T ) = A exp ( E a k B T )
A 1.74 × 1032.16 · 10-111.52 × 10111.66 × 10-11
E a 214489114600−94.6
χ 2 1.10 × 10-24.20 · 10-32.60 × 10-26.69 × 10-2
Aquilanti–Mundim (AM)
k ( T ) = A ( 1 d ε k B T ) 1 d
A 3.32 × 1061.11 · 10-101.91 × 1047.43 × 10-14
ε 318.0691702391−324.61
d −0.81−0.0860.2071.24
χ 2 3.68 × 10-26.80 · 10-42.91 × 10-22.78 × 10-3
Aquilanti–Sanchez–Coutinho–Carvalho (ASCC)
k ( T ) = A ( 1 d ε k B T + E ν ) 1 d ,   d = 1 3 ( E υ 2 ε ) 2
A 2.33 × 104---
ε 2441---
E υ 429---
χ 2 2.18 × 10-2---
Sato–Nakamura–Takayanagi (NTS)
k ( T ) = A exp [ E 0 k B ( T 2 + T 0 2 ) 1 2 ]
A 3.12 × 104---
E 0 1655---
T 0 168---
χ 2 7.38 × 10-3---
Vogel–Fulcher–Tammann (VFT)
k ( T ) = A exp ( B T T 0 )
A --1.25 × 105-
B --−1298-
T 0 --175-
χ 2 --2.16 × 10-2-

Share and Cite

MDPI and ACS Style

Machado, H.G.; Sanches-Neto, F.O.; Coutinho, N.D.; Mundim, K.C.; Palazzetti, F.; Carvalho-Silva, V.H. “Transitivity”: A Code for Computing Kinetic and Related Parameters in Chemical Transformations and Transport Phenomena. Molecules 2019, 24, 3478. https://doi.org/10.3390/molecules24193478

AMA Style

Machado HG, Sanches-Neto FO, Coutinho ND, Mundim KC, Palazzetti F, Carvalho-Silva VH. “Transitivity”: A Code for Computing Kinetic and Related Parameters in Chemical Transformations and Transport Phenomena. Molecules. 2019; 24(19):3478. https://doi.org/10.3390/molecules24193478

Chicago/Turabian Style

Machado, Hugo G., Flávio O. Sanches-Neto, Nayara D. Coutinho, Kleber C. Mundim, Federico Palazzetti, and Valter H. Carvalho-Silva. 2019. "“Transitivity”: A Code for Computing Kinetic and Related Parameters in Chemical Transformations and Transport Phenomena" Molecules 24, no. 19: 3478. https://doi.org/10.3390/molecules24193478

APA Style

Machado, H. G., Sanches-Neto, F. O., Coutinho, N. D., Mundim, K. C., Palazzetti, F., & Carvalho-Silva, V. H. (2019). “Transitivity”: A Code for Computing Kinetic and Related Parameters in Chemical Transformations and Transport Phenomena. Molecules, 24(19), 3478. https://doi.org/10.3390/molecules24193478

Article Metrics

Back to TopTop