The combustion and detonation processes of energetic materials exhibit remarkable complexity and ultra-fast transient characteristics.While reactive molecular dynamics has been extensively employed to investigate the ...The combustion and detonation processes of energetic materials exhibit remarkable complexity and ultra-fast transient characteristics.While reactive molecular dynamics has been extensively employed to investigate the reaction dynamics of energetic materials,its utility is often constrained to capturing only fundamental reaction events and species information,thereby limiting mechanistic investigations of complex reaction pathways.To elucidate the topological features of energetic material reaction networks and identify critical reaction pathways with high fidelity,this study presents ReacNetwork-an advanced large-scale reaction network analysis methodology that synergistically integrates complex network theory with molecular simulation techniques.Specifically,we have developed a multi-dimensional feature screening protocol based on node centrality metrics and K-shell decomposition algorithms.Takingα-Hexahydro-1,3,5-trinitro-1,3,5-triazine(α-RDX)as the subject,we successfully constructed a comprehensive high-temperature thermal decomposition reaction network consisting of 1,134 distinct chemical species and 3,626 elementary reactions.Through systematic application of community detection algorithms and global topological feature extraction techniques,we achieved effective dimensionality reduction and successfully identified the dominant reaction pathway within theα-RDX thermal decomposition network.The computational results not only validate the well-established initial reaction mechanism dominated by N-NO2 homolytic bond cleavage,but also provide unprecedented visualization ofα-RDX framework ring-opening dynamics and subsequent radical chain propagation networks.展开更多
Peng et al.in[Phys.Rev.Research,2020,2(3):033089,11 pp.]formulated one-way fluxes for a general chemical reaction far from equilibrium,with arbitrary complex mechanisms,multiple intermediates,and internal kinetic cycl...Peng et al.in[Phys.Rev.Research,2020,2(3):033089,11 pp.]formulated one-way fluxes for a general chemical reaction far from equilibrium,with arbitrary complex mechanisms,multiple intermediates,and internal kinetic cycles.They defined the limit of the ratio of mesoscopic one-way fluxes and the volume of the tank reactor when the volume tends to infinity as macroscopic one-way fluxes,but a rigorous proof of existence of the limit is still awaiting.In this article,we fill this gap under a mild hypothesis:the Markov chain associated with the chemical master equation has finite states and any two columns in the stoichiometric matrices are not identical.In fact,an explicit expression of the limit is obtained.展开更多
The formation of humins hampers the large-scale production of 5-hydroxymethylfurfural(HMF)in biorefinery.Here,a detailed reaction network of humin formation at the initial stage of fructose-to-HMF dehydration in water...The formation of humins hampers the large-scale production of 5-hydroxymethylfurfural(HMF)in biorefinery.Here,a detailed reaction network of humin formation at the initial stage of fructose-to-HMF dehydration in water is delineated by combined experimental,spectroscopic,and theoretical studies.Three bimolecular reaction pathways to build up soluble humins are demonstrated.That is,the intermolecular etherification of β-furanose at room temperature initiates the C12 path,whereas the C-C cleavage of a-furanose at 130-150℃ leads to C11 path,and that of open-chain fructose at 180℃ to C11' path.The successive intramolecular dehydrations and condensations of the as-formed bimolecular intermediates lead to three types of soluble humins.We show that the C12 path could be restrained by using HCl or AlCl_(3) catalyst,and both the C12 and C110 paths could be effectively inhibited by adding THF as a co-solvent or accelerating heating rate via microwave heating.展开更多
The oxidative esterification of methacrolein(MAL)is an important way to prepare high-valued methyl methacrylate(MMA),but this process is ultra-complex due to the high reactivity of both C=O and C=C bonds in MAL molecu...The oxidative esterification of methacrolein(MAL)is an important way to prepare high-valued methyl methacrylate(MMA),but this process is ultra-complex due to the high reactivity of both C=O and C=C bonds in MAL molecule.In order to further improve MMA selectivity,the reaction network and relevant mechanisms have been proposed and profoundly investigated in this paper.Five kinds of fundamental reactions are involved in this process,including(a)the acetal reaction;(b)the aerobic oxidation of hemiacetal;(c)the alkoxylation of C=C double bond;(d)the Diels-Alder reaction;and(e)the hydrogenation reaction of unsaturated double bond.Among them,the Diels-Alder reaction of MAL is non-catalyzed,and the Brönsted acid sites or the Lewis acid sites favor promoting acetal reaction of MAL with methanol,while the alkoxylation of C=C bond with methanol is enhanced under alkaline condition.In particular,by employing the Pd-based catalysts,hydrogenation products are formed in alkaline methanol solution,hence with lower than those obtained by the Au-based catalysts.Notably,it is necessary to match the hemiacetal fromation and aerobic oxidation of hemiacetal,which is relevant with the amount and strength of acid and redox sites.Consequently,this work can provide a good guidance for the further design of both catalysts and processes in future.展开更多
The treatment of a multicomponent reversible reaction network is extremely complicated because largenumber of rate constants must be precisely determined and because the calculation based on these rateconstants is ted...The treatment of a multicomponent reversible reaction network is extremely complicated because largenumber of rate constants must be precisely determined and because the calculation based on these rateconstants is tedious.In order to reduce the degrees of freedom of the process,the authors propose a methodin which the reactor and the separator are regarded as a whole.Based on this approach,an N-componentreversible reaction system can be dealt with as a two—component system.Consequently,a simple and ac-cessible way of the apparent rate determination is suggested.For fiist-order reactions,an explicit,simplifiedexpression has been derived for both lumped and distributed parameter reaction systems.展开更多
Fischer-Tropsch synthesis(FTS)wax is a mixture of linear hydrocarbons with carbon number from C7 to C70+.Converting FTS wax into high-quality diesel(no sulfur and nitrogen contents)by hydrocracking technology is attra...Fischer-Tropsch synthesis(FTS)wax is a mixture of linear hydrocarbons with carbon number from C7 to C70+.Converting FTS wax into high-quality diesel(no sulfur and nitrogen contents)by hydrocracking technology is attractive in economy and practicability.Kinetic study of the hydrocracking of FTS wax in elementary step level is very challenging because of the huge amounts of reactions and species involved.Generation of reaction networks for hydrocracking of FTS wax in which the chain length goes up to C70 is described on the basis of Boolean adjacency matrixes.Each of the species(including paraffins,olefins and carbenium ions)involved in the elementary steps is represented digitally by using a(N+3)N matrix,in which a group of standardized numbering rules are designed to guarantee the unique identity of the species.Subsequently,the elementary steps are expressed by computer-aided matrix transformations in terms of proposed reaction rules.Dynamic memory allocation is used in species storage and a characteristic vector with nine elements is designed to store the key information of a(N+3)N matrix,which obviously reduces computer memory consumption and improves computing efficiency.The detailed reaction networks of FTS wax hydrocracking can be generated smoothly and accurately by the current method.The work is the basis of advanced elementary-step-level kinetic modeling.展开更多
Biological systems use intricate networks of chemical reactions to exchange information. How to simulate complex systems with simple strand-displacement reactions is crucial to broaden the application scenario of the ...Biological systems use intricate networks of chemical reactions to exchange information. How to simulate complex systems with simple strand-displacement reactions is crucial to broaden the application scenario of the DNA reaction network. Here, we report the artificial DNA reaction network to mimic the operation and function of biological information transfer via strand-displacement reaction. DNA is used as simple artificial analogs to schematize structures and transmit information. Using chemical synapses in neural networks as an example, we show that the proposed network enables core functions of biological systems, such as the long-term potential of synapses, which underpin learning and memory. Also, we performed the simple “silicon mimetic” to link electronic circuits to chemical network-based biological structures. As such, synaptic communication simulated by the DNA reaction network provides a complete demonstration for designing artificial reaction networks based on the essence of information interaction.展开更多
This paper investigates the global exponential stability of reaction-diffusion neural networks with discrete and distributed time-varying delays. By constructing a more general type of Lyapunov-Krasovskii functional c...This paper investigates the global exponential stability of reaction-diffusion neural networks with discrete and distributed time-varying delays. By constructing a more general type of Lyapunov-Krasovskii functional combined with a free-weighting matrix approach and analysis techniques, delay-dependent exponential stability criteria are derived in the form of linear matrix inequalities. The obtained results are dependent on the size of the time-vaxying delays and the measure of the space, which are usually less conservative than delay-independent and space-independent ones. These results are easy to check, and improve upon the existing stability results. Some remarks are given to show the advantages of the obtained results over the previous results. A numerical example has been presented to show the usefulness of the derived linear matrix inequality （LMI）-based stability conditions.展开更多
In this paper, the influence of the noise and delay upon the stability property of reaction-diffusion recurrent neural networks （RNNs） with the time-varying delay is discussed. The new and easily verifiable conditio...In this paper, the influence of the noise and delay upon the stability property of reaction-diffusion recurrent neural networks （RNNs） with the time-varying delay is discussed. The new and easily verifiable conditions to guarantee the mean value exponential stability of an equilibrium solution are derived. The rate of exponential convergence can be estimated by means of a simple computation based on these criteria.展开更多
As a key component of shale oil,petroleum fractions,and chemical products,the oxidative pyrolysis behavior of paraffin directly influences energy conversion efficiency and the direction of process optimization.A deep ...As a key component of shale oil,petroleum fractions,and chemical products,the oxidative pyrolysis behavior of paraffin directly influences energy conversion efficiency and the direction of process optimization.A deep understanding of its oxidative pyrolysis mechanism is crucial for addressing wax deposition in oil and gas extraction,enhancing product selectivity in cracking processes,and advancing novel clean fuel technologies.Traditional experimental methods face challenges in capturing transient free-radical reaction pathways at high temperatures,whereas molecular dynamics simulations offer a powerful approach to bridge the research gap in elucidating atomic-scale dynamic mechanisms.This database is constructed based on high-precision molecular dynamics simulations,comprising oxidative pyrolysis trajectory data for three paraffin models featuring different straight-chain hydrocarbon distributions within the temperature range of 2100-2500 K.The COMPASS force field was employed to optimize the initial structures,and the ReaxFF reactive force field was used to simulate the oxidative pyrolysis process.The database includes atomic trajectories,species evolution information,and reaction network analysis results for both heating and isothermal cracking processes,with a total data volume of approximately 141 GB(including 150000 atomic configuration frames).The data is stored in a hierarchical directory structure,supporting multi-scale oxidative pyrolysis mechanism studies and providing atomic-scale dynamic evidence for revealing carbon chain length effects and temperature sensitivity.展开更多
Based on the experimental data relating to the reaction of FCC gasoline on acid catalyst the analysis of product distribution, and composition of gasoline and diesel fractions have been analyzed. The occurrence of dis...Based on the experimental data relating to the reaction of FCC gasoline on acid catalyst the analysis of product distribution, and composition of gasoline and diesel fractions have been analyzed. The occurrence of disproportionation reaction of FCC gasoline on acid catalyst and the network of disproportionation reaction have been identified. Study has also shown that different reaction temperatures can result in different pathways of disproportionation reactions on acid catalyst.展开更多
Inferior crude oil and fuel oil upgrading lead to escalating increase of hydrogen consumption in refineries.It is imperative to reduce the hydrogen consumption for energy-saving operations of refineries.An integration...Inferior crude oil and fuel oil upgrading lead to escalating increase of hydrogen consumption in refineries.It is imperative to reduce the hydrogen consumption for energy-saving operations of refineries.An integration strategy of hydrogen network and an operational optimization model of hydrotreating(HDT)units are proposed based on the characteristics of reaction kinetics of HDT units.By solving the proposed model,the operating conditions of HDT units are optimized,and the parameters of hydrogen sinks are determined by coupling hydrodesulfurization(HDS),hydrodenitrification(HDN)and aromatic hydrogenation(HDA)kinetics.An example case of a refinery with annual processing capacity of eight million tons is adopted to demonstrate the feasibility of the proposed optimization strategies and the model.Results show that HDS,HDN and HDA reactions are the major source of hydrogen consumption in the refinery.The total hydrogen consumption can be reduced by 18.9%by applying conventional hydrogen network optimization model.When the hydrogen network is optimized after the operational optimization of HDT units is performed,the hydrogen consumption is reduced by28.2%.When the benefit of the fuel gas recovery is further considered,the total annual cost of hydrogen network can be reduced by 3.21×10~7CNY·a^(-1),decreased by 11.9%.Therefore,the operational optimization of the HDT units in refineries should be imposed to determine the parameters of hydrogen sinks base on the characteristics of reaction kinetics of the hydrogenation processes before the optimization of the hydrogen network is performed through the source-sink matching methods.展开更多
基金support from the National Natural Science Foundation of China(Grant No.22275018).
文摘The combustion and detonation processes of energetic materials exhibit remarkable complexity and ultra-fast transient characteristics.While reactive molecular dynamics has been extensively employed to investigate the reaction dynamics of energetic materials,its utility is often constrained to capturing only fundamental reaction events and species information,thereby limiting mechanistic investigations of complex reaction pathways.To elucidate the topological features of energetic material reaction networks and identify critical reaction pathways with high fidelity,this study presents ReacNetwork-an advanced large-scale reaction network analysis methodology that synergistically integrates complex network theory with molecular simulation techniques.Specifically,we have developed a multi-dimensional feature screening protocol based on node centrality metrics and K-shell decomposition algorithms.Takingα-Hexahydro-1,3,5-trinitro-1,3,5-triazine(α-RDX)as the subject,we successfully constructed a comprehensive high-temperature thermal decomposition reaction network consisting of 1,134 distinct chemical species and 3,626 elementary reactions.Through systematic application of community detection algorithms and global topological feature extraction techniques,we achieved effective dimensionality reduction and successfully identified the dominant reaction pathway within theα-RDX thermal decomposition network.The computational results not only validate the well-established initial reaction mechanism dominated by N-NO2 homolytic bond cleavage,but also provide unprecedented visualization ofα-RDX framework ring-opening dynamics and subsequent radical chain propagation networks.
基金partially supported by NSFC(Nos.11701265,11961033).
文摘Peng et al.in[Phys.Rev.Research,2020,2(3):033089,11 pp.]formulated one-way fluxes for a general chemical reaction far from equilibrium,with arbitrary complex mechanisms,multiple intermediates,and internal kinetic cycles.They defined the limit of the ratio of mesoscopic one-way fluxes and the volume of the tank reactor when the volume tends to infinity as macroscopic one-way fluxes,but a rigorous proof of existence of the limit is still awaiting.In this article,we fill this gap under a mild hypothesis:the Markov chain associated with the chemical master equation has finite states and any two columns in the stoichiometric matrices are not identical.In fact,an explicit expression of the limit is obtained.
基金the National Natural Science Foundation of China(No.21875149)111 project(B17030)the Basal Research Fund of the Central University.
文摘The formation of humins hampers the large-scale production of 5-hydroxymethylfurfural(HMF)in biorefinery.Here,a detailed reaction network of humin formation at the initial stage of fructose-to-HMF dehydration in water is delineated by combined experimental,spectroscopic,and theoretical studies.Three bimolecular reaction pathways to build up soluble humins are demonstrated.That is,the intermolecular etherification of β-furanose at room temperature initiates the C12 path,whereas the C-C cleavage of a-furanose at 130-150℃ leads to C11 path,and that of open-chain fructose at 180℃ to C11' path.The successive intramolecular dehydrations and condensations of the as-formed bimolecular intermediates lead to three types of soluble humins.We show that the C12 path could be restrained by using HCl or AlCl_(3) catalyst,and both the C12 and C110 paths could be effectively inhibited by adding THF as a co-solvent or accelerating heating rate via microwave heating.
基金The authors acknowledge supports from the National Key Basic Research Development Plan“973”Project(2006CB202508)the SINOPEC Project(411058,413025)the National Natural Science Foundation(21808244).
文摘The oxidative esterification of methacrolein(MAL)is an important way to prepare high-valued methyl methacrylate(MMA),but this process is ultra-complex due to the high reactivity of both C=O and C=C bonds in MAL molecule.In order to further improve MMA selectivity,the reaction network and relevant mechanisms have been proposed and profoundly investigated in this paper.Five kinds of fundamental reactions are involved in this process,including(a)the acetal reaction;(b)the aerobic oxidation of hemiacetal;(c)the alkoxylation of C=C double bond;(d)the Diels-Alder reaction;and(e)the hydrogenation reaction of unsaturated double bond.Among them,the Diels-Alder reaction of MAL is non-catalyzed,and the Brönsted acid sites or the Lewis acid sites favor promoting acetal reaction of MAL with methanol,while the alkoxylation of C=C bond with methanol is enhanced under alkaline condition.In particular,by employing the Pd-based catalysts,hydrogenation products are formed in alkaline methanol solution,hence with lower than those obtained by the Au-based catalysts.Notably,it is necessary to match the hemiacetal fromation and aerobic oxidation of hemiacetal,which is relevant with the amount and strength of acid and redox sites.Consequently,this work can provide a good guidance for the further design of both catalysts and processes in future.
文摘The treatment of a multicomponent reversible reaction network is extremely complicated because largenumber of rate constants must be precisely determined and because the calculation based on these rateconstants is tedious.In order to reduce the degrees of freedom of the process,the authors propose a methodin which the reactor and the separator are regarded as a whole.Based on this approach,an N-componentreversible reaction system can be dealt with as a two—component system.Consequently,a simple and ac-cessible way of the apparent rate determination is suggested.For fiist-order reactions,an explicit,simplifiedexpression has been derived for both lumped and distributed parameter reaction systems.
基金supported by the National Key Research&Development Program of China(2020YFB0606404)National Natural Science Foundation of China(21908234)。
文摘Fischer-Tropsch synthesis(FTS)wax is a mixture of linear hydrocarbons with carbon number from C7 to C70+.Converting FTS wax into high-quality diesel(no sulfur and nitrogen contents)by hydrocracking technology is attractive in economy and practicability.Kinetic study of the hydrocracking of FTS wax in elementary step level is very challenging because of the huge amounts of reactions and species involved.Generation of reaction networks for hydrocracking of FTS wax in which the chain length goes up to C70 is described on the basis of Boolean adjacency matrixes.Each of the species(including paraffins,olefins and carbenium ions)involved in the elementary steps is represented digitally by using a(N+3)N matrix,in which a group of standardized numbering rules are designed to guarantee the unique identity of the species.Subsequently,the elementary steps are expressed by computer-aided matrix transformations in terms of proposed reaction rules.Dynamic memory allocation is used in species storage and a characteristic vector with nine elements is designed to store the key information of a(N+3)N matrix,which obviously reduces computer memory consumption and improves computing efficiency.The detailed reaction networks of FTS wax hydrocracking can be generated smoothly and accurately by the current method.The work is the basis of advanced elementary-step-level kinetic modeling.
基金financially supported by the National Natural Science Foundation of China(Nos.82172372 and 21904045)the Fundamental Research Funds for the Central Universities(2019kfyXJS169).
文摘Biological systems use intricate networks of chemical reactions to exchange information. How to simulate complex systems with simple strand-displacement reactions is crucial to broaden the application scenario of the DNA reaction network. Here, we report the artificial DNA reaction network to mimic the operation and function of biological information transfer via strand-displacement reaction. DNA is used as simple artificial analogs to schematize structures and transmit information. Using chemical synapses in neural networks as an example, we show that the proposed network enables core functions of biological systems, such as the long-term potential of synapses, which underpin learning and memory. Also, we performed the simple “silicon mimetic” to link electronic circuits to chemical network-based biological structures. As such, synaptic communication simulated by the DNA reaction network provides a complete demonstration for designing artificial reaction networks based on the essence of information interaction.
基金supported by the National Natural Science Foundation of China (Grant No. 60974139)partially supported by the Fundamental Research Funds for the Central Universities
文摘This paper investigates the global exponential stability of reaction-diffusion neural networks with discrete and distributed time-varying delays. By constructing a more general type of Lyapunov-Krasovskii functional combined with a free-weighting matrix approach and analysis techniques, delay-dependent exponential stability criteria are derived in the form of linear matrix inequalities. The obtained results are dependent on the size of the time-vaxying delays and the measure of the space, which are usually less conservative than delay-independent and space-independent ones. These results are easy to check, and improve upon the existing stability results. Some remarks are given to show the advantages of the obtained results over the previous results. A numerical example has been presented to show the usefulness of the derived linear matrix inequality （LMI）-based stability conditions.
文摘In this paper, the influence of the noise and delay upon the stability property of reaction-diffusion recurrent neural networks （RNNs） with the time-varying delay is discussed. The new and easily verifiable conditions to guarantee the mean value exponential stability of an equilibrium solution are derived. The rate of exponential convergence can be estimated by means of a simple computation based on these criteria.
基金Supported by Natural Science Foundation of Shanxi Province (202203021221219)Research on the Construction of Scientific and Technological Innovation Think Tank of Shanxi Association for Science and Technology (KXKT202542)Planning Project under Commerce Statistical Society of China (2025STY122)。
文摘As a key component of shale oil,petroleum fractions,and chemical products,the oxidative pyrolysis behavior of paraffin directly influences energy conversion efficiency and the direction of process optimization.A deep understanding of its oxidative pyrolysis mechanism is crucial for addressing wax deposition in oil and gas extraction,enhancing product selectivity in cracking processes,and advancing novel clean fuel technologies.Traditional experimental methods face challenges in capturing transient free-radical reaction pathways at high temperatures,whereas molecular dynamics simulations offer a powerful approach to bridge the research gap in elucidating atomic-scale dynamic mechanisms.This database is constructed based on high-precision molecular dynamics simulations,comprising oxidative pyrolysis trajectory data for three paraffin models featuring different straight-chain hydrocarbon distributions within the temperature range of 2100-2500 K.The COMPASS force field was employed to optimize the initial structures,and the ReaxFF reactive force field was used to simulate the oxidative pyrolysis process.The database includes atomic trajectories,species evolution information,and reaction network analysis results for both heating and isothermal cracking processes,with a total data volume of approximately 141 GB(including 150000 atomic configuration frames).The data is stored in a hierarchical directory structure,supporting multi-scale oxidative pyrolysis mechanism studies and providing atomic-scale dynamic evidence for revealing carbon chain length effects and temperature sensitivity.
文摘Based on the experimental data relating to the reaction of FCC gasoline on acid catalyst the analysis of product distribution, and composition of gasoline and diesel fractions have been analyzed. The occurrence of disproportionation reaction of FCC gasoline on acid catalyst and the network of disproportionation reaction have been identified. Study has also shown that different reaction temperatures can result in different pathways of disproportionation reactions on acid catalyst.
基金Supported by the National Natural Science Foundation of China(21376188,21676211)the Key Project of Industrial Science and Technology of Shaanxi Province(2015GY095)
文摘Inferior crude oil and fuel oil upgrading lead to escalating increase of hydrogen consumption in refineries.It is imperative to reduce the hydrogen consumption for energy-saving operations of refineries.An integration strategy of hydrogen network and an operational optimization model of hydrotreating(HDT)units are proposed based on the characteristics of reaction kinetics of HDT units.By solving the proposed model,the operating conditions of HDT units are optimized,and the parameters of hydrogen sinks are determined by coupling hydrodesulfurization(HDS),hydrodenitrification(HDN)and aromatic hydrogenation(HDA)kinetics.An example case of a refinery with annual processing capacity of eight million tons is adopted to demonstrate the feasibility of the proposed optimization strategies and the model.Results show that HDS,HDN and HDA reactions are the major source of hydrogen consumption in the refinery.The total hydrogen consumption can be reduced by 18.9%by applying conventional hydrogen network optimization model.When the hydrogen network is optimized after the operational optimization of HDT units is performed,the hydrogen consumption is reduced by28.2%.When the benefit of the fuel gas recovery is further considered,the total annual cost of hydrogen network can be reduced by 3.21×10~7CNY·a^(-1),decreased by 11.9%.Therefore,the operational optimization of the HDT units in refineries should be imposed to determine the parameters of hydrogen sinks base on the characteristics of reaction kinetics of the hydrogenation processes before the optimization of the hydrogen network is performed through the source-sink matching methods.
