Twenty-species and 15-species chemical kinetic mechanisms for cyclohexane using the local self-similarity tabulation method

The 1081 species cyclohexane-oxidation elementary reaction mechanism of Silke et al. (DOI: 10.1021/jp067592d ) is reduced in the number of species by a factor using the local self-similarity tabulation (LS2T) method. Reduced-species mechanisms of both 20 (R20) and 15 (R15) species are created in the high-pressure combustion regime typical of diesel engines. To evaluate the performance of R20 and R15 against the elementary kinetics, simulations are performed for cyclohexane/air mixtures at initial temperatures of 1150, 900, 750, and 680 K and constant pressures of 20 and 40 bar for a variety of equivalence ratios (, 1.0, and 2 for 1150 and 900 K; for 750 K; for 680 K). Very good agreement between R20 and R15 with the elementary kinetics mechanism is demonstrated at 1150 and 900 K for which the self-similarity is very well obeyed; however, only fair agreement is obtained at 750 and 680 K, a fact which is traced to the less faithful adherence to the self-similarity due to the one order of magnitude increase in ignition time over the range 750-680 K. These results are found to be quasi-independent of the tabulation grid. Future work is proposed to improve the reduction in the cold-ignition, high-pressure regime.     

Oxidation of methylamine

A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. The rate coefficients for the reactions of CHNH + O CHNH / CHNH + HO, CHNH + H CH + NH, CHNH CHNH, and CHNH + O CHNH + HO were calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CHNH and its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NH group in CHNH compared to CH. More work is desirable to confirm the products of the CHNH + NO reaction and to look into possible pathways to NH in methylamine oxidation.     

Thermal rate constant for the C(3P) + OH(X2Π) → CO(X1Σ) + H(2S) reaction using stochastic energy grained master equation method

In the present work, the kinetic mechanism of the reaction is studied. The rate constants were determined using the Master Equation Solver for Multi-Energy Well Reactions (MESMER). The master equation modeling was also employed to examine the pressure dependence for each pathway involved. The theoretical analysis shows that the overall rate coefficient is practically independent of pressure up to 100 Torr for the temperature range 125-500 K. The unusual dependence of the overall rate constant with temperature was fit with the d-Arrhenius expression , where cm³molecule⁻¹s⁻¹, , and  kJ·mol⁻¹, for 125⩽ T ⩽ 500 K. The thermal rate constant results are in relatively good agreement with other theoretical studies.     

Finding of remarkable synergistic effect on the aroxyl-radical-scavenging rates under the coexistence of α-tocopherol and catechins

Measurements of aroxyl radical (ArO^•)-scavenging rate constants () of antioxidants (AOHs) (α-tocopherol (α-TocH) and three catechins (CatHs) (ie, epicatechin (EC), epigallocatechin (EGC), and epigallocatechin gallate (EGCG)) were performed in ethanol solution, using stopped-flow spectrophotometry. values were measured not only for each AOH, but also for the mixtures of two AOHs (α-TocH and CatH). A notable synergistic effect that the value of α-TocH increases 1.29, 1.84, and 1.65 times under the coexistence of constant concentrations of EC, EGC, and EGCG, respectively, was observed for the solutions including α-TocH and CatH. Similarly, values of CatHs (EC, EGC, and EGCG) increased 1.72, 2.25, and 2.34 times under the coexistence of constant concentrations of α-TocH, respectively. UV-Vis absorption of α-tocopheroxyl radical (α-Toc^•) (λ_max = 428 nm), which had been produced by reaction of α-TocH with ArO^•, decreased remarkably under the coexistence of α-TocH and CatHs due to the fast α-TocH-regeneration reaction by CatHs. The result suggests that the prooxidant reaction due to α-Toc^• is suppressed by the coexistence of CatHs. By analyzing the formation and decay curves of α-Toc^•, it has been ascertained that one molecule of EGCG having three OH groups at B-ring may rapidly regenerate three molecules of α-Toc^• to α-TocH.     

Kinetics of two redox reactions in 2-butoxyethanol and water near its critical point of solution

The rate constants of two redox reactions and in the critical solution of 2-butoxyethanol and water have been measured by using the UV spectrophotometry at the initial reaction stage. It was found that the rate constants at various temperatures for two reactions were well described by the Arrhenius equation in the noncritical region. The critical slowing down effect was detected in the critical region. The critical slowing down exponents were determined to be 0.044 ± 0.004 and 0.046 ± 0.005 for reactions and , respectively. The values of the critical slowing down exponents showed that only dynamic critical slowing down effect, and no thermodynamic singularity could be observed for the two reactions.     

Multiconformer transition state theory rate constants for the reaction between OH and -dimethoxyfluoropolyethers

In this work, we have calculated rate constants for the tropospheric reaction between the OH radical and -dimethoxyfluoropolyethers. The latter are a specific class of the hydrofluoropolyethers family with the general formula , from which we have selected three case studies: , , and . The calculations were performed by applying a cost-effective protocol developed for bimolecular hydrogen-abstraction reactions and based on multiconformer transition state theory relying on computationally accessible M08-HX/apcseg-2//M08-HX/pcseg-1 calculations. Within the protocol's uncertainties and approximations, the results show that (1) the calculated rate constants have the same order of magnitude and (2) if observed together with previous experimental and theoretical investigations, the chain length (that varies with q and p) is seen to have a small effect on the rate constant, which is consistent with the “no discernible effect” reported in the experimental work.     

Kinetic model for Ostwald's rule of stages with applications to Boc-diphenylalanine self-assembly

A kinetic model which describes Ostwald's rule of stages, during the process of crystal growth from solution, is reported here. Reaction equations for stages are given where the stages convert from one to another. The final stage reacts to release a portion of solute back into solution, while the remainder converts to the final equilibrium form. Additionally, a remnant of the solute that was not consumed by any of the transitional stages, ultimately is converted into the final product. This particular model was motivated by a recent report for Boc-diphenylalanine self-assembly where the dissolved peptide was observed to go through two polymorphic stages before reaching the equilibrium supramolecular assembly [A. Levin et al., Nat. Commun. 5, 5219, (2014)]. Kinetic data for the concentration of solute present during the process are listed in the above-mentioned report. We show here how the model, for , describes the time-dependent behavior of the solute decay during the growth process. After comparing the model to the experimental data, we are able to report values for all of the rate constants and propose a rule whereby the relative magnitudes of these constants can be used to predict whether a supersaturated substance will noticeably pass through transitional stages or simply convert from solute to the equilibrium solid form.     

Predicting polycyclic aromatic hydrocarbon formation with an automatically generated mechanism for acetylene pyrolysis

Using Reaction Mechanism Generator (RMG), we have automatically constructed a detailed mechanism for acetylene pyrolysis, which predicts formation of polycyclic aromatic hydrocarbons (PAHs) up to pyrene. To improve the data available for formation pathways from naphthalene to pyrene, new high‐pressure limit reaction rate coefficients and species thermochemistry were calculated using a combination of electronic structure data from the literature and new quantum calculations. Pressure‐dependent kinetics for the CH potential energy surface calculated by Zádor et al. were incorporated to ensure accurate pathways for acetylene initiation reactions. After adding these new data into the RMG database, a pressure‐dependent mechanism was generated in a single RMG simulation which captures chemistry from C to C. In general, the RMG‐generated model accurately predicts major species profiles in comparison to plug‐flow reactor data from the literature. The primary shortcoming of the model is that formation of anthracene, phenanthrene, and pyrene are underpredicted, and PAHs beyond pyrene are not captured. Reaction path analysis was performed for the RMG model to identify key pathways. Notable conclusions include the importance of accounting for the acetone impurity in acetylene in accurately predicting formation of odd‐carbon species, the remarkably low contribution of acetylene dimerization to vinylacetylene or diacetylene, and the dominance of the hydrogen abstraction CH addition (HACA) mechanism in the formation pathways to all PAH species in the model. This work demonstrates the improved ability of RMG to model PAH formation, while highlighting the need for more kinetics data for elementary reaction pathways to larger PAHs.     

Kinetics of wet air oxidation of sodium sulfide over heterogeneous iron catalyst

Wet air oxidation (WAO) is an established technique for reducing the chemical oxygen demand (COD) of refinery sulfidic spent caustic waste. In the present work, the heterogeneous form of the cheap and abundant catalyst ferrous sulfate (FeSO₄) was employed for WAO of sodium sulfide. The performance of this catalyst in the oxidative destruction of this model compound is thus far unfamiliar. Kinetic data for the non-catalytic and catalytic oxidation processes was collected in a batch reactor. For the catalytic process, temperature (T), oxygen partial pressure () and catalyst concentration (ω) were varied in the ranges 80-150°C, 0.69-2.06 MPa and 0.8-2.4 g/L respectively. Around 94% COD was destroyed within 1 h when feed containing 8 g/L of sulfide was oxidized at T = 100°C, = 0.69 MPa, and ω = 0.8 g/L. First, the data on disappearance of COD were fitted to a power law model and reaction rate constants were determined. The activation energy for the non-catalytic (91 kJ/mol) and catalytic (50 kJ/mol) oxidation process was found from the temperature dependence of the rate constants. Second, hyperbolic models based on Langmuir-Hinshelwood (L-H) and Eley-Rideal (E-L) kinetics were used for fitting kinetic data. It was found that the L-H model suggesting dissociative adsorption of oxygen provided the best fit. In this way, a deep insight into oxidation kinetics of sodium sulfide was provided.     

Formation of phenanthrene via H-assisted isomerization of 2-ethynylbiphenyl produced in the reaction of phenyl with phenylacetylene

Model chemistry G3(MP2,CC)//B3LYP/6-311G(d,p) calculations of the potential energy surface for the reaction of phenyl radical (C₆H₅) with phenylacetylene (C₈H₆) have been carried out and combined with Rice-Ramsperger-Kassel-Marcus/Master Equation calculations of temperature- and pressure-dependent rate constants. The results showed that the reaction can serve as a viable source for the formation of phenanthrene via an indirect route involving a primary reaction of phenyl addition to the ortho carbon in the ring of phenylacetylene and H elimination producing 2-ethynylbiphenyl followed by secondary H-assisted isomerization of 2-ethynylbiphenyl to phenanthrene. In the secondary reaction, the H atom adds to the α carbon of the ethynyl side chain, then a six-member ring closure takes place followed by aromatization via an H loss. The channel of H addition to the side chain of 2-ethynylbiphenyl appears to be much faster than H addition to the ortho carbon in the ethynyl-substituted ring leading back to the initial C₆H₅ + C₈H₆ reactants. Rate constants for the primary C₆H₅ + C₈H₆2-ethynylbiphenyl ( p1 ) + H and secondary p1  + Hphenanthrene ( p2 ) + H reactions have been computed in the temperature range of 500-2500 K at pressures of 30 Torr, 1, 10, and 100 atm and fitted to modified Arrhenius expressions. The suggested kinetic scheme and rate constants are proposed as a prototype for the modeling of the growth of polycyclic aromatic hydrocarbons via the phenyl addition-dehydrocyclization (PAC) mechanism involving an addition of a PAH radical to an ethynyl-substituted PAH molecule.     

Resonant Fragmentation of the Water Cation by Electron Impact: a Wave-Packet Study

We have investigated the dissociation of a resonant state that can be formed in low energy electron scattering from H₂O⁺. We have chosen the second triplet resonance above the state of H₂O⁺ whose autoionization mainly produces H₂O⁺ ( ). We have considered both dissociation of the resonant state itself, dissociative recombination (DR), or the dissociation of the H₂O⁺ cation after autodetachment, dissociative excitation (DE). The time-evolution of a wave packet on the potential energy surfaces of the resonance and cationic states shows, for the initial conditions studied, that the probability for DR is about 38 % while the probability for DE is negligible.     

Fluxional bonds in quasi-planar and half-sandwich (M = K,Rb, and Cs)

Detailed molecular orbital and bonding analyses reveal the existence of both fluxional σ- and π-bonds in the global minima C_s ( 1 ) and C_s MB₁₈ ( 3 ) and transition states C_s ( 2 ) and C_s ( 4 ) of dianion and monoanions (M = K, Rb, and Cs). It is the fluxional bonds that facilitate the fluxional behaviors of the quasi-planar and half-sandwich which possess energy barriers smaller than the difference of the corresponding zero-point corrections. © 2019 Wiley Periodicals, Inc.     

On the normalized Laplacian of Möbius phenylene chain and its applications

Let denote a molecular graph of linear [n] phenylene with n hexagons and n squares, and let the Möbius phenylene chain be the graph obtained from the by identifying the opposite lateral edges in reversed way. Utilizing the decomposition theorem of the normalized Laplacian characteristic polynomial, we study the normalized Laplacian spectrum of , which consists of the eigenvalues of two symmetric matrices ℒ _R and ℒ _Q of order 3n. By investigating the relationship between the roots and coefficients of the characteristic polynomials of the two matrices above, we obtain an explicit closed-form formula of the multiplicative degree-Kirchhoff index as well as the number of spanning trees of . Furthermore, we determine the limited value for the quotient of the multiplicative degree-Kirchhoff index and the Gutman index of .     

Cover Feature: How Different are the Diamagnetic and Paramagnetic Contributions to Off-Nucleus Shielding in Aromatic and Antiaromatic Rings? (ChemPhysChem 9/2023)

The spatial variations in the diamagnetic and paramagnetic contributions to the off-nucleus isotropic shielding, , and to the zz component of the off-nucleus shielding tensor, , around benzene (C₆H₆) and cyclobutadiene (C₄H₄) are investigated using complete-active-space self-consistent field wavefunctions. Despite the substantial differences between and around the aromatic C₆H₆ and the antiaromatic C₄H₄, the diamagnetic and paramagnetic contributions to these quantities, and , and and , are found to behave similarly in the two molecules, shielding and deshielding, respectively, each ring and its surroundings. The different signs of the most popular aromaticity criterion, the nucleus-independent chemical shift (NICS), in C₆H₆ and C₄H₄ are shown to follow from a change in the balance between the respective diamagnetic and paramagnetic contributions. Thus, the different NICS values for antiaromatic and antiaromatic molecules cannot be attributed to differences in the ease of access to excited states only; differences in the electron density, which determines the overall bonding picture, also play an important role.     

Real-Space Interpretation of Interatomic Charge Transfer and Electron Exchange Effects by Combining Static and Kinetic Potentials and Associated Vector Fields**

Intricate behaviour of one-electron potentials from the Euler equation for electron density and corresponding gradient force fields in crystals was studied. Channels of locally enhanced kinetic potential and corresponding saddle Lagrange points were found between chemically bonded atoms. Superposition of electrostatic and kinetic potentials and electron density allowed partitioning any molecules and crystals into atomic - and potential-based -basins; -basins explicitly account for the electron exchange effect, which is missed for -ones. Phenomena of interatomic charge transfer and related electron exchange were explained in terms of space gaps between zero-flux surfaces of - and -basins. The gap between - and -basins represents the charge transfer, while the gap between - and -basins is a real-space manifestation of sharing the transferred electrons caused by the static exchange and kinetic effects as a response against the electron transfer. The regularity describing relative positions of -, -, and - basin boundaries between interacting atoms was proposed. The position of -boundary between - and -ones within an electron occupier atom determines the extent of transferred electron sharing. The stronger an H⋅⋅⋅O hydrogen bond is, the deeper hydrogen atom's -basin penetrates oxygen atom's -basin, while for covalent bonds a -boundary closely approaches a -one indicating almost complete sharing of the transferred electrons. In the case of ionic bonds, the same region corresponds to electron pairing within the -basin of an electron occupier atom.     

On the Mechanical Properties of Popgraphene-Based Nanotubes: a Reactive Molecular Dynamics Study

Carbon-based tubular materials have sparked a great interest in future electronics and optoelectronics device applications. In this work, we computationally studied the mechanical properties of nanotubes generated from popgraphene (PopNTs). Popgraphene is a 2D carbon allotrope composed of 5-8-5 rings. We carried out fully atomistic reactive (ReaxFF) molecular dynamics for PopNTs of different chiralities ( and ) and/or diameters and at different temperatures (from 300 up to 1200 K). Results showed that the tubes are thermally stable (at least up to 1200 K). All tubes presented stress/strain curves with a quasi-linear behavior followed by an abrupt drop of stress values. Interestingly, armchair-like PopNTs ( ) can stand a higher strain load before fracturing when contrasted to the zigzag-like ones ( ). Moreover, it was obtained that Young's modulus (Y_Mod) (750–900 GPa) and ultimate strength (σ_US) (120–150 GPa) values are similar to the ones reported for conventional armchair and zigzag carbon nanotubes. Y_Mod values obtained for PopNTs are not significantly temperature-dependent. While the σ_US values for the showed a quasi-linear dependence with the temperature, the exhibited no clear trends.     

Theoretical investigation of the formic acid decomposition kinetics

Decomposition of formic acid (HCO₂H) proceeds via three unimolecular channels: dehydration, decarboxylation, and dissociation, the latter expected to be of minor contribution to the overall kinetics. In addition, despite the similar values reported for the individual activation energies for the dehydration and decarboxylation reactions, experimental works have shown that the former is dominant in the reaction mechanism. These reactions show pressure-dependent rate coefficients, and the high-pressure condition is not yet verified at atmospheric pressure. This work aims to investigate the influence of temperature and pressure on the rate coefficients. Hence, theoretical calculations at the CCSD(T)/CBS level have been performed to accurately describe the unimolecular reaction and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficients have been calculated and integrated for the prediction of k(T,P) rate coefficients, adopting both strong and weak collision models, over the intervals 0.5-10 atm and 298-2200 K. Our results suggest that the isomerization path is important and explains the preference for the (CO + H₂O) channel. Rate coefficients for the (CO₂ + H₂) and (CO + H₂O) formations are given, in s⁻¹, as exp(−34404/T) and exp(−33785/T), respectively. The dissociation limit of 107.29 kcal mol^–1, with respect the Z-HCO₂H conformer, leading to OH + HCO, via a barrierless potential curve, with rate coefficients, in s⁻¹, expressed as k_HCO+OH(T) = 1.68 × 10¹⁷ exp(−56018/T). Temperature and pressure dependence for the HCO + OH → CO₂ + H₂ and HCO + OH → CO + H₂O reactions have also been estimated.     

The Threshold Photoelectron Spectrum of Fulvenone: A Reactive Ketene Derivative in Lignin Valorization

Unveiling reaction mechanisms by isomer-selective detection of reactive intermediates requires advanced spectroscopic knowledge. We study the photoionization of fulvenone (c-C₅H₄=C=O), a reactive ketene species relevant in catalytic pyrolysis of lignin, which was generated by pyrolysis of 2-methoxy acetophenone. The high-resolution threshold photoelectron spectrum (TPES) with vacuum ultraviolet synchrotron radiation revealed well-resolved vibrational transitions, assigned to ring deformation modes of the cyclopentadiene moiety. The adiabatic ionization energy was determined to be 8.25±0.01 eV and is assigned to the ²A₂← ¹A₁ transition. A broad and featureless band arising at 9 eV is associated with the ²B₁← ¹A₁ excitation. A conical intersection is responsible for the ultrafast relaxation of the fulvenone cation from the into the state resulting in a featureless and lifetime broadened band. These insights will increase the detection capabilities for fulvenone and thereby help to elucidate reaction mechanisms in lignin catalytic pyrolysis.     

Fluorine Substitution of TCNQ Alters the Redox-Driven Catalytic Pathway for the Ferricyanide-Thiosulfate Reaction

Mechanistic variation in catalysis through substituent-based redox tuning is well established. Fluorination of TCNQ (TCNQ=tetracyanoquinodimethane) provides ~850 mV variation in the redox potentials of the and (n=0, 2, 4) processes. With , catalysis of the kinetically very slow ferrocyanide-thiosulfate redox reaction in aqueous solution occurs via a mechanism in which the catalyst is reduced to when reacting with which is oxidised to . Subsequently, reacts with to form and reform the catalyst, in another thermodynamically favoured process. An analogous mechanism applies with as a catalyst. In contrast, since the reaction of with is thermodynamically unfavourable, an alternative mechanism is required to explain the catalytic activity observed in this non-fluorinated system. Here, upon addition of , reduction of to occurs with concomitant oxidation of to , which then acts as the catalyst for oxidation. Thermodynamic data explain the observed differences in the catalytic mechanisms. (n=0, 4) also act as catalysts for the ferricyanide-thiosulfate reaction in aqueous solution. The present study shows that homogeneous pathways are available following addition of these dissolved materials. Previously, these (n=0, 4) coordination polymers have been regarded as insoluble in water and proposed as heterogeneous catalysts for the ferricyanide-thiosulfate reaction. Details and mechanistic differences were established using UV-visible spectrophotometry and cyclic voltammetry.     

Hydrosilane σ-Adduct Intermediates in an Adaptive Zinc-Catalyzed Cross-dehydrocoupling of Si−H and O−H Bonds

Three-coordinate ^PhBOX ZnR (^PhBOX =phenyl-(4,4-dimethyl-oxazolinato; R=Me: 2 a , Et: 2 b ) catalyzes the dehydrocoupling of primary or secondary silanes and alcohols to give silyl ethers and hydrogen, with high turnover numbers (TON; up to 10⁷) under solvent-free conditions. Primary and secondary silanes react with small, medium, and large alcohols to give various degrees of substitution, from mono- to tri-alkoxylation, whereas tri-substituted silanes do not react with MeOH under these conditions. The effect of coordinative unsaturation on the behavior of the Zn catalyst is revealed through a dramatic variation of both rate law and experimental rate constants, which depend on the concentrations of both the alcohol and hydrosilane reactants. That is, the catalyst adapts its mechanism to access the most facile and efficient conversion. In particular, either alcohol or hydrosilane binds to the open coordination site on the ^PhBOX ZnOR catalyst to form a ^PhBOX ZnOR(HOR) complex under one set of conditions or an unprecedented σ-adduct ^PhBOX ZnOR(H−SiR′₃) under other conditions. Saturation kinetics provide evidence for the latter species, in support of the hypothesis that σ-bond metathesis reactions involving four-centered electrocyclic 2σ–2σ transition states are preceded by σ-adducts.