We mathematically prove that chemical reaction networks without hidden layers can solve tasks for which spiking neural networks require hidden layers. Our proof uses the deterministic mass-action kinetics formulation of chemical reaction networks. Specifically, we prove that a certain reaction network without hidden layers can learn a classification task previously proved to be achievable by a spiking neural network with hidden layers. We provide analytical regret bounds for the global behavior of the network and analyze its asymptotic behavior and Vapnik-Chervonenkis dimension. In a numerical experiment, we confirm the learning capacity of the proposed chemical reaction network for classifying handwritten digits in pixel images, and we show that it solves the task more accurately and efficiently than a spiking neural network with hidden layers. This provides a motivation for machine learning in chemical computers and a mathematical explanation for how biological cells might exhibit more efficient learning behavior within biochemical reaction networks than neuronal networks.
Reaction networks are a general framework widely used in modeling diverse phenomena in different science disciplines. The dynamical process of a reaction network endowed with mass-action kinetics is a mass-action system as an ODE defined by a directed graph, the so-called ``reaction graph''. Endotacticity is a graph property used to study persistence and permanence of mass-action systems. In this paper, we provide a detailed characterization of first order endotactic reaction graphs. Besides, we provide a sufficient condition for endotacticity of reaction networks which are not necessarily of first order. Such a characterization of a first order endotactic reaction graph yields the spectral property of the adjacency matrix of the reaction graph. As a consequence, we prove that every first order endotactic mass-action system as a linear ODE has a weakly reversible deficiency zero realization, and has a unique equilibrium which is exponentially globally asymptotically stable (and is positive) in each (positive) stoichiometric compatibility class. Using a stability result for asymptotically autonomous differential equations, examples are constructed to illustrate that the global stabi
We derive an operational-time variance kernel for a latent-order-book reaction boundary and use it to separate three objects usually collapsed in calendar-time volatility models: a structural boundary cumulant, a clock projection, and a pricing-measure choice. The reaction boundary is the zero of a bid--ask imbalance field. For a locally linear book, signed order-flow perturbations displace this zero through a damped Abel response kernel, so the variance of boundary increments is obtained as a finite-scale Green-function cumulant rather than introduced as a primitive diffusion coefficient. For long-memory forcing with exponent $0<γ<1$, the operational variance has a closed asymptotic form involving effective signed-forcing intensity, liquidity slope, resilience, memory, and operational coarse-graining scale. A deterministic activity clock gives the benchmark local-volatility projection. More general, non-unique clocks generate candidate calendar-time pricing systems. We argue that such projections are admissible only when the induced forward density operator and backward valuation operator remain adjoint on the same state space. Adjoint consistency is therefore a reality cons
Classical reaction dynamics suggests transport through an index-1 saddle is organized not just by flux, but by local symplectic width scales of bounded proxy neighborhoods near the bottleneck. We investigate if a related geometric effect appears in the quantum regime for highly squeezed Gaussian wavepackets. Building on de Gosson's symplectic approach, we analyze how transverse bath-mode squeezing modifies transmission across a quantum normal-form (QNF) bottleneck. To avoid the instability of propagating states with extreme phase-space eccentricity, we use the Weyl-symbol formulation of the QNF. For the quadratic saddle-center model, we derive an exact baseline transmission formula by convolving the bath's squeezed-state number distribution with the 1D Kemble transmission factor. For anharmonic truncated QNF models, we enforce strict algebraic energy conservation and evaluate exact Gaussian expectation-value diagnostics of the Weyl symbol via Wick-Isserlis moment formulas. Results reveal a pronounced squeeze-induced suppression of transmission. As the squeezed state's bath-plane geometric scale grows relative to the classical candidate width, the expected bath action grows rapidly.
This paper investigates a system of nonlinear reaction-diffusion equations modeling the industrial synthesis of ammonia. By applying Lie group analysis, we construct self-similar solutions and derive a reduced system of ordinary differential equations. Using comparison principles and barrier techniques, we establish sufficient conditions for the existence of global-in-time solutions in both slow-diffusion ($γ_i > 0$) and fast-diffusion ($γ_i < 0$) regimes. Detailed asymptotic analysis near the diffusion front reveals power-law behavior of concentration profiles, with explicit expressions for the decay exponents. The theoretical results are illustrated by numerical simulations, demonstrating the spatio-temporal evolution of reactant concentrations under realistic parameter values. The study provides rigorous mathematical foundations for predicting and optimizing ammonia production in catalytic reactors, with potential extensions to other chemically reacting systems.
In 1917, Marian von Smoluchowski presented a simple mathematical description of diffusion-controlled reactions on the scale of individual molecules. His model postulated that a reaction would occur when two reactants were sufficiently close and, more specifically, presented a succinct relationship between the relative proximity of two reactants at the moment of reaction and the macroscopic reaction rate. Over the last century, Smoluchowski reaction theory has been applied widely in the physical, chemical, environmental and, more recently, the biological sciences. Despite the widespread utility of the Smoluchowski theory, it only describes the rates of second order reactions and is inadequate for the description of higher order reactions for which there is no equivalent method for theoretical investigation. In this paper, we derive a generalised Smoluchowski framework in which we define what should be meant by proximity in this context when more than two reactants are involved. We derive the relationship between the macroscopic reaction rate and the critical proximity at which a reaction occurs for higher order reactions. Using this theoretical framework and using numerical experime
The origin of life is often approached through the lens of replication, heredity, or molecular specificity. This paper proposes a thermodynamic framework in which the emergence of life is driven by the persistence of reaction pathways that align energetically with fluctuating environmental inputs. We define a reaction viability inequality based on energy input, release, resilience, and expenditure, which selects for persistent chemical configurations without invoking heredity or genetic encoding. We further incorporate entropic dynamics and spatial constraints into an augmented persistence function, showing that systems far from equilibrium can simultaneously increase global entropy while supporting localized chemical order. These refinements lead to the development of the Thermodynamic Abiogenesis Likelihood Model (TALM), a probabilistic extension that estimates the likelihood of persistence-driven selection under diverse prebiotic and planetary scenarios. This framework redefines the conditions under which life-like organization may emerge and provides a testable, general theory for abiogenesis grounded in physical law.
Couplings between biochemical and mechanical processes have a profound impact on embryonic development. However, in-vitro studies capable of quantifying these interactions have remained elusive. Here, we investigate a synthetic system where a DNA reaction-diffusion (RD) front is advected by a turbulent flow generated by active matter (AM) flows in a quasi-one-dimensional geometry. Whereas the dynamics of simple RD fronts solely depend on the reaction and diffusion rates, we show that RD-AM front propagation is also influenced by the confinement geometry. We first experimentally dissected the different components of the reaction-diffusion-advection process by knocking out reaction or advection and observed how RD-AM allows for faster transport over large distances, avoiding dilution. We then show how confinement impacts active matter flow: while changes in instantaneous flow velocities are small; correlation times are dramatically increased with decreasing confinement. As a result, RD-AM front speed increased up to 8-fold compared to a RD one, in quantitative agreement with a conveyor belt RDA theoretical model. The RD-AM experimental system described here provides a framework for t
Continuous time Markov chains are commonly used as models for the stochastic behavior of chemical reaction networks. More precisely, these Stochastic Chemical Reaction Networks (SCRNs) are frequently used to gain a mechanistic understanding of how chemical reaction rate parameters impact the stochastic behavior of these systems. One property of interest is mean first passage times (MFPTs) between states. However, deriving explicit formulas for MFPTs can be highly complex. In order to address this problem, we first introduce the concept of coclique level structure and develop theorems to determine whether certain SCRNs have this feature by studying associated graphs. Additionally, we develop an algorithm to identify, under specific assumptions, all possible coclique level structures associated with a given SCRN. Finally, we demonstrate how the presence of such a structure in a SCRN allows us to derive closed form formulas for both upper and lower bounds for the MFPTs. Our methods can be applied to SCRNs taking values in a generic finite state space and can also be applied to models with non-mass-action kinetics. We illustrate our results with examples from the biological areas of ep
We present a vector-based method to balance chemical reactions. The algorithm builds candidates in a deterministic way, removes duplicates, and always prints coefficients in the lowest whole-number form. For redox cases, electrons and protons/hydroxide are treated explicitly, so both mass and charge are balanced. We also outline the basic principles of the vector formulation of stoichiometry, interpreting reactions as integer vectors in composition space, this geometric view supports compact visualizations of reagent-product interactions and helps surface distinct reaction families. The method enumerates valid balances for arbitrary user-specified species lists without special-case balancing rules or symbolic tricks, and it provides a clean foundation for developing new algorithmic variants (e.g., alternative objectives or constraints). On representative examples (neutralization, double displacement, decomposition, classical redox, small multicomponent sets) and a negative control, the method produced correct integer balances. When multiple balances exist, we report a canonical one - minimizing the total coefficient sum with a simple tie-breaker - without claiming global optimality
A stress-affected chemical reaction front propagation is considered utilizing the concept of a chemical affinity tensor. A reaction between an elastic solid and diffusing constituents, localized at the reaction front, is considered. As a result of the reaction, the elastic constituent transforms into viscoelastic one. The reaction is accompanied by volume expansion that in turn may result in stresses at the reaction front which affect the front velocity through the normal component of the chemical affinity tensor. Considering a plane strain problem with a planar chemical reaction front propagation under uniaxial deformation, we focus on the studies of the reaction front kinetics in dependence on external strains and material parameters with the use of the notion of the equilibrium concentration. Then stress relaxation behind the propagating reaction front is modelled. A standard linear solid model is used for the reaction product, and its particular cases are also considered. Analytical solutions are obtained which allow to study in explicit form the influence of strains and material parameters on the front retardation or acceleration and stress relaxation.
In this paper we prove that positive weak solutions for quasilinear parabolic equations on bounded domains subject to homogenous Neumann boundary conditions becme classical and global under the unique condition that the reaction doesn't change sign after certain positive time. We apply this result to reaction diffusion systems and prove global existence of theirs positive weak solutions under the same condition on theirs reactions. The nonlinearities growth isn't taken in consideration. The proof is based on the maximum principle.
The concept of reaction extent (the progress of a reaction, advancement of the reaction, conversion, etc.) was introduced around 100 years ago. Most of the literature provides a definition for the exceptional case of a single reaction step or gives an implicit definition that cannot be made explicit. There are views that the reaction extent somehow has to tend to 1 when the reaction goes to completion as time tends to infinity. However, there is no agreement on which function should tend to 1. Starting from the standard definition by IUPAC and following the classical works by De Donder, Aris, and Croce we extend the classic definition of the reaction extent for an arbitrary number of species and reaction steps. The new general, explicit definition is also valid for non-mass action kinetics. We also studied the mathematical properties (evolution equation, continuity, monotony, differentiability, etc.) of the defined quantity, connecting them to the formalism of modern reaction kinetics. Our approach tries to adhere to the customs of chemists and be mathematically correct simultaneously. To make the exposition easy to understand, we use simple chemical examples and many figures, thro
The past few decades have seen robust research on questions regarding the existence, form, and properties of stationary distributions of stochastically modeled reaction networks. When a stochastic model admits a stationary distribution an important practical question is: what is the rate of convergence of the distribution of the process to the stationary distribution? With the exception of \cite{XuHansenWiuf2022} pertaining to models whose state space is restricted to the non-negative integers, there has been a notable lack of results related to this rate of convergence in the reaction network literature. This paper begins the process of filling that hole in our understanding. In this paper, we characterize this rate of convergence, via the mixing times of the processes, for two classes of stochastically modeled reaction networks. Specifically, by applying a Foster-Lyapunov criteria we establish exponential ergodicity for two classes of reaction networks introduced in \cite{anderson2018some}. Moreover, we show that for one of the classes the convergence is uniform over the initial state.
Reaction-diffusion systems with mass dissipation are known to possess blow-up solutions in high dimensions when the nonlinearities have super quadratic growth rates. In dimension one, it has been shown recently that one can have global existence of bounded solutions if nonlinearities are at most cubic. For the cubic intermediate sum condition, i.e. nonlinearities might have arbitrarily high growth rates, an additional entropy inequality had to be imposed. In this article, we remove this extra entropy assumption completely and obtain global boundedness for reaction-diffusion systems with cubic intermediate sum condition. The novel idea is to show a non-concentration phenomenon for mass dissipating systems, that is the mass dissipation implies a dissipation in a Morrey space $\mathsf{M}^{1,δ}(Ω)$ for some $δ>0$. As far as we are concerned, it is the first time such a bound is derived for mass dissipating reaction-diffusion systems. The results are then applied to obtain global existence and boundedness of solutions to an oscillatory Belousov-Zhabotinsky system, which satisfies cubic intermediate sum condition but does not fulfill the entropy assumption. Extensions include global e
Studies of one-body density matrices (ODM) are performed in various correlation methods, such as the Jastrow method, the correlated basis function method, the Green's function method and the generator coordinate method aiming to extract the absolute spectroscopic factors and overlap functions (OF) for one-nucleon removal reactions from the ODM of the target nucleus. The advantage of this method is that it avoids the complicated task of calculating the total nuclear spectral function. The procedure for extracting bound-state OF's has been applied to make calculations of the cross sections of the $(e,e^{\prime}p)$ reaction on the closed-shell nuclei $^{16}$O and $^{40}$Ca as well as on the open-shell nucleus $^{32}$S consistently (using the same OF's) with the cross sections of $(p,d)$ and $(γ,p)$ reactions on the same nuclei. The analyses of the reaction cross sections and the spectroscopic factors and the comparison with the experimental data show the particular importance of these OF's, since they contain effects of nucleon correlations (short-range and/or long-range) which are accounted for to different extent in the theoretical methods considered.
Time-dependent Hartree-Fock (TDHF) theory has achieved a remarkable success in describing and understanding nuclear many-body dynamics from nucleons' degrees of freedom. We here report our investigation of multinucleon transfer (MNT) processes employing the TDHF theory. To calculate transfer probabilities for channels specified by the number of protons and neutrons included in reaction products, a particle-number projection (PNP) method has been developed. The PNP method is also used to calculate excitation energies of reaction products. Combined use of the PNP method with a statistical model, we can evaluate MNT cross sections taking account of effects of particle evaporation. Using these methods, we evaluate MNT cross sections for $^{40,48}$Ca+$^{124}$Sn, $^{40}$Ca+$^{208}$Pb, and $^{58}$Ni+$^{208}$Pb reactions. From systematic analyses, we find that cross sections for channels with a large reaction probability are in good agreement with experimental data. However, the agreement becomes less accurate as the number of transferred nucleons increases. Possible directions to improve the description are discussed.
This paper is concerned with the well-posedness and optimal control problem of a reaction-diffusion system for an epidemic Susceptible-Infected-Recovered-Susceptible (SIRS) mathematical model in which the dynamics develops in a spatially heterogeneous environment. Using as control variables the transmission rates $u_{i}$ and $u_{e}$ of contagion resulting from the contact with both asymptomatic and symptomatic persons, respectively, we optimize the number of exposed and infected individuals at a final time $T$ of the controlled evolution of the system. More precisely, we search for the optimal $u_{i}$ and $u_{e}$ such that the number of infected plus exposed does not exceed at the final time a threshold value $Λ$, fixed a priori. We prove here the existence of optimal controls in a proper functional framework and we derive the first-order necessary optimality conditions in terms of the adjoint variables.
Can a micron sized sack of interacting molecules understand, and adapt to a constantly-fluctuating environment? Cellular life provides an existence proof in the affirmative, but the principles that allow for life's existence are far from being proven. One challenge in engineering and understanding biochemical computation is the intrinsic noise due to chemical fluctuations. In this paper, we draw insights from machine learning theory, chemical reaction network theory, and statistical physics to show that the broad and biologically relevant class of detailed balanced chemical reaction networks is capable of representing and conditioning complex distributions. These results illustrate how a biochemical computer can use intrinsic chemical noise to perform complex computations. Furthermore, we use our explicit physical model to derive thermodynamic costs of inference.
This paper is concerned with the well-posedness of a diffusion-reaction system for a Susceptible-Exposed-Infected-Recovered (SEIR) mathematical model. This model is written in terms of four nonlinear partial differential equations with nonlinear diffusions, depending on the total amount of the SEIR populations. The model aims at describing the spatio-temporal spread of the COVID-19 pandemic and is a variation of the one recently introduced, discussed and tested in [A. Viguerie et al, Diffusion-reaction compartmental models formulated in a continuum mechanics framework: application to COVID-19, mathematical analysis, and numerical study, Comput. Mech. 66 (2020) 1131-1152]. Here, we deal with the mathematical analysis of the resulting Cauchy-Neumann problem: the existence of solutions is proved in a rather general setting and a suitable time discretization procedure is employed. It is worth mentioning that the uniform boundedness of the discrete solution is shown by carefully exploiting the structure of the system. Uniform estimates and passage to the limit with respect to the time step allow to complete the existence proof. Then, two uniqueness theorems are offered, one in the case