View source for Reaction-diffusion equation

A system of partial differential equations of the form

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800101.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800102.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800103.png" /> is the [[Laplace operator|Laplace operator]] in the spatial variables <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800104.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800105.png" /> is a non-negative non-zero diagonal matrix, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800106.png" /> is a function from a domain in <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800107.png" /> into <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800108.png" />. Many generalizations of these equations have also been studied, such as result when <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r0800109.png" /> depends also on the first-order <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001010.png" />-derivatives of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001011.png" />, when the operator <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001012.png" /> is replaced by other, possibly non-linear, operators, or when the matrix <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001013.png" /> is not diagonal. If extra first-order terms appear in the system as a model for convective transport effects, the system is sometimes termed reaction-advection-diffusion equation.

Such equations arise as models of diverse natural phenomena [[#References|[a1]]], but their most natural roots lie in the study of chemical systems: the components of the vector <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001014.png" /> may then represent concentrations of chemical species which are present, the term <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001015.png" /> represents the diffusive transport of those species, possibly through a chemical solution, and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001016.png" /> represents the production or destruction of species resulting from reactions among them (if the rates of all such reactions are known, as functions of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001017.png" />, then the explicit form of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001018.png" /> can be written down).

The variable <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001019.png" /> is often confined to a domain <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001020.png" /> with boundary <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001021.png" />, and then solutions are sought which satisfy specific boundary conditions on <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001022.png" />. These are generally of the form

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001023.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001024.png" /> is the derivative normal to <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001025.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001026.png" /> and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001027.png" /> are not both zero (unless <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001028.png" /> does not  "diffuse" ), and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001029.png" /> is a given function. Again, generalizations, such as to non-linear boundary conditions, abound.

The specific problems of interest are: i) the initial value problem, in which <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001030.png" /> is given and <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001031.png" /> is sought for <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001032.png" />; ii) the steady problem, in which solutions independent of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001033.png" /> are sought; and iii) the travelling-wave problem, in which <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001034.png" /> and solutions are sought of the special form <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001035.png" />.

Because of their strong connections with the applied sciences and the limited number of important properties common to all members of this unwieldy class of systems, the research impetus in this field comes more from viewing the systems as models of specific natural phenomena, rather than from interest in them for their own sake. A typical motivation, for example, may be to ask whether a certain system, in which specific natural effects are modelled, will have solutions which reflect some known natural phenomenon of interest whose causes are incompletely known. Then one looks for the existence and the stability of solutions of the system in question which have properties analogous to the phenomenon in question.

Regarding the initial value problem i), the theory of analytic semi-groups, which in this context relies on the operator <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001036.png" /> being sectorial, has developed as one of the most commonly used approaches to existence and uniqueness [[#References|[a2]]]. The study of steady solutions ii) has used a variety of methods, such as recasting the problem as a fixed-point problem for a mapping in some suitable function space and using topological-degree methods. In the case when <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001037.png" /> or the system has certain monotonicity properties, methods based on upper-and-lower solutions (cf. [[Upper-and-lower-functions method|Upper-and-lower-functions method]]) provide easier alternatives (see [[#References|[a1]]] and [[#References|[a3]]], e.g.).

The travelling-wave problem iii) can be viewed as seeking a parameter <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001038.png" /> for which there exists a connection between two critical points for the system of ordinary differential equations resulting from the substitution <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001039.png" />. One of the principal tools in this connection has been the powerful Conley index [[#References|[a4]]], [[#References|[a3]]]. Below some other methods which have been devised more recently are mentioned.

The theory of reaction-diffusion systems can be viewed as incorporating all of the theory of autonomous ordinary differential systems <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001040.png" /> (cf. [[Autonomous system|Autonomous system]]), since when homogeneous Neumann boundary conditions are imposed, solutions of the latter system automatically constitute <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001041.png" />-independent solutions of the corresponding reaction-diffusion systems. But, of course, solutions with striking spatial characteristics arise as well; and,in fact, the possible spatial structure of solutions is one of the most often investigated aspects.

Some of the best-studied examples of reaction-diffusion systems are the following.

a) The scalar Fisher equation

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001042.png" /></td> </tr></table>

[[#References|[a5]]], [[#References|[a6]]], where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001043.png" /> has exactly two zeros. This equation originally arose in connection with population genetics.

b) The scalar bistable diffusion equation, [[#References|[a7]]], [[#References|[a6]]], [[#References|[a8]]], of the same form but where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001044.png" /> has exactly three simple zeros and is negative between the first two. This equation also has connections with population genetics, but knowledge of its properties is even more important in connection with the role it plays as a component part of more complicated systems.

c) The FitzHugh–Nagumo system

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001045.png" /></td> </tr></table>

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001046.png" /></td> </tr></table>

where <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001047.png" /> has the properties given in b) (see the references in [[#References|[a3]]] and [[#References|[a9]]] for this and generalizations). It is a simplification of higher-order systems such as the Hodgkin–Huxley system, which arise as models of signal transmission on nerve axons and in cardiac tissue.

d) The thermal-diffusion model in chemical reactor theory and combustion [[#References|[a10]]], [[#References|[a11]]]. In this model, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001048.png" />, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001049.png" /> represents temperature, the other components of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001050.png" /> represent concentrations of chemical species, and the components of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001051.png" /> are given by

<table class="eq" style="width:100%;"> <tr><td valign="top" style="width:94%;text-align:center;"><img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001052.png" /></td> </tr></table>

the summation being over all reactions occurring in the material (these reactions are indexed by <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001053.png" />). Here, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001054.png" /> is a (mass action) monomial in <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001055.png" /> appropriate to the <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001056.png" />-th reaction, <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001057.png" /> is the  "reaction constant"  for that reaction, and the numbers <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001058.png" /> are  "stoichiometric parameters" , specifying the amount of species <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001059.png" /> (or heat, in case <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001060.png" />) produced or consumed in reaction <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001061.png" />.

In all of the above examples, the existence and the stability of travelling-wave solutions is of paramount importance; and in case d), other spatially or temporally ordered solutions are important as well.

In many applications, solutions with frontal or interfacial properties arise [[#References|[a9]]]. For example, a moving surface in <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001062.png" />-space may exist near which some components of <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001063.png" /> experience dramatic changes. These changes form an interior layer at the surface in question. They have been studied in the context of phase-field equations (a reaction-diffusion system with non-diagonal <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001064.png" />), in which they represent phase interfaces, of the bistable equation, and of the FitzHugh–Nagumo equations and their generalizations, in which they may represent phase changes, changes in the electrochemical properties of neural or cardiac tissue, or changes in the chemical state of a medium.

The stability of waves in one space dimension for systems with <img align="absmiddle" border="0" src="https://www.encyclopediaofmath.org/legacyimages/r/r080/r080010/r08001065.png" /> is a considerably more difficult field of investigation than is their existence. Much of the work here has been done for the FitzHugh–Nagumo equations [[#References|[a12]]]. Recently a new technique, the stability index of J. Alexander, R. Gardner and C. Jones [[#References|[a13]]], was developed and applied to a number of travelling-wave problems.

For travelling-wave and stationary solutions with interfaces (see above), a technique called the SLEP method has been developed to study stability questions (see [[#References|[a4]]] and the references therein).

For FitzHugh–Nagumo and related systems, the most important patterned solutions in two space dimensions are rotating spirals, which are extremely prevalent and apparently extremely stable structures for this and many other models for excitable media. Despite the great interest in these rotating solutions and the large number of papers the concept has generated (see the references in [[#References|[a9]]]), their mathematical foundation is still rudimentary. Analogous phenomena exist in three dimensions: structures which rotate about curves in space, called filaments, which themselves migrate according to certain approximate laws. An important challenge for the future is to better understand (and provide a firm mathematical foundation for) the connections between such dynamic spatial patterns and their laws of motion, on the one hand, and the underlying partial differential equations on the other.

====References====
<table><TR><TD valign="top">[a1]</TD> <TD valign="top">  P.C. Fife,   "Mathematical aspects of reacting and diffusing systems" , ''Lect. notes in biomath.'' , '''28''' , Springer  (1979)</TD></TR><TR><TD valign="top">[a2]</TD> <TD valign="top">  D. Henry,   "Geometric theory of semilinear parabolic equations" , ''Lect. notes in math.'' , '''840''' , Springer  (1981)</TD></TR><TR><TD valign="top">[a3]</TD> <TD valign="top">  J. Smoller,   "Shock waves and reaction-diffusion equations" , Springer  (1983)</TD></TR><TR><TD valign="top">[a4]</TD> <TD valign="top">  C. Conley,   "Isolated invariant sets and the Morse index" , Amer. Math. Soc.  (1978)</TD></TR><TR><TD valign="top">[a5]</TD> <TD valign="top">  A.N. Kolmogorov,   I.G. Petrovskii,   N.S. Piskunov,   "A study of the diffusion equation with increase in the quantity of matter, and its application to a biological problem"  ''Byull. Moskov. Gos. Univ.'' , '''17'''  (1937)  pp. 1–72  (In Russian)</TD></TR><TR><TD valign="top">[a6]</TD> <TD valign="top">  D.G. Aronson,   H.F. Weinberger,   "Multidimensional nonlinear diffusion arising in population genetics"  ''Adv. in Math.'' , '''30'''  (1978)  pp. 33–76</TD></TR><TR><TD valign="top">[a7]</TD> <TD valign="top">  Ya.I. Kanel',   "On the stability of solutions of the Cauchy problem for the equations arising in combustion theory"  ''Mat. Sb.'' , '''59'''  (1962)  pp. 245–288  (In Russian)</TD></TR><TR><TD valign="top">[a8]</TD> <TD valign="top">  P.C. Fife,   J.B. McLeod,   "The approach of solutions of nonlinear diffusion equations to travelling front solutions"  ''Arch. Rational Mech. Anal.'' , '''65'''  (1977)  pp. 335–361</TD></TR><TR><TD valign="top">[a9]</TD> <TD valign="top">  P.C. Fife,   "Dynamics of internal layers and diffusive interfaces" , ''CBMS-NSF Reg. Conf. Ser. Appl. Math.'' , '''53''' , SIAM  (1988)</TD></TR><TR><TD valign="top">[a10]</TD> <TD valign="top">  R. Aris,   "The mathematical theory of diffusion and reaction in permeable catalysts" , '''1–2''' , Clarendon Press  (1975)</TD></TR><TR><TD valign="top">[a11]</TD> <TD valign="top">  J. Buckmaster,   G.S.S. Ludford,   "Theory of laminar flow" , Cambridge Univ. Press  (1982)</TD></TR><TR><TD valign="top">[a12]</TD> <TD valign="top">  C. Jones,   "Stability of the travelling wave solution of the FitzHugh–Nagumo system"  ''Trans. Amer. Math. Soc.'' , '''286'''  (1984)  pp. 431–469</TD></TR><TR><TD valign="top">[a13]</TD> <TD valign="top">  J. Alexander,   R. Gardner,   C. Jones,   "A topological invariant arising in the stability analysis of travelling waves"  ''J. Reine Angew. Math.''  (To appear)</TD></TR><TR><TD valign="top">[a14]</TD> <TD valign="top">  Y. Nishiura,   M. Mimura,   H. Ikeda,   H. Fujii,   "Singular limit analysis of stability of travelling wave solutions in bistable reaction-diffusion systems"  ''SIAM J. Math. Anal.'' , '''21'''  (1990)  pp. 85–122</TD></TR></table>