Persistence Exponent for the Simple Diffusion Equation: The Exact Solution for any Integer Dimension

The persistence exponent $\theta_o$ for the simple diffusion equation ${\phi}_t({\it x},t) = \triangle \phi (x,t)$ , with random Gaussian initial condition {\color{red},} has been calculated exactly using a method known as selective averaging. The probability that the value of the field $\phi$ at a specified spatial coordinate remains positive throughout for a certain time $t$ behaves as $t^{-\theta_o}$ for asymptotically large time $t$. The value of $\theta_o$, calculated here for any integer dimension $d$, is $\theta_o = \frac{d}{4}$ for $d\leq 4$ and $1$ otherwise. This exact theoretical result is being reported possibly for the first time and is not in agreement with the accepted values $ \theta_o = 0.12, 0.18,0.23$ for $d=1,2,3$ respectively.


I. INTRODUCTION
The problem in the present paper is to find the persistence exponent for the simple diffusion equation φ t (x , t) = △φ(x, t). 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 φ 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 θ o called the persistence exponent. Persistence exponent for the diffusion equation has been a subject of interest to physicists [1][2][3][4][5][6][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 θ o = 1 2 . 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 φ t = △φ 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 * deva sans@yahoo.co.in may require a better understanding in the context of persistence. The article tries to find an exact solution to the problem.

II. SIMPLE DIFFUSION EQUATION , RANDOM INITIAL CONDITIONS AND PERSISTENCE EXPONENT
The diffusion equation φ t = △φ 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 θ o 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 P + (t) of φ always remaining +ve behaves in the limit of asymptotically large time as P + (t) ∼ t −θo . 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 are given by where k is the variance ofthe distribution. The solution for the diffusion equation may be written in terms of the initial condition as where G(x, t) = (4πt) −d/2 exp(−x 2 /4t). 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 x = x o starting from a definite initial value α of φ at x = x o . 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 φ(x o ) flips sign atleast 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 φ(x, 0), except when x = x o . In other words, the averaging is done over all the initial configurations such that φ In this paper the selective distribution, denoted by subscript s, is characterized by the moments, While evaluating the second order moment, we have used the relation in (3-b). Hence the mean and the variance of the distribution for φ(x o , t), represented by µ and σ 2 respectively, are In the above equation k d denotes the angular integration in d dimensional space while Γ represents the usual Gamma function. It may be mentioned that φ(x, t) in (2) is Gaussian irrespective of whether φ(x ′ , 0) , the initial Gaussian field, is correlated or not. In the present case,though, the initial field is uncorrelated and φ(x, t) 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 β + dβ at time t 1 is where µ = µ(α, t 1 ) and σ = σ(t 1 ). This probability considers all the paths that start from α to be between β and β + dβ at time t 1 including ones that flip en route β as depicted in Fig 1. Fig 1 is the projection of the trajectory of the system in the A(0, α) represents the starting point and B(t 1 , β), the destination. AB represents a path along which φ(x o ) does not flip and ADB (blue curve) is a typical path along which φ(x o ) flips. Such paths have to be excluded. The probability of reaching from A to the neighborhood B at asymptotically large time t 1 without flipping is given by, The second term represents the probability of paths such asĀDB originating fromĀ(0, −α) and terminating in the neighborhood of B at t 1 . (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 A → B to a pathĀ → B and that the probability of two such paths converge asymptotically. This part is explained in (i) in what follows. Further, to justify (8) we have to show that the "number" of paths A → B that flip and the "number" of pathsĀ → B converge asymptomatically.This is done in (ii). 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) Hence there is a one to one mapping of paths from A → B to those fromĀ → B . It may be underlined here that f → 1 as t → ∞. This implies that the probability of the two paths approach each other asymptotically. ii) 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Ā → B.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 C (t 2 , β)(not shown in the F ig 1) where t 2 > t 1 . Let Y AB = {...γ 1 , α, γ 2 , ...} be the initial configuration corresponding to path ADB ( the path in blue in Fig 1) where γ 1 ,γ 2 ... are the initial values of φ at coordinates x = x o . This path crosses zero while reaching B. It can be shown that Y AC = {...f 1 γ 1 , α, f 1 γ 2 , ...} is the corresponding initial configuration for a path A → C. The exact expression for f 1 contains a coordinate dependent term whose leading order behavior for large t is 1. Since t 2 > t 1 , we have f 1 > 1 for sufficiently large t 1 . Let the time coordinate at D be t D , then φ(x o , t D ) = 0 for the path ADB. Then one may arrive from (2) that φ(x o , t D ) < 0 for the initial configuration Y AC .Hence, one can conclude that the path corresponding to Y AC must have flipped at an earlier time than t D . Therefore, if a path from A → B flips, the corresponding path from A → C flips at an earlier time. Since t 2 > t 1 , the 'number' of paths flipping while going from A → C is more than those from A → B. Thus the 'number' of paths from A → B that flip is a fraction f 2 of those fromĀ → B where −a ) for large t 1 ,a being some positive number.
On account of (i), (ii) we say that the probability of the paths( like ADB in Fig 1) that flip while reaching B in the large time limit is given by P β| − α h correction , where h correction = 1+O(t −b ) , b = 1, being Taylor expansion in t −1 . In principle, the coefficient of t −b may be a function of β. When integrating over β -as will be done laterthe contribution to the integral comes from the vicinity of β = −µ ∼ t −d/2 which is vanishingly small in the asymptotic limit. Also, dβ ∼ σ ∼ t −d/4 . The coefficient is Taylor expanded about β = 0 and only the zeroth order term or the term independent of β is retained. So the probability of φ(x o ) not changing sign when reaching the neighborhood(dβ) of B is, for asymtotically large time, This leads us to (8). The final β may have any value as long as it remains positive. The probability of φ(x o ) starting from α and reaching a final positive value without ever changing sign is 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 µ 2 σ 2 ∼ αt −d/2 . Hence for α ≪ t d/2 , µ 2 /σ 2 ≪ 1. The expression (10) is evaluated using the identity [18] Evaluation leads to a sum of two terms -one is proportional to t −d/4 and the other is proportional to t −1 . Hence we obtain for d ≤ 4. In arriving at the above result the asymptotic expansion of 'error function' erf has been used for small argument. Finally, the expression for P + (t) is obtained by integrating α over a Gaussian distribution.
where Q(α) is the Gaussian distribution for initial φ(x o , o) = α with variance k as mentioned at the beginning. If k << t d/2 , it may be concluded from (12) and (13)  In the previous section exact calculation has been carried out to determine the probability P + (t) of the sign of the field φ remaining positive through out a asymtotically large time t. The probability is P + (t) ∼ t −d/4 . Hence the persistence exponent is θ o = d/4 valid for any arbitrary integer dimension d ≤ 4. The exponents for d = 1, 2 , 3 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 θ o 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 C(T )-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 c(T ) = sech( T 2 ) [6] of the stationary Gaussian process i.e the diffusion equation with time redefined. Similarly in [4], correlator in time F ǫ (τ − τ ′ ) 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 intial 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 the their result.