Abstract:
We consider a family of linear operators diagonalized by the Hankel transform. We express explicitly the Fredholm determinants of these operators, as restricted to $L_2[0, R]$, so that the rate of their convergence as $R\to\infty$ can be found. We use the link between these determinants and the distribution of additive functionals in a determinantal point process with Bessel kernel and estimate the distance in the Kolmogorov–Smirnov metric between the distribution of these functionals and the Gaussian distribution.
Bibliography: 27 titles.
where $J_\nu$ is the Bessel function of order $\nu$. Fix $\nu> -1$. The kernel (1.1) induces a bounded linear operator on $L_2(\mathbb{R}_+)$, which we will also denote by $B_f$. We let $B_1$ be equal to the identity operator. We refer to $B_f$ as the Bessel operator.
For any $h\in L_\infty(X)$ we also let $h$ denote the operator of pointwise multiplication by this function in $L_2(X)$. Let $\chi_A(x)$ be the characteristic function of the subset $A$.
In this work we consider the Fredholm determinant $\det(I + \chi_{[0, R]}B_{f}\chi_{[0, R]})$ for $R>0$ and a sufficiently smooth function $f$. This determinant gives an exact expression for the Laplace transform of additive functionals in the determinantal point process with Bessel kernel (see [26]). We describe the relation between the determinantal point process and the operator in § 5. Briefly, the operator $B_{\chi_{[0, 1]}}$ induces a probability measure $\mathbb{P}_{J_\nu}$ on countable subsets without accumulation points of $\mathbb{R}_+$. These subsets, called configurations, are endowed with a certain $\sigma$-algebra (see [25]). For any measurable function $b$ on $\mathbb{R}_+$ an additive functional is defined as a measurable function on configurations by $S_b(X) = \sum_{x\in X}b(x)$. The probability measure on configurations induces a random variable $S_b$. For $b\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$ the Laplace transform of $S_b$ is expressed by the formula
We conclude that the limit distribution of $S_b^R$ as $R\to\infty$ can be derived in terms of the asymptotic of the corresponding Fredholm determinant.
Recall an analogous problem for the sine process $\mathbb{P}_{\mathcal{S}}$. Unlike the determinantal point process with Bessel kernel, it is a measure on configurations on the real line $\mathbb{R}$. It is induced by the integral kernel $K_{\mathcal{S}}(x, y)=(\sin \pi(x-y))/(\pi(x-y))$. For $b\in L_1(\mathbb{R})\cap L_\infty(\mathbb{R})$ the formula for the Laplace transform is in this case
where $W_{e^{\lambda b}}$ is the Wiener–Hopf operator with symbol $e^{\lambda b}$, which we define in § 2. For these determinants we recall two results. The first is the Kac–Akhiezer formula (see [10], § 10.13), which states that under certain conditions on $b$ we have
where $Q^{\mathcal{S}}_R(\lambda b)\to 1$ as $R\to\infty$; see Theorem 3 for the values of $c_1^{\mathcal{S}}(b)$ and $c_2^{\mathcal{S}}(b)$. It follows that the distribution of additive functionals approaches here the Gaussian distribution if $c^{\mathcal{S}}_1(b)\!=\!0$. The second result is a precise formula for the remainder term $Q^{\mathcal{S}}_R(\lambda b)$; see Theorem 3 again. These results are continuous analogues of the strong Szegö limit theorem and the Borodin–Okounkov identity for Toeplitz matrices (see [10], § 10.4, and [11], respectively). Borodin and Okounkov derived an explicit expression for the determinants of Toeplitz matrices from Gessel’s theorem and Schur measures. Several different proofs under less restrictive assertions were given later [6], [9]. In a similar way a precise formula for the Wiener–Hopf operators was obtained by Basor and Chen [2]. Under weaker assumptions about the symbol this result was deduced by Bufetov [12] directly from the Borodin–Okounkov formula.
Having derived a formula for $Q^{\mathcal{S}}_R(\lambda b)$, it is possible to estimate the rate of convergence of additive functionals in the limit procedure under consideration. This was done by Bufetov in [13], Lemma 3.10. Theorem 2 gives a similar estimate for the Bessel kernel point process.
An analogue of the Kac–Akhiezer formula in case of the Bessel kernel was derived by Basor and Ehrhardt [3]. We state it in Theorem 6. One consequence of Theorem 6 is the rate at which the distribution of additive functionals approaches the Gaussian one if $c_1^{\mathcal{B}}(b)=0$. In this paper we derive a precise expression for $Q^{\mathcal{B}}_R(\lambda b)$ and estimate the rate of its convergence to 1.
This result has two special cases when $\nu=\pm 1/2$. In these cases the Bessel operator (1.1) is known to be the sum of the Wiener–Hopf and Hankel operators:
where $H_f$ is the Hankel integral operator with kernel $\chi_{[0, \infty)^2}(x, y)\widehat{f}(x+y)$. In these cases an expression for $Q^{\mathcal{B}}_R(b)$ was derived by Basor, Ehrhardt and Widom [5].
It is notable that in these cases $B_f$ has a discrete analogue, a sum of a Toeplitz and a Hankel matrix. An analogue of Szegő’s theorem for such matrices was proved by Johansson [18]. Furthermore, there is a counterpart of the Borodin–Okounkov formula for these matrices, a precise formula derived by Basor and Ehrhardt [4]. Another proof of this precise formula, which is similar to that of Borodin and Okounkov for Toeplitz matrices, was given by Betea [8], who used symplectic and orthogonal Schur measures.
Lastly, recall some results on the rate of convergence of additive functionals for matrix ensembles. For example, classical compact groups were studied by Johansson [18]. In [19] authors considered the traces of random Haar-distributed matrices multiplied by deterministic ones. Lambert, Ledoux and Webb [20] studied the rate of convergence of additive functionals in $\beta$-ensembles. Gaussian Laguerre and Jacobi ensembles were considered by Bufetov and Berezin [7], who used the Deift–Zhou asymptotic analysis of the Riemann–Hilbert problem [14]. The authors of [1] studied convergence in Wigner ensembles. However, all ensembles mentioned above are finite-dimensional. On the other hand, the Bessel kernel point process is a measure on infinite configurations. Therefore, we pursue the operator-theoretic approach used in [2] and [3].
§ 2. Statement of the result
For a function $f\in L_1(\mathbb{R})\cap L_\infty(\mathbb{R})$ we define its Fourier transform by
Here and below, if $f$ is defined on $\mathbb{R}_+$, then we extend it to an even function. In this case the Fourier transform coincides with the cosine transform. Recall that the Sobolev $p$-seminorm of a function $f$ on $\mathbb{R}_+$ is
For an integer $p$ we have $\|f\|_{\dot{H}_{p}}=(2\pi)^{-1/2}\|f^{(p)}\|_{L_2(\mathbb{R}_+)}$ by Parseval’s theorem. We set $\|f\|_{H_p} = \|f\|_{L_2(\mathbb{R}_+)} \!+\! \|f\|_{\dot{H}_{p}}$ and define the Sobolev space $H_p(\mathbb{R}_+) \!=\! \{{f\!\in\! L_2(\mathbb{R}_+)}$: $\|f\|_{H_p}\!<\!+\infty\}$ to be a Banach space with the norm $\|\,{\cdot}\,\|_{H_p}$. The Sobolev space $H_p(\mathbb{R})$ on the real line is defined similarly, and we keep the same notation for the norm.
We define the following seminorm of a function $f$ on $\mathbb{R}_+$:
Again, set $\|f\|_{\mathcal{B}}=\|f\|_{L_1(\mathbb{R}_+)} + \|f\|_{L_2(\mathbb{R}_+)}+\|xf'(x)\|_{L_\infty(\mathbb{R}_+)} +\|f\|_{\dot{\mathcal B}}$ and consider the space
where $\mathcal{F}$ is the unitary Fourier transform on $L_2(\mathbb{R})$: $\mathcal{F}h = \sqrt{2\pi}\,\widehat{h}$, $h\in L_2(\mathbb{R})$. Again, if $f$ is defined on $\mathbb{R}_+$, then we extend it to an even function.
Then $H_1(\mathbb{R}) = H_1^+(\mathbb{R})\oplus H_1^-(\mathbb{R})$. For any function $b\in H_1(\mathbb{R})$ let $b_\pm = P_\pm b$; $b = b_+ + b_-$. The above space $\mathcal{B}$ is clearly a subspace of $H_1(\mathbb{R})$ with embedding given by even continuation, so the decomposition is well defined on $\mathcal{B}$, although the components can fail to remain in $\mathcal{B}$. We refer to this decomposition as Wiener–Hopf factorization.
Our main result is then stated as follows.
Theorem 1. Let $b \in \mathcal{B}$. Then for any $R>0$
(For the first two identities, see § 6. For the last one, see [5], § II.) Indeed, the expression (2.3) for the remainder coincides with the one obtained by Basor, Ehrhardt and Widom in [5], formula (3).
We now proceed to a corollary for the determinantal point process with Bessel kernel $\mathbb{P}_{J_\nu}$ (see § 5). Let $\overline{S_f^R} = S_f^R - \mathbb{E}_{J_\nu}S_f^R$. We denote by $F_{R, f}$ and $F_{\mathcal{N}}$ the cumulative distribution functions of $\overline{S_f^R}$ and the standard Gaussian, respectively.
Theorem 2. Let $b\in\mathcal{B}$ be a real-valued function, satisfying $c_3^{\mathcal{B}}(b)=1/2$. Then there exists a constant $C=C(L(b), \|b_+\|_{L_\infty})$ providing the following estimate for the Kolmogorov–Smirnov distance for any $R\geqslant 1$
Step 2. It is shown below that the operator $W_{e^{-f_+}}W_{e^f}W_{e^{-f_-}} - I$ is trace class. We can apply the Jacobi–Dodgson identity $\det(PAP) = \det(A)\det(QA^{-1}Q)$ for orthogonal projections $P$ and $Q$ satisfying $P+Q=I$ and a determinant-class invertible operator $A$. Then we have
where $W_{e^f}^{-1} = W_{e^{-f_+}}W_{e^{-f_-}}$, as follows from the properties of Wiener–Hopf factorization.
Step 3. From the properties of the Fredholm determinant and Wiener–Hopf factorization we obtain $\det(W_{e^{-f_+}}W_{e^f}W_{e^{-f_-}}) = \det(W_{e^{f_-}}W_{e^{f_+}}W_{e^{-f_-}}W_{e^{-f_+}})$. Widom’s formula (see [27]) states that determinant on the right is well defined and its value is
By contrast to the proof of Theorem 3, it is not known if $W_{e^{-b_+}}B_{e^{b}}W_{e^{-b_-}}-I$ is trace class, so the Jacobi–Dodgson identity cannot be applied directly. However, one can observe that equality (3.2) means that the ratio of the determinants in the left-hand side does not depend on $R$. In our case we can use an analogue of the Jacobi–Dodgson identity (Corollary 2) to show the independence of this ratio of $R$.
Lemma 2. Let $b\in H_1(\mathbb{R}_+)\cap L_1(\mathbb{R}_+)$, and let $\chi_{[R,\infty)}(B_{e^{-b}} - W_{e^{-b}})\chi_{[R, \infty)}$ be a trace-class operator for each $R>0$. Then there exists $Z(b)$ such that for all $R>0$ we have
Observe that the denominator in (3.5) equals $Q_R^{\mathcal{B}}(b)$. The central step in the proof of Theorem 1 is as follows.
Lemma 3. For $b\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$ such that $\|b\|_{\dot{\mathcal B}}< \infty$ and $z\in\mathbb{C}$ the operator $\chi_{[R, \infty)}(B_{b+z} - W_{b+z})\chi_{[R, \infty)}$ is trace class for any $R>0$. There exists a constant $C$ such that for any $R\geqslant1$, $z\in\mathbb{C}$ and $b$ satisfying the conditions above we have
The operator $Q_R^{\mathcal{B}}(b)$ is related to the operator in Lemma 3 by formula (8.5) below. This allows us to obtain estimate (2.2). Further, the asymptotic for the numerator in (3.5) follows from Lemma 1 and an asymptotic result of Basor and Ehrhardt (see Theorem 6). This proves that $Z(b) = \exp(c_2^{\mathcal{B}}(b) + c_3^{\mathcal{B}}(b))$.
The rest of the paper has the following structure. In § 4 we recall some facts and notation for trace-class and Hilbert–Schmidt operators. Then in § 5 we explain the relation between the Fredholm determinant in Theorem 1 and the determinantal point process with Bessel kernel. In § 6 we establish some properties of Wiener–Hopf factorization and deduce Lemma 1. Section 7 presents the proof of Lemma 3. We complete the proof of Theorems 1 and 2 in § 8.
§ 4. Trace class and Hilbert–Schmidt operators
Here and below we denote the ideals of trace-class and Hilbert–Schmidt operators by $\mathcal{J}_1$ and $\mathcal{J}_2$, respectively. For more details of the theory, we refer the reader to [24] and [22]. We let $E$ denote $\mathbb{R}$ or $\mathbb{R}_+ = [0, \infty)$.
Definition 1. The operator $K$ on $L_2(E)$ is locally trace class if for each bounded measurable subset $B$ of $ E$ the operator $\chi_B K \chi_B$ is trace class. Let $\mathcal{J}_1^{\mathrm{loc}}(L_2(E))$ denote the space of locally trace-class operators.
Recall that for $K\in \mathcal{J}_1$ the Fredholm determinant $\det (I+K)$ is well defined, which is equal to a product, $\det(I+K)=\prod_{\lambda\in\sigma(K)}(1+\lambda)$, by Lidskii’s theorem. It is continuous with respect to the trace norm. We also define the regularized determinant by
$$
\begin{equation*}
\det{}_2 (I + K)=\exp(-\operatorname{Tr}(K))\det(I+K).
\end{equation*}
\notag
$$
It is continuous with respect to the Hilbert–Schmidt norm, so it can extended by continuity to all operators in $\mathcal{J}_2$.
The following theorem implies that for a function $f\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$ the determinants $\det (I + \chi_{[0, R]}B_f\chi_{[0, R]})$ and $\det(I + \chi_{[0, R]}W_f\chi_{[0, R]})$ are well defined.
Theorem 4 ([22], Theorem 3.11.9). Let $K$ be a continuous kernel on $[a, b]^2$ inducing a positive self-adjoint operator $K$ on $L_2[a, b]$. Then the operator $K$ is trace class. In addition, for any trace-class $K$ with continuous kernel we have
We refer to (4.1) as the calculation of the trace over the diagonal.
Recall that the kernel of an integral operator is defined almost everywhere on $[a, b]^2$. Thus, for an operator satisfying the assumptions of Theorem 4, formula (4.1) can fail to hold after we change the kernel on the diagonal, although the induced operator remains the same. This explains the assumption of continuity. Let us now extend the class of kernels for which we can calculate the trace over the diagonal.
Let $K$ be a self-adjoint trace-class operator. By the Hilbert–Schmidt theorem we have a spectral decomposition, so that its kernel can be chosen in the form
which is equal to $\displaystyle\int_a^b K(x, x)\,dx$ if we choose $K$ as in (4.2). Observe that this does not require the kernel to be continuous, but this choice is better defined: we can only change $K$ on subsets of $[a, b]^2$ having projections with zero measure. Choosing similarly the kernel of an operator $fK$:
we see that then $\displaystyle\operatorname{Tr}(fK)=\int_a^b f(x)K(x, x)\, dx$ in a similar way. We use these observations in § 5.
We will also use the following theorem.
Theorem 5 ([15], Theorem 2.2). For $A, B \in \mathfrak{B}(\mathcal{H})$ such that $[A,B]\in \mathcal{J}_1$ we have $e^Ae^Be^{-A-B}- I\in \mathcal{J}_1$. For the Fredholm determinant we have
§ 5. A link between the Bessel operator and the determinantal point process with Bessel kernel
Let $E$ be $\mathbb{R}$ or $\mathbb{R}_+ = [0, \infty)$. Recall that a configuration on $E$ is a discrete subset of $E$. We denote the set of configurations by $\operatorname{Conf}(E)$. It can be endowed with a $\sigma$-algebra of measurable subsets $\mathfrak{X}$ (for instance, see [25]). A point process $\mathbb{P}$ is defined to be a probability measure on $(\operatorname{Conf}(E), \mathfrak{X})$.
For a measurable function $f\colon E\to\mathbb{C}$ we define two measurable functions on configurations by the following formulae:
they are called additive and multiplicative functionals, respectively.
Definition 2. A point process $\mathbb{P}_K$ on $\operatorname{Conf}(E)$ is determinantal if there exists an operator $K\in \mathcal{J}_1^{\mathrm{loc}}(L_2(E))$ satisfying the following condition: for any bounded measurable function $f$ with compact support $B=\operatorname{supp} f$
$$
\begin{equation}
\mathbb{E}_K \Psi_{1+f}=\det (I + fK\chi_B).
\end{equation}
\tag{5.1}
$$
We consider the limit distribution of $S_f^R = S_{f(x/R)}$ as $R\to\infty$ for the determinantal point process with Bessel kernel $\mathbb{P}_{J_\nu}$ on $\mathbb{R}_+$ [26] induced by the operator $B_{\chi_{[0, 1]}}$.
Theorem 3 describes the limit under consideration for the sine process $\mathbb{P}_{\mathcal{S}}$ on $\mathbb{R}$ induced by the operator $\mathcal{F}^{-1}\chi_{[-\pi, \pi]} \mathcal{F}$ on $L_2(\mathbb{R})$. To be precise, we formulate the following statement.
Proposition 1. For $f\in L_1(E)\cap L_\infty(E)$ we have
By Proposition 1 it is equal to the Fredholm determinant $\det(\chi_{[0, R]}B_{e^{\lambda b}}\chi_{[0, R]})$ for $\mathbb{P}_{J_\nu}$ and to $\det(\chi_{[0, 2\pi R]} W_{e^{\lambda b}}\chi_{[0, 2\pi R]})$ for the sine process.
To prove Proposition 1 we follow Bufetov (see [13], § 2.9). We extend the class of functions for which (5.1) holds. Let $K$ be a Hilbert–Schmidt self-adjoint operator on $L_2(E)$. We choose its kernel $K$ to coincide with its spectral decomposition (see § 4). We set
if the integral of the diagonal exists. Observe, that this definition does not require $K$ to be trace class. However, if it is trace class, then we obtain the ordinary Fredholm determinant.
Let $\Pi\in \mathcal{J}_1^{\mathrm{loc}}(L_2(E))$ be an orthogonal projection. By the Macchi–Soshnikov theorem (see [21] and [25]) $\Pi$ defines a point process $\mathbb{P}_\Pi$. We choose the kernel $\Pi(x, y)$ to coincide with the spectral decomposition of $\chi_B \Pi \chi_B$ for every compact subset $B$ of $ E$. Consider the measure $d\mu_\Pi(x)=\Pi(x, x)\,dx$ on $E$.
Proposition 2. If $f\in L_1(E, d\mu_\Pi)\cap L_\infty(E, d\mu_\Pi)$, then $\Psi_{1+f}\in L_1(\operatorname{Conf}(E), \mathbb{P}_\Pi)$ and
Proof. The formula holds for $f\in L_\infty(E)$ with compact support The exponential of the integral of the diagonal in (5.4) is continuous with respect to the norm in $L_1(E, d\mu_\Pi)$. The regularized determinant is continuous with respect to the Hilbert–Schmidt norm, which is equal to
and therefore it is continuous with respect to the $L_2(E, d\mu_\Pi)$-norm. Thus, multiplicative functionals define a continuous mapping of a dense subset of $L_1(E, d\mu_\Pi)\cap L_\infty(E, d\mu_\Pi)$ and therefore of the whole space. This extension coincides by continuity with the expectation of the multiplicative functionals by Beppo Levi’s theorem, since one can approximate any multiplicative functional $\Psi_{1+|f|}$, $f\in L_1(E, d\mu_\Pi)\cap L_\infty(E, d\mu_\Pi)$ by the functionals $\Psi_{1+|f_n|}$, $|f_n|=\chi_{[-n, n]}|f|$. Hence (5.5) is an extension of (5.1) by continuity. The proposition is proved.
Proof of Proposition 1. The proof for the sine process case can be found in [13], Lemma 3.3. We proceed differently by proving the equality between the extended determinants of $I + f\mathcal{F}^{-1}\chi_{[-\pi, \pi]}\mathcal{F}$ and $I + fB_{\chi_{[0, 1]}}$ and the Fredholm determinants of $I + \chi_{[0, 2\pi]}W_f\chi_{[0, 2\pi]}$ and $I + \chi_{[0, 1]}B_f\chi_{[0, 1]}$, respectively. Note that by Theorem 4 we can take the kernel of $B_{\chi_{[0, 1]}}$ to be equal to (1.1), and the kernel of $\mathcal{F}^{-1}\chi_{[-\pi,\pi]}\mathcal{F}$ to be $\sin(\pi(x-y))/(\pi(x-y))$. It is straightforward to see in both cases that the integrals over the diagonals of the kernels coincide.
Now we show that the regularized determinants coincide. This is a consequence of unitary invariance. Recall that $B_{\chi_{[0,1]}}$ is unitarily equivalent to $\chi_{[0, 1]}$ via the Hankel transform $H_\nu$ (see § 8 and formula (8.4)). Then we have the unitary equivalence
where the last operator has the same regularized determinant as $\chi_{[0, 1]}B_f\chi_{[0, 1]}$. Notably, the operator obtained is trace class and its Fredholm determinant is well defined.
We turn to the proof for the sine kernel. We use the translation operator $T_tf(x) = f(x+t)$ to proceed as follows:
where the determinant of the obtained operator coincides, as above, with the one for $\chi_{[0, 2\pi]}W_f\chi_{[0, 2\pi]}$. Again, $\chi_{[0, 2\pi]}W_f\chi_{[0, 2\pi]}$ is trace class.
In both cases conjugation by the unitary operator $U_R f(x)=R^{-1/2}f(x/R)$ completes the proof.
The proposition is proved.
§ 6. Wiener–Hopf factorization on $H_1( {\mathbb{R}})$
Recall that the Wiener–Hopf operator as a map $\mathbf{W}\colon f\mapsto W_f$, $L_\infty(\mathbb{R})\to \mathfrak{B}(L_2(\mathbb{R}_+))$ is not a Banach algebra homomorphism. But it can be restricted to some subalgebras on which it preserves multiplication.
Recall that $H_1(\mathbb{R})$ is a Banach algebra. The spaces
are Banach subalgebras since $\operatorname{supp}(\widehat{f}*\widehat{g}) \subset \operatorname{supp}\widehat{f} + \operatorname{supp}\widehat{g}$. The operators $P_\pm = \mathcal{F}^{-1} \chi_{\mathbb{R}_\pm}\mathcal{F}$ map $H_1(\mathbb{R})$ into itself and $P_\pm H_1(\mathbb{R})=H_1^\pm(\mathbb{R})$. Then $\mathbf{W}$, as restricted to $H_1^\pm(\mathbb{R})$, is a Banach algebra homomorphism. As before, we consider $H_1(\mathbb{R}_+)$ to be embedded in $H_1(\mathbb{R})$ by an even continuation, which defines a decomposition for functions on $\mathbb{R}_+$.
be the Banach algebras $H_1^\pm(\mathbb{R})$ with adjoined unit. Then for $b\in H_1^\pm\oplus\mathbb{C}$ and an entire function $\Phi=\sum_{i=0}^\infty \Phi_iz^i$ the element
Recall that for functions in the spaces under consideration the Wiener–Hopf operators have kernels.
Lemma 4 ([3], Proposition 5.2). Let $a$ be a function on $\mathbb{R}$ such that $\|a\|_{H_{1/2}} + \|\widehat{a}\|_{L_1}<\infty$. Then $\mathcal{F}^{-1}a\mathcal{F}$ is an integral operator on $L_2(\mathbb{R})$ with kernel $k(x, y) = \widehat{a}(x-y)$.
Proposition 3. 1. The map $\boldsymbol{W}\big|_{H_1^\pm(\mathbb{R})}\colon f \mapsto W_f$ defines Banach algebra homomorphisms $H_1^\pm(\mathbb{R})\oplus \mathbb{C}\to \mathfrak{B}(L_2(\mathbb{R}_+))$.
2. For $b_{\pm}\in H_1^\pm(\mathbb{R})\oplus \mathbb{C}$ we have
3. For $b_\pm\in H_1^\pm(\mathbb{R})\oplus\mathbb{C}$ we have $W_{b_-}W_{b_+}=W_{b_-b_+}$.
In particular, $W_{e^{b_\pm}}=\exp(W_{b_\pm})$ and $W_{e^b}=\exp(W_{b_-})\exp(W_{b_+})$.
Proof. It is sufficient to prove the statements for $b_\pm\in H_1^\pm(\mathbb{R})$. We give proofs for $H_1^+(\mathbb{R})$; the case of $H_1^-(\mathbb{R})$ is proved similarly.
1. By Lemma 4 Wiener–Hopf operators have kernels. Let $a$ and $b$ be functions in $H_1^+(\mathbb{R})$. Since $\widehat{b}_\pm(x) = \chi_{\mathbb{R}_\pm}(x)\widehat{b}_\pm(x)$, we have
where the identity $\chi_{\mathbb{R}_+}(y)\chi_{\mathbb{R}_+}(t-y) = \chi_{\mathbb{R}_+}(y)\chi_{\mathbb{R}_+}(t-y)\chi_{\mathbb{R}_+}(t)$ has been used.
2. Since $\chi_{[0, R]}(x)\chi_{\mathbb{R}_+}(x-y) = \chi_{[0, R]}(x)\chi_{\mathbb{R}_+}(x-y)\chi_{[0, R]}(y)$ for all $(x, y) \in \mathbb{R}_+$, for the kernel of $(\chi_{[0, R]}W_{b_+})(x, y)$ we can write
By Theorem 4 the operator $\chi_{[0, R]}W_b\chi_{[0, R]}$ is trace class for $b\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$, so the left-hand side of (6.1) is equal to
where $\operatorname{Tr}(\chi_{[0, R]}W_b\chi_{[0, R]}) = R\widehat{b}(0)$ is calculated over the diagonal. A direct calculation shows that the operator $\chi_{[0, R]}W_f\chi_{[0, R]}$ is Hilbert–Schmidt for $f\in H_1(\mathbb{R}_+)$. Hence the commutator $[\chi_{[0, R]} W_{b_-}\chi_{[0, R]}, \chi_{[0,R]}W_{b_+}\chi_{[0, R]}]$ is trace class, and its trace is zero. Theorem 5 gives
In this section we establish Lemma 3. Our calculations follow Basor and Ehrhardt (see [3], Lemmas 2.6–2.8). The section is structured as follows. In Propositions 4 and 5 and Corollary 1 we establish estimates for the trace norms of integral operators with kernels of certain form. Then we recall the necessary properties of Bessel functions and their derivatives. We end the preparations to the proof by Lemma 5, establishing some properties of the function spaces in question. The rest of the section is devoted to the proof of Lemma 3 itself.
Proposition 4. Let $K$ be an integral operator on $L_2[R, \infty)$ with kernel
where Fubini’s theorem can be applied since the function under the integral is absolutely integrable by the assumption that the right-hand side of (7.1) is finite. Recall that $\mathfrak{B}(L_2[R, \infty)) \simeq \mathcal{J}_1^*(L_2[R, \infty))$, where $A\in\mathfrak{B}(L_2[R, \infty)$ acts on $\mathcal{J}_1^*(L_2[R, \infty))$ by the rule $X\mapsto \operatorname{Tr}(AX)$. For the norm of $K$, considering it as a functional on $\mathfrak{B}(L_2[R, \infty))$ and using the definition of a trace and the expression (7.2), we have
where $\{f_i\}_{i\in\mathbb{N}}$ is an arbitrary orthonormal basis in $L_2[R, \infty)$. Next we use the Cauchy–Bunyakovsky–Schwarz inequality to obtain
Proposition 5. Let the kernel $K(x, y)$ of an integral operator $K$ be absolutely continuous with respect to $y$ on $[a, b]$. Assume that $K = K\chi_{[a, b]}$. Consider the operator $\partial_yK$ with kernel ${\partial K(x, y)}/{\partial y}$ and assume that $\partial_y K\in\mathcal{J}_2$. Then
Proof. Let $Pf(x)=(1/(b-a))\chi_{[a, b]}(x)\langle \chi_{[a, b]}, f\rangle_{L_2}$ be a one-dimensional projector onto $\chi_{[a, b]}$, and let $Q=\chi_{[a, b]}-P$. We decompose the operator: $K = KP + KQ$. The operator $KP$ is also a one-dimensional projector with trace norm of at most
Recall that the Volterra operator $V$ on $L_2([a, b])$ is defined by the formula $\displaystyle Vf(x)=\int_a^x f(t)\,dt$. It is a Hilbert–Schmidt operator with kernel $V(x, y) = \chi_{[a, x]}(y)$, and its Hilbert–Schmidt norm is
by the definition of $V$ and $Q$. The inequality $\|\partial_yKVQ\|_{\mathcal{J}_1}\leqslant \|\partial_yK\|_{\mathcal{J}_2}\,\|V\|_{\mathcal{J}_2}$ finishes the proof.
Corollary 1. There exists a constant $C$ such that for any kernel $K$ of an integral operator on $L_2(\mathbb{R}_+)$ with the estimates
Recall several properties of Bessel functions. Set $\mathfrak{J}(x) = \sqrt{x}\,J_\nu(x)$ and $\mathfrak{D}(x) = \mathfrak{J}_\nu(x)-\sqrt{2/\pi}\cos(x-\phi_\nu)$, where $\phi_\nu={\pi}/{4} + ({\pi}/{2})\nu$. We have the following asymptotics for the Bessel function and its derivative as $x$ approaches infinity:
Lastly, we establish some properties of the function spaces we have introduced.
Lemma 5. There exists a constant $C$ such that for any $a\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$ satisfying $\|a\|_{\dot{\mathcal B}}<\infty$ we have
Since $ta''(t) = (at)'' - 2a'$, the second term is estimated by $2\|a\|_{\dot{H}_{1}} + \|ta(t)\|_{\dot{H}_{2}}$.
(2) The proof is quite similar to the previous one.
(3) The statement on $a$ follows from $\widehat{a}\in L_1(\mathbb{R})$ and the Riemann–Lebesgue lemma.
Since $a'\in H_1(\mathbb{R}_+)$ we see that $\widehat{a'}$ is absolutely integrable and the statement on $a'$ follows again from the Riemann–Lebesgue lemma.
By the first and second assertions $(ta''(t))' = a''(t) + ta'''(t)$ is absolutely integrable, so $ta''(t)$ tends to a finite limit as $t$ approaches infinity. Since $(ta(t))''$ and $a'(t)$ are square integrable, so is $ta''(t)$, which yields that the limit is zero.
(4) The expression of $a'(0)$ in terms of the cosine transform yields
We devote the rest of this section to the proof of Lemma 3.
Proof of Lemma 3. Since $B_1 - W_1 = 0$, it is sufficient to prove the statement for $b\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$ satisfying $\|b\|_{\dot{\mathcal B}}<\infty$.
Recall the formula for the kernel of the difference $\mathcal{R}_b(x, y) = B_b(x, y) - W_b(x, y)$:
Let us outline the plan of the proof. One can observe from the asymptotics (7.5) and (7.6) that the above integral contains an asymptotically small function for large $x$ and $y$. Thus, we substitute $\mathfrak{J}(x) = \sqrt{x}\,J_\nu(x) = \mathfrak{D}(x) + \sqrt{2/\pi}\cos(x-\phi_\nu)$ into the expression for $\mathcal{R}_b(x, y)$. This substitution and several integrations by parts represent the kernel as a sum of other kernels, to which Corollary 1 and Proposition 4 can be applied.
To be precise, we perform a sequence of decompositions of our kernel, for which we use the following notation:
Explicit formulae are presented below (see (7.13)–(7.15)). We prove an estimate for the trace norm for each of the above kernels separately. In particular, the estimate of the trace norm for $\chi_{[R, \infty)} \mathcal{R}_1\chi_{[R, \infty)}$ follows from Corollary 1. Estimates for the trace norms of the operators with the kernels
Next we substitute the expression (7.9) for $\sin(2\phi_\nu)/(\pi(x+y))$ into the above identity. These calculations prove the decomposition (7.10) for the following kernels $\mathcal{R}_1$ and $\mathcal{R}_2$:
The first term induces the one-dimensional projector $h_1(x)\langle h_2(y), -\rangle_{L_2}$ with trace norm bounded by $\|h_1\|_{L_2}\|h_2\|_{L_2}$. By Lemma 5, for the trace norm of the corresponding operator on $L_2[R, \infty)$ we have $|A_\nu b'(0)|/R \leqslant |A_\nu|\|b\|_{\dot{\mathcal B}}/R$.
We apply Proposition 4 to the second term. This results in the estimate
Proposition 7 ([23], Proposition 6.2.9). Let $A$ be an invertible operator on a separable Hilbert space $\mathcal{H}$ and $A-I\in \mathcal{J}_1(\mathcal{H})$. Let $P$ be an orthogonal projection in $\mathcal{H}$ and $Q=I-P$. Then
if $A$ is invertible, $A-I\in\mathcal{J}_2(\mathcal{H})$ and all Fredholm determinants are well defined.
Remark 3. It is not required that $A-I\in\mathcal{J}_1(\mathcal{H})$ in this statement.
Proof of Corollary 2. Let $A = I + K$, where $K\in\mathcal{J}_2(\mathcal{H})$. We choose a joint orthogonal basis $\{e_i\}_{i\in\mathbb{N}}$ of eigenvectors of $P_1$ and $P_2$. Let $R_n$ be the orthogonal projection onto $\operatorname{span}(\{e_i\}_{i\in 1,\dots,n})$. We regard $R_n\mathcal{H}$ as a separable Hilbert space. Since the $P_i$ and $Q_i$ commute with $R_n$, they remain orthogonal projections in $R_n\mathcal{H}$.
Since $K$ is compact, $R_nKR_n \to K$ in the operator norm. This implies that for all sufficiently large $n$ the operator $I + R_n K R_n$ is invertible on $R_n\mathcal{H}$. Applying Proposition 7 we have
where $(R_n + R_nKR_n)^{-1}\to (I+K)^{-1}$ and $R_n\to I$ strongly, and $(R_nKR_n)^2\to K$ and $R_nK^2\to K^2$ in the trace norm since $K\in\mathcal{J}_2(\mathcal{H})$. Therefore, we have
Equality (8.2) follows by taking limits in (8.3) as $n\to\infty$.
The corollary is proved.
Let us extend the definition (1.1) to all functions in $L_\infty (\mathbb{R}_+)$. Recall that the Hankel transform $H_\nu$ on $L_2(\mathbb{R}_+)$ is defined by
Similarly to the Fourier transform, it is an isometry of $L_1(\mathbb{R}_+)\cap L_2(\mathbb{R}_+)$ and extends to $L_2(\mathbb{R}_+)$ by continuity. Also recall that $H_\nu^* = H_\nu^{-1} = H_\nu$. Now for $b\in L_\infty(\mathbb{R}_+)$ we define the Bessel operator by
$$
\begin{equation}
B_b=H_\nu b H_\nu.
\end{equation}
\tag{8.4}
$$
This coincides with the definition (1.1) for $b\in L_1(\mathbb{R}_+)\cap L_\infty(\mathbb{R}_+)$, but from (8.4) it is clear that $\mathbf{B}\colon f\to B_f$, $L_\infty(\mathbb{R}_+)\to \mathfrak{B}(L_2(\mathbb{R}_+))$ is a Banach algebra homomorphism and, in particular, $B_{e^b}^{-1} = B_{e^{-b}}$.
Next we want to show that for $b\in \mathcal{B}$ the operator $B_b - W_b$ is Hilbert–Schmidt. By Lemma 3 we have $\chi_{[1, \infty)}(B_b - W_b)\chi_{[1, \infty)}\in \mathcal J_1$ and $\chi_{[0, 1]}(B_b - W_b)\chi_{[0, 1]}\in \mathcal J_1$ by Theorem 4. The following statement, completing the argument, can be verified directly.
Lemma 6 (see [3], Lemmas 3.1 and 3.2). For $b\in H_1(\mathbb{R}_+)\cap L_1(\mathbb{R}_+)$ the operators $\chi_{[0, 1]}B_b$ and $\chi_{[0, 1]}W_b$ are Hilbert–Schmidt.
Proof of Lemma 2. For arbitrary $R_1, R_2 >0$ we set $P_1 = \chi_{[0, R_1]}$, $P_2 = \chi_{[0, R_2]}$ and $A = W_{e^{-b_+}}B_{e^b}W_{e^{-b_-}}$. By the first statement of Proposition 3 and the definition (8.4) $A$ is invertible and its inverse is $W_{e^{b_-}}B_{e^{-b}}W_{e^{b_+}}$. In the notation of Corollary 2 we also have $Q_1 = \chi_{[R_1, \infty)}$ and $Q_2 = \chi_{[R_2, \infty)}$. Using the second and third assertions of Proposition 3 we can write
Also recall that for any trace-class operator $A$ and bounded $B$ we have $\|AB\|_{\mathcal{J}_1} \leqslant \|A\|_{\mathcal{J}_1}\|B\|$. By the definition of Wiener–Hopf operators $\|W_b\| = \|b\|_{L_\infty}$. These facts and Lemmas 3 and 7, in combination with the expression (8.5), prove estimate (2.2).
By the above argument the denominator in (3.5) approaches $1$. From Lemma 1 and Theorem 6, for the numerator in (3.5) we have
It is now clear that $Z(b) = \exp(c_2^{\mathcal{B}}(b) + c_3^{\mathcal{B}}(b))$. This finishes the proof.
We now proceed to the proof of Theorem 2. First we establish an estimate for the rate of convergence of $\mathbb{E}_{J_\nu}S_f^R$. The argument is based on the proof of Proposition 5.1 in [3].
Lemma 8. There exists a constant $C$ such that for any $b\in\mathcal{B}$ we have
It is clear that $F(\xi) = O(\xi^{-1})$ as $\xi\to\infty$. Therefore, we have $|F(\xi)| \leqslant C(1+\xi)^{-1}$ for some constant $C$. Integrating by parts in the formula for the expectation we obtain
which is at most ${C\|b\|_{\dot{H}_{1}}}/{\sqrt{R}}$. This completes the proof.
Proof of Theorem 2. Recall that we denote by $F_{R, b}$ and $F_{\mathcal{N}}$ the cumulative distribution functions of $\overline{S_b^R}$ and the standard Gaussian distribution $\mathcal{N}(0, 1)$, respectively. Using Feller’s smoothing estimate (see [16], Ch. XVI, § 3), Theorem 1 and Proposition 1, for any $T>0$ we have
By the second statement of Theorem 1 and Lemma 8 the integrand can be estimated by the following expression for some $C>0$ depending only on $\|b_+\|_{L_\infty}$ and $L(b)$:
Observe that if $|k| \leqslant C_1\ln R$, where $C_1C <1/2$, then there exists a constant $\widetilde{C}$ such that for any $R\geqslant 1$ the above expression is at most
Zhigang Bao and Yukun He, Quantitative CLT for linear eigenvalue statistics of Wigner matrices, arXiv: 2103.05402
2.
E. L. Basor and Yang Chen, “A note on Wiener–Hopf determinants and the Borodin–Okounkov identity”, Integral Equations Operator Theory, 45:3 (2003), 301–308
3.
E. L. Basor and T. Ehrhardt, “Asymptotics of determinants of Bessel operators”, Comm. Math. Phys., 234:3 (2003), 491–516
4.
E. L. Basor and T. Ehrhardt, Determinant computations for some classes of Toeplitz–Hankel matrices, arXiv: 0804.3073
5.
E. L. Basor, T. Ehrhardt and H. Widom, On the determinant of a certain Wiener–Hopf + Hankel operator, arXiv: math/0304002
6.
E. L. Basor and H. Widom, “On a Toeplitz determinant identity of Borodin and Okounkov”, Integral Equations Operator Theory, 37:4 (2000), 397–401
7.
S. Berezin and A. I. Bufetov, “On the rate of convergence in the central limit theorem for linear statistics of Gaussian, Laguerre, and Jacobi ensembles”, Pure Appl. Funct. Anal., 6:1 (2021), 57–99
8.
D. Betea, Correlations for symplectic and orthogonal Schur measures, arXiv: 1804.08495
9.
A. Böttcher, “On the determinant formulas by Borodin, Okounkov, Baik, Deift, and Rains”, Toeplitz matrices and singular integral equations, Oper. Theory Adv. Appl., 135, Birkhäuser Verlag, Basel, 2002, 91–99
10.
A. Böttcher and B. Silbermann, Analysis of Toeplitz operators, Springer Monogr. Math., 2nd ed., Springer-Verlag, Berlin, 2006, xiv+665 pp.
11.
A. Borodin and A. Okounkov, “A Fredholm determinant formula for Toeplitz determinants”, Integral Equations Operator Theory, 37:4 (2000), 386–396
12.
A. Bufetov, “The expectation of a multiplicative functional under the sine-process”, Funct. Anal. Appl., 58:2 (2024), 120–128
13.
A. I. Bufetov, The sine-process has excess one, arXiv: 1912.13454
14.
P. Deift and X. Zhou, “A steepest descent method for oscillatory Riemann–Hilbert problems. Asymptotics for the MKdV equation”, Ann. of Math. (2), 137:2 (1993), 295–368
15.
T. Ehrhardt, “A generalization of Pincus' formula and Toeplitz operator determinants”, Arch. Math. (Basel), 80:3 (2003), 302–309
16.
W. Feller, An introduction to probability theory and its applications, v. 2, John Wiley & Sons, Inc., New York–London–Sydney, 1966, xviii+636 pp.
17.
I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series, and products, 5th ed., Academic Press, Inc., Boston, MA, 1994, xlviii+1204 pp.
18.
K. Johansson, “On random matrices from the compact classical groups”, Ann. of Math. (2), 145:3 (1997), 519–545
19.
J. P. Keating, F. Mezzadri and B. Singphu, “Rate of convergence of linear functions on the unitary group”, J. Phys. A, 44:3 (2011), 035204, 27 pp.
20.
G. Lambert, M. Ledoux and C. Webb, “Quantitative normal approximation of linear statistics of $\beta$-ensembles”, Ann. Probab., 47:5 (2019), 2619–2685
21.
O. Macchi, “The coincidence approach to stochastic point processes”, Adv. in Appl. Probab., 7 (1975), 83–122
22.
B. Simon, Operator theory, Compr. Course Anal., 4, Amer. Math. Soc., Providence, RI, 2015, xviii+749 pp.
23.
B. Simon, Orthogonal polynomials on the unit circle, Part 1. Classical theory, Amer. Math. Soc. Colloq. Publ., 54, Part 1, Amer. Math. Soc., Providence, RI, 2005, xxvi+466 pp.
24.
B. Simon, Trace ideals and their applications, Math. Surveys Monogr., 120, 2nd ed., Amer. Math. Soc., Providence, RI, 2005, viii+150 pp.
25.
A. Soshnikov, “Determinantal random point fields”, Russian Math. Surveys, 55:5 (2000), 923–975
26.
C. A. Tracy and H. Widom, “Level spacing distributions and the Bessel kernel”, Comm. Math. Phys., 161:2 (1994), 289–309
27.
H. Widom, “A trace formula for Wiener–Hopf operators”, J. Operator Theory, 8:2 (1982), 279–298
Citation:
S. M. Gorbunov, “Rate of convergence in the central limit theorem for the determinantal point process with Bessel kernel”, Sb. Math., 215:12 (2024), 1607–1632
\Bibitem{Gor24}
\by S.~M.~Gorbunov
\paper Rate of convergence in the central limit theorem for the determinantal point process with Bessel kernel
\jour Sb. Math.
\yr 2024
\vol 215
\issue 12
\pages 1607--1632
\mathnet{http://mi.mathnet.ru/eng/sm10137}
\crossref{https://doi.org/10.4213/sm10137e}