Intrinsic Reaction Research Articles

Expounding the application of nano and micro silica as a complementary additive in metakaolin phosphate geopolymer for ceramic applications—micro and nanoscale structural investigation

Metakaolin phosphate geopolymers comprising poly-phospho-siloxo units are known for their structural performance, additionally advancing their microstructure with the transformation of crystalline berlinite phases at elevated temperatures. The intrinsic reaction of Al of metakaolin in the acid exploited, but the reaction of secondary silica phases is limitedly known. Metakaolin as a primary precursor (M) with the addition of 2% and 5% of nano silica (MS2 and MS5) and micro silica (MM2 and MM5) cast using 8-M phosphoric acid was cured at 80 °C. To enhance the utilization of geopolymer in any high-temperature applications, the structural transformations were studied after heating to various temperatures (200, 400, 600 and 800 °C) by XRD, Raman, TGA-DTA, SEM, XPS, FTIR and MAS-NMR. Sample M attained a strength of 46.2 MPa enhanced to 63.6 MPa in MS5 and 54.2 MPa in MM5. This can be ascribed from the transformation of Si–O–Al–O–Si into Si–O–Al–O–P from Raman bands. Comparing the chemical shift of Al (IV) to control, micro silica addition shifts the signal to a lower field (53 to 50 ppm) related to the increase of the number of Al-connected Si to give a tougher network. Nanoindentation is visualized from hardness and elasticity, and the corresponding values are 1.4 to 2.1 GPa and 0.8 to 1.4 GPa for loads ranging from 20 to 100 mN in silica-reinforced samples that are much higher than M. The micro and macro hardness is due to the reinforcement of quartz in micro silica around the gel. TGA-DTA showed that the reduction of mass loss is as high as 25.4% in control whereas 17.2% in MS5 and 15.8% in the MM5. Further, shrinkage rate in MS5 and MM5 was as low as − 1.1% and − 0.8% throughout the temperature range from 25 to 1000 °C and thus provides the way of use of nano and micro form of silica for better thermal resistance.Graphical

Boosting the electrocatalytic performance of CuCo2S4 via surface-state engineering for ampere current water electrolysis applications

For the efficient production of green H2 on an industrial scale, developing durable, cost-effective electrocatalysts using earth-abundant transition-metals is crucial. Herein, we report a robust, bifunctional sulfurized CuCo2O4 (CCS/Ov) catalyst with engineered oxygen-vacancies, in which the metal catalyst is tunes to high valence state, significantly altering the intrinsic reaction kinetics and generating better catalytic activity for efficient ion transport in various electrolytes (neutral, seawater, alkaline simulated-seawater (ASW), and KOH). The optimized dual-strategy synthesized CCS/Ov catalysts with hierarchical morphology exhibits modest overpotential of 383 and 355 mV at 1000 mA cm−2 for OER and HER, respectively, in 1 M KOH. Impressively, an electrolyzer cell (CCS/Ov||CCS/Ov) demonstrates low cell-voltage of 1.487 V (6 M KOH), with excellent robustness in an ASW (1 and 2 M) environment against corrosive chlorine. The bifunctional CCS/Ov||CCS/Ov catalyst outperforms a state-of-the-art paired Pt/C||RuO2 catalyst and maintains the lowest potential response at up to 2000 mA cm−2, while demonstrating remarkably stable simultaneous O2 and H2 generation over 10 days at various current rates. Additionally, DFT calculations confirmed that CCS/Ov catalysts demonstrate significantly enhanced HER and OER activities due to reduced H2O dissociation energy and strong intermediate binding energy, which is attributed to the anionic vacancy near Co3+. The excellent bifunctional performance of the hierarchical CCS/Ov highlights its potential as non-precious catalyst with facile fabrication approach.

A theoretical and kinetic study of key reactions between ammonia and fuel molecules, part III: H-atom abstraction from esters by ṄH2 radicals

Hydrogen atom abstraction reactions by ṄH2 radicals play a crucial role in determining the reactivity of ammonia/fuel binary blends. Esters are a typical component of environmentally friendly and economically promising biofuels. The feasibility of the ammonia/biofuel dual-fuel approach has been proven in practical engines. [Energy and Fuels 22 (2008) 2963] and [Int. J. Energy Res. 2023 (2023) 9920670]. ṄH2 radicals play a critical role in the combustion and pyrolysis chemistry of ammonia and N-containing-rich fuels. In ammonia/biofuels hybrid combustion, ṄH2 radicals can react with biofuel molecules in a reaction class that is particularly important especially when sufficient ammonia is blended in order to eliminate NOx emissions. To help unravel the chemistry of ammonia/biofuel blends, a systematic theoretical kinetic study of H-atom abstraction from eleven alky esters of CnH2n+1COOCH3 (n = 1–4), CH3COOCmH2m+1 (m = 1–4), and C2H5COOC2H5, by ṄH2 radicals is performed in this work. The geometry optimization, frequency, and zero-point energy calculations for all related species, as well as the hindrance potential energy surface for low frequency torsional modes in the reactants and transition states, were performed at the M06–2X/6–311++G(d,p) level of theory. Intrinsic reaction coordinate calculations were performed to validate the connections between the transition states and expected minima energy species. The energies of all of the species involved were calculated at the QCISD(T)/cc-pVXZ (X = D, T, Q) and MP2/cc-pVYZ (Y = T, Q) levels of theory and then extrapolated to the complete basis set. Rate constants of 39 reactions were calculated using the Master Equation System Solver (MESS) program in the temperature range of 500 – 2000 K. These rate constants for different H-atom abstraction sites are provided and can be extrapolated to larger esters. The kinetic effects from the functional group are also illustrated by performing detailed comparisons with the previous studies of ṄH2 radical reactions with alkanes, alcohols and ethers.

Accurately Predicting Barrier Heights for Radical Reactions in Solution Using Deep Graph Networks.

Quantitative estimates of reaction barriers and solvent effects are essential for developing kinetic mechanisms and predicting reaction outcomes. Here, we create a new data set of 5,600 unique elementary radical reactions calculated using the M06-2X/def2-QZVP//B3LYP-D3(BJ)/def2-TZVP level of theory. A conformer search is done for each species using TPSS/def2-TZVP. Gibbs free energies of activation and of reaction for these radical reactions in 40 common solvents are obtained using COSMO-RS for solvation effects. These balanced reactions involve the elements H, C, N, O, and S, contain up to 19 heavy atoms, and have atom-mapped SMILES. All transition states are verified by an intrinsic reaction coordinate calculation. We next train a deep graph network to directly estimate the Gibbs free energy of activation and of reaction in both gas and solution phases using only the atom-mapped SMILES of the reactant and product and the SMILES of the solvent. This simple input representation avoids computationally expensive optimizations for the reactant, transition state, and product structures during inference, making our model well-suited for high-throughput predictive chemistry and quickly providing information for (retro-)synthesis planning tools. To properly measure model performance, we report results on both interpolative and extrapolative data splits and also compare to several baseline models. During training and testing, the data set is augmented by including the reverse direction of each reaction and variants with different resonance structures. After data augmentation, we have around 2 million entries to train the model, which achieves a testing set mean absolute error of 1.16 kcal mol-1 for the Gibbs free energy of activation in solution. We anticipate this model will accelerate predictions for high-throughput screening to quickly identify relevant reactions in solution, and our data set will serve as a benchmark for future studies.

Built-in Electric Field in 1D/2D Heterostructure Boosts Zinc Air Battery Performance.

The realization of a rechargeable zinc-air battery (ZAB) is hindered by the low intrinsic oxygen evolution reaction (OER) and oxygen reduction reaction (ORR) activities. In this work, an abundant built-in electric field is noticed in a 1D/2D CoO/CoS2 heterostructure, triggering electron transfer from CoO to CoS2 associated with a downshifted d band center of the Co atom mitigating the strong electrochemical adsorption of *OH species on active sites; thereby, boosted OER and ORR performance are achieved. Namely, the OER specific activity of CoO/CoS2 is enhanced by 3.8- and 2.2-fold compared to the counterpart of CoO and CoS2, respectively. Furthermore, the kinetic current density of CoO/CoS2, a fingerprint of intrinsic ORR activity, is promoted by 46 and 6.6 times relative to CoO and CoS2. The rechargeable ZAB performance attains 215.6 mW cm-2, 1.6-times better than Pt/C-IrO2. Moreover, the superior performance remained for 600 h. Besides, the battery performance of the all-solid-state ZAB reaches 83.8 mW cm-2, revealing its promising application in wearable device.

A DFT Mechanistic Study on the Aza-Aldol Reaction of Boron Aza-Enolates: Relative Stability of Six-Membered Transition State and Its Relevance to the Coordination Mode of the Leaving Group.

The mechanism of the aza-aldol reaction between boron aza-enolate and benzaldehyde is investigated by using density functional theory calculations. The result shows that the syn-E isomer is preferentially formed, consistent with experimental observations. The six-membered ring transition state (TS) with the boat form leads to the E isomer, while the more unstable chair TS does to the Z isomer. The preference of the syn isomer is determined by the interactions between the substituents of aza-enolate and benzaldehyde. Structural distortion and intrinsic reaction coordinate analyses of simplified model systems provide insights into the origin of the relative stability of the rate-determining TS with boat and chair forms. The boat TS is an early TS; thus, minimal structural distortions of the reactant are required to reach this TS. The Lewis pair interactions between the boron and imine groups during B-N elimination also influenced the relative stability of the TSs. This interaction involves the nitrogen lone pair in the boat TS, while the π(N═C) orbital is involved in the chair TS. The Lewis pair with the lone pair stabilizes the TS more than that with the π orbital. The boron aza-enolate with 9-BBN generates an ate complex and forms C-C bonds sequentially, whereas that with Bpin does not generate an ate complex and exhibits the concerted formation of B-O and C-C bonds. Thus, the higher electrophilicity of boron such as 9-BBN enhances the reactivity by facilitating the formation of the ate complex. A reaction design is proposed to reverse the syn/anti selectivity. Proof-of-concept DFT calculations suggested that the modification of the imine group would change the relative stability of the boat/chair TSs and give the anti-product.

On the link between the reaction force constant and conceptual DFT.

The reaction force constant ( ), introduced by Professor Alejandro Toro-Labbé, plays a pivotal role in characterizing the reaction pathway by assessing the curvature of the potential energy profile along the intrinsic reaction coordinate. This study establishes a novel link between and the reactivity descriptors of conceptual density functional theory (c-DFT). Specifically, we derive expressions that relate the reaction force constant to nuclear softness and variations in chemical potential. Our findings indicate that regions of the reaction pathway where is negative match with significant electronic structure rearrangements, while positive regions match mostly with geometric rearrangements. This correlation between and c-DFT reactivity descriptors enhances our understanding of the underlying forces driving chemical reactions and offers new perspectives for analyzing reaction mechanisms. The internal reaction path for the proton transfer in SNOH, chemical potential, and nuclear softness were computed using DFT with B3LYP exchange-correlation functional and 6-311++G(d,2p) basis set.

Selective electro-capacitive deionization of Yb(III) by nanospherical N-doped carbon supported nickel‑iron bimetallic oxide, nitride and phosphide

The separation of rare earth elements from the mining waste water is highly important but still facing challenge. In this study, three different porous composite electrode materials were well-designed by decorating the nickel‑iron bimetallic oxide, nitride or phosphide onto nitrogen-doped carbon nanosphere, respectively (termed as NiFeO@NC, NiFeN@NC and NiFeP@NC), and utilized in capacitive deionization device for electrosorption of Yb(III) from water under the low-voltage electric field. The characterization and electrochemical results indicated that the incorporation of N or P atom lead to the reduced specific surface area and the enhanced electrochemical properties (i.e., larger specific capacitance, lower charge transfer resistance). In addition, NiFeN@NC and NiFeP@NC exhibited the excellent electrosorption capacities of Yb(III) with 105.24 mg·g−1 and 90.31 mg·g−1, respectively, under conditions of a flow rate of 20 mL·min−1, a voltage of 1.2 V, and pH of 5.0, which was extremely higher than NiFeO@NC (i.e., 50.7 mg·g−1). Moreover, all these three materials also demonstrated the remarkable electrosorption selectivity of Yb(III) from the mixed solution, and the robust durability with over 90 % of initial capacities after ten consecutive electrosorption-desorption cycles. Furthermore, the XPS characterizations and electrosorption results revealed that the electrosorption mechanism mainly depended on the synergistic effects of intrinsic Faraday redox reaction, electric double layer capacitance and active binding sites. Therefore, this work provided a feasible strategy for the separation of rare earth elements from aqueous solution, and the prospective applications related to the recovery of rare earth elements from mining wastes or spent permanent magnets in future life.

Catalytic and non-catalytic reactions of methane pyrolysis for hydrogen production in screening and electromagnetic levitation reactors using computational fluid dynamics

Molten metal (MM)-based methane (CH4) pyrolysis is a promising technology for clean hydrogen production with a minimal carbon footprint. Electromagnetic levitation (EMLR) and screening reactors (SR) were used to investigate the intrinsic catalytic reaction kinetics of MM-based CH4 pyrolysis. Heterogeneous catalytic reactions in the SR and EMLR occurred on the surfaces of the MM crucible and droplet, respectively. However, a homogenous non-catalytic reaction may occur just above the surface because the CH4 gas is preheated by the high-temperature MM. This study presents three-dimensional (3D) computational fluid dynamics (CFD) model coupled with chemical reactions to assess the extent to which non-catalytic reactions intervene in CH4 pyrolysis in both the SR and EMLR, which is difficult to measure experimentally. The CFD model validated against experimental data revealed that the non-catalytic reaction accounted for 4.0 % of CH4 pyrolysis in the CH4-Ni0.27Bi0.73 SR at 1045 °C, whereas it represented 0.02 % in the CH4-Sn EMLR at 1021 °C. This result highlights that the EMLR is suitable for studying the intrinsic catalytic reaction kinetics of MM-based CH4 pyrolysis because the non-catalytic reaction occurs to a lesser extent.

Kinetics of N2O decomposition over bulk and supported LaCoO3 perovskites

For the first time, kinetic data on the decomposition of N2O over mixed oxide LaCoO3 with a perovskite structure have been obtained. Bulk LaCoO3 synthesized using microwave activation exhibited an increased intrinsic reaction rate with a 30 kJ mol–1 lower activation energy. Perovskite samples supported on ZrO2-La demonstrated lower intrinsic reaction rates due to the lower content of the LaCoO3 phase, but the activation energies were also lower by 50–80 kJ mol–1.

Performance and mechanism investigation of Cu loading and S vacancy synergistic promotion of MoS2 catalyzed low-temperature water gas shift reaction: Experimental and theoretical calculations

In this paper, MoS2 with S vacancies and Cu loading (Cu/VS-MoS2) was synthesized and used for low-temperature water gas shift reactions (WGSR); the intrinsic reasons and reaction mechanisms for the high catalytic activity of Cu/VS-MoS2 were studied through DFT calculation. The material characterization fully confirms the successful construction of Cu loading and S vacancies on spherical cluster structures of Cu/VS-MoS2. The experimental results show that Cu/VS-MoS2 exhibits the highest CO conversion, up to 40.41%, which is 3.9 times, 1.8 times, and 1.5 times higher than MoS2, VS-MoS2, and Cu/MoS2, respectively. DFT calculations indicate that Cu/VS-MoS2 has high WGSR activity due to the superior electron transfer ability of the dual active sites, and the S vacancy is used to activate CO to generate CO2, while the Cu site activates H2O to produce hydrogen. The synergistic effect of S vacancies and Cu sites, while satisfying hydrogen production and carbon dioxide enrichment, greatly enhances the activity of MoS2. In addition, Cu/VS-MoS2 can more effectively break down the transfer barriers between surface electrons and CO, H2O molecules, allowing electrons to continuously transfer to CO and H2O using S vacancies and Cu atoms as mediators, reducing the energy barrier of the rate-determining step in association mechanism.

Topology Control of Covalent Organic Frameworks with Interlaced Unsaturated 2D and Saturated 3D Units for Boosting Electrocatalytic Hydrogen Peroxide Production.

Modulating the electronic state of multicomponent covalent organic framework (COF) electrocatalysts is crucial for enhancing catalytic activity. However, the effect of dimensionality on their physicochemical functionalities is still lacking. Herein, we report an interlaced unsaturated 2D and saturated 3D strategy to develop multicomponent-regulated COFs with tunable gradient dimensionality for high selectivity and activity electrocatalysis. Compared with the two-component 2D and 3D model COFs, the 2D/3D framework interlaced COFs with locally irregular dimensions and electronic structures are more practical in optimizing the intrinsic electrode surface reaction and mass transfer. Remarkably, the unsaturated 2D-inserted 3D TAE-COF regulates the adsorption mode of OOH* species to supply a favorable dynamic pathway for the H2O2 process, thereby achieving an excellent production rate of 8.50 mol gcat -1 h-1. Moreover, utilizing theoretical calculation and in situ ATR-FTIR experiment, we found that the central carbon atom of the tetraphenyl-based unit (site-1 and site-6) are potential active sites. This strategy of operating the adsorption ability of reactants with dimensionality-interconnected building blocks provides an idea for designing durable and efficient electrocatalysts.

Topology Control of Covalent Organic Frameworks with Interlaced Unsaturated 2D and Saturated 3D Units for Boosting Electrocatalytic Hydrogen Peroxide Production

AbstractModulating the electronic state of multicomponent covalent organic framework (COF) electrocatalysts is crucial for enhancing catalytic activity. However, the effect of dimensionality on their physicochemical functionalities is still lacking. Herein, we report an interlaced unsaturated 2D and saturated 3D strategy to develop multicomponent‐regulated COFs with tunable gradient dimensionality for high selectivity and activity electrocatalysis. Compared with the two‐component 2D and 3D model COFs, the 2D/3D framework interlaced COFs with locally irregular dimensions and electronic structures are more practical in optimizing the intrinsic electrode surface reaction and mass transfer. Remarkably, the unsaturated 2D‐inserted 3D TAE‐COF regulates the adsorption mode of OOH* species to supply a favorable dynamic pathway for the H2O2 process, thereby achieving an excellent production rate of 8.50 mol gcat−1 h−1. Moreover, utilizing theoretical calculation and in situ ATR‐FTIR experiment, we found that the central carbon atom of the tetraphenyl‐based unit (site‐1 and site‐6) are potential active sites. This strategy of operating the adsorption ability of reactants with dimensionality‐interconnected building blocks provides an idea for designing durable and efficient electrocatalysts.

Multi-Level Protocol for Mechanistic Reaction Studies Using Semi-Local Fitted Potential Energy Surfaces.

In this work, we propose a multi-level protocol for routine theoretical studies of chemical reaction mechanisms. The initial reaction paths of our investigated systems are sampled using the Nudged Elastic Band (NEB) method driven by a cheap electronic structure method. Forces recalculated at the more accurate electronic structure theory for a set of points on the path are fitted with a machine learning technique (in our case symmetric gradient domain machine learning or sGDML) to produce a semi-local reactive potential energy surface (PES), embracing reactants, products and transition state (TS) regions. This approach has been successfully applied to a unimolecular (Bergman cyclization of enediyne) and a bimolecular (SN2 substitution) reaction. In particular, we demonstrate that with only 50 to 150 energy-force evaluations with the accurate reference methods (here complete-active-space self-consistent field, CASSCF, and coupled-cluster singles and doubles, CCSD) it is possible to construct a semi-local PES giving qualitative agreement for stationary-point geometries, intrinsic reaction coordinates and barriers. Furthermore, we find a qualitative agreement in vibrational frequencies and reaction rate coefficients. The key aspect of the method's performance is its multi-level nature, which not only saves computational effort but also allows extracting meaningful information along the reaction path, characterized by zero gradients in all but one direction. Agnostic to the nature of the TS and computationally economic, the protocol can be readily automated and routinely used for mechanistic reaction studies.

An efficient technology for dye intermediate production: Process intensification of catalytic hydrogenation of aromatic nitrocompounds in a rotating packed bed reactor

AbstractThe rotating packed bed (RPB) is suitable for catalytic hydrogenation of aromatic nitrocompounds whose reaction rate is limited by the mass transfer rate. However, the matching mechanism between the mass transfer rate and intrinsic reaction rate in the RPB reactor is still unclear. In this work, the matching relationship in the RPB reactor is studied by combining the characteristics of mass transfer and reaction in detail. Effects of various operating conditions on the gas–liquid mass transfer coefficient up to 0.31/s, under working conditions were studied and the mass transfer correlation was proposed. The hydrogenation of p‐nitroanisole was used as a probe reaction to analyze the matching mechanism. Various aromatic nitrocompounds in the RPB and stirred tank reactors were investigated. The reaction efficiency of the former is 5–20 times that of the latter, which further explained the universality of the matching mechanism in the RPB reactor.

New polymorph of N-benzoyl-morpholine-4-carbothioamide with ten crystallographically independent molecules in the asymmetric unit: Crystal structure, Hirshfeld surface analysis and Density functional theory calculations

Herein, facile synthesis, single crystal X-ray diffraction (SC-XRD) structure, Hirshfeld surface (HS) analysis and Density functional theory (DFT) calculations of N-benzoyl-morpholine-4-carbothioamide 5(a) are presented. Hence, the synthesized compound, C12H14N2O2S 5(a), crystallizes in monoclinic crystal system having space group of P 21 with a = 22.2226 (9) Å, b = 11.7520 (3) Å, c = 24.9361 (8) Å, α = 90°, β = 111.054 (4)°, γ = 90°, Z = 20 and V = 6077.6 (4) Å3. The asymmetric unit of 5(a) contains ten crystallographically independent molecules, where the morpholine and benzene rings are in respective chair and planar conformations. In the crystal structure, some of the intermolecular, N—H···O, C—H···O and C—H···S, hydrogen bonds link the molecules into a 2D-network parallel to (010), enclosing R22(10) and R22(24) ring motifs, when viewed down the b-axis direction. Furthermore, contribution of the remaining hydrogen bonds consolidates a 3D-architecture. Similarly, HS analysis clarified the importance of hydrogen atom contacts with the major contributions of H…H (46.7 %), H…S/S…H (21.1 %), H…C/C…H (16.5 %), and H…O/O…H (12.5 %). Thus, hydrogen bond and van der Waals interactions play major role in the crystal packing. Finally, optimum molecular structure of compound 5(a) was related to the experimentally determined one employing DFT at B3LYP/6–311++G(d,p) level. Subsequently, the geometry of the optimized and most stable conformation turned out to be complementary with the crystal structure. The computational calculation of compound 5(a) revealed the molecule to be chemically hard in nature with high energy intrinsic reaction pathway. The associated band gap predicted compound 5(a) to be a wide band gap material. The intrinsic reaction pathway predicts hard-core behavior of the molecule.

Thermal decomposition study of 1,3-dioxolane as a representative of battery solvent

As industrial compounds and organic solvent, 1,3-dioxolane (C3H6O2) is used in lithium-ion batteries arousing extensive interest. The pyrolysis of C3H6O2 at 150 torr pressure was conducted in the temperature range of 743–1013 K. In order to quantify the major products, experimental measurements of the photoionization cross-section (PICS) and photoionization efficiency (PIE) for C3H6O2 were carried out by using synchrotron radiation photoionization mass spectrometry (SR-PIMS). Twenty-three intermediates were detected in pyrolysis of C3H6O2 including the newly observed ketene, acetaldehyde, ethenol, ethanol, dimethyl ether, diacetylene and 2-propen-1-ol which played important roles in developing mechanism. The H-abstraction of C3H6O2 by H, CH3 and OH radicals were studied via high level ab initio calculations. EStokTP code was used to find and calculate optimal geometries, vibration frequencies and intrinsic reaction coordinate (IRC) at the theory level of M06–2X/6–311+G(d,p). Conventional transition state theory (TST) was used and Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) was solved to obtain rate coefficients. For H-abstraction by H and CH3, the rates are increasing versus temperature, and rates of H-abstraction by OH show different tendency. These rate parameters were used to update the model of C3H6O2, which could well predict the ignition delay times (IDTs) in shock tube and oxidation data in jet stirred reactor (JSR) from literature. For the pyrolysis, H-abstraction reactions by H, CH3 and OH and unimolecular reactions contribute to solvent decay. Detailed pathways of C3H6O2 decomposition were also given and analyzed through this experiment and previous literature. The findings, including the experimental data and the calculations will assist to better understand pyrolysis characteristics of C3H6O2 as well.

Theoretical predictions on In-Plane aromaticity in double hydrogen transfer reactions between ethane and ethylene

Aromaticity, traditionally understood as a cornerstone concept elucidating the stability, reactivity, and structure of unsaturated organic compounds, has evolved to encompass inorganic molecules, saturated systems with delocalized electrons, and, significantly, transition structures. Our research delves into the exploration of in-plane aromaticity within the transition state (TS) of the Double Hydrogen Transfer reaction (DHTR) between ethylene and ethane, marking a novel approach to understanding cyclic reaction mechanisms. Through the application of density functional theory and verification via the Gaussian09 program’s intrinsic reaction coordinate analysis, our investigation assesses parameters such as aromatic nucleus-independent chemical shift (NICS) and Mayer bond orders to characterize the TS. The findings reveal a distinctly synchronous in-plane aromatic TS, highlighting its aromatic features. This research paves the way for future studies aimed at unravelling the influence of aromaticity on the kinetics and pathways of similar reactions, enhancing our grasp of its pivotal role in chemical processes.

Impact of mixing low-reactivity gases on the mechanism of hydrogen spontaneous combustion: A ReaxFF MD study

The phenomenon of spontaneous combustion during the release of high-pressure H2 is critical to the safe storage and transportation. Besides physical methods to suppress self-ignition, the approach of mixing low-reactivity gases has gained increasing attention. Previous work has mainly focused on the flow, physical heating effects, and boundary conditions of self-ignition in mixed gases, while the intrinsic chemical reaction mechanisms of spontaneous combustion remain underexplored. To clear this issue, this work uses the ReaxFF MD simulation method to study the spontaneous combustion process of H2/CO and H2/CH4 mixtures with various concentrations. The results indicate that: (1) in pure H2, H2/CO, and H2/CH4 systems, the primary initial reaction mechanism is H2+O2HO2+H. Additionally, in some mixed systems, CO + O2CO2+O and CH4+O2CH3+HO2 may also serve as initial reaction mechanisms. (2) Under fuel-lean combustion conditions, the RIDTs of the systems decrease, supporting the viewpoint that self-ignition typically occurs first in fuel-lean regions. The addition of CO and CH4 increases the RIDTs, which helps inhibit the development of self-ignition. (3) In H2-rich environments, CO combusts to CO2 through the formation of HOCO, while CH4 undergoes dehydrogenation and oxidation processes through intermediates such as CH3, CH3O, and HCHO, leading to the formation of CO and CO2. The combustion of H2 is also influenced by these processes. (4) The addition of CO and CH4 increases the activation energy of H2 combustion, raising the reaction energy barrier that must be overcome for spontaneous combustion. This effectively helps to suppress the occurrence and development of self-ignition during the release of high-pressure gases.

Density functional theory and ReaxFF MD study on steam-induced nitrogen migration mechanism during char gasification

The formation of nitrogen-containing species during char gasification is crucial for emission of nitrogenous pollutants. In this work, the detailed nitrogen migration mechanism during char gasification is studied. Density functional theory (DFT) coupled with ReaxFF MD methods are used to conduct an in-depth analysis on the intrinsic reaction mechanism during Char(N) and steam interaction. DFT results show that orbital electrons of carbon atoms on char surface are driven towards nitrogen atoms by nitrogen functional groups. The electron density of adjacent carbon atoms is weakened, which has a positive charge. Orbital electron properties show that the band energy gaps of three Char(N) models are 3.063 eV, 1.092 eV and 3.328 eV, indicating the reactivity order for three Char(N) models is as follows: Char(N)-2＞Char(N)-1＞Char(N)-3. DFT calculation indicates that the interaction between Char(N) and steam reduces unsaturated carbon atoms and lower the char decomposition activity, which will inhibit the yield of HCN. In contrast, through hydrogen transfer reactions, nitrogen atoms are easily combined with hydrogen atoms, which is more conducive to the formation of ammonia. ReaxFF MD modeling proves that HCN is the main product for Char(N) decomposition under inert atmosphere. Under steam atmosphere, nitrogen atoms in char are more likely to be converted into amino products. The induction of steam provides a large number of active hydrogen radicals, which is able to attack the nitrogen atoms of char and form N–H bonds, thus promoting the formation of ammonia.