arxiv: v2 [math.pr] 23 May 2014

Transcription

1 Cores of Dirichlet forms related to Random Matrix Theory Hirofumi Osada, Hideki Tanemura arxiv: v2 [math.pr] 23 May 2014 Abstract We prove the sets of polynomials on configuration spaces are cores of Dirichlet forms describing interacting Brownian motion in infinite dimensions. Typical examples of these stochastic dynamics are Dyson s Brownian motion and Airy interacting Brownian motion. Both particle systems have logarithmic interaction potentials, and naturally arise from random matrix theory. The results of the present paper will be used in a forth coming paper to prove the identity of the infinite-dimensional stochastic dynamics related to the random matrix theories constructed by apparently different methods: the method of space-time correlation functions and that of stochastic analysis. 1 Introduction. May 26, 2014 In random matrix theory, one of the main issues is to clarify the distribution of the eigenvalues and its asymptotic behavior as the size of the matrices goes to infinity. The prototypes of random matrices are Gaussian ensembles, divided into three classes and called Gaussian orthogonal/unitary/symplectic ensembles (GOE/GUE/GSE), according to their invariance under conjugation by orthogonal/unitary/symplectic groups. The eigenvalue distributions of Gaussian random matrices of N N size are then given as ˇm N β (dx N) = 1 Z h N(x N ) β e β 4 xn 2 dx N, where dx n = dx 1 dx 2 dx N, x N = (x i ) R N, and Here and after Z denotes the normalizing constant. TheGOE,GUE, and GSEcorrespondto inversetemperature β = 1,2 and 4, respectively [9, 2]. The celebrated Wigner s theorem asserts that the empirical measure (1.1) 1 N N δ xj/ N. of the eigenvalues under the distribution ˇm N β (dx N) converges to the semicircle law ς(x)dx as N [2, 9], where ς(x)dx is the probability on R such that ς(x)dx = 1 (1.2) 4 x2 1 [ 2,2] (x)dx. 2π For a countable subset {x n }, we call the σ-finite measure ξ = n δ x n a configuration if it becomes a Radon measure. The set of all configurations on R is a Polish space equipped with the vague topology, and is called the configuration space over R. We call ξ unlabeled particles, and x = (x 1,x 2,...) labeled particles. Note that for N-particle systems, there exists an obvious bijection between the distribution of the unlabeled N-particles and the symmetric distribution of the labeled N-particles, where N N. We note that this is not the case for infinite particle systems. For a given distribution µ of N-unlabeled particles, we denote by ˇµ the symmetric density of the associated N-labeled particles in the sequel. To examine the behavior of the distribution of the configuration ξ = N i=1 δ x j under ˇm N β (dx N), N there are two typical scalings, called the bulk and the soft edge. The former corresponds to the scaling suchthat y j = Nx j, andthedistributionof{y j } N under ˇm N β (dx N) is given by h N (x N ) = N x i x j. i<j Faculty of Mathematics, Kyushu University, Fukuoka, , JAPAN, osada@math.kyushu-u.ac.jp Department of Mathematics and Informatics, Faculty of Science, Chiba University, 1-33 Yayoi-cho, Inage-ku, Chiba , JAPAN, tanemura@math.s.chiba-u.ac.jp ˇµ N bulk,β(dy N ) = 1 Z h N(y N ) β e β 4N y N 2 dy N. The latter corresponds to the scaling such that y j = N 1/6 (x j 2 N), and the distribution of {y j } N under ˇm N β is given by ˇµ N soft,β (dy N) = 1 Z h N(y N ) β e β 4N 1/3 y N 2N2/3 1 N 2 dy N, 1

2 Cores of Dirichlet forms 2 where 1 N = (1,1,...,1) R N. Let β = 2. Then the limit of µ N bulk,2 is the determinantal random point field µ sin,2 with sine kernel (1.3) K sin,2 (x,y) = sin(x y) π(x y), and the limit of µ N soft,2 is the determinantal random point field µ Ai,2 with Airy kernel (1.4) K Ai,2 (x,y) = Ai(x)Ai (y) Ai (x)ai(y), x y whereai denotes the Airyfunction andai its derivative [18, 19, 9]. It is proved that these random point fields are quasi-gibbsian in [13, 14]. We consider the dynamical scaling limit corresponding to the static limit mentioned above. For this we introduce the associated stochastic dynamics describing the time evolution of N-particle systems. Let X N (t) = (Xj N(t))N be the solution of the SDE (1.5) dx N j (t) = db j (t)+ N k=1,k j dt X N j (t) XN k (t) or the SDE with Ornstein-Uhlenbeck s type drifts (1.6) dxj N (t) = db j(t) 1 2N XN j (t)dt N dt + Xj N(t) XN k (t). k=1,k j These are called Dyson s Brownian motion model with β = 2, or simply the Dyson model [3]. The solution of (1.6) is a natural reversible stochastic dynamics with respect to ˇµ N bulk,2, and that of (1.5) is also natural but has no invariant probability measures. Both have the same N-limit as we see below. Let Ξ N sin be an unlabeled process defined by Ξ N sin(t) = N δ X N j (t). Suppose that the distribution of Ξ N sin (0) is µn bulk,2. Then Ξ N sin converges in distribution to the process Ξ sin whose generating function { M Ψ t sin [exp [f] E }] f m (x)ξ sin (t m,dx), m=1 0 t 1 t 2 t M, f m C 0 (R), 1 m M, is represented by the Fredholm determinant [ ] δ st δ(x y)+k sin (s,x;t,y)χ t (y), Det (s,t) {t 1,...,t M} 2, (x,y) R 2 R with the extended sine kernel K sin [20, 7]: K sin (s,x;t,y) 1 1 due u2 (t s)/2 cos{u(y x)}, t s = π 0 1 due u2 (t s)/2 cos{u(y x)}, t < s. π 1 For the soft edge scaling, we suppose that the distribution of X N (0) is ˇµ N bulk,2, and introduce the process Y N defined by (1.7) Y N (t) = 1 (N 2/3 t) 2N 2/3 N 1/3 t+ t2 N 1/3XN 4 corresponding to (1.5), and (1.8) Y N (t) = 1 N 1/3XN (N 2/3 t) 2N 2/3 corresponding to (1.6). Then the unlabeled processes Ξ N Ai(t) = N δ Y N j (t) convergeindistributiontotheprocessξ Ai whosegenerating function is represented by the Fredholm determinant with the extended Airy kernel K Ai [17, 6, 7]: K Ai (s,x;t,y) due u(t s)/2 Ai(u+x)Ai(u+y), t s 0 = 0 due u(t s)/2 Ai(u+x)Ai(u+y), t < s. From the fact that K sin (s,x;s,y) = K sin (x,y) and K Ai (s,x;s,y) = K Ai (x,y), we see that the processes Ξ sin and Ξ Ai are reversible with respect to µ sin,2 and µ Ai,2, respectively. The scaling limits above are based on the convergence of the associated space-time correlation functions by the determinatal structures [7]. There exists another method of constructing infinite volume stochastic dynamics based on stochastic analysis. In [13, 14], unlabeled diffusion processes Ξ sin and Ξ Ai with reversible measures µ sin,2 and µ Ai,2 are constructed through the Dirichlet form technique. Let X(t) = (X j (t)) j N be a labeled process associated with Ξ sin (t) = j N δ X j(t). In [12] it is proved that the process X = (X j ) j N solves the infinitedimensional stochastic differential equation (ISDE) (sin) dx j (t) = db j (t)+ k=1,k j dt X j (t) X k (t).

3 Cores of Dirichlet forms 3 In [15] we prove that a labeled process Y = (Y j ) j N associated with Ξ Ai (t) = j N δ Y j(t) solves the ISDE (Ai) dy j (t) = db j (t) + lim { 1 r Y j (t) Y k (t) k j Y k (t) <r r r ρ(x)dx x }dt, where ρ(x) = 1 (,0)(x) x. π These two approaches are fundamentally different. Hence it is significant to prove that the resulting stochastic dynamics are the same. From the former construction we can obtain quantitative information of the limit stochastic dynamics through the calculation of space-time correlation functions; while from the latter we deduce many qualitative properties of the sample paths of the labeled diffusions through the ISDE representation of the processes. Recently, we have proved the coincidence of these pairs of stochastic dynamics Ξ sin and Ξ sin, and also Ξ Ai and Ξ Ai, throughthefollowingthreesteps: Below denotes sin or Ai for the sake of brevity. (i) Ξ has the strong Markov property. (ii) TheDirichletformsassociatedwithΞ and Ξ are both extensions of the closable form (E µ,p). Here E µ aregivenby(2.11)andp isthesetofpolynomials functions on M defined in (2.13) later. (iii) The labeled process associated with Ξ and Ξ are solutions of the ISDE ( ), and the ISDE ( ) has strong uniqueness. The claim (i) is proved in [21], and (ii) is in this article. The claim (iii) is proved in [15, 16] partly through the result in this paper. Putting these together, we will complete the proof of Ξ = Ξ for {sin,ai} in a forthcoming paper. 2 Preliminaries Let S be a closed subset of R d such that the interior S int is a connected open set and that its closure S int equals S. Let M = M(S) be the configuration space over S of unlabeled particles, the set of non-negative integer valued Radon measures on S. The space M is a Polish space endowed with the vague topology. An element ξ of M can be represented as ξ = j Λ δ x j for some countable set Λ, and the restriction of ξ on a subset A of S is denoted by ξ A = ξ( A). A function f on M is called local if f(ξ) = f(ξ K ) for some compact set K. We write ξ K = k δ x j. For a local function f with f(ξ) = f(ξ K ) we introduce the functions ˇf k on S k, k N 0 {0} N defined by ˇf 0 = f( ), where is the null configuration, and by, for k N, k ˇf k (x k ) = f δ xj for x k K k. We extend the domain of ˇfk (x k ) to S k \K k by the consistency coming from f(ξ) = f(ξ K ). Hence ˇf k, k N 0, satisfy the consistency relation (2.9) ˇfk+1 (x k,y) = ˇf k (x k ), x k S k, y / K. The infinite sequence given by (2.10) (ˇf 0, ˇf 1 (x 1 ), ˇf 2 (x 1,x 2 ),...) = (ˇf k (x k )) k 0 is a representation of the local function f. A local function f is called smooth if the ˇf k are smooth for k N 0. We denote by D the set of all local smooth functions on M. Set for x k = (x i ) k i=1 Sk, k N 0, f,g D D(f,g)(x k ) = 1 2 k d i=1 ˇf k (x k ) x ij ǧ k (x k ) x ij, where x i = (x i1,x i2,...,x id ). For given f,g D, the right hand side is a permutation invariant function, and the square field D(f,g) can be regarded as a local function with variable ξ = i N δ x i M. For a probability µ on M, we denote by L 2 (M,µ) the space of square integrable functions on M with the inner product, µ and the norm L 2 (M,µ). We consider the bilinear form (E µ,d µ ) on L2 (M,µ) defined by (2.11) E µ (f,g) = D(f, g)dµ, (2.12) M D µ = {f D : f 2 1 < }, where f 2 1 Eµ (f,f) + f 2 L 2 (M,µ). A function F on M is called a polynomial function if F is given as (2.13) F(ξ) = Q( φ 1,ξ, φ 2,ξ,..., φ l,ξ ) with φ k C 0 (R d ) and a polynomial function Q on R l, where φ,ξ = R d φ(x)dx and C 0 (Rd ) is the set of smooth functions with compact support. We denote by P the set of all polynomial functions on M, and by P 0 if we replace the set of polynomials Q on R l by C b (Rl ), the set of bounded smooth functions with bounded derivaives of any order. It is obvious that each element of P and P 0 is a local smooth function. The closability and the quasi-regularity of the bilinear form (E µ,d µ ) have been provedin [10, 13, 14], while those of (E µ,p 0 ) in [22, 1].

4 Cores of Dirichlet forms 4 Let P µ, P µ 0, and Dµ be the closures with respect to f 1 of P, P 0, and D, µ respectively. We see that P D µ and P 0 D µ under the mild assumption (A.0) below, and hence we obtain that (2.14) P µ D µ, P µ 0 Dµ. Then we deduce from Theorem 3.1 below that (2.15) P µ = P µ 0 = Dµ. The construction of unlabeled diffusions of interacting Brownian motion in infinite dimensions through the Dirichlet form approach was initiated by [10]. Later [22, 1] also used this approach but with different cores P 0 under a more restrictive assumptions of interaction potentials than [10]. The identity (2.15) above proves that these diffusions are the same. We refer to [8] and [4] for the notion of quasiregularity and Dirichlet forms. 3 Main results. We call a function ρ n the n-correlation function of µ with respect to the Lebesgue measure if ρ n : S n R is a permutation invariant function such that A k 1 1 Akm m = m M i=1 ρ n (x 1,...,x n )dx 1 dx n ξ(a i )! (ξ(a i ) k i )! dµ(ξ) for any sequence of disjoint bounded subsets A 1,...,A n S and a sequence of natural numbers {k i } with k 1 + +k m = n. We assume the following conditions on the probability measure µ on M: (A.0) The measure µ has an n-correlation function ρ n for each n N with ρ n L p (S n r,dx n) for all k,r N for some 1 < p. (A.1) (E µ,d µ ) is closable on L2 (M,µ) and its closure (E µ,d µ ) is a quasi-regular Dirichlet form. From (A.0) we easily deduce (2.14). Hence from (A.1) we see that (E µ,p) are (E µ,p 0 ) are closable on L 2 (M,µ)aswellas(E µ,d µ ). Let(E µ,p µ ), (E µ,p µ 0 ), and (E µ,d µ ) be their closures as before. The main result of this paper is the following. Theorem 3.1 Suppose that µ satisfies (A.0) and (A.1). Then (E µ,p µ ) = (E µ,p µ 0 ) = (Eµ,D µ ). In [13, 14] the sufficient conditions (A.0) and (A.1) weregiven: ifµisa(φ,ψ)-quasi-gibbsmeasure,with Borel measurable functions Φ : S R { } and Ψ : S S R { } satisfying (3.16) (3.17) c 1 Φ 0 (x) Φ(x) cφ 0 (x) c 1 Ψ 0 (x y) Ψ(x,y) cψ 0 (x y) for a positive constant c, and some upper semicontinuous functions Φ 0 and Ψ 0 being locally bounded from below and with compact core {x : Ψ 0 (x) = }, then (A.0) and (A.1) are satisfied. We next apply Theorem 3.1 to the stochastic dynamics arising from the random matrix theory. Since the processes Ξ sin and Ξ Ai have the strong Markov property [21], they are associated with the quasi-regular Dirichlet forms (E sin,d sin ) and (E Ai,D Ai ), respectively. Thus we have the desired result as a corollary of Theorem 3.1. Corollary 3.2 Let {sin, Ai}. The Dirichlet forms associated with Ξ and Ξ are both extensions of the Dirichlet form (E µ,p µ ). Furthermore, P is a core of the Dirichlet form (E µ,d µ ), and D µ D. Proof. It is proved in Proposition 7.2 of [7] that E µ (f,g) = E (f,g), f,g P. From this we deduce the first claim. The second is a direct consequence of Theorem Proof of Theorem 3.1. For simplicity we only prove the case S = R d. We set A = {a = {a r } r N : a r N,a r a r+1,r N}. For a = {a r } A, let M[a] = {ξ M : ξ(s r ) a r, for all r}, where S r = {x S : x < r}, r N. Then M[a] is compact in M endowed with the vague topology. We introduce a cut off function X of M[a] as follows X[a](ξ) = h d a (ξ), d a (ξ) = (r x j (ξ) ) 2 r=1 j J r,ξ 1/2 Here {x j (ξ)} is a sequence in S such that ξ = δxj(ξ), x j (ξ) x j+1 (ξ) for all j, and J r,ξ = {j : j > a r, x j (ξ) S r }..

5 Cores of Dirichlet forms 5 Furthermore, h : R [0,1] is a function defined by 0, t (,0), (4.18) h(t) = 1 t, t [0,1], 1, t (1, ). Note that d a (ξ) = 0 and X a (ξ) = 1 if J r,ξ = for all r N. The following is proved in [10, Lemma 2.5]. Lemma 4.1 For any f D µ and ε > 0 we can take a A such that X[a]f D and that (1 X[a])f 1 < ε. Let ψ be a smooth function on R with support in [ 1,1] such that R ψ(x)dx = 1. Then we put ψ N (x) = Nψ(Nx). For g C0 (RL ) with support in [ r,r] L we associate the following function (4.20) N N ( 2j1 r g N (x L ) = g N r,..., 2j ) Lr N r = j 1=1 j L=1 L φ r,n,jl ψ N (x l ) l=1 with x L = (x l ) L l=1 RL and ( ) N φ r,n,j (x) = r N (r +x) j (r x) N j, j where f g stands for the convolution of f and g. Then by a simple observation we have the following. Lemma 4.2 Let g C 0 (RL ) with support in [ r,r] L. Then lim N R L { g g N 2 +D(g g N,g g N )}dx L = 0. Let L = dk and x L = x dk = (x 1,...,x k ), where x i S. If a function g(x 1,...,x k ) is symmetric in (x 1,...,x k ), then we can and do regard g as a functiondefinedontheconfigurationspacemovers with support in {ξ(s) = k}. If g(x 1,...,x k ) is symmetric in (x 1,...,x k ), where x j S, then so is g N. Hence, we deduce from Lemma 4.2 the following. Lemma 4.3 Let L = dk. Let g C 0 (RL ) be symmetric in (x 1,...,x k ) with support in [ r,r] L. Let g and g N be regarded as functions on M with support in {ξ(s) = k} as above. Then (4.19) lim g g N 1 = 0. N Proof of Theorem 3.1. We prove only (E µ,p µ ) = (E µ,d µ ) because the proof of the rest is similar. Let S n off = {(x 1,x 2,...,x n ) S n ;x i x j for i j}. Letf D suchthat f(ξ) = f(ξ K ) with compactset K = [ r,r] d. For such an f we introduce continuous functions f n, n N 0 such that f 0 = ˇf 0 and that, for n N and (x 1,x 2,...,x n ) S n off, f n (x 1,x 2,...,x n ) n = ( 1) n k k=0 {i 1,...,i k } {1,2,...,n} ˇf k (x i1,...,x ik ). The values of f n on {S n off }c are defined by continuity. Then f n is a smooth symmetric function on S n vanishing out of K for n 1. Note that f n is the Möbius transformation of ˇf k, k = 0,1,...,n. Then we easily deduce that, for (x 1,x 2,...,x n ) S n off, ˇf k (x 1,x 2,...,x k ) k n=0{i 1,...,i n} {1,2,...,k} f n (x i1,x i1,...,x in ). This implies that ˇf k ˇf 0 can be represented by a linear combination of symmetric smooth functions vanishing out of K k. In (4.21) and (4.22) below, (x i1,x i2,...,x in ) and (y 1,y 2,...,y n ) are taken to be in Soff n. The equalities can be exteded to {Soff n }c by continuity of the functions. Let k = ξ K (S) and write ξ K = k i=1 δ x i. Then from (4.20) we can rewrite f as (4.21) f(ξ) = ˇf k (x 1,x 2,...,x k ) = ξ K(S) n=0 η ξ K fn (y 1,y 2,...,y n ). Here n = η(s), η = n i=1 δ y i, and η ξ K means that η(a) ξ K (A) for all A. We note that the righthand side can be regarded as a symmetric function of (x 1,...,x n ) by construction. For m N we put (4.22) f [m] (ξ) = m n=0 η ξ K fn (y 1,y 2,...,y n ). Let ε > 0 be arbitrary. Then from Lemma 4.1 and (4.22) we can take a A and a r m N such that (4.23) f f [m] 1 (1 X[a])(f f [m] ) 1 < ε. From Lemma 4.3, we approximate the symmetric function f n in (4.22) by a polynomial F n. Hence, for any ε > 0, we can take polynomials F n such that (4.24) f [m] m F n 1 < ε. n=0

6 Cores of Dirichlet forms 6 Results (4.23) and (4.24) complete the proof. Acknowledgement H.O. is supported in part by a Grant-in-Aid for Scientific Research (KIBAN-A, No ) of the Japan Society for the Promotion of Science. H.T. is supported in part by a Grant-in-Aid for Scientific Research (KIBAN-C, No ) of the Japan Society for the Promotion of Science. References [1] Albeverio, S., Kondratiev, Yu. G., Röckner, M., Analysis and geometry on configuration space: the Gibbsian case, J. Funct. Anal. 157, (1998). [2] Anderson, G.W., Guionnet, A., Zeitouni, O., An Introduction to Random Matrices, Cambridge university press, [3] Dyson, F. J., A Brownian-motion model for the eigenvalues of a random matrix, J. Math. Phys. 3, (1962). [4] Fukushima, M., Oshima, Y., Takeda, M., Dirichlet forms and symmetric Markov processes, 2nd ed., Walter de Gruyter, Berlin, [5] Honda, R., Osada, H., Infinite-dimensional stochastic differential equations related to the Bessel random point fields, preprint, arxiv: [6] Johansson, K., Discrete polynuclear growth and determinantal processes, Commun. Math. Phys. 242, (2003). [7] Katori, M., Tanemura, H., Noncolliding Brownian motion and determinantal processes, J. Stat. Phys. 129, (2007). [8] Ma, Z.-M, Röckner, M., Introduction to the theory of (non-symmetric) Dirichlet forms, springer-verlag, [9] Mehta, M. L., Random Matrices. 3rd edition, Amsterdam: Elsevier, 2004 [10] Osada, H., Dirichlet form approach to infinitedimensional Wiener processes with singular interactions, Commun. Math. Phys. 176, (1996). [11] Osada, H., Non-collision and collision properties of Dyson s model in infinite dimensions and other stochastic dynamics whose equilibrium states are determinantal random point fields, in Stochastic Analysis on Large Scale Interacting Systems, eds. T. Funaki and H. Osada, Advanced Studies in Pure Mathematics 39, 2004, [12] Osada, H., Infinite-dimensional stochastic differential equations related to random matrices, Probability Theory and Related Fields, 153, (2012). [13] Osada, H., Interacting Brownian motions in infinite dimensions with logarithmic interaction potentials, Ann. of Probab. 41, 1-49 (2013). [14] Osada, H., Interacting Brownian motions in infinite dimensions with logarithmic interaction potentials II: Airy random point field, Stochastic Processes and their applications 123, (2013) [15] Osada, H., Tanemura, H., Infinite-dimensional stochastic differential equations arising from Airy random point fields, preprint. [16] Osada, H., Tanemura, H., Strong solutions of infinite-dimensional stochastic differential equations and tail σ-fields, preprint. [17] Prähofer, M., Spohn, H., Scale invariance of the PNG droplet and the Airy process, J. Stat. Phys. 108, (2002). [18] Shirai, T., Takahashi, Y., Random point fields associated with certain Fredholm determinants I: fermion, Poisson and boson point process, J. Funct. Anal. 205, (2003). [19] Soshnikov, A., Determinantal random point fields, Russian Math. Surveys 55, (2000). [20] Spohn, H., Interacting Brownian particles: a study of Dyson s model, In: Hydrodynamic Behavior and Interacting Particle Systems, G. Papanicolaou (ed), IMA Volumes in Mathematics and its Applications, 9, pp , Springer- Verlag, Berlin (1987). [21] Tanemura, H., Strong Markov property of determinatal processes with extended kernels, preprint. [22] Yoshida, M., Construction of infinitedimensional diffusion processes through Dirichlet forms, Probability Theory and Related Fields, 106, (1996).