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

Pharmacokinetic-Pharmacodynamic Analysis of pH-Responsive Doxorubicin-Releasing Micelles with Anticancer Activity

This study aimed to investigate the effect of in vivo pH-responsive doxorubicin (DOX) release and the targetability of pilot molecules in folic acid (FA)-modified micelles using a pharmacokinetic-pharmacodynamic (PK-PD) model. The time profiles of intratumoral DOX concentrations in Walker256 tumor-bearing rats were monitored using a microdialysis probe, followed by compartmental analysis, to evaluate intratumoral tissue pharmacokinetics. Maximal DOX was released from micelles 350 min after the administration of pH-responsive DOX-releasing micelles. However, FA modification of the micelles shortened the time to peak drug concentration to 150 min. Additionally, FA modification resulted in a 27-fold increase in the tumor inflow rate constant. Walker256 tumor-bearing rats were subsequently treated with DOX, pH-responsive DOX-releasing micelles, and pH-responsive DOX-releasing FA-modified micelles to monitor the tumor growth-time profiles. An intratumoral threshold concentration of DOX (55-64 ng/g tumor) was introduced into the drug efficacy compartment to construct a PD model, followed by PK-PD analysis of the tumor growth-time profiles. Similar results of threshold concentration and drug potency of DOX were obtained across all three formulations. Cell proliferation was delayed as the drug delivery ability of DOX was improved. The PK model, which was developed using the microdialysis method, revealed the intratumoral pH-responsive DOX distribution profiles. This facilitated the estimation of intratumoral PK parameters. The PD model with threshold concentrations contributed to the estimation of PD parameters in the three formulations, with consistent mechanisms observed. We believe that our PK-PD model can objectively assess the contributions of pH-responsive release ability and pilot molecule targetability to pharmacological effects.

A Maxwell relation for dynamical timescales with application to the pressure and temperature dependence of water self-diffusion and shear viscosity

A Maxwell relation for a reaction rate constant (or other dynamical timescale) obtained under constant pressure, p, and temperature, T, is introduced and discussed. Examination of this relationship in the context of fluctuation theory provides insight into the p and T dependence of the timescale and the underlying molecular origins. This Maxwell relation motivates a suggestion for the general form of the timescale as a function of pressure and temperature. This is illustrated by accurately fitting simulation results and existing experimental data on the self-diffusion coefficient and shear viscosity of liquid water. A key advantage of this approach is that each fitting parameter is physically meaningful.

Multiscale modeling of stochastic dynamics processes with MBN Explorer

Stochastic dynamics describes processes in complex systems having the probabilistic nature. They can involve very different dynamical systems and occur on very different temporal and spatial scale. This paper discusses the concept of stochastic dynamics and its implementation in the popular program MBN Explorer. Stochastic dynamics in MBN Explorer relies on the Monte Carlo approach and permits simulations of physical, chemical, and biological processes. The paper presents the basic theoretical concepts underlying stochastic dynamics implementation and provides several examples highlighting its applicability to different systems, such as diffusing proteins seeking an anchor point on a cell membrane, deposition of nanoparticles on a surface leading to structures with fractal morphologies, and oscillations of compounds in an autocatalytic reaction. The chosen examples illustrate the diversity of applications that can be modeled by means of stochastic dynamics with MBN Explorer.

On the Dynamics of the Carbon-Bromine Bond Dissociation in the 1-Bromo-2-Methylnaphthalene Radical Anion

This paper studies the mechanism of electrochemically induced carbon-bromine dissociation in 1-Br-2-methylnaphalene in the reduction regime. In particular, the bond dissociation of the relevant radical anion is disassembled at a molecular level, exploiting quantum mechanical calculations including steady-state, equilibrium and dissociation dynamics via dynamic reaction coordinate (DRC) calculations. DRC is a molecular-dynamic-based calculation relying on an ab initio potential surface. This is to achieve a detailed picture of the dissociation process in an elementary molecular detail. From a thermodynamic point of view, all the reaction paths examined are energetically feasible. The obtained results suggest that the carbon halogen bond dissociates following the first electron uptake follow a stepwise mechanism. Indeed, the formation of the bromide anion and an organic radical occurs. The latter reacts to form a binaphthalene intrinsically chiral dimer. This paper is respectfully dedicated to Professors Anny Jutand and Christian Amatore for their outstanding contribution in the field of electrochemical catalysis and electrosynthesis.

Nanocatalyst-Mediated Chemodynamic Tumor Therapy

Traditional tumor treatments, including chemotherapy, radiotherapy, photodynamic therapy, and photothermal therapy, are developed and used to treat different types of cancer. Recently, chemodynamic therapy (CDT) has been emerged as a novel cancer therapeutic strategy. CDT utilizes Fenton or Fenton-like reaction to generate highly cytotoxic hydroxyl radicals (•OH) from endogenous hydrogen peroxide (H₂ O₂ ) to kill cancer cells, which displays promising therapeutic potentials for tumor treatment. However, the low catalytic efficiency and off-target side effects of Fenton reaction limit the biomedical application of CDT. In this regard, various strategies are implemented to potentiate CDT against tumor, including retrofitting the tumor microenvironment (e.g., increasing H₂ O₂ level, decreasing reductive substances, and reducing pH), enhancing the catalytic efficiency of nanocatalysts, and other strategies. This review aims to summarize the development of CDT and summarize these recent progresses of nanocatalyst-mediated CDT for antitumor application. The future development trend and challenges of CDT are also discussed.

"pySiRC": Machine Learning Combined with Molecular Fingerprints to Predict the Reaction Rate Constant of the Radical-Based Oxidation Processes of Aqueous Organic Contaminants

We developed a web application structured in a machine learning and molecular fingerprint algorithm for the automatic calculation of the reaction rate constant of the oxidative processes of organic pollutants by ^•OH and SO₄^•- radicals in the aqueous phase-the pySiRC platform. The model development followed the OECD principles: internal and external validation, applicability domain, and mechanistic interpretation. Three machine learning algorithms combined with molecular fingerprints were evaluated, and all the models resulted in high goodness-of-fit for the training set with R² > 0.931 for the ^•OH radical and R² > 0.916 for the SO₄^•- radical and good predictive capacity for the test set with R_ext² = Q_ext² values in the range of 0.639-0.823 and 0.767-0.824 for the ^•OH and SO₄^•- radicals. The model was interpreted using the SHAP (SHapley Additive exPlanations) method: the results showed that the model developed made the prediction based on a reasonable understanding of how electron-withdrawing and -donating groups interfere with the reactivity of the ^•OH and SO₄^•- radicals. We hope that our models and web interface can stimulate and expand the application and interpretation of kinetic research on contaminants in water treatment units based on advanced oxidative technologies.

Non-Arrhenius kinetics and slowed-diffusion mechanism of molecular aggregation of a rhodamine dye on colloidal particles

The non-covalent association is important for many fields of science, including processes in living systems. This work elucidates the mechanism of rhodamine 123 molecular aggregation in dispersions of a layered silicate and explains the mystery of the slow kinetics of this process. Chemometric analysis of thousands of spectra recorded by stopped-flow visible spectroscopy identified two parallel diffusion processes described by a two-phase exponential function. The slow and fast processes followed the super-Arrhenius kinetics and were assigned to lateral (surface) diffusion and inter-particle diffusion of dye cations, respectively. This work, supported by a large amount of data and their in-depth analysis, provides the first evidence of how these processes coexist together and provides quantitative analysis of their dependence on the reaction conditions. The implications of this work can be crucial for understanding the mechanism of the non-covalent association of adsorbed molecules in nature.

Study of deactivation in mesocellular foam carbon (MCF-C) catalyst used in gas-phase dehydrogenation of ethanol

Mesocellular foam carbon (MCF-C) is one the captivating materials for using in gas phase dehydrogenation of ethanol. Extraordinary, enlarge pore size, high surface area, high acidity, and spherical shape with interconnected pore for high diffusion. In contrary, the occurrence of the coke is a majority causes for inhibiting the active sites on catalyst surface. Thus, this study aims to investigate the occurrence of the coke to optimize the higher catalytic activity, and also to avoid the coke formation. The MCF-C was synthesized and investigated using various techniques. MCF-C was spent in gas-phase dehydrogenation of ethanol under mild conditions. The deactivation of catalyst was investigated toward different conditions. Effects of reaction condition including different reaction temperatures of 300, 350, and 400 °C on the deactivation behaviors were determined. The results indicated that the operating temperature at 400 °C significantly retained the lowest change of ethanol conversion, which favored in the higher temperature. After running reaction, the physical properties as pore size, surface area, and pore volume of spent catalysts were decreased owing to the coke formation, which possibly blocked the pore that directly affected to the difficult diffusion of reactant and caused to be lower in catalytic activity. Furthermore, a slight decrease in either acidity or basicity was observed owing to consumption of reactant at surface of catalyst or chemical change on surface caused by coke formation. Therefore, it can remarkably choose the suitable operating temperature to avoid deactivation of catalyst, and then optimize the ethanol conversion or yield of acetaldehyde.

Topography of the free energy landscape of Claisen-Schmidt condensation: solvent and temperature effects on the rate-controlling step

Recent studies have found that hydroxide elimination and the C[double bond, length as m-dash]C bond formation step in base-promoted aldol condensation have a strong influence on the overall rate of the reaction, in contrast to the well-accepted first enolization or C-C bond formation step. Here, applying theoretical models to the prototypical reaction of chalcone formation, the complete free energy profile of Claisen-Schmidt condensation is assessed, revealing how a protic solvent and a slight increase in temperature can induce the second enolization as the rate-controlling step (RCS). It is also observed: i) the nonexistence of a step with a much higher energetic barrier than the others, making the concept of RCS debatable; and ii) that the overall inverse kinetic isotopic effect does not exclude second enolization as a RCS in protic continuum medium. We expect that these results can expand the understanding of the decisive role of physicochemical factors on the choose of the RCS in the aldol condensation.

Temperature effect on water dynamics in tetramer phosphofructokinase matrix and the super-arrhenius respiration rate

Advances in understanding the temperature effect on water dynamics in cellular respiration are important for the modeling of integrated energy processes and metabolic rates. For more than half a century, experimental studies have contributed to the understanding of the catalytic role of water in respiration combustion, yet the detailed water dynamics remains elusive. We combine a super-Arrhenius model that links the temperature-dependent exponential growth rate of a population of plant cells to respiration, and an experiment on isotope labeled ¹⁸O₂ uptake to H₂¹⁸O transport role and to a rate-limiting step of cellular respiration. We use Phosphofructokinase (PFK-1) as a prototype because this enzyme is known to be a pacemaker (a rate-limiting enzyme) in the glycolysis process of respiration. The characterization shows that PFK-1 water matrix dynamics are crucial for examining how respiration (PFK-1 tetramer complex breathing) rates respond to temperature change through a water and nano-channel network created by the enzyme folding surfaces, at both short and long (evolutionary) timescales. We not only reveal the nano-channel water network of PFK-1 tetramer hydration topography but also clarify how temperature drives the underlying respiration rates by mapping the channels of water diffusion with distinct dynamics in space and time. The results show that the PFK-1 assembly tetramer possesses a sustainable capacity in the regulation of the water network toward metabolic rates. The implications and limitations of the reciprocal-activation-reciprocal-temperature relationship for interpreting PFK-1 tetramer mechanisms are briefly discussed.

Deciphering the Dynamics of Organic Nanoaggregates with AIEE Effect and Excited States: Lipophilic Benzothiadiazole Derivatives as Selective Cell Imaging Probes

An aggregation-induced emission enhancement (AIEE) effect in fluorescent lipophilic 2,1,3-benzothiadiazole (BTD) derivatives and their organic nanoaggregates were studied. A set of techniques such as single-crystal X-ray, dynamic light scattering (DLS), electron paramagnetic resonance (EPR), UV-vis, fluorescence, and density functional theory (DFT) calculations have been used to decipher the formation/break (kinetics), properties, and dynamics of the organic nanoaggregates of three BTD small organic molecules. An in-depth study of the excited-state also revealed the preferential relaxation emissive pathways for the BTD derivatives and the dynamics associated with it. The results described herein, for the first time, explain the formation of fluorescent BTD nanoaggregate derivatives and allow for the understanding of their dynamics in solution as well as the ruling forces of both aggregation and break processes along with the involved equilibrium. One of the developed dyes could be used at a nanomolar concentration to selectively stain lipid droplets emitting an intense and bright fluorescence at the red channel. The other two BTDs could also stain lipid droplets at very low concentrations and were visualized preferentially at the blue channel.

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.

Kinetics of Methyl Lactate Formation from the Transesterification of Polylactic Acid Catalyzed by Zn(II) Complexes

The kinetics of the transesterification of polylactic acid (PLA) with methanol to form methyl lactate catalyzed by Zn(II) complexes was studied experimentally and numerically. The complexes, Zn(1 ^Et )₂ and Zn(2 ^Pr )₂, were synthesized from ethylenediamine and propylenediamine Schiff bases, respectively. The temperature range covered was 313.2-383.2 K. An increase in the reaction rate with the increase in temperature was observed for the Zn(1 ^Et )₂-catalyzed reaction. The temperature relationship of the rate coefficients can be explained by a linear Arrhenius dependency with constant activation energy. The kinetics of Zn(2 ^Pr )₂, on the other hand, is only explained by non-Arrhenius kinetics with convex variable activation energy, resulting in faster methyl lactate production rates at 323.2 and 343.2 K. The formation of a new catalyst species, likely through reaction with protic reagents, appears to promote the formation of intermediate complexes, resulting in the nonlinear behavior. Stirring speed induced the stability of the intermediate complexes. Contrary to Zn(1 ^Et )₂, Zn(2 ^Pr )₂ was susceptible to the presence of air/moisture in solution. The kinetic parameters were obtained by fitting the experimental data to the mass and energy balance of a consecutive second step reversible reaction taking place in a jacketed stirred batch reactor. For the case of Zn(2 ^Pr )₂, the activation energy was fitted to a four-parameter equation. The kinetic parameters presented in this work are valuable for the design of processes involving the chemical recycling of PLA into green solvents.

Simulation Study to Investigate the Effects of Operational Conditions on Methylcyclohexane Dehydrogenation for Hydrogen Production

In the recent era, hydrogen has gained immense consideration as a clean-energy carrier. Its storage is, however, still the main hurdle in the implementation of a hydrogen-based clean economy. Liquid organic hydrogen carriers (LOHCs) are a potential option for hydrogen storage in ambient conditions, and can contribute to the clean-fuel concept in the future. In the present work, a parametric and simulation study was carried out for the storage and release of hydrogen for the methylcyclohexane toluene system. In particular, the methylcyclohexane dehydrogenation reaction is investigated over six potential catalysts for the temperature range of 300–450 °C and a pressure range of 1–3 bar to select the best catalyst under optimum operating conditions. Moreover, the effects of hydrogen addition in the feed mixture, and byproduct yield, are also studied as functions of operating conditions. The best catalyst selected for the process is 1 wt. % Pt/γ-Al2O3. The optimum operating conditions selected for the dehydrogenation process are 360 °C and 1.8 bar. Hydrogen addition in the feed reduces the percentage of methylcyclohexane conversion but is required to enhance the catalyst’s stability. Aspen HYSYS v. 9.0 (AspenTech, Lahore, Pakistan) has been used to carry out the simulation study.

Quasi-Classical Trajectory Study of NH(3∑-) + NH(3∑-) Reactive Collisions

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

Characterization of the non-Arrhenius behavior of supercooled liquids by modeling nonadditive stochastic systems

The characterization of the formation mechanisms of amorphous solids is a large avenue for research, since understanding its non-Arrhenius behavior is challenging to overcome. In this context, we present one path toward modeling the diffusive processes in supercooled liquids near glass transition through a class of nonhomogeneous continuity equations, providing a consistent theoretical basis for the physical interpretation of its non-Arrhenius behavior. More precisely, we obtain the generalized drag and diffusion coefficients that allow us to model a wide range of non-Arrhenius processes. This provides a reliable measurement of the degree of fragility of the system and an estimation of the fragile-to-strong transition in glass-forming liquids, as well as a generalized Stokes-Einstein equation, leading to a better understanding of the classical and quantum effects on the dynamics of nonadditive stochastic systems.

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

Beyond Boltzmann-Gibbs-Shannon in Physics and Elsewhere

The pillars of contemporary theoretical physics are classical mechanics, Maxwell electromagnetism, relativity, quantum mechanics, and Boltzmann–Gibbs (BG) statistical mechanics –including its connection with thermodynamics. The BG theory describes amazingly well the thermal equilibrium of a plethora of so-called simple systems. However, BG statistical mechanics and its basic additive entropy SBG started, in recent decades, to exhibit failures or inadequacies in an increasing number of complex systems. The emergence of such intriguing features became apparent in quantum systems as well, such as black holes and other area-law-like scenarios for the von Neumann entropy. In a different arena, the efficiency of the Shannon entropy—as the BG functional is currently called in engineering and communication theory—started to be perceived as not necessarily optimal in the processing of images (e.g., medical ones) and time series (e.g., economic ones). Such is the case in the presence of generic long-range space correlations, long memory, sub-exponential sensitivity to the initial conditions (hence vanishing largest Lyapunov exponents), and similar features. Finally, we witnessed, during the last two decades, an explosion of asymptotically scale-free complex networks. This wide range of important systems eventually gave support, since 1988, to the generalization of the BG theory. Nonadditive entropies generalizing the BG one and their consequences have been introduced and intensively studied worldwide. The present review focuses on these concepts and their predictions, verifications, and applications in physics and elsewhere. Some selected examples (in quantum information, high- and low-energy physics, low-dimensional nonlinear dynamical systems, earthquakes, turbulence, long-range interacting systems, and scale-free networks) illustrate successful applications. The grounding thermodynamical framework is briefly described as well.

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.

On the Theory of Electrolytic Dissociation, the Greenhouse Effect, and Activation Energy in (Electro)Catalysis: A Tribute to Svante Augustus Arrhenius

Svante Augustus Arrhenius (1859, Vik - 1927, Stockholm) received the Nobel Prize for Chemistry in 1903 "in recognition of the extraordinary services he rendered to the advancement of chemistry by his electrolytic theory of dissociation". Arrhenius was a physicist, and he received his PhD from the University of Uppsala, where he later became a professor for phyiscal chemistry, the first in the country for this subject. He was offered several positions as professor abroad, but decided to remain in Sweden and to build a Nobel Institute for physical chemistry using the Nobel funds. He remained director of the Institute until his death. There are powerful lessons to take from Svante August Arrhenius' journey leading to a Nobel laureate as there are from his tremendous contributions to chemistry and science in general, including climate science, immunochemistry and cosmology. The theory of electrolytic dissociation for which Arrhenius received the 1903 Nobel Prize in Chemistry has had a profound impact on our understanding of the chemistry of solutions, chemical reactivity, mechanisms underlying chemical transformations as well as physiological processes. As a tribute to Arrhenius, we present a brief historical perspective and present status of the theory of electrolytic dissociation, its relevance and role to the development of electrochemistry, as well as some perspectives on the possible role of the theory to future advancements in electroanalysis, electrocatalysis and electrochemical energy storage. The review briefly highlights Arrhenius' contribution to climate science owing to his studies on the potential effects of increased anthropogenic CO₂ emissions on the global climate. These studies were far ahead of their time and revealed a daunting global dilemma, global warming, that we are faced with today. Efforts to abate or reverse CO₂ accumulation constitute one of the most pressing scientific problems of our time, "man's urgent strive to save self from the adverse effects of his self-orchestrated change on the climate". Finally, we review the application of the Arrhenius equation that correlates reaction rate constants (k) and temperature (T); k = A e ( - E a / R T ) , in determining reaction barriers in catalysis with a particular focus on recent modifications of the equation to account for reactions exhibiting non-linear Arrhenius behavior with concave curvature due to prevalence of quantum mechanical tunneling, as well as infrequent convexity of Arrhenius plots due to decrease of the microcanonical rate coefficient with energy as observed for some enzyme catalyzed reactions.

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.

Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces

In this Introduction, we show the basic problems of non-statistical and non-equilibrium phenomena related to the papers collected in this themed issue. Over the past few years, significant advances in both computing power and development of theories have allowed the study of larger systems, increasing the time length of simulations and improving the quality of potential energy surfaces. In particular, the possibility of using quantum chemistry to calculate energies and forces 'on the fly' has paved the way to directly study chemical reactions. This has provided a valuable tool to explore molecular mechanisms at given temperatures and energies and to see whether these reactive trajectories follow statistical laws and/or minimum energy pathways. This themed issue collects different aspects of the problem and gives an overview of recent works and developments in different contexts, from the gas phase to the condensed phase to excited states.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'.