The Instability of Painlevé Equations in Recovering Largest Eigenvalue Distributions of GUE, LUE, JUE and an Attempt of Solution to It
Fuente:
arXiv
Guardado en:
| Autor principal: | |
|---|---|
| Formato: | Preprint |
| Publicado: |
2025
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866915684821237760 |
|---|---|
| author | Gu, Haonan |
| author_facet | Gu, Haonan |
| contents | The distribution of the largest eigenvalue for the three classical unitary ensembles -- GUE, LUE, and JUE -- admits two complementary exact descriptions: (i) as Fredholm determinants of their orthogonal polynomial correlation kernels and (ii) as isomonodromic $τ$-functions governed by Painlevé equations. For finite $n$, the associated Jimbo-Miwa-Okamoto $σ$-forms are $\PIV$ (GUE), $\mathrm{PV}$ (LUE), and $\PVI$ (JUE); under soft- or hard-edge scalings these degenerate to $\PII$ or $\PIIIp$ descriptions of the Tracy-Widom and hard-edge laws \cite{tracy1994level,forrester2003painleve,deift1999orthogonal}.
It is well known among random matrix theorists (for example Folkmar Bornemann) that the Fredholm determinant is a more numerically stable and accurate way to compute the CDF of the largest eigenvalue for GUE, LUE, JUE than direct Painlevé integration. The aim of this paper is not to improve on Fredholm methods, but to see to what extent one can numerically recover the \emph{correct} Painlevé solution from finite-$n$ data and how unstable this reconstruction is. Numerically, we verify the equality between the Fredholm- and Painlevé-based CDFs by combining (a) high-accuracy Nyström discretizations of the finite-$n$ Fredholm determinants \cite{bornemann2010numerical} with (b) an anchored, branch-locked integration of the $σ$-form ODEs, where anchors are extracted from local least-squares fits to $\log\det(I-\mathsf K)$. Our results confirm agreement across GUE/LUE/JUE with precision of $O(10^{-3})$ to $O(10^{-5})$ (occasionally $O(10^{-2})$) and illustrate the finite-$n$ to scaling-limit transition. The theoretical connections to $τ$-functions and Virasoro constraints follow the framework of \cite{adler2000random,forrester2003painleve} |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2512_08178 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | The Instability of Painlevé Equations in Recovering Largest Eigenvalue Distributions of GUE, LUE, JUE and an Attempt of Solution to It Gu, Haonan Numerical Analysis Probability The distribution of the largest eigenvalue for the three classical unitary ensembles -- GUE, LUE, and JUE -- admits two complementary exact descriptions: (i) as Fredholm determinants of their orthogonal polynomial correlation kernels and (ii) as isomonodromic $τ$-functions governed by Painlevé equations. For finite $n$, the associated Jimbo-Miwa-Okamoto $σ$-forms are $\PIV$ (GUE), $\mathrm{PV}$ (LUE), and $\PVI$ (JUE); under soft- or hard-edge scalings these degenerate to $\PII$ or $\PIIIp$ descriptions of the Tracy-Widom and hard-edge laws \cite{tracy1994level,forrester2003painleve,deift1999orthogonal}. It is well known among random matrix theorists (for example Folkmar Bornemann) that the Fredholm determinant is a more numerically stable and accurate way to compute the CDF of the largest eigenvalue for GUE, LUE, JUE than direct Painlevé integration. The aim of this paper is not to improve on Fredholm methods, but to see to what extent one can numerically recover the \emph{correct} Painlevé solution from finite-$n$ data and how unstable this reconstruction is. Numerically, we verify the equality between the Fredholm- and Painlevé-based CDFs by combining (a) high-accuracy Nyström discretizations of the finite-$n$ Fredholm determinants \cite{bornemann2010numerical} with (b) an anchored, branch-locked integration of the $σ$-form ODEs, where anchors are extracted from local least-squares fits to $\log\det(I-\mathsf K)$. Our results confirm agreement across GUE/LUE/JUE with precision of $O(10^{-3})$ to $O(10^{-5})$ (occasionally $O(10^{-2})$) and illustrate the finite-$n$ to scaling-limit transition. The theoretical connections to $τ$-functions and Virasoro constraints follow the framework of \cite{adler2000random,forrester2003painleve} |
| title | The Instability of Painlevé Equations in Recovering Largest Eigenvalue Distributions of GUE, LUE, JUE and an Attempt of Solution to It |
| topic | Numerical Analysis Probability |
| url | https://arxiv.org/abs/2512.08178 |