Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations

Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.

The Nature of the Enthalpy-Entropy Compensation and "Exotic" Arrhenius Parameters in the Denaturation Kinetics of Proteins

Protein unfolding is a ubiquitous process responsible for the loss of protein functionality (denaturation), which, in turn, can be accompanied by the death of cells and organisms. The nature of enthalpy-entropy compensation (EEC) in the kinetics of protein unfolding is a subject of debate. In order to investigate the nature of EEC, the "completely loose" transition state (TS) model has been applied to calculate the Arrhenius parameters for the unfolding of polyglycine dimers as a model process. The calculated Arrhenius parameters increase with increasing dimer length and demonstrate enthalpy-entropy compensation. It is shown that EEC results from the linear correlations of enthalpy and entropy of activation with dimer length, which are derived directly from the properties of the transition state. It is shown that EEC in solvated (hydrated, etc.) proteins is a direct consequence of EEC in proteins themselves. The suggested model allows us also to reproduce and explain "exotic" very high values of the pre-exponential factor measured for the proteins unfolding, which are drastically higher than those known for unimolecular reactions of organic molecules. A similar approach can be applied to analyzing the nature of EEC phenomena observed in other areas of chemistry.

NO Selective Catalytic Reduction over Atom-Pair Active Sites Accelerated via In Situ NO Oxidation

Selective catalytic reduction (SCR) of NO_x with NH₃ is the most efficient technology for NO_x emissions control, but the activity of catalysts decreases exponentially with the decrease in reaction temperature, hindering the application of the technology in low-temperature SCR to treat industrial stack gases. Here, we present an industrially practicable technology to significantly enhance the SCR activity at low temperatures (<250 °C). By introducing an appropriate amount of O₃ into the simulated stack gas, we find that O₃ can stoichiometrically oxidize NO to generate NO₂, which enables NO reduction to follow the fast SCR mechanism so as to accelerate SCR at low temperatures, and, in particular, an increase in SCR rate by more than four times is observed over atom-pair V₁-W₁ active sites supported on TiO₂(001) at 200 °C. Using operando SCR tests and in situ diffuse reflectance infrared Fourier transform spectra, we reveal that the introduction of O₃ allows SCR to proceed along a NH₄NO₃-mediated Langmuir-Hinshelwood model, in which the adsorbed nitrate species speed up the re-oxidation of the catalytic sites that is the rate-limiting step of SCR, thus leading to the enhancement of activity at low temperatures. This technology could be applicable in the real stack gas conditions because O₃ exclusively oxidizes NO even in the co-presence of SO₂ and H₂O, which provides a general strategy to improve low-temperature SCR efficacy from another perspective beyond designing catalysts.

Study of the Kinetics of Reduction of IrO2 on TiO2 (Anatase) by Temperature-Programmed Reduction

The interaction between IrO2 and TiO2 (anatase) in non-isothermal reduction conditions has been studied by the temperature programmed reduction technique. IrO2 clusters are of sizes between 0.5 and 0.9 nm as determined from High Resolution Transmission Electron Microscopy (HRTEM). Largely, two main regions for reduction were found and modeled at ca. 100 and 230 °C. The first region is attributed to the partial reduction of IrO2 clusters, while the second one is due to reduction of the formed crystalline (rutile IrO2), during TPR, to Ir metal. Two methods for calculating kinetic parameters were tested. First, by applying different ramping rates on a 3.5 wt.% IrO2/TiO2 using Kissinger’s method. The apparent activation energy values for the first and second reduction regions were found to be ca. 35 and 100 kJ/mol, respectively. The second method was based on fitting different kinetic models for the experimental results in order to extract qualitative information on the nature of interaction during the reduction process. It was found that the first reduction is largely due to the amount of IrO2 (reactant concentration) while the second one involved phase boundary effect as well as nucleation.

Characterization of polar surface groups on siliceous materials by inverse gas chromatography and the enthalpy-entropy compensation effect

Surface-modified porous silica is a well-established composite material. To improve its embedding and application behavior, adsorption studies of various probe molecules have been performed using the technique of inverse gas chromatography (IGC). For this purpose, IGC experiments were carried out in the infinite dilution mode on macro-porous micro glass spheres before and after surface modification with (3-mercaptopropyl)trimethoxysilane. To provide information about the polar interactions between probe molecules and the silica surface, in particular, eleven polar molecules have been injected. In summary, the free surface energy for pristine silica ( γ S t o t a l = 229 mJ/m²) and for (3-mercaptopropyl)trimethoxysilane-modified silica ( γ S t o t a l = 135 mJ/m²) indicates a reduced wettability after surface modification. This is due to the reduction of the polar component of the free surface energy ( γ S S P ) from 191 mJ/m² to 105 mJ/m². Simultaneously, with the reduction of surface silanol groups caused by surface modification of silica and, therefore, the decrease in polar interactions, a substantial loss of Lewis acidity was observed by various IGC approaches. Experiments with all silica materials have been conducted at temperatures in the range from 90°C to 120°C to determine the thermodynamic parameters, such as adsorption enthalpy ( Δ H a d s ) and adsorption entropy ( Δ S a d s ), using the Arrhenius regression procedure evaluating the IGC data. With the help of the enthalpy-entropy compensation, two types of adsorption complexes are assumed between polar probe molecules and the silica surface because of different isokinetic temperatures. Identical adsorption complexes with an isokinetic temperature of 370°C have been assigned to alkanes and weakly interacting polar probes such as benzene, toluene, dichloromethane, and chloroform. Polar probe molecules with typical functional groups such as OH, CO, and CN, having the ability to form hydrogen bonds to the silica surface, exhibit a lower isokinetic temperature of 60°C. Quantum chemical calculations of the probe molecules on a non-hydroxylated and hydroxylated silica cluster supported the formation of hydrogen bonds in the case of a strong polar adsorption complex with a bonding distance of 1.7 nm-1.9 nm to the silica surface.

An Atom-Pair Design Strategy for Optimizing the Synergistic Electron Effects of Catalytic Sites in NO Selective Reduction

Effective adsorption and speedy surface reactions are vital requirements for efficient active sites in catalysis, but it remains challenging to maximize these two functions simultaneously. We present a solution to this issue by designing a series of atom-pair catalytic sites with tunable electronic interactions. As a case study, NO selective reduction occurring on V₁ -W₁ /TiO₂ is chosen. Experimental and theoretical results reveal that the synergistic electron effect present between the paired atoms enriches high-energy spin charge around the Fermi level, simultaneously rendering reactant (NH₃ or O₂ ) adsorption more effective and subsequent surface reactions speedier as compared with single V or W atom alone, and hence higher reaction rates. This strategy enables us to rationally design a high-performance V₁ -Mo₁ /TiO₂ catalyst with optimized vanadium(IV)-molybdenum(V) electronic interactions, which has exceptional activity significantly higher than the commercial or reported catalysts.

Increasing the Distance of Adjacent Palladium Atoms for Configuration Matching in Selective Hydrogenation

Precisely tailoring the distance between adjacent metal sites to match adsorption configurations of key species for the targeted reaction pathway is a great challenge in heterogeneous catalysis. Here, we report a proof-of-concept study on the atomically sites-tailored pathway in Pd-catalyzed acetylene hydrogenation, i.e., increasing the distance of adjacent Pd atoms (d_Pd-a-Pd ) for configuration matching in acetylene semi-hydrogenation against coupling. d_Pd-a-Pd is identified as a structural descriptor for describing the competitiveness for reaction pathways, and the increased d_Pd-a-Pd prefers the semi-hydrogenation pathway due to simultaneously promoted C₂ H₄ desorption and the destabilized transition state of the C₂ H₃ * coupling. Spectroscopic, kinetics and electronic structure studies reveal that increasing d_Pd-a-Pd to 3.31 Å delivers superior selectivity and stability due to energy matching and appropriate hybridization of Pd 4d with In 2s and, especially, 2p orbitals.

MOF/POM hybrids as catalysts for organic transformations

Insertion of molecular metal oxides, e.g. polyoxometalates (POMs), into metal-organic frameworks (MOFs) opens up new research opportunities in various fields, particularly in catalysis. POM/MOF composites have strong acidity, oxygen-rich surface, and redox capacity due to typical characteristics of POMs and the large surface area, highly organized structures, tunable pore size, and shape are due to MOFs. Such hybrid materials have gained a lot of attention due to astonishing structural features, and hence have potential applications in organic catalysis, sorption and separation, proton conduction, magnetism, lithium-ion batteries, supercapacitors, electrochemistry, medicine, bio-fuel, and so on. The exceptional chemical and physical characteristics of POMOFs make them useful as catalysts in simple organic transformations with high capacity and selectivity. Here, the thorough catalytic study starts with a brief introduction related to POMs and MOFs, and is followed by the synthetic strategies and applications of these materials in several catalytic organic transformations. Furthermore, catalytic conversions like oxidation, condensation, esterification, and some other types of catalytic reactions including photocatalytic reactions are discussed in length with their plausible catalytic mechanisms. The disadvantages of the POMOFs and difficulties faced in the field have also been explored briefly from our perspectives.

Grain-growth mediated hydrogen sorption kinetics and compensation effect in single Pd nanoparticles

Grains constitute the building blocks of polycrystalline materials and their boundaries determine bulk physical properties like electrical conductivity, diffusivity and ductility. However, the structure and evolution of grains in nanostructured materials and the role of grain boundaries in reaction or phase transformation kinetics are poorly understood, despite likely importance in catalysis, batteries and hydrogen energy technology applications. Here we report an investigation of the kinetics of (de)hydriding phase transformations in individual Pd nanoparticles. We find dramatic evolution of single particle grain morphology upon cyclic exposure to hydrogen, which we identify as the reason for the observed rapidly slowing sorption kinetics, and as the origin of the observed kinetic compensation effect. These results shed light on the impact of grain growth on kinetic processes occurring inside nanoparticles, and provide mechanistic insight in the observed kinetic compensation effect.

Grains are the building blocks of crystalline solids. Here the authors show how hydrogen-sorption induced grain-growth in Pd nanoparticles slows down the hydrogen sorption kinetics and constitutes the physical origin of corresponding kinetic compensation.

Structure and Stability of Gas Adsorption Complexes in Periodic Porous Solids as Studied by VTIR Spectroscopy: An Overview

Variable-temperature infrared (VTIR) spectroscopy is an instrumental technique that enables structural characterization of gas-solid adsorption complexes by analysis of meaningful vibrational modes, and simultaneous determination of the standard enthalpy change (ΔH0) involved in the gas adsorption process, which allows one to quantify the stability of the corresponding complex. This is achieved by a van’t Hoff analysis of a set of IR spectra recorded over a sufficiently large temperature range. Herein, the use of this versatile spectroscopic technique is demonstrated by reviewing its application to the study of carbon monoxide, carbon dioxide and dinitrogen adsorption on several (alkaline) zeolites, which can be regarded as the archetype of periodic porous solids.

Single Quality Factor for Enthalpy-Entropy Compensation, Isoequilibrium and Isokinetic Relationships

Enthalpy-entropy compensation (EEC) is very often encountered in chemistry, biology and physics. Its origin is widely discussed since it would allow, for example, a very accurate tuning of the thermodynamic properties as a function of the reactants. However, EEC is often discarded as a statistical artefact, especially when only a limited temperature range is considered. We show that the likeliness of a statistical origin of an EEC can be established with a compensation quality factor (CQF) that depends only on the measured enthalpies and entropies and the experimental temperature range. This is directly derived from a comparison of the CQF with threshold values obtained from a large number of simulations with randomly generated Van 't Hoff plots. The value of CQF is furthermore a direct measure of the existence of a genuine isoequilibrium or isokinetic relationship.

Revealing Kinetics of Two-Electron Oxygen Reduction Reaction at Single-Molecule Level

By combining single-molecule fluorescence microscopy with traditional electrochemical methods, herein we report on the investigation of the electrocatalytic kinetics of two-electron (2e) pathway of oxygen reduction reaction (ORR) on a single Fe₃O₄ nanoparticle. The kinetic parameters for two-electron ORR process are successfully derived at the single-particle level, and a potential dependence of dynamic heterogeneity among individual nanoparticles is revealed. Furthermore, the performance stability of individual Fe₃O₄ nanoparticles for 2e ORR process is studied. This study deepens our understanding to the electrocatalytic ORR process, especially the 2e pathway at single-molecule and single-particle levels.

Kinetic Analysis of 4-Nitrophenol Reduction by "Water-Soluble" Palladium Nanoparticles

The most important model catalytic reaction to test the catalytic activity of metal nanoparticles is the reduction of 4-nitrophenol to 4-aminophenol by sodium borohydride as it can be precisely monitored by UV–vis spectroscopy with high accuracy. This work presents the catalytic reduction of 4-nitrophenol (4-Nip) to 4-aminophenol (4-Amp) in the presence of Pd nanoparticles and sodium borohydride as reductants in water. We first evaluate the kinetics using classical pseudo first-order kinetics. We report the effects of different initial 4-Nip and NaBH₄ concentrations, reaction temperatures, and mass of Pd nanoparticles used for catalytic reduction. The thermodynamic parameters (activation energy, enthalpy, and entropy) were also determined. Results show that the kinetics are highly dependent on the reactant ratio and that pseudo first-order simplification is not always fit to describe the kinetics of the reaction. Assuming that all steps of this reaction proceed only on the surface of Pd nanoparticles, we applied a Langmuir−Hinshelwood model to describe the kinetics of the reaction. Experimental data of the decay rate of 4-nitrophenol were successfully fitted to the theoretical values obtained from the Langmuir–Hinshelwood model and all thermodynamic parameters, the true rate constant k, as well as the adsorption constants of 4-Nip, and BH₄⁻ (K_4-Nip and K_BH4−) were determined for each temperature.

Atomically Dispersed Pt1-Polyoxometalate Catalysts: How Does Metal-Support Interaction Affect Stability and Hydrogenation Activity?

Unlike nanostructured metal catalysts, supported single-atom catalysts (SACs) contain only atomically dispersed metal atoms, hinting at much more pronounced metal-support effects. Herein, we take a series of polyoxometalate-supported Pt catalysts as examples to quantitatively investigate the stability of Pt atoms on oxide supports and how the Pt-support interaction influences the catalytic performance. For this entire series, we show that the Pt atoms prefer to stay at a 4-fold hollow site of one polyoxometalate molecule and that the least adsorption energy to obtain sintering-resistant Pt SACs is 5.50 eV, which exactly matches the cohesive energy of bulk Pt metal. Further, we compared their catalytic performance in several hydrogenation reactions and simulated the reaction pathways of propene hydrogenation by density functional theory (DFT) calculations. Both experimental and theoretical approaches suggest that despite the Pt₁-support interactions being different, the reaction pathways of various Pt₁-polyoxometalate catalysts are very similar and their effective reaction barriers are close to each other and as low as 24 kJ/mol, indicating the possibility of obtaining SACs with improved stability without compromising activity. DFT calculations show that all reaction elementary steps take place only on the Pt atom without involving neighboring O atoms and that hydrogenation proceeds from the molecularly adsorbed H₂ species. Pt SACs give a weaker H₂ adsorption energy than Pt clusters or surfaces, resulting in small adsorption equilibrium constants and small apparent activation barriers, which agree between experiment and theory.

Heterodimers for in Situ Plasmonic Spectroscopy: Cu Nanoparticle Oxidation Kinetics, Kirkendall Effect, and Compensation in the Arrhenius Parameters

The ability to study oxidation, reduction, and other chemical transformations of nanoparticles in real time and under realistic conditions is a nontrivial task due to their small dimensions and the often challenging environment in terms of temperature and pressure. For scrutinizing oxidation of metal nanoparticles, visible light optical spectroscopy based on the plasmonic properties of the metal has been established as a suitable method. However, directly relying on the plasmonic resonance of metal nanoparticles as a built-in probe to track oxidation has a number of drawbacks, including the loss of optical contrast in the late oxidation stages. To address these intrinsic limitations, we present a plasmonic heterodimer-based nanospectroscopy approach, which enables continuous self-referencing by using polarized light to eliminate parasitic signals and provides large optical contrast all the way to complete oxidation. Using Au-Cu heterodimers and combining experiments with finite-difference time-domain simulations, we quantitatively analyze the oxidation kinetics of ca. 30 nm sized Cu nanoparticles up to complete oxidation. Taking the Kirkendall effect into account, we extract the corresponding apparent Arrhenius parameters at various extents of oxidation and find that they exhibit a significant compensation effect, implying that changes in the oxidation mechanism occur as oxidation progresses and the structure of the formed oxide evolves. In a wider perspective, our work promotes the use of model-system-type in situ optical plasmonic spectroscopy experiments in combination with electrodynamics simulations to quantitatively analyze and mechanistically interpret oxidation of metal nanoparticles and the corresponding kinetics in demanding chemical environments, such as in heterogeneous catalysis.

Non-equilibrium thermodynamics as a tool to compute temperature at the catalyst surface

The NET theory predicts the coupling between reaction rates and thermal driving forces and gives new insights into why Arrhenius plots may turn out to be non-linear.

Photo-thermal reactions of ethanol over Ag/TiO2 catalysts. The role of silver plasmon resonance in the reaction kinetics

Photo-thermal catalytic reactions of ethanol over Ag/TiO₂ were conducted in order to probe into the role of plasmonic resonance response in the reaction kinetics.

Improving Kinetics of "Click-Crosslinking" for Self-Healing Nanocomposites by Graphene-Supported Cu-Nanoparticles

Investigation of the curing kinetics of crosslinking reactions and the development of optimized catalyst systems is of importance for the preparation of self-healing nanocomposites, able to significantly extend their service lifetimes. Here we study different modified low molecular weight multivalent azides for a capsule-based self-healing approach, where self-healing is mediated by graphene-supported copper-nanoparticles, able to trigger "click"-based crosslinking of trivalent azides and alkynes. When monitoring the reaction kinetics of the curing reaction via reactive dynamic scanning calorimetry (DSC), it was found that the "click-crosslinking" reactivity decreased with increasing chain length of the according azide. Additionally, we could show a remarkable "click" reactivity already at 0 °C, highlighting the potential of click-based self-healing approaches. Furthermore, we varied the reaction temperature during the preparation of our tailor-made graphene-based copper(I) catalyst to further optimize its catalytic activity. With the most active catalyst prepared at 700 °C and the optimized set-up of reactants on hand, we prepared capsule-based self-healing epoxy nanocomposites.

Probing Gas Adsorption in Zeolites by Variable-Temperature IR Spectroscopy: An Overview of Current Research

The current state of the art in the application of variable-temperature IR (VTIR) spectroscopy to the study of (i) adsorption sites in zeolites, including dual cation sites; (ii) the structure of adsorption complexes and (iii) gas-solid interaction energy is reviewed. The main focus is placed on the potential use of zeolites for gas separation, purification and transport, but possible extension to the field of heterogeneous catalysis is also envisaged. A critical comparison with classical IR spectroscopy and adsorption calorimetry shows that the main merits of VTIR spectroscopy are (i) its ability to provide simultaneously the spectroscopic signature of the adsorption complex and the standard enthalpy change involved in the adsorption process; and (ii) the enhanced potential of VTIR to be site specific in favorable cases.

Surface single-molecule dynamics controlled by entropy at low temperatures

Configuration transitions of individual molecules and atoms on surfaces are traditionally described using an Arrhenius equation with energy barrier and pre-exponential factor (attempt rate) parameters. Characteristic parameters can vary even for identical systems, and pre-exponential factors sometimes differ by orders of magnitude. Using low-temperature scanning tunnelling microscopy (STM) to measure an individual dibutyl sulfide molecule on Au(111), we show that the differences arise when the relative position of tip apex and molecule changes by a fraction of the molecule size. Altering the tip position on that scale modifies the transition's barrier and attempt rate in a highly correlated fashion, which results in a single-molecular enthalpy-entropy compensation. Conversely, appropriately positioning the STM tip allows selecting the operating point on the compensation line and modifying the transition rates. The results highlight the need to consider entropy in transition rates of single molecules, even at low temperatures.

STM is capable of imaging the configurations of molecules on surfaces and measuring the rate of transitions between them. Here the authors demonstrate that, controlled by the STM tip position, the entropic and conservative forces on the molecule can modify the rate by orders of magnitude.

Hydrogenolysis of alkanes: reactions of n-butane on Ru/zeolite catalysts

Ru/13X zeolite catalysts containing 0.6, 2 and 5.7% Ru prepared by reduction or decomposition of the hexammine complex are active and stable for n-butane hydrogenolysis between 400 and 460 K; oxidation leads to greatly increased activity.

Improving the catalytic efficiency of carbon-based active sites by trace oxide promoters for highly productive olefin synthesis

The catalytic ODH activation energy on CNTs significantly decreased by the introduction of metal oxides. The propane turn-over efficiency of each original carbon-based active site consequently increased exponentially.

Direct inference of site strength in basic solids upon CO2 adsorption: enthalpy-entropy compensation effects

The adsorption of CO2 coupled to calorimetry is a state-of-the-art technique for characterizing the basic properties of solids. In this paper, we show that the differential heat and entropy curves measured upon CO2 adsorption on a basic solid can be reasonably estimated from a single CO2 isotherm with no need for any independent heat (calorimetric) measurement. Our method relies on two important observations: (1) formulation of generalized F-H-TS thermodynamic isotherms, the former (F) being directly generated from the raw CO2 isotherms, and (2) the presence of unexpected enthalpy-entropy compensation effects upon CO2 adsorption linking the integral enthalpy and entropy of adsorption until saturation for different solids. Our thermodynamic method has been validated using a broad library of basic solids with variable site strength and heterogeneity. Finally, a new scale of basicity is proposed using the parameters fitted from the thermodynamic isotherm (free energy basis) as descriptors of basic strength. This method opens an avenue to the inference of site strength of basic solids without the need for expensive calorimeters.

Nickel–antimony nanoparticles confined in SBA-15 as highly efficient catalysts for the hydrogenation of nitroarenes

Nickel–antimony nanoparticles confined in SBA-15 are found to be highly efficient catalysts for nitroarene reduction reactions.