Modern Stochastics: Theory and Applications logo


  • Help
Login Register

  1. Home
  2. Issues
  3. Volume 9, Issue 1 (2022)
  4. Interacting Brownian motions in infinite ...

Interacting Brownian motions in infinite dimensions related to the origin of the spectrum of random matrices
Volume 9, Issue 1 (2022), pp. 89–122
Yosuke Kawamoto ORCID icon link to view author Yosuke Kawamoto details  

Authors

 
Placeholder
https://doi.org/10.15559/21-VMSTA193
Pub. online: 10 January 2022      Type: Research Article      Open accessOpen Access

Received
11 August 2021
Revised
8 November 2021
Accepted
8 November 2021
Published
10 January 2022

Abstract

The generalised sine random point field arises from the scaling limit at the origin of the eigenvalues of the generalised Gaussian ensembles. We solve an infinite-dimensional stochastic differential equation (ISDE) describing an infinite number of interacting Brownian particles which is reversible with respect to the generalised sine random point field. Moreover, finite particle approximation of the ISDE is shown, that is, a solution to the ISDE is approximated by solutions to finite-dimensional SDEs describing finite-particle systems related to the generalised Gaussian ensembles.

1 Introduction

Consider the unitary ensembles of random matrices whose density is given by
(1.1)
\[ \frac{1}{\mathcal{Z}}|\mathrm{det}M{|^{2\alpha }}{e^{-\mathrm{Tr}{M^{\ast }}M}}dM\hspace{1em}\hspace{2.5pt}\text{for}\hspace{2.5pt}\alpha >-\frac{1}{2},\]
where $dM$ is the usual flat Lebesgue measure on the space of $N\times N$ Hermite matrices, and $\mathcal{Z}$ is the normalising constant. Hereafter, by abuse of notation, we use the same letter $\mathcal{Z}$ for normalising constants of several ensembles. Let ${\mathbf{x}_{N}}=({x_{1}},\dots ,{x_{N}})\in {\mathbb{R}^{N}}$ be the eigenvalues of an $N\times N$ matrix M. When M is distributed as (1.1), the probability density function of its eigenvalues is written as the following ${m_{G,\alpha }^{N}}$, which is called generalised Gaussian ensembles [20]:
(1.2)
\[ {m_{G,\alpha }^{N}}(d{\mathbf{x}_{N}})=\frac{1}{\mathcal{Z}}\prod \limits_{1\le i<j\le N}|{x_{i}}-{x_{j}}{|^{2}}\prod \limits_{1\le k\le N}|{x_{k}}{|^{2\alpha }}{e^{-{x_{k}^{2}}}}d{\mathbf{x}_{N}}.\]
Then, ${m_{G,\alpha }^{N}}$ naturally gives a random point field ${\mu _{G,\alpha }^{N}}$ in the sense that the labelled density of ${\mu _{G,\alpha }^{N}}$ with respect to the Lebesgue measure is given by ${m_{G,\alpha }^{N}}$.
Note that ${\mu _{G,\alpha }^{N}}$ is a determinantal random point field. Let ${\{{p_{\alpha }^{n}}\}_{n\in \mathbb{N}}}$ be the monic orthogonal polynomials with respect to $|x{|^{2\alpha }}{e^{-{x^{2}}}}dx$, and we set ${h_{n}}={\textstyle\int _{\mathbb{R}}}{({p_{\alpha }^{n}}(x))^{2}}|x{|^{2\alpha }}{e^{-{x^{2}}}}dx$. Let ${K_{G,\alpha }^{N}}:\mathbb{R}\times \mathbb{R}\to \mathbb{R}$ be the determinantal kernel defined as
(1.3)
\[ {K_{G,\alpha }^{N}}(x,y)=|x{|^{\alpha }}|y{|^{\alpha }}{e^{-\frac{{x^{2}}}{2}-\frac{{y^{2}}}{2}}}{\sum \limits_{k=0}^{N-1}}\frac{{p_{\alpha }^{k}}(x){p_{\alpha }^{k}}(y)}{{h_{k}}}.\]
Then, for each n, the n-correlation function ${\rho _{G,\alpha }^{N,n}}$ of ${\mu _{G,\alpha }^{N}}$ is given by
\[ {\rho _{G,\alpha }^{N,n}}({x_{1}},\dots ,{x_{n}})=\underset{1\le i,j\le n}{\det }[{K_{G,\alpha }^{N}}({x_{i}},{x_{j}})].\]
From the Cristoffel–Darboux formula [6, Proposition 5.1.3], we have
\[ {K_{G,\alpha }^{N}}(x,y)=\frac{|x{|^{\alpha }}|y{|^{\alpha }}{e^{-\frac{{x^{2}}}{2}-\frac{{y^{2}}}{2}}}}{{h_{N-1}}}\frac{{p_{\alpha }^{N}}(x){p_{\alpha }^{N-1}}(y)-{p_{\alpha }^{N-1}}(x){p_{\alpha }^{N}}(y)}{x-y}.\]
To focus on fluctuation of the generalised Gaussian ensembles around the origin, we take a scaling. We set
(1.4)
\[ {K_{\alpha }^{N}}(x,y)=\frac{1}{\sqrt{N}}{K_{G,\alpha }^{N}}\Big(\frac{x}{\sqrt{N}},\frac{y}{\sqrt{N}}\Big).\]
Let ${\mu _{\alpha }^{N}}$ be the determinantal random point field with kernel ${K_{\alpha }^{N}}$. Clearly, the labelled density ${m_{\alpha }^{N}}$ of ${\mu _{\alpha }^{N}}$ is written in the form
(1.5)
\[ {m_{\alpha }^{N}}(d{\mathbf{x}_{N}})=\frac{1}{\mathcal{Z}}\prod \limits_{1\le i<j\le N}|{x_{i}}-{x_{j}}{|^{2}}\prod \limits_{1\le k\le N}|{x_{k}}{|^{2\alpha }}{e^{-\frac{{x_{k}^{2}}}{N}}}d{\mathbf{x}_{N}}.\]
Then, the scaled kernel ${K_{\alpha }^{N}}$ has a nontrivial limit
(1.6)
\[ \underset{N\to \infty }{\lim }{K_{\alpha }^{N}}(x,y)={K_{\alpha }}(x,y)\hspace{1em}\hspace{2.5pt}\text{compact uniformly}.\]
Here, the limit kernel ${K_{\alpha }}(x,y):\mathbb{R}\times \mathbb{R}\to \mathbb{R}$ is given by, for $x\ne y$,
(1.7)
\[\begin{aligned}{}& {K_{\alpha }}(x,y)\\ {} & =\frac{x\sqrt{|{x^{-1}}y|}{J_{\alpha +\frac{1}{2}}}(\sqrt{2}|x|){J_{\alpha -\frac{1}{2}}}(\sqrt{2}|y|)-y\sqrt{|x{y^{-1}}|}{J_{\alpha -\frac{1}{2}}}(\sqrt{2}|x|){J_{\alpha +\frac{1}{2}}}(\sqrt{2}|y|)}{\sqrt{2}(x-y)},\end{aligned}\]
and, for $x=y$,
(1.8)
\[\begin{aligned}{}{K_{\alpha }}(x,x)& =\frac{|x|}{2}\big\{{J_{\alpha +\frac{1}{2}}^{2}}(\sqrt{2}|x|)+{J_{\alpha -\frac{1}{2}}^{2}}(\sqrt{2}|x|)\\ {} & \hspace{2em}-{J_{\alpha +\frac{3}{2}}}(\sqrt{2}|x|){J_{\alpha -\frac{1}{2}}}(\sqrt{2}|x|)-{J_{\alpha +\frac{1}{2}}}(\sqrt{2}|x|){J_{\alpha -\frac{3}{2}}}(\sqrt{2}|x|)\big\},\end{aligned}\]
where ${J_{\nu }}$ denotes the Bessel function of the first kind of order ν. Remark that (1.6) yields weak convergence of random point fields. Let ${\mu _{\alpha }}$ be the determinantal random point field whose kernel is ${K_{\alpha }}$. More precisely, the n-correlation function ${\rho _{\alpha }^{n}}$ of ${\mu _{\alpha }}$ is given by
(1.9)
\[ {\rho _{\alpha }^{n}}({x_{1}},\dots ,{x_{n}})=\underset{1\le i,j\le n}{\det }[{K_{\alpha }}({x_{i}},{x_{j}})].\]
Then, we obtain
\[ \underset{N\to \infty }{\lim }{\mu _{\alpha }^{N}}={\mu _{\alpha }}\hspace{2.5pt}\text{weakly}.\]
The kernel ${K_{\alpha }}$ is called the generalised sine kernel of order α [1]: when $\alpha =0$, ${K_{\alpha }}$ becomes the classical sine kernel $\sin (x-y)/(x-y)$. Note that, for $x,y>0$ or $\alpha \in \mathbb{N}$, the kernel (1.7) takes a simpler form
\[ {K_{\alpha }}(x,y)=\sqrt{xy}\frac{{J_{\alpha +\frac{1}{2}}}(\sqrt{2}x){J_{\alpha -\frac{1}{2}}}(\sqrt{2}y)-{J_{\alpha -\frac{1}{2}}}(\sqrt{2}x){J_{\alpha +\frac{1}{2}}}(\sqrt{2}y)}{\sqrt{2}(x-y)}.\]
As with other random matrix models, the universality of the limit kernel has been studied. Let $V:\mathbb{R}\to \mathbb{R}$ be a function that increases fast enough at infinity, for instance, an even degree polynomial with a positive leading coefficient. Then, we consider the following ${m_{V,\alpha }^{N}}$, which is a generalisation of (1.2):
\[ {m_{V,\alpha }^{N}}(d{\mathbf{x}_{N}})=\frac{1}{\mathcal{Z}}\prod \limits_{1\le i<j\le N}|{x_{i}}-{x_{j}}{|^{2}}\prod \limits_{1\le k\le N}|{x_{k}}{|^{2\alpha }}{e^{-V({x_{k}})}}d{\mathbf{x}_{N}}.\]
Then, a natural problem is the universality of random matrices. More precisely, it is of interest to show that, for a potential V of a quite wide class, local fluctuation at the origin for ${m_{V,\alpha }^{N}}$ yields the universal random point field ${\mu _{\alpha }}$ under suitable scaling. Akemann et al. showed the universality under the assumptions that α is a nonnegative integer, V is an even degree polynomial, and recurrence coefficients of associated orthogonal polynomials satisfy certain condition [1]. The assumption that α is a nonnegative integer was excepted by Kanzieper and Freilikher [10]. Moreover, the restriction that V is even degree polynomial was removed by Kuijlaars and Vanlessen. They showed universality for real analytic potential with mild assumptions [15].
It is natural to try to study stochastic dynamics with infinitely many particles related to the universal random point field ${\mu _{\alpha }}$. More specifically, our purpose is to construct reversible diffusion whose equilibrium measure is ${\mu _{\alpha }}$, and describe its infinite-dimensional stochastic differential equation (ISDE). In order to derive such an ISDE, we consider a stochastic process related to ${\mu _{\alpha }^{N}}$, which is given by
(1.10)
\[\begin{aligned}{}d{X_{t}^{N,i}}& =d{B_{t}^{i}}+\Big\{\frac{\alpha }{{X_{t}^{N,i}}}-\frac{{X_{t}^{N,i}}}{N}+{\sum \limits_{j\ne i}^{N}}\frac{1}{{X_{t}^{N,i}}-{X_{t}^{N,j}}}\Big\}dt,\hspace{1em}1\le i\le N,\end{aligned}\]
where ${({B^{i}})_{i=1}^{N}}$ is the N-dimensional Brownian motion. We derive (1.10) as follows. Consider a Dirichlet form
\[ \mathcal{E}(f,g)={\int _{{\mathbb{R}^{N}}}}\frac{1}{2}{\sum \limits_{i=1}^{N}}\frac{\partial f}{\partial {x_{i}}}\frac{\partial g}{\partial {x_{i}}}{m_{\alpha }^{N}}(d{\mathbf{x}_{N}})\]
on ${L^{2}}({\mathbb{R}^{N}},{m_{\alpha }^{N}}(d{\mathbf{x}_{N}}))$. Then, by using integration by parts and (1.5), the generator $-{A^{N}}$ of $\mathcal{E}$ on ${L^{2}}({\mathbb{R}^{N}},{m_{\alpha }^{N}}(d{\mathbf{x}_{N}}))$ is given by
\[ {A^{N}}=\frac{1}{2}\Delta +{\sum \limits_{i=1}^{N}}\Big\{\frac{\alpha }{{x_{i}}}-\frac{{x_{i}}}{N}+{\sum \limits_{j\ne i}^{N}}\frac{1}{{x_{i}}-{x_{j}}}\Big\}\frac{\partial }{\partial {x_{i}}},\]
which corresponds to (1.10).
Taking N to infinity in (1.10), it is expected that the limit ISDE is given by
(1.11)
\[\begin{aligned}{}d{X_{t}^{i}}& =d{B_{t}^{i}}+\Big\{\frac{\alpha }{{X_{t}^{i}}}+\underset{r\to \infty }{\lim }\sum \limits_{|{X_{t}^{i}}-{X_{t}^{j}}|<r,j\ne i}\frac{1}{{X_{t}^{i}}-{X_{t}^{j}}}\Big\}dt,\hspace{1em}i\in \mathbb{N}.\end{aligned}\]
Here, ${({B^{i}})_{i\in \mathbb{N}}}$ is the collection of independent copies of the one-dimensional Brownian motion. Clearly, (1.11) becomes the Dyson model in the case $\alpha =0$. Our purpose is to construct a solution ${({X^{i}})_{i\in \mathbb{N}}}$ to (1.11).
It is important to point out that (1.11) is essentially different from the Bessel interacting ISDE, and deriving (1.11) from (1.10) is nontrivial for the long-range effect of logarithmic interaction. Honda and Osada derived the following Bessel interacting ISDE for $\alpha >1/2$ [9]:
(1.12)
\[\begin{aligned}{}d{X_{t}^{i}}& =d{B_{t}^{i}}+\Big\{\frac{\alpha }{{X_{t}^{i}}}+\sum \limits_{j\in \mathbb{N},j\ne i}\frac{1}{{X_{t}^{i}}-{X_{t}^{j}}}\Big\}dt,\hspace{1em}i\in \mathbb{N}.\end{aligned}\]
Note that the (unlabelled) solution to (1.12) is reversible with respect to the Bessel random point field. Seemingly, our model (1.11) resembles (1.12). However, the sum of the drift term in (1.12) converges absolutely unlike (1.11). The Bessel random point field is a random point field on the positive line, and its one-correlation function, which stands for the mean density of particles, decreases with order ${x^{-1/2}}$ as $x\to \infty $; therefore, the sum of the interaction term in (1.12) converges absolutely. In contrast, the generalised sine random point field is on $\mathbb{R}$, and its one-correlation function is of constant order as $|x|\to \infty $; this implies that the drift term in (1.11) does not converge absolutely. Thus, we need careful computation to show the conditional convergence of the drift term as opposed to (1.12), and (1.11) is more difficult to study than (1.12) in this respect.
We refer to historical remarks on interacting Brownian motions with infinitely many particles. For a free potential Φ and an interaction potential Ψ, interacting Brownian motions in infinite-dimensions are given by the ISDE of the form
(1.13)
\[\begin{aligned}{}d{X_{t}^{i}}& =d{B_{t}^{i}}-\frac{\beta }{2}{\nabla _{x}}\Phi ({X_{t}^{i}})dt-\frac{\beta }{2}{\sum \limits_{j\ne i}^{\infty }}{\nabla _{x}}\Psi ({X_{t}^{i}},{X_{t}^{j}})dt,\hspace{1em}i\in \mathbb{N}.\end{aligned}\]
Here, ${({B^{i}})_{i\in \mathbb{N}}}$ is the collection of independent copies of d-dimensional Brownian motion and $\beta >0$ is the inverse temperature. Lang derived general solutions to ISDE (1.13) under the condition $\Phi =0$ and $\Psi \in {C_{0}^{3}}({\mathbb{R}^{d}})$ [16, 17]. Fritz explicitly described the set of starting points for up to four dimensions [7], and Tanemura solved equations for hardcore Brownian balls [31]. These results were achieved through Itô’s method, and required the coefficients to be smooth and have compact support. These conditions exclude physically interesting examples of long-range interaction potentials, such as the Lennard-Jones 6-12 potential and Riesz potentials. In particular, the logarithmic potential, which appears in random matrix theory, is also excluded. Osada established a Dirichlet form approach to solve ISDEs with a long-range potential, which can be applied to the logarithmic interacting particle systems [22, 24–27]. Osada and Tanemura also showed strong uniqueness of solutions [29]. Furthermore, Tsai constructed nonequilibrium solutions to the Dyson model for $1\le \beta <\infty $ [32].
We use the Dirichlet form approach to construct the unique strong solution to (1.11). More precisely, we construct a reversible diffusion on the configuration space with respect to ${\mu _{\alpha }}$, and show that the process satisfies (1.11). To show ISDEs by the Dirichlet form approach, expression of the logarithmic derivatives is crucial, because the logarithmic derivatives correspond the drift terms of ISDEs. Bufetov, Dymov, and Osada introduced a method to compute the logarithmic derivatives for determinantal random point fields [2]. Since ${\mu _{\alpha }}$ is determinantal, their result seems to be applicable to our case. In spite of this, we use finite-particle approximation method to find the logarithmic derivative of ${\mu _{\alpha }}$, which was introduced in [25]. This is because finite-particle approximation of the logarithmic derivative implies that of dynamics [12]. More accurately, let ${({X^{N,i}})_{1\le i\le N}}$ and ${({X^{i}})_{i\in \mathbb{N}}}$ be the solutions to (1.10) and (1.11), respectively. Then, provided suitable labelling, the approximation of the logarithmic derivatives imply that the first m-particles of ${({X^{N,i}})_{1\le i\le N}}$ converge to that of ${({X^{i}})_{i\in \mathbb{N}}}$ weakly in the path space. We will show such convergence in Theorem 2.4. This is an advantage of our method.
This paper is organised as follows. In Section 2, main results are shown. In Section 3, we prepare estimates of determinantal kernels. In Section 4 and Section 5, the logarithmic derivative and quasi-Gibbs property of ${\mu _{\alpha }}$ are shown, respectively. In Section 6, the strong uniqueness of our solution is proved.

2 Set up and main results

2.1 The strong uniqueness of solutions to ISDEs and dynamical convergence

Let $S=\mathbb{R}$. Let $\mathsf{S}$ be the configuration space over S given by
\[\begin{aligned}{}& \mathsf{S}=\{\mathsf{s}=\sum \limits_{i}{\delta _{{s_{i}}}}\hspace{0.1667em};\hspace{0.1667em}{s_{i}}\in S,\mathsf{s}(K)<\infty \hspace{2.5pt}\text{for any compact set}\hspace{2.5pt}K\subset S\}.\end{aligned}\]
We regard the zero measure as an element of $\mathsf{S}$ by convention. We equip $\mathsf{S}$ with the vague topology, which makes $\mathsf{S}$ to be a Polish space. A probability measure μ on $(\mathsf{S},\mathfrak{B}(\mathsf{S}))$ is called a random point field on S. Let
(2.1)
\[\begin{aligned}{}& {\mathsf{S}_{\mathrm{s}}}=\{\mathsf{s}\in \mathsf{S}\hspace{0.1667em};\hspace{0.1667em}\mathsf{s}(\{x\})\le 1\hspace{2.5pt}\text{for any}\hspace{2.5pt}x\in S\},\hspace{2em}{\mathsf{S}_{\mathrm{i}}}=\{\mathsf{s}\in \mathsf{S}\hspace{0.1667em};\hspace{0.1667em}\mathsf{s}(S)=\infty \},\\ {} & {\mathsf{S}_{\mathrm{s},\mathrm{i}}}={\mathsf{S}_{\mathrm{s}}}\cap {\mathsf{S}_{\mathrm{i}}}.\end{aligned}\]
For each $n\in \mathbb{N}$, a symmetric and locally integrable function ${\rho ^{n}}:{S^{n}}\to [0,\infty )$ is called the n-correlation function of μ with respect to the Lebesgue measure if ${\rho ^{n}}$ satisfies
\[\begin{aligned}{}& {\int _{{A_{1}^{{k_{1}}}}\times \cdots \times {A_{m}^{{k_{m}}}}}}{\rho ^{n}}({x_{1}},\dots ,{x_{n}})d{x_{1}}\cdots d{x_{n}}={\int _{\mathsf{S}}}{\prod \limits_{i=1}^{m}}\frac{\mathsf{s}({A_{i}})!}{(\mathsf{s}({A_{i}})-{k_{i}})!}d\mu (\mathsf{s})\end{aligned}\]
for any sequence of disjoint bounded sets ${A_{1}},\dots ,{A_{m}}\in \mathfrak{B}(S)$ and any sequence of natural numbers ${k_{1}},\dots ,{k_{m}}$ satisfying ${k_{1}}+\cdots +{k_{m}}=n$.
We define an unlabelling map $\mathfrak{u}:\{{\textstyle\bigcup _{k=0}^{\infty }}{S^{k}}\}\cup {S^{\mathbb{N}}}\to \mathsf{S}$ as $\mathfrak{u}({({s_{i}})_{i}})={\textstyle\sum _{i}}{\delta _{{s_{i}}}}$. Here, ${S^{0}}=\{\varnothing \}$ and $\mathfrak{u}(\varnothing )=\mathfrak{o}$, where $\mathfrak{o}$ is the zero measure. Furthermore, a measurable function $\mathfrak{l}:{\mathsf{S}_{\mathrm{s},\mathrm{i}}}\to {S^{\mathbb{N}}}$ is called a labelling map if $\mathfrak{u}\circ \mathfrak{l}$ is the identity map.
Recall that ${\mu _{\alpha }}$ is the random point field whose correlation functions are given by (1.9). The next theorem is the main result of the present paper.
Theorem 2.1.
Assume $\alpha >1/2$. Then, there exists a set ${\mathsf{S}_{\alpha }}$ satisfying
\[ {\mu _{\alpha }}({\mathsf{S}_{\alpha }})=1,\hspace{1em}{\mathsf{S}_{\alpha }}\subset {\mathsf{S}_{\mathrm{s},\mathrm{i}}},\]
such that the following holds: for any $\mathbf{s}\in {\mathfrak{u}^{-1}}({\mathsf{S}_{\alpha }})$, there exists an ${S^{\infty }}$-valued continuous process $\mathbf{X}={({X^{i}})_{i\in \mathbb{N}}}$ and an ${\mathbb{R}^{\mathbb{N}}}$-valued Brownian motion $\mathbf{B}={({B^{i}})_{i\in \mathbb{N}}}$ satisfying
(2.2)
\[\begin{aligned}{}d{X_{t}^{i}}& =d{B_{t}^{i}}+\Big\{\frac{\alpha }{{X_{t}^{i}}}+\underset{r\to \infty }{\lim }\sum \limits_{|{X_{t}^{i}}-{X_{t}^{j}}|<r}\frac{1}{{X_{t}^{i}}-{X_{t}^{j}}}\Big\}dt,\end{aligned}\]
(2.3)
\[\begin{aligned}{}{\mathbf{X}_{0}}& =\mathbf{s},\end{aligned}\]
and
\[ P(\mathfrak{u}({\mathbf{X}_{t}})\in {\mathsf{S}_{\alpha }},0{\le ^{\forall }}t<\infty )=1.\]
The solutions obtained in Theorem 2.1 are just weak solutions. We next show the strong uniqueness of solutions to (2.2). We shall give the definitions of (IFC), (AC), (SIN), (NBJ), and (MF) in Section 6.
Theorem 2.2.
For ${\mu _{\alpha }}\circ {\mathfrak{l}^{-1}}$-a.s. s, there exists a strong solution $(\mathbf{X},\mathbf{B})$ to (2.2)–(2.3) satisfying (IFC), (${\mu _{\alpha }}$-AC), (SIN), and (NBJ). Furthermore, the strong uniqueness of solutions to (2.2)–(2.3) holds under the constraints of (MF), (IFC), (${\mu _{\alpha }}$-AC), (SIN), and (NBJ).
Remark 2.3.
Five assumptions (IFC), (AC), (SIN), (NBJ), and (MF) are required for the uniqueness of solutions in Theorem 2.2. It is not difficult to check these assumptions in practice. In fact, we showed that solutions to ISDEs satisfy (IFC) under mild assumptions, and furthermore, if a solution is associated with a Dirichlet form, we can check the five assumptions [13, 14]. Thus, as an application of the uniqueness of solutions in Theorem 2.2, we can show the uniqueness of Dirichlet forms associated with the generalised sine random point field [13], but we do not pursue here.
Let $(\mathbf{X},\mathbf{B})$ be the unique solution obtained in Theorem 2.2 and we write $\mathbf{X}={({X^{i}})_{i\in \mathbb{N}}}$. Let ${\mathbf{X}^{N}}={({X^{N,i}})_{1\le i\le N}}$ be the (unique) solution to (1.10). In analogy to the labelling map $\mathfrak{l}$, let ${\mathfrak{l}^{N}}:\{\mathsf{s}\in \mathsf{S}\hspace{0.1667em};\hspace{0.1667em}\mathsf{s}(S)=N\}\to {S^{N}}$ be a labelling map for N-particles. For labelling maps $\mathfrak{l}={({l^{i}})_{i\in \mathbb{N}}}$ and ${\mathfrak{l}^{N}}={({l^{N,i}})_{1\le i\le N}}$, we write ${\mathfrak{l}_{m}}={({l^{i}})_{1\le i\le m}}$ and ${\mathfrak{l}_{m}^{N}}={({l^{N,i}})_{1\le i\le m}}$, respectively. Here and subsequently, $W(A)$ stands for the set of all continuous paths $w:[0,\infty )\to A$. Then, we see a weak convergence from ${\mathbf{X}^{N}}$ to X in a path space.
Theorem 2.4.
Assume $\alpha >1/2$. Suppose that, for each $m\in \mathbb{N}$,
(2.4)
\[ \underset{N\to \infty }{\lim }{\mu _{\alpha }^{N}}\circ {({\mathfrak{l}_{m}^{N}})^{-1}}={\mu _{\alpha }}\circ {({\mathfrak{l}_{m}})^{-1}}\hspace{2.5pt}\textit{weakly}.\]
Assume that ${\mathbf{X}_{0}^{N}}={\mu _{\alpha }^{N}}\circ {({\mathfrak{l}^{N}})^{-1}}$ and ${\mathbf{X}_{0}}={\mu _{\alpha }}\circ {\mathfrak{l}^{-1}}$ in distribution. Then, for each $m\in \mathbb{N}$,
\[ \underset{N\to \infty }{\lim }({X^{N,1}},{X^{N,2}},\dots ,{X^{N,m}})=({X^{1}},{X^{2}},\dots ,{X^{m}})\]
weakly in $W({S^{m}})$.

2.2 Quasi-Gibbs property and logarithmic derivative

Theorem 2.1 is proved by making use of the Dirichlet form approach. In this framework, quasi-Gibbs property and logarithmic derivatives play important roles. We first introduce the concept of quasi-Gibbs measures.
We set ${S_{r}}=\{|x|\le r\}$ and ${\mathsf{S}_{r}^{m}}=\{\mathsf{s}\in \mathsf{S}\hspace{0.1667em};\hspace{0.1667em}\mathsf{s}({S_{r}})=m\}$. Let ${\Lambda _{r}}$ be the Poisson random point field whose intensity is the Lebesgue measure on ${S_{r}}$ and set ${\Lambda _{r}^{m}}={\Lambda _{r}}(\cdot \cap {\mathsf{S}_{r}^{m}})$. Let ${\pi _{r}},{\pi _{r}^{c}}:\mathsf{S}\to \mathsf{S}$ be such that ${\pi _{r}}(\mathsf{s})=\mathsf{s}(\cdot \cap {S_{r}})$ and ${\pi _{r}^{c}}(\mathsf{s})=\mathsf{s}(\cdot \cap {S_{r}^{c}})$, respectively.
Definition 2.5 ([26, 27]).
A random point field μ is said to be a $(\Phi ,\Psi )$-quasi-Gibbs measure if, for each $r,m\in \mathbb{N}$ and μ-a.s. $\mathsf{s}$, the regular conditional probability
\[\begin{aligned}{}& {\mu _{r,\mathsf{s}}^{m}}:=\mu (\hspace{0.1667em}{\pi _{r}}(\mathsf{x})\in \cdot \hspace{0.1667em}|\hspace{0.1667em}{\pi _{r}^{c}}(\mathsf{x})={\pi _{r}^{c}}(\mathsf{s}),\hspace{0.1667em}\mathsf{x}({S_{r}})=m)\end{aligned}\]
satisfies
\[ {c_{\text{1}}^{-1}}{e^{-{\mathcal{H}_{r}}(\mathsf{x})}}{\Lambda _{r}^{m}}(d\mathsf{x})\le {\mu _{r,\mathsf{s}}^{m}}(d\mathsf{x})\le {c_{\text{1}}}{e^{-{\mathcal{H}_{r}}(\mathsf{x})}}{\Lambda _{r}^{m}}(d\mathsf{x}).\]
Here, ${c_{1}}={c_{\text{1}}}(r,m,\mathsf{s})$ is a positive constant depending on r, m, and $\mathsf{s}$, and ${\mathcal{H}_{r}}$ is the Hamiltonian on ${S_{r}}$ defined as
\[\begin{aligned}{}& {\mathcal{H}_{r}}(\mathsf{x})=\sum \limits_{{x_{i}}\in {S_{r}}}\Phi ({x_{i}})+\sum \limits_{{x_{j}},{x_{k}}\in {S_{r}}}\Psi ({x_{j}},{x_{k}})\hspace{1em}\hspace{2.5pt}\text{for}\hspace{2.5pt}\mathsf{x}=\sum \limits_{i}{\delta _{{x_{i}}}}.\end{aligned}\]
Moreover, for two measures $\mu ,\nu $ on a σ-field $\mathcal{F}$, we write $\mu \le \nu $ if $\mu (A)\le \nu (A)$ for all $A\in \mathcal{F}$.
Theorem 2.6.
For any $\alpha >1/2$, ${\mu _{\alpha }}$ is a $(-2\alpha \log |x|,-2\log |x-y|)$-quasi-Gibbs measure.
From Theorem 2.6, a reversible diffusion with respect to ${\mu _{\alpha }}$ is constructed by the Dirichlet form theory. Let ${\mathcal{E}^{{\mu _{\alpha }}}}$ and ${\mathcal{D}_{\circ }^{{\mu _{\alpha }}}}$ be as in (2.8) and (2.9) with $k=0$ and $\mu ={\mu ^{\alpha }}$, respectively. See, e.g., [8, 19] for notions of the general theory of Dirichlet forms such as locality, quasi-regularity, associated diffusions, and capacity.
Corollary 2.7.
Assume $\alpha >1/2$. Then, $({\mathcal{E}^{{\mu _{\alpha }}}},{\mathcal{D}_{\circ }^{{\mu _{\alpha }}}})$ is closable on ${L^{2}}({\mu _{\alpha }})$, and its closure $({\mathcal{E}^{{\mu _{\alpha }}}},{\mathcal{D}^{{\mu _{\alpha }}}})$ is a local and quasi-regular Dirichlet form. Thus, there exists an $\mathsf{S}$-valued diffusion $(\mathsf{X},{\{{\mathsf{P}_{\mathsf{s}}}\}_{\mathsf{s}\in \mathsf{S}}})$ associated with $({\mathcal{E}^{{\mu _{\alpha }}}},{\mathcal{D}^{{\mu _{\alpha }}}})$.
Proof.
This corollary immediately follows from Theorem 2.6 with [26, Lemma 2.1, Corollary 2.1].  □
We remark that each particle does not hit the origin. Let ${\mathrm{Cap}^{{\mu _{\alpha }}}}$ be the one-capacity with respect to $({\mathcal{E}^{{\mu _{\alpha }}}},{\mathcal{D}^{{\mu _{\alpha }}}},{L^{2}}({\mu _{\alpha }}))$.
Lemma 2.8.
Assume $\alpha >1/2$. Let $\mathsf{A}=\{\mathsf{s}\hspace{0.1667em};\hspace{0.1667em}\mathsf{s}(\{0\})\ge 1\}$. Then,
\[ {\mathrm{Cap}^{{\mu _{\alpha }}}}(\mathsf{A})=0.\]
Proof.
From (1.8), we see that, for each r, there exists a positive constant ${c_{2}}$ such that
(2.5)
\[ {\rho _{\alpha }^{1}}(x)\le {c_{\text{2}}}|x{|^{2\alpha }}\hspace{1em}\hspace{2.5pt}\text{for}\hspace{2.5pt}x\in {S_{r}}.\]
Once we have obtained (2.5), the lemma is proved by the same argument as in [9, Lemma B.1].  □
We next prepare some quantities to introduce the logarithmic derivative, which is a crucial quantity for the representation of ISDE. For a random point field μ and $\mathbf{x}=({x_{1}},\dots ,{x_{k}})\in S$, for some $k\in \mathbb{N}$, we call ${\mu _{\mathbf{x}}}$ the (reduced) Palm measure of μ conditioned at x if ${\mu _{\mathbf{x}}}$ is the regular conditional probability defined as
\[\begin{aligned}{}& {\mu _{\mathbf{x}}}=\mu \Big(\hspace{0.1667em}\cdot \hspace{0.1667em}-{\sum \limits_{i=1}^{k}}{\delta _{{x_{i}}}}\hspace{0.1667em}|\hspace{0.1667em}\mathsf{s}(\{{x_{i}}\})\ge 1\hspace{2.5pt}\text{for}\hspace{2.5pt}i=1,\dots ,k\Big).\end{aligned}\]
Recall that ${\rho ^{k}}$ denotes the k-correlation function of μ. A Radon measure ${\mu ^{[k]}}$ on ${S^{k}}\times \mathsf{S}$ is called the k-Campbell measure of μ if ${\mu ^{[k]}}$ is given by
\[\begin{aligned}{}& {\mu ^{[k]}}(d\mathbf{x}d\mathsf{s})={\rho ^{k}}(\mathbf{x}){\mu _{\mathbf{x}}}(d\mathsf{s})d\mathbf{x}.\end{aligned}\]
Let ${L_{\mathrm{loc}}^{p}}(S\times \mathsf{S},{\mu ^{[1]}})={\textstyle\bigcap _{r=1}^{\infty }}{L^{p}}({S_{r}}\times \mathsf{S},{\mu ^{[1]}})$.
A function $f:\mathsf{S}\to \mathbb{R}$ is said to be local if f is $\sigma [{\pi _{r}}]$-measurable for some $r\in \mathbb{N}$. Moreover, we say that f is smooth if $\check{f}$ is smooth, where $\check{f}$ is the permutation invariant function in ${({s_{i}})_{i}}$ such that $f(\mathsf{s})=\check{f}({({s_{i}})_{i}})$ for $\mathsf{s}={\textstyle\sum _{i}}{\delta _{{s_{i}}}}$. Let ${\mathcal{D}_{\circ }}$ be the set of all local smooth functions on $\mathsf{S}$, and we write ${\mathcal{D}_{\circ ,b}}=\{f\in {\mathcal{D}_{\circ }}\hspace{0.1667em};\hspace{0.1667em}f\hspace{2.5pt}\text{is bounded}\}$. We set
\[\begin{aligned}{}& {C_{0}^{\infty }}(\mathbb{R}\setminus \{0\})\otimes {\mathcal{D}_{\circ ,b}}=\big\{{\sum \limits_{i=1}^{m}}{f_{i}}(x){g_{i}}(\mathsf{y})\hspace{0.1667em};\hspace{0.1667em}{f_{i}}\in {C_{0}^{\infty }}(\mathbb{R}\setminus \{0\}),\hspace{2.5pt}{g_{i}}\in {\mathcal{D}_{\circ ,b}},\hspace{0.1667em}m\in \mathbb{N}\big\}.\end{aligned}\]
Definition 2.9 ([25]).
An $\mathbb{R}$-valued function ${\mathsf{d}^{\mu }}\in {L_{\mathrm{loc}}^{1}}(S\times \mathsf{S},{\mu ^{[1]}})$ is called the logarithmic derivative of μ if, for all $\varphi \in {C_{0}^{\infty }}(\mathbb{R}\setminus \{0\})\otimes {\mathcal{D}_{\circ ,b}}$,
\[\begin{aligned}{}& {\int _{S\times \mathsf{S}}}{\mathsf{d}^{\mu }}(x,\mathsf{y})\varphi (x,\mathsf{y}){\mu ^{[1]}}(dxd\mathsf{y})=-{\int _{S\times \mathsf{S}}}{\nabla _{x}}\varphi (x,\mathsf{y}){\mu ^{[1]}}(dxd\mathsf{y}).\end{aligned}\]
Remark 2.10.
The space of test functions is not ${C_{0}^{\infty }}(\mathbb{R})\otimes {\mathcal{D}_{\circ ,b}}$, but ${C_{0}^{\infty }}(\mathbb{R}\setminus \{0\})\otimes {\mathcal{D}_{\circ ,b}}$, since particles never hit the origin for the case $\mu ={\mu _{\alpha }}$ with $\alpha >1/2$ from Lemma 2.8.
The next claim is the key theorem in the present paper.
Theorem 2.11.
For $\alpha >1/2$, the logarithmic derivative of ${\mu _{\alpha }}$ exists in ${L_{\mathrm{loc}}^{p}}(S\times \mathsf{S},{\mu _{\alpha }^{[1]}})$ for $1\le p<2$, and it is given by
(2.6)
\[ {\mathsf{d}^{{\mu _{\alpha }}}}(x,\mathsf{y})=\frac{2\alpha }{x}+\underset{r\to \infty }{\lim }\sum \limits_{|x-{y_{j}}|<r}\frac{2}{x-{y_{j}}}\]
for $\mathsf{y}={\textstyle\sum _{j}}{\delta _{{y_{j}}}}$. Here, the limit in the right hand side of (2.6) is taken in ${L_{\mathrm{loc}}^{p}}(S\times \mathsf{S},{\mu _{\alpha }^{[1]}})$.

2.3 The general theory for ISDEs

This subsection is devoted to the general framework – the Dirichlet form approach for infinite particle systems. We first introduce Dirichlet forms describing k-labelled processes.
For $f,g\in {\mathcal{D}_{\circ }}$, we define $\mathbb{D}[f,g]:\mathsf{S}\to \mathbb{R}$ as
(2.7)
\[ \mathbb{D}[f,g](\mathsf{s})=\frac{1}{2}\sum \limits_{i}\frac{\partial \check{f}(\mathbf{s})}{\partial {s_{i}}}\frac{\partial \check{g}(\mathbf{s})}{\partial {s_{i}}}.\]
Here $\mathsf{s}={\textstyle\sum _{i}}{\delta _{{s_{i}}}}$ and $\mathbf{s}={({s_{i}})_{i}}$. Remark that $\mathbb{D}[f,g]$ is well defined, because the right hand side of (2.7) depends only on $\mathsf{s}$. For $f,g\in {C_{0}^{\infty }}({S^{k}})\otimes {\mathcal{D}_{\circ }}$, let ${\nabla ^{k}}[f,g]$ be the function on ${S^{k}}\times \mathsf{S}$ defined as
\[ {\nabla ^{k}}[f,g](\mathbf{x},\mathsf{s})=\frac{1}{2}{\sum \limits_{i=1}^{k}}\frac{\partial f(\mathbf{x},\mathsf{s})}{\partial {x_{i}}}\frac{\partial g(\mathbf{x},\mathsf{s})}{\partial {x_{i}}},\]
where $\mathbf{x}={({x_{i}})_{i=1}^{k}}\in {S^{k}}$. For $f,g\in {C_{0}^{\infty }}({S^{k}})\otimes {\mathcal{D}_{\circ }}$, we set
\[ {\mathbb{D}^{k}}[f,g](\mathbf{x},\mathsf{s})={\nabla ^{k}}[f,g](\mathbf{x},\mathsf{s})+\mathbb{D}[f(\mathbf{x},\cdot ),g(\mathbf{x},\cdot )](\mathsf{s}).\]
Then, we set the bilinear form $({\mathcal{E}^{{\mu ^{[k]}}}},{\mathcal{D}_{\circ }^{{\mu ^{[k]}}}})$ as
(2.8)
\[\begin{aligned}{}{\mathcal{E}^{{\mu ^{[k]}}}}(f,g)& ={\int _{{S^{k}}\times \mathsf{S}}}{\mathbb{D}^{k}}[f,g]d{\mu ^{[k]}},\end{aligned}\]
(2.9)
\[\begin{aligned}{}{\mathcal{D}_{\circ }^{{\mu ^{[k]}}}}& =\{f\in ({C_{0}^{\infty }}({S^{k}})\otimes {\mathcal{D}_{\circ }})\cap {L^{2}}({S^{k}}\times \mathsf{S},{\mu ^{[k]}})\hspace{0.1667em};\hspace{0.1667em}{\mathcal{E}^{{\mu ^{[k]}}}}(f,f)<\infty \}.\end{aligned}\]
When $k=0$, we interpret ${\mathbb{D}^{0}}=\mathbb{D}$ and ${\mu ^{[0]}}=\mu $. For simplicity, we write ${L^{2}}({\mu ^{[k]}})={L^{2}}({S^{k}}\times \mathsf{S},{\mu ^{[k]}})$.
Let $\mathrm{Erf}(t)=(1/\sqrt{2\pi }){\textstyle\int _{t}^{\infty }}{e^{-{x^{2}}/2}}dx$ be the complementary error function. Suppose that a random point field μ satisfies the following (A1)–(A5).
  • (A1) For each n, there exists n-correlation functions ${\rho ^{n}}$ of μ and ${\rho ^{n}}$ is locally bounded.
  • (A2) There exists the logarithmic derivative ${\mathsf{d}^{\mu }}$.
  • (A3) For each $k\in \{0\}\cup \mathbb{N}$, $({\mathcal{E}^{{\mu ^{[k]}}}},{\mathcal{D}_{\circ }^{{\mu ^{[k]}}}})$ is closable on ${L^{2}}({\mu ^{[k]}})$.
  • (A4) ${\mathrm{Cap}^{\mu }}({\mathsf{S}_{\mathrm{s}}^{c}})=0$.
  • (A5) There exists $T>0$ such that, for each $R>0$,
    \[ \underset{r\to \infty }{\liminf }\Big\{{\int _{|x|\le r+R}}{\rho ^{1}}(x)dx\Big\}\mathrm{Erf}\Big(\frac{r}{\sqrt{(r+R)T}}\Big)=0.\]
From (A3), let $({\mathcal{E}^{{\mu ^{[k]}}}},{\mathcal{D}^{{\mu ^{[k]}}}})$ be the closure of $({\mathcal{E}^{{\mu ^{[k]}}}},{\mathcal{D}_{\circ }^{{\mu ^{[k]}}}})$ on ${L^{2}}({\mu ^{[k]}})$. It is known that $({\mathcal{E}^{{\mu ^{[k]}}}},{\mathcal{D}^{{\mu ^{[k]}}}})$ is a quasi-regular Dirichlet form [24]. Then, ${\mathrm{Cap}^{\mu }}$ in (A4) denotes the one-capacity associated with $({\mathcal{E}^{{\mu ^{}}}},{\mathcal{D}^{{\mu ^{}}}},{L^{2}}({\mu ^{}}))$. Assumption (A4) means that particles never collide. Furthermore, (A5) implies that each tagged particle does not explode [24].
Recall that $W(A)=C([0,\infty ),A)$. Each $\mathsf{w}\in W({\mathsf{S}_{\mathrm{s},\mathrm{i}}})$ can be written as ${\mathsf{w}_{t}}={\textstyle\sum _{i}}{\delta _{{w_{t}^{i}}}}$, where ${w^{i}}$ is an S-valued continuous path defined on an interval ${I_{i}}$ of the form $[0,{b_{i}})$ or $({a_{i}},{b_{i}})$, where $0\le {a_{i}}<{b_{i}}\le \infty $. Taking maximal intervals of this form, we can choose $[0,{b_{i}})$ and $({a_{i}},{b_{i}})$ uniquely up to labelling. We remark that ${\lim \nolimits_{t\downarrow {a_{i}}}}|{w_{t}^{i}}|=\infty $ and ${\lim \nolimits_{t\uparrow {b_{i}}}}|{w_{t}^{i}}|=\infty $ for ${b_{i}}<\infty $ for all i. We call ${w^{i}}$ a tagged path of $\mathsf{w}$ and ${I_{i}}$ the defining interval of ${w^{i}}$. Let
\[ {W_{\mathrm{NE}}}({\mathsf{S}_{\mathrm{s},\mathrm{i}}})=\{\mathsf{w}\in W({\mathsf{S}_{\mathrm{s},\mathrm{i}}})\hspace{0.1667em};\hspace{0.1667em}{I_{i}}=[0,\infty )\hspace{2.5pt}\text{for all}\hspace{2.5pt}i\}.\]
It is said that the tagged path ${w^{i}}$ of $\mathsf{w}$ does not explode if ${b_{i}}=\infty $, and does not enter if ${I_{i}}=[0,{b_{i}})$, where ${b_{i}}$ is the right end of the defining interval of ${w^{i}}$. Thus, ${W_{\mathrm{NE}}}({\mathsf{S}_{\mathrm{s},\mathrm{i}}})$ is the set consisting of all nonexploding and nonentering paths in $W({\mathsf{S}_{\mathrm{s},\mathrm{i}}})$.
We can naturally lift each $\mathsf{w}={\{{\textstyle\sum _{i}}{\delta _{{w_{t}^{i}}}}\}_{t\in [0,\infty )}}\in {W_{\mathrm{NE}}}({\mathsf{S}_{\mathrm{s},\mathrm{i}}})$ to the labelled path
\[\begin{aligned}{}& \mathbf{w}={({w^{i}})_{i\in \mathbb{N}}}={\{{\mathbf{w}_{t}}\}_{t\in [0,\infty )}}={\{{({w_{t}^{i}})_{i\in \mathbb{N}}}\}_{t\in [0,\infty )}}\in W({S^{\mathbb{N}}})\end{aligned}\]
using a label $\mathfrak{l}={({\mathfrak{l}^{i}})_{i\in \mathbb{N}}}$. Indeed, for each $\mathsf{w}\in {W_{\mathrm{NE}}}({\mathsf{S}_{\mathrm{s},\mathrm{i}}})$, we can construct the labelled process $\mathbf{w}={\{{({w_{t}^{i}})_{i\in \mathbb{N}}}\}_{t\in [0,\infty )}}$ such that ${\mathbf{w}_{0}}=\mathfrak{l}({\mathsf{w}_{0}})$, because each tagged particle can carry the initial label i from the noncollision and nonexplosion properties of $\mathsf{w}$. We write this correspondence as
\[ {\mathfrak{l}_{\mathrm{path}}}(\mathsf{w})={({\mathfrak{l}_{\mathrm{path}}^{i}}(\mathsf{w}))_{i\in \mathbb{N}}}.\]
From (A1) and (A3), there exists a diffusion $\mathsf{X}$ on $\mathsf{S}$ associated with $({\mathcal{E}^{{\mu ^{}}}},{\mathcal{D}^{{\mu ^{}}}},{L^{2}}({\mu ^{}}))$ [22]. Then, we set $\mathbf{X}={\mathfrak{l}_{\mathrm{path}}}(\mathsf{X})$. The labelled process X gives a weak solution to an ISDE as follows.
Lemma 2.12 ([25, Theorem 26]).
Assume (A1)–(A5). Then, there exists a set $\mathsf{H}$ satisfying
\[ \mu (\mathsf{H})=1,\hspace{1em}\mathsf{H}\subset {\mathsf{S}_{\mathrm{s},\mathrm{i}}},\]
such that the following holds: for all $\mathbf{s}\in {\mathfrak{u}^{-1}}(\mathsf{H})$, there exists an ${S^{\mathbb{N}}}$-valued Brownian motion $\mathbf{B}={({B^{i}})_{i\in \mathbb{N}}}$ such that $(\mathbf{X},\mathbf{B})$ is a weak solution to
\[\begin{aligned}{}d{X_{t}^{i}}& =d{B_{t}^{i}}+\frac{1}{2}{\mathsf{d}^{\mu }}({X_{t}^{i}},{\mathsf{X}_{t}^{i\lozenge }})dt,\\ {} {\mathbf{X}_{0}}& =\mathbf{s},\end{aligned}\]
where ${\mathsf{X}_{t}^{i\lozenge }}={\textstyle\sum _{j\in \mathbb{N},j\ne i}}{\delta _{{X_{t}^{j}}}}$. Moreover, it holds that
\[ P({\mathsf{X}_{t}}\in \mathsf{H},0{\le ^{\forall }}t<\infty )=1.\]

2.4 Derivation of Theorem 2.1 from Theorem 2.6 and Theorem 2.11

We shall apply Lemma 2.12 for ${\mu _{\alpha }}$ to show Theorem 2.1.
Lemma 2.13.
For any $\alpha >1/2$, ${\mu _{\alpha }}$ satisfies (A1), (A4) and (A5).
Proof.
It is easy to see that ${K_{\alpha }}$ is bounded. Then, we have (A1) and (A5). Furthermore, because ${K_{\alpha }}$ is locally Lipschitz continuous, (A4) follows from [23, Theorem 2.1].  □
We derive Theorem 2.1 from Theorem 2.6 and Theorem 2.11.
Proof of Theorem 2.1.
Assumption (A2) holds from Theorem 2.11. Furthermore, we obtain (A3) from Theorem 2.6. Actually, (A3) for $k=0$ follows from the quasi-Gibbs property of ${\mu _{\alpha }}$ [26, Lemma 3.6]. For general $k\ge 1$, (A3) also follows from the quasi-Gibbs property of ${\mu _{\alpha }}$ in a similar manner.
Combining these with Lemma 2.13, we have checked (A1)–(A5) for ${\mu _{\alpha }}$. Therefore, Lemma 2.12 completes the proof of Theorem 2.1.  □
With the above argument, Theorem 2.1 is reduced to Theorem 2.6 and Theorem 2.11. The next section is devoted to preparing key tools to show Theorem 2.6 and Theorem 2.11.

3 Key tools for the proofs of main results

3.1 Determinantal kernels

For $\nu >-1$, let ${\{{L_{\nu }^{n}}(x)\}_{n\in \mathbb{N}}}$ be the classical Laguerre polynomials, that is,
\[ {L_{\nu }^{n}}(x)=\frac{{e^{x}}{x^{-\nu }}}{n!}\frac{{d^{n}}}{d{x^{n}}}({e^{-x}}{x^{n+\nu }}).\]
For simplicity, we write
\[ {\lambda _{\nu }^{n}}(x)=|x{|^{\nu }}{e^{-\frac{{x^{2}}}{2}}}{L_{\nu }^{n}}({x^{2}}).\]
The orthogonal polynomials ${\{{p_{\alpha }^{n}}\}_{n\in \mathbb{N}}}$, which appear in (1.3), are represented by the Laguerre polynomials. Accordingly, ${K_{G,\alpha }^{N}}$ can be rewritten in terms of ${\{{\lambda _{\nu }^{n}}\}_{n\in \mathbb{N}}}$ as follows.
Lemma 3.1 ([20]).
For any $k\in \{0\}\cup \mathbb{N}$,
\[\begin{aligned}{}& {K_{G,\alpha }^{N}}(x,y)\\ {} & =\left\{\begin{array}{l@{\hskip10.0pt}l}\frac{\Gamma (k+1)}{\Gamma (k+\alpha +\frac{1}{2})}\frac{x\sqrt{|{x^{-1}}y|}{\lambda _{\alpha +\frac{1}{2}}^{k}}(x){\lambda _{\alpha -\frac{1}{2}}^{k}}(y)-y\sqrt{|x{y^{-1}}|}{\lambda _{\alpha -\frac{1}{2}}^{k}}(x){\lambda _{\alpha +\frac{1}{2}}^{k}}(y)}{x-y}\hspace{1em}& \textit{for}\hspace{2.5pt}N=2k+1,\\ {} \frac{\Gamma (k+1)}{\Gamma (k+\alpha +\frac{1}{2})}\frac{x\sqrt{|{x^{-1}}y|}{\lambda _{\alpha +\frac{1}{2}}^{k-1}}(x){\lambda _{\alpha -\frac{1}{2}}^{k}}(y)-y\sqrt{|x{y^{-1}}|}{\lambda _{\alpha -\frac{1}{2}}^{k}}(x){\lambda _{\alpha +\frac{1}{2}}^{k-1}}(y)}{x-y}\hspace{1em}& \textit{for}\hspace{2.5pt}N=2k.\end{array}\right.\end{aligned}\]
This lemma makes an expression of the one-correlation function of ${\mu _{\mathrm{G},\alpha }^{N}}$. We note the symmetry ${\rho _{\mathrm{G},\alpha }^{N,1}}(-x)={\rho _{\mathrm{G},\alpha }^{N,1}}(x)$ and ${\rho _{\alpha }^{1}}(-x)={\rho _{\alpha }^{1}}(x)$. Then, for simplicity, we consider the one-correlation functions for nonnegative x in Lemma 3.2 and Lemma 3.3.
Lemma 3.2.
The following holds for $x\ge 0$: when $N=2k+1$,
\[\begin{aligned}{}& {\rho _{\mathrm{G},\alpha }^{N,1}}(x)={k^{\frac{1}{2}-\alpha }}\Big(2x\big\{{\lambda _{\alpha +\frac{1}{2}}^{k}}(x){\lambda _{\alpha +\frac{1}{2}}^{k-1}}(x)-{\lambda _{\alpha -\frac{1}{2}}^{k}}(x){\lambda _{\alpha +\frac{3}{2}}^{k-1}}(x)\big\}+{\lambda _{\alpha +\frac{1}{2}}^{k}}(x){\lambda _{\alpha -\frac{1}{2}}^{k}}(x)\Big)\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\times (1+\mathcal{O}({N^{-1}})),\end{aligned}\]
and when $N=2k$,
\[\begin{aligned}{}& {\rho _{\mathrm{G},\alpha }^{N,1}}(x)={k^{\frac{1}{2}-\alpha }}\Big(2x\big\{{\big({\lambda _{\alpha +\frac{1}{2}}^{k-1}}(x)\big)^{2}}-{\lambda _{\alpha -\frac{1}{2}}^{k}}(x){\lambda _{\alpha +\frac{3}{2}}^{k-2}}(x)\big\}+{\lambda _{\alpha +\frac{1}{2}}^{k-1}}(x){\lambda _{\alpha -\frac{1}{2}}^{k}}(x)\Big)\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\times (1+\mathcal{O}({N^{-1}}))\end{aligned}\]
as $N\to \infty $. Here, $\mathcal{O}({N^{-1}})$ terms are independent of $x\in \mathbb{R}$.
Proof.
Noting that $(d/dx){L_{\nu }^{k}}(x)=-{L_{\nu +1}^{k-1}}(x)$, we have the equality $(d/dx){L_{\nu }^{k}}({x^{2}})=-2x{L_{\nu +1}^{k-1}}({x^{2}})$. Combining this with Lemma 3.1 and the fact that
\[ \frac{\Gamma (k+a)}{\Gamma (k+b)}={k^{a-b}}(1+\mathcal{O}({k^{-1}}))\hspace{1em}\hspace{2.5pt}\text{as}\hspace{2.5pt}k\to \infty ,\]
we have the desired result.  □
Lemma 3.3.
For $x>0$, we have ${\rho _{\alpha }^{1}}(x)>0$.
Proof.
From (1.8) and the fact that ${J_{\nu +1}}(x)=2\nu {x^{-1}}{J_{\nu }}(x)-{J_{\nu -1}}(x)$, we see
\[\begin{aligned}{}{\rho _{\alpha }^{1}}(x)& =x\Big\{{\Big({J_{\alpha +\frac{1}{2}}}(\sqrt{2}x)-\frac{\alpha }{\sqrt{2}x}{J_{\alpha -\frac{1}{2}}}(\sqrt{2}x)\Big)^{2}}+\Big(1-\frac{{\alpha ^{2}}}{2{x^{2}}}\Big){J_{\alpha -\frac{1}{2}}^{2}}(\sqrt{2}x)\Big\}.\end{aligned}\]
Combining this with the fact that ${J_{\alpha +\frac{1}{2}}}(\sqrt{2}x)$ and ${J_{\alpha -\frac{1}{2}}}(\sqrt{2}x)$ have no common zero points, we see ${\rho _{\alpha }^{1}}(x)>0$ for $\sqrt{2}x>\alpha $. Recall that ${J_{\nu -1}}-{J_{\nu +1}}=2{J^{\prime }_{\nu }}$. Then, (1.8) yields that
\[\begin{aligned}{}{\rho _{\alpha }^{1}}(x)& =x\{{J^{\prime }_{\alpha +\frac{1}{2}}}(\sqrt{2}x){J_{\alpha -\frac{1}{2}}}(\sqrt{2}x)-{J^{\prime }_{\alpha -\frac{1}{2}}}(\sqrt{2}x){J_{\alpha +\frac{1}{2}}}(\sqrt{2}x)\}\\ {} & =\sqrt{2}\alpha {\int _{0}^{\sqrt{2}x}}\frac{{J_{\alpha +\frac{1}{2}}}(t){J_{\alpha -\frac{1}{2}}}(t)}{t}dt\\ {} & =\sqrt{2}\alpha {\int _{0}^{\sqrt{2}x}}\frac{{J_{\alpha +\frac{1}{2}}}(t)\{{J_{\alpha +\frac{3}{2}}}(t)+2{J^{\prime }_{\alpha +\frac{1}{2}}}(t)\}}{t}dt.\end{aligned}\]
We remark that the second equality is shown by differentiation of both sides with differential equation ${x^{2}}{J^{\prime\prime }_{\nu }}(x)+x{J^{\prime }_{\nu }}(x)+({x^{2}}-{\nu ^{2}}){J_{\nu }}(x)=0$. Because ${J_{\nu }}(z)>0$ and ${J^{\prime }_{\nu }}(z)>0$ for $z\in (0,\nu )$ (see [21, p. 246]), we have ${\rho _{\alpha }^{1}}(x)>0$ for $0<\sqrt{2}x\le \alpha $.  □

3.2 Asymptotic behaviour of the Laguerre polynomials

Because of Lemma 3.1, asymptotic analysis of the correlation functions boils down to that of Laguerre polynomials. We quote the Plancherel–Rotach type asymptotic results for the Laguerre polynomials by Erdélyi. We introduce the quantity
\[ {A_{\nu }^{n}}=n+\frac{\nu +1}{2}.\]
Lemma 3.4 ([5]).
Let $\nu >0$. Then, we have the following asymptotics:
  • (i) For any $0<\rho <\pi /2$,
    (3.1)
    \[\begin{aligned}{}& {L_{\nu }^{n}}(4{A_{\nu }^{n}}{\cos ^{2}}\rho )=\frac{{(-1)^{n}}{e^{2{A_{\nu }^{n}}{\cos ^{2}}\rho }}}{{(2\cos \rho )^{\nu }}{(\pi n\sin 2\rho )^{1/2}}}\\ {} & \hspace{2em}\times \Big\{\cos \left({A_{\nu }^{n}}(\sin 2\rho -2\rho )+\frac{\pi }{4}\right)+\mathcal{O}\Big(\frac{1}{n{\rho ^{3}}}\Big)+\mathcal{O}\Big(\frac{1}{n(\pi /2-\rho )}\Big)\Big\}\end{aligned}\]
    as $n{\rho ^{3}}\to \infty $ and $n(\pi /2-\rho )\to \infty $.
  • (ii) For any $0<\rho $,
    (3.2)
    \[ {L_{\nu }^{n}}(4{A_{\nu }^{n}}{\cosh ^{2}}\rho )=\frac{{(-1)^{n}}{e^{{A_{\nu }^{n}}(1+2\rho +{e^{-2\rho }})}}}{{(2\cosh \rho )^{\nu +1}}{(2\pi {A_{\nu }^{n}}\tanh \rho )^{1/2}}}\Big\{1+\mathcal{O}\Big(\frac{1+{\rho ^{3}}}{n{\rho ^{3}}}\Big)\Big\}\]
    as $n{\rho ^{3}}\to \infty $.
  • (iii) For any x satisfying $x-4{A_{\nu }^{n}}=o({n^{\frac{3}{5}}})$,
    (3.3)
    \[\begin{aligned}{}{L_{\nu }^{n}}(x)& =\frac{{(-1)^{n}}{e^{x/2}}}{{2^{\nu +1/3}}{n^{1/3}}}\Big\{\mathrm{Ai}\Big(\frac{x-4{A_{\nu }^{n}}}{{(16{A_{\nu }^{n}})^{1/3}}}\Big)\\ {} & +\tilde{\mathrm{Ai}}\Big(\frac{x-4{A_{\nu }^{n}}}{{(16{A_{\nu }^{n}})^{1/3}}}\Big)\Big[\mathcal{O}\Big(\frac{1}{n}\Big)+\mathcal{O}\Big(\frac{x-4{A_{\nu }^{n}}}{n}\Big)+\mathcal{O}\Big(\frac{{(x-4{A_{\nu }^{n}})^{5/2}}}{{n^{3/2}}}\Big)\Big]\Big\}\end{aligned}\]
    as $n\to \infty $. Here, $\mathrm{Ai}$ is the Airy function of the first kind, and $\tilde{\mathrm{Ai}}$ is defined as
    \[ \tilde{\mathrm{Ai}}(x)=\left\{\begin{array}{l@{\hskip10.0pt}l}\mathrm{Ai}(x)\hspace{1em}& \textit{for}\hspace{2.5pt}x\ge 0,\\ {} {\{|\mathrm{Ai}(x){|^{2}}+|\mathrm{Bi}(x){|^{2}}\}^{\frac{1}{2}}}\hspace{1em}& \hspace{5pt}\textit{for}\hspace{2.5pt}x<0,\end{array}\right.\]
    where $\mathrm{Bi}$ is the Airy function of the second kind.
Remark that several $\mathcal{O}$-terms in Lemma 3.4 contain two variables. For example, $f(n,\rho )\in \mathcal{O}({(n{\rho ^{3}})^{-1}})$ as $n{\rho ^{3}}\to \infty $ means that there exist positive constants C and L such that $|f(n,\rho )|\le C{(n{\rho ^{3}})^{-1}}$ holds for any $n,\rho $ satisfying $n{\rho ^{3}}>L$. Other terms have similar meaning.
With a computation similar to that used in [11], Lemma 3.4 yields the following asymptotics.
Lemma 3.5.
For $\nu >0$ and $m\in \{-2,-1,0\}$, the following hold:
  • (i) For any $0<\rho <\pi /2$, we have
    (3.4)
    \[\begin{aligned}{}{\lambda _{\nu }^{n+m}}(2{n^{\frac{1}{2}}}\cos \rho )& =\frac{{(-1)^{n+m}}{n^{\frac{\nu -1}{2}}}}{{(\pi \sin 2\rho )^{1/2}}}\Big\{\cos \left(n(\sin 2\rho -2\rho )-2\rho {A_{\nu }^{m}}+\frac{\pi }{4}\right)\\ {} & \hspace{2em}+\mathcal{O}\Big(\frac{1}{n{\rho ^{3}}}\Big)+\mathcal{O}\Big(\frac{1}{n(\pi /2-\rho )}\Big)\Big\}\end{aligned}\]
    as $n{\rho ^{3}}\to \infty $ and $n(\pi /2-\rho )\to \infty $.
  • (ii) For any $0<\rho $, we have
    (3.5)
    \[\begin{aligned}{}{\lambda _{\nu }^{n+m}}(2{n^{\frac{1}{2}}}\cosh \rho )& =\frac{{(-1)^{n+m}}{n^{\frac{\nu -1}{2}}}\exp \{n(2\rho -\sinh 2\rho )+2\rho {A_{\nu }^{m}}\}}{{2^{3/2}}\cosh \rho {(\pi \tanh \rho )^{1/2}}}\\ {} & \hspace{2em}\times \Big\{1+\mathcal{O}\Big(\frac{1+{\rho ^{3}}}{n{\rho ^{3}}}\Big)\Big\}\end{aligned}\]
    as $n{\rho ^{3}}\to \infty $.
Proof.
For n large enough, there exists $\theta =\theta (n,\rho )$ such that
(3.6)
\[ \cos \rho ={\Big(\frac{{A_{\nu }^{n+m}}}{n}\Big)^{\frac{1}{2}}}\cos (\rho +\theta ),\]
where ${A_{\nu }^{n+m}}=n+m+(\nu +1)/2$. Then, it is easy to see that
(3.7)
\[\begin{aligned}{}\theta =& \frac{\mathcal{O}({n^{-1}})}{\sin \rho },\end{aligned}\]
(3.8)
\[\begin{aligned}{}\cos (\rho +\theta )=& \Big\{1-\frac{{A_{\nu }^{m}}}{2n}+\mathcal{O}({n^{-2}})\Big\}\cos \rho \end{aligned}\]
as $n\to \infty $. Combining these with $\cos (\rho +\theta )=\cos \rho -\theta \sin \rho +\mathcal{O}({\theta ^{2}})$ as $\theta \to 0$, we obtain
(3.9)
\[\begin{aligned}{}\theta & =\frac{{A_{\nu }^{m}}}{2n}\frac{\cos \rho }{\sin \rho }+\frac{\mathcal{O}({n^{-2}})}{{\sin ^{3}}\rho }\end{aligned}\]
as $n\to \infty $. Furthermore, the Taylor expansion with (3.7) yields
(3.10)
\[ \sin 2(\rho +\theta )-2(\rho +\theta )=\sin 2\rho -2\rho -4\theta {\sin ^{2}}\rho +\frac{\mathcal{O}({n^{-2}})}{{\sin ^{2}}\rho }\]
as $n\to \infty $. Then, by straightforward computation with (3.10) and (3.9), we have
(3.11)
\[\begin{aligned}{}& {A_{\nu }^{n+m}}\{\sin 2(\rho +\theta )-2(\rho +\theta )\}=n(\sin 2\rho -2\rho )-2\rho {A_{\nu }^{m}}+\frac{\mathcal{O}({n^{-1}})}{{\sin ^{2}}\rho }\end{aligned}\]
as $n\to \infty $.
From (3.1) and (3.6), we see that
(3.12)
\[\begin{aligned}{}& {L_{\nu }^{n+m}}(4n{\cos ^{2}}\rho )={L_{\nu }^{n+m}}(4{A_{\nu }^{n+m}}{\cos ^{2}}(\rho +\theta ))\\ {} & =\frac{{(-1)^{n+m}}{e^{2n{\cos ^{2}}\rho }}}{{(2\cos \rho )^{\nu }}{(\pi n\sin 2(\rho +\theta ))^{1/2}}}\\ {} & \times \Big\{\cos \left({A_{\nu }^{n+m}}\{\sin 2(\rho +\theta )-2(\rho +\theta )\}+\frac{\pi }{4}\right)+\mathcal{O}\Big(\frac{1}{n{\rho ^{3}}}\Big)+\mathcal{O}\Big(\frac{1}{n(\pi /2-\rho )}\Big)\Big\}\end{aligned}\]
as $n{\rho ^{3}}\to \infty $ and $n(\pi /2-\rho )\to \infty $. Therefore, substituting (3.11) for (3.12), we finally derive (3.4).
To show (ii), let $\theta =\theta (n,\rho )$ be such that
(3.13)
\[ \cosh \rho ={\Big(\frac{{A_{\nu }^{n+m}}}{n}\Big)^{\frac{1}{2}}}\cosh (\rho +\theta ).\]
From this, we have that
(3.14)
\[\begin{aligned}{}\theta & =\frac{\mathcal{O}({n^{-1}})}{\sinh \rho },\end{aligned}\]
(3.15)
\[\begin{aligned}{}\cosh (\rho +\theta )& =\Big\{1-\frac{{A_{\nu }^{m}}}{2n}+\mathcal{O}({n^{-2}})\Big\}\cosh \rho \end{aligned}\]
as $n\to \infty $. Combining this with $\cosh (\rho +\theta )=\cosh \rho +\theta \sinh \rho +\mathcal{O}({\theta ^{2}})$ as $\theta \to 0$, we obtain
(3.16)
\[\begin{aligned}{}\theta & =-\frac{{A_{\nu }^{m}}}{2n}\frac{\cosh \rho }{\sinh \rho }+\frac{\mathcal{O}({n^{-2}})}{{\sinh ^{3}}\rho }\end{aligned}\]
as $n\to \infty $. By the Taylor expansion, (3.16) and (3.14) yield
\[\begin{aligned}{}{A_{\nu }^{n+m}}\{2\theta +{e^{-2\rho }}({e^{-2\theta }}-1)\}& ={A_{\nu }^{n+m}}\Big\{2\theta (1-{e^{-2\rho }})+\frac{\mathcal{O}({n^{-2}})}{{\sinh ^{2}}\rho }\Big\}\\ {} & =-{A_{\nu }^{m}}(1+{e^{-2\rho }})+\frac{\mathcal{O}({n^{-1}})}{{\sinh ^{2}}\rho }\end{aligned}\]
as $n\to \infty $. Thereby, combining this with the fact that $1+{e^{-2\rho }}=2{\cosh ^{2}}\rho -\sinh 2\rho $, we get
(3.17)
\[\begin{aligned}{}& {A_{\nu }^{n+m}}(1+2(\rho +\theta )+{e^{-2(\rho +\theta )}})\\ {} & ={A_{\nu }^{n+m}}(1+2\rho +{e^{-2\rho }})+{A_{\nu }^{n+m}}(2\theta +{e^{-2\rho }}({e^{-2\theta }}-1))\\ {} & =2n{\cosh ^{2}}\rho +n(2\rho -\sinh 2\rho )+2\rho {A_{\nu }^{m}}+\frac{\mathcal{O}({n^{-1}})}{{\sinh ^{2}}\rho }\end{aligned}\]
as $n\to \infty $.
From (3.2) and (3.13), we have
(3.18)
\[\begin{aligned}{}& {L_{\nu }^{n+m}}(4n{\cosh ^{2}}\rho )\\ {} & ={L_{\nu }^{n+m}}(4{A_{\nu }^{n+m}}{\cosh ^{2}}(\rho +\theta ))\\ {} & =\frac{{(-1)^{n+m}}\exp \{{A_{\nu }^{n+m}}(1+2(\rho +\theta )+{e^{-2(\rho +\theta )}})\}}{{(2\cosh (\rho +\theta ))^{\nu +1}}{(2\pi {A_{\nu }^{n+m}}\tanh (\rho +\theta ))^{1/2}}}\Big\{1+\mathcal{O}\Big(\frac{1+{\rho ^{3}}}{n{\rho ^{3}}}\Big)\Big\}\end{aligned}\]
as $n{\rho ^{3}}\to \infty $. Finally, substituting (3.17) for (3.18), we obtain (3.5).  □

3.3 Asymptotic estimates of determinantal kernels

Using the results in Section 3.2, we derive asymptotic estimates of determinantal kernels. For convenience, we set
\[ {M_{\alpha }^{N}}(x,y)=\frac{1}{\sqrt{N}}{K_{G,\alpha }^{N}}(\sqrt{N}x,\sqrt{N}y).\]
Note that (1.4) implies
(3.19)
\[ {M_{\alpha }^{N}}(x,y)={K_{\alpha }^{N}}(Nx,Ny).\]
Fix a constant γ satisfying
(3.20)
\[ -\frac{1}{2}<\gamma <-\frac{2}{5},\]
and we set
\[\begin{aligned}{}& {U_{1}^{N}}=[-\sqrt{2}+{N^{\gamma }},0)\cup (0,\sqrt{2}-{N^{\gamma }}],\\ {} & {U_{2}^{N}}=(-\infty ,-\sqrt{2}-{N^{\gamma }}]\cup [\sqrt{2}+{N^{\gamma }},\infty ),\\ {} & {T^{N}}=[-\sqrt{2}-{N^{\gamma }},-\sqrt{2}+{N^{\gamma }}]\cup [\sqrt{2}-{N^{\gamma }},\sqrt{2}+{N^{\gamma }}].\end{aligned}\]
Remark that $\mathbb{R}\setminus \{0\}={U_{1}^{N}}\cup {U_{2}^{N}}\cup {T^{N}}$. Furthermore, we make the definition
\[ {U^{N}}={U_{1}^{N}}\cup {U_{2}^{N}}.\]
Hereafter, we always suppose $\alpha >1/2$.
Lemma 3.6.
  • (i) For any $x\in {U_{1}^{N}}$, we have
    (3.21)
    \[ {M_{\alpha }^{N}}(x,x)=\frac{\sqrt{2-{x^{2}}}}{\pi }+\mathcal{O}\Big(\frac{1}{{N^{1+2\gamma }}}\Big)+\mathcal{O}\Big(\frac{1}{N|x|}\Big)\]
    as $N\to \infty $ and $N|x|\to \infty $.
  • (ii) There exist positive constants ${c_{3}}$ and L such that, for any $x\in {U_{2}^{N}}$ and $N\in \mathbb{N}$ satisfying $N{({x^{2}}-2)^{2}}>L$, we have
    (3.22)
    \[ {M_{\alpha }^{N}}(x,x)\le \frac{{c_{\textit{3}}}}{N{({x^{2}}-2)^{2}}}.\]
Proof.
We prove for the case $N=2k$ for $k\in \mathbb{N}$; the case of odd N can be proved in the same way. It is sufficient to see for $x>0$. Put $\sqrt{N}x=2\sqrt{k}\cos \rho $. Considering $|2-{x^{2}}|\ge {N^{\gamma }}$ for $x\in {U_{1}^{N}}$, we see ${\rho ^{-1}}=\mathcal{O}({N^{-\gamma /2}})$ as $N\to \infty $. Set $\tilde{\rho }=\sin 2\rho -2\rho $ for simplicity, and (3.4) becomes
\[\begin{aligned}{}& {\lambda _{\nu }^{k+m}}(\sqrt{N}x)=\frac{{(-1)^{k+m}}{k^{\frac{\nu -1}{2}}}}{{(\pi \sin 2\rho )^{1/2}}}\Big\{\cos \Big(k\tilde{\rho }-2\rho {A_{\nu }^{m}}+\frac{\pi }{4}\Big)\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}+\mathcal{O}\Big(\frac{1}{{N^{1+3\gamma /2}}}\Big)+\mathcal{O}\Big(\frac{1}{N|x|}\Big)\Big\}\end{aligned}\]
as $N\to \infty $ and $N|x|\to \infty $. Then, this yields
(3.23)
\[\begin{aligned}{}& 2\sqrt{N}x\big\{{({\lambda _{\alpha +\frac{1}{2}}^{k-1}}(\sqrt{N}x))^{2}}-{\lambda _{\alpha -\frac{1}{2}}^{k}}(\sqrt{N}x){\lambda _{\alpha +\frac{3}{2}}^{k-2}}(\sqrt{N}x)\big\}\\ {} & \hspace{2em}=\frac{2{k^{\alpha }}}{\pi }\Big\{\sin \rho +\mathcal{O}\Big(\frac{1}{{N^{1+2\gamma }}}\Big)+\mathcal{O}\Big(\frac{1}{N|x|}\Big)\Big\}\end{aligned}\]
as $N\to \infty $ and $N|x|\to \infty $. Here, we used the fact that ${\cos ^{2}}A-\cos (A+B)\cos (A-B)={\sin ^{2}}B$ with $A=k\tilde{\rho }-(\alpha -\frac{1}{2})\rho +\frac{\pi }{4}$ and $B=\rho $. Furthermore, we have
(3.24)
\[\begin{aligned}{}& {\lambda _{\alpha +\frac{1}{2}}^{k-1}}(\sqrt{N}x){\lambda _{\alpha -\frac{1}{2}}^{k}}(\sqrt{N}x)\\ {} & \hspace{2em}=\frac{{k^{\alpha -1}}}{\pi \sin 2\rho }\Big\{\mathcal{O}(1)+\mathcal{O}\Big(\frac{1}{{N^{1+3\gamma /2}}}\Big)+\mathcal{O}\Big(\frac{1}{N|x|}\Big)\Big\}\end{aligned}\]
as $N\to \infty $ and $N|x|\to \infty $. Then, (3.23) and (3.24) with Lemma 3.2 yield (3.21).
Next, we shall show (ii). Using (3.5) for $\sqrt{N}x=2\sqrt{k}\cosh \rho $, we get
\[\begin{aligned}{}& 2\sqrt{N}x\{{({\lambda _{\alpha +\frac{1}{2}}^{k-1}}(\sqrt{N}x))^{2}}-{\lambda _{\alpha -\frac{1}{2}}^{k}}(\sqrt{N}x){\lambda _{\alpha +\frac{3}{2}}^{k-2}}(\sqrt{N}x)\}\\ {} & \hspace{2em}\hspace{2em}=\frac{{k^{\alpha }}{e^{2k(2\rho -\sinh 2\rho )+(2\alpha -1)\rho }}}{2\pi \sinh \rho }\times \mathcal{O}\Big(\frac{1+{\rho ^{3}}}{N{\rho ^{3}}}\Big)\end{aligned}\]
and
\[\begin{aligned}{}{\lambda _{\alpha +\frac{1}{2}}^{k-1}}(\sqrt{N}x){\lambda _{\alpha -\frac{1}{2}}^{k}}(\sqrt{N}x)& =\frac{-{k^{\alpha -1}}{e^{2k(2\rho -\sinh 2\rho )+2\alpha \rho }}}{{2^{3}}\pi \sinh \rho \cosh \rho }\Big\{1+\mathcal{O}\Big(\frac{1+{\rho ^{3}}}{N{\rho ^{3}}}\Big)\Big\}\end{aligned}\]
as $N{\rho ^{3}}\to \infty $. Then, combining these with Lemma 3.2, we obtain (3.22).  □
Lemma 3.7.
  • (i) There exist a positive constant ${c_{4}}$ and $f(N,x)\in \mathcal{O}({(N|x|)^{-1}})$ as $N|x|\to \infty $ such that, for all $N\in \mathbb{N}$ and $x\in {U^{N}}$,
    (3.25)
    \[ {M_{\alpha }^{N}}(x,x)\le {c_{\textit{4}}}\big(1+f(N,x)\big)\]
  • (ii) There exist a positive constant ${c_{5}}$ such that, for all $N\in \mathbb{N}$ and $x\in {T^{N}}$,
    (3.26)
    \[ {M_{\alpha }^{N}}(x,x)\le {c_{\textit{5}}}{N^{\frac{1}{3}}}.\]
  • (iii) There exists a positive constant ${c_{6}}$ and $f(N,x)\in \mathcal{O}({(N|x|)^{-1}})$ as $N|x|\to \infty $ such that, for all $N\in \mathbb{N}$ and $x,y\in {U^{N}}$,
    (3.27)
    \[ |{M_{\alpha }^{N}}(x,y)|\le \frac{{c_{\textit{6}}}(1+f(N,x))(1+f(N,y))}{N|x-y||2-{x^{2}}{|^{1/4}}|2-{y^{2}}{|^{1/4}}}.\]
Proof.
Inequality (3.25) follows from (3.21) and (3.22). Let $N=2k$ (the case $N=2k+1$ is proven in the same manner). For $\nu >0$, we see $N{x^{2}}-4{A_{\nu }^{k}}=o({k^{\frac{3}{5}}})$ for $x\in {T^{N}}$ since $\gamma <-2/5$. Therefore, (3.3) yields
\[\begin{aligned}{}|{\lambda _{\nu }^{k}}(\sqrt{N}x)|& ={N^{\nu /2}}|x{|^{\nu }}{e^{-\frac{N{x^{2}}}{2}}}|{L_{\nu }^{k}}(N{x^{2}})|\le C{N^{\frac{\nu }{2}-\frac{1}{3}}}\end{aligned}\]
for some constant C. Here, we used the fact that $\mathrm{Ai}$ is bounded on $\mathbb{R}$ and $\mathrm{Bi}$ is bounded on $(-\infty ,0)$ (see [4, Section 9], for example). This with Lemma 3.2 yields (3.26). Furthermore, combining Lemma 3.1 with (3.4) and (3.5), we get (3.27).  □

4 The logarithmic derivative

4.1 Finite particle approximation of the logarithmic derivative

The logarithmic derivative ${\mathsf{d}^{\mu }}$ of μ can be approximated by that of finite particle systems. Let ${\{{\mu ^{N}}\}_{N\in \mathbb{N}}}$ be a sequence of random point fields such that ${\mu ^{N}}(\mathsf{s}(S)=N)=1$. Moreover, ${\rho ^{n}}$ and ${\rho ^{N,n}}$ stand for the n-correlation function of μ and ${\mu ^{N}}$, respectively. Then, we assume the following:
  • (B1) For each $R\in \mathbb{N}$,
    (4.1)
    \[\begin{aligned}{}& \underset{N\to \infty }{\lim }{\rho ^{N,n}}={\rho ^{n}}\hspace{2.5pt}\text{uniformly on}\hspace{2.5pt}{S_{R}^{n}}\hspace{2.5pt}\text{for each}\hspace{2.5pt}n\in \mathbb{N},\end{aligned}\]
    (4.2)
    \[\begin{aligned}{}& \underset{N\in \mathbb{N}}{\sup }\underset{{\mathbf{x}_{n}}\in {S_{R}^{n}}}{\sup }{\rho ^{N,n}}({\mathbf{x}_{n}})\le {c_{\text{7}}^{n}}{n^{{c_{\text{8}}}n}}\hspace{2.5pt}\text{for any}\hspace{2.5pt}n\in \mathbb{N},\end{aligned}\]
    where $0<{c_{7}}(R)<\infty $ and $0<{c_{8}}(R)<1$ are constants independent of n.
Note that (B1) implies ${\lim \nolimits_{N\to \infty }}{\mu ^{N}}=\mu $ weakly.
Let $u,{u^{N}}:S\to \mathbb{R}$ and $g:{S^{2}}\to \mathbb{R}$ be measurable functions. For $r>0$, we set
\[ {\mathsf{g}_{r}}(x,\mathsf{y})=\sum \limits_{i}{\chi _{r}}(x-{y_{i}})g(x,{y_{i}}),\hspace{2em}{\mathsf{w}_{r}}(x,\mathsf{y})=\sum \limits_{i}(1-{\chi _{r}}(x-{y_{i}}))g(x,{y_{i}}),\]
where $\mathsf{y}={\textstyle\sum _{i}}{\delta _{{y_{i}}}}$ and ${\chi _{r}}\in {C_{0}^{\infty }}(S)$ is a cut-off function such that $0\le {\chi _{r}}\le 1$, ${\chi _{r}}(x)=0$, for $|x|\ge r+1$ and ${\chi _{r}}(x)=1$ for $|x|\le r$. We set ${\tilde{S}_{R}}=\{x\hspace{0.1667em};\hspace{0.1667em}{R^{-1}}\le |x|\le R\}$. Following [12], we make assumption (B2).
  • (B2) For each N, ${\mu ^{N}}$ has the logarithmic derivative ${\mathsf{d}^{{\mu ^{N}}}}$ such that
    \[\begin{aligned}{}& {\mathsf{d}^{{\mu ^{N}}}}(x,\mathsf{y})={u^{N}}(x)+{\mathsf{g}_{r}}(x,\mathsf{y})+{\mathsf{w}_{r}}(x,\mathsf{y}).\end{aligned}\]
    Furthermore, ${u^{N}}$, ${\mathsf{g}_{r}}$, and ${\mathsf{w}_{r}}$ satisfy the following (i)–(iii) for some $\hat{p}>1$:
    • (i) It holds that ${u^{N}}\in {C^{1}}(S\setminus \{0\})$. Furthermore, ${u^{N}}$ and $\nabla {u^{N}}$ converge uniformly to u and $\nabla u$ on each compact set in $S\setminus \{0\}$, respectively.
    • (ii) It holds that $g\in {C^{1}}({S^{2}}\cap \{x\ne y\})$. In addition, for each $R\in \mathbb{N}$,
      (4.3)
      \[\begin{aligned}{}& \underset{\mathsf{p}\to \infty }{\lim }\underset{N\to \infty }{\limsup }{\int _{x\in {\tilde{S}_{R}},|x-y|\le {2^{-\mathsf{p}}}}}{\chi _{r}}(x-y)|g(x,y){|^{\hat{p}}}\hspace{0.1667em}{\rho _{x}^{N,1}}(y)dxdy=0,\end{aligned}\]
      where ${\rho _{x}^{N,1}}$ is the one-correlation function of the reduced Palm measure ${\mu _{x}^{N}}$.
    • (iii) For each $R\in \mathbb{N}$,
      (4.4)
      \[ \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }{\int _{{\tilde{S}_{R}}\times \mathsf{S}}}|{\mathsf{w}_{r}}(x,\mathsf{y}){|^{\hat{p}}}d{\mu ^{N,[1]}}=0.\]
Then, we obtain an explicit expression of the logarithmic derivative of μ by finite particle approximation.
Lemma 4.1 ([25, Theorem 45]).
Assume (B1) and (B2). Then, the logarithmic derivative ${\mathsf{d}^{\mu }}$ of μ exists in ${L_{\mathrm{loc}}^{p}}({\mu ^{[1]}})$ for $1\le p<\hat{p}$, and it is represented as
\[ {\mathsf{d}^{\mu }}(x,\mathsf{y})=u(x)+\underset{r\to \infty }{\lim }{\mathsf{g}_{r}}(x,\mathsf{y}).\]
Here, the convergence ${\lim \nolimits_{r\to \infty }}{\mathsf{g}_{r}}$ takes place in ${L_{\mathrm{loc}}^{p}}({\mu ^{[1]}})$.

4.2 The logarithmic derivative of the generalised sine random point field

Let ${\mu _{\alpha }^{N}}$ be the determinantal random point field whose kernel is ${K_{\alpha }^{N}}$ as in (1.4). In other words, the (labelled) density of ${\mu _{\alpha }^{N}}$ is given by
(4.5)
\[ {m_{\alpha }^{N}}(d{\mathbf{x}_{N}})=\frac{1}{\mathcal{Z}}\prod \limits_{1\le i<j\le N}|{x_{i}}-{x_{j}}{|^{2}}{\prod \limits_{k=1}^{N}}|{x_{k}}{|^{2\alpha }}{e^{-\frac{{x_{k}^{2}}}{N}}}d{\mathbf{x}_{N}}.\]
Let ${\rho _{\alpha }^{N,n}}$ be the n-correlation function of ${\mu _{\alpha }^{N}}$. Furthermore, let ${\mu _{\alpha ,x}^{N}}$ be the reduced Palm measure of ${\mu _{\alpha }^{N}}$ conditioned at x, and ${\rho _{\alpha ,x}^{N,n}}$ denotes its n-correlation function. We shall apply the general result in Section 4.1 taking $\mu ={\mu _{\alpha }}$ and ${\mu ^{N}}={\mu _{\alpha }^{N}}$.
From (4.5), the logarithmic derivative of ${\mu _{\alpha }^{N}}$ is given by
\[ {\mathsf{d}^{{\mu _{\alpha }^{N}}}}(x,\mathsf{y})=\sum \limits_{j}\frac{2}{x-{y_{j}}}+\frac{2\alpha }{x}-\frac{2x}{N}\]
for $\mathsf{y}={\textstyle\sum _{i=1}^{N-1}}{\delta _{{y_{j}}}}$. Then, we shall show (B2) under the setting of
(4.6)
\[ {u^{N}}(x)=\frac{2\alpha }{x}-\frac{2x}{N},\hspace{1em}g(x,y)=\frac{2}{x-y}.\]
The most tough assumption in Lemma 4.1 is (4.4). We introduce a sufficient condition for (4.4). We set $||\cdot |{|_{R}}={\sup _{x\in {\tilde{S}_{R}}}}|\cdot |$.
Lemma 4.2.
We set ${S_{x,r}}=\{y\in S\hspace{0.1667em};\hspace{0.1667em}r\le |x-y|<\infty \}$. The following conditions (4.7)–(4.10) imply (4.4) with $\hat{p}=2$:
(4.7)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{S_{x,r}}}}\frac{{\rho _{\alpha }^{N,1}}(y)}{x-y}dy{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.8)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{S_{x,r}}}}\frac{{\rho _{\alpha ,x}^{N,1}}(y)-{\rho _{\alpha }^{N,1}}(y)}{x-y}dy{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.9)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{S_{x,r}}}}\frac{{\rho _{\alpha }^{N,1}}(y)}{{(x-y)^{2}}}dy\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}+{\int _{{({S_{x,r}})^{2}}}}\frac{{\rho _{\alpha }^{N,2}}(y,z)-{\rho _{\alpha }^{N,1}}(y){\rho _{\alpha }^{N,1}}(z)}{(x-y)(x-z)}dydz{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.10)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{S_{x,r}}}}\frac{{\rho _{\alpha ,x}^{N,1}}(y)-{\rho _{\alpha }^{N,1}}(y)}{{(x-y)^{2}}}dy+\\ {} & {\int _{{({S_{x,r}})^{2}}}}\frac{{\rho _{\alpha ,x}^{N,2}}(y,z)-{\rho _{\alpha ,x}^{N,1}}(y){\rho _{\alpha ,x}^{N,1}}(z)-{\rho _{\alpha }^{N,2}}(y,z)+{\rho _{\alpha }^{N,1}}(y){\rho _{\alpha }^{N,1}}(z)}{(x-y)(x-z)}dydz{\bigg|\bigg|_{R}}=0.\end{aligned}\]
Proof.
This lemma follows from [25, Lemma 52, Lemma 53]. Remark that, in [25], the sup norm $||\cdot |{|_{R}}$ is defined as ${\sup _{x\in {S_{R}}}}||\cdot ||$. Here, it is taken over ${\tilde{S}_{R}}$ instead of ${S_{R}}$, considering the space of test functions of the logarithmic derivative in Definition 2.9.  □

4.3 Proof of Theorem 2.11

We begin by rewriting the conditions in Lemma 4.2 in terms of determinantal kernels. We recall that reduced Palm measures of determinantal random point fields are also determinantal. From [30], ${\mu _{\alpha ,x}^{N}}$ is determinantal, and its kernel ${K_{\alpha ,x}^{N}}$ is given by
(4.11)
\[ {K_{\alpha ,x}^{N}}(y,z)={K_{\alpha }^{N}}(y,z)-\frac{{K_{\alpha }^{N}}(x,y){K_{\alpha }^{N}}(x,z)}{{K_{\alpha }^{N}}(x,x)}.\]
With the aid of (4.11), we see that
(4.12)
\[\begin{aligned}{}& {\rho _{\alpha }^{N,1}}(y)={K_{\alpha }^{N}}(y,y)\end{aligned}\]
(4.13)
\[\begin{aligned}{}& {\rho _{\alpha ,x}^{N,1}}(y)-{\rho _{\alpha }^{N,1}}(y)=\frac{{K_{\alpha }^{N}}{(x,y)^{2}}}{{K_{\alpha }^{N}}(x,x)}\end{aligned}\]
(4.14)
\[\begin{aligned}{}& {\rho _{\alpha }^{N,2}}(y,z)-{\rho _{\alpha }^{N,1}}(y){\rho _{\alpha }^{N,1}}(z)=-{K_{\alpha }^{N}}{(y,z)^{2}}\end{aligned}\]
(4.15)
\[\begin{aligned}{}& {\rho _{\alpha ,x}^{N,2}}(y,z)-{\rho _{\alpha ,x}^{N,1}}(y){\rho _{\alpha ,x}^{N,1}}(z)-\{{\rho _{\alpha }^{N,2}}(y,z)-{\rho _{\alpha }^{N,1}}(y){\rho _{\alpha }^{N,1}}(z)\}\\ {} & \hspace{2em}=2\frac{{K_{\alpha }^{N}}(y,z){K_{\alpha }^{N}}(x,y){K_{\alpha }^{N}}(x,z)}{{K_{\alpha }^{N}}(x,x)}-\frac{{K_{\alpha }^{N}}{(x,y)^{2}}{K_{\alpha }^{N}}{(x,z)^{2}}}{{K_{\alpha }^{N}}{(x,x)^{2}}}.\end{aligned}\]
Thereby, (4.7)–(4.10) in Lemma 4.2 are rewritten as follows.
Lemma 4.3.
To simplify notation, we set
\[ {x_{N}}=\frac{x}{N},\hspace{2em}{T_{x,r}^{N}}=\{y\in \mathbb{R}\hspace{0.1667em};\hspace{0.1667em}\frac{r}{N}\le |{x_{N}}-y|<\infty \}.\]
Then, (4.7)–(4.10) are equivalent to (4.16)–(4.19), respectively.
(4.16)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{{x_{N}}-y}dy{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.17)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{1}{{x_{N}}-y}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.18)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{N|{x_{N}}-y{|^{2}}}dy-{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{({x_{N}}-y)({x_{N}}-z)}dydz{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.19)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{1}{N|{x_{N}}-y{|^{2}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}+{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{1}{({x_{N}}-y)({x_{N}}-z)}\Big(\frac{2{M_{\alpha }^{N}}(y,z){M_{\alpha }^{N}}({x_{N}},y){M_{\alpha }^{N}}({x_{N}},z)}{{M_{\alpha }^{N}}({x_{N}},{x_{N}})}\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}\hspace{2em}-\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}{M_{\alpha }^{N}}{({x_{N}},z)^{2}}}{{M_{\alpha }^{N}}{({x_{N}},{x_{N}})^{2}}}\Big)dydz{\bigg|\bigg|_{R}}=0.\end{aligned}\]
Proof.
Recall that ${M_{\alpha }^{N}}(x,y)={K_{\alpha }^{N}}(Nx,Ny)$ in (3.19). Then, by (4.12) and the change of variables $y\mapsto Ny$, we have that
\[ {\int _{{S_{x,r}}}}\frac{{\rho _{\alpha }^{N,1}}(y)}{x-y}dy={\int _{{T_{x,r}^{N}}}}\frac{{K_{\alpha }^{N}}(Ny,Ny)}{x-Ny}Ndy={\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{{x_{N}}-y}dy.\]
Therefore, we obtain the equivalence between (4.7) and (4.16). By the same argument with (4.13)–(4.15), it holds that (4.8)–(4.10) are equivalent to (4.17)–(4.19), respectively.  □
We shall prove (4.16)–(4.19) in order, in the rest of this section.
Lemma 4.4.
For $0<q\le 1$, we have
(4.20)
\[ \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T^{N}}\cup {U_{2}^{N}}}}\frac{{M_{\alpha }^{N}}{(y,y)^{q}}}{|{x_{N}}-y|}dy{\bigg|\bigg|_{R}}=0.\]
Proof.
From (3.26), we have
(4.21)
\[\begin{aligned}{}& \bigg|\bigg|{\int _{{T^{N}}}}\frac{{M_{\alpha }^{N}}{(y,y)^{q}}}{|{x_{N}}-y|}dy{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{{T^{N}}}}\frac{{c_{\text{4}}^{q}}{N^{q/3}}}{|{x_{N}}-y|}dy{\bigg|\bigg|_{R}}\\ {} & ={c_{\text{4}}^{q}}{N^{q/3}}\big(\log |{x_{N}}-(\sqrt{2}-{N^{\gamma }})|-\log |{x_{N}}-(\sqrt{2}+{N^{\gamma }})|\\ {} & \hspace{2em}\hspace{2em}\hspace{2em}+\log |{x_{N}}-(-\sqrt{2}-{N^{\gamma }})|-\log |{x_{N}}-(-\sqrt{2}+{N^{\gamma }})|\big)\\ {} & =\mathcal{O}({N^{\frac{q}{3}+\gamma }})\to 0\hspace{2em}\hspace{2.5pt}\text{as}\hspace{2.5pt}N\to \infty .\end{aligned}\]
Here, we used the fact that $q\le 1$ and $\gamma <-2/5$ in (3.20) in the last line. Moreover, from (3.22) and $-1/2<\gamma $, we have
(4.22)
\[\begin{aligned}{}\bigg|\bigg|{\int _{{U_{2}^{N}}}}\frac{{M_{\alpha }^{N}}{(y,y)^{q}}}{|{x_{N}}-y|}dy{\bigg|\bigg|_{R}}& \le \bigg|\bigg|{\int _{{U_{2}^{N}}}}\frac{{c_{\text{3}}^{q}}}{|{x_{N}}-y|{N^{q}}|{y^{2}}-2{|^{2q}}}dy{\bigg|\bigg|_{R}}\\ {} & \to 0\hspace{2em}\hspace{2.5pt}\text{as}\hspace{2.5pt}N\to \infty .\end{aligned}\]
Hence, from (4.21) and (4.22), we have (4.20).  □
Lemma 4.5.
Equation (4.16) holds.
Proof.
Our proof starts with the observation that (3.21) yields
(4.23)
\[ \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{{x_{N}}-y}dy{\bigg|\bigg|_{R}}=0.\]
Indeed, we see that
(4.24)
\[\begin{aligned}{}\underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{1}{{x_{N}}-y}\frac{\sqrt{2-{y^{2}}}}{\pi }dy{\bigg|\bigg|_{R}}& =\bigg|\mathrm{P}.\mathrm{V}.{\int _{-\sqrt{2}}^{\sqrt{2}}}\frac{1}{y}\frac{\sqrt{2-{y^{2}}}}{\pi }dy\bigg|\\ {} & =0,\end{aligned}\]
and, from $-1/2<\gamma $,
(4.25)
\[ \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{1}{|{x_{N}}-y|}\frac{1}{{N^{1+2\gamma }}}dy{\bigg|\bigg|_{R}}=0.\]
Moreover,
(4.26)
\[ \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{1}{|{x_{N}}-y|}\frac{1}{N|y|}dy{\bigg|\bigg|_{R}}=\bigg|\bigg|\frac{1}{x}\log \frac{x+r}{r-x}{\bigg|\bigg|_{R}}.\]
From (4.24), (4.25), and (4.26) with (3.21), we have (4.23). Therefore, (4.16) follows from (4.20) with $q=1$ and (4.23).  □
Applying the Schwartz inequality to (1.3), we have that ${K_{G,\alpha }^{N}}{(x,y)^{2}}\le {K_{G,\alpha }^{N}}(x,x){K_{G,\alpha }^{N}}(y,y)$, which yields
(4.27)
\[ {M_{\alpha }^{N}}{(x,y)^{2}}\le {M_{\alpha }^{N}}(x,x){M_{\alpha }^{N}}(y,y).\]
From ${\tilde{S}_{R}}=\{{R^{-1}}\le |x|\le R\}$, Lemma 3.3 yields
\[ 0<\underset{x\in {\tilde{S}_{R}}}{\inf }{K_{\alpha }}(x,x).\]
Then, since ${M_{\alpha }^{N}}({x_{N}},{x_{N}})={K_{\alpha }^{N}}(x,x)$ and ${\lim \nolimits_{N\to \infty }}{K_{\alpha }^{N}}(x,x)={K_{\alpha }}(x,x)$ uniformly for $x\in {\tilde{S}_{R}}$, we see that ${c_{9}}$ below is finite for each $R\in \mathbb{N}$:
(4.28)
\[ {c_{\text{9}}}=\underset{N\in \mathbb{N}}{\sup }\underset{x\in {\tilde{S}_{R}}}{\sup }\frac{1}{{M_{\alpha }^{N}}({x_{N}},{x_{N}})}<\infty .\]
Let $f(N,y)\in \mathcal{O}({(N|y|)^{-1}})$ as $N|y|\to \infty $. Then, there exist positive constants C and L such that $|f(N,y)|\le C{(N|y|)^{-1}}$ for $N,y$ satisfying $N|y|>L$. From this, for sufficiently large r, we have
\[ \underset{y\in {T_{x,r}^{N}},N\in \mathbb{N},x\in {\tilde{S}_{R}}}{\sup }|f(N,y)|<\infty .\]
Thus, from (3.25) and (3.27), there exists a positive constant ${c_{10}}$ independent of N and r such that
(4.29)
\[\begin{array}{r@{\hskip0pt}l@{\hskip0pt}r@{\hskip0pt}l}\displaystyle {M_{\alpha }^{N}}(y,y)& \displaystyle \le {c_{\text{10}}}& & \displaystyle \hspace{2.5pt}\text{for any}\hspace{2.5pt}y\in {T_{x,r}^{N}}\cap {U^{N}},\end{array}\]
(4.30)
\[\begin{array}{r@{\hskip0pt}l@{\hskip0pt}r@{\hskip0pt}l}\displaystyle |{M_{\alpha }^{N}}({x_{N}},y)|& \displaystyle \le \frac{{c_{\text{10}}}}{N|{x_{N}}-y||2-{y^{2}}{|^{1/4}}}& & \displaystyle \hspace{2.5pt}\text{for any}\hspace{2.5pt}x\in {\tilde{S}_{R}},\hspace{0.1667em}\hspace{0.1667em}y\in {T_{x,r}^{N}}\cap {U^{N}},\end{array}\]
(4.31)
\[\begin{array}{r@{\hskip0pt}l@{\hskip0pt}r@{\hskip0pt}l}\displaystyle |{M_{\alpha }^{N}}(y,z)|& \displaystyle \le \frac{{c_{\text{10}}}}{N|y-z||2-{y^{2}}{|^{1/4}}|2-{z^{2}}{|^{1/4}}}\hspace{1em}& & \displaystyle \hspace{2.5pt}\text{for any}\hspace{2.5pt}y,z\in {T_{x,r}^{N}}\cap {U^{N}}.\end{array}\]
Lemma 4.6.
The following (4.32) and (4.33) hold.
(4.32)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.33)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}({x_{N}},y)}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}=0.\end{aligned}\]
In particular, (4.17) holds.
Proof.
Equation (4.17) follows from (4.32) immediately. Then, we first prove (4.32). From (4.27) and (4.20) with $q=1$, we deduce that
(4.34)
\[\begin{aligned}{}& \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T^{N}}\cup {U_{2}^{N}}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}\\ {} & \hspace{2em}\le \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T^{N}}\cup {U_{2}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{|{x_{N}}-y|}dy{\bigg|\bigg|_{R}}=0.\end{aligned}\]
From (4.28) and (4.30), we have, for large r,
(4.35)
\[\begin{aligned}{}& \underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}\\ {} & \le \underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{{c_{\text{9}}}{c_{\text{10}}^{2}}}{{N^{2}}|{x_{N}}-y{|^{3}}{(2-{y^{2}})^{1/2}}}dy{\bigg|\bigg|_{R}}\le \frac{{c_{\text{9}}}{c_{\text{10}}^{2}}}{{r^{2}}}.\end{aligned}\]
Here, we used the fact that $|2-{y^{2}}|\ge {N^{\gamma }}$ for $y\in {U_{1}^{N}}$. Combining (4.34) with (4.35), we get (4.32).
By the same argument as above, (4.28) and (4.30) yield that
(4.36)
\[\begin{aligned}{}& \underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{{M_{\alpha }^{N}}({x_{N}},y)}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}\\ {} & \le \underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U_{1}^{N}}}}\frac{{c_{\text{9}}}{c_{\text{10}}}}{N|{x_{N}}-y|{(2-{y^{2}})^{1/4}}}dy{\bigg|\bigg|_{R}}\le \frac{2{c_{\text{9}}}{c_{\text{10}}}}{r}.\end{aligned}\]
Furthermore, from (4.27), (4.28), and (4.20) with $q=1/2$, we see that
(4.37)
\[\begin{aligned}{}& \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T^{N}}\cup {U_{2}^{N}}}}\frac{{M_{\alpha }^{N}}({x_{N}},y)}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}\\ {} & \hspace{2em}\hspace{2em}\le \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T^{N}}\cup {U_{2}^{N}}}}\frac{{c_{\text{9}}^{\frac{1}{2}}}{M_{\alpha }^{N}}{(y,y)^{\frac{1}{2}}}}{|{x_{N}}-y|}dy{\bigg|\bigg|_{R}}=0.\end{aligned}\]
Hence, (4.37) and (4.36) yield (4.33).  □
Lemma 4.7.
The following (4.38) and (4.39) hold:
(4.38)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{N|{x_{N}}-y{|^{2}}}dy{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.39)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}=0.\end{aligned}\]
In particular, (4.18) holds.
Proof.
From (4.29), we see that
(4.40)
\[ \underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{N|{x_{N}}-y{|^{2}}}dy{\bigg|\bigg|_{R}}\le \frac{2{c_{\text{10}}}}{r}.\]
Define a positive constant ${c_{11}}$ as
\[ {c_{\text{11}}}=\underset{N\to \infty }{\limsup }\underset{y\in {T^{N}}}{\sup }||\hspace{0.1667em}|{x_{N}}-y{|^{-1}}|{|_{R}}.\]
Clearly, ${c_{\text{11}}}<\infty $. From (3.26), we have
(4.41)
\[ \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {T^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{N|{x_{N}}-y{|^{2}}}dy{\bigg|\bigg|_{R}}\le \underset{N\to \infty }{\lim }\frac{4{c_{\text{4}}}{c_{\text{11}}^{2}}{N^{1/3}}{N^{\gamma }}}{N}=0.\]
Here, we used the fact that $|{T^{N}}|=4{N^{\gamma }}$ and $\gamma <-2/5$. Therefore, (4.38) follows from (4.40) and (4.41).
Next we prove (4.39). First, we estimate the integration on ${U^{N}}\times {U^{N}}$. We can see
(4.42)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\ge {N^{-1}}\}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}=0,\end{aligned}\]
(4.43)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\le {N^{-1}}\}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}=0.\end{aligned}\]
Actually, from (4.31) and the Schwartz inequality, we have
\[\begin{aligned}{}& \bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\ge {N^{-1}}\}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\ge {N^{-1}}\}}}\\ {} & \hspace{1em}\times \frac{{c_{\text{10}}^{2}}}{{N^{2}}|y-z{|^{2}}|2-{y^{2}}{|^{1/2}}|2-{z^{2}}{|^{1/2}}|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\ge {N^{-1}}\}}}\frac{{c_{\text{10}}^{2}}}{{N^{2}}|y-z{|^{2}}|2-{y^{2}}||{x_{N}}-y{|^{2}}}dydz{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U^{N}}}}\frac{{c_{\text{10}}^{2}}}{{N^{2}}|2-{y^{2}}||{x_{N}}-y{|^{2}}}\Bigg\{{\int _{\{|y-z|\ge {N^{-1}}\}}}\frac{1}{|y-z{|^{2}}}dz\Bigg\}dy{\bigg|\bigg|_{R}}\\ {} & =\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U^{N}}}}\frac{2{c_{\text{10}}^{2}}}{N|2-{y^{2}}||{x_{N}}-y{|^{2}}}dy{\bigg|\bigg|_{R}}.\end{aligned}\]
Hence, computation similar to (4.35) derives
\[ \underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\ge {N^{-1}}\}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\le \frac{4{c_{\text{10}}^{2}}}{r},\]
which implies (4.42). We next show (4.43). From (4.27) and (4.29), it holds that
\[\begin{aligned}{}& \bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\le {N^{-1}}\}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\le {N^{-1}}\}}}\frac{{c_{\text{10}}^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}\cap \{|y-z|\le {N^{-1}}\}}}\frac{{c_{\text{10}}^{2}}}{2}\Big\{\frac{1}{|{x_{N}}-y{|^{2}}}+\frac{1}{|{x_{N}}-z{|^{2}}}\Big\}dydz{\bigg|\bigg|_{R}}\\ {} & \le \frac{2{c_{\text{10}}^{2}}}{N}\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{1}{|{x_{N}}-y{|^{2}}}dy{\bigg|\bigg|_{R}}=\frac{4{c_{\text{10}}^{2}}}{r},\end{aligned}\]
which implies (4.43). Then, (4.42) and (4.43) gives
(4.44)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {U^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}=0.\end{aligned}\]
We next consider the integration on ${T^{N}}\times {T^{N}}$. From (3.26) and (4.27), we obtain
(4.45)
\[\begin{aligned}{}& \underset{N\to \infty }{\lim }\bigg|\bigg|{\int _{{({T_{x,r}^{N}}\cap {T^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & \hspace{2em}\hspace{2em}\le \underset{N\to \infty }{\lim }{c_{\text{4}}^{2}}{c_{\text{11}}^{2}}{N^{2/3}}{(4{N^{\gamma }})^{2}}=0.\end{aligned}\]
Here, we used $|{T^{N}}|=4{N^{\gamma }}$ and $\gamma <-2/5$.
Lastly, we consider the case ${U^{N}}\times {T^{N}}$. A similar argument deduces
(4.46)
\[\begin{aligned}{}& \bigg|\bigg|{\int _{({T_{x,r}^{N}}\cap {U^{N}})\times ({T_{x,r}^{N}}\cap {T^{N}})}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & \le \bigg|\bigg|{\int _{({T_{x,r}^{N}}\cap {U^{N}})\times ({T_{x,r}^{N}}\cap {T^{N}})}}\frac{{M_{\alpha }^{N}}(y,y){M_{\alpha }^{N}}(z,z)}{|{x_{N}}-y||{x_{N}}-z|}dydz{\bigg|\bigg|_{R}}\\ {} & =\bigg|\bigg|{\int _{{T_{x,r}^{N}}\cap {U^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{|{x_{N}}-y|}dy{\int _{{T_{x,r}^{N}}\cap {T^{N}}}}\frac{{M_{\alpha }^{N}}(z,z)}{|{x_{N}}-z|}dz{\bigg|\bigg|_{R}}\\ {} & =\mathcal{O}(\log N)\mathcal{O}({N^{1/3+\gamma }})\to 0\hspace{2.5pt}\text{as}\hspace{2.5pt}N\to \infty .\end{aligned}\]
Collecting (4.44), (4.45), and (4.46), we finally obtain (4.39). Clearly, (4.18) follows from (4.38) and (4.39).  □
Lemma 4.8.
Equation (4.19) holds.
Proof.
We estimate three terms in (4.19).
The first term in (4.19) vanishes. Indeed, from (4.27) and (4.38), we obtain
(4.47)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{1}{N|{x_{N}}-y{|^{2}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}\\ {} & \hspace{2em}\hspace{2em}\le \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}(y,y)}{N|{x_{N}}-y{|^{2}}}dy{\bigg|\bigg|_{R}}=0.\end{aligned}\]
For the second term in (4.19), from the Schwartz inequality, we have
\[\begin{aligned}{}& \bigg|\bigg|{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}(y,z){M_{\alpha }^{N}}({x_{N}},y){M_{\alpha }^{N}}({x_{N}},z)}{|{x_{N}}-y||{x_{N}}-z|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dydz{\bigg|\bigg|_{R}}\\ {} & =\bigg|\bigg|{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}{(y,z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|}dydz\bigg|{\bigg|_{R}^{\frac{1}{2}}}\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy{\bigg|\bigg|_{R}}.\end{aligned}\]
Thereby, (4.32) and (4.39) implies
(4.48)
\[ \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}(y,z){M_{\alpha }^{N}}({x_{N}},y){M_{\alpha }^{N}}({x_{N}},z)}{|{x_{N}}-y||{x_{N}}-z|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dydz{\bigg|\bigg|_{R}}=0.\]
Lastly, the third term in (4.19) converges to zero. Indeed, from (4.32), we see
(4.49)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{({T_{x,r}^{N}})^{2}}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}{M_{\alpha }^{N}}{({x_{N}},z)^{2}}}{|{x_{N}}-y||{x_{N}}-z|{M_{\alpha }^{N}}{({x_{N}},{x_{N}})^{2}}}dydz{\bigg|\bigg|_{R}}\\ {} & =\underset{r\to \infty }{\lim }\underset{N\to \infty }{\limsup }\bigg|\bigg|{\int _{{T_{x,r}^{N}}}}\frac{{M_{\alpha }^{N}}{({x_{N}},y)^{2}}}{|{x_{N}}-y|{M_{\alpha }^{N}}({x_{N}},{x_{N}})}dy\bigg|{\bigg|_{R}^{2}}=0.\end{aligned}\]
Therefore, we conclude (4.19) from (4.47), (4.48), and (4.49).  □
Proof of Theorem 2.11.
It is easy to see (B1). Indeed, (4.1) is verified by (1.6), and (4.2) follows from (1.9) and the Hadamard inequality. Additionally, (B2) holds for $\hat{p}=2$ under the setting of (4.6). (B2) (i) is clear, and direct computation with (4.11) yields (4.3) in (ii). Furthermore, we have (B2) (iii), since (4.16)–(4.19) in Lemma 4.3 are followed from Lemma 4.5, Lemma 4.6, Lemma 4.7, and Lemma 4.8, respectively. Then, we have checked the all conditions in Lemma 4.1, which completes the proof of Theorem 2.11.  □

5 Quasi-Gibbs property

5.1 Sufficient conditions for the quasi-Gibbs property

To show Theorem 2.6, we use sufficient conditions for the quasi-Gibbs property of a random point field μ with logarithmic interaction which is derived in [27]. We set the following conditions (QG1)–(QG2):
  • (QG1) There exists a sequence of random point fields ${\{{\mu ^{N}}\}_{N\in \mathbb{N}}}$ satisfying (B1) and the following (i)–(iii):
    • (i) ${\mu ^{N}}(\mathsf{s}(S)=N)=1$ for each $N\in \mathbb{N}$.
    • (ii) ${\mu ^{N}}$ is a $({\Phi ^{N}},-\beta \log |x-y|)$-canonical Gibbs measure for each $N\in \mathbb{N}$.
    • (iii) For each R, self-potential ${\Phi ^{N}}$ satisfies the following:
      \[ \underset{N\to \infty }{\lim }{\Phi ^{N}}(x)=\Phi (x)\hspace{2.5pt}\text{for a.e.}\hspace{2.5pt}x,\hspace{2em}\underset{n\in \mathbb{N}}{\inf }\underset{x\in {S_{R}}}{\inf }{\Phi ^{N}}(x)>-\infty .\]
Let $\mathsf{x}={\textstyle\sum _{i}}{\delta _{{x_{i}}}}$. For $l,r\in \mathbb{N}$, we define ${\mathsf{v}_{l,r}}:\mathsf{S}\to \mathbb{R}$ as
\[ {\mathsf{v}_{l,r}}(\mathsf{x})=\beta \sum \limits_{{x_{i}}\in {S_{r}^{c}}}\frac{1}{{x_{i}^{l}}}.\]
  • (QG2) There exists ${l_{0}}\in \mathbb{N}$ such that
    \[ \underset{N\in \mathbb{N}}{\sup }\Big\{{\int _{1\le |x|<\infty }}\frac{1}{|x{|^{{l_{0}}}}}{\rho ^{N,1}}(x)dx\Big\}<\infty \]
    and that, for each $1\le l<{l_{0}}$,
    \[ \underset{r\to \infty }{\lim }\underset{N\in \mathbb{N}}{\sup }||{\mathsf{v}_{l,r}}|{|_{{L^{1}}(\mathsf{S},{\mu ^{N}})}}=0.\]
Lemma 5.1.
[27, Theorem 2.2] Assume (QG1) and (QG2). Then, μ is a $(\Phi ,-\beta \log |x-y|)$-quasi-Gibbs measure.

5.2 Proof of Theorem 2.6

We first check (QG2) to use Lemma 5.1.
Lemma 5.2.
It holds that
(5.1)
\[\begin{aligned}{}& \underset{N\in \mathbb{N}}{\sup }\Big\{{\int _{1\le |x|<\infty }}\frac{1}{{x^{2}}}{\rho _{\alpha }^{N,1}}(x)dx\Big\}<\infty ,\end{aligned}\]
(5.2)
\[\begin{aligned}{}& \underset{r\to \infty }{\lim }\underset{N\in \mathbb{N}}{\sup }||{\mathsf{v}_{1,r}}|{|_{{L^{1}}(\mathsf{S},{\mu _{\alpha }^{N}})}}=0.\end{aligned}\]
In particular, we have (QG2) for ${l_{0}}=2$.
Proof.
Recall ${\rho _{\alpha }^{N,1}}(x)={K_{\alpha }^{N}}(x,x)={M_{\alpha }^{N}}(x/N,x/N)$. Then, from (3.25), there exists a positive constant ${c_{12}}$ such that
\[\begin{aligned}{}& \underset{N\in \mathbb{N}}{\sup }\underset{1\le |x|\le N}{\sup }{\rho _{\alpha }^{N,1}}(x)\le {c_{\text{12}}}.\end{aligned}\]
Therefore, we have
(5.3)
\[ {\int _{1\le |x|<\infty }}\frac{1}{|x{|^{2}}}{\rho _{\alpha }^{N,1}}(x)dx\le 2{c_{\text{12}}}{\int _{1}^{N}}\frac{1}{{x^{2}}}dx+2{\int _{N}^{\infty }}\frac{{\rho _{\alpha }^{N,1}}(x)}{{x^{2}}}dx.\]
Using ${\textstyle\int _{S}}{\rho _{\alpha }^{N,1}}(x)dx=N$, we see that
\[ {\int _{N}^{\infty }}\frac{{\rho _{\alpha }^{N,1}}(x)}{{x^{2}}}dx\le \frac{1}{{N^{2}}}{\int _{N}^{\infty }}{\rho _{\alpha }^{N,1}}(x)dx\le \frac{1}{N}.\]
Thereby, combining this with (5.3), we obtain (5.1). Moreover, we see that
\[ ||{\mathsf{v}_{1,r}}|{|_{{L^{1}}(\mathsf{S},{\mu _{\alpha }^{N}})}}={\int _{|x|>r}}\frac{2{\rho _{\alpha }^{N,1}}(x)}{x}dx={\int _{|x|>\frac{r}{N}}}\frac{2{M_{\alpha }^{N}}(x,x)}{x}dx.\]
Then, by arguments similar to that yield Lemma 4.5, we prove (5.2). Therefore, we obtain (QG2) for ${l_{0}}=2$.  □
Proof of Theorem 2.6.
Taking ${\mu ^{N}}={\mu _{\alpha }^{N}}$ and $\mu ={\mu _{\alpha }}$, we check assumption (QG1) and (QG2). It is easily to see (QG1) with ${\Phi ^{N}}(x)={x^{2}}{N^{-1}}-2\alpha \log |x|$. Furthermore, (QG2) holds for ${l_{0}}=2$ from Lemma 5.2. Hence, we conclude Theorem 2.6 from Lemma 5.1.  □

6 Strong uniqueness

6.1 General framework for strong uniqueness

Following [29], we introduce a framework of strong uniqueness of ISDEs in the case when $S=\mathbb{R}$ and the diffusion coefficient $\sigma =1$. In order not to make this paper too long, several terminologies and symbols are not defined here. See [29] for precise definitions.
Recall that ${\mathsf{S}_{\mathrm{s},\mathrm{i}}}$ is defined in (2.1). Let $\mathsf{H}$ and ${\mathsf{S}_{\mathrm{sde}}}$ be Borel subsets of $\mathsf{S}$ such that
\[ \mathsf{H}\subset {\mathsf{S}_{\mathrm{sde}}}\subset {\mathsf{S}_{\mathrm{s},\mathrm{i}}}.\]
Define ${\mathbf{S}_{\mathrm{sde}}}\subset {S^{\mathbb{N}}}$ and ${\mathsf{S}_{\mathrm{sde}}^{[1]}}\subset S\times \mathsf{S}$ as
\[ {\mathbf{S}_{\mathrm{sde}}}={\mathfrak{u}^{-1}}({\mathsf{S}_{\mathrm{sde}}}),\hspace{2em}{\mathsf{S}_{\mathrm{sde}}^{[1]}}={\mathfrak{u}_{[1]}^{-1}}({\mathsf{S}_{\mathrm{sde}}}),\]
where ${\mathfrak{u}_{[1]}}:S\times \mathsf{S}\to \mathsf{S}$ is given by ${\mathfrak{u}_{[1]}}(x,\mathsf{s})={\delta _{x}}+\mathsf{s}$.
Let $(\mathbf{X},\mathbf{B})$ be an ${S^{\mathbb{N}}}\times {\mathbb{R}^{\mathbb{N}}}$-valued continuous process defined on a filtered space $(\Omega ,\mathfrak{F},P,\{{\mathfrak{F}_{t}}\})$. We assume that $(\Omega ,\mathfrak{F},P)$ is a standard probability space. Then, the regular conditional probability ${P_{\mathbf{s}}}=P(\cdot |{\mathbf{X}_{0}}=\mathbf{s})$ exists for s, $P\circ {\mathbf{X}_{0}^{-1}}$-almost everywhere.
Let $b:{\mathsf{S}_{\mathrm{sde}}^{[1]}}\to \mathbb{R}$ be a Borel measurable function. Then, consider an ISDE of $\mathbf{X}={({X^{i}})_{i\in \mathbb{N}}}$ starting from $\mathfrak{l}(\mathsf{H})$ with state space ${\mathbf{S}_{\mathrm{sde}}}$:
(6.1)
\[\begin{aligned}{}& d{X_{t}^{i}}=d{B_{t}^{i}}+b({X_{t}^{i}},{\mathsf{X}_{t}^{i\lozenge }})dt,\end{aligned}\]
(6.2)
\[\begin{aligned}{}& \mathbf{X}\in W({\mathbf{S}_{\mathrm{sde}}}),\\ {} & {\mathbf{X}_{0}}\in \mathfrak{l}(\mathsf{H}).\end{aligned}\]
For $\mathbf{X}={({X^{i}})_{i\in \mathbb{N}}}$, we set ${\mathsf{X}_{t}^{m\ast }}={\textstyle\sum _{i=m+1}^{\infty }}{\delta _{{X_{t}^{i}}}}$. For $(u,\mathbf{v})\in {S^{m}}$ and $\mathsf{v}={\textstyle\sum _{i=1}^{m-1}}{\delta _{{v_{i}}}}$, where $\mathbf{v}=({v_{1}},\dots ,{v_{m-1}})$, we define ${b_{\mathsf{X}}^{m}}:[0,\infty )\times {S^{m}}\to \mathbb{R}$ as
\[ {b_{\mathsf{X}}^{m}}(t,(u,\mathbf{v}))=b(u,\mathsf{v}+{\mathsf{X}_{t}^{m\ast }}).\]
Let
\[ {\mathbf{S}_{\mathrm{sde}}^{m}}(t,\mathsf{w})=\{{\mathbf{s}_{m}}=({s_{1}},\dots ,{s_{m}})\in {S^{m}}\hspace{0.1667em};\hspace{0.1667em}\mathfrak{u}({\mathbf{s}_{m}})+{\mathsf{w}_{t}^{m\ast }}\in {\mathsf{S}_{\mathrm{sde}}}\},\]
where ${\mathsf{w}_{t}^{m\ast }}={\textstyle\sum _{i=m+1}^{\infty }}{\delta _{{w_{t}^{i}}}}$ for ${\mathsf{w}_{t}}={\textstyle\sum _{i=1}^{\infty }}{\delta _{{w_{t}^{i}}}}$. Let ${\mathbf{Y}^{m}}=({Y^{1}},\dots ,{Y^{m}})$ be a solution to the following SDE with random environment $\mathsf{X}$ defined on $(\Omega ,\mathfrak{F},{P_{\mathbf{s}}},\{{\mathfrak{F}_{t}}\})$:
(6.3)
\[\begin{aligned}{}& d{Y_{t}^{m,i}}=d{B_{t}^{i}}+{b_{\mathsf{X}}^{m}}(t,({Y_{t}^{m,i}},{\mathbf{Y}_{t}^{m,i\lozenge }}))dt,\end{aligned}\]
(6.4)
\[\begin{aligned}{}& {\mathbf{Y}_{t}^{m}}\in {\mathbf{S}_{\mathrm{sde}}^{m}}(t,\mathsf{X})\hspace{2.5pt}\text{for all}\hspace{2.5pt}t,\end{aligned}\]
(6.5)
\[\begin{aligned}{}& {\mathbf{Y}_{t}^{m}}={\mathbf{s}_{m}}.\end{aligned}\]
Here, we set ${\mathbf{Y}_{t}^{i\lozenge }}={({Y^{m,j}})_{j\ne i}^{m}}$ and ${\mathbf{s}_{m}}=({s_{1}},\dots ,{s_{m}})$ for $\mathbf{s}=({s_{i}})\in {S^{\mathbb{N}}}$.
We formulate strong solution to (6.3)–(6.5) and its uniqueness. We set ${W_{0}}({\mathbb{R}^{m}})=\{\mathbf{w}\in W({\mathbb{R}^{m}})\hspace{0.1667em};\hspace{0.1667em}{\mathbf{w}_{0}}=\mathbf{0}\}$. Let ${\mathfrak{C}^{m}}$, ${\mathfrak{C}_{t}^{m}}$ be σ-fields on ${W_{0}}({\mathbb{R}^{m}})\times W({S^{\mathbb{N}}})$ defined as in [29, 1154 p.]. Let ${\mathfrak{B}_{t}^{m}}=\sigma [{\mathbf{w}_{s}}\hspace{0.1667em};\hspace{0.1667em}0\le s\le t,\mathbf{w}\in W({\mathbb{R}^{m}})]$.
Definition 6.1 ([29, Definition 3.9]).
We say that ${\mathbf{Y}^{m}}$ a strong solution to (6.3)–(6.5) for $(\mathbf{X},\mathbf{B})$ under ${P_{\mathbf{s}}}$ if $({\mathbf{Y}^{m}},{\mathbf{B}^{m}},{\mathbf{X}^{m\ast }})$ satisfies (6.3)–(6.5) and there exists a function
\[ {F_{\mathbf{s}}^{m}}:{W_{0}}({\mathbb{R}^{m}})\times W({S^{\mathbb{N}}})\to W({S^{m}})\]
such that ${\mathfrak{C}^{m}}$-measurable and ${\mathfrak{C}_{t}^{m}}/{\mathfrak{B}_{t}^{m}}$-measurable for any t, and it holds that
\[ {\mathbf{Y}^{m}}={F_{\mathbf{s}}^{m}}({\mathbf{B}^{m}},{\mathbf{X}^{m\ast }})\hspace{2em}\text{for}\hspace{2.5pt}{P_{\mathbf{s}}}\text{-a.s.}\]
Here, ${\mathbf{B}^{m}}=({B^{1}},\dots ,{B^{m}})$ is the first m-components of the Brownian motion B.
Definition 6.2 ([29, Definition 3.10]).
SDE (6.3)–(6.5) is said to have a unique strong solution for $(\mathbf{X},\mathbf{B})$ under ${P_{\mathbf{s}}}$ if there exists a function ${F_{\mathbf{s}}^{m}}$ satisfying the same condition as in Definition 6.1 and, for any weak solution $({\hat{\mathbf{Y}}^{m}},{\mathbf{B}^{m}},{\mathbf{X}^{m\ast }})$ to (6.3)–(6.5) under ${P_{\mathbf{s}}}$, we have
\[ {\hat{\mathbf{Y}}^{m}}={F_{\mathbf{s}}^{m}}({\mathbf{B}^{m}},{\mathbf{X}^{m\ast }}).\]
Then, we introduce the IFC condition for $(\mathbf{X},\mathbf{B})$ defined on $(\Omega ,\mathfrak{F},P,\{{\mathfrak{F}_{t}}\})$.
  • (IFC) For each $m\in \mathbb{N}$, SDE (6.3)–(6.5) has a unique strong solution ${F_{\mathbf{s}}^{m}}({\mathbf{B}^{m}},{\mathbf{X}^{m\ast }})$ for $(\mathbf{X},\mathbf{B})$ under ${P_{\mathbf{s}}}$ for $P\circ {\mathbf{X}_{0}^{-1}}$-a.s. s.
We next introduce several conditions for the strong uniqueness of solutions to ISDEs. See [29, Section 3] for the definitions of solutions to ISDEs such as strong solutions and unique strong solutions under constraints. For a process $\mathbf{X}={({X^{i}})_{i\in \mathbb{N}}}\in W({S^{\mathbb{N}}})$, let ${\mathsf{X}_{t}}={\textstyle\sum _{i}}{\delta _{{X_{t}^{i}}}}$. We make assumptions for μ and X under P.
  • (TT) μ is tail trivial.
  • (μ-AC) $P\circ {\mathsf{X}_{t}^{-1}}\prec \mu $ for all $0<t<\infty $.
  • (SIN) $P(\mathsf{X}\in {W_{\mathrm{NE}}}({\mathsf{S}_{\mathrm{s},\mathrm{i}}}))=1$.
  • (NBJ) $P({\mathsf{m}_{r,T}}(\mathbf{X})<\infty )=1$ for each $r,T\in \mathbb{N}$, where, for $\mathbf{w}\in W({S^{\mathbb{N}}})$,
    \[ {\mathsf{m}_{r,T}}(\mathbf{w})=\inf \{m\in \mathbb{N}\hspace{0.1667em};\hspace{0.1667em}\underset{t\in [0,T]}{\min }|{w_{t}^{n}}|>r\hspace{2.5pt}\text{for all}\hspace{2.5pt}n\in \mathbb{N}\hspace{2.5pt}\text{such that}\hspace{2.5pt}n>m\}.\]
Let ${F_{\mathbf{s}}}:{W_{0}}({\mathbb{R}^{\mathbb{N}}})\to W({\mathbb{R}^{\mathbb{N}}})$ be a strong solution to (6.1)–(6.2) starting at s for $P\circ {\mathbf{X}_{0}^{-1}}$-a.s. s.
  • (MF) $P({F_{\mathbf{s}}}(\mathbf{B})\in A)$ is ${\overline{\mathfrak{B}({S^{\mathbb{N}}})}^{P\circ {\mathbf{X}_{0}^{-1}}}}$-measurable in s for any $A\in \mathfrak{B}(W({S^{\mathbb{N}}}))$.
Lemma 6.3 ([29, Theorem 3.1]).
Assume (TT). Assume that (6.1)–(6.2) has a weak solution under P satisfying (IFC), (μ-AC), (SIN), and (NBJ). Then, (6.1)–(6.2) has a family of unique strong solutions $\{{F_{\mathbf{s}}}\}$ starting at s for $P\circ {\mathbf{X}_{0}^{-1}}$-a.s. s under the constraints of (MF), (IFC), (μ-AC), (SIN), and (NBJ).

6.2 Proof of Theorem 2.2

Let $(\mathsf{X},\{{\mathsf{P}_{\mathsf{s}}}\})$ be the diffusion associated with $({\mathcal{E}^{{\mu _{\alpha }}}},{\mathcal{D}^{{\mu _{\alpha }}}})$, which is obtained in Corollary 2.7. We take X and P in Section 6.1 as $\mathbf{X}={\mathfrak{l}_{\mathrm{path}}}(\mathsf{X})$ and $P={\mathsf{P}_{{\mu _{\alpha }}}}:=\textstyle\int {\mathsf{P}_{\mathsf{s}}}d{\mu _{\alpha }}$.
Proof of Theorem 2.2.
We shall check assumptions of Lemma 6.3. Since ${\mu _{\alpha }}$ is determinantal, (TT) for ${\mu _{\alpha }}$ is known [3, 18, 28]. From the fact that $\mathsf{X}$ is reversible with respect to ${\mu _{\alpha }}$, (${\mu _{\alpha }}$-AC) holds. Furthermore, (SIN) and (NBJ) hold from the same argument as in [29, Section 10.].
We can check (IFC) for $(\mathbf{X},\mathbf{B})$ using a result of [14]. We take ${a_{q}}(r)=C(q){r^{1+\varepsilon }}$, where ε is a sufficiently small positive constant and ${\{C(q)\}_{q\in \mathbb{N}}}$ is a increasing sequence. It is easy to see that there exists a sequence $\{C(q)\}$ such that (5.17) and (5.41) in [14] holds. Using [14, Theorem 6.1], $\{\mathbf{B}\mathbf{1}\}$ in [14] holds for $(\mathbf{X},\mathbf{B})$. Remark that the assumption $\textstyle\int \tilde{\chi }d{\mu _{\alpha }}<\infty $ in [14, Theorem 6.1] is proved by [14, Lemma 5.4]. Moreover, we see that $\{\mathbf{C}\mathbf{1}\}$ and $\{\mathbf{C}\mathbf{2}\}$ in [14] holds from Theorem 2.11. Then, [14, Theorem 3.3] implies that $(\mathbf{X},\mathbf{B})$ satisfies $\{\mathbf{B}\mathbf{2}\}$ and $\{\mathbf{B}\mathbf{3}\}$ in [14]. Therefore, from [14, Theorem 3.2], $(\mathbf{X},\mathbf{B})$ satisfies (IFC).
Finally, we have checked all assumptions in Lemma 6.3, which completes the proof of Theorem 2.2.  □

7 Dynamical convergence

7.1 General theory of finite particle approximation to ISDEs

This subsection is devoted to prepare a general result for finite particle approximations to ISDEs. Recall that ${\mu ^{N}}$ and μ satisfy (B1). Additionally, we suppose the following:
  • (B3) For each $m\in \mathbb{N}$,
    \[ \underset{N\to \infty }{\lim }{\mu ^{N}}\circ {({\mathfrak{l}_{m}^{N}})^{-1}}=\mu \circ {({\mathfrak{l}_{m}})^{-1}}\hspace{1em}\hspace{2.5pt}\text{weakly in}\hspace{2.5pt}{S^{m}}.\]
We shall later take ${\mu ^{N}}\circ {({\mathfrak{l}^{N}})^{-1}}$ as an initial distribution of labelled processes for finite particle systems. Hence, (B3) means weak convergence of the initial distribution of the labelled dynamics.
For a labelled process ${\mathbf{X}^{N}}={({X^{N,i}})_{i=1}^{N}}\in W({S^{N}})$, we set
\[\begin{aligned}{}& {\mathsf{X}_{t}^{N,\lozenge i}}={\sum \limits_{j\ne i}^{N}}{\delta _{{X_{t}^{N,j}}}}.\end{aligned}\]
Recall that ${\mathsf{d}^{{\mu ^{N}}}}$ denotes the logarithmic derivative of ${\mu ^{N}}$. We introduce a finite-dimensional SDE of ${\mathbf{X}^{N}}={({X^{N,i}})_{i=1}^{N}}$: for $1\le i\le N$,
(7.1)
\[\begin{aligned}{}d{X_{t}^{N,i}}& =d{B_{t}^{i}}+\frac{1}{2}{\mathsf{d}^{{\mu ^{N}}}}({X_{t}^{N,i}},{\mathsf{X}_{t}^{N,\lozenge i}})dt.\end{aligned}\]
Then, we assume the following:
  • (B4) For each N, there exists the logarithmic derivative ${\mathsf{d}^{{\mu ^{N}}}}$ of ${\mu ^{N}}$ and SDE (7.1) with initial condition ${\mathbf{X}_{0}^{N}}={\mathbf{s}_{N}}$ has a unique solution for ${\mu ^{N}}\circ {({\mathfrak{l}^{N}})^{-1}}$-a.s. ${\mathbf{s}_{N}}$. Furthermore, this solution does not explode.
  • (B5) There exists $T>0$ such that, for each $R>0$,
    \[ \underset{r\to \infty }{\liminf }\underset{N\in \mathbb{N}}{\sup }\Big\{{\int _{|x|\le r+R}}{\rho ^{N,1}}(x)dx\Big\}\mathrm{Erf}\Big(\frac{r}{\sqrt{(r+R)T}}\Big)=0.\]
    Moreover, we write ${s_{i}}={\mathfrak{l}^{N,i}}(\mathsf{s})$ and, for each $R,T\in \mathbb{N}$,
    \[ \underset{L\to \infty }{\lim }\underset{N\to \infty }{\limsup }\sum \limits_{i>L}{\int _{\mathsf{S}}}\mathrm{Erf}\Big(\frac{|{s_{i}}|-R}{T}\Big){\mu ^{N}}(d\mathsf{s})=0.\]
We say ${\{{\mathbf{X}^{N}}\}_{N\in \mathbb{N}}}$ is tight in $W({S^{\mathbb{N}}})$ if each subsequence ${\mathbf{X}^{{N^{\prime }}}}$ of ${\mathbf{X}^{N}}$ contains a subsequence ${\mathbf{X}^{{N^{\prime\prime }}}}$ such that, for each $m\in \mathbb{N}$, ${\mathbf{X}^{{N^{\prime\prime }},m}}={({X^{{N^{\prime\prime }},i}})_{i=1}^{m}}$ is convergent weakly in $W({S^{m}})$. The following theorem is a special case of [12, Theorem 2.2].
Lemma 7.1.
Assume (A3), (A4), and (B1)–(B5). Assume that ${\mathbf{X}_{0}^{N}}={\mu ^{N}}\circ {({\mathfrak{l}^{N}})^{-1}}$ in distribution. Then, the following (1) and (2) hold:
  • (i) The family of processes ${\{{\mathbf{X}^{N}}\}_{N\in \mathbb{N}}}$ is tight in $W({S^{\mathbb{N}}})$, and each limit point X of ${\{{\mathbf{X}^{N}}\}_{N\in \mathbb{N}}}$ is a solution to
    \[\begin{aligned}{}& d{X_{t}^{i}}=d{B_{t}^{i}}+\frac{1}{2}{\mathsf{d}^{\mu }}({X_{t}^{i}},{\mathsf{X}_{t}^{i\lozenge }})dt\end{aligned}\]
    with initial distribution $\mu \circ {\mathfrak{l}^{-1}}$.
  • (ii) If each limit point X satisfies (IFC), (SIN), and (NBJ), then, for each $m\in \mathbb{N}$,
    (7.2)
    \[ \underset{N\to \infty }{\lim }({X^{N,1}},{X^{N,2}},\dots ,{X^{N,m}})=({X^{1}},{X^{2}},\dots ,{X^{m}})\]
    weakly in $W({S^{m}})$.
Proof.
Claim (i) follows from [12, Theorem 2.2]. Since X is a limit of $\{{\mathbf{X}^{N}}\}$ and ${\textstyle\sum _{1\le i\le N}}{\delta _{{X^{N,i}}}}$ starts from reversible measure ${\mu ^{N}}$, each limit point X satisfies (μ-AC). Then, from [29, Corollary 3.2], the distribution of each limit point X is unique. Therefore, combining the uniqueness with the tightness proved in (i), we obtain (7.2).  □

7.2 Proof of Theorem 2.4

We use Lemma 7.1 (ii) for $\mu ={\mu _{\alpha }}$ and ${\mu ^{N}}={\mu _{\alpha }^{N}}$ to prove Theorem 2.4.
Proof of Theorem 2.4.
We have already shown (B1) and (B2). Assumption (B3) corresponds (2.4). It is easy to see that (B4) holds. With the aid of [12, Lemma 4.6], we can check (B5).
Let X be an arbitrary limit point of ${\mathbf{X}^{N}}$. From Lemma 7.1 (1), X is a weak solution to (2.2) with initial distribution ${\mu _{\alpha }}\circ {\mathfrak{l}^{-1}}$. Since $\mathfrak{u}(\mathbf{X})$ is reversible with respect to ${\mu _{\alpha }}$, X has the Lyons–Zheng type decomposition (see [14, Section 9]). Then, the noncollision property holds from [14, Theorem 6.1]. It is easy to see the nonexplosion of tagged particles. Therefore, we have (SIN) for X. Furthermore, the Lyons–Zheng type decomposition makes it possible to use the same argument as in [29, Lemma 10.3], which proves (NBJ). From the same argument as the proof of Theorem 2.2, X satisfies (IFC).
Thus, Lemma 7.1 (ii) completes the proof of Theorem 2.4.  □

References

[1] 
Akemann, G., Damgaard, P.H., Magnea, U., Nishigaki, S.: Universality of random matrices in the microscopic limit and the Dirac operator spectrum. Nucl. Phys. B 487(3), 721–738 (1997). MR1432823. https://doi.org/10.1016/S0550-3213(96)00713-4
[2] 
Bufetov, A.I., Dymov, A.V., Osada, H.: The logarithmic derivative for point processes with equivalent Palm measures. J. Math. Soc. Jpn. 71(2), 451–469 (2019). MR3943446. https://doi.org/10.2969/jmsj/78397839
[3] 
Bufetov, A.I., Qiu, Y., Shamov, A.: Kernels of conditional determinantal measures and the Lyons–Peres completeness conjecture. J. Eur. Math. Soc. 23(5), 1477–1519 (2021). MR4244512. https://doi.org/10.4171/JEMS/1038
[4] 
Olver, F.W.J., Lozier, D.W., Boisvert, R.F., Clark, C.W. (eds.): NIST Handbook of Mathematical Functions Cambridge University Press, Cambridge (2010). MR2723248
[5] 
Erdélyi, A.: Asymptotic forms for Laguerre polynomials. J. Indian Math. Soc. (N.S.) 24(1960), 235–250 (1961). MR0123751
[6] 
Forrester, P.J.: Log-Gases and Random Matrices. London Mathematical Society Monographs. Princeton University Press, Princeton (2010). MR2641363. https://doi.org/10.1515/9781400835416
[7] 
Fritz, J.: Gradient dynamics of infinite point systems. Ann. Probab. 15, 478–514 (1987). MR0885128
[8] 
Fukushima, M., Oshima, Y., Takeda, M.: Dirichlet Forms and Symmetric Markov Processes, 2nd edn. Walter de Gruyter, Berlin (2011). MR2778606
[9] 
Honda, R., Osada, H.: Infinite-dimensional stochastic differential equations related to the Bessel random point fields. Stoch. Process. Appl. 125(10), 3801–3822 (2015). MR3373304. https://doi.org/10.1016/j.spa.2015.05.005
[10] 
Kanzieper, E., Freilikher, V.: Random matrix models with log-singular level confinement: method of fictitious fermions. Philos. Mag. B 77(5), 1161–1172 (1998). https://doi.org/10.1080/13642819808205006
[11] 
Katori, M., Tanemura, H.: Markov property of determinantal processes with extended sine, Airy, and Bessel kernels. Markov Process. Relat. Fields 17(4), 541–580 (2011). MR2918121
[12] 
Kawamoto, Y., Osada, H.: Finite-particle approximations for interacting Brownian particles with logarithmic potentials. J. Math. Soc. Jpn. 70(3), 921–952 (2018). MR3830792. https://doi.org/10.2969/jmsj/75717571
[13] 
Kawamoto, Y., Osada, H., Tanemura, H.: Uniqueness of Dirichlet forms related to infinite systems of interacting Brownian motions. Potential Anal. (2020). Published online. MR4341065. https://doi.org/10.1007/s11118-020-09872-2
[14] 
Kawamoto, Y., Osada, H., Tanemura, H.: Infinite-dimensional stochastic differential equations and tail σ-fields II: the IFC condition. J. Math. Soc. Japan. to appear. MR4126938. https://doi.org/10.1007/s00440-020-00981-y
[15] 
Kuijlaars, A.B.J., Vanlessen, M.: Universality for eigenvalue correlations at the origin of the spectrum. Commun. Math. Phys. 243(1), 163–191 (2003). MR2020225. https://doi.org/10.1007/s00220-003-0960-z
[16] 
Lang, R.: Unendlich-dimensionale Wienerprocesse mit Wechselwirkung I. Z. Wahrscheinlichkeitstheor. Verw. Geb. 38, 55–72 (1977). MR0431435. https://doi.org/10.1007/BF00534170
[17] 
Lang, R.: Unendlich-dimensionale Wienerprocesse mit Wechselwirkung II. Z. Wahrscheinlichkeitstheor. Verw. Geb. 39, 277–299 (1978). MR0431435. https://doi.org/10.1007/BF01877496
[18] 
Lyons, R.: A note on tail triviality for determinantal point processes. Electron. Commun. Probab. 23, 1–3 (2018). MR3866045. https://doi.org/10.1214/18-ECP175
[19] 
Ma, Z.-M., Röckner, M.: Introduction to the Theory of (Non-symmetric) Dirichlet Forms. Springer, Berlin (1992). MR1214375. https://doi.org/10.1007/978-3-642-77739-4
[20] 
Nagao, T., Slevin, K.: Nonuniversal correlations for random matrix ensembles. J. Math. Phys. 34(5), 2075–2085 (1993). MR1214509. https://doi.org/10.1063/1.530157
[21] 
Olver, F.W.J.: Asymptotics and Special Functions. A K Peters, Ltd., Wellesley, MA (1997). Reprint of the 1974 original, AKP Classics MR1429619
[22] 
Osada, H.: Dirichlet form approach to infinite-dimensional Wiener processes with singular interactions. Commun. Math. Phys. 176, 117–131 (1996). MR1372820. https://doi.org/10.1007/BF02099365
[23] 
Osada, H.: Non-collision and collision properties of Dyson’s model in infinite dimensions and other stochastic dynamics whose equilibrium states are determinantal random point fields. In: Funaki, T., Osada, H. (eds.) Stochastic Analysis on Large Scale Interacting Systems. Advanced Studies in Pure Mathematics, vol. 39, pp. 325–343 (2004). MR2073339. https://doi.org/10.2969/aspm/03910325
[24] 
Osada, H.: Tagged particle processes and their non-explosion criteria. J. Math. Soc. Jpn. 62(3), 867–894 (2010). MR2648065. https://doi.org/10.2969/jmsj/06230867
[25] 
Osada, H.: Infinite-dimensional stochastic differential equations related to random matrices. Probab. Theory Relat. Fields 153, 471–509 (2012). MR2948684. https://doi.org/10.1007/s00440-011-0352-9
[26] 
Osada, H.: Interacting Brownian motions in infinite dimensions with logarithmic interaction potentials. Ann. Probab. 41, 1–49 (2013). MR3059192. https://doi.org/10.1214/11-AOP736
[27] 
Osada, H.: Interacting Brownian motions in infinite dimensions with logarithmic interaction potentials II: airy random point field. Stoch. Process. Appl. 123, 813–838 (2013). MR3005007. https://doi.org/10.1016/j.spa.2012.11.002
[28] 
Osada, H., Osada, S.: Discrete approximations of determinantal point processes on continuous spaces: tree representations and tail triviality. J. Stat. Phys. 170, 421–435 (2018). MR3744393. https://doi.org/10.1007/s10955-017-1928-2
[29] 
Osada, H., Tanemura, H.: Infinite-dimensional stochastic differential equations and tail σ-fields. Probab. Theory Relat. Fields 177, 1137–1242 (2020). MR4126938. https://doi.org/10.1007/s00440-020-00981-y
[30] 
Shirai, T., Takahashi, Y.: Random point fields associated with certain Fredholm determinants I: fermion, Poisson and boson point process. J. Funct. Anal. 205, 414–463 (2003). MR2018415. https://doi.org/10.1016/S0022-1236(03)00171-X
[31] 
Tanemura, H.: A system of infinitely many mutually reflecting Brownian balls in ${\mathbb{R}^{d}}$. Probab. Theory Relat. Fields 104, 399–426 (1996). MR1376344. https://doi.org/10.1007/BF01213687
[32] 
Tsai, L.-C.: Infinite dimensional stochastic differential equations for Dyson’s model. Probab. Theory Relat. Fields 166, 801–850 (2016). MR3568040. https://doi.org/10.1007/s00440-015-0672-2
Exit Reading PDF XML


Table of contents
  • 1 Introduction
  • 2 Set up and main results
  • 3 Key tools for the proofs of main results
  • 4 The logarithmic derivative
  • 5 Quasi-Gibbs property
  • 6 Strong uniqueness
  • 7 Dynamical convergence
  • References

Export citation

Copy and paste formatted citation
Placeholder

Download citation in file


Share


RSS

MSTA

MSTA

  • Online ISSN: 2351-6054
  • Print ISSN: 2351-6046
  • Copyright © 2018 VTeX

About

  • About journal
  • Indexed in
  • Editors-in-Chief

For contributors

  • Submit
  • OA Policy
  • Become a Peer-reviewer

Contact us

  • ejournals-vmsta@vtex.lt
  • Mokslininkų 2A
  • LT-08412 Vilnius
  • Lithuania
Powered by PubliMill  •  Privacy policy