Persistence Exponent for the Simple Diffusion Equation: The Exact Solution for Any Integer Dimension ()
1. Introduction
The problem in the present paper is to find the persistence exponent for the simple diffusion equation
. The diffusion equation is an equation that has no stochasticity. In the present problem, the stochasticity is introduced through the random initial conditions. The problem is about evaluating the probability of a certain event. The event is that
at a specified location remains positive throughout the time evolution till a certain time t i.e. the
at the location does not change sign even once. This probability for asymptotically large time is characterised by an exponent
called the persistence exponent. Persistence exponent for the diffusion equation has been a subject of interest to physicists [1] - [7] etc., researchers in mathematics and statistics [8] [9] etc. as well as experimentalists [10] . The interest in the persistence exponent is just not confined to the diffusion equation but to other areas of non-equilibrium physics. Among them random walk [11] , walk in a random environment with or without bias [12] , surface growth [13] , diffusing particle in a random potential with a small concentration of absorbers [14] , behaviour of financial markets [15] etc. are worth mentioning. There are few exact calculations for the persistence exponent in the literature. The case of a simple random walk in one dimension
gives the exponent
. Even the calculation of persistence exponents for Gaussian processes may not be straight forward.
We revisit the problem of simple diffusion. It is strongly non-Markovian in nature. The problem involves the partial differential equation
with random Gaussian initial conditions. It appears to remain an unsolved problem even though results [1] [2] and several others have been reported. The problem of diffusion may require a better understanding in the context of persistence. The article tries to find an exact solution to the problem.
2. Simple Diffusion Equation, Random Initial Conditions and Persistence Exponent
The diffusion equation
is a coarse grained differential equation whose solution is uniquely determined by the initial condition. In the present problem, the initial condition is not fixed but is chosen from a distribution. The initial value of
at every coordinate is chosen from a Gaussian distribution with mean 0, variance k and the initial values of
at any two coordinates are statistically independent.
In order to calculate persistence exponent
we have to calculate the probability that the field
at a specified coordinate does not flip sign even once throughout a time t. This probability
of
always remaining +ve behaves in the limit of asymptotically large time as
. This is true for a non-stationary process like in the present case. In this article any position x coordinate is a vector quantity in a d dimensional space. The moments of the initial condition distribution described above are given by
(1-a)
(1-b)
where k is the variance of the distribution. The solution for the diffusion equation may be written in terms of the initial condition as
(2)
where
. The plan for the evaluation of the exponent is as follows. First, we have to calculate the probability of
attaining a specific final value
at a certain
starting from a definite initial value
of
at
. In order to evaluate it we use the method of selective averaging. The paths that take the initial
to the final value
also comprise those where
flips sign at least once during time evolution. The probability of such paths is to be subtracted out. Finally, there has to be an integration over the final
from 0 to
, followed by an integration over
from 0 to
.
Selective averaging means averaging over the initial field
, except when
. In other words, the averaging is done over all the initial configurations such that
at
is kept fixed at
(say) i.e.
while for
varies according to Gaussian distribution. In this paper the selective distribution, denoted by subscript s, is characterized by the moments,
(3-a)
(3-b)
It may be verified from (3-a), (3-b) that if
,
,
, we get (1-a), (1-b) and for
, (3-a), (3-b) give
,
as expected. Using (3-a), (3-b), we can calculate the moments of the random variable
,
(4)
(5)
While evaluating the second order moment, we have used the relation in (3-b). Hence the mean and the variance of the distribution for
, represented by
and
respectively, are
(6-a)
(6-b)
In the above equation
denotes the angular integration in d dimensional space while
represents the usual Gamma function. It may be mentioned that
in (2) is Gaussian irrespective of whether
, the initial Gaussian field, is correlated or not. In the present case, though, the initial field is uncorrelated and
can be proved to be Gaussian using characteristic functions in probability theory [16] . It may be noted that the
function distribution is the limiting case of a Gaussian distribution. The expression for the conditional probability for starting at
and being between
and
at time
is
(7)
where
and
. This probability considers all the paths that start from
to be between
and
at time
including ones that flip en route
as depicted in Figure 1. Figure 1 is the projection of the trajectory of the system in the infinite dimensional
space on to the
plane.
represents the starting point and
, the destination. AB represents a path along which
does not flip and ADB (blue curve) is a typical path along which
flips. Such paths have to be excluded. The probability of reaching from A to the neighborhood B at asymptotically large time
without flipping is given by,
(8)
The second term represents the probability of paths such as
originating from
and terminating in the neighborhood of B at
. (8) is not be confused with the method of images in [17] . (8) follows a very different logic in the present case and holds good asymptotically. To prove (8) we will show that there is a one to one mapping from a path
to a path
and that the probability of two such paths converges asymptotically. This part is explained in 1) in what follows. Further, to justify (8) we have to show that the “number” of paths
that flip and the “number” of paths
converge asymptomatically. This is done in 2). In the subsequent analysis we will consider a d dimensional lattice—lattice spacing being infinitesimally small— instead of continuum for the sake of notational convenience only. The reason for (8) follows.
Figure 1. Projection of
trajectory onto the
plane.
1) An initial configuration at A of Figure 1 given by
, is considered, where
are the initial values of
at coordinates
. The corresponding path takes initial
to B, then it may be concluded from (2) that
(
) is a
configuration at
which takes initial
to B. Hence there is a one to one mapping of paths from
to those from
. It may be underlined here that
as
. This implies that the probability of the two paths approach each other asymptotically.
2) In this part we will address the fact that in the asymptotically large time limit it is a very good approximation to say that there is a one to one correspondence between the paths from A that flip to those from
. This may be used as it is a controlled approximation for it improves with increasing t. In order to see this point let us consider a point
(not shown in the Figure 1) where
. Let
be the initial configuration corresponding to path ADB ( the path in blue in Figure 1) where
are the initial values of
at coordinates
. This path crosses zero while reaching B. It can be
shown that
(
) is the
corresponding initial configuration for a path
. The exact expression for
contains a coordinate dependent term whose leading order behavior for large t is 1. Since
, we have
for sufficiently large
. Let the time coordinate at D be
, then
for the path ADB. Then one may arrive from (2) that
for the initial configuration
. Hence, one can conclude that the path corresponding to
must have flipped at an earlier time than
. Therefore, if a path from
flips, the corresponding path from
flips at an earlier time. Since
, the “number” of paths flipping while going from
is more than those from
. Thus the “number” of paths from
that flip is a fraction
of those from
where
for large
, a being some positive number.
On account of 1), 2) we say that the probability of the paths (like ADB in Figure 1) that flip while reaching B in the large time limit is given by
, where
,
, being Taylor expansion in
. In principle, the coefficient of
may be a function of
. When integrating over
—as will be done later—the contribution to the integral comes from the vicinity of
which is vanishingly small in the asymptotic limit. Also,
. The coefficient is Taylor expanded about
and only the zeroth order term or the term independent of
is retained. So the probability of
not changing sign when reaching the neighborhood (
) of B is, for asymtotically large time,
(9)
This leads us to (8). The final
may have any value as long as it remains positive. The probability of
starting from
and reaching a final positive value without ever changing sign is
(10)
We would now calculate (10) for asymptotically large value of t. Under the circumstances the second term on the R.H.S of (6-b) can be neglected. Further
. Hence for
,
. The expression (10) is evaluated using the identity [18]
(11)
Evaluation leads to a sum of two terms—one is proportional to
and the other is proportional to
. Hence we obtain
(12)
for
. In arriving at the above result the asymptotic expansion of “error function” erf has been used for small argument. Finally, the expression for
is obtained by integrating
over a Gaussian distribution.
(13)
where
is the Gaussian distribution for initial
with variance k as mentioned at the beginning. If
, it may be concluded from (12) and (13) that
or
depending on whether
or not. This gives
or 1.
3. Result and Conclusion
In the previous section, exact calculation has been carried out to determine the probability
of the sign of the field
remaining positive throughout an asymtotically large time t. The probability is
. Hence, the persistence exponent is
valid for any arbitrary integer dimension
. The exponents for
are 0.25, 0.50, 0.75 respectively.
The result may be experimentally verified for a system initially at thermal equilibrium defined by a temperature T. The equilibrium is then disturbed in a suitable manner. The time evolution of the coarse grained temperature at any point satisfies the simple diffusion equation; hence, this time evolution can be studied to find the persistence exponent.
The answer for the exponent
obtained in this paper is in disagreement with all the papers cited in the beginning. The first results for the persistence exponent in the case of the diffusion problem were published in [1] [2] back to back. The papers used a two time correlation function and explicitly applied the approximation (IIA), Independent Interval Approximation, to get to the answer. Application of the two time correlation function is not suitable here and so is IIA which is a Markovian approximation. Further, the papers use Monte Carlo simulation to confirm the result. Monte Carlo method appears to be unsuitable for this problem. Hence all the papers that reproduce the results of [1] [2] are not expected to give the correct answer. In [3] , the authors have defined a correlation function
, just like [1] [2] , to carry the calculation forward. Let us also consider [6] where the authors have used Kac Polynomials [19] to obtain the “exact exponent” in 2d. The answer obtained agrees perfectly with [1] [2] . In the course of the calculation, they have used that the zero crossing property is
governed by the covariance
[6] of the stationary Gaussian
process i.e. the diffusion equation with time redefined. Similarly in [4] , correlator in time
has been used in the calculation. The point is that the covariance/correlator/correlation function is a misleading quantity for the problem for reasons mentioned below. The model presented in the paper has randomness only in the initial condition. Once the system starts evolving, there is no further randomness. It evolves in accordance with the kernel in (2). It is encoded in the initial condition when and where the
will flip. The probability of each path is uniquely determined by the probability of initial condition, hence the problem with covariance/correlator. The covariance function imposes stochasticity on the present problem throughout the entire time evolution. We now have a different model with the same correlation function but no unique dependence of the probability of the path on initial condition. It also makes the problem Markovian. Hence, all the previous results are in perfect agreement though the calculated exponent will be different from the actual value. The value of the exponent does not depend on only correlation function, but it depends on other details of the model too. Further, there even appears to be experimental proof [10] for the results of [1] [2] . The experimental setup of [10] does not represent the diffusion model described in this paper. The setup satisfies the approximations of the previous papers and hence the agreement with their result.
Acknowledgements
The author would like to thank CSIR, India for Fellowship during the course of the work (2004) at IACS, India.