SEA s workshop- MIT - July 10-14
|
|
- Posy Weaver
- 5 years ago
- Views:
Transcription
1 Matrix-valued Stochastic Processes- Eigenvalues Processes and Free Probability SEA s workshop- MIT - July July 13, 2006
2 Outline Matrix-valued stochastic processes. 1- definition and examples. 2- Multivariate statistics and multivariate functions. 3- Eigenvalues process : non-colliding particles. 4- The Hermitian complex case : determinantal processes. Root systems 1- β-processes (β > 0). 2- Radial Dunkl processes. Free probability : free processes.
3 Matrix-valued processes : definitions Let (Ω, F, (F t ) t 0, P) be a filtered probability space. A matrix-valued process is defined by X : R Ω M m (C) (t, ω) X t (ω) where M m (C) is the set of square complex matrices. Let S m, H m denote the spaces of symmetric and Hermitian matrices.
4 Examples Dyson model : B 1 = (Bij 1) i,j, B 2 = (Bij 2) i,j : two independent m m real Brownian matrices. { B 1 ii (t) i = j X ij (t) = i < j B 1 ij (t)+ 1B 2 ij (t) 2
5 Examples Dyson model : B 1 = (Bij 1) i,j, B 2 = (Bij 2) i,j : two independent m m real Brownian matrices. { B 1 ii (t) i = j X ij (t) = i < j B 1 ij (t)+ 1B 2 ij (t) 2 Wishart process (Bru, Donati-Doumerc-Matsumoto-Yor,): X (t) = B T (t)b(t), B : n m real Brownian matrix. Complex version: Laguerre process. n is the dimension and m is the size.
6 Examples Dyson model : B 1 = (Bij 1) i,j, B 2 = (Bij 2) i,j : two independent m m real Brownian matrices. { B 1 ii (t) i = j X ij (t) = i < j B 1 ij (t)+ 1B 2 ij (t) 2 Wishart process (Bru, Donati-Doumerc-Matsumoto-Yor,): X (t) = B T (t)b(t), B : n m real Brownian matrix. Complex version: Laguerre process. n is the dimension and m is the size. Other models: matrix Jacobi process (Y. Doumerc), Hermitian model of Katori and Tanemura (Brownian bridges).
7 Some SDE When it makes sense, one has for : Wishart (W (n, m, X 0 )), Laguerre (L(n, m, X 0 )) processes : dx t = dnt Xt + X t dn t + βni m dt, (β = 1, 2) Real and complex matrix Jacobi processes J(p, q, m, X 0 ): dx t = X t dnt Im X t + I m X t dn t Xt + β(pi m (p + q)x t )dt, (β = 1, 2) (N t ) t 0 is a square real Brownian matrix of size m.
8 Multivariate statistics t = 1, X 0 = 0, Dyson model, symmetric BM GUE and GOE. t = 1, X 0 = 0, Wishart, Laguerre LOE, LUE. Non-central Wishart and complex Wishart distributions with parameters M = X 0, Σ = βti m (James, Muirhead, Chikuze). Stationnary Jacobi matrix MANOVA.
9 Key tools : multivariate functions Hypergeometric function of matrix argument : pf β q ((a i ) 1 i q, (b j ) 1 j q, X ) = (a 1 ) τ (a p ) τ Jτ β (X ) (b 1 ) τ (b q ) τ k! k=0 τ β = 1 : zonal polynomial (Muirhead). β = 2 : Schur function (Macdonald) : J 2 τ (x 1,, x m ) = det(x k j +m j i ) det(x m j i ) Hypergeometric function of two matrix arguments : pf β q ((a i ) 1 i q, (b j ) 1 j q, X, Y ) = (a 1 ) τ (a p ) τ Jτ β (X )Jτ β (Y ) (b 1 ) τ (b q ) τ Jτ β (I m )k! k=0 τ
10 Some expressions Wishart and Laguerre semi-groups : 0 F β 1, β = 1, 2. Generalized Hartman-Watson Law : 0 F β 1, β = 1, 2. Tail distribution of T 0 := inf{t, det(x t ) = 0} : 1 F β 1, β = 1, 2. Real symmetric case : quite complicated. More precise results in the complex Hermitian case : determinantal representations of multivariate functions (Gross and Richards, Demni, Lassalle for orthogonal polynomials).
11 determinantal formulae pf 2 q ((a) 1 i p, (b j ) 1 j q, X ) = (1) det(x m j i pf q ((a j + 1) 1 i p, (b l j + 1) 1 l q, x i ) det(x m j i ) (Gross et Richards) pfq 2 ((m + µ i ) 1 i p, (m + φ j ) 1 j q ; X, Y ) = Γ m (m) p π m(m 1) (p q 1) (Γ(µ i + 1)) m q Γ m (m + φ j ) 2 Γ m (m + µ i ) (Γ(φ j + 1)) m (2) i=1 j=1 det ( p F q ((µ i + 1) 1 i p, (1 + φ j ) 1 j q ; x l y f ) l,f h(x)h(y) for (µ i ), (φ j ) > 1 (Gross and Richards, Demni).
12 The complex case : determinantal processes A determinantal process : R m -valued process with semi group writable as determinant : q t (x, y) = det(k t (x i, y j )) i,j, x, y R m In the complex Hermitian case, the eigenvalue process is determinantal : Weyl integration formula + determinantal representation of the two matrix arguments functions. Probability technics (Doob h-transform). q t (x, y) = V (y) ( 1 V (x) det exp (y j x i ) 2 ) 2πt 2t Laguerre process : K t = the squared Bessel process semi-group. i,j
13 The eigenvalues process : SDE Let X be a matrix-valued process and (λ 1, λ 2,..., λ m ) its eigenvalues process with starting point (λ 1 (0) > > λ m (0)). Then dλ i (t) = d X i (t) + SD(t), 1 i m X i L = Xii, independent. SD : singular drift showing interaction between particles.
14 Examples Symmetric and Hermitian Brownian matrix : dλ i (t) = db i (t) + β 2 j i dt λ i (t) λ j (t)
15 Examples Symmetric and Hermitian Brownian matrix : dλ i (t) = db i (t) + β 2 j i dt λ i (t) λ j (t) Wishart and Laguerre processes δ > m 1, m. dλ i (t) = 2 λ i (t)db i (t) + βδdt + β j i λ i (t) + λ j (t) λ i (t) λ j (t) dt
16 Examples Symmetric and Hermitian Brownian matrix : dλ i (t) = db i (t) + β 2 j i dt λ i (t) λ j (t) Wishart and Laguerre processes δ > m 1, m. Jacobi dλ i (t) = 2 λ i (t)db i (t) + βδdt + β j i λ i (t) + λ j (t) λ i (t) λ j (t) dt dλ i (t) = 2 λ i (t)(1 λ i (t)db i (t) + (p (p + q)λ i (t))dt + β j i λ i (t)(1 λ j (t)) + λ j (t)(1 λ i (t)) dt λ i (t) λ j (t) for β = 1, 2. What happens for arbitrary β > 0?
17 Root systems Let α R m \ {0} and let σ α denotes the reflection with respect to the hyperplane H α orthogonal to α : σ α (x) = x 2 < x, α > < α, α > α. A root system R is a non-empty subset of non-null vectors of R m satisfying : R Rα = {±α}, A simple system is a basis of R such that each α R is either a positive or negative linear combination of vectors of. The first kind of roots constitute the positive subsystem, denoted by R +, and are called by the way positive roots.
18 Root systems Let α R m \ {0} and let σ α denotes the reflection with respect to the hyperplane H α orthogonal to α : σ α (x) = x 2 < x, α > < α, α > α. A root system R is a non-empty subset of non-null vectors of R m satisfying : R Rα = {±α}, σ α (R) = R, α R A simple system is a basis of R such that each α R is either a positive or negative linear combination of vectors of. The first kind of roots constitute the positive subsystem, denoted by R +, and are called by the way positive roots.
19 Weyl Group Weyl group W : Multiplicity function k : constant on each orbit. W := span{σ α, α R} O(R m ) Positive Weyl chamber : k : {orbits of W } R α k(α) C := {x R m, < α, x >> 0 α R + }
20 Radial Dunkl process Let C be the closure of C. The radial Dunkl process X W is a paths- continuous C-valued Markov process with extended generator given by : L u(x) = 1 2 u(x) + k(α) < u(x), α > < x, α >, α R + where u C 2 0 (C) such that < u(x), α >= 0 for x H α, α R +.
21 A m 1 -type R = {±(e i ± e j ), 1 i < j m}. R + = {e i e j, 1 i < j m}. = {e i e i+1, 1 i m 1}. C = {x R m, x 1 > > x m }. k = k 0 > 0, (e i ) 1 i m is the standard basis of R m.
22 A m 1 -type R = {±(e i ± e j ), 1 i < j m}. R + = {e i e j, 1 i < j m}. = {e i e i+1, 1 i m 1}. C = {x R m, x 1 > > x m }. k = k 0 > 0, (e i ) 1 i m is the standard basis of R m. Set k 0 = β/2, β > 0. dλ i (t) = db i (t) + β 2 j i dt λ i (t) λ j (t), 1 i m. with λ 1 (0) >... λ m (0) (Cépa and Lépingle).
23 A m 1 -type R = {±(e i ± e j ), 1 i < j m}. R + = {e i e j, 1 i < j m}. = {e i e i+1, 1 i m 1}. C = {x R m, x 1 > > x m }. k = k 0 > 0, (e i ) 1 i m is the standard basis of R m. Set k 0 = β/2, β > 0. dλ i (t) = db i (t) + β 2 j i dt λ i (t) λ j (t), 1 i m. with λ 1 (0) >... λ m (0) (Cépa and Lépingle). X 0 = 0, t = 1 β-hermite ensemble (Edelman-Dumitriu).
24 The B m -type R = {±e i, 1 i m, ±(e i ± e j ), 1 i < j m}. R + = {e i, 1 i m, e i ± e j, 1 i < j m}. = {e i e i+1, 1 i m 1, e m }. C = {x R m, x 1 > x 2... x m > 0}. Two conjugacy classes k = (k 0, k 1 ).
25 β-laguerre processes (Demni) Let β, δ > 0. A β-laguerre process (λ(t)) t 0 starting at (λ 1 (0) > > λ m (0)) is a solution when it exists of dλ i (t) = 2 λ i (t)dν i (t)+β δ + λ i (t) + λ j (t) dt, λ i (t) λ j (t) i j 1 i m for t < τ, the first collision time, where (ν i ) are independent BM. Let R 0 := inf{t, λ m (t) = 0 for some i} then T 0 = τ R 0, moreover :
26 A B m - Radial Dunkl process (r(t)) t T0 = ( λ(t)) t T0 satisfies : dr i (t) = dν i (t)+ β 2 + [ j i 1 r i (t) r j (t) + 1 r i (t) + r j (t) β(δ m + 1) 1 dt, 2r i (t) 1 i m (r t ) t T0 is a B m -radial Dunkl process with multiplicity function given by 2k 0 = β(δ m + 1) 1 > 0 and 2k 1 = β > 0. ] dt
27 A B m - Radial Dunkl process (r(t)) t T0 = ( λ(t)) t T0 satisfies : dr i (t) = dν i (t)+ β 2 + [ j i 1 r i (t) r j (t) + 1 r i (t) + r j (t) β(δ m + 1) 1 dt, 2r i (t) 1 i m (r t ) t T0 is a B m -radial Dunkl process with multiplicity function given by 2k 0 = β(δ m + 1) 1 > 0 and 2k 1 = β > 0. ] dt X 0 = 0, t = 1 β-laguerre ensemble (Edelman-Dumitriu).
28 A unique strong solution for all t 0 Theorem 1 (Demni): Let B be a m-dimensional BM. Then, the radial Dunkl process (X W t ) t 0 is the unique strong solution of the SDE dy t = db t Φ(Y t )dt, Y 0 C where Φ(x) = α R + k(α) ln(< α, x >) for k(α) > 0 α R +.
29 A β-matrix model (Demni) Q : Is there a matrix-valued process corresponding to : dλ i (t) = db i (t) + β 2 j i dt λ i (t) λ j (t) 1 i m
30 A β-matrix model (Demni) Q : Is there a matrix-valued process corresponding to : dλ i (t) = db i (t) + β 2 j i A : β-hermitian model, 0 < β 2 : Bii 1(t) X ij (t) = β 2 dt λ i (t) λ j (t) ( ) Bij 1(t)+B2 ij (t) 2 1 i m i = j i < j where B 1 ii, B2 jj t = (1 β/2)t and B 1, B 2 are two ind. m m Brownian matrices.
31 Free probability : free processes Non-commutative probability space : unital algebra A + state A C, Φ(1) = 1. Examples A m = p>0 L p (Ω, F, (F t ) t 0, P) M m (C) the set of m m random matrices with finite moments, and the normalized trace expectation : Φ m := 1 m E(tr) := E(tr m) B(H) : the set of bounded linear operators on a Hilbert space H with the pure state Φ(a) =< ax, x >, a B(H), where x H is a unit element.
32 Other properties - algebra C or W - non commutative probability space. involutive Banach algebra : norm s.t a = a, a A + completion. C -algebra involutive Banach algebra + aa = a 2 for all a A. Φ can be : 1. tracial : Φ(ab) = Φ(ba). In the matrix example, involution has to be the usual adjonction and conditions are obviously fulfilled. As in classical probability, we endow our space with a family (A t ) t 0 of increasing C -subalgebras called filtration conditional expectation Φ(a t /A s ), s t.
33 Other properties - algebra C or W - non commutative probability space. involutive Banach algebra : norm s.t a = a, a A + completion. C -algebra involutive Banach algebra + aa = a 2 for all a A. Φ can be : 1. tracial : Φ(ab) = Φ(ba). 2. faithful : Φ(aa ) = 0 a = 0. In the matrix example, involution has to be the usual adjonction and conditions are obviously fulfilled. As in classical probability, we endow our space with a family (A t ) t 0 of increasing C -subalgebras called filtration conditional expectation Φ(a t /A s ), s t.
34 Other properties - algebra C or W - non commutative probability space. involutive Banach algebra : norm s.t a = a, a A + completion. C -algebra involutive Banach algebra + aa = a 2 for all a A. Φ can be : 1. tracial : Φ(ab) = Φ(ba). 2. faithful : Φ(aa ) = 0 a = normal In the matrix example, involution has to be the usual adjonction and conditions are obviously fulfilled. As in classical probability, we endow our space with a family (A t ) t 0 of increasing C -subalgebras called filtration conditional expectation Φ(a t /A s ), s t.
35 Examples A free process is a family (a t ) t 0 of free variables. It is said to be adapted if a t A t for all t 0. The free additive and free multiplicative Brownian motion.
36 Examples A free process is a family (a t ) t 0 of free variables. It is said to be adapted if a t A t for all t 0. The free additive and free multiplicative Brownian motion. The free Wishart process (Donati-Capitaine).
37 Examples A free process is a family (a t ) t 0 of free variables. It is said to be adapted if a t A t for all t 0. The free additive and free multiplicative Brownian motion. The free Wishart process (Donati-Capitaine). The free Jacobi process (Demni).
38 Asymptotic freeness and random matrices Definitions: A family (U i (m)) i I of random matrices converges in distribution to (U i ) i I in (A, φ) if lim m E(tr m(u i1 (m)... U ip (m))) = φ(u i1... U ip ) for any collection i 1,..., i p I. A family is said to be asymptotically free if it cv in distribution to free variables. Connection with random matrices : Voiculescu result on GUE.
39 The A t -free additive Brownian motion X It is a adapted and selfadjoint process such that : X 0 = 0.
40 The A t -free additive Brownian motion X It is a adapted and selfadjoint process such that : X 0 = 0. X t has the semi-circle law of mean 0 and variance t given by : σ t (dy) = 1 4t y 2πt 2 1 { 2 t,2 t} (y)dy
41 The A t -free additive Brownian motion X It is a adapted and selfadjoint process such that : X 0 = 0. X t has the semi-circle law of mean 0 and variance t given by : σ t (dy) = 1 4t y 2πt 2 1 { 2 t,2 t} (y)dy For any collection t 0 < t 1... < t k, X t0, X t1 X t0,..., X tk X tk 1 are free and stationnary.
42 Matrix-valued and free processes Let X m be a normalized Hermitian Brownian matrix : Xij m (t) = B m ii (t) m if i = j B m ij (t)+ 1 B m ij (t) 2m if i < j where B m, B m are two independent m m real Brownian matrices. Voiculescu result : X m additive free Brownian motion X. Corollary: Let Z m = (B m + 1B m )/ 2m be a non-selfadjoint process. Z m complex free Brownian motion Z defined by Z = (X 1 + 1X 2 )/ 2 where X 1, X 2 are free A t -free Brownian motions.
43 The A t -free multiplicative Brownian motion Y It is a adapted unitary process such that : Y 0 = I Biane showed that this process is the limit in distribution of the m m unitary Brownian motion Y (m).
44 The A t -free multiplicative Brownian motion Y It is a adapted unitary process such that : Y 0 = I ν t, the law of Y t is supported in the unit circle and is given by its Σ-transform (Bercovici et Voiculescu) : Σ νt (z) = e t 2 1+z 1 z ν t+s = ν t ν s Biane showed that this process is the limit in distribution of the m m unitary Brownian motion Y (m).
45 The A t -free multiplicative Brownian motion Y It is a adapted unitary process such that : Y 0 = I ν t, the law of Y t is supported in the unit circle and is given by its Σ-transform (Bercovici et Voiculescu) : Σ νt (z) = e t 2 1+z 1 z ν t+s = ν t ν s For any collection t 0 < t 1... < t k, Y t0, Y t1 Yt 1 0,..., Y tk Yt 1 k 1 are free and stationnary. Biane showed that this process is the limit in distribution of the m m unitary Brownian motion Y (m).
46 Free SDE For suitable parameters : Free Wishart process (Capitaine-Donati) : dw t = W t dz t + dzt Wt + λpdt Thanks.
47 Free SDE For suitable parameters : Free Wishart process (Capitaine-Donati) : dw t = W t dz t + dz t Wt + λpdt Free multiplicative BM (Kummerer-Speicher) Let X be a free additive BM. Then dy t = i dx t Y t 1 2 Y tdt, Y 0 = 1 Thanks.
48 Free SDE For suitable parameters : Free Wishart process (Capitaine-Donati) : dw t = W t dz t + dz t Wt + λpdt Free multiplicative BM (Kummerer-Speicher) Let X be a free additive BM. Then dy t = i dx t Y t 1 2 Y tdt, Y 0 = 1 Free Jacobi process (Demni) : dj t = λθ P J t dz t Jt + λθ J t dz t P Jt +(θp J t ) dt where Z is a complex free BM. Thanks.
49 Free SDE For suitable parameters : Free Wishart process (Capitaine-Donati) : dw t = W t dz t + dz t Wt + λpdt Free multiplicative BM (Kummerer-Speicher) Let X be a free additive BM. Then dy t = i dx t Y t 1 2 Y tdt, Y 0 = 1 Free Jacobi process (Demni) : dj t = λθ P J t dz t Jt + λθ J t dz t P Jt +(θp J t ) dt where Z is a complex free BM. condition : Injectivity. Thanks.
NONCOLLIDING BROWNIAN MOTIONS AND HARISH-CHANDRA FORMULA
Elect. Comm. in Probab. 8 (003) ELECTRONIC COMMUNICATIONS in PROBABILITY NONCOLLIDING BROWNIAN MOTIONS AND HARISH-CHANDRA FORMULA MAKOTO KATORI Department of Physics, Faculty of Science and Engineering,
More informationThe norm of polynomials in large random matrices
The norm of polynomials in large random matrices Camille Mâle, École Normale Supérieure de Lyon, Ph.D. Student under the direction of Alice Guionnet. with a significant contribution by Dimitri Shlyakhtenko.
More informationLecture III: Applications of Voiculescu s random matrix model to operator algebras
Lecture III: Applications of Voiculescu s random matrix model to operator algebras Steen Thorbjørnsen, University of Aarhus Voiculescu s Random Matrix Model Theorem [Voiculescu]. For each n in N, let X
More informationFree Probability Theory and Random Matrices. Roland Speicher Queen s University Kingston, Canada
Free Probability Theory and Random Matrices Roland Speicher Queen s University Kingston, Canada We are interested in the limiting eigenvalue distribution of N N random matrices for N. Usually, large N
More informationStrong Markov property of determinantal processes associated with extended kernels
Strong Markov property of determinantal processes associated with extended kernels Hideki Tanemura Chiba university (Chiba, Japan) (November 22, 2013) Hideki Tanemura (Chiba univ.) () Markov process (November
More informationSecond Order Freeness and Random Orthogonal Matrices
Second Order Freeness and Random Orthogonal Matrices Jamie Mingo (Queen s University) (joint work with Mihai Popa and Emily Redelmeier) AMS San Diego Meeting, January 11, 2013 1 / 15 Random Matrices X
More informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. MATRIX ANAL. APPL. Vol. XX, No. X, pp. XX XX c 005 Society for Industrial and Applied Mathematics DISTRIBUTIONS OF THE EXTREME EIGENVALUES OF THE COMPLEX JACOBI RANDOM MATRIX ENSEMBLE PLAMEN KOEV
More informationProgress in the method of Ghosts and Shadows for Beta Ensembles
Progress in the method of Ghosts and Shadows for Beta Ensembles Alan Edelman (MIT) Alex Dubbs (MIT) and Plamen Koev (SJS) Oct 8, 2012 1/47 Wishart Matrices (arbitrary covariance) G=mxn matrix of Gaussians
More informationHERMITE AND LAGUERRE POLYNOMIALS AND MATRIX- VALUED STOCHASTIC PROCESSES
Elect. Comm. in Probab. 13 8, 67 8 ELECTRONIC COMMUNICATIONS in PROBABILITY HERMITE AND LAGUERRE POLYNOMIALS AND MATRIX- VALUED STOCHASTIC PROCESSES STEPHAN LAWI 1 KBC Financial Products Hong Kong Limited
More informationMOMENTS OF THE HERMITIAN JACOBI PROCESS
MOMENTS OF THE HERMITIAN JACOBI PROCESS LUC DELEAVAL AND NIZAR DEMNI Abstract. In this paper, we compute the expectation of (Jt n), n, where (Jt) t 0 is the hermitian Jacobi process of size m n. To proceed,
More informationEigenvalues of the Laguerre Process as Non-colliding Squared Bessel Processes
Eigenvalues of the Laguerre Process as Non-colliding Squared Bessel Processes Wolfgang Konig, Neil O Connell Basic Research Institute in the Mathematical Science HP Laboratories Bristol HPL-BRIMS-200-07
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationStochastic Differential Equations Related to Soft-Edge Scaling Limit
Stochastic Differential Equations Related to Soft-Edge Scaling Limit Hideki Tanemura Chiba univ. (Japan) joint work with Hirofumi Osada (Kyushu Unv.) 2012 March 29 Hideki Tanemura (Chiba univ.) () SDEs
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationOperator norm convergence for sequence of matrices and application to QIT
Operator norm convergence for sequence of matrices and application to QIT Benoît Collins University of Ottawa & AIMR, Tohoku University Cambridge, INI, October 15, 2013 Overview Overview Plan: 1. Norm
More informationEigenvalues and Singular Values of Random Matrices: A Tutorial Introduction
Random Matrix Theory and its applications to Statistics and Wireless Communications Eigenvalues and Singular Values of Random Matrices: A Tutorial Introduction Sergio Verdú Princeton University National
More informationMatrix-valued stochastic processes
Matrix-valued stochastic processes and applications Małgorzata Snarska (Cracow University of Economics) Grodek, February 2017 M. Snarska (UEK) Matrix Diffusion Grodek, February 2017 1 / 18 Outline 1 Introduction
More informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationUpper triangular forms for some classes of infinite dimensional operators
Upper triangular forms for some classes of infinite dimensional operators Ken Dykema, 1 Fedor Sukochev, 2 Dmitriy Zanin 2 1 Department of Mathematics Texas A&M University College Station, TX, USA. 2 School
More informationFrom random matrices to free groups, through non-crossing partitions. Michael Anshelevich
From random matrices to free groups, through non-crossing partitions Michael Anshelevich March 4, 22 RANDOM MATRICES For each N, A (N), B (N) = independent N N symmetric Gaussian random matrices, i.e.
More informationLimit Laws for Random Matrices from Traffic Probability
Limit Laws for Random Matrices from Traffic Probability arxiv:1601.02188 Slides available at math.berkeley.edu/ bensonau Benson Au UC Berkeley May 9th, 2016 Benson Au (UC Berkeley) Random Matrices from
More informationPainlevé Representations for Distribution Functions for Next-Largest, Next-Next-Largest, etc., Eigenvalues of GOE, GUE and GSE
Painlevé Representations for Distribution Functions for Next-Largest, Next-Next-Largest, etc., Eigenvalues of GOE, GUE and GSE Craig A. Tracy UC Davis RHPIA 2005 SISSA, Trieste 1 Figure 1: Paul Painlevé,
More informationFree Probability Theory and Non-crossing Partitions. Roland Speicher Queen s University Kingston, Canada
Free Probability Theory and Non-crossing Partitions Roland Speicher Queen s University Kingston, Canada Freeness Definition [Voiculescu 1985]: Let (A, ϕ) be a non-commutative probability space, i.e. A
More informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationBurgers equation in the complex plane. Govind Menon Division of Applied Mathematics Brown University
Burgers equation in the complex plane Govind Menon Division of Applied Mathematics Brown University What this talk contains Interesting instances of the appearance of Burgers equation in the complex plane
More informationMOMENTS OF THE HERMITIAN MATRIX JACOBI PROCESS
MOMENTS OF THE HERMITIAN MATRIX JACOBI PROCESS LUC DELEAVAL AND NIZAR DEMNI Abstract In this paper, we compute the expectation of traces of powers of the hermitian matrix Jacobi process for a large enough
More informationProperties of an infinite dimensional EDS system : the Muller s ratchet
Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population
More informationStrong Markov property of determinatal processes
Strong Markov property of determinatal processes Hideki Tanemura Chiba university (Chiba, Japan) (August 2, 2013) Hideki Tanemura (Chiba univ.) () Markov process (August 2, 2013) 1 / 27 Introduction The
More informationFluctuations of Random Matrices and Second Order Freeness
Fluctuations of Random Matrices and Second Order Freeness james mingo with b. collins p. śniady r. speicher SEA 06 Workshop Massachusetts Institute of Technology July 9-14, 2006 1 0.4 0.2 0-2 -1 0 1 2-2
More informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
More informationDifferential Equations for Dyson Processes
Differential Equations for Dyson Processes Joint work with Harold Widom I. Overview We call Dyson process any invariant process on ensembles of matrices in which the entries undergo diffusion. Dyson Brownian
More informationFREE PROBABILITY THEORY
FREE PROBABILITY THEORY ROLAND SPEICHER Lecture 4 Applications of Freeness to Operator Algebras Now we want to see what kind of information the idea can yield that free group factors can be realized by
More informationLINEAR PRESERVER PROBLEMS: generalized inverse
LINEAR PRESERVER PROBLEMS: generalized inverse Université Lille 1, France Banach Algebras 2011, Waterloo August 3-10, 2011 I. Introduction Linear preserver problems is an active research area in Matrix,
More informationRECTANGULAR RANDOM MATRICES, ENTROPY, AND FISHER S INFORMATION
J. OPERATOR THEORY 00:0(XXXX), 00 00 Copyright by THETA, XXXX RECTANGULAR RANDOM MATRICES, ENTROPY, AND FISHER S INFORMATION FLORENT BENAYCH-GEORGES Communicated by Serban Stratila ABSTRACT. We prove that
More informationFredholm determinant with the confluent hypergeometric kernel
Fredholm determinant with the confluent hypergeometric kernel J. Vasylevska joint work with I. Krasovsky Brunel University Dec 19, 2008 Problem Statement Let K s be the integral operator acting on L 2
More informationGeometric Dyson Brownian motion and May Wigner stability
Geometric Dyson Brownian motion and May Wigner stability Jesper R. Ipsen University of Melbourne Summer school: Randomness in Physics and Mathematics 2 Bielefeld 2016 joint work with Henning Schomerus
More information1 Intro to RMT (Gene)
M705 Spring 2013 Summary for Week 2 1 Intro to RMT (Gene) (Also see the Anderson - Guionnet - Zeitouni book, pp.6-11(?) ) We start with two independent families of R.V.s, {Z i,j } 1 i
More informationFree probability and quantum information
Free probability and quantum information Benoît Collins WPI-AIMR, Tohoku University & University of Ottawa Tokyo, Nov 8, 2013 Overview Overview Plan: 1. Quantum Information theory: the additivity problem
More informationLecture I: Asymptotics for large GUE random matrices
Lecture I: Asymptotics for large GUE random matrices Steen Thorbjørnsen, University of Aarhus andom Matrices Definition. Let (Ω, F, P) be a probability space and let n be a positive integer. Then a random
More informationMaximal height of non-intersecting Brownian motions
Maximal height of non-intersecting Brownian motions G. Schehr Laboratoire de Physique Théorique et Modèles Statistiques CNRS-Université Paris Sud-XI, Orsay Collaborators: A. Comtet (LPTMS, Orsay) P. J.
More informationFree Meixner distributions and random matrices
Free Meixner distributions and random matrices Michael Anshelevich July 13, 2006 Some common distributions first... 1 Gaussian Negative binomial Gamma Pascal chi-square geometric exponential 1 2πt e x2
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationAiry and Pearcey Processes
Airy and Pearcey Processes Craig A. Tracy UC Davis Probability, Geometry and Integrable Systems MSRI December 2005 1 Probability Space: (Ω, Pr, F): Random Matrix Models Gaussian Orthogonal Ensemble (GOE,
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationVerona Course April Lecture 1. Review of probability
Verona Course April 215. Lecture 1. Review of probability Viorel Barbu Al.I. Cuza University of Iaşi and the Romanian Academy A probability space is a triple (Ω, F, P) where Ω is an abstract set, F is
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationCombinatorial Aspects of Free Probability and Free Stochastic Calculus
Combinatorial Aspects of Free Probability and Free Stochastic Calculus Roland Speicher Saarland University Saarbrücken, Germany supported by ERC Advanced Grant Non-Commutative Distributions in Free Probability
More informationLectures on the Combinatorics of Free Probability Theory. Alexandru Nica Roland Speicher
Lectures on the Combinatorics of Free Probability Theory Alexandru Nica Roland Speicher Department of Pure Mathematics, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada E-mail address: anica@math.uwaterloo.ca
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationHomogenization of the Dyson Brownian Motion
Homogenization of the Dyson Brownian Motion P. Bourgade, joint work with L. Erdős, J. Yin, H.-T. Yau Cincinnati symposium on probability theory and applications, September 2014 Introduction...........
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationGaussian Free Field in beta ensembles and random surfaces. Alexei Borodin
Gaussian Free Field in beta ensembles and random surfaces Alexei Borodin Corners of random matrices and GFF spectra height function liquid region Theorem As Unscaled fluctuations Gaussian (massless) Free
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationBulk scaling limits, open questions
Bulk scaling limits, open questions Based on: Continuum limits of random matrices and the Brownian carousel B. Valkó, B. Virág. Inventiones (2009). Eigenvalue statistics for CMV matrices: from Poisson
More informationQuantum Symmetric States
Quantum Symmetric States Ken Dykema Department of Mathematics Texas A&M University College Station, TX, USA. Free Probability and the Large N limit IV, Berkeley, March 2014 [DK] K. Dykema, C. Köstler,
More informationAn Overview of the Martingale Representation Theorem
An Overview of the Martingale Representation Theorem Nuno Azevedo CEMAPRE - ISEG - UTL nazevedo@iseg.utl.pt September 3, 21 Nuno Azevedo (CEMAPRE - ISEG - UTL) LXDS Seminar September 3, 21 1 / 25 Background
More informationEigenvalue PDFs. Peter Forrester, M&S, University of Melbourne
Outline Eigenvalue PDFs Peter Forrester, M&S, University of Melbourne Hermitian matrices with real, complex or real quaternion elements Circular ensembles and classical groups Products of random matrices
More informationTHE NONCOMMUTATIVE TORUS
THE NONCOMMUTATIVE TORUS The noncommutative torus as a twisted convolution An ordinary two-torus T 2 with coordinate functions given by where x 1, x 2 [0, 1]. U 1 = e 2πix 1, U 2 = e 2πix 2, (1) By Fourier
More informationStochastic Calculus. Kevin Sinclair. August 2, 2016
Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationGeometric measure on quasi-fuchsian groups via singular traces. Part I.
Geometric measure on quasi-fuchsian groups via singular traces. Part I. Dmitriy Zanin (in collaboration with A. Connes and F. Sukochev) University of New South Wales September 19, 2018 Dmitriy Zanin Geometric
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationFree Probability Theory
Noname manuscript No. (will be inserted by the editor) Free Probability Theory and its avatars in representation theory, random matrices, and operator algebras; also featuring: non-commutative distributions
More informationGaussian Free Field in (self-adjoint) random matrices and random surfaces. Alexei Borodin
Gaussian Free Field in (self-adjoint) random matrices and random surfaces Alexei Borodin Corners of random matrices and GFF spectra height function liquid region Theorem As Unscaled fluctuations Gaussian
More informationTriangular matrices and biorthogonal ensembles
/26 Triangular matrices and biorthogonal ensembles Dimitris Cheliotis Department of Mathematics University of Athens UK Easter Probability Meeting April 8, 206 2/26 Special densities on R n Example. n
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationLog-Convexity Properties of Schur Functions and Generalized Hypergeometric Functions of Matrix Argument. Donald St. P. Richards.
Log-Convexity Properties of Schur Functions and Generalized Hypergeometric Functions of Matrix Argument Donald St. P. Richards August 22, 2009 Abstract We establish a positivity property for the difference
More informationFree Probability, Sample Covariance Matrices and Stochastic Eigen-Inference
Free Probability, Sample Covariance Matrices and Stochastic Eigen-Inference Alan Edelman Department of Mathematics, Computer Science and AI Laboratories. E-mail: edelman@math.mit.edu N. Raj Rao Deparment
More informationConvolutions and fluctuations: free, finite, quantized.
Convolutions and fluctuations: free, finite, quantized. Vadim Gorin MIT (Cambridge) and IITP (Moscow) February 2018 Hermitian matrix operations a 1 0 b 1 0 0 a 2 0 0 b 2 0 A =. 0.. B = 0. 0.. 0 0 a N 0
More informationMultiple orthogonal polynomials. Bessel weights
for modified Bessel weights KU Leuven, Belgium Madison WI, December 7, 2013 Classical orthogonal polynomials The (very) classical orthogonal polynomials are those of Jacobi, Laguerre and Hermite. Classical
More information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationVector valued Jack polynomials. Jean-Gabriel Luque
Vector valued Jack polynomials Jean-Gabriel Luque What are Jack polynomials? The story begins in 1960's... Inner product: Orthogonal polynomials How to define multivariate orthogonal polynomials? Interpretation
More informationExtreme eigenvalue fluctutations for GUE
Extreme eigenvalue fluctutations for GUE C. Donati-Martin 204 Program Women and Mathematics, IAS Introduction andom matrices were introduced in multivariate statistics, in the thirties by Wishart [Wis]
More informationCAAM 335 Matrix Analysis
CAAM 335 Matrix Analysis Solutions to Homework 8 Problem (5+5+5=5 points The partial fraction expansion of the resolvent for the matrix B = is given by (si B = s } {{ } =P + s + } {{ } =P + (s (5 points
More informationFree Probability and Random Matrices: from isomorphisms to universality
Free Probability and Random Matrices: from isomorphisms to universality Alice Guionnet MIT TexAMP, November 21, 2014 Joint works with F. Bekerman, Y. Dabrowski, A.Figalli, E. Maurel-Segala, J. Novak, D.
More informationAppearance of determinants for stochastic growth models
Appearance of determinants for stochastic growth models T. Sasamoto 19 Jun 2015 @GGI 1 0. Free fermion and non free fermion models From discussions yesterday after the talk by Sanjay Ramassamy Non free
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationApplications and fundamental results on random Vandermon
Applications and fundamental results on random Vandermonde matrices May 2008 Some important concepts from classical probability Random variables are functions (i.e. they commute w.r.t. multiplication)
More informationLax Solution Part 4. October 27, 2016
Lax Solution Part 4 www.mathtuition88.com October 27, 2016 Textbook: Functional Analysis by Peter D. Lax Exercises: Ch 16: Q2 4. Ch 21: Q1, 2, 9, 10. Ch 28: 1, 5, 9, 10. 1 Chapter 16 Exercise 2 Let h =
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part II Manuela Girotti based on M. Girotti s PhD thesis, A. Kuijlaars notes from Les Houches Winter School 202 and B. Eynard s notes
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationConvergence of spectral measures and eigenvalue rigidity
Convergence of spectral measures and eigenvalue rigidity Elizabeth Meckes Case Western Reserve University ICERM, March 1, 2018 Macroscopic scale: the empirical spectral measure Macroscopic scale: the empirical
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationPath Decomposition of Markov Processes. Götz Kersting. University of Frankfurt/Main
Path Decomposition of Markov Processes Götz Kersting University of Frankfurt/Main joint work with Kaya Memisoglu, Jim Pitman 1 A Brownian path with positive drift 50 40 30 20 10 0 0 200 400 600 800 1000-10
More informationNumerical analysis and random matrix theory. Tom Trogdon UC Irvine
Numerical analysis and random matrix theory Tom Trogdon ttrogdon@math.uci.edu UC Irvine Acknowledgements This is joint work with: Percy Deift Govind Menon Sheehan Olver Raj Rao Numerical analysis and random
More informationAlgebraic reformulation of Connes embedding problem and the free group algebra.
Algebraic reformulation of Connes embedding problem and the free group algebra. Kate Juschenko, Stanislav Popovych Abstract We give a modification of I. Klep and M. Schweighofer algebraic reformulation
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationLyapunov exponents of free operators
Journal of Functional Analysis 255 (28) 1874 1888 www.elsevier.com/locate/jfa Lyapunov exponents of free operators Vladislav Kargin Department of Mathematics, Stanford University, Stanford, CA 9435, USA
More informationLecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 15 Review of Matrix Theory III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Matrix An m n matrix is a rectangular or square array of
More informationGeometric RSK, Whittaker functions and random polymers
Geometric RSK, Whittaker functions and random polymers Neil O Connell University of Warwick Advances in Probability: Integrability, Universality and Beyond Oxford, September 29, 2014 Collaborators: I.
More informationCommutative Banach algebras 79
8. Commutative Banach algebras In this chapter, we analyze commutative Banach algebras in greater detail. So we always assume that xy = yx for all x, y A here. Definition 8.1. Let A be a (commutative)
More informationParticle models for Wasserstein type diffusion
Particle models for Wasserstein type diffusion Vitalii Konarovskyi University of Leipzig Bonn, 216 Joint work with Max von Renesse konarovskyi@gmail.com Wasserstein diffusion and Varadhan s formula (von
More information