A convergence framework for Airy_β line ensemble
via pole evolution

Jiaoyang Huang Department of Statistics and Data Science, University of Pennsylvania, Pennsylvania, PA, USA [email protected] and Lingfu Zhang Division of Physics, Mathematics and Astronomy, California Institute of Technology, Pasadena, CA, USA. and Department of Statistics, University of California, Berkeley, CA, USA. [email protected]

Abstract.

The Airy_β line ensemble is an infinite sequence of random curves. It is a natural extension of the Tracy-Widom_β distributions, and is expected to be the universal edge scaling limit of a range of models in random matrix theory and statistical mechanics. In this work, we provide a framework of proving convergence to the Airy_β line ensemble, via a characterization through the pole evolution of meromorphic functions satisfying certain stochastic differential equations. Our framework is then applied to prove the universality of the Airy_β line ensemble as the edge limit of various continuous time processes, including Dyson Brownian motions with general $\beta$ and potentials, Laguerre processes and Jacobi processes.

1. Introduction

During the 18th century, De Moivre established the Gaussian distribution for sums of independent binomial variables, a concept later generalized by Laplace. Gauss popularized the central limit theorem, used for error evaluation in systems characterized by independence. In recent decades, there has been increasing interest in understanding fluctuations in highly correlated systems, leading to the emergence of a different family of distributions known as Tracy-Widom_β, which are indexed by a positive parameter $\beta$ . Historically, such Tracy-Widom_β distributions for $\beta=1,2,4$ were first observed in Random Matrix Theory, as the scaling limit of the extreme eigenvalues of the classical matrix ensembles [121, 122, 58, 82, 84, 86]. Later, such extreme eigenvalue statistics are proven to be universal, in the sense that Tracy-Widom_β limits hold for a wide range of random matrix models, including adjacency matrices of random graphs, which are usually sparse. See e.g., [118, 120, 98, 54, 90, 78, 71, 72, 79, 16]. Beyond matrix models, Tracy-Widom_β distributions also appear in lots of different random systems, such as random tilings, 1+1 dimensional random growth models, exclusion processes, planar random geometry such as first/last passage percolation models, and square ice models (six-vertex models). See e.g., [11, 82, 8, 111, 84, 83, 104, 5, 12, 115, 21, 30, 127, 70, 80, 13, 14, 10, 28, 4, 126]. Many of these models are related to each other through their underlying integrable structures, and are in the so called Kardar-Parisi-Zhang (KPZ) universality class, a topic in probability theory that has been intensively studied in recent years.

A motivation of the current paper is to better understand the mathematical structures behind Tracy-Widom_β distributions, and to develop new methods of proving convergence to them. The main object studied here is the Airy_β line ensemble (ALE_β) $\{{\mathcal{A}}^{\beta}_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ , which can be defined as a random process on $\mathbb{N}\times{\mathbb{R}}$ or a family of continuous random processes on ${\mathbb{R}}$ , with $\beta>0$ being a parameter, and is ordered, i.e., ${\mathcal{A}}^{\beta}_{1}(t)\geq{\mathcal{A}}^{\beta}_{2}(t)\geq{\mathcal{A}}^% {\beta}_{3}(t)\geq\cdots$ for any $t$ . They are natural generalizations of the Tracy-Widom_β distributions, akin to Brownian motions being generaliztions of Gaussian distributions, and have Tracy-Widom_β as the one-point distribution of the top line ${\mathcal{A}}^{\beta}_{1}$ . These ALE_β are believed to be universal objects, in the sense of being the scaling limit of many random matrix models and interacting particle systems. However, basic properties of these processes as well as these convergences remain quite mysterious so far, except for the special setting of $\beta=2$ , where a determinantal structure is present and has been largely exploited using algebraic methods.

In this paper, we take a new perspective, and our main result is a characterization of ALE_β in terms of its Stieltjes transform and a system of stochastic differential equations (SDE). Our result provides a new framework to prove convergence. As some examples, we prove convergence to ALE_β from various random processes, including the classical $\beta$ -Hermite/Laguerre/Jacobi processes and their generalizations. We note that some of these were previously unknown even in the $\beta=1,4$ cases, where such convergences can be interpreted as the joint convergence of extreme eigenvalues of correlated real/quaternion random matrices. Beyond these, our framework should also be applicable to prove convergence to ALE_β for many other models; and for some of them, even the Tracy-Widom_β limits were previously unknown. Moreover, our characterization also reveals some useful information for ALE_β, such as Hölder properties and collision of adjacent lines.

1.1. Background

We next provide some setup, starting with some more classical processes.

1.1.1. Edge limit of general $\beta$ -ensembles

Tracy-Widom_β distributions for general $\beta>0$ was introduced and studied in [53, 119] by Edelman and Sutton, and [113] by Ramirez, Rider and Virág, as the scaling limit of the extreme eigenvalue of Gaussian $\beta$ -ensembles. More generally, $\beta$ -ensemble is a probability distribution on $n$ particle system $x_{1}\geq x_{2}\geq\cdots\geq x_{n}$ , with probability density:

(1.1)

\frac{1}{Z_{n,\beta,W}}\prod_{i<j}|x_{i}-x_{j}|^{\beta}\prod_{i=1}^{n}W(x_{i}),

where $Z_{n,\beta,W}$ is a renormalization constant, and $W\geq 0$ is the weight function. There are three special cases referred to as the classical ensembles, which are defined by

(1.5)

\displaystyle W(x)=\left\{\begin{array}[]{cc}e^{-\beta x^{2}/4},&\text{Hermite% /Gaussian ensemble}\\ x^{\beta(m-n+1)/2-1}e^{-\beta x/2},&\text{Laguerre ensemble}\\ x^{\beta(p-n+1)/2-1}(1-x)^{\beta(q-n+1)/2-1},&\text{Jacobi ensemble}\\ \end{array}\right.

These classical ensembles for $\beta=1,2,4$ originated from the study of eigenvalue distributions of random matrices. They represent the joint distributions of the eigenvalues of size $n$ Gaussian, Wishart, and Jacobi matrices. These matrices and their extreme eigenvalues, with $\beta=1$ corresponding to the real case, have been widely used in high-dimensional statistical inference (see the survey by Johnstone [87]). Beyond random matrix theory, $\beta$ -ensembles also describe the one-dimensional Coulomb gas in physics [58], and are connected to orthogonal polynomial systems [91].

As mentioned earlier, as $n\to\infty$ , the distribution of the largest eigenvalues converges to the Tracy-Widom_β distribution for $\beta=1,2,4$ , respectively. More generally, one can consider the edge limit, i.e., the joint scaling limit of the top $k$ eigenvalues for any arbitrary $k$ , as a point process. It has been shown [41, 42, 43] that this edge limit does not depend on the potential function $W$ , but it varies for each of $\beta=1,2,4$ . The $\beta=2$ edge limit is also known as the Airy point process.

For $\beta$ other than $1,2,4$ , obtaining such edge limit was a challenging problem, partly due to the relative lack of exact-solvable structures. Based on a tri-diagonal random matrix model discovered by Dumitriu and Edelman [49], this was resolved in [113], where the edge limit of Gaussian and Laguerre $\beta$ -ensemble is shown to be the eigenvalues of the $\beta$ stochastic Airy operator (SAO_β), which is also called the Airy_β point process. In particular, for each $\beta>0$ , the law of the largest eigenvalue of SAO_β, i.e., the top particle in the Airy_β point process, is then defined as the Tracy-Widom_β distribution. Later, such Airy_β point process limit is also extended to more general $W$ [95, 25, 18, 17]. Analogous edge limits for discrete $\beta$ -ensembles have been derived in [67].

1.1.2. Airy line ensemble

In another direction, the Tracy-Widom₂ distribution and the ( $\beta=2$ ) Airy point process are extended to the Airy line ensemble (ALE), an ordered family of random processes $\{{\mathcal{A}}_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ , where each ${\mathcal{A}}_{i}$ is continuous, and they are jointly stationary in the ${\mathbb{R}}$ direction (see Figure 1). ALE was introduced by Prahöfer and Spohn [109], as the scaling limit for the multi-layer polynuclear growth (PNG) model from the KPZ universality class. The top line ${\mathcal{A}}_{1}$ is known as the stationary Airy₂ process, whose one-point marginal is the Tracy-Widom₂ distribution; and for any $t\in{\mathbb{R}}$ , $\{{\mathcal{A}}_{i}(t)\}_{i\in\mathbb{N}}$ is the Airy point process. ALE plays a central role in KPZ, in particular through the construction of the directed landscape [37]. (See also [38] for computing passage times in the directed landscape from ALE.)

ALE is particularly useful in KPZ, partly due to its Brownian Gibbs property, which was recognized by Corwin and Hammond [35]. Specifically, for ALE minus a parabola, it inside any domain, conditional on the boundary, is given by non-intersecting Brownian bridges. This fact is later widely used as a powerful tool to study ALE and many KPZ class models. Aggarwal and the first-named author provided a strong characterization of ALE, demonstrating that ALE (minus a parabola) is the only random process on $\mathbb{N}\times{\mathbb{R}}$ with the Brownian Gibbs property as well as approximating a parabola [6]. Such a strong characterization would be a powerful tool to prove convergence to ALE and establishing KPZ universality for various models; see e.g. [4].

Figure 1. An illustration of ALE

1.2. Airy_β line ensemble

From the success of ALE, the next question is to construct a time dependent evolution for Tracy-Widom_β (and more generally Airy_β point process) for any $\beta>0$ . Following [64] where this is formally introduced, we call it the Airy_β line ensemble (ALE_β). There are several problems in this program:

•

Construction What should it be? How to construct it?
•

Description What are its properties? Ideally, can some precise information be given?
•

Universality Why is ALE_β natural and interesting? Can it be shown to be the universal scaling limit of many natural random processes, as the $\beta=2$ case?

Towards these questions, there have been many results focusing on different aspects of ALE_β (some tracing back to the studies of Airy_β point process, or for special $\beta$ ): infinite-dimensional SDE [97, 106, 107, 88], correlation function [104, 11, 85, 109], and Laplace transform [118, 102, 103, 117, 81, 64]. In this paper, we provide a different perspective using Stieltjes transform, tailored to the universality problem. Moreover, we view results presented here and in the concurrent paper [64] (by Gorin, Xu, and the second-named author, to be discussed shortly) complement to each other.

We now give a more detailed account on the state of the art, and further motivate our results.

1.2.1. The edge limit construction and convergence

Beyond ALE where $\beta=2$ , ALE_β was also constructed and studied for several other special values: in [117] for $\beta=1$ (and the arguments there should also go through for $\beta=4$ ), and in [60] for $\beta=\infty$ . For general $\beta>0$ , even its construction is relatively recent. One potential way, as inspired by the fact (from [35]) that ALE is the edge scaling limit of the $\beta=2$ Dyson Brownian motion (DBM), is to consider general $\beta$ DBM:

(1.6)

{\rm d}\lambda_{i}(t)=\sqrt{\frac{2}{\beta}}{\rm d}B_{i}(t)+\sum_{\begin{% subarray}{c}1\leq j\leq n\\ j\neq i\end{subarray}}\frac{{\rm d}t}{\lambda_{i}(t)-\lambda_{j}(t)}-\frac{1}{% 2}\lambda_{i}(t){\rm d}t,

where $n\in\mathbb{N}$ and $\{B_{i}\}_{i=1}^{n}$ are independent two-sided Brownian motions. There is a solution to this SDE, such that $\{\lambda_{i}(t)\}_{i=1}^{n}$ for any $t$ is a Gaussian $\beta$ -ensemble; and this is known as the (stationary) DBM of size $n$ with parameter $\beta$ . (See e.g., [9] for some more backgrounds on DBM.) One can then define ALE_β $\{{\mathcal{A}}^{\beta}_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ as its edge limit, i.e., the limit of $(i,t)\mapsto n^{1/6}\lambda_{i}(tn^{-1/3})-2n^{2/3}$ as $n\to\infty$ . Note that for any fixed $t$ , $\{{\mathcal{A}}^{\beta}_{i}(t)\}_{i\in\mathbb{N}}$ should be the Airy_β point process, and ${\mathcal{A}}^{\beta}_{1}(t)$ should follow the Tracy-Widom_β distribution.

Beyond DBM, another potential way of constructing ALE_β is via the edge limit of Gaussian corners processes, which are random Gelfand-Tsetlin patterns, and can be viewed as multi-level generalizations of Gaussian $\beta$ ensembles. For $\beta=1,2,4$ , they can be defined as eigenvalues of the top-left corners of different sizes, for Wigner matrices with Gaussian real, or complex, or quaternionic entries (see e.g. [59, 15, 101]).

To justify the above definition of ALE_β, the above stated convergence should be proved. In [97], this is achieved for stationary DBM and any $\beta\geq 1$ via a coupling argument. In [64] the convergences for both stationary DBM and Gaussian corners process are established, for any $\beta>0$ . The proofs are via explicit formulas, also showing that both limits are the same. Very recently, in [52] a tridiagonal model for DBM is proposed, which may also be used to demonstrate convergence to ALE_β. Additionally, the arguments in this paper provide an alternative proof of DBM convergence.

1.2.2. Description via explicit expressions

Historically, the theory of Tracy-Widom_β distributions used to largely rely on explicit formulas, based on determinantal/Pfaffian structures or matrix models (see e.g., [121, 123, 116, 118, 113, 63]). As for ALE_β, formulas used to be available only for $\beta=2$ , i.e., ALE, in the construction by Prahöfer and Spohn [109]; and for $\beta=\infty$ in [60]. There are various challenges in obtaining precise formulas for general $\beta$ , primarily due to the lack of structures in this generality. As above mentioned, the first construction of ALE_β for $\beta\geq 1$ in [97] was via the edge limit of DBM, using a more abstract convergence argument. Then in [64], precise formulas for the Laplace transform of ALE_β are obtained. Such Laplace transforms also determine the law of ALE_β, thereby [64] gives a direct and explicit definition of ALE_β for any $\beta>0$ .

1.2.3. Uniqueness and universality

As indicated by the convergence of both DBM and Gaussian corners process to the same limit, i.e., ALE_β, in [64], it is natural to expect that ALE_β is also the scaling limit of many other well-known processes. Some examples include DBM with general potentials, non-intersecting random walks [61, 93, 74], various other $\beta$ –corners processes [62, 22, 65] and measures on Gelfand-Tsetlin patterns [29, 108], Macdonald processes [20], and $(q,\kappa)$ -distributions on lozenge tilings [23, 46].

A robust and general approach to establishing convergence involves a suitable characterization of ALE_β. Specifically, this means identifying easily verifiable properties of ALE_β and demonstrating that these properties uniquely determine ALE_β. To prove convergence, it would then suffice to establish tightness and confirm that any subsequential limit satisfies these properties.

For ALE_β with $\beta=2$ , an elegant characterization is the Brownian Gibbs property [6] as mentioned above. However, this does not hold for any $\beta\neq 2$ . The next natural candidate of characterization would be an ‘infinite dimensional DBM’, by taking $n\to\infty$ in (1.6). For example, [97] shows that ALE_β (for $\beta\geq 1$ ), i.e., the edge limit of finite dimension DBM, is a solution to such an infinite dimensional DBM in a weak sense.

However, there are several technical challenges in developing a characterization and convergence framework along this direction. First, particles (i.e., those $\lambda_{i}$ in (1.6)) may collide or adjacent particles may get too close, leading to singularities in the drift $1/(\lambda_{i}(t)-\lambda_{j}(t))$ term. In [97], for $\beta\geq 1$ such singularities are ruled out using existing level replusion estimates (see [97, Theorem 2.2] and [25, Theorem 3.2]), which are known only for stationary DBM whose laws are given by $\beta$ -ensembles. Deriving such estimates for other models could be difficult. Moreover, for $\beta<1$ collisions are inevitable. Second, the long-range interactions introduce additional complications when analyzing infinitely many particles. In fact, even the well-posedness of the infinite-dimensional DBM starting from a fixed intial condition appears to be absent from the literature, except for the cases $\beta=1,2,4$ [106, 107, 88] where the underlying algebraic structure is used¹¹1Note that for an analogous infinite-dimensional SDE corresponding to the bulk limit of DBM, such well-posedness has been achieved for $\beta\geq 1$ (see [89, 105, 106, 124]). A key property used in the bulk setting is the cancellation of particle interactions from left and right, and that is absent at the edge.. As a result, for the infinite-dimensional DBM, both proving the uniqueness of solution and verifying it for any subsequential limit face various barriers.

To overcome these difficulties, in this paper we take an alternative approach, and characterize ALE_β as the pole dynamics of meromorphic functions, satisfying a funciton-valued SDE. In other words, we characterize ALE_β via the dynamics of its Stieltjes transform. This method completely circumvent the issue of long-range interactions and collisions, and is applicable for any $\beta>0$ .

1.3. Main characterization result

In the rest of this paper, we fix $\beta>0$ . We study (infinite) line ensembles, by which we mean ordered families of continuous random processes, denoted by $\{{\bm{x}}(t)\}_{t\in I}=\{x_{i}(t)\}_{i\in\mathbb{N},t\in I}$ , for $I={\mathbb{R}}$ or any interval, satisfying $x_{1}(t)\geq x_{2}(t)\geq x_{3}(t)\geq\cdots$ .

As already alluded to, we shall characterize ALE_β using its Stieltjes transform. If $\{{\bm{x}}(t)\}_{t\in{\mathbb{R}}}$ were ALE_β, it would be imperative that (for each $t\in{\mathbb{R}}$ ) the Stieltjes transform of $\{x_{i}(t)\}_{i\in\mathbb{N}}$ (with complex variable $w$ ) asymptotically behaves like $\sqrt{w}$ ²²2Here and throughout this paper, $\sqrt{w}$ , or any rational power of $w\in{\mathbb{C}}$ , is taken to be the branch on ${\mathbb{C}}\setminus{\mathbb{R}}_{-}$ as $w\rightarrow\infty$ in the complex plane. This is because we expect $\{x_{i}(t)\}_{i\in\mathbb{N}}$ , as a one time slice of ALE_β, to be the Airy_β point process; hence, the particles should be close to the zeros of the Airy function $\operatorname{Ai}(w)$ . Therefore, the Stieltjes transform of $\{x_{i}(t)\}_{i\in\mathbb{N}}$ should exhibit similar asymptotic behaviors as $-\operatorname{Ai}^{\prime}(w)/\operatorname{Ai}(w)$ , which is known to behave like $\sqrt{w}$ as $w\rightarrow\infty$ .

We next give a more precise formulation of such asymptotic behaviors. Below we use $\mathbb{H}$ to denote the open upper half complex plane. For Stieltjes transforms, we shall use the notion of Nevanlinna functions, i.e., functions from $\mathbb{H}$ to $\mathbb{H}\cup{\mathbb{R}}$ that are holomorphic.

Definition 1.1.

A measure $\mu$ on ${\mathbb{R}}$ is particle-generated, if it is locally finite, and is in the form of $\sum_{x\in P}\delta_{x}$ , where $P$ (the particles) is an at most countable multiset of real numbers. A Nevanlinna function $Y$ has a Nevanlinna representation of the form

(1.7)

\displaystyle Y(w)=b+cw+\int_{\mathbb{R}}\left(\frac{1}{x-w}-\frac{x}{1+x^{2}}% \right){\rm d}\mu(x),

where $b,c\in{\mathbb{R}}$ , $c\geq 0$ . We say that $Y$ is particle-generated, if the measure $\mu$ in its Nevanlinna representation is particle-generated, namely $\mu=\sum_{x\in P}\delta_{x}$ .

We note that such $Y$ can be extended to a meromorphic function on ${\mathbb{C}}$ with $Y(\overline{w})=\overline{Y(w)}$ , with $P$ being all the poles, and each residue equals $1$ ³³3More precisely, for each $x\in{\mathbb{R}}$ , the residue at $x$ equals the multiplicity of $x$ in $P$ ..

Definition 1.2.

For ${\mathfrak{d}},C_{*}>0$ , a Nevanlinna function $Y$ is $({\mathfrak{d}},C_{*})$ -Airy-like, or simply Airy-like, if

(1)

it is particle-generated, and all the poles are $\leq C_{*}$ ;

(2)

for all $w$ with ${\rm Im}[w]\geq C_{*}\sqrt{{\rm Re}[w]\vee 0+1}$ ,

(1.8)

\displaystyle|Y(w)-\sqrt{w}|\leq\frac{C_{*}{\rm Im}[\sqrt{w}]^{1-{\mathfrak{d}% }}}{{\rm Im}[w]}.

Remark 1.3.

The condition (2) implies the existence of infinitely many poles $x_{1}\geq x_{2}\geq\cdots$ . As we will see shortly (2.3), bounds similar to (1) and (2) (but may with a different domain for $w$ ) would imply that the density of these poles would be close to $\sqrt{-x}$ in ${\mathbb{R}}_{-}$ . Such density closeness can be phrased as quantitative bounds on the distances between the poles and zeros of the Airy function ${\mathfrak{a}}_{1}>{\mathfrak{a}}_{2}>\cdots$ , as stated in 2.3, and will be frequently used in our proofs. Therefore, as a slight misuse of notions, we will refer to such density closeness as Airy-zero approximation.

As another comment on (2): again note that the domain of $w$ is different from that in 2.3, since the domain here is taken to be easily verifiable for the sub-sequential limit of various models, as will be evident in Section 7. In fact, bounds on $|Y(w)-\sqrt{w}|$ for $w$ closer to ${\mathbb{R}}_{+}$ can be readily deduced for Airy-like Nevanlinna functions (see 2.5).

Take a family of random Nevanlinna functions $\{Y_{t}\}_{t\in{\mathbb{R}}}$ . We next state two assumptions.

Assumption 1.4.

For any $t\in{\mathbb{R}}$ , $Y_{t}$ is particle-generated. Moreover, there exists a (deterministic) ${\mathfrak{d}}>0$ , a sequence $t_{1},t_{2},\cdots\to-\infty$ , and a tight family of random numbers $\{C_{*,j}\}_{j\in\mathbb{N}}$ , so that each $Y_{t_{j}}$ is $({\mathfrak{d}},C_{*,j})$ -Airy-like.

Assumption 1.5.

Such $Y_{t}$ is continuous in $t$ , and satisfies the following SDE:

(1.9)

\displaystyle{\rm d}Y_{t}(w)={\rm d}M_{t}(w)+\left(\frac{2-\beta}{2\beta}% \operatorname{\partial}^{2}_{w}Y_{t}(w)+\frac{1}{2}\operatorname{\partial}_{w}% Y_{t}(w)^{2}-\frac{1}{2}\right){\rm d}t,

where $M_{t}(w)$ are complex valued Martingales, with quadratic variation given by

(1.10)

\displaystyle\frac{{\rm d}}{{\rm d}t}\langle M_{t}(w)\rangle=\frac{1}{3\beta}% \operatorname{\partial}_{w}^{3}Y_{t}(w),

and

(1.11)

\displaystyle\frac{{\rm d}}{{\rm d}t}\langle M_{t}(w),M_{t}(w^{\prime})\rangle% =\frac{2}{\beta}\operatorname{\partial}_{w}\operatorname{\partial}_{w^{\prime}% }\frac{Y_{t}(w)-Y_{t}(w^{\prime})}{w-w^{\prime}},

for $w\neq w^{\prime}$ .

We now explain where this SDE comes from. Take the $n$ dimensional stationary DBM $\{\lambda_{i}\}_{i=1}^{n}$ with parameter $\beta$ , i.e., the stationary solution to (1.6). Let $m_{t}$ be the Stieltjes transform of $\{\lambda_{i}(t)\}_{i=1}^{n}$ , i.e.,

m_{t}(z)=\sum_{i=1}^{n}\frac{1}{\lambda_{i}(t)-z},\quad z\in{\mathbb{C}}.

Then one can use Ito’s formula to write out an SDE satisfied by $m_{t}$ ; by taking an appropriate scaling limit from there, one gets the SDE (1.9). More details on this derivation can be found in Section 7.

Our main result states that these two assumptions are sufficient to determine ALE_β uniquely.

Theorem 1.6.

For any $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.4 and 1.5, its poles give a line ensemble, which has the same law as ALE_β.

Several remarks are in line.

(i) Essentiality of 1.4 (in characterizing ALE_β). Both (1) and (2) in Definition 1.2 are necessary: without (1), the line ensemble may be ALE_β plus some additional lines (see [2, 45] for an example in the $\beta=2$ setting); while (2) rules out the possibility that the line ensemble is ALE_β shifted by a (deterministic or random) constant. As already mentioned in Remark 1.3, we’ve aimed to make 1.4 as minimal as possible to ensure broad applicability of our convergence framework. As will be seen in Section 7, Assumption 1.4 is straightforward to verify in these examples.

(ii) DBM convergence. Our proof of 1.6 does not a priori assume the convergence at the edge of stationary DBM. Instead, in Section 7, we show the tightness at the edge of stationary DBM, and that any subsequential limit satisfies 1.4 and 1.5. Therefore, we essentially provide an alternative construction of ALE_β.

(iii) Stationarity. We also note that in 1.6, we do not assume that $Y_{t}$ is stationary. Rather, it is a consequence of the theorem that the poles of $Y_{t}$ converge to ALE_β, which is stationary, and hence $Y_{t}$ is stationary as well. Our proof of 1.6 in fact establishes a natural relaxation for the SDE (1.9). Specifically, for a family of random particle-generated Nevanlinna functions $\{Y_{t}\}_{t\geq 0}$ , if there is a (deterministic) ${\mathfrak{d}}>0$ and a random number $C$ such that $Y_{0}$ is $({\mathfrak{d}},C)$ -Airy-like, and $\{Y_{t}\}_{t\geq 0}$ satisfies 1.5, then for $T\rightarrow\infty$ , the poles of $\{Y_{t}\}_{t\geq T}$ converge to ALE_β, under the uniform in compact topology.

(iv) Stieltjes transform and poles dynamics. Stieltjes transforms and Nevanlinna functions have been widely used to investigate and characterize eigenvalue distributions of random matrix ensembles. See [57] for studies on eigenvalue rigidity, [56, 7] for bulk limits, and [98, 90, 31] for edge limits.

The concept of characterizing the evolution of interacting particle systems through the pole dynamics of meromorphic functions has been explored previously. In integrable systems, it has been demonstrated that the movement of poles in certain solutions of various nonlinear PDEs can be formally linked to the dynamics of particle systems interacting through simple two-body potentials. This discovery was initially made in [34] for equations such as the Korteweg-de Vries and Burgers-Hopf equations, and in [99] for specific integrable Hamiltonian systems. Subsequently, these observations were extended to include elliptic solutions of equations such as the Kadomtsev-Petviashvili equation [94], the Korteweg–de Vries equation [40], the Kadomtsev-Petviashvili hierarchy [125], and the Toda lattice hierarchy [110]. Our results can be interpreted as a stochastic counterpart to these findings, wherein ALE_β is characterized as the pole evolution of the SDE (1.9).

1.4. Convergence framework

Given the characterization presented in 1.6, to prove convergence to ALE_β, it suffices to

(1)

establish the tightness of the Stieltjes transforms of the empirical particle density at the microscopic scale, and
(2)

verify that the scaling limit is Airy-like, and satisfies the SDE (1.9).

As a demonstration of this approach, we prove the convergence to ALE_β for several continuous interacting particle systems. We next give the formal statement of our result.

We use a strong topology of uniform in compact convergence for line ensembles. More precisely, for a sequence of ordered families of functions $\{f_{i}^{(1)}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}},\{f_{i}^{(2)}(t)\}_{i\in% \mathbb{N},t\in{\mathbb{R}}},\ldots$ , they converge to $\{f_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ under the uniform in compact topology, if for each $i\in\mathbb{N}$ , $\lim_{n\to\infty}f_{i}^{(n)}=f_{i}$ uniformly in any compact interval.

Theorem 1.7.

ALE_β is the edge scaling limit of stationary DBM with certain general potentials (satisfying 7.1 below), stationary Laguerre process, and stationary Jacobi process, all with parameter $\beta$ , under the uniform in compact topology. We refer to 7.2 for a more detailed statement.

The definitions and background of these processes, as well as the precise statement and proof, will be given in Section 7. We emphasize that these convergence results are new even for the classical cases of $\beta=1,4$ (except for the DBM with $\beta=1$ ), which can be viewed as the joint convergence of eigenvalues of time-evolved classical ensembles with real or quaternion entries.

We remark that the developed framework can also be applied to prove convergence to ALE_β for the other mentioned models. The main remaining task is to establish the desired tightness given in (1). While such tightness are not available from [26], a plausible way is to utilize the dynamical loop equation, as in [75] where local laws down to any mesoscopic scale have been proven for random tilings. We leave this for future works.

1.5. Other properties

In addition to proving convergence to ALE_β, our new characterization can be leveraged to further investigate its properties. First, we can study the regularity of ALE_β. The Brownian regularity for the ALE has been intensively studied in [35, 69, 68, 36]. For ALE_β with $\beta\geq 1$ , it has been established in [97] that the lines of ALE_β are locally Brownian. In Section 4 we show that the lines of ALE_β are Hölder continuous with an exponent $1/2$ for any $\beta>0$ . The second property we study is the collision of lines. For $\beta\geq 1$ , it has been established in [97] that the lines of ALE_β do not collide. Conversely, for $\beta\in(0,1)$ , collisions among lines are anticipated. We prove in Section 5 that the occurrence of collisions is almost surely of measure zero.

1.6. Proof ideas

We give an outline of our proofs, highlighting the main difficulties and ideas.

To prove the uniqueness in law as stated in 1.6, the overall strategy is to establish a certain sense of ‘mixing in time’ of the dynamics (given by 1.5). More precisely, we take two families of random particle-generated Nevanlinna functions, both satisfying the two assumptions. Using the SDE (1.9) we reconstruct the dynamics of the poles, which are ‘infinite dimensional DBM’ in a certain sense. We couple the two ‘infinite dimensional DBM’ obtained from both functions, by coupling the driven Brownian motions. Then we show that the poles get closer in time under this coupling. Thus since both dynamics start from time $-\infty$ , necessarily they are the same.

For both the reconstruction of DBM and the coupling, an essential input is that the poles have Airy-zero approximation, uniformly in time. This is implied by the uniform in time Airy-like property, as explained in Remark 1.3. Then under 1.4, it remains to show that such an approximation propagates in time, for which we again resort to the SDE (1.9).

In summary, three tasks are inline: propagation of Airy-zero approximation, reconstruction of DBM, and coupling. We next explain each of them in more details.

1.6.1. Propagation of Airy-zero approximation

Our 1.4 concerns specific times, implying that the $i$ -th pole remains constant away from the $i$ -th zero of the Airy function. Utilizing the SDE (1.9), we manage to get refined estimates: the $i$ -th pole approximates the $i$ -th zero with a polynomially small error over arbitrarily long time intervals with high probability, as demanded for later steps. To achieve this, in Section 3, we analyze (1.9) along certain characteristics which offset the singularity of the nonlinear term. This idea has previously been used (see e.g., [76, 1, 24]) to study DBM down to any mesoscopic scale, where the distance from the spectral parameter $w$ to the particles is much bigger than particle fluctuations. However, in our analysis of (1.9), we operate at a microscopic scale, where the distance from $w$ to the poles is of the same order as their fluctuation size. While a straightforward union bound over characteristic flows from polynomially many points suffices at the mesoscopic scale, our case demands careful selection of characteristic flows and precise estimation of error terms in the SDE, tailored to their initial positions.

1.6.2. DBM reconstruction

As already mentioned, there are significant challenges in analyzing DBM due to the singular repulsive interaction and possible particle collisions, in particular when $\beta\in(0,1)$ . Even for finite dimensional DBM, establishing the existence and uniqueness of a strong solution requires the theory of multivalued SDE, see [33, 32]. Our approach through pole evolutions circumvents these issues entirely. Notably, there are no singularities even when poles collide.

On the other hand, a key challenge of our method lies in reconstructing the dynamics of each pole, which requires ruling out the possibility of poles adhering to each other for prolonged periods. To address these, in Section 4 we first establish that the trajectory of each pole is Hölder continuous solely utilizing (1.9). Together with the Airy-zero approximation, for any short time interval, we can identify a large index $k$ such that the $k$ -th and $(k+1)$ -th poles remain bounded away from each other. This enables us to localize the system and study the evolution of the first $k$ poles, treating the remaining poles’ influence as an additional potential. For this $k$ poles system, in Section 5, by employing classical Itô calculus on certain elementary symmetric polynomials, we show that the time of collisions almost surely has measure zero. We note that similar ideas have been employed to show the instant diffraction of the particles for DBM [66]. Subsequently, the evolution of each pole can be reconstructed by performing a contour integral of (1.9). As a result, the $k$ poles system can be interpreted as a $k$ -dimensional DBM with a time-dependent random drift, which exhibits a monotonicity property.

1.6.3. Uniqueness via coupling

In Section 6, we take two solutions to (1.9), and design a coupling where the poles get closer in time. Consider the $k$ -dimensional DBMs with random drifts reconstructed in the previous step, for these two solutions respectively. Our coupling is by using the same set of driven Brownian motions for both. Note that such $k$ -dimensional DBMs with random drifts are constructed with random $k$ , and only for a short time interval; but we need a coupling for a long time (to let the poles get closer). A trick here is to concatenate these short intervals, and allow for different $k$ in each of them, as long as $k$ is always large enough.

There is a monotonicity property: if the $i$ -th pole of the initial data for the first solution dominates that of the second solution for each $i$ , then at any time after the $i$ -th pole of the first solution dominates that of the second solution for each $i$ . Then we can sandwich one solution between affine shifts of the other, while keeping the error arising from the affine shifts arbitrarily small. Such sandwiching forces the poles of the two solutions to get closer in time. By taking long enough time intervals, they must coincide, establishing the uniqueness as desired. Such coupling and sandwiching strategies have been used to establish local statistics universality in random lozenge tilings [3, 5, 77].

Notations

In the rest of this paper, for any $a\leq b\in{\mathbb{R}}$ , we let $\llbracket a,b\rrbracket=[a,b]\cap\mathbb{Z}$ . For any $w\in{\mathbb{C}}$ , we use $\operatorname{\mathcal{O}}(w)$ to denote some $w^{\prime}\in{\mathbb{C}}$ , satisfying $|w^{\prime}|<C|w|$ for some universal constant $C>0$ . We also write $w^{\prime}\lesssim w$ for $w^{\prime}=\operatorname{\mathcal{O}}(w)$ .

Acknowledgement

The research of J.H. is supported by NSF grant DMS-2331096 and DMS-2337795, and the Sloan research award. The research of L.Z. is supported by NSF grant DMS-2246664 and partially by the Miller Institute for Basic Research in Science. Part of this project was done when L.Z. was visiting University of Pennsylvania in the spring of 2023, and he thanks them for their hospitality. The authors would like to thank Paul Bourgade, Vadim Gorin, Benjamin Landon, and Bálint Virág for helpful discussions.

2. Preliminaries and decomposition

In this section, we set up some preliminaries of our arguments.

We start with an explicit expression for any Airy-like Nevanlinna $Y$ from Definition 1.2. The expression involves the Airy function, which is usually denoted by $\operatorname{Ai}$ , and is a special function that appears in various areas of mathematics and physics. It can be defined as an entire function, and the solution to the Airy equation: $\operatorname{Ai}^{\prime\prime}(w)-w\operatorname{Ai}(w)=0$ with $\operatorname{Ai}(w)\to 0$ as $w\to\infty$ along ${\mathbb{R}}_{+}$ . All the zeros of $\operatorname{Ai}$ are on the real line, and are all negative, and we denote them as $0>{\mathfrak{a}}_{1}>{\mathfrak{a}}_{2}>{\mathfrak{a}}_{3}>\cdots$ .

We now give the expression.

Proposition 2.1.

For any particle-generated Nevanlinna function $Y:\mathbb{H}\to\mathbb{H}\cup{\mathbb{R}}$ with infinitely many poles $x_{1}\geq x_{2}\geq\cdots$ , if i) $\sup_{i\in\mathbb{N}}|x_{i}-{\mathfrak{a}}_{i}|<\infty$ ; and ii) there exists a sequence of complex numbers $w_{n}\to\infty$ along any direction in $(0,3\pi/4)$ , such that $Y(w_{n})-\sqrt{w_{n}}\rightarrow 0$ , then

(2.1)

\displaystyle Y(w)=\sum_{i=1}^{\infty}\frac{1}{x_{i}-w}-\frac{1}{{\mathfrak{a}% }_{i}}-\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}.

Moreover, (2.1) holds if $Y$ is Airy-like.

The Nevanlinna representation (1.7), for any particle-generated Nevanlinna function $Y$ with poles $P$ (a multiset), can be written as

(2.2)

\displaystyle Y(w)=b+cw+\sum_{x\in P}\frac{1}{x-w}-\frac{x}{1+x^{2}},

where $b,c\in{\mathbb{R}}$ , $c\geq 0$ . We remark that it is possible that $P$ contains only finitely many numbers, and the summation in (2.2) is finite. Then to prove 2.1, it remains to determine $b$ and $c$ in (2.2) for $Y_{t}$ , and establish that the sum $\sum_{i=1}^{\infty}\frac{1}{{\mathfrak{a}}_{i}}-\frac{x_{i}}{1+x_{i}^{2}}$ converges.

To start with, we first collect some basic estimates on the Airy function $\operatorname{Ai}$ and Nevanlinna functions, which will also be used repeatedly in the rest of this paper.

2.1. Airy function

The Airy function has the following asymptotic formula. For $|\arg(w)|<\pi$ ,

\operatorname{Ai}(w)\sim\frac{\exp(-\zeta)}{w^{1/4}}\sum_{n=0}^{\infty}\frac{% \Gamma(n+5/6)\Gamma(n+1/6)}{4\pi^{3/2}n!(-2\zeta)^{n}},

\operatorname{Ai}^{\prime}(w)\sim-w^{1/4}\exp(-\zeta)\sum_{n=0}^{\infty}\frac{% 1+6n}{1-6n}\cdot\frac{\Gamma(n+5/6)\Gamma(n+1/6)}{4\pi^{3/2}n!(-2\zeta)^{n}},

where $\zeta=\frac{2}{3}w^{3/2}$ . In particular, there is

\left|w^{1/4}\operatorname{Ai}(w)-\frac{\exp(-2w^{3/2}/3)}{2\sqrt{\pi}}\left(1% -\frac{5}{48w^{3/2}}\right)\right|\leq D(|\arg(w)|)\frac{|\exp(-2w^{3/2}/3)|}{% |w^{3}|},

\left|w^{-1/4}\operatorname{Ai}^{\prime}(w)+\frac{\exp(-2w^{3/2}/3)}{2\sqrt{% \pi}}\left(1+\frac{7}{48w^{3/2}}\right)\right|\leq D(|\arg(w)|)\frac{|\exp(-2w% ^{3/2}/3)|}{|w^{3}|},

for any $w\in{\mathbb{C}}\setminus{\mathbb{R}}_{\leq 0}$ , where $D:[0,\pi)\to{\mathbb{R}}$ is a continuous function. See e.g., [47, Subsection 9.7.iv] and [100, Appendix B]. It follows that

(2.3)

\displaystyle\left|\frac{\operatorname{Ai}^{\prime}(w)}{\operatorname{Ai}(w)}+% \sqrt{w}\right|\lesssim|w|^{-1},

for any $w\in{\mathbb{C}}$ with $|w|$ large enough and $|\arg(w)|<3\pi/4$ . The Weierstrass representation gives

(2.4)

\displaystyle-\frac{\operatorname{Ai}^{\prime}(w)}{\operatorname{Ai}(w)}=\sum_% {i=1}^{\infty}\frac{1}{{\mathfrak{a}}_{i}-w}-\frac{1}{{\mathfrak{a}}_{i}}-% \frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}

It is also known that ${\mathfrak{a}}_{i}$ is around $-(3i\pi/2)^{2/3}$ . More precisely, for any $i\in\mathbb{N}$ we have

(2.5)

\displaystyle\left|{\mathfrak{a}}_{k}+\left(\frac{3\pi i}{2}\right)^{2/3}% \right|\lesssim i^{-1/3}.

2.2. Estimates on Nevanlinna functions

We now present some estimates on (particle-generated) Nevanlinna functions, which will be used in the Airy-like function part of 2.1. We note that some of them are also used repeatedly in subsequent sections.

For any particle-generated Nevanlinna function $Y:\mathbb{H}\to\mathbb{H}\cup{\mathbb{R}}$ , from (2.2) we have

(2.6)

\displaystyle{\rm Im}[Y(w)]=c{\rm Im}[w]+\sum_{x\in P}\frac{{\rm Im}[w]}{|x-w|% ^{2}}.

Lemma 2.2.

The quantity ${\rm Im}[w]{\rm Im}[Y(w)]$ is monotone in ${\rm Im}[w]$ ; the derivatives of $Y$ satisfy

(2.7)

\displaystyle|Y^{(k)}(w)|\leq\frac{k!{\rm Im}[Y(w)]}{{\rm Im}[w]^{k}}\leq\frac% {k!|Y(w)|}{{\rm Im}[w]^{k}},

where $Y^{(k)}$ is the $k$ -th derivative of $Y$ , for any integer $k\geq 2$ .

Proof.

We first consider $k=1$ . Denote $w=E+\mathrm{i}\eta$ , then (2.6) gives

(2.8)

\displaystyle{\rm Im}[w]{\rm Im}[Y(w)]=c\eta^{2}+\sum_{x\in P}\frac{\eta^{2}}{% |E-x|^{2}+\eta^{2}},

which is increasing in $\eta\geq 0$ . Using (2.2), the derivative of $Y(w)$ satisfies

(2.9)

\displaystyle|Y^{\prime}(w)|=\left|c+\sum_{x\in P}\frac{1}{(x-w)^{2}}\right|% \leq c+\sum_{x\in P}\frac{1}{|x-w|^{2}}=\frac{{\rm Im}[Y(w)]}{{\rm Im}[w]}.

And by

|Y^{(k)}(w)|\leq\sum_{x\in P}\frac{k!}{|x-w|^{k+1}}\leq\frac{k!|Y^{\prime}(w)|% }{{\rm Im}[w]^{k-1}},

the $k\geq 2$ case follows. ∎

As already alluded to, if a particle-generated Nevanlinna function $Y$ is close to the function $\sqrt{w}$ , its poles would be close to the Airy function zeros. More precisely, we have the following estimate.

Lemma 2.3.

Take any parameters $K>100$ and $0\leq\delta<1$ . Suppose that a particle-generated Nevanlinna function $Y$ satisfies the following conditions:

•

there is no pole of $Y$ in $(K,\infty)$ ;

•

for any $w=x+\mathrm{i}y$ with $x\leq K^{2}$ and $y\geq 4K^{2}/(1+|x|^{\delta/2})$ , we have

\left|Y(w)-\sqrt{w}\right|\leq\frac{{\rm Im}[\sqrt{w}]^{1-\delta}}{{\rm Im}[w]}.

Then $Y$ has infinitely many poles $x_{1}\geq x_{2}\geq\cdots$ , and $|x_{i}-{\mathfrak{a}}_{i}|<CK^{4}i^{-\delta/6}$ for any $i\in\mathbb{N}$ , where $C>0$ is a universal constant.

In particular, these conditions are satisfied by Airy-like Nevanlinna functions (with $K$ large and $\delta={\mathfrak{d}}$ ).

Corollary 2.4.

For any $({\mathfrak{d}},C_{*})$ -Airy-like Nevanlinna function $Y$ (with poles $x_{1}\geq x_{2}\geq\cdots$ ), there exists $B>0$ depending only on $C_{*}$ , such that $|x_{i}-{\mathfrak{a}}_{i}|\leq B$ for each $i\in\mathbb{N}$ .

The proof of 2.3 relies on the Helffer-Sjöstrand formula, which has become standard in random matrix theory. Therefore, we defer it to Appendix A.

2.3. Proof of 2.1

Thanks to (2.3) and (2.4), we have

(2.10)

\displaystyle\sum_{i=1}^{\infty}\frac{1}{{\mathfrak{a}}_{i}-w}-\frac{1}{{% \mathfrak{a}}_{i}}-\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}-% \sqrt{w}\rightarrow 0,

when $w\to\infty$ along any direction in $(-3\pi/4,3\pi/4)$ . By $|{\mathfrak{a}}_{i}|\sim(3\pi i/2)^{2/3}$ from (2.5), we have

(2.11)

\displaystyle|{\mathfrak{a}}_{i}+(3\pi i/2)^{2/3}|,|x_{i}+(3\pi i/2)^{2/3}|\leq B

for a large $B>0$ ; in particular $x_{1}\leq B$ . If we take $w$ with $\arg(w)\in(-3\pi/4,3\pi/4)$ and $|w|>2B$ ,

(2.12)

\displaystyle\begin{split}&\phantom{{}={}}\left|\sum_{i=1}^{\infty}\frac{1}{x_% {i}-w}-\frac{1}{{\mathfrak{a}}_{i}}-\sum_{i=1}^{\infty}\frac{1}{{\mathfrak{a}}% _{i}-w}-\frac{1}{{\mathfrak{a}}_{i}}\right|\\ &\leq\sum_{i=1}^{\infty}\frac{|{\mathfrak{a}}_{i}-x_{i}|}{|x_{i}-w||{\mathfrak% {a}}_{i}-w|}\leq\sum_{i=1}^{\infty}\frac{B}{|x_{i}-w||{\mathfrak{a}}_{i}-w|}\\ &=\sum_{i>|w|^{3/2}}\frac{B}{|x_{i}-w||{\mathfrak{a}}_{i}-w|}+\sum_{i=1}^{% \lfloor|w|^{3/2}\rfloor}\frac{B}{|x_{i}-w||{\mathfrak{a}}_{i}-w|}\\ &\lesssim\sum_{i>|w|^{3/2}}\frac{B}{i^{4/3}}+\sum_{i=1}^{\lfloor|w|^{3/2}% \rfloor}\frac{B}{|w|^{2}}\lesssim\frac{1}{|w|^{1/2}},\end{split}

where in the last line we used that, when $i>|w|^{3/2}$ , $|x_{i}-w||{\mathfrak{a}}_{i}-w|\gtrsim i^{4/3}$ for the first term, and $|x_{i}-w||{\mathfrak{a}}_{i}-w|\gtrsim|w|^{2}$ for the second term (since $|w|>2B$ ).

Therefore, (2.10) and (LABEL:e:replace) together give that

(2.13)

\displaystyle\sum_{i=1}^{\infty}\frac{1}{x_{i}-w}-\frac{1}{{\mathfrak{a}}_{i}}% -\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}-\sqrt{w}% \rightarrow 0,

for $w\to\infty$ along any direction in $(-3\pi/4,3\pi/4)$ . Also note that, thanks to (2.11), we have

(2.14)

\displaystyle\left|\sum_{i=1}^{\infty}\frac{1}{{\mathfrak{a}}_{i}}-\frac{x_{i}% }{1+x_{i}^{2}}\right|\leq\sum_{i=1}^{\infty}\left|\frac{1+x_{i}(x_{i}-{% \mathfrak{a}}_{i})}{{\mathfrak{a}}_{i}(1+x_{i}^{2})}\right|\leq\sum_{i=1}^{% \infty}\frac{1+|{\mathfrak{a}}_{i}|B+B^{2}}{|{\mathfrak{a}}_{i}|(1+(({% \mathfrak{a}}_{i}+B)\wedge 0)^{2})}<\infty.

By plugging (2.14) into the representation (2.2) for $Y$ , we can rewrite $Y$ as (for some $b^{\prime}\in{\mathbb{R}}$ )

(2.15)

\displaystyle Y(w)=b+cw+\sum_{i=1}^{\infty}\frac{1}{x_{i}-w}-\frac{x_{i}}{1+x_% {i}^{2}}=b^{\prime}+cw+\sum_{i=1}^{\infty}\frac{1}{x_{i}-w}-\frac{1}{{% \mathfrak{a}}_{i}}-\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}.

By our assumption, $Y(w_{n})-\sqrt{w}_{n}\rightarrow 0$ as $n\to\infty$ . Taking $w=w_{n}$ in (2.15), comparing with (2.13), we conclude that $b^{\prime}+cw_{n}\rightarrow 0$ when $n\to\infty$ . It follows that $b^{\prime}=c=0$ , and (2.1) holds.

Finally, if $Y$ is Airy-like, then 2.4 implies that $\sup_{i\in\mathbb{N}}|x_{i}-{\mathfrak{a}}_{i}|<\infty$ , and $|Y(n\mathrm{i})-\sqrt{n\mathrm{i}}|\to 0$ as $n\to\infty$ . These verify the assumptions in 2.1, and (2.1) holds. ∎

2.4. Domain extension

As stated in Remark 1.3, for an Airy-like Nevanlinna function $Y$ , we also provide a bound of $|Y(w)-\sqrt{w}|$ for $w$ close to ${\mathbb{R}}_{+}$ , which will be useful later.

Lemma 2.5.

For $Y:\mathbb{H}\to\mathbb{H}\cup{\mathbb{R}}$ being any $({\mathfrak{d}},C_{*})$ -Airy-like Nevanlinna function, we have

|Y(w)-\sqrt{w}|\leq B|w|^{-1/2},\quad\forall\arg(w)\in(0,3\pi/4),|w|>B.

for $B>0$ depending only on ${\mathfrak{d}}$ and $C_{*}$ .

Proof.

By 2.4, there is $B^{\prime}>0$ with each $|x_{i}-{\mathfrak{a}}_{i}|\leq B^{\prime}$ . Then, for any $w\in\mathbb{H}$ with $\arg(w)\in(0,3\pi/4),|w|>2B^{\prime}$ , we have $|{\mathfrak{a}}_{i}-w|\leq 2|x_{i}-w|$ .

By 2.1 and (2.4), we have

\left|Y(w)+\frac{\operatorname{Ai}^{\prime}(w)}{\operatorname{Ai}(w)}\right|% \leq\sum_{i=1}^{\infty}\left|\frac{1}{{\mathfrak{a}}_{i}-w}-\frac{1}{x_{i}-w}% \right|=\sum_{i=1}^{\infty}\left|\frac{x_{i}-{\mathfrak{a}}_{i}}{(x_{i}-w)({% \mathfrak{a}}_{i}-w)}\right|\leq\sum_{i=1}^{\infty}\frac{2B^{\prime}}{|{% \mathfrak{a}}_{i}-w|^{2}}.

Using (2.5), we have that

\sum_{i=1}^{\infty}\frac{1}{|{\mathfrak{a}}_{i}-w|^{2}}=\sum_{i=1}^{\lfloor|w|% ^{3/2}\rfloor}\frac{1}{|{\mathfrak{a}}_{i}-w|^{2}}+\sum_{i=\lfloor|w|^{3/2}% \rfloor+1}^{\infty}\frac{1}{|{\mathfrak{a}}_{i}-w|^{2}}\lesssim\frac{|w|^{3/2}% }{|w|^{2}}+\sum_{i=\lfloor|w|^{3/2}\rfloor+1}^{\infty}\frac{1}{i^{4/3}}% \lesssim|w|^{-1/2}.

Thus with (2.3), the conclusion follows. ∎

2.5. Topological statements

As we shall derive convergence to ALE_β from Stieltjes transforms, we will need several statements on the functional spaces, which we provide here.

Definition 2.6.

For any locally finite measures $\mu_{1},\mu_{2},\cdots$ on ${\mathbb{R}}$ , we say that they converge in the vague topology to another locally finite $\mu$ , if $\mu_{n}(f)\to\mu(f)$ , for any $f:{\mathbb{R}}\to{\mathbb{R}}$ that is compactly supported and smooth.

Such vague topology arises naturally from Nevanlinna function convergence.

Lemma 2.7.

Take Nevanlinna functions $Y_{1},Y_{2},\cdots$ and $Y$ such that $Y_{n}\to Y$ as $n\to\infty$ , uniformly in any compact subset of $\mathbb{H}$ . Suppose the corresponding measures (in their Nevanlinna representation) are $\mu_{1},\mu_{2},\cdots$ and $\mu$ , respectively, then $\mu_{n}\to\mu$ in the vague topology.

Proof.

Take any $f$ that is compactly supported and smooth, and let $K$ be a large enough number such that $f=0$ outside $[-K,K]$ . Take a smooth function $\chi:{\mathbb{R}}_{+}\to{\mathbb{R}}$ , such that $\chi=1$ on $(0,1)$ , and $\chi=0$ on $(2,\infty)$ . By A.1, we have

\mu_{n}(f)=\frac{1}{\pi}\int_{x+\mathrm{i}y\in\mathbb{H}}-{\rm Re}[Y_{n}(x+% \mathrm{i}y)]yf^{\prime}(x)\chi^{\prime}(y)-{\rm Im}[Y_{n}(x+\mathrm{i}y)](yf^% {\prime\prime}(x)\chi(y)+f(x)\chi^{\prime}(y)){\rm d}x{\rm d}y.

We note that $yf^{\prime}(x)\chi^{\prime}(y)=f(x)\chi^{\prime}(y)=0$ whenever $(x,y)\not\in[-K,K]\times[1,2]$ . Also, $yf^{\prime\prime}(x)\chi(y)=0$ whenever $(x,y)\not\in[-K,K]\times(0,2]$ ; and for $y\leq 2$ , from Nevalinna representation we have that ${\rm Im}[Y_{n}(x+\mathrm{i}y)]\leq\frac{2}{y}{\rm Im}[Y_{n}(x+2\mathrm{i})]$ . Therefore, we have that the integrand in the above integral is non-zero only in $[-K,K]\times(0,2]$ ; and it is bounded by a constant, which is independent of $n$ by the uniform convergence of $Y_{n}$ in $[-K,K]\times[1,2]$ . Thus we can apply dominated convergence theorem to deduce that the above integral converges to

\mu(f)=\frac{1}{\pi}\int_{x+\mathrm{i}y\in\mathbb{H}}-{\rm Re}[Y(x+\mathrm{i}y% )]yf^{\prime}(x)\chi^{\prime}(y)-{\rm Im}[Y(x+\mathrm{i}y)](yf^{\prime\prime}(% x)\chi(y)+f(x)\chi^{\prime}(y)){\rm d}x{\rm d}y.

So the conclusion follows. ∎

On the other hand, in the setting of particle-generated measures, vague topology convergence can imply pole convergence.

Lemma 2.8.

For a sequence of particle-generated measures $\mu_{1},\mu_{2},\ldots$ , such that $\mu_{k}\to\mu$ as $k\to\infty$ in the vague topology, the limit $\mu$ must also be particle-generated. Moreover, if there is some $K>0$ such that $\mu_{k}([K,\infty))=0$ for each $k$ , then the followings are true. We denote by $x^{k}_{i}$ the $i$ -th largest pole of $\mu_{k}$ (with the convention that $x^{k}_{i}=-\infty$ if there are less than $i$ poles). For each $i\in\mathbb{N}$ , either $x^{k}_{i}\to-\infty$ as $k\to\infty$ , or $\lim_{k\to\infty}x^{k}_{i}$ exists and is a pole of $\mu$ . Also, all the poles of $\mu$ are given by such limits.

Proof.

By vague topology convergence, for any $a<b$ , we have that $\limsup_{k\to\infty}\mu_{k}([a,b])\leq\mu([a,b])$ , and $\liminf_{k\to\infty}\mu_{k}((a,b))\geq\mu((a,b))$ . Then for any $a<b$ with $\mu([a,b])<1$ , we must have $\mu_{k}([a,b])=0$ for $k$ large enough, since each $\mu_{k}([a,b])$ must be an integer. Therefore, $\mu_{k}((a,b))=0$ for $k$ large enough, and $\mu((a,b))=0$ . These imply that $\mu$ in any open interval is either $\geq 1$ or zero. Therefore, in any compact interval, $\mu$ is supported on finitely many points.

Now take any $x$ where $\mu(\{x\})>0$ . Take any $a<x<b$ such that $\mu((a,b))=\mu([a,b])=\mu(\{x\})$ . Then for $k$ large enough, $\mu_{k}((a,b))=\mu_{k}([a,b])=\mu(\{x\})$ ; so $\mu(\{x\})\in\mathbb{N}$ . These imply that $\mu$ is particle-generated.

Under the additional assumption, there is $\mu((K,\infty))=0$ . So we can write the poles of $\mu$ as $x_{1}\geq x_{2}\geq\cdots$ (with the convention that $x_{i}=-\infty$ if there are less than $i$ poles). Then we can show $x^{k}_{i}\to x_{i}$ via induction in $i\in\mathbb{N}$ , using that for each $a<b$ with $a,b\not\in\{x_{i}\}_{i=1}^{\infty}$ , $\mu_{k}((a,b))=\mu((a,b))$ for any $k$ large enough. ∎

3. Pole evolution: uniform rigidity in time

In this section, we prove a uniform in time estimate for the poles. More precisely, the following proposition states that for the SDE (1.9), with high probability, its pole evolution gives a line ensemble (i.e., all the poles are bounded from above, and the trajectories are continuous), and the poles are close to the zeros of the Airy function, uniformly in time.

Proposition 3.1.

For any ${\mathfrak{d}},C_{*}>0$ , there exist small $\delta,c>0$ and large $C>0$ , such that the following holds. Take any particle-generated $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.5, and large $T>0$ . Conditional on the event that $Y_{0}$ is $({\mathfrak{d}},C_{*})$ -Airy-like, with probability at least $1-e^{-c(\log T)^{2}}$ , we have

(1)

The poles of $\{Y_{t}\}_{t\in[T,2T]}$ give a line ensemble $\{x_{i}(t)\}_{i\in\mathbb{N},t\in[T,2T]}$ , and for each $t\in[T,2T]$ and $w\in\mathbb{H}$ ,

(3.1)

\displaystyle Y_{t}(w)=\sum_{i=1}^{\infty}\frac{1}{x_{i}(t)-w}-\frac{1}{{% \mathfrak{a}}_{i}}-\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}.

(2)

For each $t\in[T,2T]$ , $i\in\mathbb{N}$ , $y\leq C$ , we have that

(3.2)

|x_{i}(t)-{\mathfrak{a}}_{i}|\leq\frac{C(\log T)^{40}}{i^{\delta}},

and

(3.3)

|\{i\in\mathbb{N}:x_{i}(t)\in[y-1,y+1]\}|\leq C\sqrt{|y|+1}.

This immediately implies the first part of 1.6.

Corollary 3.2.

For any $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.4 and 1.5, its poles give a line ensemble $\{x_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ , and $Y_{t}(w)=\sum_{i=1}^{\infty}\frac{1}{x_{i}(t)-w}-\frac{1}{{\mathfrak{a}}_{i}}-% \frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}$ for any $t\in{\mathbb{R}}$ .

We note that (3.3) (which is a Wegner estimate) follows easily from (3.2), plus

|\{i\in\mathbb{N}:{\mathfrak{a}}_{i}(t)\in[y-1,y+1]\}|\leq C\sqrt{|y|+1},

which directly follows from (2.5).

Our general strategy is to obtain uniform in time estimates for $Y_{t}-\sqrt{w}$ , and to apply 2.3. The main tasks are (1) to estimate bulk pole densities, via bounding $|Y_{t}(w)-\sqrt{w}|$ for $w$ in a reasonable domain contained in $\mathbb{H}$ (in particular, allowing for polynomially close to the real axis, when ${\rm Re}[w]\to-\infty$ ); (2) to bound the first pole $x_{1}$ . Both these are to be achieved through analyzing the SDE (1.9).

3.1. Characteristic flow

We consider the characteristic flow,

(3.4)

\displaystyle\operatorname{\partial}_{t}w_{t}=-\sqrt{w_{t}},\quad w_{0}\in% \mathbb{H},

which can be solved as

(3.5)

\displaystyle\operatorname{\partial}_{t}\sqrt{w_{t}}=-\frac{1}{2},\quad w_{t}=% (w_{0}-t/2)^{2}.

By plugging this characteristic flow into (1.9), itô’s formula gives the following semi-martingale decomposition for $Y_{t}(w_{t})-\sqrt{w_{t}}$

(3.6)

\displaystyle{\rm d}(Y_{t}(w_{t})-\sqrt{w_{t}})=({\rm d}M_{t})(w_{t})+\left(% \frac{2-\beta}{2\beta}\operatorname{\partial}^{2}_{w}Y_{t}(w_{t})+(Y_{t}(w_{t}% )-\sqrt{w_{t}})\operatorname{\partial}_{w}Y_{t}(w_{t})\right){\rm d}t,

with the Martingale term $({\rm d}M_{t})(w_{t})={\rm d}(M_{t}(w_{t}))-(\operatorname{\partial}_{w}M_{t})% (w_{t}){\rm d}t$ , whose quadratic variations are given by (using (1.10)):

(3.7)

\frac{{\rm d}}{{\rm d}t}\left\langle\int_{0}^{t}({\rm d}M_{s})(w_{s})\right% \rangle=\left.\frac{1}{3\beta}\operatorname{\partial}_{w}^{3}Y_{t}(w)\right|_{% w=w_{t}}.

In the rest of this section, we fix ${\mathfrak{d}},C_{*}>0$ , and take $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.5, and (unless otherwise noted) conditional on $Y_{0}$ which is $({\mathfrak{d}},C_{*})$ -Airy-like. All the constants below (including those in $\lesssim$ and $\operatorname{\mathcal{O}}(\cdot)$ ) can depend on ${\mathfrak{d}}$ and $C_{*}$ . We take $T$ to be a large number, and set $K=(\log T)^{8}$ .

3.2. Estimates for the bulk

We next prove the following bound of $|Y_{t}(w)-\sqrt{w}|$ for $w$ in a domain contained in $\mathbb{H}$ . It will be used to bound the bulk pole densities.

Proposition 3.3.

There exist small $\delta,c>0$ such that the following holds. Define the spectral domain

(3.8)

\displaystyle{\mathcal{D}}=\{\kappa+\mathrm{i}\eta:\eta\geq K,\eta^{-1-\delta}% K\leq\kappa\leq\sqrt{K^{2}+\eta^{2}}\}.

With probability $1-e^{-c(\log T)^{2}}$ , for any $t\in[T,2T]$ , and $\sqrt{w}\in{\mathcal{D}}$ , it holds

(3.9)

\displaystyle\left|Y_{t}(w)-\sqrt{w}\right|\leq\frac{{\rm Im}[\sqrt{w}]^{1-% \delta}}{{\rm Im}[w]}.

Thanks to the (away from the real axis) Lipschitz property of $Y_{t}(w)-\sqrt{w}$ in 2.2, we only need to prove (3.9) for a set of carefully chosen mesh points. Namely, we consider the following mesh of points in the upper half plane:

(3.10)

\displaystyle{\mathcal{L}}=\left\{\kappa+\mathrm{i}\eta:\eta^{3}\in\mathbb{Z},% \eta\geq K,\kappa=\frac{\mathbb{Z}}{\eta^{2}},K\leq\kappa\leq T+K+\eta\right\}.

Lemma 3.4.

For any $\kappa^{\prime}+\mathrm{i}\eta^{\prime}\in{\mathcal{D}}$ as defined in (3.8) and $t\in[T,2T]$ , there exists $\kappa+\mathrm{i}\eta\in{\mathcal{L}}$ such that $t\leq 2\kappa-2\eta^{-1-\delta}K$ , and

(3.11)

\displaystyle|\kappa^{\prime}-\kappa+t/2|,|\eta^{\prime}-\eta|\leq\frac{1}{% \eta^{2}}.

Proof.

Suppose that $\eta^{\prime}\in[i^{1/3},(i+1)^{1/3}]$ for some integer $i\geq K^{3}$ , we can take $\eta=i^{1/3}$ , therefore $|\eta-\eta^{\prime}|\leq(i+1)^{1/3}-i^{1/3}\leq i^{-2/3}/3=\eta^{-2}/3$ . Suppose that $\kappa^{\prime}+t/2\in[j\eta^{-2},(j+1)\eta^{-2}]$ for some $j\in\mathbb{Z}$ , we can take $\kappa=(j+1)\eta^{-2}$ . Then $\kappa\geq\kappa^{\prime}+T/2>K$ , and $\kappa\leq\kappa^{\prime}+T+\eta^{-2}\leq T+K+\eta$ , and $|\kappa^{\prime}-\kappa+t/2|\leq\eta^{-2}$ . Also $\kappa-t/2\geq\kappa^{\prime}\geq\eta^{-1-\delta}K$ , thus $t\leq 2\kappa-2\eta^{-1-\delta}K$ . ∎

Now 3.3 follows from the following estimate on one point.

Lemma 3.5.

The following holds true for small enough $\delta,c>0$ . Take any $u=\kappa_{0}+\eta\mathrm{i}\in{\mathcal{L}}$ , and let $\sqrt{w_{t}}=u-t/2=:\kappa_{t}+\eta\mathrm{i}$ for each $t<2\kappa_{0}$ . Conditional on $Y_{0}$ with $|Y_{0}(w_{0})-\sqrt{w_{0}}|\leq{{\rm Im}[\sqrt{w_{0}}]^{1-{\mathfrak{d}}}}/{{% \rm Im}[w_{0}]}$ , with probability $1-e^{-c\sqrt{\eta}}$ ,

(3.12)

\displaystyle|Y_{t}(w_{t})-\sqrt{w_{t}}|\leq\frac{1}{\kappa_{t}\eta^{\delta}},% \quad\forall\;0\leq t\leq 2T\wedge(2\kappa_{0}-2\eta^{-1-\delta}K).

Proof of 3.3.

By a union bound over all the points in ${\mathcal{L}}$ , 3.5 implies that

(3.13)

\displaystyle|Y_{t}(w_{t})-\sqrt{w_{t}}|\leq\frac{1}{\kappa_{t}\eta^{\delta}},% \quad\forall\;0\leq t\leq 2T\wedge(2\kappa_{0}-2\eta^{-1-\delta}K),\;\sqrt{w_{% 0}}\in{\mathcal{L}},

with probability at least

(3.14)

\displaystyle 1-\sum_{\eta\geq K,\eta^{3}\in\mathbb{Z}}Ce^{-c\sqrt{\eta}}\eta^% {2}(2T+2\eta)\geq 1-e^{-c(\log T)^{2}},

where we used that $K=(\log T)^{8}$ and $T$ is large enough, and $C>0$ is a large constant.

For any $t\in[T,2T]$ and $\sqrt{w^{\prime}}=\kappa^{\prime}+\mathrm{i}\eta^{\prime}\in{\mathcal{D}}$ , thanks to 3.4, there exists $\sqrt{w_{0}}=\kappa+\mathrm{i}\eta\in{\mathcal{L}}$ such that $|\kappa^{\prime}-\kappa_{t}|,|\eta^{\prime}-\eta|\leq\eta^{-2}$ . Then we have $\eta^{-1-\delta}K/2<\kappa_{t}<2\eta$ . It also follows that

(3.15)

|\sqrt{w_{t}}-\sqrt{w^{\prime}}|\leq\sqrt{|\kappa^{\prime}-\kappa_{t}|^{2}+|% \eta^{\prime}-\eta|^{2}}\leq\sqrt{2}\eta^{-2},

and

(3.16)

\displaystyle|w_{t}-w^{\prime}|\leq|\sqrt{w_{t}}+\sqrt{w^{\prime}}||\sqrt{w_{t% }}-\sqrt{w^{\prime}}|\leq(2\sqrt{\kappa^{\prime 2}+\eta^{2}}+2\eta^{-2})\sqrt{% 2}\eta^{-2}\lesssim\eta^{-1}.

In particular, this implies that $|w_{t}-w^{\prime}|$ is much smaller than $2\kappa_{t}\eta={\rm Im}[w_{t}]$ (which is at least $\eta^{-\delta}K$ ). It then follows from 2.2 that

(3.17)

\displaystyle|Y_{t}(w_{t})-Y_{t}(w^{\prime})|\lesssim\frac{{\rm Im}\sqrt{w_{t}% }}{{\rm Im}[w_{t}]}|w_{t}-w^{\prime}|=\frac{|w_{t}-w^{\prime}|}{2\kappa_{t}}% \lesssim\frac{1}{\kappa_{t}\eta}.

Combining this with (3.14) and (3.15), and using that $\eta^{-2}$ is much smaller than $\frac{1}{\kappa_{t}\eta^{\delta}}$ (which is at least $\eta^{-1-\delta}/2$ ), we conclude that

|Y_{t}(w^{\prime})-\sqrt{w^{\prime}}|\lesssim\frac{1}{\kappa_{t}\eta^{\delta}}% \lesssim\frac{1}{2\kappa^{\prime}{\eta^{\prime}}^{\delta}}=\frac{{\rm Im}[% \sqrt{w}]^{1-\delta}}{{\rm Im}[w]}.

Then the proof finishes by taking a slightly smaller $\delta$ (to remove the constant factor). ∎

We now prove the one point estimate, using (3.6).

Proof of 3.5.

We introduce the following stopping time,

(3.18)

\displaystyle\sigma=\inf\left\{0\leq t\leq 2T:|Y_{t}(w_{t})-\sqrt{w_{t}}|\geq% \frac{1}{\kappa_{t}\eta^{\delta}},\text{ or }\kappa_{t}\leq\eta^{-1-\delta}K,% \text{ or }t=2T\right\}.

We now bound the terms in the RHS of (3.6). For $0\leq t\leq\sigma$ , we have ${\rm Im}[Y_{t}(w_{t})]\leq 2{\rm Im}[\sqrt{w_{t}}]=2\eta$ , since

(3.19)

\displaystyle|{\rm Im}[Y_{t}(w_{t})]-\eta|\leq\frac{1}{\kappa_{t}\eta^{\delta}% }\leq\frac{\eta}{K}.

Therefore (using 2.2) we get

(3.20)

\displaystyle|\operatorname{\partial}_{w}^{2}Y_{t}(w_{t})|\leq\frac{2{\rm Im}[% Y_{t}(w_{t})]}{{\rm Im}[w_{t}]^{2}}\leq\frac{4{\rm Im}[\sqrt{w_{t}}]}{{\rm Im}% [w_{t}]^{2}}=\frac{1}{\kappa_{t}^{2}\eta}.

The quadratic variation of the martingale term is given by (3.7); then (again using 2.2)

\left|\frac{{\rm d}}{{\rm d}t}\left\langle\int_{0}^{t}({\rm d}M_{s})(w_{s})% \right\rangle\right|\leq\frac{2{\rm Im}[Y_{t}(w_{t})]}{\beta{\rm Im}[w_{t}]^{3% }}\leq\frac{4{\rm Im}[\sqrt{w_{t}}]}{\beta{\rm Im}[w_{t}]^{3}}=\frac{1}{2\beta% \kappa_{t}^{3}\eta^{2}}.

By integrating in time, we have

\int_{0}^{t\wedge\sigma}\frac{{\rm d}s}{2\beta\kappa_{t}^{3}\eta^{2}}\leq\frac% {1}{2\beta\kappa_{t}^{2}\eta^{2}}.

Take $\theta=\sqrt{\eta}$ . By the Burkholder-Davis-Gundy inequality, for any $0\leq t\leq 2T$ the following holds with probability⁴⁴4Here and in the rest of this proof, $C$ is used to denote a large constant, whose value may change from line to line. $\geq 1-Ce^{-c\theta}$ :

\sup_{0\leq s\leq t\wedge\sigma}\left|\int_{0}^{s}({\rm d}M_{u})(w_{u})\right|% \leq\frac{\theta}{2\kappa_{t}\eta}.

Then we can take a union bound over times

t=\left(2-\frac{1}{2^{i-1}}\right)\kappa_{0},\quad\kappa_{t}=\frac{\kappa_{0}}% {2^{i}},

for $i\in\llbracket 1,3\log_{2}(T)\rrbracket$ , and get that with probability $1-C\log(T)e^{-c\theta}$ ,

(3.21)

\displaystyle\left|\int_{0}^{t}({\rm d}M_{s})(w_{s})\right|\leq\frac{\theta}{% \kappa_{t}\eta},\quad\forall\;0\leq t\leq\sigma.

Now by (3.6), and using (3.20) and (3.21), we have that for any $0\leq t\leq\sigma$ ,

(3.22)

\displaystyle\left|Y_{t}(w_{t})-\sqrt{w_{t}})\right|\leq\int_{0}^{t}\left|Y_{s% }(w_{s})-\sqrt{w_{s}}\right|\left|\partial_{w}Y_{s}(w_{s})\right|{\rm d}s+% \frac{2\theta}{\kappa_{t}\eta}+\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|.

Here we used that $\int_{0}^{t}\frac{2}{\kappa_{s}^{2}\eta}{\rm d}s<\frac{4}{\kappa_{t}\eta}$ , which is much smaller than $\frac{\theta}{\kappa_{t}\eta}$ .

For $0\leq t\leq\sigma$ , (using 2.2) we have

|\operatorname{\partial}_{w}Y_{t}(w_{t})|\leq\frac{{\rm Im}[Y_{t}(w_{t})]}{{% \rm Im}[w_{t}]}\leq\frac{1}{2\kappa_{t}}+\frac{|Y_{t}(w_{t})-\sqrt{w_{t}}|}{2% \kappa_{t}\eta}\leq\frac{1}{2\kappa_{t}}+\frac{1}{2\kappa_{t}^{2}\eta^{1+% \delta}}=:\gamma(t)\leq\frac{1}{\kappa_{t}}.

Then for any $0\leq s<t\leq\sigma$ ,

(3.23)

\displaystyle\int_{s}^{t}\gamma(u){\rm d}u\leq\log(\kappa_{s}/\kappa_{t})+% \frac{1}{\kappa_{t}\eta^{1+\delta}}\leq\log(\kappa_{s}/\kappa_{t})+\frac{1}{K}.

By Grönwall’s inequality and (3.22), for any $0\leq t\leq\sigma$ , we can bound $|Y_{t}(w_{t})-\sqrt{w_{t}}|$ by:

(3.24)

\displaystyle\frac{2\theta}{\kappa_{t}\eta}+\left|Y_{0}(w_{0})-\sqrt{w_{0}}% \right|+\int_{0}^{t}\gamma(s)\left(\frac{2\theta}{\kappa_{s}\eta}+\left|Y_{0}(% w_{0})-\sqrt{w_{0}}\right|\right)\exp\left(\int_{s}^{t}\gamma(u){\rm d}u\right% ){\rm d}s.

We note that by (3.23), $\exp(\int_{s}^{t}\gamma(u){\rm d}u)\leq\frac{2\kappa_{s}}{\kappa_{t}}$ . Also

\int_{0}^{t}\gamma(s)\frac{\theta}{\kappa_{s}\eta}\cdot\frac{\kappa_{s}}{% \kappa_{t}}{\rm d}s\leq\frac{\theta}{\kappa_{t}\eta}\int_{0}^{t}\frac{1}{% \kappa_{s}}{\rm d}s=\frac{2\log(\kappa_{0}/\kappa_{t})\theta}{\kappa_{t}\eta},

and

\int_{0}^{t}\gamma(s)\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|\frac{\kappa_{s}}{% \kappa_{t}}{\rm d}s\leq\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|\int_{0}^{t}\frac% {1}{\kappa_{t}}{\rm d}s=\frac{t\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|}{\kappa_% {t}}.

Therefore we have

\left|Y_{t}(w_{t})-\sqrt{w_{t}}\right|\leq\frac{2(1+4\log(\kappa_{0}/\kappa_{t% }))\theta}{\kappa_{t}\eta}+\left(1+\frac{2t}{\kappa_{t}}\right)\left|Y_{0}(w_{% 0})-\sqrt{w_{0}}\right|\leq\frac{1}{2\kappa_{t}\eta^{\delta}},

where in the last inequality we used $\log(\kappa_{0}/\kappa_{t})\leq\log(T\eta^{2})$ , which is much smaller than $\eta^{1/2-\delta}$ provided $K$ is large enough; and

(3.25)

\displaystyle\left(1+\frac{2t}{\kappa_{t}}\right)\left|Y_{0}(w_{0})-\sqrt{w_{0% }}\right|\leq\left(1+\frac{2t}{\kappa_{t}}\right)\frac{{\rm Im}[\sqrt{w_{0}}]^% {1-{\mathfrak{d}}}}{{\rm Im}[w_{0}]}=\left(1+\frac{2t}{\kappa_{t}}\right)\frac% {1}{(2\kappa_{t}+t)\eta^{{\mathfrak{d}}}}\leq\frac{1}{3\kappa_{t}\eta^{\delta}}.

Therefore, since $Y_{t}(w_{t})-\sqrt{w_{t}}$ is continuous in $t$ , we conclude that $\sigma=2T\wedge(2\kappa-2\eta^{-1-\delta}K)$ with probability $1-Ce^{-c\sqrt{\eta}}$ . ∎

3.3. Estimates for the first pole

Under the same setup as the previous subsection (in particular, $T$ is taken to be any large number, and $K=(\log T)^{8}$ ), we next prove the following proposition, which states that with high probability all the poles are bounded by $K$ .

Proposition 3.6.

There exists a small number $c>0$ such that the following holds. With probability $1-e^{-c(\log T)^{4}}$ , for any $t\in[0,2T]$ , $Y_{t}$ has no pole in $(K,\infty)$ .

We introduce the following stopping time, which is the first time the largest pole exceeds $K$ ,

(3.26)

\displaystyle\sigma_{0}=\inf\left\{0\leq t\leq 2T:Y_{t}\text{ has a pole in }(% K,\infty)\text{, or }t=2T\right\}.

We now denote $K_{j}=jK$ for $j\in\mathbb{N}$ (in particular $K_{1}=K$ ), and consider a mesh of points:

{\mathcal{L}}=\left\{\kappa+\mathrm{i}\eta:\eta=(400K_{j})^{-1/4},\kappa=% \mathbb{Z}\eta,\sqrt{K_{j}+\eta^{2}}\leq\kappa\leq\sqrt{K_{j+1}+\eta^{2}}+2T,j% \in\mathbb{N}\right\}.

Then for any $\sqrt{w}=\kappa+\mathrm{i}\eta\in{\mathcal{L}}$ , it holds $\kappa\eta^{2}\geq 1/4$ , and ${\rm Im}[\sqrt{w}]=\eta\geq 1/{\kappa\eta^{1-\delta}}$ , provided $K$ is large enough.

Now 3.6 follows from the following estimate on one point.

Lemma 3.7.

For any $B>0$ , the following holds true for small enough $\delta,c>0$ . Take any $j\geq 0$ , $u=\kappa_{0}+\eta\mathrm{i}\in{\mathcal{L}}$ with $\eta=(400K_{j})^{-1/4}$ , and let $\sqrt{w_{t}}=u-t/2=:\kappa_{t}+\eta\mathrm{i}$ for each $t$ . Conditional on $Y_{0}$ with $|Y_{0}(w_{0})-\sqrt{w_{0}}|\leq B|w_{0}|^{-1/2}$ , with probability $\geq 1-e^{-c\sqrt{K_{j}}}$ ,

(3.27)

\displaystyle|Y_{t}(w_{t})-\sqrt{w_{t}}|\leq\frac{1}{\kappa_{t}\eta^{1-\delta}% },\quad\forall\;0\leq t\leq\sigma_{0}\wedge(2\kappa_{0}-2\sqrt{K_{j}+\eta^{2}}).

Proof of 3.6.

As $Y_{0}$ is $({\mathfrak{d}},C_{*})$ -Airy-like, 2.5 implies that (for some $B>0$ )

(3.28)

\displaystyle|Y_{0}(w)-\sqrt{w}|\leq B|w|^{-1/2}\text{ for all }\arg(w)\in(0,3% \pi/4),|w|>B.

By an union bound over all the points in ${\mathcal{L}}$ , 3.7 implies that conditional on (3.28), with probability at least (for some $c^{\prime}>0$ , and taking $c<c^{\prime}$ and $T$ large)

(3.29)

\displaystyle 1-\sum_{j=1}^{\infty}3T(400K_{j})^{1/4}e^{-c^{\prime}\sqrt{K_{j}% }}=1-\sum_{j=1}^{\infty}3T(400jK)^{1/4}e^{-c^{\prime}\sqrt{jK}}\geq 1-e^{-c% \sqrt{K}},

it holds that for any $j\in\mathbb{N}$ and $\sqrt{w_{0}}=\kappa_{0}+(400K_{j})^{-1/4}\mathrm{i}\in{\mathcal{L}}$ ,

(3.30)

\displaystyle|Y_{t}(w_{t})-\sqrt{w_{t}}|\leq\frac{1}{\kappa_{t}\eta^{1-\delta}% },\quad\forall\;0\leq t\leq\sigma_{0}\wedge(2\kappa_{0}-2\sqrt{K_{j}+\eta^{2}}).

Next we show that (3.30) implies that $\sigma_{0}=2T$ . Take any $t\in[0,2T]$ , and $x\geq K$ . Then $K_{j}\leq x<K_{j+1}$ for some $j\in\mathbb{N}$ . Take $\eta=(400K_{j})^{-1/4}$ . Then $\sqrt{x+\eta^{2}}+t/2\in[i\eta,(i+1)\eta]$ for some $i\in\mathbb{N}$ . We can take $\kappa_{0}=(i+1)\eta$ and $\sqrt{w_{0}}=\kappa_{0}+\mathrm{i}\eta\in{\mathcal{L}}$ , so that $\kappa_{t}\in[\sqrt{x+\eta^{2}},\sqrt{x+\eta^{2}}+\eta]$ . Then we have

	$\displaystyle{\rm Re}[w_{t}]$	$\displaystyle=\kappa_{t}^{2}-\eta^{2}\in[x,x+3\sqrt{x}\eta],$
	$\displaystyle{\rm Im}[w_{t}]$	$\displaystyle=2\kappa_{t}\eta\geq 2\sqrt{x}\eta.$

If there is a pole of $Y_{t}$ at $x$ , (by Nevanlinna representation (2.2)) necessarily,

{\rm Im}[Y_{t}(w_{t})]\geq\frac{{\rm Im}[w_{t}]}{|x-w_{t}|^{2}}\geq\frac{{\rm Im% }[w_{t}]}{(3\sqrt{x}\eta)^{2}+{\rm Im}[w_{t}]^{2}}=\frac{2\kappa_{t}}{9x\eta+4% \kappa_{t}^{2}\eta^{2}}\geq\frac{1}{7\sqrt{x}\eta},

However, (3.30) (note that $\kappa_{t}\geq\sqrt{K_{j}+\eta^{2}}$ , so $t\leq 2\kappa_{0}-2\sqrt{K_{j}+\eta^{2}}$ ) implies

{\rm Im}[Y_{t}(w_{t})]\leq{\rm Im}[\sqrt{w_{t}}]+\frac{1}{\kappa_{t}\eta^{1-% \delta}}=\eta+\frac{1}{\kappa_{t}\eta^{1-\delta}}\leq\eta+\frac{1}{\sqrt{x}% \eta^{1-\delta}}<\frac{1}{7\sqrt{x}\eta}

This leads to a contradiction. Therefore (with probability $\geq 1-e^{-c\sqrt{K}}$ ) there is no pole in $(K,\infty)$ at any time in $[0,\sigma_{0}]$ . Recall the definition of $\sigma_{0}$ from (3.26). Using that $Y_{t}$ is continuous in $t$ , and 2.7 and 2.8, we conclude that $\sigma_{0}=2T$ , and 3.6 follows. ∎

Proof of 3.7.

We introduce the following stopping time

(3.31)

\displaystyle\sigma=\inf\left\{t\leq\sigma_{0}:|Y_{t}(w_{t})-\sqrt{w_{t}}|\geq% \frac{1}{\kappa_{t}\eta^{1-\delta}},\text{ or }\kappa_{t}^{2}-\eta^{2}\leq K_{% j},\text{ or }t=\sigma_{0}\right\}.

For $t\leq\sigma$ , it is necessary that $t\leq\sigma_{0}$ , and $Y_{t}$ has no pole in $(K,\infty)$ . Moreover, we also have that $\kappa_{t}\geq\sqrt{K_{j}+\eta^{2}}\geq\sqrt{K_{j}}$ , and $\kappa_{t}\eta^{2-\delta}\geq 1$ , provided that $K$ is large enough.

We now consider the terms in the RHS of (3.6). For $0\leq t\leq\sigma$ , using Nevanlinna representation (2.2) we have

|\operatorname{\partial}_{w}^{2}Y_{t}(w_{t})|\leq\sum_{x\in P}\frac{1}{|x-w_{t% }|^{3}}\leq\frac{1}{|w_{t}-K|}\sum_{x\in P}\frac{1}{|x-w_{t}|^{2}}\leq\frac{{% \rm Im}[Y_{t}(w_{t})]}{|w_{t}-K|{\rm Im}[w_{t}]},

where the second inequality follows from that ${\rm Re}[w_{t}]=\kappa_{t}^{2}-\eta^{2}\geq K_{j}\geq K\geq x$ for any $x\in P$ . As $|{\rm Im}[Y_{t}(w_{t})]-{\rm Im}[\sqrt{w_{t}}]|\leq 1/(\kappa_{t}\eta^{1-% \delta})\leq\eta$ , we have ${\rm Im}[Y_{t}(w_{t})]\leq 2\eta$ , so

\displaystyle|\operatorname{\partial}_{w}^{2}Y_{t}(w_{t})|\leq\frac{2\eta}{|w_% {t}-K|{\rm Im}[w_{t}]}\leq\frac{\sqrt{2}}{(\kappa_{t}^{2}-\eta^{2}-K+2\kappa_{% t}\eta)\kappa_{t}},

where we used ${\rm Im}[w_{t}]=2\kappa_{t}\eta$ , and $|w_{t}-K|=|\kappa_{t}^{2}-\eta^{2}-K+2\mathrm{i}\kappa_{t}\eta|$ for the second inequality. By integrating in time, we have

(3.32)

\displaystyle\begin{split}&\int_{0}^{t}|\operatorname{\partial}_{w}^{2}Y_{s}(w% _{s})|{\rm d}s\leq\int_{0}^{t}\frac{\sqrt{2}}{(\kappa_{t}\kappa_{s}-\eta^{2}-K% +2\kappa_{t}\eta)\kappa_{t}}{\rm d}s\\ \leq&\frac{2\sqrt{2}(\log(\kappa_{t}\kappa_{0}-\eta^{2}-K+2\kappa_{t}\eta)-% \log(\kappa_{t}^{2}-\eta^{2}-K+2\kappa_{t}\eta))}{\kappa_{t}^{2}}\leq\frac{2% \sqrt{2}\log(1+\kappa_{0}/(2\eta))}{\kappa_{t}^{2}}.\end{split}

The quadratic variation of the martingale term is given by (3.7). By using Nevanlinna representation (2.2) for $Y_{t}$ , we have

\left|\frac{{\rm d}}{{\rm d}t}\left\langle\int_{0}^{t}({\rm d}M_{s})(w_{s})% \right\rangle\right|\leq\sum_{x\in P}\frac{2}{\beta|x-w_{t}|^{4}}.

Then using ${\rm Re}(w_{t})\geq K_{j}\geq K\geq x$ for any $x\in P$ , we can bound the above by

\frac{2}{\beta|w_{t}-K|^{2}}\sum_{x\in P}\frac{1}{|x-w_{t}|^{2}}\leq\frac{2{% \rm Im}[Y_{t}(w_{t})]}{\beta|w_{t}-K|^{2}{\rm Im}[w_{t}]}\leq\frac{4}{\beta(% \kappa_{t}^{2}-\eta^{2}-K+2\kappa_{t}\eta)^{2}\kappa_{t}},

where we used ${\rm Im}[Y_{t}(w_{t})]\leq 2\eta$ , ${\rm Im}[w_{t}]=2\kappa_{t}\eta$ , and $|w_{t}-K|=|\kappa_{t}^{2}-\eta^{2}-K+2\mathrm{i}\kappa_{t}\eta|$ for the last inequality. By integrating in time, we have

\int_{0}^{t}\frac{4{\rm d}s}{\beta(\kappa_{s}^{2}-\eta^{2}-K+2\kappa_{s}\eta)^% {2}\kappa_{s}}\leq\int_{0}^{t}\frac{4{\rm d}s}{\beta(\kappa_{t}\kappa_{s}-\eta% ^{2}-K+2\kappa_{t}\eta)^{2}\kappa_{t}}\\ \leq\frac{8}{\beta(\kappa_{t}^{2}-\eta^{2}-K+2\kappa_{t}\eta)\kappa_{t}^{2}}% \leq\frac{4}{\beta\kappa_{t}^{3}\eta}.

Similarly to (3.21), taking $\theta=K_{j}^{1/4}$ , a union bound implies that with probability $1-\log(T)e^{-c\theta}$ ,

(3.33)

\left|\int_{0}^{t}({\rm d}M_{s})(w_{s})\right|\leq\frac{\theta}{\sqrt{\kappa_{% t}^{3}\eta}},\quad\forall\;0\leq t\leq\sigma.

Now by (3.6), and using (LABEL:eq:lasder) and (3.33), we have that for any $0\leq t\leq\sigma$ ,

(3.34)

\displaystyle\left|Y_{t}(w_{t})-\sqrt{w_{t}})\right|\leq\int_{0}^{t}\left|Y_{s% }(w_{s})-\sqrt{w_{s}}\right|\left|\partial_{w}Y_{s}(w_{s})\right|{\rm d}s+% \frac{2\theta}{\sqrt{\kappa_{t}^{3}\eta}}+\left|Y_{0}(w_{0})-\sqrt{w_{0}}% \right|.

Here we used that $\frac{2\sqrt{2}\log(1+\kappa_{0}/(2\eta))}{\kappa_{t}^{2}}<\frac{\theta}{\sqrt% {\kappa_{t}^{3}\eta}}$ , provided $K$ is large enough.

As for $\operatorname{\partial}_{w}Y_{t}(w_{t})$ , for $0\leq t\leq\sigma$ (using 2.2) we have

\displaystyle|\operatorname{\partial}_{w}Y_{t}(w_{t})|\leq\frac{{\rm Im}[Y_{t}% (w_{t})]}{{\rm Im}[w_{t}]}\leq\frac{1}{2\kappa_{t}}+\frac{|Y_{t}(w_{t})-\sqrt{% w_{t}}|}{2\kappa_{t}\eta}\leq\frac{1}{2\kappa_{t}}+\frac{1}{2\kappa_{t}^{2}% \eta^{2-\delta}}=:\gamma(t)\leq\frac{1}{\kappa_{t}}.

Then for any $0\leq s<t\leq\sigma$ (noticing that $\kappa_{t}\eta^{2}\geq 1/20$ ),

(3.35)

\displaystyle\int_{s}^{t}\gamma(u){\rm d}u\leq\log(\kappa_{s}/\kappa_{t})+% \frac{1}{\kappa_{t}\eta^{2-\delta}}\leq\log(\kappa_{s}/\kappa_{t})+20\eta^{% \delta}.

By Grönwall’s inequality and (3.34), for any $0\leq t\leq\sigma$ , we can bound $|Y_{t}(w_{t})-\sqrt{w_{t}}|$ by:

\frac{2\theta}{\sqrt{\kappa_{t}^{3}\eta}}+\left|Y_{0}(w_{0})-\sqrt{w_{0}}% \right|+\int_{0}^{t}\gamma(s)\left(\frac{2\theta}{\sqrt{\kappa_{s}^{3}\eta}}+% \left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|\right)\exp\left(\int_{s}^{t}\gamma(u){% \rm d}u\right){\rm d}s.

By (3.35), $\exp(\int_{s}^{t}\gamma(u){\rm d}u)\leq\frac{2\kappa_{s}}{\kappa_{t}}$ . Also

\int_{0}^{t}\gamma(s)\frac{\theta}{\sqrt{\kappa_{s}^{3}\eta}}\cdot\frac{\kappa% _{s}}{\kappa_{t}}{\rm d}s\leq\frac{\theta}{\kappa_{t}}\int_{0}^{t}\frac{1}{% \sqrt{\kappa_{s}^{3}\eta}}{\rm d}s\leq\frac{\theta}{\sqrt{\kappa_{t}^{3}\eta}},

and

\int_{0}^{t}\gamma(s)\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|\frac{\kappa_{s}}{% \kappa_{t}}{\rm d}s\leq\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|\int_{0}^{t}\frac% {1}{\kappa_{t}}{\rm d}s=\frac{t\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|}{\kappa_% {t}}.

Therefore for $\theta=K_{j}^{1/4}$ , it holds

\left|Y_{t}(w_{t})-\sqrt{w_{t}}\right|\leq\frac{6\theta}{\sqrt{\kappa_{t}^{3}% \eta}}+\left(1+\frac{2t}{\kappa_{t}}\right)\left|Y_{0}(w_{0})-\sqrt{w_{0}}% \right|\leq\frac{1}{2\kappa_{t}\eta^{1-\delta}},

where in the last inequality we used that

\left(1+\frac{2t}{\kappa_{t}}\right)\left|Y_{0}(w_{0})-\sqrt{w_{0}}\right|\leq B% \left(1+\frac{2t}{\kappa_{t}}\right)|w_{0}|^{-1/2}\leq B\left(1+\frac{2t}{% \kappa_{t}}\right)\frac{1}{\kappa_{t}+t/2}\leq\frac{1}{3\kappa_{t}\eta^{1-% \delta}}.

Therefore, since $Y_{t}(w_{t})-\sqrt{w_{t}}$ is continuous in $t$ , we conclude that $\sigma=\sigma_{0}\wedge(2\kappa_{0}-2\sqrt{K_{j}+\eta^{2}})$ with the desired probability. ∎

3.4. Proof of 3.1

As before, take $T$ large enough and $K=(\log T)^{8}$ . 3.6 and 3.3 verify the two assumptions in 2.3, respectively. Thus with probability $1-e^{-c(\log T)^{2}}$ , for any $t\in[T,2T]$ , $Y_{t}$ has infinitely many poles $x_{1}(t)\geq x_{2}(t)\geq\cdots$ , satisfying

(3.36)

\displaystyle|x_{i}(t)-{\mathfrak{a}}_{i}|<CK^{4}i^{-\delta/6}.

This gives the second statement in 3.1, after replacing $\delta/6$ by $\delta$ .

The statement (3.36) also verifies the first assumption in 2.1. Moreover, in (3.9) we can take a sequence of complex numbers $w_{n}=n\mathrm{i}$ , so that $|Y_{t}(w_{n})-\sqrt{w_{n}}|\to 0$ as $n\to\infty$ . This verifies the second assumption in 2.1, from which (3.1) holds.

Finally, for each $i\in\mathbb{N}$ , the continuity of $x_{i}(t)$ in $t\in[T,2T]$ follows from the continuity of $Y_{t}$ in $t$ , and 2.7, 2.8. ∎

4. Hölder Regularity

In this section, we upgrade the trajectory continuity into Hölder regularity.

Proposition 4.1.

For any $B>0$ , there exist large and small $C,c>0$ such that the following holds. Take any $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.5, and that its poles are given by a line ensemble $\{x_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ , and that $Y_{0}(w)=\sum_{i=1}^{\infty}\frac{1}{x_{i}(0)-w}-\frac{1}{{\mathfrak{a}}_{i}}-% \frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}$ for any $w\in\mathbb{H}$ . Take $k\in\mathbb{N}$ large enough and any $0<\xi<ck^{-4/3}$ . Then conditional on the event that

(4.1)

|x_{i}(0)-{\mathfrak{a}}_{i}|\leq B,\quad\forall i\in\llbracket k/2,2k\rrbracket,

with probability $>1-e^{-ck^{1/6}}$ we have

\max_{0\leq s\leq\xi}|x_{k}(s)-x_{k}(0)|\leq Ck^{2/3}\xi^{1/2}.

Proof.

We first prove the (with conditional probability $>1-e^{-ck^{1/6}}$ ) upper bound

(4.2)

\displaystyle\max_{0\leq s\leq\xi}x_{k}(s)\leq x_{k}(0)+Ck^{2/3}\xi^{1/2}.

The lower bound $\min_{0\leq s\leq\xi}x_{k}(s)\geq x_{k}(0)-Ck^{2/3}\xi^{1/2}$ can be proven in the same way.

By (4.1) with (2.5), we conclude that there exists a large constant $C_{1}$ such that

(4.3)

\displaystyle|\{i:x_{i}(0)\in[x_{k}(0)-1,x_{k}(0)+1]\}|\leq C_{1}k^{1/3}.

Take small $\delta\leq 1/(C_{1}k^{1/3})$ . Then (4.3) implies that there exists some $1\leq\ell\leq C_{1}k^{1/3}$ , such that $x_{k-\ell-1}(0)-x_{k-\ell}(0)\geq\delta$ . We take the smallest such $\ell$ , then $|x_{k-\ell}(0)-x_{k}(0)|\leq\delta\ell$ .

Let $E=x_{k-\ell}(0)+\delta/2$ , and take a small $b>0$ and $w=E+\mathrm{i}b$ . Next we show that, provided $\delta\geq 8bC_{1}^{1/2}k^{1/6}$ ,

(4.4)

\displaystyle{\rm Im}[Y_{0}(w)]={\rm Im}[Y_{0}(E+\mathrm{i}b)]=\sum_{i=1}^{% \infty}\frac{b}{|E+\mathrm{i}b-x_{i}(0)|^{2}}\leq\frac{1}{4b}.

For this, from (4.1) and (2.5), we have $|\{i:x_{i}(0)\in[E-\delta/2-m,E-\delta/2-(m-1)]\}|\leq C_{1}m^{1/2}k^{1/3}$ for each $m\in\mathbb{N}$ . It follows that

	$\displaystyle\sum_{i:x_{i}(0)\leq E-\delta/2}\frac{b}{\|E+\mathrm{i}b-x_{i}(0)\|% ^{2}}$	$\displaystyle\leq\sum_{m=1}^{\infty}\sum_{i:x_{i}(0)\in[E-\delta/2-m,E-\delta/% 2-(m-1)]}\frac{b}{\|E+\mathrm{i}b-x_{i}(0)\|^{2}}$
		$\displaystyle\leq\sum_{m=1}^{\infty}\frac{C_{1}bm^{1/2}k^{1/3}}{(\delta/2+m-1)% ^{2}}\leq\frac{8C_{1}bk^{1/3}}{\delta^{2}}\leq\frac{1}{8b},$

using that $\delta\geq 8bC_{1}^{1/2}k^{1/6}$ . By a similar argument, we can also upper bound the summation over $x_{i}(0)\geq E+\delta/2$ , and (4.4) follows.

We now introduce a stopping time $\tau$ :

\tau=\inf\left\{s\geq 0:{\rm Im}[Y_{s}(w)]\geq\frac{1}{2b}\right\}\wedge\xi.

For $0\leq s\leq\tau$ , it follows from 2.2

|\partial_{w}^{2}Y_{s}(w)|\leq\frac{2{\rm Im}[Y_{s}(w)]}{b^{2}}\leq\frac{1}{b^% {3}},\quad|\operatorname{\partial}_{w}Y_{s}(w)^{2}|=2|\partial_{w}Y_{s}(w)|% \cdot|Y_{s}(w)|\leq\frac{|Y_{s}(w)|}{b^{2}}.

Thus by 1.5,

(4.5)

\displaystyle|Y_{s}(w)-Y_{0}(w)|\leq\left|\int_{0}^{s}{\rm d}M_{u}(w)\right|+% \frac{\int_{0}^{s}|Y_{u}(w)|{\rm d}u}{2b^{2}}+\operatorname{\mathcal{O}}\left(% \frac{s}{b^{3}}\right).

For the Martingale term, using (1.10) and 2.2, it follows that (for $0\leq s\leq\tau$ )

\left|\frac{{\rm d}}{{\rm d}s}\langle M_{s}(w)\rangle\right|\leq\frac{2{\rm Im% }[Y_{s}(w)]}{\beta b^{3}}\leq\frac{1}{\beta b^{4}}.

Therefore, by the Burkholder-Davis-Gundy inequality, there exists a small constant $c_{1}>0$ ,

(4.6)

\displaystyle\mathbb{P}\left(\sup_{0\leq s\leq\tau}\left|\int_{0}^{s}{\rm d}M_% {u}(w)\right|\geq\frac{k^{1/6}\xi^{1/2}}{b^{2}}\right)\leq e^{-c_{1}k^{1/6}}

By plugging (4.6) into (4.5), it follows that with probability $e^{-c_{1}k^{1/6}}$ ,

(4.7)

\displaystyle\left|Y_{s}(w)-Y_{0}(w)\right|=\operatorname{\mathcal{O}}\left(% \frac{k^{1/6}\xi^{1/2}}{b^{2}}+\frac{\xi}{b^{3}}\right),\quad\forall s\in[0,% \tau],

provided that $\xi\leq b^{2}$ . Now we take a large $C_{2}>0$ , and assume that $\xi\leq b^{2}/(C_{2}k^{1/3})$ . Then for any $s\in[0,\tau]$ , and using (4.4), we have

{\rm Im}[Y_{s}(w)]\leq|Y_{s}(w)-Y_{0}(w)|+{\rm Im}[Y_{0}(w)]<\frac{1}{2b},

which, in particular, implies that $\tau=\xi$ . Since ${\rm Im}[Y_{s}(w)]=\sum_{i=1}^{\infty}\frac{b}{|w-x_{i}(s)|^{2}}$ , this further implies that $\{x_{i}(s)\}_{i\in\mathbb{N}}\cap[E-b,E+b]=\emptyset$ , for any $s\in[0,\xi]$ .

In summary, we conclude that with probability $1-e^{-c_{1}k^{1/6}}$ , $\{x_{i}(s)\}_{i\in\mathbb{N}}\cap[E-b,E+b]=\emptyset$ , for any $s\in[0,\xi]$ . Since $x_{k}(0)<E$ , it follows that (with probability $1-e^{-c_{1}k^{1/6}}$ ) for any $s\in[0,\xi]$ ,

x_{k}(s)\leq E-b=x_{k-\ell}(0)+\delta/2-b\leq x_{k}(0)+\delta\ell+\delta/2\leq x% _{k}(0)+2C_{1}k^{1/3}\delta.

Finally, we choose the parameter $\delta$ and $b$ , satisfying all the above constraints:

\delta\leq 1/(C_{1}k^{1/3}),\quad\delta\geq 8bC_{1}^{1/2}k^{1/6},\quad\xi\leq b% ^{2}/(C_{2}k^{1/3}).

Then we can take $b=C_{2}^{1/2}k^{1/6}\xi^{1/2}$ and $\delta=8C_{1}^{1/2}C_{2}^{1/2}k^{1/3}\xi^{1/2}$ . By taking $C$ large enough and $c$ small enough (depending on $C_{1}$ , $C_{2}$ , and $c_{1}$ ), the constraint $\delta\leq 1/(C_{1}k^{1/3})$ is also satisfied since $\xi<ck^{-4/3}$ ; and (4.2) follows. ∎

5. Recover Dyson Brownian Motion

In this section, we localize any line ensemble given by the poles of the SDE (1.9), by deriving another SDE satisfied by the first finitely many poles (in the sense of weak solution).

More precisely, for the line ensemble $\{x_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ of poles, we prove that, if for some large $k\in\mathbb{N}$ , $x_{k}(t)$ and $x_{k+1}(t)$ are bounded away from each other for certain amount of time, the evolution of $x_{1}(t)\geq x_{2}(t)\geq\cdots\geq x_{k}(t)$ is then described by DBM plus a drift term, describing the effect of $\{x_{i}(t)\}_{i=k+1}^{\infty}$ .

Proposition 5.1.

For any $C>0$ the following is true. Take any $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.5, and that its poles are given by a line ensemble $\{x_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R}}}$ . Fix a large $k\in\mathbb{N}$ , and denote the stopping time

(5.1)

\displaystyle\tau=\inf\{t\geq 0:x_{k}(t)-x_{k+1}(t)\leq 1/(Ck^{1/3})\}\cup\{1\}.

Conditional on the event $\tau>0$ , there exist independent Brownian motions $B_{1},\cdots,B_{k}$ adapted to the filtration ${\mathcal{F}}_{t}=\sigma(\{{\bm{x}}(u)\}_{u\leq t})$ , satisfying

(5.2)

{\rm d}x_{i}(t)=\sqrt{\frac{2}{\beta}}{\rm d}B_{i}(t)+\sum_{1\leq j\leq k\atop j% \neq i}\frac{{\rm d}t}{x_{i}(t)-x_{j}(t)}+W_{t}(x_{i}(t)){\rm d}t,\quad\forall i% \in\llbracket 1,k\rrbracket,\;t\in[0,\tau],

where $W_{t}$ is a random meromorphic function, defined as

W_{t}(w)=-Y_{t}(w)+\sum_{i=1}^{k}\frac{1}{x_{i}(t)-w}.

Moreover, almost surely the following holds:

(5.3)

\displaystyle\int_{0}^{\tau}\mathds{1}(\exists 1\leq i<j\leq k:x_{i}(t)=x_{j}(% t)){\rm d}t=0.

The rest of this section is devoted to proving 5.1. Using that $\{x_{i}(t)\}_{i\in\mathbb{N}}$ are poles of $Y_{t}$ , we derive (5.2) from 1.5, using a contour integral. For this, we need first establish that the poles do not collide at almost every time (i.e., (5.3)). The idea to establish the collision time estimate is to consider the process $(x_{i}(t)-x_{j}(t))^{2}$ for some $i<j$ , showing that its level- $0$ local time equals $0$ . We mainly follow the standard argument used to study the Bessel process, see [114, Chapter XI, Section 1]. To analyze such processes we again resort to contour integrals, and therefore an induction will be used.

We next give the semi-martingale decomposition of a process, which is the sum of $(x_{i}(t)-x_{j}(t))^{2}$ for $i,j$ in an interval.

For any line ensemble $\{{\bm{x}}(t)\}_{t\in{\mathbb{R}}}$ satisfying 1.5, and any $a,\alpha\in\mathbb{N}$ with $\alpha\geq 2$ , denote

(5.4)

W^{a,\alpha}_{t}(w)=-Y_{t}(w)+\sum_{i=a}^{a+\alpha-1}\frac{1}{x_{i}(t)-w},

and

Z^{a,\alpha}(t)=\sum_{a\leq i<j\leq a+\alpha-1}(x_{i}(t)-x_{j}(t))^{2}=\alpha% \sum_{i=1}^{a+\alpha-1}x_{i}(t)^{2}-\left(\sum_{i=1}^{a+\alpha-1}x_{i}(t)% \right)^{2}.

Lemma 5.2.

In the above setup, take any $t_{0}\in{\mathbb{R}}$ , and denote the stopping time

\sigma=\inf\{t\geq t_{0}:x_{a-1}(t)=x_{a}(t)\text{ or }x_{a+\alpha-1}(t)=x_{a+% \alpha}(t)\text{ or }t=t_{0}+1\},

where we use the convention that $x_{0}(t)=\infty$ , when $a=1$ . Then $Z^{a,\alpha}$ in $[t_{0},\sigma)$ satisfies ${\rm d}Z^{a,\alpha}(t)={\rm d}M^{a,\alpha}(t)+V^{a,\alpha}(t){\rm d}t$ , where

V^{a,\alpha}(t)=2\sum_{1\leq i<j\leq\alpha}(x_{i}(t)-x_{j}(t))(W^{a,\alpha}_{t% }(x_{i})-W^{a,\alpha}_{t}(x_{j}))+\alpha^{2}(\alpha-2)+\frac{2\alpha(\alpha-1)% }{\beta},

and ${\rm d}M^{a,\alpha}(t)$ is the Martingale term, with quadratic variation

(5.5)

\frac{{\rm d}}{{\rm d}t}\langle M^{a,\alpha}(t)\rangle=\frac{8\alpha}{\beta}Z^% {a,\alpha}(t).

Proof.

For simplicity of notations, we fix $a,\alpha$ , and write $Z(t)=Z(t)^{a,\alpha}$ , $V(t)=V(t)^{a,\alpha}$ , an ${\rm d}M(t)={\rm d}M(t)^{a,\alpha}$ within this proof.

For $t$ with $x_{a-1}(t)>x_{a}(t)$ and $x_{a+\alpha-1}(t)>x_{a+\alpha}(t)$ , we take a contour ${\mathcal{C}}={\mathcal{C}}_{t}$ enclosing $x_{a}(t),\cdots,x_{a+\alpha-1}(t)$ , but not any $x_{i}(t)$ for $i<a$ or $i\geq a+\alpha$ . Then by 1.5, we have

(5.6)

\displaystyle\begin{split}{\rm d}\sum_{i=a}^{a+\alpha-1}x_{i}(t)&=-\frac{1}{2% \pi\mathrm{i}}\oint_{\mathcal{C}}w{\rm d}Y_{t}(w){\rm d}w\\ &=-\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}w{\rm d}w\left({\rm d}M_{t}(w)+% \left(\frac{2-\beta}{2\beta}\operatorname{\partial}^{2}_{w}Y_{t}(w)+\frac{1}{2% }\operatorname{\partial}_{w}Y_{t}(w)^{2}-\frac{1}{2}\right){\rm d}t\right).% \end{split}

Note that

\oint_{\mathcal{C}}w{\rm d}w=0,\quad\oint_{\mathcal{C}}w\partial_{w}^{2}Y_{t}(% w){\rm d}w=\oint_{\mathcal{C}}\sum_{i=1}^{\infty}\frac{2w{\rm d}w}{(x_{i}(t)-w% )^{3}}=0,

and

-\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}\frac{w}{2}\partial_{w}Y_{t}(w)^{2% }{\rm d}w=\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}\frac{Y_{t}(w)^{2}}{2}{% \rm d}w=\sum_{i=a}^{a+\alpha-1}W_{t}^{i,1}(x_{i}(t)).

Note that the poles of $W_{t}^{i,1}$ are $x_{1}(t),x_{2}(t),\ldots$ , except for $x_{i}(t)$ . Then we have

(5.7)

{\rm d}\sum_{i=a}^{a+\alpha-1}x_{i}(t)=-\frac{1}{2\pi\mathrm{i}}\oint_{% \mathcal{C}}w{\rm d}M_{t}(w){\rm d}w+\sum_{i=a}^{a+\alpha-1}W^{a,\alpha}_{t}(x% _{i}(t)){\rm d}t,

using that $\sum_{i=a}^{a+\alpha-1}W_{t}^{i,1}(x_{i}(t))=\sum_{i=a}^{a+\alpha-1}W^{a,% \alpha}_{t}(x_{i}(t))$ .

Similarly, we have

(5.8)

\displaystyle\begin{split}{\rm d}\sum_{i=a}^{a+\alpha-1}x^{2}_{i}(t)&=-\frac{1% }{2\pi\mathrm{i}}\oint_{\mathcal{C}}w^{2}{\rm d}Y_{t}(w){\rm d}w\\ &=-\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}w^{2}{\rm d}w\left({\rm d}M_{t}(% w)+\left(\frac{2-\beta}{2\beta}\operatorname{\partial}^{2}_{w}Y_{t}(w)+\frac{1% }{2}\operatorname{\partial}_{w}Y_{t}(w)^{2}-\frac{1}{2}\right){\rm d}t\right)% \\ &=-\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}w^{2}{\rm d}M_{t}(w){\rm d}w+% \left(2\sum_{i=a}^{a+\alpha-1}W_{t}^{a,\alpha}(x_{i}(t))x_{i}(t)+\alpha(\alpha% -2)+\frac{2\alpha}{\beta}\right){\rm d}t.\end{split}

Here for the last equality, we used that

\oint_{\mathcal{C}}w^{2}{\rm d}w=0,\quad-\frac{1}{2\pi\mathrm{i}}\oint_{% \mathcal{C}}w^{2}\partial_{w}^{2}Y_{t}(w){\rm d}w=-\frac{1}{2\pi\mathrm{i}}% \oint_{\mathcal{C}}\sum_{i=1}^{\infty}\frac{2w^{2}{\rm d}w}{(x_{i}(t)-w)^{3}}=% 2\alpha,

-\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}\frac{w^{2}}{2}\partial_{w}Y_{t}(w% )^{2}{\rm d}w=\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}wY_{t}(w)^{2}{\rm d}w% =2\sum_{i=a}^{a+\alpha-1}W_{t}^{i,1}(x_{i}(t))x_{i}(t),

and that $\sum_{i=a}^{a+\alpha-1}W_{t}^{i,1}(x_{i}(t))x_{i}(t)=\sum_{i=a}^{a+\alpha-1}W^% {a,\alpha}_{t}(x_{i}(t))x_{i}(t)+\frac{\alpha(\alpha-1)}{2}$ .

Now using (5.6) and (5.8), and Ito’s formula, we can write ${\rm d}Z(t)={\rm d}M(t)+V(t){\rm d}t$ , with

\displaystyle\begin{split}V(t)=&2\alpha\sum_{i=a}^{a+\alpha-1}W_{t}^{a,\alpha}% (x_{i}(t))x_{i}(t)+\alpha^{2}(\alpha-2)+\frac{2\alpha^{2}}{\beta}\\ &-2\left(\sum_{i=a}^{a+\alpha-1}x_{i}(t)\right)\left(\sum_{i=a}^{a+\alpha-1}W_% {t}^{i,1}(x_{i}(t))\right)-\frac{2}{(2\pi\mathrm{i})^{2}\beta}\oiint_{{% \mathcal{C}}^{2}}ww^{\prime}\operatorname{\partial}_{w}\operatorname{\partial}% _{w^{\prime}}\frac{Y_{t}(w)-Y_{t}(w^{\prime})}{w-w^{\prime}}{\rm d}w{\rm d}w^{% \prime}\\ =&2\sum_{1\leq i<j\leq\alpha}(x_{i}(t)-x_{j}(t))(W^{a,\alpha}_{t}(x_{i})-W^{a,% \alpha}_{t}(x_{j}))+\alpha^{2}(\alpha-2)+\frac{2\alpha(\alpha-1)}{\beta}.\end{split}

Here we used (1.11) in the first equality. For the second equality, it is by evaluating the contour integral in $w$ and $w^{\prime}$ , via integration by parts.

As for ${\rm d}M(t)$ , we have

{\rm d}M(t)=\frac{1}{2\pi\mathrm{i}}\oint_{\mathcal{C}}\left(2w\sum_{i=a}^{a+% \alpha-1}x_{i}(t)-\alpha w^{2}\right){\rm d}M_{t}(w){\rm d}w.

By (1.11), the quadratic variation ${\rm d}\langle M(t)\rangle/{\rm d}t$ therefore equals

		$\displaystyle\frac{2}{(2\pi\mathrm{i})^{2}\beta}\oiint_{{\mathcal{C}}^{2}}% \left(2w\sum_{i=a}^{a+\alpha-1}x_{i}(t)-\alpha w^{2}\right)\left(2w^{\prime}% \sum_{i=a}^{a+\alpha-1}x_{i}(t)-\alpha w^{\prime 2}\right)\operatorname{% \partial}_{w}\operatorname{\partial}_{w^{\prime}}\frac{Y_{t}(w)-Y_{t}(w^{% \prime})}{w-w^{\prime}}{\rm d}w{\rm d}w^{\prime}$
	$\displaystyle=$	$\displaystyle\frac{8}{(2\pi\mathrm{i})^{2}\beta}\oiint_{{\mathcal{C}}^{2}}% \left(\sum_{i=a}^{a+\alpha-1}x_{i}(t)-\alpha w\right)\left(\sum_{i=a}^{a+% \alpha-1}x_{i}(t)-\alpha w^{\prime}\right)\frac{Y_{t}(w)-Y_{t}(w^{\prime})}{w-% w^{\prime}}{\rm d}w{\rm d}w^{\prime}.$

By taking the $w^{\prime}$ residues at $x_{a}(t),\ldots,x_{a+\alpha-1}(t)$ , we get

\frac{8}{2\pi\mathrm{i}\beta}\oint_{\mathcal{C}}\left(\sum_{i=a}^{a+\alpha-1}x% _{i}(t)-\alpha w\right)\left(\sum_{j=a}^{a+\alpha-1}\frac{\sum_{i=a}^{a+\alpha% -1}x_{i}(t)-\alpha x_{j}(t)}{w-x_{j}(t)}\right){\rm d}w.

By further taking the $w$ residues at $x_{1}(t),\ldots,x_{\alpha}(t)$ , this equals

\frac{8\alpha}{\beta}\left(\alpha\sum_{i=a}^{a+\alpha-1}x_{i}^{2}(t)-\left(% \sum_{i=a}^{a+\alpha-1}x_{i}(t)\right)^{2}\right)=\frac{8\alpha}{\beta}Z(t),

and the conclusion follows. ∎

We next establish the collision time estimate, for poles whose indices are in an interval.

Lemma 5.3.

Under the same setup as 5.2, almost surely

\int_{t_{0}}^{\sigma}\mathds{1}[x_{a}(t)=x_{a+1}(t)=\cdots=x_{a+\alpha-1}(t)]{% \rm d}t=0.

Proof.

Again, we write $Z(t)=Z(t)^{a,\alpha}$ and $V(t)=V(t)^{a,\alpha}$ in this proof.

We use the local time of $Z(t)$ to analyze its boundary behavior at zero. According to 5.2, $Z(t)$ for $t\in[t_{0},\sigma)$ is a semi-martingale. We let $L_{t}^{h}$ be the level $h$ local time in $[t_{0},\sigma)$ (with $L_{t_{0}}^{h}=0$ for each $h\in{\mathbb{R}}$ ). Then by [114, Chapter VI, Theorem 1.7], almost surely $L_{t}^{h}$ is continuous in $t$ and cadlag (right continuous) in $h$ , and $L_{t}^{0}=L_{t}^{0}-L_{t}^{0-}$ satisfies

(5.9)

L_{\sigma}^{0}=2\int_{t_{0}}^{\sigma}\mathds{1}(Z(t)=0)V(t){\rm d}t=\left(% \alpha^{2}(\alpha-2)+\frac{2\alpha(\alpha-1)}{\beta}\right)\int_{t_{0}}^{% \sigma}\mathds{1}(Z(t)=0){\rm d}t,

where for the second equality, we used that if $Z(t)=0$ , then $x_{a}(t)=\cdots=x_{a+\alpha-1}(t)$ and $V(t)=\alpha^{2}(\alpha-2)+\frac{2\alpha(\alpha-1)}{\beta}$ .

Thanks to the occupation time formula [114, Chapter VI, Corollary 1.6], we have

\displaystyle\int_{0}^{\infty}h^{-1}L_{\sigma}^{h}{\rm d}h=\int_{t_{0}}^{% \sigma}Z(t)^{-1}{\rm d}\langle Z(t)\rangle\leq\frac{8\alpha}{\beta}<\infty,

where we used (5.5) which gives ${\rm d}\langle Z(t)\rangle=8\alpha Z(t){\rm d}t/\beta$ , and $\sigma\leq t_{0}+1$ . Since $L^{h}_{\sigma}$ is right continuous in $h$ , it follows that $L^{0}_{\sigma}=0$ , and hence (5.9) implies that almost surely $\int_{t_{0}}^{\sigma}\mathds{1}(Z(t)=0){\rm d}t=0$ . Thus the conclusion follows. ∎

Proof of 5.1.

We take the following two steps.

Step 1: Non-collision. We will first show (5.3), i.e., ${\rm d}t$ almost everywhere poles do not collide. More precisely, we will prove inductively on $\ell=1,2,3,\cdots,k-1$

(5.10)

\displaystyle\int_{0}^{\tau}\mathds{1}(x_{1}(t),x_{2}(t),\cdots,x_{k}(t)\text{% take at most $\ell$ distinct values}){\rm d}t=0.

The claim of 5.1 follows from the case of $\ell=k-1$ in (5.10).

For the base case where $\ell=1$ , it follows from 5.3 with $a=1$ and $\alpha=k$ , and $t_{0}=0$ . (Note that in this case, we always have $\sigma\geq\tau$ )

We next give the induction step: if (5.10) holds for some $2\leq\ell<k-1$ , then it holds for $\ell+1$ .

Under the induction hypothesis, ${\rm d}t$ almost everywhere, there exist $0<\alpha_{1}<\alpha_{2}<\cdots<\alpha_{\ell}<k$ such that

(5.11)

\displaystyle x_{\alpha_{1}}(t)>x_{\alpha_{1}+1}(t),\quad x_{\alpha_{2}}(t)>x_% {\alpha_{2}+1}(t),\quad\cdots,\quad x_{\alpha_{\ell}}(t)>x_{\alpha_{\ell}+1}(t),

Fix the indices $0<\alpha_{1}<\alpha_{2}<\cdots<\alpha_{\ell}<k$ , the set of time $t\in[0,\tau]$ such that (5.11) holds is a random open set, and we denote it by $I\subset[0,\tau]$ .

Next we show that almost surely, for almost every $t\in I$ , $x_{1}(t),x_{2}(t),\cdots,x_{k}(t)$ take at least $\ell+2$ distinct values. This implies that (5.10) holds for $\ell+1$ . For the convenience of notations, we denote $\alpha_{0}=0$ and $\alpha_{\ell+1}=k$ . Since $\ell\leq k-2$ , there exists some $\nu\in\llbracket 0,\ell\rrbracket$ such that $\alpha_{\nu+1}-\alpha_{\nu}\geq 2$ . We then apply 5.3 with $a=\alpha_{\nu}+1$ and $\alpha=\alpha_{\nu+1}-\alpha_{\nu}$ , and $t_{0}$ taking any rational numbers. We note that the union of all such $[t_{0},\sigma)$ would cover $I$ , therefore

\int_{I}\mathds{1}(x_{a}(t)=\cdots=x_{a+\alpha-1}(t)){\rm d}t=0.

Thus for almost every $t\in I$ , $x_{a}(t),\cdots,x_{a+\alpha-1}(t)$ would take at least two distinct values, so we finish the induction step.

Then by induction principle, we finish the proof of (5.3).

Step 2: Dyson Brownian motion. We next prove (5.2). For that we need to construct the Brownian motions $B_{i}(t)$ in (5.2). By (5.3) and that each $x_{i}(t)$ is continuous, for any $t\in[0,\tau]$ outside a closed measure zero set (i.e., in a countable union of open intervals, whose closure is $[0,\tau]$ ), we have $x_{i}(t)>x_{i+1}(t)$ for each $i\in\llbracket 1,k\rrbracket$ . Then we can take a small contour ${\mathcal{C}}_{i}={\mathcal{C}}_{i,t}$ enclosing $x_{i}(t)$ but not any other poles. From (5.7) in the proof of 5.2, we have

(5.12)

{\rm d}x_{i}(t)=\sqrt{\frac{2}{\beta}}{\rm d}B_{i}(t)+W_{t}^{i,1}(x_{i}(t)){% \rm d}t,

where

{\rm d}B_{i}(t)=-\sqrt{\frac{\beta}{2}}\frac{1}{2\pi\mathrm{i}}\oint_{{% \mathcal{C}}_{i}}w{\rm d}M_{t}(w){\rm d}w.

We can then further extend $B_{i}(t)$ to all of $[0,\tau]$ as a continuous process. The quadratic variations are given by

\frac{{\rm d}}{{\rm d}t}\langle B_{i}(t),B_{j}(t)\rangle=\frac{1}{(2\pi\mathrm% {i})^{2}}\oiint_{{\mathcal{C}}_{i}\times{\mathcal{C}}_{j}}ww^{\prime}% \operatorname{\partial}_{w}\operatorname{\partial}_{w^{\prime}}\frac{Y_{t}(w)-% Y_{t}(w^{\prime})}{w-w^{\prime}}=\frac{1}{(2\pi\mathrm{i})^{2}}\oiint_{{% \mathcal{C}}_{i}\times{\mathcal{C}}_{j}}\frac{Y_{t}(w)-Y_{t}(w^{\prime})}{w-w^% {\prime}},

which equals $\mathds{1}(i=j)$ . Then it follows that $\{B_{i}(t)\}_{i=1}^{k}$ are independent Brownian motions.

Noting that $W_{t}=W_{t}^{1,k}$ , we get (5.2) from (5.12). ∎

Remark 5.4.

Another approach to study DBM developed in [66] is based on applying Ito’s formula to the elementary symmetric functions $\sum_{1\leq j_{1}<j_{2}\cdots<j_{n}}x_{j_{1}}x_{j_{2}}\cdots x_{j_{n}}$ . There a large family of Dyson type interacting particle systems are considered. For $\beta\geq 1$ , they show that if the initial data has some particles at the same location, they will separate instantly. The method there could potentially be adapted and derive non-collision in the above proof as well.

6. Coupling and Uniqueness

In this section we prove the uniqueness part of 1.6.

Take any $\{Y_{t}\}_{t\in{\mathbb{R}}}$ satisfying 1.4 and 1.5. Let $\{{\bm{x}}(t)\}_{t\in{\mathbb{R}}}=\{x_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R% }}}$ be the line ensemble given by its poles (from 3.2). We also take another line ensemble $\{{\bm{y}}(t)\}_{t\in{\mathbb{R}}}=\{y_{i}(t)\}_{i\in\mathbb{N},t\in{\mathbb{R% }}}$ through the same way.

Proposition 6.1.

The two line ensembles $\{{\bm{x}}(t)\}_{t\in{\mathbb{R}}}$ and $\{{\bm{y}}(t)\}_{t\in{\mathbb{R}}}$ have the same law.

Our general strategy is to construct a coupling of the dynamics in $t$ , where ${\bm{x}}(t)$ and ${\bm{y}}(t)$ would get closer as $t$ increases. Then by sending the starting time of the dynamics to $-\infty$ , one concludes that these two line ensembles must equal in law.

The coupling. There are four parameters $\delta,C,T,n$ in the definition of this coupling. Here $\delta,C>0$ are small and large real numbers; $-T$ is among the sequence $t_{1},t_{2},\cdots\to-\infty$ in 1.4, and $T$ is large enough depending on $\delta,C$ ; and we let $n=\lfloor T\rfloor$ . We shall mainly consider the dynamics of the first order $n$ many paths, for $t\in[-T,T]$ .

For each $t\in[-T,T]$ , let ${\mathcal{E}}[t]$ be the event where

|x_{i}(t)-{\mathfrak{a}}_{i}|,|y_{i}(t)-{\mathfrak{a}}_{i}|\leq\frac{C(\log T)% ^{40}}{i^{\delta}},

for each $i\in\mathbb{N}$ .

For each $\ell\in\llbracket 0,2Tn^{3}\rrbracket$ , denote $t_{\ell}=-T+\ell n^{-3}$ . Under ${\mathcal{E}}[t_{\ell}]$ , by (2.5), we let $k_{\ell}$ and $k^{\prime}_{\ell}$ be the smallest numbers in $[\![n,2n-1$ ]], such that

|x_{k_{\ell}}(t_{\ell})-x_{k_{\ell}+1}(t_{\ell})|\geq\frac{2}{C^{\prime}n^{1/3% }},\quad|y_{k_{\ell}^{\prime}}(t_{\ell})-y_{k_{\ell}^{\prime}+1}(t_{\ell})|% \geq\frac{2}{C^{\prime}n^{1/3}},

where $C^{\prime}$ is a large enough universal constant. Note that when $n$ is large enough depending on $C^{\prime}$ , such $k_{\ell}$ and $k^{\prime}_{\ell}$ exist.

We introduce a stopping time $\tau$ (with respect to the filtration ${\mathcal{F}}_{t}=\sigma(\{{\bm{x}}(u)\}_{u\leq t},\{{\bm{y}}(u)\}_{u\leq t})$ ), as follows. If there exists any $t\in[-T,T]$ such that ${\mathcal{E}}[t]$ does not hold, or if there is any $\ell\in\llbracket 0,2Tn^{3}-1\rrbracket$ and $t\in[t_{\ell},t_{\ell+1}]$ , such that

|x_{k_{\ell}}(t)-x_{k_{\ell}+1}(t)|\wedge|y_{k_{\ell}^{\prime}}(t)-y_{k_{\ell}% ^{\prime}+1}(t)|\leq\frac{1}{C^{\prime}n^{1/3}},

we let $\tau$ be the smallest such $t$ . Otherwise, we let $\tau=\infty$ .

Lemma 6.2.

For any $\varepsilon>0$ , there exist $n,\delta,C$ , such that $\mathbb{P}[\tau=\infty]\geq 1-\varepsilon$ .

Proof.

By 3.1, for small enough $\delta$ , large enough $C$ and $n$ , we have $\mathbb{P}\left[\bigcap_{t\in[-T,T]}{\mathcal{E}}[t]\right]\geq 1-\varepsilon/2$ . Then by the Hölder continuity estimate 4.1,

\mathbb{P}[\tau=\infty]\geq 1-\varepsilon/2-2Tn^{4}e^{-cn^{1/6}}\geq 1-\varepsilon,

where $c>0$ is small enough depending on $C$ , and the second inequality is by taking $n$ large. ∎

By 5.1, we can find a family of independent Brownian motions $\{B_{i}\}_{i=1}^{2n}$ , such that for each $\ell\in\llbracket 0,2Tn^{3}-1\rrbracket$ and $t\in[t_{\ell}\wedge\tau,t_{\ell+1}\wedge\tau]$ , $i\in\llbracket 1,k_{\ell}\rrbracket$ ,

(6.1)

\displaystyle\begin{split}{\rm d}x_{i}(t)&=\sqrt{\frac{2}{\beta}}{\rm d}B_{i}(% t)+\left(\sum_{1\leq j\leq k_{\ell}\atop j\neq i}\frac{1}{x_{i}(t)-x_{j}(t)}+W% _{t}(x_{i}(t))\right){\rm d}t,\\ W_{t}(w)&=\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)}+\sum_{i=1% }^{k_{\ell}}\frac{1}{{\mathfrak{a}}_{i}}+\sum_{i=k_{\ell}+1}^{\infty}\frac{1}{% w-x_{i}(t)}+\frac{1}{{\mathfrak{a}}_{i}},\end{split}

where we used 3.2 for the expression of $W_{t}$ . We can similarly find a family of independent Brownian motions $\{\overline{B}_{i}\}_{i=1}^{2n}$ , such that for each $\ell\in\llbracket 0,2Tn^{3}-1\rrbracket$ and $t\in[t_{\ell}\wedge\tau,t_{\ell+1}\wedge\tau]$ , $i\in\llbracket 1,k_{\ell}^{\prime}\rrbracket$ ,

(6.2)

\displaystyle\begin{split}{\rm d}y_{i}(t)&=\sqrt{\frac{2}{\beta}}{\rm d}% \overline{B}_{i}(t)+\left(\sum_{1\leq j\leq k_{\ell}^{\prime}\atop j\neq i}% \frac{1}{y_{i}(t)-y_{j}(t)}+\overline{W}_{t}(y_{i}(t))\right){\rm d}t,\\ \overline{W}_{t}(w)&=\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0)% }+\sum_{i=1}^{k_{\ell}^{\prime}}\frac{1}{{\mathfrak{a}}_{i}}+\sum_{i=k_{\ell}^% {\prime}+1}^{\infty}\frac{1}{w-y_{i}(t)}+\frac{1}{{\mathfrak{a}}_{i}}.\end{split}

We now couple $\{B_{i}\}_{i=1}^{2n}$ and $\{\overline{B}_{i}\}_{i=1}^{2n}$ so that they equal almost surely. Thereby, we get a coupling between $\{{\bm{x}}(t)\}_{t\in{\mathbb{R}}}$ and $\{{\bm{y}}(t)\}_{t\in{\mathbb{R}}}$ .

The following proposition states that under this coupling, these two line ensembles are close to each other with high probability.

Proposition 6.3.

Fix any $\varepsilon,\theta>0$ and $S>0$ . Then there exist $n,\delta,C$ , such that under the above coupling with probability $1-\varepsilon$ ,

|x_{i}(t)-y_{i}(t)|\leq\theta,\quad\forall t\in[-S,S],\;i\in\llbracket 1,n\rrbracket.

In the following, we will prove that with probability $1-\varepsilon$ ,

(6.3)

\displaystyle x_{i}(t)\leq y_{i}(t)+\theta,\quad\forall t\in[-S,S],\;i\in% \llbracket 1,n\rrbracket.

The lower bound that $x_{i}(t)\geq y_{i}(t)-\theta$ can be proven in the same way.

Our strategy is to consider a shifted version of $\{{\bm{y}}(t)\}_{t\in{\mathbb{R}}}$ , which at $t=-T$ is much larger than ${\bm{x}}(-T)$ ; then we show that it is larger than $\{{\bm{y}}(t)\}_{t\in{\mathbb{R}}}$ for $t\in[-S,S]$ (under the coupling), while the amount of shift is $\leq\theta$ in $[-S,S]$ .

We now define the shifted version of $\{{\bm{y}}(t)\}_{t\in{\mathbb{R}}}$ . For any $t\in{\mathbb{R}}$ and $i\in\mathbb{N}$ , we let

\widetilde{y}_{i}(t)=y_{i}(t)+M-\kappa(t+T),

where $M$ and $\kappa$ are taken as follows. By 1.4 and 2.4, we take $M$ taken large enough (depending only on $\varepsilon$ ) such that with probability $1-\varepsilon/2$ ,

(6.4)

\displaystyle y_{i}(-T)+M>x_{i}(-T),\quad\forall i\in\mathbb{N}.

We then take $\kappa$ such that

M-\kappa(T-S)=\theta.

Then for $n=\lfloor T\rfloor$ large enough (depending on $M,S,\theta$ ), the above choice of parameters imply

(6.5)

\displaystyle\kappa<\frac{2M}{n},\quad M-\kappa(S+T)=\theta-2S\kappa\geq\frac{% \theta}{2}.

We can now rewrite (6.2) in terms of $\{\widetilde{\bm{y}}(t)\}_{t\in{\mathbb{R}}}$ . For each $\ell\in\llbracket 0,2Tn^{3}-1\rrbracket$ and $t\in[t_{\ell}\wedge\tau,t_{\ell+1}\wedge\tau]$ , $i\in\llbracket 1,k_{\ell}^{\prime}\rrbracket$ , we have

(6.6)

\displaystyle\begin{split}{\rm d}\widetilde{y}_{i}(t)&=\sqrt{\frac{2}{\beta}}{% \rm d}B_{i}(t)+\left(\sum_{1\leq j\leq k_{\ell}^{\prime}\atop j\neq i}\frac{1}% {\widetilde{y}_{i}(t)-\widetilde{y}_{j}(t)}+\widetilde{W}_{t}(\widetilde{y}_{i% }(t))\right){\rm d}t,\\ \widetilde{W}_{t}(w)&=\frac{\operatorname{Ai}^{\prime}(0)}{\operatorname{Ai}(0% )}-\kappa+\sum_{i=1}^{k_{\ell}^{\prime}}\frac{1}{{\mathfrak{a}}_{i}}+\sum_{i=k% _{\ell}^{\prime}+1}^{\infty}\frac{1}{w-\widetilde{y}_{i}(t)}+\frac{1}{{% \mathfrak{a}}_{i}}.\end{split}

Lemma 6.4.

There exist $n,\delta,C$ , such under the above coupling the following hodls. Take any $\ell\in\llbracket 0,(S+T)n^{3}\rrbracket$ . Assuming that

\widetilde{y}_{i}(t_{\ell}\wedge\tau)>x_{i}(t_{\ell}\wedge\tau),\quad\forall i% \in\llbracket 1,n\rrbracket,

then

(6.7)

\displaystyle\widetilde{y}_{i}(t)>x_{i}(t),\quad\forall t\in[t_{\ell}\wedge% \tau,t_{\ell+1}\wedge\tau],\;i\in\llbracket 1,n\rrbracket.

Assuming this lemma, we can now finish proving the uniqueness in law of line ensembles.

Proof of 6.3.

As already alluded to, it suffices to prove (6.3). From our choice of $M$ (see (6.4)), we have that $\widetilde{y}_{i}(-T)>x_{i}(-T)$ for all $i\geq 1$ . Then by repeatedly applying 6.4 for $\ell\in\llbracket 0,(S+T)n^{3}\rrbracket$ , and 6.2, we conclude that with probability $1-\varepsilon$ , we have $\tau=\infty$ and $\widetilde{y}_{i}(t)>x_{i}(t)$ for all $t\in[-T,S]$ and $i\in\llbracket 1,n\rrbracket$ . In particular for $t\in[-S,S]$ , this gives

\displaystyle y_{i}(t)+\theta=y_{i}(t)+(M-\kappa(T-S))\geq y_{i}(t)+(M-\kappa(% t+T))>x_{i}(t).

This finishes the proof of (6.3). ∎

Proof of 6.1.

The conclusion follows from taking $\theta,\varepsilon$ to zero and $S$ to infinity in 6.3. ∎

The rest of this section is devoted to proving 6.4. The idea is straightforward: from the coupling we take the difference between (6.1) and (6.6), to cancel out the Brownian motions; and the rest are deterministic arguments.

Proof of 6.4.

For simplicity of notation, in this proof we fix $\ell$ , and write $k=k_{\ell}$ and $k^{\prime}=k_{\ell}^{\prime}$ . Recall that $k,k^{\prime}\in[\![n,2n-1$ ]]. We can take the difference between (6.1) and (6.6), so that for any $i\in\llbracket 1,n\rrbracket$ ,

(6.8)

\displaystyle\begin{split}&\phantom{{}={}}{\rm d}(\widetilde{y}_{i}(t)-x_{i}(t% ))=\sum_{j\in\llbracket 1,n\rrbracket,j\neq i}\frac{(x_{i}(t)-\widetilde{y}_{i% }(t))-(x_{j}(t)-\widetilde{y}_{j}(t))}{(\widetilde{y}_{i}(t)-\widetilde{y}_{j}% (t))(x_{i}(t)-x_{j}(t))}{\rm d}t\\ &+\left(\widetilde{W}_{t}(\widetilde{y}_{i}(t))+\sum_{j=n+1}^{k^{\prime}}\frac% {1}{\widetilde{y}_{i}(t)-\widetilde{y}_{j}(t)}-W_{t}(x_{i}(t))-\sum_{j=n+1}^{k% }\frac{1}{x_{i}(t)-x_{j}(t)}\right){\rm d}t.\end{split}

Denote the stopping time $\sigma$ to be the first time after $t_{\ell}\wedge\tau$ , such that there exists at least one index $i_{*}\in\llbracket 1,n\rrbracket$ with $x_{i_{*}}(\sigma)=\widetilde{y}_{i_{*}}(\sigma)$ (if there were multiple such indices, take $i_{*}$ to be the smallest one). We will prove that $\sigma\geq t_{\ell+1}\wedge\tau$ then (6.7) holds.

We prove by contradiction, and assume that $\sigma<t_{\ell+1}\wedge\tau$ . By the definition of the stopping time $\tau$ , for each $i\in\mathbb{N}$ , $i\geq n/3$ , and $t\in[-T,S\wedge\tau]$ , we have

(6.9)

x_{i}(t)\leq{\mathfrak{a}}_{i}+\frac{C(\log T)^{40}}{i^{\delta}}\leq y_{i}(t)+% \frac{2C(\log T)^{40}}{i^{\delta}}<y_{i}(t)+M-\kappa(t+T)=\widetilde{y}_{i}(t).

We let $a$ (resp. $b$ ) be the smallest (resp. largest) index with $x_{a}(\sigma)=x_{i_{*}}(\sigma)$ (resp. $x_{b}(\sigma)=x_{i_{*}}(\sigma)$ ); and we let $a^{\prime},b^{\prime}\in\mathbb{N}$ be the corresponding indices for $y_{i_{*}}(\sigma)$ . By (6.9), and that $\widetilde{y}_{i}(\sigma)\geq x_{i}(\sigma)$ for each $i\in\llbracket 1,n\rrbracket$ , necessarily $1\leq a^{\prime}\leq a\leq b^{\prime}\leq b<n/2$ . Now for (LABEL:e:yxdiff), by summing over $i\in\llbracket a,b^{\prime}\rrbracket$ , and integrating from $\sigma-\iota$ to $\sigma$ for a sufficiently small $\iota$ , we have

(6.10)

\displaystyle\begin{split}0>&\sum_{i=a}^{b^{\prime}}(\widetilde{y}_{i}(t)-x_{i% }(t))\Big{|}^{\sigma}_{\sigma-\iota}=\int_{\sigma-\iota}^{\sigma}\sum_{i=a}^{b% ^{\prime}}\sum_{j\in\llbracket 1,n\rrbracket\setminus\llbracket a,b^{\prime}% \rrbracket}\frac{(x_{i}(t)-\widetilde{y}_{i}(t))-(x_{j}(t)-\widetilde{y}_{j}(t% ))}{(\widetilde{y}_{i}(t)-\widetilde{y}_{j}(t))(x_{i}(t)-x_{j}(t))}{\rm d}t\\ &+\int_{\sigma-\iota}^{\sigma}\sum_{i=a}^{b^{\prime}}\widetilde{W}_{t}(% \widetilde{y}_{i}(t))+\sum_{j=n+1}^{k^{\prime}}\frac{1}{\widetilde{y}_{i}(t)-% \widetilde{y}_{j}(t)}-W_{t}(x_{i}(t))-\sum_{j=n+1}^{k}\frac{1}{x_{i}(t)-x_{j}(% t)}{\rm d}t.\end{split}

Consider the first term on the RHS of (6.10).

A convergence framework for Airyβ line ensemble via pole evolution

Abstract.

1. Introduction

1.1. Background

1.1.1. Edge limit of general β𝛽\betaitalic_β-ensembles

1.1.2. Airy line ensemble

1.2. Airyβ line ensemble

1.2.1. The edge limit construction and convergence

1.2.2. Description via explicit expressions

1.2.3. Uniqueness and universality

1.3. Main characterization result

Definition 1.1.

Definition 1.2.

Remark 1.3.

Assumption 1.4.

Assumption 1.5.

Theorem 1.6.

1.4. Convergence framework

Theorem 1.7.

1.5. Other properties

1.6. Proof ideas

1.6.1. Propagation of Airy-zero approximation

1.6.2. DBM reconstruction

1.6.3. Uniqueness via coupling

Notations

Acknowledgement

2. Preliminaries and decomposition

Proposition 2.1.

2.1. Airy function

2.2. Estimates on Nevanlinna functions

Lemma 2.2.

Proof.

Lemma 2.3.

Corollary 2.4.

2.3. Proof of 2.1

2.4. Domain extension

Lemma 2.5.

Proof.

2.5. Topological statements

Definition 2.6.

Lemma 2.7.

Proof.

Lemma 2.8.

Proof.

3. Pole evolution: uniform rigidity in time

Proposition 3.1.

Corollary 3.2.

3.1. Characteristic flow

3.2. Estimates for the bulk

Proposition 3.3.

Lemma 3.4.

Proof.

Lemma 3.5.

Proof of 3.3.

Proof of 3.5.

3.3. Estimates for the first pole

Proposition 3.6.

Lemma 3.7.

Proof of 3.6.

Proof of 3.7.

3.4. Proof of 3.1

4. Hölder Regularity

Proposition 4.1.

Proof.

5. Recover Dyson Brownian Motion

Proposition 5.1.

Lemma 5.2.

Proof.

Lemma 5.3.

Proof.

Proof of 5.1.

Remark 5.4.

6. Coupling and Uniqueness

Proposition 6.1.

Lemma 6.2.

Proof.

Proposition 6.3.

Lemma 6.4.

Proof of 6.3.

Proof of 6.1.

A convergence framework for Airy_β line ensemble
via pole evolution

1.1.1. Edge limit of general $\beta$ -ensembles

1.2. Airy_β line ensemble