Simplified approach for calculating moments of action for linear reaction-diffusion equations

• Published in: Physical review. E, Statistical, nonlinear, and soft matter physics Vol. 88; no. 5; p. 054102
• Main Authors: Ellery, Adam J, Simpson, Matthew J, McCue, Scott W, Baker, Ruth E
• Format: Journal Article
• Language: English
• Published: United States 01.11.2013
• ISSN: 1550-2376
• DOI: 10.1103/PhysRevE.88.054102

The mean action time is the mean of a probability density function that can be interpreted as a critical time, which is a finite estimate of the time taken for the transient solution of a reaction-diffusion equation to effectively reach steady state. For high-variance distributions, the mean action time underapproximates the critical time since it neglects to account for the spread about the mean. We can improve our estimate of the critical time by calculating the higher moments of the probability density function, called the moments of action, which provide additional information regarding the spread about the mean. Existing methods for calculating the nth moment of action require the solution of n nonhomogeneous boundary value problems which can be difficult and tedious to solve exactly. Here we present a simplified approach using Laplace transforms which allows us to calculate the nth moment of action without solving this family of boundary value problems and also without solving for the transient solution of the underlying reaction-diffusion problem. We demonstrate the generality of our method by calculating exact expressions for the moments of action for three problems from the biophysics literature. While the first problem we consider can be solved using existing methods, the second problem, which is readily solved using our approach, is intractable using previous techniques. The third problem illustrates how the Laplace transform approach can be used to study coupled linear reaction-diffusion equations.