We review the authors’ recent work where we obtain the uniform large s asymptotics for the Fredholm determinant , 0 ≤ γ ≤ 1. The operator Ks acts with kernel Ks(x, y) = sin(s(x − y))/(π(x − y)), and D(s, γ) appears for instance in Dyson’s model of a Coulomb log-gas with varying external potential or in the bulk scaling analysis of the thinned Gaussian unitary ensemble.
Dedicated to the memory of Ludwig Faddeev
I. INTRODUCTION
Consider the determinantal point process {xj} on the real line with a correlation kernel
It is well known that this process appears in the bulk scaling limit for random matrices (in particular, in the Gaussian Unitary Ensemble, GUE) where the average distance between particles xj is unity. See Ref. 19 and the references therein to Lenard’s work, that for any whose support is inside a bounded measurable set B,
where the determinant on the right-hand side of (1.1) is the Fredholm determinant of the trace-class integral operator ϕK on L2 (B) with kernel ϕ(x)K(x, y). Provided we define
then applying (1.1) to ϕ(x) = 1 − e−W(x), we find
where
We can interpret > 0 above as an external potential applied to the Coulomb log-gas particle system {xj} whose effect is to push particles out of the interval . We are interested in D(s, γ) for large s and 0 ≤ ≤ +∞. This has been a long standing problem, and we shall first review relevant classical results.
II. ASYMPTOTIC RESULTS
A. Fixed and slowly growing
If = +∞, i.e., γ = 1, then D(s, 1) is the probability that there are no sine process particles in the interval . In this case,
where and ζ′(z) is the derivative of Riemann’s zeta function. The main term, , in (2.1) was first conjectured by Dyson, followed by a conjecture for by des Cloizeaux and Mehta in Ref. 12. A full asymptotic expansion for D(s, 1), including the numerical value of c0, was then identified by Dyson15 in 1976 who used inverse scattering techniques for Schrödinger operators and an earlier study of Widom21 on Toeplitz determinants. Dyson’s calculations were not fully rigorous, and the first proof of the main term in (2.1) was given by Widom22 in 1994. A proof of the full expansion (except for the value of c0) was carried out in Ref. 14, and c0, the so-called Widom-Dyson constant, was finally proved in three different ways in Refs. 13, 17, and 18.
If (which includes the case of fixed finite v), there exist constants s0, cj > 0 such that for s > s0,
where G(z) is the Barnes G-function
with Euler’s constant γE, and
Expansion (2.2) in the case of fixed finite as s → +∞ and with r(s, ) replaced by was first established by Basor and Widom in 19831 and later, independently, by Budylin and Buslaev.10 It was extended to the presented range with varying in Ref. 7. At this point, it is worthwhile to contrast the Gaussian decay in the leading order of (2.1) with the leading exponential decay (for fixed ) in (2.2). The underlying non-trivial transition (as worked out in the authors’ studies, Refs. 6–8) is summarized in Subsection II B.
B. Faster growing
We shall now describe the asymptotic transition between (2.1) and (2.2) which takes place when is allowed to grow faster with s than . In more detail, we describe the large s asymptotics of D(s, γ) in the quarter-plane (s, ), s, > 0; see Fig. 1. Heuristically, this problem was first addressed by Dyson in Ref. 16 who discovered an oscillatory behavior in the asymptotics and the appearance of Jacobi-theta functions (see Ref. 6 for a discussion of Dyson’s calculations). Let
so that γ = 1 − = 1 − e−2κs. First, we have
The quarter plane (s, v) with s, v, q sufficiently large. We indicate the major diagonal v = s as a solid black line and a few dashed-dotted Stokes curves . The validity of Theorems 2.1 (with q = 3) and 2.2 is roughly indicated with two circular arcs. The dotted green curve is .
The quarter plane (s, v) with s, v, q sufficiently large. We indicate the major diagonal v = s as a solid black line and a few dashed-dotted Stokes curves . The validity of Theorems 2.1 (with q = 3) and 2.2 is roughly indicated with two circular arcs. The dotted green curve is .
Observe that can be written as which implies that = su. Thus, if , only the factors indexed with j = 0, 1, …, q − 1 contribute to the product in (2.4), while terms with larger indices j yield a decaying with s contribution. In geometric terms, if one moves to the right in Fig. 1, each time a Stokes curve is crossed, a factor is added to (2.4). One can roughly interpret this phenomenon as a particle jumping to the center of the interval (−1, 1), i.e., the rescaled original interval . Note that for , estimate (2.3) shows that the leading order asymptotics of D(s, γ) are the same as for the case = ∞. The asymptotic expansions (2.3) and (2.4) in the case of fixed q (rather than growing with s, as stated in Theorem 2.1) were proven in Ref. 6 using general operator theoretical arguments and a result of Slepian20 on the asymptotics of the eigenvalues of the sine-kernel integral operator. Each factor in the product (2.4) is related to an eigenvalue of the associated integral operator. In full generality, Theorem 2.1 is proven in Ref. 8 using Riemann-Hilbert nonlinear steepest descent techniques.
Theorems 2.1 and 2.2 provide a full analytic description of the transition between (2.1) and (2.2). In terms of the initial log-gas particle system, the following rough picture for the particles’ behaviors follows from our analysis: First, for , no particles are expected in the interval (−1, 1). As the rate of increase of with s slows down, particles begin jumping one by one to the center of the interval: particles at the interval center correspond to the asymptotics (2.4), and their behavior is described by classical Hermite polynomials; see Ref. 8. As q increases, a transition takes place to the theta-function behavior of (2.5), at first with a close to zero. The particles then tend to fill in the centermost subinterval (−a, a) ⊂ (−1, 1), where the soft edges ±a are characterized by an approximate Airy function behavior; see Ref. 6. As the growth of with s slows down further, the subinterval (−a, a) increases and eventually a approaches 1. At this point, a transition occurs to the limit s → ∞ with fixed.
III. THINNED SINE PROCESS
In closing, we recall the following additional interpretation of D(s, γ), noticed by Bohigas and Pato:3 (see, the discussion in Ref. 6) The Fredholm determinant D(s, γ) equals the probability that the interval contains no particles of the thinned process obtained from {xj} by removing each particle xj independently with likelihood 1 − γ. Indeed, using the well-known Fredholm determinant expression for the probability of n sine-kernel process particles xj being in a given interval, we obtain via inclusion-exclusion principle
In short, D(s, γ) is simply the gap-probability of in the thinned sine-process and thus also the central building block in the evaluation of the thinned Gaudin distribution, i.e., the limiting distribution of gaps between consecutive thinned bulk eigenvalues. Following the arguments of Bohigas and Pato, we contract the s-scale in order to keep to unity the average distance between the thinned particles; i.e., instead of D(s, γ), we consider
It follows from (2.2) that at leading order, as γ = 1 − ↓ 0,
which is the gap probability of in the particle system obtained from a Poisson point process by removing each particle with probability 1 − γ [in order to obtain (3.1), Bohigas and Pato used the asymptotics for fixed , but a rigorous analysis requires the result for varying γ, hence varying ; see Ref. 6]. Indeed, the gap probability in this thinned Poissonian system is exactly given by
On the other hand, for γ = 1, we obtain , i.e., by thinning bulk eigenvalues in the GUE, we can interpolate between a Poissonian and a random matrix theory system. For other matrix ensembles, such as the Circular Unitary Ensemble, CUE, corresponding gap probability phenomena have been studied in Ref. 11. For interpolation mechanisms between random matrix theory systems and Weibull modelled particle systems, we refer the interested reader to Refs. 2, 5, 4, and 9.
ACKNOWLEDGMENTS
T.B. acknowledges support of the AMS and the Simons Foundation through a travel grant. The work of P.D. is supported under NSF Grant No. DMS-1300965, and A.I. acknowledges support of the NSF Grant No. DMS-1700261 and the Russian Science Foundation Grant No. 17-11-01126. I.K.’s work is supported by the Leverhulme Trust research project Grant No. RPG-2018-260.