Tracy-Widom Law: Fredholm determinants and Airy functions

This is a continuation of the previous post available here. This section is quite different from the previous post. Here we present a series of propositions. The proofs of the propositions are not present here. I may include them later. However the proofs of the propositions (except proposition 4) are not necessary in our proof.

In the previous post, we introduced Hermite polynomials and wave functions. We noticed the Hermite polynomials are orthogonal polynomials. There are a plethora of results involving orthonormal polynomials. We merely state only few of them that we will be needing in the our proof.

Proposition 1: Let H_n(x) be the Hermite polynomial of degree n. Let \psi_n(x) be the corresponding wave function. Let \displaystyle K^{(n)}(x,y)=\sum_{k=0}^{n-1} \psi_k(x)\psi_k(y). We have the following results.

(a) H_n(x+t)=\sum_{k=0}^{n} \binom{n}{k}H_k(x)t^{n-k}

(b) \displaystyle \frac1{\sqrt{n}}K^{(n)}(x,y)= \frac{\psi_n(x)\psi_{n-1}(y)-\psi_{n-1}(x)\psi_n(y)}{x-y}

(c) \displaystyle K^{(n)}(x,y)=\frac{\psi_n(x)\psi_n'(y)-\psi_n'(x)\psi_n(y)}{x-y}-\frac12\psi_n(x)\psi_n(y)

(d) \displaystyle \frac{d}{dx}K^{(n)}(x,x)=-\sqrt{n}\psi_n(x)\psi_{n-1}(x)

Proof: Not included

Let us recall Theorem 3 once again.

Theorem 3: For any measurable subset A of \mathbb{R},

\displaystyle P\left(\bigcap_{i=1}^n \{\lambda_i \in A\}\right) =1+\sum_{k=1}^{\infty} \frac{(-1)^k}{k!}\int_{A^c}\cdots\int_{A^c} \det_{i,j=1}^k K^{(n)}(x_i,x_j)\prod_{i=1}^k dx_i

Let us replace A by A^c and rewrite the above result as

\displaystyle P\left(\bigcap_{i=1}^n \{\lambda_i \not\in A\}\right) =1+\sum_{k=1}^{\infty} \frac{(-1)^k}{k!}\int_{A}\cdots\int_{A} \det_{i,j=1}^k K^{(n)}(x_i,x_j)\prod_{i=1}^k dx_i

Now it is right time to see what exactly are we trying to show if we want to prove weak convergence. For that we need some basic idea about Airy function.

Continue reading


Tracy-Widom Law: The starting point

Hi there!

Introduction: RMT deals with matrices with random variables as its entries. RMT is a fascinating field which is expanding rapidly in recent years. In this series of posts, I will be proving the vague convergence of the celebrated Tracy-Widom Law. I will be mostly following the book ‘An Introduction to Random Matrices’ by Anderson


Let T,W be independent N(0,1/2) variables. Define Z=T+iW to be a complex random variable. We say that Z follows standard complex normal distribution and denote it as Z \sim CN(0,1), since E(Z)=0 and E(Z\bar{Z})=1. Let U be a n\times n random matrix with iid CN(0,1)  entries and V be a n\times n random matrix with iid N(0,1)  entries. Define

\displaystyle  X=\frac{V+V^T}{\sqrt{2}},\ \ Y=\frac{U+U^*}{\sqrt2}.

Let P_n^{(1)} denote the law of X.  The random matrix X is said to belong to the Gaussian orthogonal ensemble [\mathcal{H}_n^{(1)}] (GOE(n)). Observe that X is symmetric and X_{1,1} is a N(0,2) random variable while X_{1,2} is N(0,1) random variable.

Let P_n^{(2)} denote the law of Y. Then the random matrix Y is said to belong to the Gaussian unitary ensemble [\mathcal{H}_n^{(2)}] (GUE(n)). Observe that Y is hermitian and Y_{1,1} is a N(0,1) random variable while Y_{1,2} is CN(0,1) random variable.

Why the law P_n^{(\beta)} special?

Because the density at H \in \mathcal{H}_N^{(\beta)} is proportional to \exp\left(-\mbox{tr}\frac{\beta H^2}{4}\right). Thus for the law P_n^{(1)}, the density at H and OHO^T is same for any orthogonal matrix O and for law P_n^{(2)}, the density at H and UHU^{*} is same for any unitary matrix U (and hence the names GOE and GUE). This also enables us to find the joint distribution of the eigenvalues exactly!

Target: Let \lambda_1^n\le \lambda_2^n \le \cdots \le \lambda_n^n be the eigenvalues of a matrix X from GUE/GOE. Then it is known that the empirical distribution of the eigenvalues L_n:=\frac1n\sum_{i=1}^n \delta_{\lambda_i^n/sqrt{n}} converges weakly, in probability to the semicircle distribution \sigma(x)\,dx with density \displaystyle \sigma(x) =\frac1{2\pi}\sqrt{4-x^2}\mathbb{I}_{|x|\le 2}. We also have that \displaystyle \frac1{\sqrt{n}}\lambda_n^n \stackrel{a.e.}{\to} 2. What we are interested in the limiting distribution of this largest eigenvalue. In fact what we intend to show is that \displaystyle n^{2/3}\left(\frac{\lambda_n^n}{\sqrt{n}}-2\right) converges weakly. In the initial series of posts we will be only showing the vague convergence of above expression.

Let us adapt this convention: If we are using the notation \lambda_1^n, \ldots \lambda_n^n for eigenvalues then they are ordered. If we use \lambda_1, \ldots, \lambda_n for eigenvalues then they are unordered.

Notations: Let A=((a_{ij}))_{i,j=1}^n be a square matrix of order n. The determinant of A will be denoted as

|A|=\det(A)=\det\limits_{i,j=1}^n a_{ij}

By laplace expansion we also have the following formula for the determinants.

\det\limits_{i,j=1}^n a_{i,j}=\sum\limits_{\sigma \in S_n} \epsilon_{\sigma} \prod\limits_{i=1}^{n} a_{i,\sigma(i)}

Theorem 1: Let X\in \mathcal{H}_n^{(\beta)} be a random matrix with law P_n^{(\beta)} (\beta=1,2). Let \lambda_1^n \le \lambda_2^n \le \ldots \le \lambda_n^n denotes the ordered eigenvalues of the random matrix X. The joint distribution of the ordered eigenvalues has density with respect to Lebesgue measure which equals

\displaystyle \hat\rho_{n}^{(\beta)}(x_1,x_2,\ldots,x_n):=n!C_n^{(\beta)}|\Delta(x)|^{\beta}\prod_{i=1}^n e^{-\beta x_i^2/4}\mathbb{I}_{x_1\le x_2\le \ldots \le x_n}

where C_n^{(\beta)} is the normalizing constant and \Delta(x) is Vandermonde determinant associated with x=(x_1,x_2,\ldots,x_n) which equals to \displaystyle \det_{i,j=1}^n x_i^{j-1}=\prod_{1\le i<j\le n} (x_j-x_i).

Proof: The proof will be skipped. Essentially the proof involves the change of variable formula. It takes a certain suitable transformation with a smooth inverse. Then we obtain the required density by integrating out the extra variables. The \Delta(x) part emerges as  the jacobian of the transformation and e^{-\beta x_i^2/4} appears by making the change of variables. However to justify it clearly, lot more mathematics is needed.

Remark: What we are interested in, is the distribution of the unordered eigenvalues which is obtained by symmetrizing the density \hat\rho_n^{(\beta)}. We say the joint density of the unordered eigenvalues is given by

\displaystyle \rho_n^{(\beta)}:=C_n^{(\beta)}|\Delta(x)|^{\beta}\prod_{i=1}^n e^{-\beta x_i^2/4}

Technically speaking, this is just a symmetrized valid density. This should not be called as the joint density of the unordered eigenvalues as integrating this joint density n-1 times gives us a marginal density. But that marginal density is the density of which eigenvalue? The best we can say is that if we choose one eigenvalue randomly, it is the density of that randomly chosen eigenvalue.

Continue reading