Tracy–Widom Distribution
Tracy–Widom distribution is the limiting distribution of the scaled largest eigenvalue of a random matrix belonging to Gaussian ensembles. It also appears in various different disciplines, such as counting problems, random growth models, phase transitions, etc., and gives accurate predictions.
Tracy–Widom distribution consists of three classes ,
, and
. Each corresponds to different Gaussian ensembles; see the corresponding PDFs.
![](assets.en/tracy-widom-distribution/O_38.png)
Use MatrixPropertyDistribution to represent the scaled largest eigenvalue of a matrix from GUE.
![Click for copyable input](assets.en/tracy-widom-distribution/In_80.png)
ev\[ScriptCapitalD][2, n_] :=
MatrixPropertyDistribution[(Max[Eigenvalues[x]] - 2 Sqrt[n]) n^(1/6),
x \[Distributed] GaussianUnitaryMatrixDistribution[n]]
Sample from the distribution and compare the histogram with the PDF.
![Click for copyable input](assets.en/tracy-widom-distribution/In_81.png)
sample = RandomVariate[ev\[ScriptCapitalD][2, 250], 2000];
![](assets.en/tracy-widom-distribution/O_39.png)
Tracy–Widom distribution can be well approximated by gamma distribution in the central region.
![Click for copyable input](assets.en/tracy-widom-distribution/In_83.png)
gdist = GammaDistribution[k, s, 1, a];
PDF[gdist, x]
![](assets.en/tracy-widom-distribution/O_40.png)
Fit gamma distribution with Tracy–Widom distribution of by matching the first three moments.
![Click for copyable input](assets.en/tracy-widom-distribution/In_84.png)
moments = Through[{Mean, Variance, Skewness}[gdist]];
nmoments =
N[Through[{Mean, Variance, Skewness}[TracyWidomDistribution[1]]]];
sol = FindRoot[Thread[moments == nmoments], {{k, 1}, {s, 1}, {a, 1}}]
![](assets.en/tracy-widom-distribution/O_41.png)
Compare the PDFs.
![](assets.en/tracy-widom-distribution/O_42.png)