Harold Widom’s work in random matrix theory

Abstract

This is a survey of Harold Widom’s work in random matrices. We start with his pioneering papers on the sine-kernel determinant, continue with his and Craig Tracy’s groundbreaking results concerning the distribution functions of random matrix theory, touch on the remarkable universality of the Tracy–Widom distributions in mathematics and physics, and close with Tracy and Widom’s remarkable work on the asymmetric simple exclusion process.

1. Introduction

The distributions of random matrix theory govern the statistical properties of a wide variety of large systems which do not obey the usual laws of classical probability. Such systems appear in many different areas of applied science and technology, including heavy nuclei, polymer growth, high-dimensional data analysis, and certain percolation processes. The four distribution functions play a particularly important role in the mathematical apparatus of random matrices. The first one describes the emptiness formation probability in the bulk of the spectrum of a large random matrix, and it is explicitly given in terms of the sine-kernel Fredholm determinant. The second, the third, and the forth distributions are given in terms of the Airy kernel Fredholm determinant, and they describe the edge fluctuations of the eigenvalues in the large size limit of the matrices taken from the three classical Gaussian ensembles, unitary (GUE), orthogonal (GOE), and symplectic (GSE). The last three of these distributions are known now as the Tracy–Widom distribution functions.

A key analytical observation concerning the distribution functions of random matrix theory, which was made on many occasions in the papers 567101112, is that they satisfy certain nonlinear integrable PDEs. This property, which generalizes the first result of this type established in JMMS for the sine-kernel determinant, follows in turn from a remarkable Fredholm determinant representation for the random matrix distributions. The existence of such representations in several important examples beyond the sine-kernel case was also first shown in the above-mentioned papers.

What follows is an overview of the principal contributions of Harold Widom to random matrix theory. We start with his papers 3 and 9 on the sine-kernel determinant, which brought mathematical rigor to the theory of this, the first universal distribution function of random matrix theory.

2. The sine-kernel determinant

Let $J = \bigcup _{k=0}^{n}J_k = \bigcup _{k=0}^{n}(a_k, b_k)$ be a union of $n+1$ disjoint intervals in $\mathbb{R}$. Consider the Fredholm determinant

$$\begin{equation*} P_x= \det \left( 1 - K_{\text{sine}}\right), \end{equation*}$$

where $K_{\text{sine}}$ is the trace class operator in $L_2\Bigl (J; dz\Bigr )$ with kernel

$$\begin{equation*} K_{\text{sine}}(z,z') = \frac{\sin x(z-z')}{\pi (z-z')}. \end{equation*}$$

The determinant $P_x$ plays a central role in random matrix theory. Indeed, it is the probability of finding no eigenvalues in the union of intervals $\frac{xJ}{\pi } \coloneq \bigcup _{k=0}^{n}\left(\frac{xa_k}{\pi }, \frac{xb_k}{\pi }\right)$ for a random Hermitian matrix chosen from the Gaussian unitary ensemble (GUE), in the bulk scaling limit with mean spacing 1 (see Meh and Figure 1). Moreover, in the one interval case, the second derivative of $P_x$ describes the distribution of normalized spacings of eigenvalues of large random GUE matrices (see, e.g., Dei1). The determinant $P_x$ also appears in quantum and statistical mechanics. For instance, it describes the emptiness formation probability in the one-dimensional impenetrable Bose gas Len and the gap probability in the one-dimensional Coulomb gas, at inverse temperature $\beta = 2$ (see, e.g., Dys1). The key analytical issue related to $P_x$ is its large $x$ behavior, i.e., the large gap asymptotics. Harold Widom made major contributions to the resolution of this question.

In the one interval case, $n =0$, after rescaling and translation, we may assume $J = (-1, 1)$. For this case, in 1973, des Cloizeaux and Mehta dCM showed that as $x\to \infty$,

$$\begin{equation} \ln P_x=-{x^2\over 2} -{1\over 4}\ln x+ c_0 +o(1), {\tag{2.1}} \end{equation}$$

for some constant $c_0$. In 1976 Dyson Dys2 showed that $P_x$ in fact has a full asymptotic expansion of the form

$$\begin{equation} \ln P_x=-{x^2\over 2} -{1\over 4}\ln x+ c_0 +{a_1\over x}+ {a_2\over x^2}+\cdots . {\tag{2.2}} \end{equation}$$

Dyson identified all the constants $c_0$, $a_1$, $a_2$, $\dots$. Of particular interest is the constant $c_0$, which he found to be

$$\begin{equation} c_0={1\over 12}\ln 2+ 3\zeta '(-1), {\tag{2.3}} \end{equation}$$

where $\zeta (z)$ is the Riemann zeta-function. It should be noted that Dyson obtained this result using one of the early results of Widom 1 on the asymptotics of Toeplitz determinants with symbols supported on circular arcs.

The results in dCM and Dys2 were not fully rigorous. In 3, using an adaptation of Szegő’s classical method to the continuous analogues of orthogonal polynomials (the so-called Krein functions), Widom gave the first rigorous proof of the leading asymptotics in 2.1 in the form

$$\begin{equation} \ln P_x=-{x^2\over 2}(1+o(1)) \tag{2.4} \end{equation}$$

as $x \to \infty$. Actually, Widom proved a slightly stronger result,

$$\begin{equation*} \frac{d}{dx}\ln P_x = -x + O(1) \end{equation*}$$

as $x \to \infty$. In addition, Widom computed the leading asymptotics of the quantity $\operatorname {tr}\left( (1- K_{\text{sine}})^{-1}K_{\text{sine}}\right)$, which is the ratio of the probability that there is at most one eigenvalue in the interval $\left(\frac{-x}{\pi }, \frac{x}{\pi }\right)$ to $P_x$. In the subsequent paper 9, Widom considered the multi-interval case, $n > 0$, and showed that as $x \to \infty$,

$$\begin{equation*} \frac{d}{dx}\ln P_x = c_1 x + c_2(x) + o(1), \end{equation*}$$

where $c_1$ is a negative constant and $c_2(x)$ is a certain bounded oscillatory function of $x$, which was described up to the solution of a Jacobi inversion problem.⁠¹ The method of 9 is a further development of the approach in 3. As in 3, Widom also computes the leading asymptotics of $\operatorname {tr}\left( (1- K_{\text{sine}})^{-1}K_{\text{sine}}\right)$.

¹

An explicit formula for $c_2(x)$, in terms of the Riemann theta-function, was obtained later in DIZ.

✖

The formula 2.3 for $c_0$ in Dys2 was in the form of conjecture: A rigorous proof of 2.3 was only given in 2004. In fact, two proofs of this formula were given independently in Kra and Ehr. It is remarkable that the proof of 2.3 in Kra again uses (in conjunction with the nonlinear steepest descent method) Widom’s computation in 1.

3. The Tracy–Widom distribution functions

The joint probability densities of the eigenvalues of random matrices from the GOE, GUE, and GSE ensembles are given by

$$\begin{equation} P_{N\beta }(\lambda _1, \dots , \lambda _N) = Z^{-1}_{N\beta }\,e^{-\frac{1}{2}\beta \sum _{j=1}^{N}\lambda ^{2}_j} \prod _{j<k}|\lambda _j - \lambda _k|^{\beta }, {\tag{3.1}} \end{equation}$$

where $Z_{N\beta }$ is a normalization constant (or partition function) and

$$\begin{equation*} \beta \coloneq \begin{cases} 1 & \text{for GOE},\\ 2 & \text{for GUE},\\ 4 & \text{for GSE}. \end{cases} \end{equation*}$$

The famous Tracy–Widom distribution functions, commonly denoted as $F_{\beta }(t)$, describe the edge fluctuation of the eigenvalues in the large $N$ limit. They are defined via the scaling limits,

$$\begin{equation} F_{\beta }(t) = \lim _{N\to \infty }\operatorname {Prob} \bigl ((\lambda _{\max } -\sqrt {2N})2^{1/2}N^{1/6} \leq t\bigl ), {\tag{3.2}} \end{equation}$$

where $\lambda _{\max }$ is the largest eigenvalue drawn from the ensembles with density 3.1. The central theme of the series of papers 5671012 is the following analytical description of these distributions, which links them to the theory of integrable systems.

Let $K_{\text{Airy}}$ be the trace-class operator in $L_2\bigl ((t, \infty ); dz\bigr )$ with kernel

$$\begin{equation} K_{\text{Airy}}(z, z') = \frac{\operatorname {Ai}(z)\operatorname {Ai}'(z') - \operatorname {Ai}'(z)\operatorname {Ai}(z')}{z - z'}, {\tag{3.3}} \end{equation}$$

where $\operatorname {Ai}(z)$ is the Airy function,

$$\begin{equation*} \operatorname {Ai}(z) = \frac{1}{2\pi }\int _{-\infty }^{\infty }e^{-\frac{i}{3}s^3 - izs}ds. \end{equation*}$$

Then, the $\beta =2$ Tracy–Widom distribution, $F_{2}(t)$, is given by the Airy-kernel Fredholm determinant (see For),

$$\begin{equation} F_{2}(t) = \det \left(1 - K_{\text{Airy}}\right). {\tag{3.4}} \end{equation}$$

Moreover, as shown in 6, the following formula is valid for the Fredholm determinant on the right-hand side of 3.4:

$$\begin{equation} \det (1 - K_{\text{Airy}}) = \exp \left\{ - \int _{t}^{\infty }(s-t)u^2(s)ds\right\}, {\tag{3.5}} \end{equation}$$

where $u(x)$ is the Hastings–McLeod solution of the second Painlevé equation, i.e., the solution of the ODE

$$\begin{equation} u_{tt} = tu +2u^3, {\tag{3.6}} \end{equation}$$

uniquely determined by the boundary condition,

$$\begin{equation} u(t) = \frac{1}{2\sqrt {\pi }}t^{-1/4}e^{-\frac{2}{3}t^{3/2}}(1 + o(1)), \quad t \to +\infty . {\tag{3.7}} \end{equation}$$

Similar Painlevé representations for the other two Tracy–Widom distribution functions (see 10) have the form

$$\begin{equation} F_{1}(t) = F_{2}^{1/2}(t)E(t), \quad F_{4}(t) = \frac{1}{2}\left\{E(t) + \frac{1}{E(t)}\right\}F_{2}^{1/2}(t), {\tag{3.8}} \end{equation}$$

where

$$\begin{equation*} E(t) = \exp \left\{ - \frac{1}{2}\int _{t}^{\infty }u(s)ds\right\}, \end{equation*}$$

and $u(x)$ is the same Hastings–McLeod solution to the second Painlevé transcendent.⁠² Painlevé equations are integrable in the sense of Lax pairs, which, in particular, implies that their solutions admit a Riemann–Hilbert (RH) representation. A Riemann–Hilbert representation can be viewed as the nonabelian analogue of the familiar integral representations of the classical special functions, such as the Bessel functions, the Airy function, etc. A key consequence of the RH representation is that the Painlevé functions possess one of the principal features of a classical special function—a mechanism, viz. the nonlinear steepest descent method, to evaluate explicitly relevant asymptotic connection formulae (see, e.g., FIKN). Specifically, in the case of the Hastings–McLeod solution, the integrability of the second Painlevé equation underlies the fact that, in addition to the asymptotic behavior 3.7 at $t = +\infty$, one also knows the asymptotic behavior of $u(t)$ at $t = -\infty$, which is described in HaMc⁠³ as

²

A similar representation for the sine-kernel determinant (involving, this time, a special solution of the fifth Painlevé equation) was given earlier in JMMS.

✖

³

Hastings and McLeod derived 3.9 using the inverse scattering transform which was a precursor of the Riemann–Hilbert method.

✖

$$\begin{equation} u(t) = \sqrt {-\frac{t}{2}} + O\Bigl ((-t)^{-5/2}\Bigr ), \quad t \to -\infty . {\tag{3.9}} \end{equation}$$

Herein lies the importance of the Tracy–Widom formulae 3.5 and 3.8—they provide the key distribution functions of random matrix theory with explicit representations that are amenable to detailed asymptotic analysis.

An immediate corollary of formula 3.5 (and the known asymptotics 3.9 of the Painlevé function) is an explicit formula for the large negative $t$ behavior of the distribution function $F_{2}(t)$:

$$\begin{equation} \ln F_2(t) = \frac{t^3}{12} - \frac{1}{8}\ln (- t )+ c_0 + o(1), \quad t \to -\infty . {\tag{3.10}} \end{equation}$$

The value of the constant $c_0$ was conjectured by Tracy and Widom in the paper 6 to be the same as for the sine-kernel determinant, which is given by

$$\begin{equation} c_0 = \frac{1}{24}\ln 2 + \zeta '(-1), {\tag{3.11}} \end{equation}$$

where $\zeta (s)$ is the Riemann zeta-function.⁠⁴

⁴

This conjecture was independently proved in DIK and BBD. Also, in BBD, similar asymptotic results were established for the other two Tracy–Widom distributions.

✖

Another important advantage of the Tracy–Widom formula 3.8 is that it makes it possible to design, using the exact connection formula 3.7, 3.8, 3.9, a very efficient scheme (see Die, DBT) for the numerical evaluation⁠⁵ of the distribution functions $F_{\beta }$.

⁵

It should be noted that the Hastings–McLeod solution is very unstable; indeed, a small change in the pre-exponential numerical factor in the normalization condition 3.7 yields a completely different type of behavior of $u(x)$ for negative $x$, including the appearance of singularities (e.g., see again FIKN). Herein lies the importance of the knowledge of the behavior 3.9 in order to adjust the numerical procedure appropriately. Said differently, knowledge of the connection formulae makes it possible to transform an unstable ODE initial value problem into a stable ODE boundary value problem.

✖

The appearance of the Painlevé equations in the Tracy–Widom formulae is not accidental. As follows from results in 7, the connection to integrable systems is already encoded in the determinant formula 3.4. Indeed, the Airy-kernel integral operator belongs to a class of integral operators with kernels of the form

$$\begin{equation} K(z,z') = \frac{\phi (z)\psi (z') - \psi (z)\phi (z') }{z-z'}, {\tag{3.12}} \end{equation}$$

acting in $L_{2}(J)$, where $J$ is a union of intervals,

$$\begin{equation*} J = \bigcup _{k=1}^{l}[a_{2k-1}, a_{2k}], \end{equation*}$$

and $\phi (z)$, $\psi (z)$ are $C^{1}$-functions. This class of integral operators has appeared frequently in many applications related to random matrices and statistical mechanics.⁠⁶ In 7, it is shown that if the functions $\phi (z)$ and $\psi (z)$ satisfy a linear differential equation,

⁶

An integral operator with kernel 3.12 is a special case of a so-called integrable Fredholm operator, i.e., an integral operator whose kernel is of the form

$$\begin{equation*} K(z,z') = \frac{\sum _{j=1}^{m}f_j(z)h_{j}(z')}{z-z'}, \quad z, z' \in C, \end{equation*}$$

with some functions $f_{j}(z)$ and $h_{j}(z)$ defined on a contour $C$. This type of integral operator was singled out as a distinguished class in IIKS (see also Dei2; in a different context, unrelated to integrable systems, these operators were also studied in the earlier work Sak). A crucial property of a kernel 3.13 is that the associated resolvent kernel is again an integrable kernel. Moreover, the functions “$f$” and “$h$” corresponding to the resolvent are determined via an auxiliary matrix Riemann–Hilbert problem whose jump-matrix is explicitly constructed in terms of the original $f, h$-functions.

✖

$$\begin{equation*} \frac{d}{dz}\begin{pmatrix} \phi (z)\\ \psi (z)\end{pmatrix} =\Omega (z)\begin{pmatrix} \phi (z)\\ \psi (z)\end{pmatrix}, \end{equation*}$$

where

$$\begin{equation*} \Omega (z) = \frac{1}{m(z)}\begin{pmatrix} A(z)&B(z)\\ C(z)&-A(z)\end{pmatrix}, \end{equation*}$$

and $m(z)$, $A(z)$, $B(z)$, and $C(z)$ are polynomials, then the Fredholm determinant, $\det (1 - K)$, can be expressed in terms of the solution to a certain system of nonlinear partial differential equations with the end points $a_j$ as independent variables. This system is integrable in the sense of Lax,⁠⁷ and in the case of the Airy kernel it reduces to a single ODE—the second Painlevé equation.

⁷

The Lax-integrability of the Tracy–Widom system associated with the kernel 3.12 was proven in Pal, where this system was identified as a special case of isomonodromy deformation equations and the Fredholm determinant $\det (1-K)$ was identified as the corresponding $\tau$-function.

✖

In 8, the results of 7 were extended to the general $m\times m$ case, that is, to the kernels of the form

$$\begin{equation} K(z,z') = \frac{\sum _{j,i=1}^{m}c_{ij}\phi _j(z)\phi _{i}(z')}{z-z'}, {\tag{3.13}} \end{equation}$$

acting in $L_{2}(J)$. Here, $J$ is, as before, a union of intervals, the constant matrix $C \equiv \{c_{ij}\}$ is antisymmetric, and $\phi _j(z)$ are $C^{1}$-functions satisfying the linear differential equation

$$\begin{equation} m(z)\frac{d\phi (z)}{dz} = \Omega (z)\phi (z), \quad \phi (z) \coloneq \begin{pmatrix} \phi _1(z)\\ \phi _2(z)\\ .\\ .\\ .\\ \phi _m(z)\end{pmatrix}, {\tag{3.14}} \end{equation}$$

where $m(z)$ is a scalar polynomial while $\Omega (z)$ is a matrix with polynomial entries connected to the matrix $C$ by

$$\begin{equation*} \Omega ^{T}(z)C + C\Omega (z) =0. \end{equation*}$$

The main result of 8 is the derivation, for such operators $K$, of a system of PDEs, with the $a_j$ as independent variables, whose solution determines the logarithmic derivatives of $\det (1-K)$ with respect to $a_j$.

Figure 2.

Harold and his wife, Linda Larkin, at Lake Sonoma in 2013. Photo courtesy of Craig Tracy.

4. Universality of the Tracy–Widom distributions

A remarkable fact is that the Tracy–Widom distribution functions $F_{\beta }(z)$ appear in a large and growing number of applications in random matrix theory as well as in areas beyond random matrices. In this section we will describe some of these applications.

4.1. General invariant ensembles and the Wigner ensembles

The GUE, GOE, and GSE ensembles are, respectively, particular examples of unitary (UE), orthogonal (OE), and symplectic (SE) ensembles of random matrix theory. The joint probability density of the eigenvalues for these ensembles are given by (cf. 3.1)

$$\begin{equation} P_{N \beta }(\lambda _1, \dots , \lambda _N) = Z^{-1}_{N\beta }\,e^{-\eta _{\beta }\sum _{j=1}^{N}V(\lambda _j)} \prod _{j<k}|\lambda _j - \lambda _k|^{\beta }, {\tag{4.1}} \end{equation}$$

where $\eta _{1}=1$, $\eta _{2} = 1$, $\eta _4 = 2$ (this reflects the fact that the eigenvalues of self-dual Hermitian matrices come in pairs), and the potential $V(\lambda )$ is a polynomial with even leading term,

$$\begin{equation} V(\lambda ) = \sum _{k=1}^{2m}t_k\lambda ^k, \quad t_{2m} > 0. {\tag{4.2}} \end{equation}$$

(More general potential functions are also allowed.) GUE, GOE, and GSE correspond to the choice $V(\lambda ) = \lambda ^2$. Consider again the upper edge fluctuations of the spectrum. Edge-universality states that for a given polynomial potential $V(\lambda )$ there exist scalars $z_N^{(\beta )}$ and $s_{N}^{(\beta )}$ such that

$$\begin{equation} \lim _{N\to \infty }\operatorname {Prob}\left(\frac{\lambda _{\text{max}} -z^{(\beta )}_N}{s^{(N)}_{N}} \leq t\right) = F_{\beta }(t), {\tag{4.3}} \end{equation}$$

where $F_{\beta }(t)$, $\beta = 1, 2, 4$ are exactly the same Tracy–Widom distribution functions as in 3.2. For $\beta = 2$ this result was obtained in DKMVZ1DKMVZ2 and for $\beta = 1, 4$ in DeGi1DeGi2. It should be noted that the authors of DeGi1DeGi2 based their analysis on Widom’s papers 1112 where important relations between symplectic, orthogonal, and unitary ensembles were established for finite $N$. These relations allowed the authors of DeGi1DeGi2 to use, in the symplectic and orthogonal cases, the asymptotic estimates previously obtained in DKMVZ2 for the unitary ensembles. There are also a variety of universality results for eigenvalues in the bulk of the spectrum (see, e.g., Dei3 and references therein) for UE, OE, and SE.

Complementary to the invariant ensembles (UE, OE, and SE) are the so-called Wigner matrix ensembles, i.e., random matrices with independent identically distributed entries. Tracy–Widom edge universality for these ensembles is also valid, and this important fact was proved by Soshnikov Sos. In addition, universality in the bulk of the spectrum for the Wigner ensembles is now established; see, e.g., Erd for a comprehensive survey of the results and a historical review.

The Tracy–Widom distribution appears also in Hermitian matrix models with varying weights. These are unitary ensembles with the potential $V(\lambda )$ in 4.1 replaced by $NV(\lambda )$. Universality results of the form 4.3 are also true for such scaled potentials (see, e.g., DKMVZ1, Dei3 and the references therein). For potentials $NV(\lambda )$ a key issue is the large $N$ asymptotics of the partition functions

$$\begin{equation} Z_N = \int _{-\infty }^{\infty }\cdots \int _{-\infty }^{\infty }e^{-N\sum _{j=1}^{N}V(\lambda _j)} \prod _{j<k}(\lambda _j - \lambda _k)^{2}d\lambda _1\cdots d\lambda _N. {\tag{4.4}} \end{equation}$$

Under certain regularity conditions, the leading large $N$ behavior of the partition function is described by the relation Joh1

$$\begin{equation} Z_N \sim e^{-N^2F_0},\quad N \to \infty , {\tag{4.5}} \end{equation}$$

where $F_0$ is the free energy of the model, and is given explicitly in terms of the equilibrium measure corresponding to the potential $V(\lambda )$. The equilibrium measure is the unique probability measure $\nu _{\text{eq}}$ which minimizes the energy functional

$$\begin{equation*} I_{V}(\nu ) = -\int \int _{\mathbb{R}^2}\ln |\lambda -\lambda '|d\nu (\lambda )d\nu (\lambda ') + \int _{\mathbb{R}}V(\lambda )d\lambda . \end{equation*}$$

The formula for the free energy reads

$$\begin{equation*} F_0 = I_{V}(\nu _{\text{eq}}). \end{equation*}$$

As shown in DKM, in the case that $V(\lambda )$ is real analytic and $\frac{V(\lambda )}{\ln |\lambda |} \to + \infty$ as $|\lambda | \to \infty$, the equilibrium measure is absolutely continuous and is supported on a finite number of intervals. The derivative of $\nu _{\text{eq}}$ coincides with the mean limiting density of eigenvalues of the associated UE. If the support of the equilibrium measure consists of one interval, then one can prove (see EM and also BI) the existence of a full asymptotic expansion,

$$\begin{equation} \ln \frac{Z_{N}}{Z_{N}^{\text{Gauss}}} \sim -N^2e_0 + \sum _{g=1}^{\infty }\frac{e_{g}}{N^{2g-2}}, \quad N \to \infty , {\tag{4.6}} \end{equation}$$

where $Z_{N}^{\text{Gauss}}$ stands for the partition function of GUE ( $V(\lambda ) = \lambda ^2$). Since 1980 the expansion 4.6 has played a prominent role in the application of random matrices to enumerative topology. This role was first recognized in BIZ together with the discovery of its deep connections to the counting of graphs on Riemann surfaces (see EM for more details).

The number of intervals in the support of the equilibrium measure depends on the values of the parameters $t_1, \dots , t_{2m}$ in 4.2. It is of great interest to study the transition regimes; i.e., regimes involving values of $t_j$ near critical values when the number of the intervals in the support of the equilibrium measure changes. The first example where one can observe this important critical phenomenon is the even quartic potential $V(\lambda ) = t_4\lambda ^4 + t_2\lambda ^2$. Up to a trivial renormalization, the potential $V$ can be chosen in the form

$$\begin{equation} V(\lambda ) = \frac{1}{4t^2}\lambda ^4 + \left(1-\frac{2}{t}\right)\lambda ^2. {\tag{4.7}} \end{equation}$$

For $t>1$ the support of the equilibrium measure corresponding to this potential consists of one interval. Moreover, as shown in BI, all the coefficients $e_{g}(t)$ of the asymptotic series 4.6, as functions of $t$, are real analytic on $[1, \infty )$. If $t <1$, the support becomes two intervals. At the critical value, $t = 1$, the support consists of one interval but the density function has a double zero inside the support at $\lambda =0$. The relevant double scaling limit near the critical point $t =1$ is prescribed by the scaling condition

$$\begin{equation} \Bigl |(t-1)N^{2/3}\Bigr | < C. {\tag{4.8}} \end{equation}$$

As $N \to \infty$ and $t \to 1$ satisfying 4.8, the partition function admits the following asymptotic representation BI: