Quantum Mechanical Calculations Research Articles

Highly efficient upconversion photodynamic performance of rare-earth-coupled dual-photosensitizers: ultrafast experiments and excited-state calculations

Abstract Upconversion photodynamic therapy (UC-PDT), which integrates upconversion nanoparticles (UCNPs) with photosensitizers (PSs), presents a promising advancement in the field of phototherapy. However, despite the extensive studies focused on the design and synthesis of UCNPs, there is a paucity of systematic research on the mechanisms underlying the synergistic upconversion photodynamic effects. Here we have synthesized upconversion core@dotted-shell nanoparticles (CDSNPs) and covalently tethered them with two distinct PSs, thereby constructing a dual-PS UC-PDT system with high synergistic photodynamic performance. To unravel the mechanism underlying the synergism, we employed a combination of quantum mechanical calculations and ultrafast time-resolved spectroscopy techniques. The results indicate that rare earth oxides play a pivotal role in enhancing the intersystem crossing processes of PSs through modulating their excited electronic states. Additionally, Förster resonance energy transfer between two distinct PSs contributes to the amplification of triplet state populations, thus further enhancing the photodynamic effect. In vitro experiments demonstrate that the prepared CDSNPs based dual-PS system exhibits excellent biocompatibility with normal cells and exceptional synergistic photodynamic efficacy against tumor cells upon near-infrared excitation. This research contributes theoretical insights into the design and application of multi-photosensitizer UC-PDT systems, laying the groundwork for more efficient preclinical implementations in the future.

DESIGN AND IN SILICO VALIDATION OF VITAMIN B5 AND ITS DERIVATIVES AS A POTENTIAL TARGET AGAINST CYCLOPHILIN A (CyPA)

The title molecule (N-(2,4-dihydroxy-3,3-dimethylbutyryl)b-alanine), DDBBA and its derivatives were selected for theoretical investigations viz geometry optimization, ADME profiling, binding affinity using quantum-mechanical calculations and modelling simulation tools. Geometry optimization by Gaussian 09 program revealed the stability and electrophilic nature of the investigated molecules. In order to depict the charge density distributions that may be related to biological activity, the contour maps of HOMO-LUMO as well as the associated chemical descriptors such as chemical potential (µ), electronegativity (χ), electrophilicity (ω), hardness (η) and softness (σ) were explored. Utilizing molecular docking, the antiviral, antibacterial, and anticancer activities was examined. The docked molecules showed strong propensity for binding to 2HQ6 cancer protein active sites. Vit B5-CH=CF2 and VitB5-CCl3 showed lowest binding energies (-5.861 and -5.478 kcal/mol) and low inhibition constant values (1.43 M). Studies on the (NBO) natural bond orbital, the Mulliken population, and the Fukui function were all analyzed. Further, the interactions between the derivatives and other molecules were studied using Hirshfeld surface analysis.

C–H bond dissociation enthalpy prediction with machine learning reinforced semi-empirical quantum mechanical calculations

Abstract We introduce a combined fast semi-empirical quantum mechanical and machine learning (SQM/ML) approach capable of matching the C–H bond dissociation enthalpies (BDEs) computed with the highly accurate (RO)CBS-QB3 method. The usefulness of our proposed SQM/ML model is corroborated by the fact that a single C–H BDE of a molecule is calculated in seconds and the mean absolute error amounts to only 1 to 2 kcal/mol.

Chlorine degradation of semi-aromatic polypiperazine-amide membranes and the mechanisms

Semi-aromatic polypiperazine-amide (PPA) nanofiltration (NF) membranes are widely used in water treatment applications, but the resistance to free chlorine and other biocidal agents need to be improved. Understanding the chlorination mechanisms of PPA NF membranes is of great significance for avoiding membrane degradation and developing new chlorine-resistant membranes. In this paper, the effects of chlorination on the rejection selectivity as well as the surface charge, hydrophilicity and roughness, and pore size distribution of PPA NF membranes are presented. In addition, the chlorination mechanisms of the membranes are proposed through functional groups characterization of the stripped PPA active layer, element composition analysis, quantum mechanical calculations based on density-functional theory, and use of model compound ((3,5-dimethylphenyl)(4-methylpiperazin-1-yl)methanone) for PPA material. The molecular active chlorine species, HOCl, predominantly induces the cleavage of C–N bonds, whereas the primary reaction pathway caused by the ionic active chlorine species, ClO−, involves the ring-opening of benzene. Under the same chlorination conditions, HClO has a more severe impact on membrane properties and performance than ClO−. NF membranes with a looser polyamide layer are more susceptible to the attack by active chlorine species.

Efficient enzymatic synthesis of theaflavin and its production mechanism.

In this study, a novel preparation method of theaflavin (TF) has been established. Our findings indicated that the formation of TF was significantly enhanced by using an ice bath (2-3°C). Additionally, increasing the ratio of (-)-epigallocatechin (EGC) under the ice bath could further improve its yield. This approach prevented the appearance of a dark solution within 3h, effectively protecting TF from oxidation. Our study on the generation mechanism of TF suggested that EGC-quinone I (EGC-Q-I) with two carbanions could potentially serve as one of synthons based on the retrosynthetic analysis of the bicyclo[3.2.1]octane-type intermediate. Subsequently, quantum mechanical calculations further supported this hypothesis. Practical Application: In this study, we have developed a novel method for the synthesis of theaflavin (TF), demonstrating that the use of ice bath significantly enhanced its yield. Increasing the ratio of (-)-epigallocatechin (EGC) under the ice bath further improved TF yields and prevented darkening of the solution for at least 3h, thereby protecting TF from oxidation. Our study suggested that EGC-quinone I is a potential synthon based on the retrosynthetic analysis of the bicyclo[3.2.1]octane-type intermediate (BOI). This hypothesis is supported by QM calculations.

Elastic properties and constitutive behaviour of graphene at finite temperature and large deformation

In the earlier reported theoretical studies, the bending modulus of the graphene sheet predicted with the inclusion of the dihedral energy term is observed to be almost double as compared to those reported without considering the dihedral energy term. The bending modulus reported considering dihedral energy term is observed to be more accurate when compared with the quantum mechanical calculations. The findings of these studies for bending modulus (with/without considering dihedral energy terms) are reported without considering the effect of finite temperature. In the present work, a multiscale computational framework incorporating the dihedral energy term and finite temperature in the constitutive model is presented for the first time. A temperature-dependent bond length is expressed in terms of strain and curvature tensor through the temperature-related quadratic-type Cauchy-Born rule. The total Helmholtz free energy of the unit cell is expressed as a sum of the potential energy of the unit cell and thermal energy due to the thermal vibrations of atoms arising from the finite temperature. The atomic interactions are modelled through the Tersoff-Brenner potential considering different sets of empirical parameters, including the authors’ recently proposed set of empirical parameters. The stress/moment tensors and the tangent constitutive matrix are obtained from the derivatives of the Helmholtz free energy density with strains and curvatures. The thermal properties, elastic properties, and tangent stiffness coefficients for armchair and zigzag graphene with/without dihedral energy term at finite strain, curvature, and temperature are investigated and discussed in detail. The proposed multiscale model facilitates the accurate prediction of the constitutive behaviour of the graphene in a thermal environment with remarkable improvements for the prediction of bending modulus due to the inclusion of the dihedral energy term at finite temperature.

A systematic study on interplay between intermolecular hydrogen bonding and aromaticity

In this study, the interplay between aromaticity and intermolecular hydrogen bonding was investigated in some heterocyclic compounds using quantum mechanical calculations. The NH … O intermolecular hydrogen bonds between formamide and NH and CO functional groups of the selected compounds enhance the delocalisation of the six-membered ring electron cloud, resulting in an increase in the aromaticity of the ring. The contribution of various factors to the binding energy, including the stabilisation energy of the interplay between aromaticity and H-bonding, was estimated using two model systems. The relative binding energies and the estimated stabilisation energies of the aforementioned interplay can be explained by the weights of resonance structures calculated using the natural resonance theory (NRT) analysis with the natural bond orbitals (NBO) method. These weights reflect the tendency of ring nitrogen lone pairs to contribute to the π electron cloud of the ring. The results of the topological analysis of electron charge densities using the atoms in molecules (AIM) method were used to estimate the aromaticity of the ring. Aromaticity indices, particularly nucleus-independent chemical shifts (NICS) calculated one angstrom above the ring plane, are in good linear correlation with the estimated stabilisation energies and electron density at the H-bond critical points. The sum of the total nuclear spin–spin coupling 1J of adjacent atoms of the ring (Σ1J), and the J values of H-bond donor and H-bond acceptor atoms, 1HJNO, which are in good agreement with the energy data and aromaticity indices, can be used as measures of the interplay between aromaticity and H-bonding. The interplay between H-bonding and the aromaticity of the ring is clearly evident from the changes in the pattern of the molecular electrostatic potential (MEP) contour maps at 1.7 Å above the plane of the molecule upon H-bond complexation.

Mass spectrometry-guided isolation of thiodiketopiperazines from an EtOAc-extract of Setosphaeria rostrata culture medium and their anti-skin aging effects on TNF-α-induced human dermal fibroblasts.

Using mass spectrometry (MS)-guided isolation methods, a new thiodiketopiperazine derivative (1) and exserohilone (2) were isolated from an EtOAc-extract of Setosphaeria rostrata culture medium. The chemical structure of the new compound was elucidated by MS and NMR spectroscopy, and the absolute configurations were established by the quantum mechanical calculations of electronic circular dichroism. All isolated compounds were examined for their effects on reactive oxygen species (ROS) production, matrix metalloproteinase 1 (MMP-1) secretion, and procollagen type I α1 secretion in tumor necrosis factor (TNF)-α-induced human dermal fibroblasts. Compound 1 and exserohilone (2) exhibited the inhibition of TNF-α-induced ROS generation and MMP-1 secretion. Additionally, compound 1 and exserohilone (2) increased the procollagen type I α1 secretion. Compound 1 docked computationally into the active site of MMP-1 (-6.0 kcal/mol).

Unveiling Stereo‐Electronic Effects in Homogeneous Catalysis Integrating Theory and Experiments: the Potential of Dimeric Iron(III) Salen Complexes in Methyl Levulinate Transesterification

AbstractBiomass represents an abundant and renewable source of raw materials that can be effectively converted into intermediates and base chemicals. Among these, levulinic acid has garnered significant attention due to its direct derivation from lignocellulosic biomass and the presence of two functional groups, a carboxylic acid, and a ketone, enabling diverse transformations into valuable products. Levulinic acid esters have already found applications in various industrial sectors such as cosmetics, plasticizers, and solvents. In this study, a broad array of alkyl levulinates has been successfully synthesized through homogeneous transesterification catalyzed by Lewis acids. A systematic screening of a library of μ‐oxo‐dinuclear Fe(III) salen complexes has been conducted, aiming to elucidate the impact of steric and electronic properties of their substituents. A complementary theoretical‐computational study based on quantum mechanics calculations and machine learning regression techniques has been employed to unravel the relationship between the structure of the catalysts and their corresponding properties, shedding light on the key factors influencing their catalytic performance.

Aluminium 8-Hydroxyquinolinate N-Oxide as a Precursor to Heterometallic Aluminium-Lanthanide Complexes.

A reaction in anhydrous toluene between the formally unsaturated fragment [Ln(hfac)3] (Ln3+ = Eu3+, Gd3+ and Er3+; Hhfac = hexafluoroacetylacetone) and [Al(qNO)3] (HqNO = 8-hydroxyquinoline N-oxide), here prepared for the first time from [Al(OtBu)3] and HqNO, affords the dinuclear heterometallic compounds [Ln(hfac)3Al(qNO)3] (Ln3+ = Eu3+, Gd3+ and Er3+) in high yields. The molecular structures of these new compounds revealed a dinuclear species with three phenolic oxygen atoms bridging the two metal atoms. While the europium and gadolinium complexes show the coordination number (CN) 9 for the lanthanide centre, in the complex featuring the smaller erbium ion, only two oxygens bridge the two metal atoms for a resulting CN of 8. The reaction of [Eu(hfac)3] with [Alq3] (Hq = 8-hydroxyquinoline) in the same conditions yields a heterometallic product of composition [Eu(hfac)3Alq3]. A recrystallization attempt from hot heptane in air produced single crystals of two different morphologies and compositions: [Eu2(hfac)6Al2q4(OH)2] and [Eu2(hfac)6(µ-Hq)2]. The latter compound can be directly prepared from [Eu(hfac)3] and Hq at room temperature. Quantum mechanical calculations confirm (i) the higher stability of [Eu(hfac)3Al(qNO)3] vs. the corresponding [Eu(hfac)3Alq3] and (ii) the preference of the Er complexes for the CN 8, justifying the different behaviour in terms of the Lewis acidity of the metal centre.

Structural, Biochemical, and Computational Characterization of Sulfamides as Bimetallic Peptidase Inhibitors.

The sulfonamide function is used extensively as a general building block in various inhibitory scaffolds and, more specifically, as a zinc-binding group (ZBG) of metalloenzyme inhibitors. Here, we provide biochemical, structural, and computational characterization of a metallopeptidase in complex with inhibitors, where the mono- and bisubstituted sulfamide functions are designed to directly engage zinc ions of a bimetallic enzyme site. Structural data showed that while monosubstituted sulfamides coordinate active-site zinc ions via the free negatively charged amino group in a canonical manner, their bisubstituted counterparts adopt an atypical binding pattern divergent from expected positioning of corresponding tetrahedral reaction intermediates. Accompanying quantum mechanics calculations revealed that electroneutrality of the sulfamide function is a major factor contributing to the markedly lower potency of bisubstituted compounds by considerably lowering their interaction energy with the enzyme. Overall, while bisubstituted uncharged sulfamide functions can bolster favorable pharmacological properties of a given inhibitor, their use as ZBGs in metalloenzyme inhibitors might be less advantageous due to their suboptimal metal-ligand properties.

Triplet Excimer Formation in a DNA Duplex with Silver Ion-Mediated Base Pairs.

The dynamics of excited electronic states in self-assembled structures formed between silver(I) ions and cytosine-containing DNA strands or monomeric cytosine derivatives were investigated by time-resolved infrared (TRIR) spectroscopy and quantum mechanical calculations. The steady-state and time-resolved spectra depend sensitively on the underlying structures, which change with pH and the nucleobase and silver ion concentrations. At pH ∼ 4 and low dC20 strand concentration, an intramolecularly folded i-motif is observed, in which protons, and not silver ions, mediate C-C base pairing. However, at the higher strand concentrations used in the TRIR measurements, dC20 strands associate pairwise to yield duplex structures containing C-Ag+-C base pairs with a high degree of propeller twisting. UV excitation of the silver ion-mediated duplex produces a long-lived excited state, which we assign to a triplet excimer state localized on a pair of stacked cytosines. The computational results indicate that the propeller-twisted motifs induced by metal-ion binding are responsible for the enhanced intersystem crossing that populates the triplet state and not a generic heavy atom effect. Although triplet excimer states have been discussed frequently as intermediates in the formation of cyclobutane pyrimidine dimers, we find neither computational nor experimental evidence for cytosine-cytosine photoproduct formation in the systems studied. These findings provide a rare demonstration of a long-lived triplet excited state that is formed in a significant yield in a DNA duplex, demonstrating that supramolecular structural changes induced by metal ion binding profoundly affect DNA photophysics.

Quantum biochemical analysis of the binding interactions between a potential inhibitory drug and the Ebola viral glycoprotein

Ebola virus disease (EVD) causes outbreaks and epidemics in West Africa that persist until today. The envelope glycoprotein of Ebola virus (GP) consists of two subunits, GP1 and GP2, and plays a key role in anchoring or fusing the virus to the host cell in its active form on the virion surface. Toremifene (TOR) is a ligand that mainly acts as an estrogen receptor antagonist; however, a recent study showed a strong and efficient interaction with GP. In this context, we aimed to evaluate the energetic affinity features involved in the interaction between GP and toremifene by computer simulation techniques using the Molecular Fractionation Method with Conjugate Caps (MFCC) scheme and quantum-mechanical (QM) calculations, as well as missense mutations to assess protein stability. We identified ASP522, GLU100, TYR517, THR519, LEU186, LEU515 as the most attractive residues in the EBOV glycoprotein structure that form the binding pocket. We divided toremifene into three regions and evaluated that region i was more important than region iii and region ii for the formation of the TOR-GP1/GP2 complex, which might control the molecular remodeling process of TOR. The mutations that caused more destabilization were ARG134, LEU515, TYR517 and ARG559, while those that caused stabilization were GLU523 and ASP522. TYR517 is a critical residue for the binding of TOR, and is highly conserved among EBOV species. Our results may help to elucidate the mechanism of drug action on the GP protein of the Ebola virus and subsequently develop new pharmacological approaches against EVD. Communicated by Ramaswamy H. Sarma

Study on photochemistry and component transport coupling processes of azobenzene molecular solar thermal system

Molecular solar thermal (MOST) system is a new technology to collect, transform and store solar energy. The mechanism of photochemical reaction and component transport of molecular photoswitches in the reactor is still unclear. In this paper, a simulated theoretical framework for energy conversion rate and energy storage efficiency of solar energy collection and conversion devices based on a molecular solar thermal system was presented. The mathematical model of the photochemical reaction studied was based on an azobenzene derivative, which acted as a molecular photoswitch that can store energy by the switch of molecular structures. The microscopic parameters of reactants were obtained by quantum mechanics calculation and applied to the macroscopic energy conversion analysis of reaction processes. The theoretical framework can predict the conversion and energy storage efficiency of photoisomerization reactions, and can be used for component transport analysis, energy conversion analysis and reaction rate contribution analysis in reactors. We focused on the influence of macroscopic factors, including the flow velocity, inlet reactants concentration and reactor size on the photoisomerization energy storage process. The results of the theoretical analysis were in good agreement with those of existing laboratory equipment. The relative error of the conversion rate is only 0.71% in the photostable state, and the maximum relative error of the conversion rate is 10.99% in the whole process. In addition, in the photochemical reaction, we introduced the concept of absorbance field for the absorbance attenuation caused by the absorption of photons in the medium, so that the local material conversion, absorbance and energy conversion of this type of photochemical reaction can be more clearly understood. This simulated theoretical framework allows the study of molecules in a macroscopic setting, thereby establishing a link between practical engineering and basic chemistry, and providing guidance for the design of molecular solar thermal systems in the future.

Novichok Nerve Agents as Inhibitors of Acetylcholinesterase-In Silico Study of Their Non-Covalent Binding Affinity.

In silico studies were performed to assess the binding affinity of selected organophosphorus compounds toward the acetylcholinesterase enzyme (AChE). Quantum mechanical calculations, molecular docking, and molecular dynamics (MD) with molecular mechanics Generalized-Born surface area (MM/GBSA) were applied to assess quantitatively differences between the binding energies of acetylcholine (ACh; the natural agonist of AChE) and neurotoxic, synthetic correlatives (so-called "Novichoks", and selected compounds from the G- and V-series). Several additional quantitative descriptors like root-mean-square fluctuation (RMSF) and the solvent accessible surface area (SASA) were briefly discussed to give-to the best of our knowledge-the first quantitative in silico description of AChE-Novichok non-covalent binding process and thus facilitate the search for an efficient and effective treatment for Novichok intoxication and in a broader sense-intoxication with other warfare nerve agents as well.

Myoglobin-Catalyzed Azide Reduction Proceeds via an Anionic Metal Amide Intermediate.

Nitrene transfer reactions catalyzed by heme proteins have broad potential for the stereoselective formation of carbon-nitrogen bonds. However, competition between productive nitrene transfer and the undesirable reduction of nitrene precursors limits the broad implementation of such biocatalytic methods. Here, we investigated the reduction of azides by the model heme protein myoglobin to gain mechanistic insights into the factors that control the fate of key reaction intermediates. In this system, the reaction proceeds via a proposed nitrene intermediate that is rapidly reduced and protonated to give a reactive ferrous amide species, which we characterized by UV/vis and Mössbauer spectroscopies, quantum mechanical calculations, and X-ray crystallography. Rate-limiting protonation of the ferrous amide to produce the corresponding amine is the final step in the catalytic cycle. These findings contribute to our understanding of the heme protein-catalyzed reduction of azides and provide a guide for future enzyme engineering campaigns to create more efficient nitrene transferases. Moreover, harnessing the reduction reaction in a chemoenzymatic cascade provided a potentially practical route to substituted pyrroles.

Calcium bioaccessibility as affected by strontium. Temperature effect on citrate binding to strontium and calcium alone and in combination

Strontium as a bone health promoter was found to increase the bioaccessibility of essential calcium from solid calcium citrate tetrahydrate in aqueous suspensions. This is unexpected since calcium citrate binds calcium stronger than strontium with ΔH0ass,c = 5.07 kJ⋅mol−1, ΔS0ass,c = 78.0 J⋅mol−1⋅K−1 for 1:1 Ca-Citrate complex formation and ΔH0ass,c = 1.67 kJ⋅mol−1, ΔS0ass,c = 58.3 J⋅mol−1⋅K−1 for 1:1 Sr-Citrate complex formation as determined both electrochemically and by isothermal titration calorimetry. Calcium citrate tetrahydrate has reverse solubility around physiological temperature in contrast to strontium citrate pentahydrate. This difference results in opposite effect of strontium upon calcium ion liberation from calcium citrate suspensions at both low and high temperature compared to physiological temperature, for which strontium was found to stabilize solid calcium citrate tetrahydrate. This effect, confirmed by quantum mechanical calculation and by comparison of spontaneous aqueous supersaturation of calcium citrate and strontium citrate, is suggested to be involved in the positive effect of strontium on bone and teeth health. The conclusions provide a theoretical foundation for the research of molecular mechanisms among calcium, strontium, and citrate ions in relation to bone and teeth health.

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors.

The accurate prediction of binding free energy is a major challenge in structure-based drug design. Quantum mechanics (QM)-based approaches show promising potential in predicting ligand-protein binding affinity by accurately describing the behavior and structure of electrons. However, traditional QM calculations face computational limitations, hindering their practical application in drug design. Nevertheless, the fragment molecular orbital (FMO) method has gained widespread application in drug design due to its ability to reduce computational costs and achieve efficient ab initio QM calculations. Although the FMO method has demonstrated its reliability in calculating the gas phase potential energy, the binding of proteins and ligands also involves other contributing energy terms, such as solvent effects, the 'deformation energy' of a ligand's bioactive conformations, and entropy. Particularly in cases involving ionized fragments, the calculation of solvation free energy becomes particularly crucial. We conducted an evaluation of some previously reported implicit solvent methods on the same data set to assess their potential for improving the performance of the FMO method. Herein, we develop a new QM-based binding free energy calculation method called FMOScore, which enhances the performance of the FMO method. The FMOScore method incorporates linear fitting of various terms, including gas-phase potential energy, deformation energy, and solvation free energy. Compared to other widely used traditional prediction methods such as FEP+, MM/PBSA, MM/GBSA, and Autodock vina, FMOScore showed good performance in prediction accuracies. By constructing a retrospective case study, it was observed that incorporating calculations for solvation free energy and deformation energy can further enhance the precision of FMO predictions for binding affinity. Furthermore, using FMOScore-guided lead optimization against Src homology-2-containing protein tyrosine phosphatase 2 (SHP-2), we discovered a novel and potent allosteric SHP-2 inhibitor (compound 8).

New fusidane-type nortriterpenoids from the Arctic marine-derived fungus Simplicillium lamellicola culture medium with their inhibitory effect on benign prostatic hyperplasia

Three new fusidane-type nortriterpenoids, simplifusinolide A, 24-epi simplifusinolide A, and simplifusidic acid L (1–3), were isolated from the EtOAc extract of the Arctic marine-derived fungus Simplicillium lamellicola culture medium, together with fusidic acid (4) and 16-O-deacetylfusicid acid (5). The structures of the isolated compounds were elucidated by NMR and MS analyses. The absolute configurations of compounds 1–3 were established by the quantum mechanical calculations of electronic circular dichroism and gauge-including atomic orbital NMR chemical shifts, followed by DP4 + analysis.Benign prostatic hyperplasia (BPH) is a major urological disorder in men worldwide. The anti-BPH potentials of the isolated compounds were evaluated using BPH-1 and WPMY-1 cells. Treatment with simplifusidic acid L (3) and fusidic acid (4) significantly downregulated the mRNA levels of the androgen receptor (AR) and its downstream effectors, inhibiting the proliferation of BPH-1 cells. Specifically, treatment with 24-epi simplifusinolide A (2) significantly suppressed the cell proliferation of both BPH-1 and DHT-stimulated WPMY-1 cells by inhibiting AR signaling. These results suggest the potential of 24-epi simplifusinolide A (2), simplifusidic acid L (3) and fusidic acid (4) as alternative agents for BPH treatment by targeting AR signaling.

Effect of mechanical energy input during mechanical crushing on the macrokinetics of the coal–oxygen reaction: A laboratory–scale study

Underground mining is usually simplified as a process of coal crushing and pulverization. However, this process is not purely physical. The mechanical energy input to this process also induces a mechanochemical effect. The specific impact of this effect on the macrokinetics of the coal–oxygen reaction (COR) has not been independently investigated, which motivates this work. Thermal analysis tests showed that mechanical energy input contributed to the strength of the COR and exhibited regularities in its impact on COR's kinetic parameters. Explanations for this were given based on mechanochemistry, free radicals chain reaction and quantum chemistry theories. A coal macromolecule model, influenced by mechanical energy input, was developed based on the typical Wiser model, and subsequent quantum mechanical calculations were conducted. It was found that mechanical energy input facilitated the breaking of specific chemical bonds in coal, redistributing the electrostatic potential and increasing the number of active sites within the coal molecule. Moreover, both the free radicals from bonds breaking and the remaining active structural sites exhibited an increasing trend in their electrified capacities. Different from the traditional physical pathway, this work refined the mechanism by which mechanical forces affect COR through the mechanochemical pathway.