Maximal distance travelled by N vicious walkers till their survival arxiv: v1 [cond-mat.stat-mech] 17 Feb 2014

Transcription

1 Noname manuscript No. will be inserted by the editor Anupam Kundu Satya N. Majumdar Grégory Schehr Maximal distance travelled by N vicious walkers till their survival arxiv:4.3964v [cond-mat.stat-mech] 7 Feb 4 October, 8 Abstract We consider N Brownian particles moving on a line starting from initial positions u {u, u,... u N } such that < u < u < < u N. Their motion gets stopped at time t s when either two of them collide or when the particle closest to the origin hits the origin for the first time. For N =, we study the probability distribution function p m u and p m u of the maximal distance travelled by the st and nd walker till t s. For general N particles with identical diffusion constants D, we show that the probability distribution p N m u of the global maximum m N, has a power law tail p i m u N B N F N u/m ν N with exponent ν N = N +. We obtain explicit expressions of the function F N u and of the N dependent amplitude B N which we also analyze for large N using techniques from random matrix theory. We verify our analytical results through direct numerical simulations. Introduction Extreme value statistics EVS is by now a major issue with a variety of applications in several areas of sciences including physics, statistics or finance, to name just a few []. For N independent and identically distributed i.i.d. random variables y,, y N, the distribution of the maximum y max = max i N y i or the minimum y min is well understood with the identification, in the large N thermodynamical limit, of three distinct universality classes, depending on the parent distribution of the y i s []. However, these results for i.i.d. random variables do not apply when the random variables are correlated [, 3]. Recently, there has been a surge of interest in EVS of strongly correlated random variables, which is very often the interesting case in statistical physics. Physically relevant examples include for instance the extreme statistics of a stochastic process ut, with strong temporal correlations, like Brownian motion or its variants. Many studies in this context are focused on extremal properties, like the maximum of ut, over a fixed time interval, t [, T ] [4,5,6,7,8,9,]. However in many cases the length of this time interval is itself a random variable t s, which can thus vary from one realization of the stochastic process to another. This time t s is usually strongly A. Kundu Laboratoire de Physique Théorique et Modèles Statistiques, Université Paris-Sud, Bât., 945 Orsay Cedex, France. anupamoraja@gmail.com S. N. Majumdar Laboratoire de Physique Théorique et Modèles Statistiques, Université Paris-Sud, Bât., 945 Orsay Cedex, France. satya.majumdar@u-psud.fr G. Schehr Laboratoire de Physique Théorique et Modèles Statistiques, Université Paris-Sud, Bât., 945 Orsay Cedex, France. gregory.schehr@u-psud.fr

2 correlated to the process ut itself. An interesting situation is the case where t s is a stopping time [], i.e. when it is associated to the stopping of the process ut if a certain event occurs for the first time. For example, in a queuing process starting from an initial queue length l >, t s is the time when the queue length l t becomes zero for the first time also called the busy period [,3]. In finance t s might correspond to the time when a stock price S t reaches some specified level for the first time [5, 4, 5]. Stopping times also naturally arise in various statistical physics models ranging from capture processes [6, 7] or target annihilating problems [8] all the way to reaction-diffusion kinetics [9,, ] or coarsening dynamics of domain walls in Ising model []. In the context of stochastic control theory, stochastic processes with stopping time have been widely studied [3]. The simplest example of a stopped stochastic process is the motion of a single Brownian particle starting form u = u > which is observed till time t s when the walker crosses the origin for the first time. This time is called the first passage time [4,5]. A natural extreme value question is then: what is the distribution pm u of the maximal displacement m = max t ts ut travelled by the walker till its first passage time t s? It can be shown [6] that the cumulative probability Q u L = Prob[m L u ] that the maximum stays below L till the first passage time is given by Q u L = u /L, hence pm u = u /m. In the context of polymer translocation through a small pore, the quantity Q u L is precisely the probability of complete translocation of a polymer of length L. For generic subdiffusive process, this translocation probability is shown to scale as u /L φ for large L with φ = θ p /H where θ p is the persistence exponent [5,7] and H is the Hurst exponent [8]. Other related questions like the statistics of the time when the walker reaches the maximal displacement before its first passage time t s or the fluctuations of the area enclosed under the Brownian motion till t s, have also been studied in connection with several applications including queuing theory or lattice polygon models [,3,6,9,3,3]. Stopped processes involving N > particles are also interesting and have been considered in the literature. For instance, the maximal displacement between the leader" and the laggard" among N particles has been studied for N = 3 particles in Ref. [3]. Very recently the authors of Ref. [33] have studied the probability distribution function PDF pm u of the global maximum m N of N non-interacting and identical Brownian walkers i.e. with the same diffusion constant before their first exit from the positive half-line, given that they had started from positions u {u, u,... u N }. They showed that the tail of the PDF pm u of the global maximum m N till the time t s when any one of the N walkers crosses the origin for the first time, is given by pm u [ N b N N i= u i ] m N+, m u N, where b N is an N-dependent constant that behaves for large N as, b N exp [ N loglog N]. This result holds for non-interacting particles and it is natural to wonder about the effects of interactions on the statistics of the global maximum till the stopping time of this multi-particle process. This is precisely the question which we address in this article, by considering non-intersecting Brownian motions, which is one of the simplest though non trivial interacting particles system. More precisely, we consider N Brownian walkers moving on a line with position u i t at time t for i =,,..., N. They evolve with time according to the Langevin equations d dt u it = η i t, with η i t = and η i tη j t = D i δ ij δt t, where D i is the diffusion constant of the ith particle and η i s are independent Gaussian white noises. The initial positions of these particles are u i = u i such that < u < u < < u N. The process gets stopped at a random time t = t s when a specific event occurs. In this paper we consider two different mechanisms of stopping event called process [see Fig. i] and process [see Fig. ii]: In process, we consider the evolution of the N Brownian walkers till time t s when either the first particle crosses the origin for the first time before any two walkers meet each other or any two particles meet for the first time before the first particle crosses the origin [see Fig. i]. In process, t s is the time when the first particle crosses the origin for the first time before any two walkers meet [see Fig. ii].

3 3 i t s t Process L L L L t s u t ii t u t Process m m m m u t u t, u u space, u u space Fig. Color online Schematic representation of the process and process with N = particles: u t represents the trajectory of the st particle the leftmost one and u t represents the trajectory of the nd particle the rightmost one. In both cases, the trajectories of the particles are non-intersecting. In the physics literature, such nonintersecting Brownian motions are called vicious walkers [34, 35] and have been recently studied in various contexts [36,37,38,39,4,4,4]. Let p i m udm, with i =,,, N, denote the probability that m i [m, m + dm], where m i = max t ts u i t is the maximal distance travelled by the ith walker till the stopping time t s. Here we mainly focus on the PDFs p m u and p N m u because m and m N provide characterization of certain geometrical properties of the Brownian walker trajectories. For instance, one may think about m as an estimate of the common region visited by all the N walkers till the process stops given that all the particles initially started very close to the origin. Similarly, m N characterizes the number of distinct sites visited by the N walkers till t s. Recently we have studied the distributions of the number of distinct sites and common sites visited by N independent walkers over a fixed time interval [, t] [43]. Our initial motivation was to generalize this case to interacting walkers over a fixed time interval. But it is a harder problem to solve. However we show in this paper that the problem with a stopping time is solvable even in the presence of interactions. It is also interesting to note that introducing an extra random variable namely the stopping time t s renders the problem analytically tractable. Before presenting the details of our calculations, it is useful to give a summary of our results. We first study the N = particle problem because it is fully solvable even when the diffusion constants of the two particles are different i.e D D and also because the basic concepts are easy to present in this case. Solving a backward Fokker-Planck BFP equation we are able to find the full distributions p m u, u and p m u, u corresponding to the maximal displacements m and m of the right and left particle, respectively see Fig.. We show explicitly for both processes and that the PDFs p m u, u and p m u, u have power law tails valid for m u, as p m u, u A u, u, D, D m ν, 3 p m u, u A u, u, D, D, 4 m ν D 3 arctan D with exponents ν = ν =. 5 D arctan D The functions A i u, u, D, D are the amplitudes associated to the algebraic tails of the PDF p i m u, u with i =,. While these amplitudes differ from process to process, the exponents

4 4 ν = ν for the right and left walkers are process independent. The amplitudes depend explicitly on the initial positions u, u as well as on the diffusion constants D, D. Explicit expressions of A i u, u, D, D for both processes and are given in Eqs. 3 to 6. Next we consider the general N-particle problem. In this case, based on the results for the noninteracting case [Eq. ] as well as on the results of the N = -particle problem, one generally expects that the PDF p i m u of the maximal distance m i of the i th particle till the stopping time t s, has an algebraic tail: p i m u A i u, D m νi, m u N, i =,, N. 6 The exponents ν i s and the amplitudes A i s are, in general, different for the two processes for N > note that for N =, while the exponents are same, the amplitudes are different. They also depend explicitly on the number of particles N and on the diffusion constants D = D, D,, D N. Proving the result in Eq. 6 for any i =,,..., N and general N is a hard task. However, one can make some progress for i = N i.e for the maximal distance m N travelled by the rightmost walker. When the walkers are identical i.e. when they have identical diffusion constants D = D =... = D N = D, we estimate the tail of the PDF p N m u using a heuristic scaling argument based on the distribution f N t s u of the stopping time t s. This argument, for both processes and, yields : ν N = N +, when D = D =... = D N = D. 7 We also obtain an explicit expression of the prefactor A N u, D = D,, D N = D in Eq. 6. We observe that for identical diffusion constants this prefactor does not depend on D explicitly. Hence suppressing D from the argument, we denote A N u, D = D,, D N = D = A N u and show that it is given by A N u N B N Fu, where Y N u = N i= u i with Fu = Y Nu S N u, 8 i<j N u j u i, 9 and S N u is an exit probability whose value is for process and smaller than for process [given in Eq. 86]. We also present a formal exact expression of the N dependent constant B N, which for large N, is shown to grow asymptotically as [ ] N B N exp log N + olog N, where olog N represents terms smaller than log N. This large N asymptotic form of B N should be compared with the corresponding behavior b N exp [ N loglog N] in the non-interacting case in Eq.. The paper is organized as follows. In section we consider the two walkers problem where we evaluate the PDFs p m u, u and p m u, u corresponding to m and m respectively. In this section we solve a BFP equation, which under the stopping time framework becomes a Laplace s equation. From the solution of the BFP equation we find the distributions of the individual maximal distances of the first and second particles. In section 3 we consider the general N-particle problem. This section is divided into two subsections. In the first subsection 3., we give a heuristic scaling argument based on the distribution of the stopping time t s, to find the exponent ν N of the power law tail of the PDF p N m u corresponding to the global maximum m N. In the next subsection 3., we present a more rigorous calculation based on N-particle Green s function to establish the power law obtained in the previous section 3.. This calculation also provides exact expressions for the amplitudes associated to the tail of p N m u. Some technical details have been left in Appendices A, B, C and D.

5 5 Two walkers problem N = : exact solution Let us consider the motion of two non-identical Brownian walkers u t and u t given by u i t = η i t, with η i t =, for i =,, and η tη t = D δt t, η tη t = D δt t, η tη t =, where D and D are the diffusion constants of the first left and second right particle respectively. To compute the PDFs of the individual maximum displacements m and m, respectively, of the first and second particle, we start by defining the joint cumulative distribution function QL u QL, L u, u = Prob.[m L, m L < u < L ; u < u < L ], 3 given that the initial positional order is maintained till t s and L = L, L. The marginal cumulative distribution Q L u, u = Prob.[m L < u < L; u < u < ] is obtained by taking the limits L L and L in QL u whereas the marginal cumulative distribution Q L u, u = Prob.[m L < u < L; u < u < L] is obtained by taking the limit L L = L of QL u. To find QL u we consider a different problem. We consider the first exit problem of a single Brownian walker ut = u t, u t moving in two dimensions inside the region W = OBCD described in Fig. i. We are interested in the probability with which the d-walker exits from W = OBCD through specific boundaries for the first time. We denote this first exit probability by F u; L for both processes and. For process, the exit probability F u; L represents the probability that the d-walker, starting from position u, u, exits from the region W through boundary OD or OB for the first time. When the d-walker exits through OB, it corresponds, in the original two-particle picture Fig., to the first left particle meeting the second right particle before it hits the origin for the first time at t = t s while keeping m L and m L over [, t s ]. In contrast, first exit of the d-walker through OD corresponds to the first particle hitting the origin before meeting the second particle for the first time at t = t s while maintaining m L and m L over [, t s ]. On the other hand, for process the function F u; L represents the probability that the d-walker exits from the region W only through boundary OD for the first time. This exit event, in the two-particle picture, corresponds to the first particle hitting the origin for the first time before the two particles meet each other while keeping m L and m L. In the limit L and L, we get the ultimate exit probability S u, u = lim lim L L F u; L, 4 which, for process, represents the the probability that the first particle hits either the origin or the second particle ultimately. Of course this occurs with probability S u, u = in this case. On the other hand, for process, S u, u represents the probability that the first particle hits the origin for the first time before it collides with the second particle. This exit probability S u, u, in case of process, is precisely the survival probability of a lamb in the so-called lamb-lion problem where it is being chased by a single diffusing lion in the presence of a refuge. If we identify the first particle as the lamb, the second particle as the lion and the origin as the refuge [4,44] then S u, u is the probability that the lamb survives i.e. reaches the refuge before being caught by the lion. This probability is smaller than one since there is a finite probability that the lion catches the lamb i.e. the second particle hits the first particle before the later hits the origin. In particular for process, one can show that [44] u S u, u = arctan u D D [ arctan D D ], 5 which for D = D = D becomes independent of D and is given by S u, u = 4 arctanu /u. The quantity S u, u has nice interpretations in terms of the trajectories of the two walkers. It represents the volume of a set of trajectories which contains all pairs of such trajectories which, starting

6 6 from positions u, u, stay non-intersecting till t s, whereas the quantity F u; L represents the volume of a subset, which contains such pairs of non-intersecting trajectories that are constrained by m L and m L. Hence the ratio F u;l S u,u gives the fraction of such pairs of vicious trajectories which have m L and m L. This fraction precisely represents the cumulative probability QL u defined in Eq. 3. Hence, if we know the exit probability F u; L, the cumulative probability QL u is obtained from QL u = F u; L S u, u, where S u, u = lim lim L L F u; L. 6 The next question is then how to compute this exit probability F u; L in Eq. 6. In the next subsection we show that the probability F u; L satisfies a Laplace s equation which we solve with boundary conditions specified for both process and process.. Backward Fokker-Planck equation for F u; L A powerful tool to study the PDF of first passage times, like t s in our problem [see Fig. ], is the backward Fokker-Planck equation [4,5]. Here we are actually dealing with functionals of t s, m i = max{x i t, t t s }. For such functional also, it is possible to use an approach based on BFP equation see Ref. [9] for a review. Here we write down a BFP equation for the quantity F u; L treating the initial coordinates u i as independent variables. To do this, we consider trajectories u t and u t of the two Brownian particles over the interval [, t s ], which evolve according to the Langevin Eqs.. We first split the time interval [, t s ] into two parts: [, t] and [ t, t s ]. In the first infinitesimal time window [, t] the two Brownian particles will move from their initial positions {u, u } to new positions {u + u, u + u }, where u = t η t dt and u = t η t dt. 7 These two new positions are considered as new" initial positions of the two Brownian particles, respectively, for the evolution in the subsequent time interval [ t, t s ]. Since the evolution of the positions of the particles are Markovian, we have F u; L = F u + u; L. 8 u By Taylor expanding the right hand side of the above equation in u, u we have F u; L = F u; L + i= From the Langevin Eqs. one can easily show that [ F u i + F ] u u i u i + F u u i u u u i =, u i u j = δ ij D i t for i =,. Using these relations in Eq. 9 and keeping only terms of O t we obtain the following partial differential equation D F u + D F u =, with boundary conditions BCs determined by the stopping rules, which are thus different for process and process. The Eq. is valid over the region W = { u L ; u u L }. It is convenient to perform the following rescaling v i = u i Di, l i = L i Di,

7 7 u i v Process ii v Process L D C l d F= c= l, l l d F= c= l, l W B F= T F= F= T F= O /4 L u o F= β b = l, l tan β v l o F= β b = l, l tan β v l Fig. Color online i Motion of a d-walker inside the region W, ii Boundaries and boundary conditions associated to the Laplace s equation in 3 for process and process. Note that l i = L i Di and tanβ = D/D. which transforms the trapezium OBCD in the u, u plane to the trapezium obcd in the v, v plane where tanβ = D /D. Under this transformation the Fokker-Planck equation in becomes the Laplace s equation F v + F v =, 3 which holds over the region T = obcd = { v l ; v tanβ v l } [see Fig. ii] with appropriate BCs. We give the BCs in Table, which can be understood from the following arguments: Boundary conditions with tanβ = D D Boundary process process v = [od] F v =, v ; l, l = F v =, v ; l, l = v = tan βv [ob] F v, v = tanβv ; l, l = F v, v = tanβv ; l, l = v = l [bc] F v = l, v ; l, l = F v = l, v ; l, l = v = l [cd] F v, v = l ; l, l = F v, v = l ; l, l = Table Table of boundary conditions associated to the Laplace s equation in 3 for process and process. BC on the segment [od] v = : If the first particle starts with u = and u u L i.e. v = and v l, then the first particle immediately crosses the origin, which implies t s =, for both process and process. Clearly then the maximal displacement of the two particles remains m = < L and m = u L implying the BC F v =, v ; l, l = on [od] [see Fig. ii] for both processes and. BC on the segment [ob] [v = tanβv ]: If both particles start from the same position i.e. u = u < L < L, then they immediately collide implying t s = and hence, for process, F v, v = tanβ v ; l, l =. On the other hand, process excludes the possibility of any collision between the two particles even at time t = t s. This implies F v, v = tanβ v ; l, l = instead of [see Fig. ii]. BC on the segment [bc] v = l : When the initial positions are say u = L and L u L i.e. v = l and l v l, then clearly m = L at t = and it will definitely become larger than L in the next subsequent instant. Hence the BC on the segment [bc] is F v = l, v ; l, l = [see Fig. ii] for both processes and. BC on the segment [cd] v = l : When the second particle starts from its initial position u = L i.e. v = l then m = L right at the beginning and m will definitely become larger than L in the next subsequent instant implying F v, v = l ; l, l = on the segment [cd] [see Fig. ii] for both processes and.

8 8 To summarize, we finally have to solve a Laplace s equation in 3, which holds inside the polygon T in the plane v, v shown in Fig. with BCs specified in Table for both process and process.. Solution of the Laplace s equation via conformal mapping Solving the Laplace s equation in 3 for any given BC is not a priori an easy task. However using a conformal transformation of the variables, one can transform boundaries of the domain T to a much simpler geometry, while leaving the Laplace s equation itself invariant. Following Ref. [3], we here use the Schwarz-Christoffel S-C transformation which operates as follows: For a polygon P see Fig. 3 in the W plane having n vertices {w, w,..., w n } with corresponding interior angles {α, α,..., α n }, there exists a transformation W = W z from complex z-plane to W plane such that the upper half R of the z-plane gets mapped onto the interior region R of the polygon in the W plane. Under this transformation W = W z, the real axis in the z-plane gets mapped onto the boundary of the polygon P with the n vertices {w, w,..., w n } being images of the n specific points {x, x,..., x n } on the real axis. As a result, solving the Laplace s equation with complicated boundaries reduces to finding the electrostatic potential on the upper half of the complex z-plane when the potential is given on the real axis: the electrostatic potential can then be obtained explicitly from the Poisson s integral formula. The S-C transformation reads as [45] W z = B z t x α t x α... t x n αn dt + C, 4 where B and C are arbitrary constants. It is convenient to choose one point, say x n, at, such that the last factor t x n αn present in the integrand of Eq. 4 is absent. In our problem, we v W plane z plane y α 5 w α α R α 3 α 4 w 5 w 4 R w w 3 v, x x x, 3 x 4 x 5 Fig. 3 Color online Schematic representation of Schwarz-Christoffel transformation W z in Eq. 4, such that w i = W x i. x have a trapezium obcd as shown in Fig. 4 for both processes and. We chose a point b on the real line of the z-plane at, which corresponds to the image of vertex b on the W plane see Fig. 4. Moreover, let us consider that the points c, d, o on the real line with coordinates x = a, x = and x = are mapped onto the vertices c, d, o see Fig. 4 under the transformation W z. One thus has W z = B z t β + t + t + a dt + C, where β = arctan D D, 5 where a, B and C are unknown constants to be determined. Since in our case the origin is mapped onto itself under the transformation W z, i.e. W =, we have C =. Hence, from Eq. 5 and Eq. we have u + i u z dt = W z = B D D t θ + t a + t where θ = arctan D. 6 D

9 9 d v W plane c = l, l y z plane o, β +β β b = l, l tanβ v b a, c,, d o Fig. 4 Color online Schematic representation of Schwarz-Christoffel transformation in Eq. 5. x Here we note that for D = D the exponent θ = 4. The unknown constants B and a in Eq. 6 are determined as follows. The points c a, and d, on the real axis get mapped onto the points c l = L D, l = L D and d, l on the W plane respectively, which implies d = W d = c = W c = where the variable α is defined as i l α tanβ = B l + i α tanβ = B a t θ + t a + t dt, 7 t θ + t a + t dt, 8 α = L /L. 9 Simplifying Eqs. 7 and 8 one obtains the following two expressions, B = l α tanβ k θ a h θ a k θ a = α tanβ = L D, 3 L D e iβ = L k θ a D exp [ i arctan D D ], 3 with h θ a = a t θ t θ dt and k θ a = dt, 3 t a t t a t which determine the two constants a and B. The solution of the Eq. 3 gives the value of a for given α and β whereas using this solution for a in Eq. 3 we get B. When α, i.e. L L, the vertices c and b of the trapezium obcd approach to each other. This means that the point c on the real axis of the z-plane Fig. 4 should approach b implying a as L L. On the other hand, when α, i.e. L L, the point c should approach d implying a in this limit. Hence we expect that the value of a should lie in the interval [, + for < α < where both the integrals h θ a and k θ a are smooth real valued functions of a for given θ. In Fig. 5 i and ii we show how the integrals h θ a and k θ a behave as a function of a for θ = 3. Solving the Eq. 3 we obtain the value of a for given α and β and using this value of a in Eq. 3 we can get B. In Fig. 6 we plot a as a function of α for D =.3 and D =.5 and see that a diverges when α goes to whereas a approaches when α as expected from the above arguments. Let us analyze the integrals h θ a and k θ a in detail to see how a behaves as a function of α when α and α separately.

10 h θ a i. 5 5 a k θ a ii 5 5 a h θ a / k θ a..5 iii k θ 3 iv a θ Fig. 5 Color online Properties of the integrals h θ a and k θ a given in Eq. 3. In figures i, ii and iii we have used the value θ = /3. We first consider the case when α approaches from above i.e. a. Expanding the two functions entering the left hand side l.h.s. of Eq. 3 for large a we get [ ] Γ [θ] h θ a a tan β Γ [ + θ] + θ Γ [ 3 + θ] a + Oa, 33 [ ] Γ [θ] k θ a a Γ [ + θ] + θ Γ [ + θ] θ 4Γ [ 3 + θ] a + Oa, 34 where Γ [x] is the Gamma function. Using the above expressions in Eq. 3 we see that a diverges as 3 a α Fig. 6 Color online Plot of a as a function of α obtained by numerically solving Eq. 3 for a given value of β see Eq. 5 corresponding to D =.3 and D =.5. a /α as α approaches. Next we consider the case α where we expect a to approach. Expanding the functions h θ a and k θ a around a = we find h θ a + θ tanβa + O a, 35 k θ a loga + O a loga, 36

11 which with the help of Eq. 3 yields a e α tanβ. In Table we summarize the values of a and B for different α = L /L. Once the values of B and a are determined, the conformal transformation α a /α B = LΓ [ + θ] tan β ae iβ DΓ [θ] < α < < a < B = Lα tanβ Dk θ a e iβ α a e α tanβ B = L D e iβ Table Table of values of a and B for given L, α = L /L and β = arctan D D. in Eq. 6 is uniquely defined. Under this transformation 6 the Laplace s equation 3 remains invariant i.e we still have F x, y x + F x, y y =, 37 in the new variables x, y which holds over the upper half complex plane. The BCs on the real axis are i F x, = for x < and F x, = for x process 38 and ii F x, = for x and F x, = otherwise process. 39 The solution of the Laplace s equation in the upper half complex plane can be written explicitly in terms of the values at the boundary by using Poisson s integral formula F x, y = y F x, y + x x dx. 4 Using the BCs in Eqs. 38, 39 and performing the integral in both cases we get the following explicit solutions { F x, y = arctan y +x for process [ arctan y x arctan y ] 4 for process, expressed in terms of x and y. +x.3 Results and discussions In the previous section., we have solved the Laplace s equation in 3 using S-C conformal mapping which provides the first exit probability F in terms of the x, y coordinates i.e. in the z-plane. Then, to obtain the cumulative distribution QL, L u, u defined in Eq. 3, we first need to express the solution F x, y in terms of our original coordinates v = u D, v = u D which can, in principle, be done by inverting the conformal transformation W z = W in Eq. 6. Once this inversion is performed, the marginal cumulative distribution Q L u, u = Prob.[m L < u < L ; u < u < L ] is obtained by taking L L and L limits of QL, L u, u whereas the marginal cumulative distribution Q L u, u = Prob.[m L < u < L ; u < u < L ] is obtained by taking L L = L limit of QL, L u, u. The inversion of the transformation W z = W for any given L and L can not be done analytically in terms of elementary functions but one can do this inversion numerically. In the asymptotic limit L and L, as we will show below, the

12 S-C transformation gets simplified and there one can invert the transformation analytically. In section.3. we evaluate Q L u, u and Q L u, u obtained by numerical inversion of the conformal transformation. These expressions can then be compared to their numerical estimation obtained by direct simulation of Langevin equations. In section.3. we present large L asymptotics which allows to study the tails of the marginal distributions p i m u, u. Finally, in section.3.3 we study the correlations between m and m..3. Evaluation of Q L u, u and Q L u, u through numerical inversion We first compute Q L u, u, i.e. the probability that the maximum m of the st particle remains below L till the stopping of the two-particle process. It is obtained from QL, L u, u in the L limit keeping L = L fixed. This corresponds to α = L /L which, with the help of table, implies a and B = L e iβ where β is given in Eq. 5. Using this value of B D and taking the limit a of the S-C transformation in Eq. 6 we have D u + D u z L e iψ β t θ = D + t where we have written the complex coordinate W = W = u D dt, 4 + i u D on the left hand side as D u + D u e iψ D u, with ψ = arctan. 43 D D D u For given u, u and L, we numerically solve the above equation 4 for z = x + iy. Plugging this solution into Eq. 4 first and then using Eq. 6 we obtain Q L u, u. In Fig. 7 i and ii we compare the value of Q L u, u obtained from numerical inversion of the transformation W z [i.e. solving Eq. 4] to the value of Q L u, u obtained from direct simulation of the Langevin Eqs. for both process and process, and for u =.7, u = 3.5, D =.3 and D =.5. We observe a very good agreement between the analytical and numerical curves. Similarly, to evaluate Q L u, u, i.e. the probability that the maximum m of the nd particle stays below L till t s, we first take the limit L L = L i.e. α of QL, L u, u. We know from Table that, when α the coordinate a goes to. As a result, the S-C transformation in Eq. 6 now reads Γ [θ]γ [ θ] D u + D u z L e iψ β = D + D t θ + t dt. 44 For given u, u and L, we solve the above equation for z = x + iy numerically and plug the solution into Eq. 4 to obtain Q L u, u from Eq. 6. In Fig. 7 iii and iv we compare the value of Q L u, u, obtained by numerically inverting the transformation W z [i.e. solving Eq. 44] to the value of Q L u, u obtained from direct simulation of the Langevin Eqs. for both process and process, and for u =.7, u = 3.5, D =.3 and D =.5: here also one observes a very good agreement between the analytical and numerical curves..3. Large L and L limits We now focus on the limit where both L and L are large, keeping the ratio α = L /L fixed. In this limit, the S-C transformation gets simplified which makes it possible to invert the transformation W = W z analytically. For simplicity we present here our calculation assuming D = D = D. The calculation for D D can be done similarly. For D = D = D we have β = 4, θ = 4, W = u + i u u = + u e iψ u where ψ = arctan. 45 D D From Eq. 3, we see that B diverges linearly with L as k θ a is finite see Fig. 5 ii. Dividing both sides of Eq. 6 by B and taking L and L limit while keeping α = L /L fixed, we see u

13 3 Q L u, u Q L u, u Process Process. Simulation Theory i iii Simulation Theory L Q L u, u Q L u, u Simulation Theory ii iv Simulation Theory L Fig. 7 Color online Plots of marginal cumulative distributions. The open circles represent the data obtained from numerical simulations while the solid line corresponds to a numerical evaluation of our exact formula. In figures i and ii we show a plot of Q L u, u as a function of L while in figures iii and iv we show a plot of Q L u, u again as a function L for both processes. The parameters used for this plot are u =.7, u = 3.5,D =.3 and D =.5. the left hand side of Eq. 6 decreases to zero. This suggests us to expand the integral on the right hand side of Eq. 6 around z = to get W = u + iu B B D = z [ 4 4 a ] + a z + Oz 5a. 46 Following Ref. [3], we now invert the transformation W = W z from z-plane to u, u plane to obtain z Re iψ where, denoting L = L and L = αl, we have R = [X a u + u + ax a u L ] u 4 cos 4ψ 5a L 8 + OL, 47 + ax a u Ψ = 4ψ + + u sin 4ψ 5a L 4 + OL, 48 h/4 a 4 a with X a =, 49 4 and h /4 a is defined in Eq. 3. This large L expansion can in principle be carried out systematically to arbitrary order. We now take the small z limit of the explicit solutions F x, y in Eq. 4 and then inject the above large L expansion of z Re iψ into it to get the exit probability, mentioned in Eq. 6, as F u, u ; L, αl X a u + u sin 4ψ L 4 + OL 8 for process, 4ψ 3a 5a X au + u sin 4ψ L 4 + OL for process, 5

14 4 where ψ and X a are given in Eq. 45 and 49 respectively, with a implicitly determined from Eq. 3. Putting L = in the above equation 5 we get the probability S u, u defined in Eq. 4, which for process is equal to and for process is equal to 4ψ. Using the expression of ψ from Eq. 45, we get explicit expression of S u, u for process, as announced in Eq. 5 with D = D. One can follow the same calculation to get S u, u for D D. After few simplifications, one can rewrite the exit probability F u, u ; L, αl in Eq. 5 in the following form : F u, u ; L, αl S u, u cay u, u L 4 for L u, 5 for both processes and, where the constant ca is given by ca = 4X a, for process 3a 4X a 5a, for process, 5 and Y u, u = u u u u for both processes and. 53 Plugging the expression of F u, u ; L, αl from Eq. 5 into Eq. 6, we get the joint cumulative distribution QL, αl u, u ca Y u, u S u, u L 4, for L u. 54 Note that the α dependence in the above expression comes only through a since it is a function of α = L L [see Eq. 3]. Taking the limit α i.e. a see Table in the above expression, we get the marginal cumulative distribution Q L u, u [defined below Eq. 3] of the maximum m given that the st particle always stayed below the nd particle till the stopping time t s. Similarly, if we take the limit α i.e. a limit, we get Q L u, u, the cumulative distribution of the maximum m. After taking the derivative of Q i L u, u with respect to L and putting L = m we get the marginal PDFs, p i m u, u for i =, which behave like p m u, u k Y u, u m 5, for m u, 55 p m u, u k Y u, u S u, u m 5, for m u, 56 where the numerical constants k and k are obtained by taking, respectively, a and a limits of ca and finally multiplying it by 4 [coming from the derivative of L 4 w. r. t. L in 5]. The constants k and k are explicitly given by k = 3, for process 6 3, for process 8 and k = Γ [ 4 ]8, for process Γ [ 4 ]8, for process We can easily see that the tails of the PDFs p i m u, u in Eq. 56 are of the form announced in Eqs. 3 and 4 with D = D. The Vandermonde determinant in the expression of Y u, u 53 reflects the fact that the two walkers are non-intersecting and is reminiscent of the connection between vicious walkers and random matrix theory [36]. A similar calculation can be performed in the case of different diffusion constants D D to obtain QL, αl u, u Bu, u, α, µ L µ with µ = arctan 57 D D for L u 58

15 5 which finally provides the PDFs p i m u, u A iu, u, µ m ν with ν = µ + for m u. 59 The explicit expressions of the functions Bu, u, α, µ and A i u, u, µ have been left in Appendix B. Here we present the results of p m u, u and p m u, u obtained by simulating directly the Langevin equations in Eq. and compare them with the analytical prediction in Eq. 59. In Fig. Process Process p m u, u - - i Simulation Theory p m u, u - - ii Simulation Theory m m p m u, u - - iii Simulation Theory p m u, u - - iv Simulation Theory m m Fig. 8 Color online Plots of the PDFs of m, the maximum of the first left particle and of m, the maximum of the second right particle up to the stopping time t s. The open circles correspond to numerical simulations while the solid lines represent the exact asymptotic behaviors the red ones for the large argument asymptotics and the green one for the small argument asymptotics as explained in the text. There is no fitting parameters. In figures i and ii we show a plot of p m u, u as a function of m whereas, in figures iii and iv we show a plot of p m u, u as a function of m for both processes. The parameters used for these plots are u =.7, u = 3.5, D =.3 and D =.5. The red solid lines have a slope ν = 4.86 as expected from Eqs. 58 and 59 while the green solid lines have a slope. 8 i and ii we show a plot p m u, u with open circles obtained from simulation, respectively for process and process, with u =.7, u = 3.5, D =.3 and D =.5. We have also plotted the large m asymptotic behavior obtained from our analytical prediction 59 with the same set of parameters for which, one expects from Eqs. 58 and 59 that ν 4.86 and, from Eqs. 3 and 4, A u, u, µ for process and A u, u, µ 5.3 for process. We see that the agreement between our numerical simulations and our analytical results is very good. Notice also that for m u, the first particle does not feel the presence of the second particle and therefore one expects that in this limit p m u, u u /m for process and p m u, u u /[S u, u m ] for process, with S u, u.549 in the later case [see Eq. 5]. These asymptotic behaviors for small m are shown as green line in Fig. 8 i and ii: the agreement with the numerical data is rather good. Finally, in Fig. 8 iii and iv we show a comparison between p m u, u evaluated numerically open circles and the analytical predictions for the tails 59 [see also Eqs. 5 and 6]. Here also the agreement is very good.

16 H β α α Fig. 9 Color online Plot of H β α, which gives a measure of correlation of the two maxima m and m for D =.3 and D = Correlation between the maxima m and m We end up this section by considering the correlations between the maximal displacements m and m of the st and nd particle, respectively. To characterize these correlations, we define the following quantity C θ α, L; u, u = [QL, αl u, u Q L u, u Q αl u, u ]L θ, 6 with α = L /L and θ is given in Eq. 6. This quantity measures the difference between the joint cumulative probability QL, αl u, u of m, m and the product of their individual marginal cumulative probabilities Q L u, u and Q αl u, u. If the two maxima m and m are independent of each other then the quantity defined above would be identically zero for any α and L. We plug the large L expression of QL, αl u, u from Eq. 58 and the large L of Q L u, u and Q αl u, u obtained from Eq. 58 by taking α and α limit, respectively, to see that the function C θ α, L; u, u becomes independent of L [the factor L /θ in 6 is chosen for this purpose] and takes the following form C θ α, L; u, u L = H β αl β u, u, 6 for both processes, where the function H β α carries the information on the correlations the function L β u, u can also be computed explicitly but we do not discuss it here. This function H β α is given by H β α = [ tanβ θ + α θ Γ [θ] θ kθ a a Γ [ + θ] α ] θ, with θ = β and β = arctan D. 6 D Here we should keep in mind that according to Eq. 3, a is a function of α and β. In Fig. 9, we plot H β α as a function of α for D =.3 and D =.5. This shows that even when both m, m are large, m and m are strongly correlated as long as they are of the same order of magnitude. As expected, these correlations vanish when the two particles are very far away from each other. 3 Multi-particle problem: N > In this section we generalize the two vicious walkers problem to N vicious walkers problem. We focus on the PDF p N m u of the global maximum m N maximal distance travelled by the rightmost particle.

17 7 Moreover, we assume that the N walkers are identical i..e. they have the same diffusion constant D = D =... = D N = D. In this case, we expect that the PDF p N m u will not depend on D and will have the following power law tail p N m u A N u m ν N with ν N = N +, 63 for both processes and. We first show from a heuristic scaling argument that one can predict the power law in the above equation. The scaling argument is based on the large time tail of the distribution f N t s u of the stopping time t s. This argument provides the N-dependence of the exponent ν N accurately but it does not predict the prefactor precisely. For this we study the N-particle problem rigorously in the next subsection, where we follow an approach different from what we have done for the N = case. In particular, we have used the Green s function approach directly rather than solving a N-dimensional Laplace s equation inside an N-dimensional complicated Weyl chamber because the later approach becomes difficult as we do not have at our disposal any generalized Schwarz-Christoffel transformation valid in dimensions d >. 3. A heuristic argument for N particles To justify the power law for p N m u given in Eq. 63, we present a simple scaling argument which is based on the power law tail of the PDF of the stopping time t s itself. This argument is valid for both process and process and it goes as follows. Let Q N L u = Prob.[m N L u] be the probability that the global maximum m N stays below the level L given that the N walkers starting from positions u stay non-intersecting till time t s. For large L, the dominant contribution to Q N L u = Prob.[m N L u] comes from trajectories which typically have large global maximum m N. On the other hand, using the connection between non-intersecting Brownian motions and random matrix theory RMT one can argue that the average value of the global maximum m N over the time interval [, t s ] grows as m N t s N, whereas the fluctuations of m N around this mean value decays as N /6 [36,39]. As a result, the distribution of m N over the time interval [, t s ] will be highly peaked around m N t s N for very large N. Therefore one expects that for large t s and N, the random variable m N will be typically of the order of t s N. Hence, for large N and L, the tail of the cumulative distribution Q N L u is obtained from : Q N L u = Prob.[m N > L u] Prob. with, Prob.[t s > t u] = t ] [t s > c L u, 64 N f N t udt, 65 where f N t s u is the PDF of the stopping time t s and c is an undetermined constant, irrelevant for the present argument. To find the PDF of the stopping time f N t u for identical walkers i.e. for D = D =... = D N = D, we start with the Green s function G N y, t; u, of N non-intersecting Brownian walkers with an absorbing wall at the origin. This Green s function represents the probability density of the positions y = y, y,..., y N of the N walkers at time t given that they had started from positions u = u, u,..., u N initially. Using the Karlin-McGregor formula [46], this N-particle Green s function can be expressed as the determinant of a N N matrix [G i,j ] [gu i, y j, t] where gu, y, t = 4Dt exp [ y u 4Dt ] exp [ y + u 4Dt is the single particle Green s function with an absorbing wall at the origin. One can also map the problem of finding G N y, t; u, to the problem of finding the wave function of N free fermions with an infinite wall at the origin and this mapping allows us to write [36] G N y, t u, = Dt N N! N... dq dq... dq N Φ N q ] u Dt Φ N q 66 y Dt e q, 67

18 8 where q {q, q,..., q N }, q = N Φ N q u = i= q i and sin q u sin q u sin q u N sin q u sin q u sin q u N sin q N u sin q N u sin q N u N In case of process, the motion of the N walkers stop when either any two particles meet each other for the first time before the leftmost particle hits the origin or the first particle crosses the origin for the first time before any two particles meet each other. The survival probability S N t u of such N-particle process, is given by S N t u = d N y G N y, t u,. 69 W The Weyl chamber W in the above expression is defined as W = [y R N + y y y N ] where R + is the set of non-negative real numbers. Hence, for process, the PDF of the stopping time is f N t u = S N t u = dy dy... dy N G N y, t u, for process. 7 t t W On the other hand, for process the reasoning is a bit different as the process gets stopped only when the first particle hits the origin for the first time before any two other particles collide. This implies that the PDF f N t x of the stopping time is obtained from the outward flux through the y = hyperplane W = [y W y = ] of the Weyl chamber W. Hence integrating the outward probability current density D GN y,t u, y y = over the hyperplane W we get GN y, t u, f N t u = D dy dy 3... dy N W y y = for process, 7 where the Green s function is explicitly given in Eq. 67. To obtain the tail of the PDF f N t u, we need to find the large t behavior of G N y, t u,, which can be obtained by expanding the function Φ N q u t in Eq. 68 for large t and finite u. One can show [37] that for large t and finite u, [ ] u Φ N q = det sin u j q i γ N Y N q Y N u Dt i,j N Dt N + O t N +/,7 Dt where, γ N = NN N i! i= and Y N u = N i= u i i<j N u j u i. 73 Plugging the above large t approximation 7 into Eq. 67 and performing the integrations over the variables q i, we get G N y, t u, Y Nu N ] Dt N [ N Dt exp y y Y N i! N 4Dt Dt. 74 i= Finally, putting this large t form of G N y, t u, into Eqs. 7 and 7 and performing the rest of the integrations over the variables y i, we obtain Y N u f N t s u Dδ N Dt s ζ, with ζ = N +, 75

19 9 for both process and process, where δ N is an N dependent constant different for process and process. The explicit expressions of δ N for both processes are given in Appendix A. The above result for f N t s u has also been proved in [47,48] for process. Plugging the large t s behavior of f N t s u from Eq. 75 into Eqs. 64 and 65 we get, Q N L u B N Y N u/l N, for large L 76 where, the function Y N u is given in Eq. 73 and B N is an N dependent constant. Upon deriving both sides of the Eq. 76 with respect to L at L = m, we get the tail of the PDF p N m u as given in Eq. 63 with ν N = N + and A N = N B N Y N u. This heuristic argument also provides a rough estimate of B N for large N and that is B N exp N log N. Similar scaling arguments have been successfully used to study the distribution of the global maximum m N of N non-interacting particles till their first exit from the half space [33]. This result 76 is in line with the following general results valid for a generic self-affine process. For such a process xt starting from x >, the cumulative distribution QL x = Prob.m L x of the maximum m till the stopping time t s time of first passage through x =, or equivalently the exit probability QL x from the box [, L] through the origin, has been recently studied in [8] where it was shown that QL x x/l φ in the x/l limit. The exponent φ is related to the persistence exponent θ p and the Hurst exponent H via the scaling relation φ = θ p /H [8]. The persistence exponent θ p characterizes the late time power-law decay of the survival probability, i.e. the probability that the process stays on the positive half-axis up to time t [5,7], whereas the Hurst exponent characterizes the typical growth of xt t H with time t. Thus, the PDF of the maximum decays for large m as P m x m φ with φ = θ p /H. From Eq. 75 we see that for the multiparticle process the corresponding persistence exponent is θ p = N. If we consider this N-particle process as a single self-affine process in N-dimensional space with H = /, then the general argument from [8] suggests that QL x /L φ, with φ = θ p /H = N, which is in accordance with Eq. 76, although this argument can not predict the precise dependence on the initial positions. In the next subsection we prove ν N = N + on firmer grounds and compute the amplitude exactly. 3. The distribution p N m u of the global maximum for N > Here we study the distribution p N m u of the global maximum m N of N identical i.e. with identical diffusion constant D = D =... = D N = D vicious walkers using the N-particle Green s function. We first compute the cumulative probability Q N L u = Prob.m N L u which represents the probability that the global maximum m N of the rightmost particle is less or equal to L given that the walkers, starting from positions u = u, u,..., u N stayed non-intersecting till the stopping time t s. Upon taking the derivative of Q N L u with respect to L at L = m we get the PDF p N m u. To compute this cumulative probability Q N L u we consider the first exit problem of a single N- dimensional Brownian walker ut = u t, u t,..., u N t from the region T N L = {u R N + < u < u <... < u N < L}, as done for the N = -particle case see Fig.. For process we consider the first exit probability of the walker through any of the boundaries u = or u i+ = u i with i =,,..., N. These exit events correspond, in the original N-particle problem, to the following events: a leftmost particle crossing the origin for the first time before any two particles meet or b any two particle meet for the first time before the leftmost particle hits the origin. On the other hand for process, we consider the first exit probability of the N-dimensional walker only through the boundary u =. This event corresponds to the leftmost particle hitting the origin for the first time before any two particles meet in the N-particle picture. We denote this first exit probability for both process and process by F N u, L. For process, one can see that F N u, L is equal to the time integration from t = to t = of the total outward probability flux through all the boundaries of T N L which is equal to minus time integration of the outward flux through the boundary u N = L whereas for process F N u, L is equal to the time integration of the outward flux only through the boundary u =. Hence, the probability F N u, L can be expressed in terms of the Green s function