Search

Multiple time scales problems are investigated by combining geometrical and analytical approaches. More precisely, for fast-slow reaction-diffusion systems, we first prove the existence of slow manifolds for the abstract problem under the assumption that cross diffusion is small. This is done by extending the Fenichel theory to an infinite dimensional setting, where a main idea is to introduce a suitable space splitting corresponding to a small parameter, which controls additional fast contributions of the slow variable. These results require a strong convergence in $L^\infty(0,T;H^2(\Omega))$, which is the subtle analytical issue in fast reaction problems in comparison to previous works. By considering a nonlinear fast reversible reaction in one dimension, we successfully prove this convergence and therefore obtain the slow manifold for a PDE with fast reactions. Moreover, the obtained convergence also shows the influence of the cross diffusion term and illuminates the role of the initial layer. Explicit approximations of the slow manifold are also carried out in the case of linear systems.

The global existence of bounded solutions to reaction-diffusion systems with fractional diffusion in the whole space $\mathbb R^N$ is investigated. The systems are assumed to preserve the non-negativity of initial data and to dissipate total mass. We first show that if the nonlinearities are at most quadratic then there exists a unique global bounded solution regardless of the fractional order. This is done by combining a regularizing effect of the fractional diffusion operator and the H\"older continuity of a non-local inhomogeneous parabolic equation. When the nonlinearities might be super-quadratic, but satisfy some intermediate sum conditions, we prove the global existence of bounded solutions by adapting the well-known duality methods to the case of fractional diffusion. In this case, the order of the intermediate sum conditions depends on the fractional order. These results extend the existing theory for mass dissipated reaction-diffusion systems to the case of non-local diffusion and unbounded domains., Comment: 49 pages

The global existence and boundedness of solutions to volume-surface reaction diffusion systems with a mass control condition are investigated. Such systems arise typically in e.g. cell biology, ecology or fluid mechanics, when some concentrations or densities are inside a domain and some others are on its boundary. Comparing to previous works, the difficulty of systems under consideration here is that the nonlinearities on the surface can have a sub-quadratic growth rates in all dimensions. To overcome this, we first use the Moser iteration to get some uniform bounds of the time integration of the solutions. Then by combining these bounds with an $L^p$-energy method and a duality argument, we obtain the global existence of solutions. Moreover, under mass dissipation conditions, the solution is shown to be bounded uniformly in time., Comment: Comments are welcome

Large time dynamics of reaction-diffusion systems modeling some irreversible reaction networks are investigated. Depending on initial masses, these networks possibly possess boundary equilibria, where some of the chemical concentrations are completely used up. In the absence of these equilibria, we show an explicit convergence to equilibrium by a modified entropy method, where it is shown that reactions in a measurable set with positive measure is sufficient to combine with diffusion and to drive the system towards equilibrium. When the boundary equilibria are present, we show that they are unstable (in Lyapunov sense) using some bootstrap instability technique from fluid mechanics, while the nonlinear stability of the positive equilibrium is proved by exploiting a spectral gap of the linearized operator and the uniform-in-time boundedness of solutions., Comment: Title slightly modified. Comments are welcome

We consider hard-potential cutoff multi-species Boltzmann operators modeling microscopic binary elastic collisions and bimolecular reversible chemical reactions inside a gaseous mixture. We prove that the spectral gap estimate derived for the linearized elastic collision operator can be exploited to deduce an explicit negative upper bound for the Dirichlet form of the linearized chemical Boltzmann operator. Such estimate may be used to quantify explicitly the rate of convergence of close-to-equilibrium solutions to the reactive Boltzmann equation toward the global chemical equilibrium of the mixture.

We investigate a fast-reaction--diffusion system modelling the effect of autotoxicity on plant-growth dynamics, in which the fast-reaction terms are based on the dichotomy between healthy and exposed roots depending on the toxicity. The model was proposed in [Giannino, Iuorio, Soresina, forthcoming] to account for stable stationary spacial patterns considering only biomass and toxicity, and its fast-reaction (cross-diffusion) limit was formally derived and numerically investigated. In this paper, the cross-diffusion limiting system is rigorously obtained as the fast-reaction limit of the reaction-diffusion system with fast-reaction terms by performing a bootstrap argument involving energies. Then, a thorough well-posedness analysis of the cross-diffusion system is presented, including a $L^\infty_{x,t}$-bound, uniqueness, stability, and regularity of weak solutions. This analysis, in turn, becomes crucial to establish the convergence rate for the fast reaction limit, thanks to the key idea of using an inverse Neumann Laplacian operator. Finally, numerical experiments illustrate the analytical findings on the convergence rate.

In this paper, a reaction-diffusion system modeling injection of a chemotherapeutic drug on the surface of a living tissue during a treatment for cancer patients is studied. The system describes the interaction of the chemotherapeutic drug and the normal, tumor and immune cells. We first establish well-posedness for the nonlinear reaction-diffusion system, then investigate the long-time behavior of solutions. Particularly, it is shown that the cancer cells will be eliminated assuming that its reproduction rate is sufficiently small in a short time period in each treatment interval. The analysis is then essentially exploited to study an optimal drug injection rate problem during a chemotherapeutic drug treatment for tumor cells, which is formulated as an optimal boundary control problem with constraints. For this, we show that the existence of an optimal drug injection rate through the boundary, and derive the first-order optimality condition., Comment: 20 pages, 1 figure. We added some new references. Comments are very welcome

The trend to equilibrium for reaction-diffusion systems modelling chemical reaction networks is investigated, in the case when reaction processes happen on subsets of the domain. We prove the convergence to equilibrium by directly showing functional inequalities in terms of entropy method. Our approach allows us to deal with nonlinearities of arbitrary orders, for which only global renormalised solutions are known to globally exist. For bounded solutions, we also prove the convergence to equilibrium when the diffusion as well as the reaction are degenerate, that is both diffusion and reaction processes only act on specific subsets of the domain., Comment: Title and abstract are slightly changed, and some typos are corrected. Comments are welcome

The global existence and boundedness of solutions to quasi-linear reaction-diffusion systems are investigated. The system arises from compartmental models describing the spread of infectious diseases proposed in [Viguerie et al, Appl. Math. Lett. (2021); Viguerie et al, Comput. Mech. (2020)], where the diffusion rate is assumed to depend on the total population, leading to quasilinear diffusion with possible degeneracy. The mathematical analysis of this model has been addressed recently in [Auricchio et al, Math. Method Appl. Sci. (2023] where it was essentially assumed that all sub-populations diffuse at the same rate, which yields a positive lower bound of the total population, thus removing the degeneracy. In this work, we remove this assumption completely and show the global existence and boundedness of solutions by exploiting a recently developed $L^p$-energy method. Our approach is applicable to a larger class of systems and is sufficiently robust to allow model variants and different boundary conditions., Comment: Comments are very welcome! arXiv admin note: text overlap with arXiv:2103.16863

The stabilization by noise for parabolic equations in perforated domains, i.e. domains with small holes, is investigated. We show that when the holes are small enough, one can stabilize the unstable equations using suitable multiplicative It\^o noise. The results are quantitative, in the sense that we can explicitly estimate the size of the holes and diffusion coefficients for which stabilization by noise takes place. This is done by using the asymptotic behaviour of the first eigenvalue of the Laplacian as the hole shrinks to a point.

To model bio-chemical reaction systems with diffusion one can either use stochastic, microscopic reaction-diffusion master equations or deterministic, macroscopic reaction-diffusion system. The connection between these two models is not only theoretically important but also plays an essential role in applications. This paper considers the macroscopic limits of the chemical reaction-diffusion master equation for first-order chemical reaction systems in highly heterogeneous environments. More precisely, the diffusion coefficients as well as the reaction rates are spatially inhomogeneous and the reaction rates may also be discontinuous. By carefully discretizing these heterogeneities within a reaction-diffusion master equation model, we show that in the limit we recover the macroscopic reaction-diffusion system with inhomogeneous diffusion and reaction rates.

Reactions with enzymes are critical in biochemistry, where the enzymes act as catalysis in the process. One of the most used mechanisms for modeling enzyme-catalyzed reactions is the Michaelis-Menten (MM) kinetic. In the ODE level, i.e. concentrations are only on time-dependent, this kinetic can be rigorously derived from mass action law using quasi-steady-state approximation. This issue in the PDE setting, for instance when molecular diffusion is taken into account, is considerably more challenging and only formal derivations have been established. In this paper, we prove this derivation rigorously and obtain MM kinetic in the presence of spatial diffusion. In particular, we show that, in general, the reduced problem is a cross-diffusion-reaction system. Our proof is based on improved duality method, heat regularisation and a suitable modified energy function. To the best of our knowledge, this work provides the first rigorous derivation of MM kinetic from mass action kinetic in the PDE setting., Comment: Comments are welcome!

Bulk-surface systems on evolving domains are studied. Such problems appear typically from modelling receptor-ligand dynamics in biological cells. Our first main result is the global existence and boundedness of solutions in all dimensions. This is achieved by proving $L^p$-maximal regularity of parabolic equations and duality methods in moving surfaces, which are of independent interest. The second main result is the large time dynamics where we show, under the assumption that the volume/area of the moving domain/surface is unchanged and that the material velocities are decaying for large time, that the solution converges to a unique spatially homogeneous equilibrium. The result is proved by extending the entropy method to bulk-surface systems in evolving domains., Comment: 44 pages. Some typos are corrected. Comments are very welcome

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,\delta}(\Omega)$ for some $\delta>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 existence mass controlled systems with slightly-super cubic intermediate sum condition., Comment: Title is slightly changed. We've added a subsection in the Introduction to show the application of our results to a realistic oscillatory Belousov-Zhabotinsky reaction system. Comments are welcome!

We analyze semilinear reaction-diffusion systems that are mass controlled, and have nonlinearities that satisfy critical growth rates. The systems under consideration are only assumed to satisfy natural assumptions, namely the preservation of non-negativity and a control of the total mass. It is proved in dimension one that if nonlinearities have (slightly super-) cubic growth rates then the system has a unique global classical solutions. Moreover, in the case of mass dissipation, the solution is bounded uniformly in time in sup-norm. One key idea in the proof is the H\"older continuity of gradient of solutions to parabolic equation with possibly discontinuous diffusion coefficients and low regular forcing terms. When the system possesses additionally an entropy inequality, the global existence and boundedness of a unique classical solution is shown for nonlinearities satisfying a cubic intermediate sum condition, which is a significant generalization of cubic growth rates. The main idea in this case is to combine a modified Gagliardo-Nirenberg inequality and the newly developed $L^p$-energy method in \cite{morgan2021global,fitzgibbon2021reaction}. This idea also allows us to deal with the case of discontinuous diffusion coefficients in higher dimensions, which has \blue{only recently been touched} in the context of mass controlled reaction-diffusion systems., Comment: Title is changed, minor typos and misprints are corrected. Accepted in Journal of Evolution Equations

We propose a mathematical model, namely a reaction-diffusion system, to describe social behaviour of cockroaches. An essential new aspect in our model is that the dispersion behaviour due to overcrowding effect is taken into account {as a counterpart to commonly studied aggregation}. This consideration leads to an intriguing new phenomenon which has not been observed in the literature. Namely, due to the competition between aggregation towards areas of higher concentration of pheromone and dispersion avoiding overcrowded areas, the cockroaches aggregate more at the transition area of pheromone. Moreover, we also consider the fast reaction limit where the switching rate between active and inactive subpopulations tends to infinity. By utilising improved duality and energy methods, together with the regularisation of heat operator, we prove that the weak solution of the reaction-diffusion system converges to that of a reaction-cross-diffusion system., Comment: 41 pages, 14 figures

The global existence of renormalised solutions and convergence to equilibrium for reaction-diffusion systems with non-linear diffusion are investigated. The system is assumed to have quasi-positive non-linearities and to satisfy an entropy inequality. The difficulties in establishing global renormalised solutions caused by possibly degenerate diffusion are overcome by introducing a new class of weighted truncation functions. By means of the obtained global renormalised solutions, we study the large-time behaviour of complex balanced systems arising from chemical reaction network theory with non-linear diffusion. When the reaction network does not admit boundary equilibria, the complex balanced equilibrium is shown, by using the entropy method, to exponentially attract all renormalised solutions in the same compatibility class. This convergence extends even to a range of non-linear diffusion, where global existence is an open problem, yet we are able to show that solutions to approximate systems converge exponentially to equilibrium uniformly in the regularisation parameter., Comment: 36 pages

In this work we investigate the convergence to equilibrium for mass action reaction-diffusion systems which model irreversible enzyme reactions. Using the standard entropy method in this situation is not feasible as the irreversibility of the system implies that the concentrations of the substrate and the complex decay to zero. The key idea we utilise in this work to circumvent this issue is to introduce a family of cut-off partial entropy functions which, when combined with the dissipation of a mass like term of the substrate and the complex, yield an explicit exponential convergence to equilibrium. This method is also applicable in the case where the enzyme and complex molecules do not diffuse, corresponding to chemically relevant situation where these molecules are large in size., Comment: 45 pages, 1 figure. Minor typos are corrected and some references are added

In this paper, we study unique, globally defined uniformly bounded weak solutions for a class of semilinear reaction-diffusion-advection systems. The coefficients of the differential operators and the initial data are only required to be measurable and uniformly bounded. The nonlinearities are quasi-positive and satisfy a commonly called mass control or dissipation of mass property. Moreover, we assume the intermediate sum condition of a certain order. The key feature of this work is the surprising discovery that quasi-positive systems that satisfy an intermediate sum condition automatically give rise to a new class of $L^p$-energy type functionals that allow us to obtain requisite uniform a priori bounds. Our methods are sufficiently robust to extend to different boundary conditions, or to certain quasi-linear systems. We also show that in case of mass dissipation, the solution is bounded in sup-norm uniformly in time. We illustrate the applicability of results by showing global existence and large time behavior of models arising from a spatio-temporal spread of infectious disease., Comment: 36 pages. Accepted in SIAM Journal on Mathematical Analysis. The section of applications is significantly shortened where the last two examples are removed. Some typos are corrected

We study the global existence of classical solutions to volume-surface reaction-diffusion systems with control of mass. Such systems appear naturally from modeling evolution of concentrations or densities appearing both in a volume domain and its surface, and therefore have attracted considerable attention. Due to the characteristic volume-surface coupling, global existence of solutions to general systems is a challenging issue. In particular, the duality method, which is powerful in dealing with mass conserved systems in domains, is not applicable on its own. In this paper, we introduce a new family of $L^p$-energy functions and combine them with a suitable duality method for volume-surface systems, to ultimately obtain global existence of classical solutions under a general assumption called the \textit{intermediate sum condition}. For systems that conserve mass, but do not satisfy this condition, global solutions are shown under a quasi-uniform condition, that is, under the assumption that the diffusion coefficients are close to each other. In the case of mass dissipation, we also show that the solution is bounded uniformly in time by studying the system on each time-space cylinder of unit size, and showing that the solution is sup-norm bounded independently of the cylinder. Applications of our results include global existence and boundedness for systems arising from membrane protein clustering or activation of Cdc42 in cell polarization., Comment: 57 pages

We consider a system of reaction-diffusion equations including chemotaxis terms and coming out of the modeling of multiple sclerosis. The global existence of strong solutions to this system in any dimension is proved, and it is also shown that the solution is bounded uniformly in time. Finally, a nonlinear stability result is obtained when the chemotaxis term is not too big. We also perform numerical simulations to show the appearance of Turing patterns when the chemotaxis term is large., Comment: 3 figures

The global existence of weak solutions to a class of quasilinear parabolic equations with nonlinearities depending on first order terms and integrable data in a moving domain is investigated. The class includes the $p$-Laplace equation as a special case. Weak solutions are shown to be global by obtaining appropriate estimates on the gradient as well as a suitable version of Aubin-Lions lemma in moving domains., Comment: accepted in Nonlinear Analysis (2020)

In this work we study global well-posedness and large time behaviour for a typical reaction--diffusion system, which include degenerate diffusion, and whose non-linearities arise from chemical reactions. We show that there is an {\it indirect diffusion effect}, i.e. an effective diffusion for the non-diffusive species which is incurred by a combination of diffusion from diffusive species and reversible reactions between the species. Deriving new estimates for such degenerate reaction-diffusion system, we show, by applying the entropy method, that the solution converges exponentially to equilibrium, and provide explicit convergence rates and associated constants., Comment: Some typos are corrected

We study the uniform boundedness of solutions to reaction-diffusion systems possessing a Lyapunov-like function and satisfying an {\it intermediate sum condition}. This significantly generalizes the mass dissipation condition in the literature and thus allows the nonlinearities to have arbitrary polynomial growth. We show that two dimensional reaction-diffusion systems, with quadratic intermediate sum conditions, have global solutions which are bounded uniformly in time. In higher dimension, bounded solutions are obtained under the condition that the diffusion coefficients are {\it quasi-uniform}, i.e. they are close to each other. Applications include boundedness of solutions to chemical reaction networks with diffusion.

We consider a chemotaxis-Navier-Stokes system modelling cellular swimming in fluid drops where an exchange of oxygen between the drop and its environment is taken into account. This phenomenon results in an inhomogeneous Robin-type boundary condition. Moreover, the system is studied without the logistic growth of the bacteria population. We prove that in two dimensions, the system has a unique global classical solution, while the existence of a global weak solution is shown in three dimensions. In the latter case, we show that the energy is bounded uniformly in time. A key idea is to utilise a boundary energy to derive suitable {\it a priori} estimates. Moreover, we are able to remove the convexity assumption on the domain.

Uniform-in-time bounds of nonnegative classical solutions to reaction-diffusion systems in all space dimension are proved. The systems are assumed to dissipate the total mass and to have locally Lipschitz nonlinearities of at most (slightly super-) quadratic growth. This pushes forward the recent advances concerning global existence of reaction-diffusion systems dissipating mass in which a uniform-in-time bound has been known only in space dimension one or two. As an application, skew-symmetric Lotka-Volterra systems are shown to have unique classical solutions which are uniformly bounded in time in all dimensions with relatively compact trajectories in $C(\overline{\Omega})^m$., Comment: 17 pages

We study the global existence and uniform-in-time bounds of classical solutions in all dimensions to reaction-diffusion systems dissipating mass. By utilizing the duality method and the regularization of the heat operator, we show that if the diffusion coefficients are close to each other, or if the diffusion coefficients are large enough compared to initial data, then the local classical solution exists globally and is bounded uniformly in time. Applications of the results include the validity of the Global Attractor Conjecture for complex balanced reaction systems with large diffusion.

Energy conservations are studied for inhomogeneous incompressible and compressible Euler equations with general pressure law in a torus or a bounded domain. We provide sufficient conditions for a weak solution to conserve the energy. By exploiting a suitable test function, the spatial regularity for the density is only required to be of order $2/3$ in the incompressible case, and of order $1/3$ in the compressible case. When the density is constant, we recover the existing results for classical incompressible Euler equation., Comment: The new version deals with general pressure law for compressible equations. By using a new test function, we allow the density to be zero instead of excluding completely vacuum. Required regularities for the density in the compressible case are also reduced

The global existence of classical solutions to reaction-diffusion systems in arbitrary space dimensions is studied. The nonlinearities are assumed to be quasi-positive, to have (slightly super-) quadratic growth, and to possess a mass control, which includes the important cases as mass conservation and mass dissipation. Under these assumptions, the local classical solution is shown to be global and, in case of mass conservation or mass dissipation, to have $L^{\infty}$-norm growing at most polynomially in time. Applications include skew-symmetric Lotka-Volterra systems and quadratic reversible chemical reactions., Comment: The results are improved from mass conservation to mass control, which includes the case of mass dissipation. A new author is added

The convergence to equilibrium of renormalized solutions to reaction-cross-diffusion systems in a bounded domain under no-flux boundary conditions is studied. The reactions model complex balanced chemical reaction networks coming from mass-action kinetics and thus do not obey any growth condition, while the diffusion matrix is of cross-diffusion type and hence nondiagonal and neither symmetric nor positive semi-definite, but the system admits a formal gradient-flow or entropy structure. The diffusion term generalizes the population model of Shigesada, Kawasaki and Teramoto to an arbitrary number of species. By showing that any renormalized solution satisfies the conservation of masses and a weak entropy-entropy production inequality, it can be proved under the assumption of no boundary equilibria that {\it all} renormalized solutions converge exponentially to the complex balanced equilibrium with a rate which is explicit up to a finite dimensional inequality., Comment: The convergence of $I_6(M)$ is corrected

The problem of recovering acoustic sources, more specifically monopoles, from point-wise measurements of the corresponding acoustic pressure at a limited number of frequencies is addressed. To this purpose, a family of sparse optimization problems in measure space in combination with the Helmholtz equation on a bounded domain is considered. A weighted norm with unbounded weight near the observation points is incorporated into the formulation. Optimality conditions and conditions for recovery in the small noise case are discussed, which motivates concrete choices of the weight. The numerical realization is based on an accelerated conditional gradient method in measure space and a finite element discretization.

The stabilisation by noise on the boundary of the Chafee-Infante equation with dynamical boundary conditions subject to a multiplicative It\^o noise is studied. In particular, we show that there exists a finite range of noise intensities that imply the exponential stability of the trivial steady state. This differs from previous works on the stabilisation by noise of parabolic PDEs, where the noise acts inside the domain and stabilisation typically occurs for an infinite range of noise intensities. To the best of our knowledge, this is the first result on the stabilisation of PDEs by boundary noise., Comment: to appear in Discrete and Continuous Dynamical Systems - Series B

The large-time asymptotics of weak solutions to Maxwell--Stefan diffusion systems for chemically reacting fluids with different molar masses and reversible reactions are investigated. The diffusion matrix of the system is generally neither symmetric nor positive definite, but the equations admit a formal gradient-flow structure which provides entropy (free energy) estimates. The main result is the exponential decay to the unique equilibrium with a rate that is constructive up to a finite-dimensional inequality. The key elements of the proof are the existence of a unique detailed-balanced equilibrium and the derivation of an inequality relating the entropy and the entropy production. The main difficulty comes from the fact that the reactions are represented by molar fractions while the conservation laws hold for the concentrations. The idea is to enlarge the space of $n$ partial concentrations by adding the total concentration, viewed as an independent variable, thus working with $n+1$ variables. Further results concern the existence of global bounded weak solutions to the parabolic system and an extension of the results to complex-balanced systems.

We study the boundedness and convergence to equilibrium of weak solutions to reaction-diffusion systems with nonlinear diffusion. The nonlinear diffusion is of porous medium type and the nonlinear reaction terms are assumed to grow polynomially and to dissipate (or conserve) the total mass. By utilising duality estimates, the dissipation of the total mass and the smoothing effect of the porous medium equation, we prove that if the exponents of the nonlinear diffusion terms are high enough, then weak solutions are bounded, locally H\"older continuous and their $L^{\infty}(\Omega)$-norm grows in time at most polynomially. In order to show convergence to equilibrium, we consider a specific class of nonlinear reaction-diffusion models, which describe a single reversible reaction with arbitrarily many chemical substances. By exploiting a generalised Logarithmic Sobolev Inequality, an indirect diffusion effect and the polynomial in time growth of the $L^{\infty}(\Omega)$-norm, we show an entropy entropy-production inequality which implies exponential convergence to equilibrium in $L^p(\Omega)$-norm, for any $1\leq p < \infty$, with explicit rates and constants., Comment: 29 pages

The convergence to equilibrium for renormalised solutions to nonlinear reaction-diffusion systems is studied. The considered reaction-diffusion systems arise from chemical reaction networks with mass action kinetics and satisfy the complex balanced condition. By applying the so-called entropy method, we show that if the system does not have boundary equilibria, then any renormalised solution converges exponentially to the complex balanced equilibrium with a rate, which can be computed explicitly up to a finite dimensional inequality. This inequality is proven via a contradiction argument and thus not explicitly. An explicit method of proof, however, is provided for a specific application modelling a reversible enzyme reaction by exploiting the specific structure of the conservation laws. Our approach is also useful to study the trend to equilibrium for systems possessing boundary equilibria. More precisely, to show the convergence to equilibrium for systems with boundary equilibria, we establish a sufficient condition in terms of a modified finite dimensional inequality along trajectories of the system. By assuming this condition, which roughly means that the system produces too much entropy to stay close to a boundary equilibrium for infinite time, the entropy method shows exponential convergence to equilibrium for renormalised solutions to complex balanced systems with boundary equilibria., Comment: 25 pages

The close-to-equilibrium regularity of solutions to a class of reaction-diffusion systems is investigated. The considered systems typically arise from chemical reaction networks and satisfy a complex balanced condition. Under some restrictions on spatial dimensions ($d\leq 4$) and order of nonlinearities ($\mu = 1 + 4/d$), we show that if the initial data is close to a complex balanced equilibrium in $L^2$-norm, then classical solutions are shown global and converging exponentially to equilibrium in $L^{\infty}$-norm. Possible extensions to higher dimensions and order of nonlinearities are also discussed. The results of this paper improve the recent work [M.J. C\'aceres and J.A. Ca\~nizo, Nonlinear Analysis: TMA 159 (2017): 62-84].