arxiv: v1 [cond-mat.stat-mech] 15 Mar 2019

Transcription

1 arxiv: v1 [cond-mat.stat-mech] 15 Mar 19 Distribution of Brownian coincidences Alexandre Krajenbrink Laboratoire de Physique de l Ecole Normale Supérieure, PSL University, CNS, Sorbonne Universités, 4 rue Lhomond, 7531 Paris Cedex 5, France Bertrand Lacroix-A-Chez-Toine LPTMS, CNS, Univ. Paris-Sud, Université Paris-Saclay, 9145 Orsay, France Pierre Le Doussal Laboratoire de Physique de l Ecole Normale Supérieure, PSL University CNS, Sorbonne Universités, 4 rue Lhomond, 7531 Paris Cedex 5, France Abstract. We study the probability distribution, P N (T ), of the coincidence time T, i.e. the total local time of all pairwise coincidences of N independent Brownian walkers. We consider in details two geometries: Brownian motions all starting from, and Brownian bridges. Using a Feynman-Kač representation for the moment generating function of this coincidence time, we map this problem onto some observables in three related models (i) the propagator of the Lieb Liniger model of quantum particles with pairwise delta function interactions (ii) the moments of the partition function of a directed polymer in a random medium (iii) the exponential moments of the solution of the Kardar-Parisi-Zhang equation. Using these mappings, we obtain closed formulae for the probability distribution of the coincidence time, its tails and some of its moments. Its asymptotics at large and small coincidence time are also obtained for arbitrary fixed endpoints. The universal large T tail, P N (T ) exp( 3T /(N 3 N)) is obtained, and is independent of the geometry. We investigate the large deviations in the limit of a large number of walkers through a Coulomb gas approach. Some of our analytical results are compared with numerical simulations.

2 CONTENTS Contents 1 Introduction and main results Introduction Summary of the main results Coincidence time for Brownian bridges and KPZ equation with droplet initial conditions 8.1 Bethe ansatz solution of the Lieb-Liniger model for droplet initial conditions Probability Distribution Function of the coincidence time Full distribution of the coincidence time for N =, 3 Brownian bridges 1.4 Asymptotic behaviours of P N,BB (T ) for arbitrary N Small T limit of the Probability Distribution Function Large T limit of the Probability Distribution Function Mean and variance of the coincidence time for N Brownian bridges Coulomb gas method for N Coincidence time for Brownian motions and KPZ equation with flat initial condition Bethe ansatz solution of the Lieb-Liniger model for half-flat initial conditions Flat initial conditions as a limit of the half-flat Moment generating function and probability distribution function for N =, Moment generating function for arbitrary N Probability Distribution Function of the coincidence time Asymptotic behaviour of P N,BM (T ) for arbitrary N Small T limit of the Probability Distribution Function Large T limit of the Probability Distribution Function Mean and variance of the coincidence time for N Brownian motions. 5 4 Coincidence time for arbitrary fixed final points Exact distribution for N =, Small T asymptotics of the PDF for any N and for any fixed final points Large T asymptotics of the PDF for any N and arbitrary final and initial points Mean coincidence time for arbitrary final points Conclusion 9 A Expected value of the coincidence time for non identical diffusing particles 31 B Details of the calculations of the distribution for N = 3 Brownian bridges 31 C Details of the numerical simulations 3

3 CONTENTS 3 1. Introduction and main results 1.1. Introduction In this article, we study independent and identical diffusive Brownian processes. We want to characterise in details the statistics of the time spent by identical diffusing particles in the vicinity of each other, which we call the coincidence time. This problem is relevant for instance to analyse reaction networks of the type A + A B. (1) Indeed, for chemical species to react, they first have to encounter each other and then overcome the potential barriers. This takes a finite amount of time and this coincidence time plays a crucial role in the kinetics of the reaction. The behaviour of this coincidence time is highly dependent on the spatial dimension, as it is clear that the higher the dimension, the harder it is for diffusing particles to encounter one another. We limit our study of this problem to the case of dimension d = 1, where the number of encounters is maximum. We model the identical diffusing particles via independent one dimensional Brownian motions x i (τ) on the time interval τ [, t] ẋ i (τ) = η i (τ) with η i (τ) = η i (τ)η j (τ ) = Dδ i,j δ(τ τ ). () where below we choose units such that D = 1. We consider the case where the N diffusing particles are emitted at t = from a single point-like source x i () = x for i = 1,..., N. We define the total time spent by these particles in a close vicinity of length l to each other as T N (l; t) = N i j t [ ] l Θ x i(τ) x j (τ) dτ, (3) where Θ(x) is the Heaviside step-function. We are particularly interested in the limit where the length l is small compared to all the other scales of the problem and therefore define in the limit l, 1 T N (t) = lim l l T N(l; t) = N i j t δ[x i (τ) x j (τ)]dτ. (4) We will refer to this observable as the coincidence time of N diffusing particles, i.e. the amount of time that these independent particles spend crossing each other. Note however that T N (t) does not have the dimension of a time, which is a usual property of the local time of a stochastic process [1,, 3]. Using the Brownian scaling, the process has a trivial rescaling to the time interval τ [, 1], and we obtain the equality T N (t = 1) = T N = T N (t) t. (5) We will see in the following that the N-dependence is however highly non-trivial. The case N = can be solved quite easily by considering the process z(τ) = for calculational convenience below, each pair appears twice in the sum in our definition (4).

4 CONTENTS 4 1 (x 1 (τ) x (τ)) which is also a simple diffusive process of same diffusion coefficient D = 1, ż(τ) = ξ(τ), with ξ(τ) = ξ(τ)ξ(τ ) = Dδ(τ τ ). (6) The coincidence time T of the two Brownian particles can be expressed in terms of the local time L t (x) of the process z(τ) defined as Setting x =, we obtain L t () = t = L t (x) = t dτδ(z(τ)) = t dτδ(z(τ) x). (7) t ( ) 1 dτδ (x 1 (τ) x (τ)) dτδ (x 1 (τ) x (τ)) = T (t). (8) The joint PDF of the local time L t () = L and the final position z(t) = x f obtained by Borodin et al. [] and exploited by Pitman [3] and reads was P joint (L, x f ) = 1 1 Dπt 3 ( x f + DL)e ( x f +DL) 4Dt. (9) We set t = 1 and remind that D = 1 in our convention. Using now the identities L 1 () = T / together with z(1) = (x 1 (1) x (1))/, we obtain the joint PDF of the coincidence time T = T and the final algebraic distance between the diffusive particles d = x 1 (1) x (1) as P joint (T, d) = 1 ( ) T P joint, d = ( d + T ) ( d +T ) 8π e 8. (1) However, as soon as N = 3, one cannot define independent processes for which our coincidence time is a simple observable. Here we will study the Probability Distribution Function (PDF) of T N for two cases, Brownian motions and Brownian bridges. In the first case (Brownian motion) one defines the Moment Generating Function (MGF) e ct N (t) x, where... x denotes the expectation value with respect to the PDF of T N (t) for given initial conditions x(τ = ) = x = (x,..., x ). This MGF can be expressed in the Feynman- Kač framework as an N-dimensional path integral e ct N (t) x = dyz N (y, t x ; c), with, (11) N t N ẋ i (τ) N + c δ[x i (τ) x j (τ)] dτ x(t)=y 4 i=1 i<j Z N (y, t x ; c) = Dx(τ) e, x()=x (1) where x(τ) = (x 1 (τ), x (τ),, x N (τ)). In this expression we wrote the sum i j as i<j to represent the coincidence as a delta interaction between each pair of

5 CONTENTS 5 Brownian motions. As the final point y is arbitrary, we integrated here over all the possible realisations. This situation corresponds to the right panel in Fig. 1. In the second case (Brownian bridges) we define the MGF as the following expectation value with fixed initial and final positions e ct N (t) x,y = Z N(y, t x ; c) Z N (y, t x ; c = ) (13) where Z N (y, t x ; c = ) = (4πt) N/ e (y x) /(4t) is the free propagator. This situation corresponds to the left panel in Fig Figure 1. Plot of a simulation of N = 3 diffusive particles for a diffusion coefficient D = 1 starting from x = (,, ). Left: Brownian bridges with endpoints in y = (,, ). ight: Brownian motions with arbitrary final points. The function Z N (y, t x ; c) can be interpreted as the imaginary time N-body propagator of the bosonic Lieb-Liniger Hamiltonian [4] N N t Z N = ĤN(c)Z N with ĤN (c) = ˆp i + c δ(ˆx i ˆx j ). (14) We therefore see that this MGF, eq. (13), for non-interacting diffusing classical particles is mapped onto the propagator of a quantum problem of identical bosons with a contact interaction. This problem is also related to other well-studied models: the directed polymer in a random potential, equivalent to the Kardar-Parisi-Zhang (KPZ) stochastic growth equation. Indeed, consider an elastic polymer of length t and fixed endpoints x and y, in a random potential ξ(x, t). Its partition function reads t [ẋ(τ) ] x(t)=y + ξ(x(τ), τ) dτ Z(y, x, t) = 4. (15) x()=x Dx(τ)e One may then compute the N-point correlation function of the partition function, averaged with respect to the random potential (denoted as ) and obtain N Z(y i, x, t) = i=1 x(t)=y Dx(τ)e x()=x t i=1 N [ẋi (τ) i=1 4 i<j + ξ(x i (τ), τ) ] dτ. (16)

6 CONTENTS 6 Taking now ξ(x, t) as a Gaussian white noise of variance c >, ξ(x, t) = and ξ(x, t)ξ(x, t ) = cδ(x x )δ(t t ), (17) the replica average in Eq. (16) reduces to the Lieb-Liniger real time propagator in Eq. (1) [5, 6, 7] with the value of the parameter c = c, and satisfies the same evolution equation (14). The directed polymer partition function Z(y, x, t) is also equal to the droplet solution of the KPZ equation, through a Cole-Hopf transformation, log Z(y, x, t) = h(y, x, t), where t h = yh + ( y h) + c ξ(y, t). (18) Both the directed polymer and the KPZ equation correspond to the choice c = c <, i.e. to the attractive Lieb-Liniger model. The moments (16) of the directed polymer problem, which correspond to exponential moments for the KPZ height field, have been studied extensively recently and exact formulae have been obtained for the droplet solution at arbitrary time [7, 8, 9]. The MGF of our diffusion problem, for N independent Brownian bridges, and negative value of the parameter c = c is thus related to these moments as follows e ct N (t), = Z(,, t)n Z (,, t) N (19) where Z (x,, t) = 1 4πt e x /(4t) is the free Brownian propagator (for c = ). Explicit expressions of these moments were derived from the Bethe ansatz. Note that all initial and final points here are set to. In this paper, we obtain complementary exact formula for the MGF for positive value c >, i.e. for the Laplace transform of the coincidence time e ct N (t). We study in particular the case of N independent Brownian bridges denoted e ct N (t) BB e ct N (t),, and the case of N Brownian motions starting from the same point, x =, with free final points, denoted e ct N (t) BM e ct N (t). We also use the formula (19) for negative c = c < to obtain the large T asymptotics of the PDF of T N (t). For that particular asymptotic behavior we obtain some more results for arbitrary fixed final points, using a formula for e ct N (t),y. We now summarize our main results. 1.. Summary of the main results Brownian bridges We obtain an exact formula for the PDF of the rescaled coincidence time T N for N Brownian bridges as ( P N,BB (T ) = N+1 T dr Θ T ) ( ) r l det e (r i r j ) 4 () N 1 i,j N + l The constraint i j in (1) originates from the Itô prescription in the (Feynman-Kač) stochastic heat equation satistied by Z, equivalently in the definition of the path integral.

7 CONTENTS 7 For N =, 3, we obtain explicit formulae for P N,BB (T ) respectively in Eq. (4) and (46). For general N, we extract the small and large T asymptotic behaviours as N(N 1) G(N + ) Γ(N(N 1)) T N(N 1) 1 + O(T N(N 1)+1 ), T, P N,BB (T ) = N! N 1 π N 1 N 3/ α N N H N 1 ( α N T ) e α N T + O(e β N T ), T, (1) where G(n) = n k=1 k! is the Barnes-G function, Γ is the usual Gamma function, H p (x) = e x ( x ) p e x is the Hermite polynomial of degree p and the exponential factors are 3 α N = (N 1)N(N + 1), β 3 N = (N )(N 1)N. () We also obtain the mean coincidence time in Eq. (56) and the variance in Eq. (6). Brownian motions We obtain an exact formula for the MGF of the rescaled coincidence time T N for N Brownian motions (i.e. with free final points), which depend on the parity of N. It is given in Eq. (95) for even values of N and in Eq. (96) for odd values. For N =, 3, we obtain explicit formulae for the PDF P N,BM (T ), P,BM (T ) = π e T, P 3,BM (T ) = π e T 3 (e For arbitrary values of N, we obtained the expression of the PDF as: N odd P N,BM (T ) = ( ) N dr Θ T r l N N+1 T N + l=1 1 k<l N sign(r k r l ) 8 T 1). (3) [ ( )] rk r l Pf erf 1 k,l N 1 (4) N even P N,BM (T ) = ( ) N dr Θ T r l N+1 T N + l=1 1 k<l N sign(r k r l ) [ ( )] rk r l (5) Pf erf 1 k,l N We further extracted the small and large T behaviours as I N T N(N 1) 1 + O(T N(N 1) +1 ), T, P N,BM (T ) = N 1 3 πn(n 1) e α N T + O(e β N T ), T, (6) where I N is given in Eq. (135). The coefficients α N and β N are the same as for the Brownian bridge in Eq. (). We also obtain the mean coincidence time in Eq. (119) and the variance in Eq. (11).

8 CONTENTS 8 Finally, we extended our study to N Brownian walkers with arbitrary fixed endpoints. The small T asymptotics of the PDF for arbitrary final points y (all initial points being at ) is given in (133). The large T asymptotics of the PDF for any choice of initial and final points is displayed in (139). A trivial consequence of our work is the PDF PN c (T ) of the coincidence time for N Brownian walkers interacting pairwise via a delta interaction of strength c (of any sign). It is given in terms of the PDF for non-interacting Brownian P N (T ) calculated here, simply as PN c e ct P N (T ) (T ) = + (7) dt e ct P N (T ) The remaining of the paper is organised as follows. Section, is dedicated to the results for the case of Brownian bridges. In subsection.1, we describe how the Bethe ansatz applied to the Lieb-Liniger model allows to obtain a formula for the MGF.In subsection., we derive the exact formula in Eq. () for the PDF of the coincidence time T N = T for arbitrary N. In subsection.3, we analyse this formula for N =, 3. In subsection.4, we derive the asymptotic behaviours for T and T for arbitrary N. In subsection.5, we obtain the mean and variance for arbitrary N. Finally, in subsection.6, we analyse the large N limit via Coulomb gas techniques. In Section 3, we extend all previous results - except the Coulomb gas - to the Brownian motions in the same fashion and derive the exact PDF for arbitrary N Brownian motions given in Eqs. (5) and (4). In Section 4, we study some properties of the coincidence time for N Brownian with arbitrary fixed endpoints. We obtain the exact expression of its PDF for N = and N = 3, the small T asymptotics for any N and the large T asymptotics for arbitrary initial and final points. The latter is obtained using the properties of the ground state of the Lieb-Liniger model.. Coincidence time for Brownian bridges and KPZ equation with droplet initial conditions.1. Bethe ansatz solution of the Lieb-Liniger model for droplet initial conditions The Lieb-Liniger Hamiltonian in Eq. (14) is an exactly solvable model using the Bethe ansatz. Here we need only the symmetric eigenstates of ĤN (c), which in the repulsive case c >, are single particle states. The eigen-energies E k of this system defined on a circle of radius L form, in the large L limit, a continuum indexed by N real momenta {k i } s such that Ĥ N (c) Ψ c k = E k Ψ c k, with E k = k = x Ψ c k = Cc k N! σ S N i<j Ck c = (k i k j ) (k i<j i k j ) + c N ki, (8) i=1 ( 1 ic sign(x j x i ) k σ(j) k σ(i) 1/ ) e i N j=1 xjk σ(j), (9). (3) where the Ψ c k denote the eigenstates of ĤN(c). The calculation is similar to the one in refs. [7, 8] except that we retain here for c > only the particle states. Consider

9 CONTENTS 9 the imaginary time propagator at coinciding endpoints Z N (, t ; c) = e ĤN (c)t dk = (π) N Ψc k Ψ c k e Ekt. (31) We now use that Ψ c k = Cc k and we arrive at the following Laplace transform for the PDF of the coincidence time for N Brownian bridges e ct N (t) BB = (4πt) N dk (π) N e k t (k i k j ) (k N i<j i k j ) + c. (3) where we have divided by the result at c = which coincides with the free propagator Z N (, t ; ) = (4πt) N. We note that there is an alternative formula for the moments of the directed polymer problem, which was derived using Macdonald processes [9, 1]. Using this formula we can write the MGF of N Brownian bridges equivalently as e ct N (t) BB = (4πt) N dz z i z j i (iπ) N etz z N i<j i z j + c. (33) Note that this formula is valid only for c >, otherwise the contours are different. Using the symmetrization identity z σ(i) z σ(j) + c = n! (34) z σ(i) z σ(j) σ S n i>j one can show that this formula is in agreement with (3)... Probability Distribution Function of the coincidence time In this section we obtain from Eq. (3) an expression for the PDF of the coincidence time for the N Brownian bridges. First, we use the Cauchy identity i<j N (k i k j ) ( ) (k i k j ) + c = det c 1 l,m N c + i(k l k m ) = σ S N sign(σ) We then introduce the auxiliary variables r = (r 1,, r N ) such that i<j (k i k j ) (k i k j ) + c = cn sign(σ) σ S N = c N sign(σ) σ S N N + N + dr e c l r l dr e c l r l N l=1 N l=1 c c + i(k l k σ(l) ). e ir l(k l k σ(l) ) (35) (36) N e ik l(r l r σ(l) ). (37) eplacing in Eq. (3), we are now able to compute the integrals over k, yielding e ct N (t) BB = c N = c N N + N + l=1 dr e c N l r l sign(σ) σ S N l=1 ( dr e c l r l 4t det 1 i,j N e (r i r j ) e (r l r σ(l) ) 4t (38) ) (39)

10 CONTENTS 1 Inverting the Laplace transform of this expression, we obtain the PDF P N,BB (T ) of the rescaled coincidence time T N = T N (t = 1) as ( P N,BB (T ) = T N dr δ T ) ( ) r l det e (r i r j ) 4 N 1 i,j N + l ( = N+1 T dr Θ T ) ( ) r l det e (r i r j ) 4 (4) N 1 i,j N + l where Θ(x) is the Heaviside step function. We start by analysing the case of N =, 3, where the PDF can be computed explicitly before considering the general behaviour for arbitrary values of N. This formula (4) is the most compact general expression that we could find. We show in the next section that it can be used to derive explicitly the PDF for small values of N =, Full distribution of the coincidence time for N =, 3 Brownian bridges Distribution for N =. In the case of N =, we will see how to extract the full PDF from Eq. (4). Setting N =, we may compute exactly the determinant in the integrand, yielding T T r1 [ ] T T r1 P,BB (T ) = T 3 dr 1 dr 1 e (r 1 r ) = T 3 dr 1 dr e (r 1 r ). (41) The first term in the integrand vanishes when deriving with respect to T. The PDF can then be derived in a few steps as P,BB (T ) = T T (r T ) dre = T T The asymptotic behaviours of this PDF are simple to extract as T + O(T 3 ), T, P,BB (T ) = T e T, T. due u = T e T. (4) In Fig., we compare our analytical formula of Eq. (4) with the numerical simulations of Brownian bridges (see Appendix C for the details of the simulations), showing an excellent agreement. Note that in the case N =, this result can be obtained directly from Eq. (1) by setting the final distance between the diffusive particles d = and ensuring the normalisation of the probability P,BB (T ) = (43) P joint (T, d = ) P joint (T, d = )dt = T e T. (44) As noted previously this method cannot be generalised to larger values of N, a contrario with our exact formula in Eq. (4). Distribution for N = 3. Starting from N = 3, the distribution cannot be obtained from a simple rescaling of the local time of a single Brownian motion,and we

11 CONTENTS 11.6 numerical data P, BB (T) T numerical data P, BB (T) T Figure. Comparison between the PDF of the coincidence time obtained numerically for N = Brownian bridges on the time interval τ [, 1] and diffusion coefficient D = 1 and the analytical result in Eq. (4) (some details on the simulations are provided in Appendix C). Left: Linear scale, ight: Logarithmic scale. need to use our results in Eq. (4) for the Brownian bridge. In the case of N = 3, the determinant appearing in the integrand can still be computed quite easily, yielding T T r1 [ ] T r1 r 3 P 3,BB (T ) = T 4 dr 1 dr dr 3 1 3e (r 1 r ) + e (r i r i+1 ) 4, where we used the symmetry between the r i s. Appendix B for details), we finally obtain the PDF P 3,BB (T ) = T e T 1 π + 8 Its asymptotic behaviours read T 5 3 e T i=1 (45) After some computations (See 8 (T 4) erf ( ) 3 8 T. (46) 8 + O(T 7 ), T, P 3,BB (T ) = (47) 1 π ( T 4 ) e T 8 + O(e T ), T. 8 3 In Fig. 3, we compare the analytical formula of Eq. (46) with the numerical simulations of Brownian bridges (see Appendix C for the details of the simulations), showing an excellent agreement..4. Asymptotic behaviours of P N,BB (T ) for arbitrary N We now come back to the case of N Brownian bridges and analyse the tails of the PDF P N,BB (T ) of the rescaled coincidence time T N = T N (t)/ t respectively for T and T Small T limit of the Probability Distribution Function To obtain the small T limit of P N,BB (T ), it is convenient to consider the c + limit of e ct N BB,

12 CONTENTS 1.35 numerical data.3 P 3, BB (T) T numerical data P 3, BB (T) T (T 4)e T T Figure 3. Comparison between the coincidence time obtained numerically for N = 3 Brownian bridges on the time interval τ [, 1] and diffusion coefficient D = 1 and the analytical result in Eq. (46). Left: Linear scale, ight: Logarithmic scale. which is given in Eq. (3). The trajectories which contribute to a small coincidence time are repelled from each others. In the Lieb-Liniger picture, this is consistent with the repulsive case where the states are described as single particle states rather than bound state. Taking the large c + limit of Eq. (3), we obtain e ct N BB 1 c N(N 1) dke k π N N i<j (k i k j ) = N(N 1) G(N + ) c N(N 1), (48) where G(n) = n k=1 k! is the Barnes G function. We have used the Mehta integral formula e.g. (1.5)-(1.6) in [11]. Inverting the Laplace transform, we obtain the small T behavior P N,BB (T ) = N(N 1) G(N + ) T Γ(N(N 1)) T N(N 1) 1 + O(T N(N 1)+1 ). (49) Note that setting N =, 3, using G(4) =, G(5) = 1 and Γ(6) = 1, one recovers explicitly the result of the first line of Eqs. (43) and (47)..4.. Large T limit of the Probability Distribution Function To study the large T limit of P N,BB (T ), one needs to investigate the c limit of the expectation e ct N BB which is dominated by large values of T N. The trajectories which contribute to a large coincidence time are attracted to each others. In the Lieb-Liniger picture, this corresponds to the attractive case, where particles form bound states called strings. Upon increasing c = c, the potential becomes more and more attractive and the trajectories become dominated by the configuration where all the bosons are bounded into a single string. Mathematically, for c = c < the expression of the moments of the directed polymer partition function is more involved as it involves a sum over string states [7, 8, 9, 1]. Nonetheless, in the limit c +, the energy spectrum of the Lieb- Liniger model is dominated by its ground state which is a single string containing all particles with energy E (N) = c N(N 1) 1. It follows from the contribution of the

13 CONTENTS 13 ground state to the moment (see e.g. (65-66) in Supp. Mat. of [13] replacing t by c t and multiplying by c N ) that e ct N BB = N 1 N!(4π) N 3/ N(N N 1 c e c 1) 1 [1 + O(e 1 4 N(N 1) c )] (5) One can check that this behavior is consistent with the following T asymptotics P N,BB (T ) = N!N 1 π N 1 αn ( N 3/ T ) N 1 e α N T + O(e β N T ) (51) = N!N 1 π N 1 N 3/ α N N H N 1 ( α N T ) e α N T + O(e β N T ), (5) where H p (x) = e x ( x ) p e x is the Hermite polynomial of degree p. The exponential factors are equal to 3 α N = N(N 1) 3 (53) β N = N 3 3N + N This is checked by calculating the Laplace transform of (51) using a saddle point approximation (with a saddle point for T = cα 1 N > for c = c > ). For completeness, we provide this asymptotics explicitly for N =, 3. N =, P,BB (T ) = T e T. This matches eq. (4). N = 3, P 3,BB (T ) = 1 8 π 3 e T 8 ( T 4 ) + O(e T ). This matches eq. (47)..5. Mean and variance of the coincidence time for N Brownian bridges In this section, we compute the two first cumulant of the Brownian bridge coincidence time T N for arbitrary values of N. Mean value of T N,BB To compute the first moment T N,BB, there are two alternative methods. We may obtain it by expanding Eq. (33) for small c or we may use the Brownian bridge propagator to compute it directly. In this section, we present both of these methods. We start by considering Eq. (33). We take a derivative ofthis equation with respect to c to obtain T N e ct N BB = (4π) N dz (z i z j ) 1 i (iπ) N ez (z N i<j i z j ) + c (z i<j i z j ) + c. (54) We introduce the auxiliary variables r and take the limit c + to recast this integral as T N BB = (4π) N dz (iπ) N ez dre (zi zj++ )r. (55) i N where the + term is explicitly written for convergence. We may then change variables i<j

14 CONTENTS 14 from z l ik l and compute the integrals for l i, j, yielding N(N 1) dk 1 dk T N BB = 4π dr (π) = N(N 1) dre r = N(N 1) e k 1 k e i(k 1 k )r + r π. (56) Alternatively, this computation can be done by introducing the Brownian bridge propagator (x x ) e 4Dτ(1 τ) P BB (x, τ x,, 1) =. (57) 4πDτ(1 τ) where the symbol P BB (x, τ x,, 1) denotes the propagator of the Brownian from x at time t = to x at time τ conditioned on being x at time t = 1. The mean value of the coincidence time of our process can simply be computed as (now taking D = 1) T N BB = 1 dτ i j δ(x i (τ) x j (τ)) = = N(N 1) 1 dτ τ(1 τ) 1 dx dτ i j 4π e x N(N 1) = P BB (x i, τ x,, 1) dx i π, (58) which coincides with the above result. Variance of T N,BB The second moment of the distribution can be obtained in a similar fashion either using the cumulant generating function in Eq. (33) or the Brownian bridge propagator in Eq. (57). Using the propagator we obtain that there are three different kinds of contributions 1 1 TN = dτ 1 dτ δ(x i (τ 1 ) x j (τ 1 ))δ(x m (τ ) x l (τ )) (59) = [ N! 1 (N 4)! N! + 8 (N 3)! N! + 4 (N )! 1 dτ i j l m P BB (x, τ, ) dx 1 dτ 1 dτ P BB (x, τ 1,, 1) τ 1 P (y x, τ τ 1, 1)P BB (y, τ,, 1)dxdy 1 dτ 1 1 τ 1 dτ ] P (x, τ 1, ) P (y x, τ τ 1, 1) dxdy, where P (x, τ, ) = (4πτ) 1/ e x /(4τ) is the free Brownian propagator. The first N! term in this expansion gives a disconnected contribution 4(N 4)! T. After a careful computation, we obtain the variance for the Brownian bridge [ Var (T N ) BB = [T N T N BB ] N! 8π BB = (N 3)! 9 3 π ] N! [ + 1 π ]. (N )! 4 (6) This is the last result that we obtain for arbitrary N Brownian bridges. In the limit of a large number of walkers, N + limit we find T N BB 1 ( π 8π N, Var (T N ) BB 9 3 π ) N 3 (61)

15 CONTENTS 15 This scaling with N is characteristic of the regime of typical fluctuations. We will now explore another regime, the large deviation regime, dominated by rare fluctuations..6. Coulomb gas method for N We now study the limit of a large number of walkers N. We will study the moment generating function of the coincidence time T N, in that limit for c >. As we show below, there is a natural scaling regime, c N, which allows to use a Coulomb gas approach. This regime corresponds to large positive c, hence to events where the coincidence time T N is much smaller than its average T N. Exponentiating the integrand of eq. (3), we get e ct N (t) BB = π N/ (6) where the action is given by S N (k) = i k i + i<j N dke S N (k) [ c ] log 1 + (k i k j ) For the two terms to be of the same order we must scale k i = Np i and c = N c. Introducing the empirical density ρ(p) = 1 N i δ(p p i) we rewrite the action as (63) S N (k) = N S c [ρ] + o(n ) (64) where we have defined the Coulomb gas energy functional S c [ρ] = dp p ρ(p) + 1 dpdp ρ(p)ρ(p ) log [ 1 + c ] (p p ). (65) We can evaluate the integral in (6) using a saddle point approximation and obtain [ ( )] e ct N c (t) BB exp N Ψ, (66) N where from its definition the rate function Ψ( c) is defined on c [, + [ and must be positive, increasing Ψ ( c) and concave Ψ ( c), with Ψ() =. It is determined by minimizing the Coulomb gas energy [ ( )] Ψ( c) = min ρ S c [ρ] µ ρ(p)dp 1 with respect to the density ρ(p), where we have introduce a Lagrange multiplier µ to enforce the normalization constraint dpρ(p) = 1. We denote ρ ( c, p) the minimizer of the action S c [ρ] under the normalisation constraint. Note that for c =, there is no normalised density allowing to minimise S[ρ]. Computing the functional derivative, we obtain an integral equation for ρ ( c, p), valid for p in the support of ρ ( c; p ) [ p + dp ρ ( c; p c ] ) log 1 + (p p ) = µ( c). (68) (67)

16 CONTENTS 16 Multiplying this equation with ρ ( c; p) and integrating with respect to p, we obtain [ dpdp ρ ( c; p)ρ ( c; p c ] ) log 1 + (p p ) = µ( c) dpρ ( c; p)p. (69) eplacing in Eq. (65), this yields a simpler expression, which can be interpreted as a virial theorem, for the rate function Ψ( c) for general value of c Ψ( c) = dp p ρ ( c; p) + µ( c). (7) Note also that taking a derivative of Eq. (68) with respect to p, we obtain an equation that does not depend on the Lagrange multiplier µ and reads, for any p in the support of ρ ( c; p ) ρ ( c; p )dp c p = p p c + (p p ). (71) Solving this integral equation for arbitrary values of c is quite non-trivial and left for future studies. Here we will only solve this equation perturbatively in the regime c. The density ρ ( c; p) can be expressed as a perturbative series of c ρ ( c; p) = ρ sc (p) + 1 c δρ 1(p) + O( ), (7) 1 c 4 where the leading order of the density ρ sc (p) does not depend on c. It is solution of the integral equation ρ sc (p )dp p p = p. (73) b ρ(p The general solution of the equation )dp a p p = g(p) (i.e. assuming a single interval support) can be obtained by the Tricomi formula [14] ρ(p) = 1 π p a b p [ b a b ρ(p )dp dp p a ] b p a π p p g(p ). (74) We expect in this case that a = b and that the density ρ sc (p) is normalised. Inserting these two conditions, we obtain ρ(p) = + a p π a p, a p a. (75) Imposing that this density vanishes at the edges in ±a, we obtain a =. solution of Eq. (73) is then the Wigner semi-circle law The ρ sc (p) = p π, p. (76) Computing the next order term, we obtain for p [, ] δρ 1 (p )dp p p = ρ sc (p )(p p )dp = p. (77)

17 CONTENTS 17 The solution at first order of perturbation must have a vanishing integral δρ 1(p)dp for the density to be normalised. in Eq. (74), this yields Using again the Tricomi formula δρ 1 (p) = 1 p π p, p. (78) We may then compute the rate function at order c. Evaluating the Lagrange multiplier from Eq. (68) at, e.g. p =, we obtain µ( c) = log c dpρ sc (p) log p + 1 c [ ] dpρ sc (p)p dpδρ 1 (p) log p + O( c 4 ) = log c log + 1 c (1 + 1) + O( c 4 ). where we used the large c expansion [ c ] log 1 + (p p ) = log c log p p + 1 c (p p ) + O( ) (8) 1 c 4 Gathering the different terms, we obtain from (7) the large c behavior Ψ( c) = log c log O( ) (81) c 1 c 4 which is consistent, to leading order with Ψ( c) being increasing and concave. This behaviour, obtained here in the regime N, c + with large c = c/ N, can be compared with the large c expansion at fixed N that we obtained in Eq. (48), 1 N log e ct N log G(N + ) BB c + N + N + log (79) N(N 1) N (log c + 1 log ) (8) [ ] c + 3 N log (83) where we used the asymptotic behaviour of the Barnes function log G(z + ) = (z /) log z 3z /4 + O(z). This estimate coincides precisely with the first three leading terms at large c in (81). The matching of these two regimes is illustrated in the Figure 4. If we assume that the above variational problem has a well behaved solution for any c, then the MGF has the large deviation form Eq. (66). That is then consistent with the PDF of T having the following large deviation tail in the regime T N 3/ P N,BB (T ) e N G( T ) T, T =. (84) N 3/ Indeed substituting this form in the definition of the Laplace transform, the integral is dominated at large N by the maximum of its integrand + dt P N,BB (T ) e ct + dt e N (G( T )+ c T ) e N min T [G( T )+ c T ] (85)

18 CONTENTS 18 Hence the rate function G( T ) is then determined via a Legendre transform Ψ( c) = min [G( T ) + c T ], G( T ) = max [Ψ( c) c T ] (86) T c Denoting the saddle points T ( c) and c ( T ) we have T ( c) c ( T ) 1 T T at small T. We obtain the small T expansion 1 c 1 c 3 at large c and G( T ) = log( 1 T ) T log + + o( T ) (87) Again, the first three terms match exactly the large N limit of (49). The large deviation regime studied above thus corresponds to events of probability log P = O(N ) in the double limit N, T + with T N 3/, i.e. to coincidence times T much smaller than the average (and typical) one T = T N BB. It equivalently corresponds to configurations of Brownian paths which strongly repel each others, i.e. in the Laplace variable this is the regime c N >. This is illustrated in Figure 4. It is then clear from the figure that the result (49) for the PDF, which corresponds to small T = O(1) and fixed N, should match for large N the large c limit. One cannot exclude however that there exists intermediate deviation regimes, e.g. for T T N BB N. One notes for instance that the upper tail at fixed N and large T behaves as P N (T ) e T /N 3 (see formula (51)). Inserting T T N BB we see that this factor in the upper tail can be of order e O(N) at most. Studying this regime at large N would require to extend the above Coulomb gas on the c < side, which is left for future studies. 3. Coincidence time for Brownian motions and KPZ equation with flat initial condition Let us now turn to the case of the Brownian motion, where the final points are not fixed (i.e. they are integrated upon). This is connected to the directed polymer with one free endpoint, equivalently with the Kardar-Parisi-Zhang equation with a flat initial condition. A solution for the latter was given in [15, 16]. We will use this solution here both for c > and c <. Since it was obtained by first calculating the moments for the half-flat initial condition we start by the latter, which corresponds to N Brownians all starting at and ending on a given half line Bethe ansatz solution of the Lieb-Liniger model for half-flat initial conditions Consider the partition sum of the directed polymer with half-flat initial condition Z w (x, t) = dye wy Z(x, y, t) (88) From Eq (5) in [17], or Eq. (88) in [16], we have the following: restoring the factors of c (by a change of units), with c = c, restricting to c and retaining only the single particle states (as they are the only states of the repulsive Lieb-Liniger model),

19 CONTENTS 19 P N (T) typical fluctuations e N G ( T N 3/ ) N 3/ e T N e N F( T N ) T N T 1 T N 3/ T=< T N N > T Figure 4. Schematic representation of the different regimes for the PDF of the coincidence time T, P N (T ), in the limit of large N. The average T is T N N, and the typical fluctuations live in a window T T N N 3/ around it. The large deviation regime studied in the text via the Coulomb gas is represented in blue and correspond to T N 3/ T N and has rate N. We have shown that its left part, which corresponds to T N 3/ (and c 1 for the corresponding MGF) matches the regime T 1 which is also indicated. The matching of its right part, i.e. N 3/ T N (and c 1 for the corresponding MGF), towards the typical regime remains to be studied. We note that the right tail of the regime of typical fluctuations is consistent with the existence of a right large deviation regime for T T N N of rapidity N, which remains to be studied. we obtain Z w (x, t) N = N j=1 dk j π e k j t ixkj ik j + w 1 i<j N (k i k j ) ik i + ik j + w + c (k i k j ) + c ik i + ik j + w From this we obtain the MGF for the coincidence time of N Brownian walkers which start at x and reach the half line ], ], as e ct N half BM = Z w(x, t) N Zw(x, t) N, dk e k t ixk Z w(x, t) = (9) π ik + w It is written here in presence of an extra weight e wx, which can then be considered in the limit w = Flat initial conditions as a limit of the half-flat The moments of the partition sum of the directed polymer with one free endpoint i.e. Z flat (t) = dyz(x, y, t), associated to the Kardar-Parisi-Zhang equation with flat initial condition, can be obtained from the ones for the half-flat initial condition in the double limit x and w +, as performed in [16]. In that limit only paired strings and strings with zero momenta remain: this provides the MGF for the coincidence time of unconstrained Brownians walkers all starting at. (89)

20 CONTENTS Moment generating function and probability distribution function for N =, 3 We start with the lowest moments, from eqs. (5) and (56) in [16], and consider c >, in which case we keep only the contribution of single particle states. The first moment is simply Z flat (t) = 1 for all times t. We consequently focus on the moments at time t = 1. N = From the second moment we obtain the result for two walkers e ct BM = Z flat (1) dk e k = 4c π 4k + c = e c c (1 erf( )) (91) The inverse Laplace transform then yields P,BM (T ) = L 1 e ct BM = dk π e k 4 cos(kt ) = π e T This result compares very well with the numerical simulations in Fig. 5. (9) numerical data P, BM (T) T numerical data P, BM (T) T Figure 5. Comparison between the coincidence time obtained numerically for N = Brownian motions on the time interval τ [, 1] and diffusion coefficient D = 1 and the analytical result in Eq. (9) (some details on the simulations are provided in Appendix C). Left: Linear scale, ight: Logarithmic scale. N = 3 The third moment gives the result for three walkers e ct3 BM = Z flat (1) 3 dk k e k = 1c π (k + c )(4k + c ) The inverse Laplace transform then yields P 3,BM (T ) = L 1 dk e ct3 BM = = π 4e k π e T (e 3 (cos(kt ) cos(kt )) 8 T 1) This result compares very well with the numerical simulations in Fig. 6. (93) (94) 3... Moment generating function for arbitrary N In the case of Brownian motions, the expression of the MGF of the coincidence time of N walkers depends on the parity of N:

21 CONTENTS numerical data P 3, BM (T) T numerical data P 3, BM (T) 3 4 T e T T Figure 6. Comparison between the coincidence time obtained numerically for N = 3 Brownian motions on the time interval τ [, 1] and diffusion coefficient D = 1 and the analytical result in Eq. (94). Left: Linear scale, ight: Logarithmic scale. For N even For N odd e ct N BM = (c) N N! ( N )! e ct N BM = 1 p<q N N p=1 N 1 (c) N! ( N 1 )! 1 p<q N 1 dk p π e k p 4k p + c (95) (k p k q ) (k p + k q ) (k p k q ) + c (k p + k q ) + c N 1 p=1 dk p π k pe k p (k p + c )(4k p + c ) (k p k q ) (k p + k q ) (k p k q ) + c (k p + k q ) + c These results are obtained from Eq. (18) in [16] by noting that for c > and N even only particle states with paired momenta (k 1, k 1,..., k N/, k N/ ) contribute, while for N odd there are N 1 paired momenta and one zero momentum Probability Distribution Function of the coincidence time In this section, we obtain from Eqs. (95) and (96) an expression for the PDF of the coincidence time for N Brownian motions. For N even For p [1, N/], define the conjugate variables (96) X p = c ik p, X p 1 = c + ik p (97) and observe that (k p k q ) (k p + k q ) (k p k q ) + c (k p + k q ) + c = (X q X p )(X q 1 X p 1 )(X q 1 X p )(X q X p 1 ) (X p + X q )(X p 1 + X q 1 )(X p 1 + X q )(X p + X q 1 ) (98)

22 CONTENTS Then one rewrites Eq. (95) as e ct N BM = (c) N N! ( N )! N p=1 dk p π e k p N p=1 1 X p X p 1 X p + X p 1 X p X p 1 1 p<q N Using the explicit values of X p and X p 1 we simplify the expression as e ct N BM = cn N! ( N )! N p=1 dk p π ie k p k p N 1 X p=1 p 1 p<q N X q X p X q + X p (99) X q X p X q + X p (1) Using formulae (9.1) and (9.4) of Bruijn [18], we write the last part of (1) as N 1 X p=1 p 1 p<q N X q X p = dr e c l r l X q + X p N + 1 k<l N sign(r k r l ) N p=1 e ikp(rp 1 rp) (11) From this identity, we obtain after anti-symmetrizing the last exponential w.r.t r p and r p 1 e ct N BM = cn N! ( N )! N + dr e c l r l = cn N! N ( N )! dr e c l r l N + 1 k<l N 1 k<l N sign(r k r l ) N p=1 sign(r k r l ) N p=1 dk p π e k p sin(k p(r p 1 r p )) ( ) rp 1 r p erf k p (1) where we computed the integral w.r.t {k p } s. The final trick consists in (i) changing the labels in the r variables in the error function from r p to r σ(p) where σ belongs to the symmetric group S N, (ii) using the fact that there are N!/( N/ (N/)!) ways of pairing N objects (N is even) and (iii) using the definition of the Pfaffian of an antisymmetric matrix Pf(A) = σ S N,σ(p 1)<σ(p) sign(σ) N/ p=1 A σ(p 1),σ(p) to finally obtain e ct N BM = c N N + dr e c l r l 1 k<l N sign(r k r l ) [ ( )] rk r l Pf erf 1 k,l N (13) Inverting the Laplace transform of this expression, we obtain the PDF P N,BM (T ) of the rescaled coincidence time T N = T N (t = 1) as P N,BM (T ) = N T N + dr δ ( = N+1 T dr Θ T N + ( T ) N r l l=1 ) N r l l=1 1 k<l N 1 k<l N sign(r k r l ) sign(r k r l ) [ ( )] rk r l Pf erf 1 k,l N [ ( )] rk r l Pf erf 1 k,l N (14)

23 CONTENTS 3 where Θ(x) is the Heaviside step function. Even though the product of sign functions can itself be written as a Pfaffian, we have not tried to simplify further Eq. (14). For N odd For p [1, (N 1)/], define the conjugate variables X p = c ik p, X p 1 = c + ik p (15) and define X N = c. Then, in a similar fashion as the even N case, one rewrites Eq. (96) as N 1 e ct N BM = cn N! dk p ie k p ( N 1 )! p=1 π k p N 1 X p=1 p 1 p<q N X q X p X q + X p (16) By the same argument as in the even N case, using Bruijn s results [18], the moment generating function reads e ct N BM = c N N! N 1 ( N 1 )! N + dr e c l r l 1 k<l N sign(r k r l ) N 1 p=1 ( ) rp 1 r p erf (17) We have to be careful in the odd case as there are N variables r but only N 1 of them are involved in the error functions whereas all N variables are involved in the product of sign functions. We now wish to employ the same Pfaffian trick as the even N case and hence have to consider all (N 1)!/( N 1 ( N 1 )!) ways of pairings N 1 objects. Hence, we obtain our final expression for the moment generating function in the odd N case e ct N BM = Nc N dr e c [ ( )] l r rk r l l sign(r k r l ) Pf erf 1 k,l N 1 N + 1 k<l N (18) Inverting the Laplace transform of this expression, we obtain the PDF P N,BM (T ) of the rescaled coincidence time T N = T N (t = 1) as P N,BM (T ) (19) ( ) N [ ( )] = N T N rk r l dr δ T r l sign(r k r l ) erf N + l=1 ( = N N+1 T dr Θ T N + ) N r l l=1 1 k<l N 1 k<l N sign(r k r l ) Pf 1 k,l N 1 [ ( )] rk r l Pf erf 1 k,l N 1 (11) where Θ(x) is the Heaviside step function. We have verified that for N =, 3, these formulae give back the results for the PDF and the MGF for the coincidence time of Brownian motions. For N =, 3 we only have to use a Pfaffian, Pf A ij = A 1, 1 i,j hence the calculation is elementary.

24 CONTENTS Asymptotic behaviour of P N,BM (T ) for arbitrary N Small T limit of the Probability Distribution Function We now discuss the small T behavior of the PDF P N,BM (T ) of T N = T N (t)/ t. It can be extracted from the c + limit of e ct N BM. When c is increased, the interaction in the Lieb- Liniger model becomes more repulsive and the corresponding Brownian trajectories with small coincidence time are those where none of the Brownian are bounded together. Taking the large c + limit of Eqs. (95) and (96) we see that the leading term is in all cases c N(N 1)/. Inverting the Laplace transform, we obtain the small T behavior of P N,BM (T ) as where I N depends on the parity of N. For N even P N,BM (T ) = I N T N(N 1) 1 + O(T N(N 1) +1 ), (111) T I N = For N odd = N N! ( N )!Γ ( N(N 1) N(N ) ) N p=1 dk p π e k p N 1 k=1 1 p<q N (k p k q) Γ(N + 1) ( ) Γ(k + 1) (11) (π) N 4 Γ N(N 1) I N = N 1 N! ( ( N 1 )!Γ N(N 1) ) N 1 p=1 Γ( N = + 1) ( ) N(N 3) (π) N+1 4 Γ N(N 1) dk p π k pe k N 1 k=1 p 1 p<q N 1 (k p k q) Γ(k + 1). (113) The k integrals are typical examples of Selberg integrals, and further details can be found in ef. [11] in Section 1.4 and after changing variables from {x i } s to {k i = x i / }. One can check that these results agree with the small T expansion of the formula for the PDF in the cases N =, 3 given above, i.e. I = I 3 = 3 4 π Large T limit of the Probability Distribution Function Similarly to Section.4., we determine the large T tail of the PDF P N,BB (T ) of T N = T N (t)/ t by investigating the c limit of the Lieb-Liniger model. It again determined by the same ground state as in Section.4.. From Eqs. (65-66) in ef. [13] it follows that the MGF in the c limit reads π and e ct N BM = N 1 e c N(N 1) 1 [1 + O(e 1 4 N(N 1) c )] (114)

25 CONTENTS 5 As in Section.4. we find the T asymptotics P N,BM (T ) = N 1 3 πn(n 1) e α N T + O(e β N T ) (115) 3 with the exponential factors again given by eqs. (53), i.e. α N = N(N 1) and 3 β N = N 3 3N +N. For completeness, we compute this asymptotics explicitly for N =, 3. N = This matches eq. (9). N = 3 This matches eq. (94). P,BM (T ) = P 3,BM (T ) = π e T π e T (116) 8 + O(e T ) (117) 3.5. Mean and variance of the coincidence time for N Brownian motions Mean value of T N,BM To compute the first moment, we may use the propagator for a single Brownian motion with diffusion coefficient D (x x ) P (x, τ x, ) = e 4Dτ, (118) 4πDτ together with the fact that we consider i.i.d. variables. The mean value of the coincidence time of the Brownian motion can be obtained (setting D = 1) as T N BM = = 1 dτ i j δ(x i (τ) x j (τ)) = N(N 1) 1 dτ τ As expected, we have T N BM < T N BB. dx 4π e x 1 dτ i j N(N 1) = P (x i, τ x, ) dx i π. (119) Variance of T N,BM We may still use the free propagator to compute the second moment of the distribution. Using again the Brownian propagator, we obtain T N = = 1 1 dτ 1 dτ δ(x i (τ 1 ) x j (τ 1 ))δ(x m (τ ) x l (τ )) (1) [ N! 1 (N 4)! N! + 8 (N 3)! N! + 4x (N )! 1 1 dτ i j l m dτ 1 1 τ 1 dτ dτ 1 1 P (x, τ, ) dx τ 1 dτ ] P (x, τ 1, ) P (y x, τ τ 1, ) P (y, τ, )dxdy P (x, τ 1, ) P (y x, τ τ 1, ) dxdy.

26 CONTENTS 6 After a careful computation, we finally obtain [ Var (T N ) BM = [T N T N BM ] N! BM = (N 3)! 3 ] [ N! 1 + π (N )! 1 ]. π (11) We note that the variance ratio Var (T N ) BM /Var (T N ) BB decreases monotonically from for N = to for infinite N. The ratio of relative variances T N BB Var (T N) BM /(Var (T N ) BB T N BM ) decreases from for N = to for infinite N. 4. Coincidence time for arbitrary fixed final points As mentionned in Section.1 there is an alternative formula (33) to the one from the Bethe ansatz, for the moments of the directed polymer problem, obtained from the study of Macdonald processes. It turns out that this formula can be extended to arbitrary final points y [9, 1] (while the starting points are still all at zero, i.e. x = ) and for arbitrary values of c. From this formula, specialized to the case c >, we immediately obtain that for y 1 y y N, the Laplace transform of the PDF of the coincidence time for N such Brownian walkers is e ct N (t),y = (4πt) N dz y i (iπ) N et(z+ t ) N i<j z i z j z i z j + c. (1) We note that the r.h.s. of (1) is invariant by a global shift of the endpoints y i y i + w. This global shift corresponds to adding a linear drift to each of the Brownians which has no effect on the coincidence time Exact distribution for N =, 3 Distribution for N =. We start from Eq. (1) where we have set N = and t = 1. This yields dz e ct 1 dz y 1,y = 4π i (iπ) e(z1+ ) +(z + y ) z 1 z z 1 z + c, y 1 y. (13) Taking the inverse Laplace transform and introducing the change of variables z j ik j one obtains P (T, y 1, y ) = 4π = T dk 1 dk (π) dk 1 dk π (ik 1 ik )e (k1 iy 1 ) (k iy ) T (ik 1 ik ) e (k1 iy 1 ) (k iy ) it (k 1 k ). (14) The two Gaussian integrals can now be computed and one obtains P (T, y 1, y ) = 1 (T + y y 1 )e T (T +y y1), y 1 y. (15) Note that setting y = y 1, we recover the result for the Brownian bridge (4). Alternatively, the result for the Brownian motion (9) is obtained by computing P,BM (T ) = dy 1 y 1 e y 1 +y 4 dy P (T, y 1, y ). (16) 4π

27 CONTENTS 7 Finally, symmetrising the problem by releasing the constraint y 1 y, by defining P,sym (T, y 1, y ) = 1 (T + y y 1 )e T (T + y y1 ), (17) we may obtain the joint PDF of the coincidence time, the algebraic distance d between the Brownian defined as d = y 1 y (with d ) and their center of mass position u = (y 1 + y )/. It reads +u ) + ( d +u ) ( P (T, d, u) = e (d 4 P,sym T, d 4π + u, d ) + u = e u 8π (T + d )e 1 8 (T + d ). (18) We may then recover the joint PDF of d and T + in Eq. (1) by integrating over the center of mass position u, P joint (T, d) = P (T, d, u)du = (T + d ) 8π e 1 8 (T + d ) (19) Distribution for N = 3. The distribution for N = 3 can be obtained from Eq. (1) but its expression is quite cumbersome and not very enlightening. Therefore, we only reproduce here its asymptotic behaviours for y 1 < y < y 3 P 3 (T, y 1, y, y 3 ) = 1 π 4 6 e (y1+y3 y) (y y 1 )(y 3 y )(y 3 y 1 ) T + O(T 3 ), T 16 [ 4 T 4 (y3 y1) T 8 ] (T + y 3 y 1 ) 4, T. (13) Note that taking two equal points y = y 3, the PDF behaves for small T as P 3 (T, y 1, y, y ) = (y 1 y ) T 3 /4 + O(T 4 ), while for all equal points y 1 = y = y 3, we recover the expression in Eq. (47). 4.. Small T asymptotics of the PDF for any N and for any fixed final points From (1) one can extract the small T asymptotics for any N and any fixed final points y. In the large c limit it becomes e ct N,y (4π) N c N(N 1) dz y ) (z i z j ). (131) = (4π) N c N(N 1) = (c) N(N 1) i (iπ) N e(z+ N i<j dz i (iπ) N ez ( y i + y j ). N i<j (y j y i ) (13) i<j The second line follows from the fact that the first line is antisymmetric in the y i and upon the shift z i z i y i /, that is is a polynomial of degree N(N 1)/ in the