Quantitative uniform resolvent estimates

Piero D’Ancona,^a Jérémy Faupin,^b and David Krejčiřík ^c

(

$a)$

Department of Mathematics “Guido Castelnuovo”, University of Rome “La Sapienza”, Piazzale Aldo Moro 5, 00185 Rome, Italy; dancona@mat.uniroma1.it.

$b)$

Université de Lorraine, CNRS, IECL, F-57000 Metz, France; jeremy.faupin@univ-lorraine.fr.

$c)$

Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Trojanova 13, 12000 Prague 2, Czechia; david.krejcirik@fjfi.cvut.cz.

24 June 2026)

Abstract

We derive quantitative uniform resolvent estimates for Schrödinger operators on the half-line with inverse-square potentials, which provide a sharp behaviour in the limit of large coupling. Our approach is based on a matrix representation of the boundary value of a weighted resolvent. The partial wave decomposition then turns these one-dimensional channel estimates into explicit weighted resolvent estimates for the Laplacian, its inverse-square potential perturbations and for the magnetic Laplacian with an Aharonov–Bohm potential. We also obtain exact Simon-type identities for the imaginary parts of the weighted resolvents of these operators.

Keywords:

Uniform resolvent estimates; Bessel operators; Hankel transform; inverse square potentials; Aharonov–Bohm potential; Kato smoothness.
MSC (2020):

Primary 35J10, 35P25; Secondary 33C10, 47A10, 47A55, 81Q10.

1 Introduction

Let $-\Delta$ be the self-adjoint realisation of the Laplacian in $L^{2}(\mathbb{R}^{d})$ with $d\geq 3$ . The resolvent cannot be bounded uniformly near its spectrum $[0,\infty)$ , since

\|(-\Delta-z)^{-1}\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}=\frac{1}{\mathrm{dist}(z,[0,\infty))}.

(1.1)

However, it is also well known that a uniform resolvent estimate can be achieved when the Laplacian is reconsidered as an operator between weighted spaces. For instance, one has

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}<\infty.

(1.2)

Here the weight $|x|^{-1}$ is scale invariant for the Laplacian and it is naturally associated with the classical Hardy inequality. Estimate (1.2) is behind the limiting absorption principle and Kato smoothness and smoothing estimates [19, 20, 31]. It has a perturbative interpretation as well: Kato’s abstract theory and the Birman–Schwinger principle convert such bounds into spectral stability statements for small (possibly non-self-adjoint) perturbations subordinated to the same critical weight [19, 16, 14, 18].

More specifically, potentials of the form $c\,|x|^{-2}$ with $c\in\mathbb{C}$ represent critical perturbations of the Laplacian, in view of the optimal Hardy inequality $-\Delta\geq c_{d}\,|x|^{-2}$ with $c_{d}:=(d-2)^{2}/4$ . The related critical electromagnetic Helmholtz, Strichartz and smoothing estimates, as well as spectral bounds, have been intensively studied in recent years [4, 5, 1, 3, 2, 24]. What is more, even if the classical Hardy inequality fails in two dimensions, there exist magnetic Hardy inequalities [23, 6] and the inverse-square potential is particularly relevant for the singular Aharonov–Bohm field [27, 22, 15].

The purpose of this paper is to make (1.2) quantitative and to identify the mechanism behind the constants. Let us now explain the main idea of this paper and present our main results.

As a first step, we note that a standard maximum principle reduction and scaling invariance (see Section 3.2 below for the detailed argument) yield

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}=\left\||x|^{-1}(-\Delta-(1\pm i0))^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}.

(1.3)

In this way, the problem is reduced to limiting values at a single point of the positive real axis.

1.1 Reduction to one-dimensional Bessel operators

The main insight of this paper is that the mechanism to compute the right-hand side of (1.3) is one-dimensional. Indeed, after a decomposition into spherical harmonics, every angular momentum channel is governed by a Bessel operator

h_{\nu}:=-\partial_{r}^{2}+\frac{\nu^{2}-1/4}{r^{2}}\qquad\hbox{in}\qquad L^{2}((0,\infty))=:L^{2}\,,

(1.4)

introduced as the Friedrichs extension of the operator initially defined on $C_{0}^{\infty}((0,\infty))$ . The central object of this paper is therefore the weighted boundary value

A_{\nu}^{\pm}:=r^{-1}(h_{\nu}-(1\pm i0))^{-1}r^{-1}.

(1.5)

The spectral theory of $h_{\nu}$ , including the Hankel transform to diagonalise it, is classical; see, for example, [11]. The new point here is different. We compute $A_{\nu}^{\pm}$ as an explicit operator matrix in a basis adapted to the Hankel transform. This matrix separates the threshold vector, the low energy part, and the high energy part. It gives lower bounds, the large order asymptotic law, and the sharp imaginary part discovered by Simon in the Laplacian setting [31].

1.2 Multidimensional examples beyond the Laplacian

In this picture, the relevant parameter $\nu$ becomes the smallest effective Bessel order, which depends on the particular model considered.

The canonical example from (1.3) is the Laplacian $-\Delta$ in $L^{2}(\mathbb{R}^{d})$ with $d\geq 3$ defined on $H^{2}(\mathbb{R}^{d})$ . In this case, the orders are

\nu_{\ell,d}:=\ell+\frac{d}{2}-1,\qquad\ell\in\mathbb{N}:=\{0,1,2,\dots\}.

Another example is the inverse square Hamiltonian

H_{d,c}:=-\Delta+c\,|x|^{-2}\qquad\mbox{in}\qquad L^{2}(\mathbb{R}^{d})

(1.6)

with $d\geq 2$ and real $c>-c_{d}=-(d-2)^{2}/4$ . We introduce it as the Friedrichs extension of this operator initially defined on $C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ . In this case, replacing $-\Delta$ by $H_{d,c}$ in (1.3), the orders become

\nu_{\ell,d,c}:=\sqrt{\nu_{\ell,d}^{2}+c}.

Finally, we consider the two-dimensional magnetic Hamiltonian

-\Delta_{\alpha}:=(-i\nabla-A_{\alpha})^{2}\qquad\mbox{in}\qquad L^{2}(\mathbb{R}^{2})

(1.7)

with the Aharonov–Bohm potential $A_{\alpha}(x):=\alpha\,(-x_{2},x_{1})\,|x|^{-2}$ , where $\alpha\in\mathbb{R}$ , introduced as the Friedrichs extension of this operator initially defined on $C_{0}^{\infty}(\mathbb{R}^{2}\setminus\{0\})$ . In this case,

\nu_{m,\alpha}:=|m-\alpha|,\qquad m\in\mathbb{Z}.

The applications below are consequences of this single channel principle.

1.3 Main results for the one-dimensional Bessel operators

The main methodological novelty of this paper is Theorem 2.1 below, in which we reveal a suitable matrix representation of $A_{\nu}^{\pm}$ introduced in (1.5). It would be too cumbersome to present it in the introduction. Let us only mention that it is based on an orthonormal basis adapted to the Hankel transform. Instead, we collect here the main results about the norm of $A_{\nu}^{\pm}$ , which we establish in this paper based particularly on the crucial Theorem 2.1.

Theorem 1.1 (Bessel operators).

Given any positive $\nu$ , let $A_{\nu}^{\pm}$ be as in (1.5). Then

(i)

$\displaystyle\|\mathrm{Im}\,A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}=\frac{\pi}{4\nu}$ ,
(ii)

$\displaystyle\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}=\frac{1}{\nu^{2}}$ if $0<\nu\leq\frac{1}{2}$ ,
(iii)

$\displaystyle\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}=\frac{1}{\nu}+O(\nu^{-3/2})$ as $\nu\to\infty$ ,

(iv)

$\displaystyle\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}\leq\mathcal{C}(\nu)$ , where

\mathcal{C}(\nu):=\begin{cases}1/\nu^{2},&0<\nu\leq\frac{1}{2},\\ 2/\nu,&\frac{1}{2}<\nu\leq\frac{1}{\sqrt{3}},\\ 1/\sqrt{\nu^{2}-\frac{1}{4}},&\frac{1}{\sqrt{3}}<\nu,\end{cases}

(1.8)

Part (i) is an analogue of the Simon-type identity [31] for the Laplacian, see (1.11) below. Part (ii) demonstrates an interesting fact that an exact identity can be derived also for the full complex value of $A_{\nu}^{\pm}$ , provided that the argument $\nu$ is sufficiently small. For all the other values of $\nu$ , we still have explicit lower bounds (see Theorem 2.4 and Proposition 3.10). For large values of $\nu$ , part (iii) gives a sharp lower bound. Part (iv) is the main result in the spirit of this paper to derive quantitative uniform resolvent estimates.

1.4 Main results for the multidimensional examples

We next state the consequences for the full weighted resolvent. For the free Laplacian, the Bessel channel estimates imply the following explicit bound.

Theorem 1.2 (Laplacian).

For every $d\geq 3$ ,

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|\leq\begin{cases}4,&d=3,\\ \displaystyle\frac{2}{\sqrt{(d-1)(d-3)}},&d\geq 4.\end{cases}

(1.9)

As a byproduct, we obtain the following monotonicity property.

Proposition 1.3 (Dimensional monotonicity).

For each $z\in\mathbb{C}\setminus[0,\infty)$ , the maps

d\mapsto\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{2d+1})\to L^{2}(\mathbb{R}^{2d+1})}

and

d\mapsto\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{2d+2})\to L^{2}(\mathbb{R}^{2d+2})}

are nonincreasing on $\{1,2,\dots\}$ .

We emphasise that our approach provides not only an explicit constant in (1.9), but also a new proof of the very fact (1.2). Actually, we give two arguments, the alternative being the method of multipliers, which has been used by other authors, too.

Except for the well-known three-dimensional case, we are not aware of attempts to quantify the uniform resolvent estimate (1.2). (An exception is the crude estimate obtained by Cossetti and one of the present authors in [9].) Be all as it may be, we note that the same bound (1.9) can be deciphered from the proof of [5, Thm. 3] (after quantifying the generic constants appearing there), which is close to our multiplier approach (see also [4, Thm. 1]). We also note that a similar bound with a worse constant appears in [25, Thm. 4] (see also [26, Thm. 1]).

The constant of Theorem 1.2 is sharp in three dimensions:

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{3})\to L^{2}(\mathbb{R}^{3})}=4.

(1.10)

This follows from an explicit lower bound that we also derive (Proposition 4.2). We leave as an open problem whether our estimate (1.9) is sharp for $d\geq 4$ , but we prove that it is sharp in the limit of high dimensions (Proposition 4.5):

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}=\frac{2}{d}+O(d^{-3/2})\qquad\mbox{as}\qquad d\to\infty.

Sharp results can be obtained also for the imaginary part of the operators. In particular, we recover Simon’s identity [31]

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}\mathrm{Im}\,(-\Delta-z)^{-1}|x|^{-1}\right\|=\frac{\pi}{2(d-2)}.

(1.11)

In our approach, Simon’s identity is not isolated but just a special case of analogous exact identities for the other critical models. Indeed, (1.11) is the first member of a family governed only by the lowest effective angular momentum, as the following results indicate.

Theorem 1.4 (Exact imaginary parts).

Given $d\geq 2$ and $c>-c_{d}$ , let $H_{d,c}$ be the inverse square Hamiltonian defined in (1.6). Then

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}\mathrm{Im}\,(H_{d,c}-z)^{-1}|x|^{-1}\right\|=\frac{\pi}{4\nu_{0,d,c}}.

(1.12)

Moreover, for the two-dimensional Aharonov–Bohm Hamiltonian $-\Delta_{\alpha}$ defined in (1.7), if $\alpha\notin\mathbb{Z}$ , then

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}\mathrm{Im}\,(-\Delta_{\alpha}-z)^{-1}|x|^{-1}\right\|=\frac{\pi}{4\,\mathrm{dist}(\alpha,\mathbb{Z})}.

(1.13)

Coming back to the uniform resolvent estimates, for the inverse square family we obtain the following result.

Theorem 1.5 (Inverse square potentials).

Given $d\geq 2$ and $c>-c_{d}$ , let $H_{d,c}$ be as in (1.6). Then

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(H_{d,c}-z)^{-1}|x|^{-1}\right\|\leq\mathcal{C}(\nu_{0,d,c}),

(1.14)

where $\mathcal{C}(\nu)$ is defined in (1.8). If $0<\nu_{0,d,c}\leq 1/2$ , then equality holds in (1.14).

Finally, the same channel analysis yields the Aharonov–Bohm estimate in two dimensions. Here the lowest effective order is the distance from the magnetic flux to the integers.

Theorem 1.6 (Aharonov–Bohm potentials).

Given $\alpha\notin\mathbb{Z}$ , let $-\Delta_{\alpha}$ be as in (1.7). Then

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta_{\alpha}-z)^{-1}|x|^{-1}\right\|=\frac{1}{\mathrm{dist}(\alpha,\mathbb{Z})^{2}}.

(1.15)

For $\alpha\in\mathbb{Z}$ , $-\Delta_{\alpha}$ is unitarily equivalent to $-\Delta$ in $L^{2}(\mathbb{R}^{2})$ and there is no uniform resolvent estimate for the Laplacian in two dimensions.

1.5 Applications

Kato smoothing and evolution control

Let $H=-\Delta$ or $H_{d,c}$ or $-\Delta_{\alpha}$ . The exact imaginary part identities (1.11), (1.12) and (1.13) give sharp smoothness constants in Kato’s theorem [19, 10]. Indeed, denoting

B_{H}:=\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}\mathrm{Im}\,(H-z)^{-1}|x|^{-1}\right\|

we have, for every $u\in L^{2}(\mathbb{R}^{d})$ ,

\int_{\mathbb{R}}\left\||x|^{-1}e^{-itH}u\right\|_{L^{2}(\mathbb{R}^{d})}^{2}\,\mathrm{d}t\leq 2B_{H}\|u\|_{L^{2}(\mathbb{R}^{d})}^{2}.

In the cases covered by (1.11) and Theorem 1.4, this gives the explicit constants

2B_{-\Delta}=\frac{\pi}{d-2},\qquad 2B_{H_{d,c}}=\frac{\pi}{2\nu_{0,d,c}},\qquad 2B_{-\Delta_{\alpha}}=\frac{\pi}{2\,\mathrm{dist}(\alpha,\mathbb{Z})}.

The full weighted resolvent bounds of Theorems 1.2, 1.5 and 1.6 give the corresponding supersmooth estimates, which are equivalent to the inhomogeneous smoothing estimate below [10]. Namely, with

C_{H}:=\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(H-z)^{-1}|x|^{-1}\right\|,

(1.16)

Kato’s retarded estimate yields

\left\||x|^{-1}\int_{0}^{t}e^{-i(t-s)H}F(s)\,\mathrm{d}s\right\|_{L^{2}(\mathbb{R};L^{2}(\mathbb{R}^{d})))}\leq 2\pi C_{H}\||x|F\|_{L^{2}(\mathbb{R};L^{2}(\mathbb{R}^{d})))}.

Here $C_{H}$ is bounded by the constants in Theorems 1.2, 1.5, and 1.6; in the Aharonov–Bohm case (Theorem 1.6), $C_{H}$ is exactly $\mathrm{dist}(\alpha,\mathbb{Z})^{-2}$ .

Spectral enclosures and stability

Let again $H$ be one of the self-adjoint model operators above. Consider a possibly complex-valued perturbation $V:\mathbb{R}^{d}\to\mathbb{C}$ such that $|x|^{2}\,V\in L^{\infty}(\mathbb{R}^{d})$ . The Birman–Schwinger principle [18] gives the following necessary condition for an eigenvalue $z\in\mathbb{C}\setminus[0,\infty)$ of a closed realisation $H_{V}$ of $H+V$ to exist:

1\leq\||x|^{2}V\|_{L^{\infty}(\mathbb{R}^{d})}\left\||x|^{-1}(H-z)^{-1}|x|^{-1}\right\|.

Equivalently,

\sigma_{\rm p}(H_{V})\setminus[0,\infty)\subset\left\{z\in\mathbb{C}\setminus[0,\infty):\left\||x|^{-1}(H-z)^{-1}|x|^{-1}\right\|\geq\||x|^{2}V\|_{L^{\infty}(\mathbb{R}^{d})}^{-1}\right\}.

In particular, a sufficient condition for the absence of eigenvalues of $H$ in $\mathbb{C}\setminus[0,\infty)$ reads

C_{H}\,\||x|^{2}V\|_{L^{\infty}(\mathbb{R}^{d})}<1,

(1.17)

where $C_{H}$ is defined in (1.16). What is more, the same smallness condition (1.17) ensures the total absence of eigenvalues (including eigenvalues in $[0,\infty)$ ) [18, Corol. 2]. In fact, $H$ is similar to $-\Delta$ via bounded and boundedly invertible transformation [19].

Theorems 1.2, 1.5, and 1.6 make the sufficient condition (1.17) quantitative. For instance, in the Aharonov–Bohm case this smallness condition reads $\||x|^{2}V\|_{L^{\infty}(\mathbb{R}^{d})}<\mathrm{dist}(\alpha,\mathbb{Z})^{2}$ .

1.6 Structure of the paper

The paper is organised as follows. In Section 2, we begin with the matrix representation $A_{\nu}^{\pm}$ and its immediate consequences; in particular, we establish parts (i) and (iii) of Theorem 1.1. More subtle one-channel estimates, including parts (ii) and (iv) of Theorem 1.1, are established in Section 3. The applications to the Laplacian, inverse square potentials and the Aharonov–Bohm magnetic fields are performed in Section 4. In the final Section 5, we present an alternative way to produce the uniform resolvent estimate for the Laplacian (Theorem 1.2) via the method of multipliers.

2 The matrix computation for one Bessel channel

Given any $\nu>0$ , recall the definition (1.4) of the Bessel operator $h_{\nu}$ from the introduction. We are interested in the boundary values $A_{\nu}^{\pm}$ of the weighted resolvent at the positive energy $1$ introduced in (1.5). By a scaling invariance and the maximum principle (see Section 3.2 below), the norm of $r^{-1}(h_{\nu}-(\lambda\pm i0))^{-1}r^{-1}$ is independent of $\lambda>0$ and equals $\sup_{z\in\mathbb{C}\setminus[0,\infty)}\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\|_{L^{2}\to L^{2}}$ . So it is indeed enough to understand (1.5).

In this section, we decompose $A^{\pm}_{\nu}$ in a suitable orthonormal basis of $L^{2}$ , and then deduce from it both a lower bound and the asymptotic behaviour as $\nu\to\infty$ of $\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}$ . More specifically, the main result of this section is the following theorem.

Theorem 2.1.

Let $\nu>0$ . There is an orthonormal basis $(e_{\nu,\ell})_{\ell\in\mathbb{Z}}$ of $L^{2}$ such that $A^{\pm}_{\nu}$ decomposes in this basis as follows: Denoting $a_{\ell\ell^{\prime}}:=\langle e_{\nu,\ell},A^{\pm}_{\nu}e_{\nu,\ell^{\prime}}\rangle_{L^{2}}$ ,

\displaystyle(a_{\ell\ell^{\prime}})_{\ell,\ell^{\prime}\in\mathbb{Z}}=\begin{pmatrix}\ddots&\vdots&\vdots&\vdots&\vdots&\vdots&\reflectbox{$\ddots$}\\ \dots&a_{-2,-2}&0&a_{-2,0}&0&0&\dots\\ \dots&0&a_{-1,-1}&a_{-1,0}&0&0&\dots\\ \dots&a_{0,-2}&a_{0,-1}&a_{0,0}&a_{0,1}&a_{0,2}&\dots\\ \dots&0&0&a_{1,0}&a_{1,1}&a_{1,2}&\dots\\ \dots&0&0&a_{2,0}&a_{2,1}&a_{2,2}&\dots\\ \reflectbox{$\ddots$}&\vdots&\vdots&\vdots&\vdots&\vdots&\ddots\end{pmatrix},

where for $\ell,\ell^{\prime}\geq 1$ ,

	$\displaystyle a_{0,0}=\pm i\frac{\pi}{4\nu}+\frac{1}{4\nu^{2}},$
	$\displaystyle a_{-\ell,\ell^{\prime}}=a_{\ell^{\prime},-\ell}=0,$
	$\displaystyle a_{-\ell,-\ell^{\prime}}=-\frac{\delta_{\ell\ell^{\prime}}}{4\ell(\nu+\ell)},$
	$\displaystyle a_{0,-\ell}=a_{-\ell,0}=\frac{(-1)^{\ell}}{4\ell(\nu+\ell)}\sqrt{1+\frac{2\ell}{\nu}},$
	$\displaystyle a_{0,\ell}=a_{\ell,0}=-\frac{1}{4\ell(\nu+\ell)}\left[(-1)^{\ell}-\frac{\ell!}{(\nu)_{\ell}}\right]\sqrt{1+\frac{2\ell}{\nu}},$
	$\displaystyle a_{\ell,\ell^{\prime}}=\frac{1}{4\ell\ell^{\prime}}\frac{\max(\ell,\ell^{\prime})!\,\Gamma(\nu+\min(\ell,\ell^{\prime}))}{(\min(\ell,\ell^{\prime})-1)!\,\Gamma(\nu+\max(\ell,\ell^{\prime})+1)}\sqrt{1+\frac{2\ell}{\nu}}\sqrt{1+\frac{2\ell^{\prime}}{\nu}},$

where $(\nu)_{\ell}=\nu(\nu+1)\cdots(\nu+\ell-1)$ is the Pochammer symbol.

Theorem 2.1 will be proven below. Before, we discuss its consequences.

2.1 A lower bound to the norm of $A_{\nu}^{\pm}$

In what follows, $J_{\nu}$ , $Y_{\nu}$ are the usual Bessel functions of the first and second kind, respectively. Moreover, $H_{\nu}^{(1)}=J_{\nu}+iY_{\nu}$ and $H_{\nu}^{(2)}=J_{\nu}-iY_{\nu}$ are the Hankel function of the first and second kind, respectively. To shorten the notations, we denote $p_{\nu}(r):=r^{-1/2}J_{\nu}(r),$ and define the distinguished vector

e_{\nu,0}:=(2\nu)^{1/2}\,p_{\nu}=(2\nu)^{1/2}r^{-1/2}J_{\nu}(r)\,.

(2.1)

First, we are concerned with the imaginary part of $A_{\nu}^{\pm}$ , establishing part (i) of Theorem 1.1.

Proposition 2.2.

Let $\nu>0$ . One has

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}\mathrm{Im}\,(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}=\frac{\pi}{4\nu}.

(2.2)

Proof.

By the maximum principle and scaling (see Section 3.2 below) it is enough to compute the norm of the boundary values at $1\pm i0$ . The integral kernel of $h_{\nu}$ at $1\pm i0$ reads

	$\displaystyle(h_{\nu}-(1+i0))^{-1}(r,R)$	$\displaystyle=\frac{i\pi}{2}(rR)^{1/2}$		(2.3)
	$\displaystyle(h_{\nu}-(1-i0))^{-1}(r,R)$	$\displaystyle=-\frac{i\pi}{2}(rR)^{1/2}$		(2.3)

Since $J_{\nu}$ , $Y_{\nu}$ are real-valued, one has

\mathrm{Im}\,A_{\nu}^{\pm}(r,R)=\pm\frac{\pi}{2}(rR)^{-1/2}J_{\nu}(r)J_{\nu}(R).

Thus the imaginary part is the rank-one operator

\pm\mathrm{Im}\,A_{\nu}^{\pm}=\frac{\pi}{2}|r^{-1/2}J_{\nu}\rangle\langle r^{-1/2}J_{\nu}|=\frac{\pi}{2}|p_{\nu}\rangle\langle p_{\nu}|=\frac{\pi}{4\nu}|e_{\nu,0}\rangle\langle e_{\nu,0}|,

(2.4)

where we have used $\|p_{\nu}\|_{L^{2}}^{2}=\|r^{-1/2}J_{\nu}\|_{L^{2}}^{2}=(2\nu)^{-1}$ . This proves (2.2). ∎

As for the real part of $A_{\nu}^{\pm}$ , we have the following result.

Proposition 2.3.

Let $\nu>0$ . One has

\left\langle e_{\nu,0},\mathrm{Re}\,A_{\nu}^{\pm}e_{\nu,0}\right\rangle_{L^{2}}=\frac{1}{4\nu^{2}}.

(2.5)

Proof.

From (2.3) we deduce

\mathrm{Re}\,A_{\nu}^{\pm}(r,R)=-\frac{\pi}{2}(rR)^{-1/2}\begin{cases}J_{\nu}(r)Y_{\nu}(R)&\mbox{if}\quad r<R,\\ J_{\nu}(R)Y_{\nu}(r)&\mbox{if}\quad R<r.\end{cases}

(2.6)

Hence

\displaystyle\left\langle p_{\nu},A_{\nu}^{\pm}p_{\nu}\right\rangle_{L^{2}}=-\pi\int_{0}^{\infty}\frac{J_{\nu}(r)^{2}}{r}\int_{r}^{\infty}\frac{Y_{\nu}(R)J_{\nu}(R)}{R}\,\mathrm{d}R\,\mathrm{d}r.

We compute the last integral by the Mellin transform. First,

\displaystyle\mathcal{M}\{r^{-1}J_{\nu}(r)^{2}\}(s)=\int_{0}^{\infty}r^{s-2}J_{\nu}(r)^{2}\,\mathrm{d}r2^{s-2}\frac{\Gamma(2-s)\Gamma\big(\nu-\frac{1}{2}+\frac{s}{2}\big)}{\Gamma\big(\frac{3}{2}-\frac{s}{2}\big)^{2}\Gamma\big(\nu+\frac{3}{2}-\frac{s}{2}\big)}

for $1-2\nu<\mathrm{Re}(s)<2$ , see [17, 6.574.2]. Likewise

	$\displaystyle\mathcal{M}\{r^{-1}J_{\nu}(r)Y_{\nu}(r)\}(s)$	$\displaystyle=\int_{0}^{\infty}r^{s-2}J_{\nu}(r)Y_{\nu}(r)\,\mathrm{d}r$
		$\displaystyle=-\frac{2^{s-1}}{\pi}\cos\left(\frac{(s-1)\pi}{2}\right)\frac{\Gamma(\nu-\frac{1}{2}+\frac{s}{2})\Gamma(-\frac{1}{2}+\frac{s}{2})\Gamma(2-s)}{\Gamma(\nu+\frac{3}{2}-\frac{s}{2})\Gamma(\frac{3}{2}-\frac{s}{2})}$

for $0<\mathrm{Re}\,s<1$ , see [28, 2.13.15.5]. Therefore

\displaystyle\mathcal{M}\left\{\int_{r}^{\infty}\frac{Y_{\nu}(R)J_{\nu}(R)}{R}\,\mathrm{d}R\right\}(s)=-\frac{2^{s}}{\pi s}\cos\left(\frac{s\pi}{2}\right)\frac{\Gamma(\nu+\frac{s}{2})\Gamma(\frac{s}{2})\Gamma(1-s)}{\Gamma(\nu+1-\frac{s}{2})\Gamma(1-\frac{s}{2})}

for $-1<\mathrm{Re}\,s<0$ . Parseval’s formula gives, for $0<c<1$ ,

\displaystyle-\pi\int_{0}^{\infty}\frac{J_{\nu}(r)^{2}}{r}\int_{r}^{\infty}\frac{Y_{\nu}(R)J_{\nu}(R)}{R}\,\mathrm{d}R\,\mathrm{d}r=\frac{1}{2i\pi}\int_{c-i\infty}^{c+i\infty}\frac{1}{(1-s)^{2}}\frac{1}{\big(\nu+\frac{1}{2}-\frac{s}{2}\big)\big(\nu-\frac{1}{2}+\frac{s}{2}\big)}\,\mathrm{d}s.

The residue at $s=-2\nu+1$ equals $1/(8\nu^{3})$ . This proves the desired identity after recalling the relationship (2.1). ∎

The imaginary part identity from Proposition 2.2 gives a lower bound for $\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}$ . Indeed, for all $\varphi\in L^{2}$ , one has

\displaystyle\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}\geq\big|\langle\varphi,A_{\nu}^{\pm}\varphi\rangle_{L^{2}}\big|

\displaystyle\geq\big|\mathrm{Im}\,\langle\varphi,A_{\nu}^{\pm}\varphi\rangle_{L^{2}}\big|=\big|\langle\varphi,\big[\mathrm{Im}\,A_{\nu}^{\pm}\big]\varphi\rangle_{L^{2}}\big|.

Taking the supremum over $\varphi\in L^{2}$ such that $\|\varphi\|_{L^{2}}=1$ and using Proposition 2.2 gives

\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}\geq\frac{\pi}{4\nu}.

Now we are in a position to improve this bound. Indeed, from the matrix representation of Theorem 2.1 one can compute $\|A_{\nu}^{\pm}e_{\nu,0}\|_{L^{2}}$ , as the following theorem shows.

Theorem 2.4.

Let $\nu>0$ . One has

\displaystyle\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}^{2}\geq\big\|A_{\nu}^{\pm}e_{\nu,0}\big\|_{L^{2}}^{2}=\frac{\pi^{2}}{12\nu^{2}}+\frac{1}{8\nu^{4}}.

Proof.

Since $e_{\nu,0}(r)=(2\nu)^{\frac{1}{2}}r^{-\frac{1}{2}}J_{\nu}(r)$ is real-valued, we can write $A_{\nu}^{\pm}e_{\nu,0}=(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}+i(\mathrm{Im}\,A_{\nu}^{\pm})e_{\nu,0}$ . This implies

\displaystyle\big\|A_{\nu}^{\pm}\big\|_{L^{2}\to L^{2}}^{2}\geq\big\|A_{\nu}^{\pm}e_{\nu,0}\big\|_{L^{2}}^{2}=\big\|(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}\big\|_{L^{2}}^{2}+\big\|(\mathrm{Im}\,A_{\nu}^{\pm})e_{\nu,0}\big\|_{L^{2}}^{2}=\big\|(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}\big\|_{L^{2}}^{2}+\frac{\pi^{2}}{16\nu^{2}},

where we have used (2.4) in the last equality. It remains to compute $\big\|(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}\big\|_{L^{2}}^{2}$ . From the matrix representation of Theorem 2.1 and Parseval’s identity, it follows that

$\displaystyle\big\\|(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}\big\\|_{L^{2}}^{2}$	$\displaystyle=\frac{1}{16\nu^{4}}+\frac{1}{16}\sum_{\ell\geq 1}\frac{1}{\ell^{2}(\nu+\ell)^{2}}\Big(1+\frac{2\ell}{\nu}\Big)+\frac{1}{16}\sum_{\ell\geq 1}\frac{1}{\ell^{2}(\nu+\ell)^{2}}\Big((-1)^{\ell}-\frac{l!}{(\nu)_{\ell}}\Big)^{2}\Big(1+\frac{2\ell}{\nu}\Big)$
	$\displaystyle=\frac{1}{16\nu^{4}}+\frac{1}{8\nu^{2}}\sum_{\ell\geq 1}\Big(\frac{1}{\ell^{2}}-\frac{1}{(\nu+\ell)^{2}}\Big)-\frac{1}{8\nu^{2}}\sum_{\ell\geq 1}(-1)^{\ell}\Big(\frac{1}{\ell^{2}}-\frac{1}{(\nu+\ell)^{2}}\Big)\frac{\ell!}{(\nu)_{\ell}}$
	$\displaystyle\quad+\frac{1}{16\nu^{2}}\sum_{\ell\geq 1}\Big(\frac{1}{\ell^{2}}-\frac{1}{(\nu+\ell)^{2}}\Big)\Big(\frac{\ell!}{(\nu)_{\ell}}\Big)^{2}.$	(2.7)

We compute the three sums on the right-hand side of (2.7). For the first sum, by definition of the trigamma function, we have

\displaystyle\sum_{\ell\geq 1}\Big(\frac{1}{\ell^{2}}-\frac{1}{(\nu+\ell)^{2}}\Big)=\frac{\pi^{2}}{6}-\psi^{(1)}(\nu)+\frac{1}{\nu^{2}}.

(2.8)

To compute the second sum, we define

\displaystyle S_{\nu}:=\sum_{\ell\geq 1}(-1)^{\ell}\Big(\frac{1}{\ell^{2}}-\frac{1}{(\nu+\ell)^{2}}\Big)\frac{\ell!}{(\nu)_{\ell}},

and then compute

\displaystyle S_{\nu+1}-S_{\nu}

\displaystyle=\sum_{\ell=1}^{\infty}(-1)^{\ell}\frac{\ell!}{(\nu)_{\ell}}\Big[-\frac{1}{\ell(\nu+\ell)}+\frac{\ell+1}{(\nu+\ell)(\nu+\ell+1)^{2}}+\frac{1}{(\nu+\ell)^{2}(\nu+\ell+1)}\Big]=\frac{1}{(\nu+1)^{2}},

which yields

\displaystyle S_{\nu}=\sum_{k\geq 0}(S_{\nu+k}-S_{\nu+k+1})=-\sum_{k\geq 0}\frac{1}{(\nu+k+1)^{2}}=\frac{1}{\nu^{2}}-\psi^{(1)}(\nu).

(2.9)

The third sum in (2.7) simplifies to

\displaystyle\sum_{\ell\geq 1}\Big(\frac{1}{\ell^{2}}-\frac{1}{(\nu+\ell)^{2}}\Big)\Big(\frac{\ell!}{(\nu)_{\ell}}\Big)^{2}=\sum_{\ell\geq 1}\Big(\Big(\frac{(\ell-1)!}{(\nu)_{\ell}}\Big)^{2}-\Big(\frac{\ell!}{(\nu)_{\ell+1}}\Big)^{2}\Big)=\frac{1}{\nu^{2}}.

(2.10)

Inserting (2.8)–(2.10) into (2.7), the result follows. ∎

2.2 Asymptotics of the norm of $A_{\nu}^{\pm}$ as $\nu\to\infty$

It is not clear to us how to compute an exact expression of $\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}$ from the matrix representation of Theorem 2.1. However, we can obtain the following asymptotics as $\nu\to\infty$ , establishing part (iii) of Theorem 1.1.

Theorem 2.5.

Let $\nu>0$ . One has

\displaystyle\big\|A_{\nu}^{\pm}\big\|_{L^{2}\to L^{2}}=\frac{1}{\nu}+O(\nu^{-3/2})\qquad\mbox{as}\qquad\quad\nu\to\infty.

Remark 2.6.

We will show more precisely that, for all $\nu>0$ ,

\displaystyle\Big|\big\|A_{\nu}^{\pm}\big\|_{L^{2}\to L^{2}}-\frac{1}{\nu}\Big|\leq\frac{\sqrt{3\psi^{(1)}(\nu+1)}}{2\sqrt{2}\nu}+\frac{1}{\nu^{2}},

(2.11)

where $\psi^{(1)}$ is the trigamma function (in particular $\psi^{(1)}(\nu)\sim\frac{1}{\nu}$ as $\nu\to\infty$ ), without trying to optimise the remainder term.

We prove Theorem 2.5 thanks to a suitable matrix representation of the auxiliary operator

B_{\nu}^{\pm}:=A_{\nu}^{\pm}-r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}.

(2.12)

Remark 2.7.

The operator $r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}$ is self-adjoint and its spectrum reads

\displaystyle\sigma\big(r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}\big)=\sigma_{\mathrm{ess}}\big(r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}\big)=[0,\nu^{-2}].

(2.13)

Consequently,

\|r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}\|_{L^{2}\to L^{2}}=\nu^{-2}.

(2.14)

This follows from the integral kernel

h_{\nu}^{-1}(r,R)=\frac{1}{2\nu}\begin{cases}r^{\nu+\frac{1}{2}}R^{-\nu+\frac{1}{2}}&\mbox{if}\quad\quad r<R\,,\\ r^{-\nu+\frac{1}{2}}R^{\nu+\frac{1}{2}}&\mbox{if}\quad\quad R<r\,.\end{cases}

In more detail, let $\mathcal{D}:L^{2}((0,\infty))\to L^{2}(\mathbb{R})$ be the unitary operator $(\mathcal{D}\varphi)(x)=e^{x/2}\varphi(e^{x}).$ Under this transform, $r^{-1}h_{\nu}^{-1}r^{-1}$ is unitarily equivalent to the operator with kernel $(2\nu)^{-1}e^{-\nu|x-y|}$ . After Fourier transform, this becomes the multiplier $(\xi^{2}+\nu^{2})^{-1}$ . This implies that

\displaystyle\sigma\big(r^{-1}h_{\nu}^{-1}r^{-1}\big)=\sigma_{\mathrm{ess}}\big(r^{-1}h_{\nu}^{-1}r^{-1}\big)=[0,\nu^{-2}],

and that $\|r^{-1}h_{\nu}^{-1}r^{-1}\|=\nu^{-2}$ . By scaling, the cutoff operator with $\mathbf{1}_{h_{\nu}\geq 1}$ is unitarily equivalent to the cutoff at $\varepsilon>0$ . Letting $\varepsilon\to 0$ and using strong convergence, we arrive at (2.13) by a standard spectral convergence results (see, e.g., [29, Thm. VIII.24]).

Now we are in a position to present the matrix representation of $B_{\nu}^{\pm}$ .

Proposition 2.8.

Let $\nu>0$ . In the orthonormal basis $(e_{\nu,\ell})_{\ell\in\mathbb{Z}}$ of $L^{2}$ from Theorem 2.1, the operator $B^{\pm}_{\nu}$ decomposes as follows: Denoting $b_{\ell\ell^{\prime}}:=\langle e_{\nu,\ell},B_{\nu}^{\pm}e_{\nu,\ell^{\prime}}\rangle$ ,

\displaystyle(b_{\ell\ell^{\prime}})_{\ell,\ell^{\prime}\in\mathbb{Z}}=\begin{pmatrix}\ddots&\vdots&\vdots&\vdots&\vdots&\vdots&\reflectbox{$\ddots$}\\ \dots&b_{-2,-2}&0&b_{-2,0}&0&0&\dots\\ \dots&0&b_{-1,-1}&b_{-1,0}&0&0&\dots\\ \dots&b_{0,-2}&b_{0,-1}&b_{0,0}&b_{0,1}&b_{0,2}&\dots\\ \dots&0&0&b_{1,0}&b_{1,1}&0&\dots\\ \dots&0&0&b_{2,0}&0&b_{2,2}&\dots\\ \reflectbox{$\ddots$}&\vdots&\vdots&\vdots&\vdots&\vdots&\ddots\end{pmatrix},

where for $\ell,\ell^{\prime}\geq 1$ ,

	$\displaystyle b_{0,0}=\pm i\frac{\pi}{4\nu},$
	$\displaystyle b_{-\ell,\ell^{\prime}}=b_{\ell^{\prime},-\ell}=0,$
	$\displaystyle b_{-\ell,-\ell^{\prime}}=-\frac{\delta_{\ell\ell^{\prime}}}{4\ell(\nu+\ell)},$
	$\displaystyle b_{0,-\ell}=b_{-\ell,0}=\frac{(-1)^{\ell}}{4\ell(\nu+\ell)}\sqrt{1+\frac{2\ell}{\nu}},$
	$\displaystyle b_{0,\ell}=b_{\ell,0}=\frac{(-1)^{\ell+1}}{4\ell(\nu+\ell)}\sqrt{1+\frac{2\ell}{\nu}},$
	$\displaystyle b_{\ell,\ell^{\prime}}=\frac{\delta_{\ell\ell^{\prime}}}{4\ell(\nu+\ell)}.$

Remark 2.9.

Proposition 2.8 implies by direct computation that, for all $\nu>0$ , the operator $B_{\nu}^{\pm}$ is Hilbert–Schmidt with norm

	$\displaystyle\big\\|B_{\nu}^{\pm}\big\\|_{\mathrm{HS}}^{2}$	$\displaystyle=\sum_{\ell\in\mathbb{Z}}\big\\|B_{\nu}^{\pm}e_{\nu,\ell}\big\\|_{L^{2}}^{2}$
		$\displaystyle=\frac{\pi^{2}}{16\nu^{2}}+\sum_{\ell\geq 1}\frac{1+\frac{2\ell}{\nu}}{8\ell^{2}(\nu+\ell)^{2}}+\sum_{\ell\geq 1}\frac{1+\frac{\ell}{\nu}}{4\ell^{2}(\nu+\ell)^{2}}=\frac{\pi^{2}}{8\nu^{2}}-\frac{\psi^{(1)}(\nu+1)}{8\nu^{2}}-\frac{\psi(\nu+1)+\gamma}{4\nu^{3}},$

where as before, $\psi^{(1)}$ is the trigamma function, $\psi$ is the digamma function and $\gamma$ is Euler’s constant. Combined with Remark 2.7, this in turn implies that

\displaystyle\sigma_{\mathrm{ess}}\big(A_{\nu}^{\pm}\big)=[0,\nu^{-2}],

for all $\nu>0$ .

Proposition 2.8 will be proven below, after the proof of Theorem 2.1. In order to deduce Theorem 2.5 from the former, we first identify a matrix $C^{\pm}$ , independent of $\nu$ , such that $\|B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\|_{L^{2}\to L^{2}}$ is negligible compared to $\nu^{-1}$ as $\nu\to\infty$ . Let $C^{\pm}$ be the operator in $L^{2}$ , whose matrix $(c_{\ell,\ell^{\prime}})_{\ell,\ell^{\prime}\in\mathbb{Z}}$ with respect to the orthonormal basis $(e_{\nu,\ell})_{\ell\in\mathbb{Z}}$ reads as follows:

	$\displaystyle c_{0,0}=\pm i\frac{\pi}{4},$
	$\displaystyle c_{-\ell,\ell^{\prime}}=c_{\ell^{\prime},-\ell}=0,$
	$\displaystyle c_{-\ell,-\ell^{\prime}}=-\frac{\delta_{\ell\ell^{\prime}}}{4\ell},$
	$\displaystyle c_{0,-\ell}=c_{-\ell,0}=\frac{(-1)^{\ell}}{4\ell},$
	$\displaystyle c_{0,\ell}=c_{\ell,0}=\frac{(-1)^{\ell+1}}{4\ell},$
	$\displaystyle c_{\ell,\ell^{\prime}}=\frac{\delta_{\ell\ell^{\prime}}}{4\ell}.$

Lemma 2.10.

Let $\nu>0$ . One has

\displaystyle\big\|B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\big\|^{2}_{\mathrm{HS}}\leq\frac{3\psi^{(1)}(\nu+1)}{8\nu^{2}}.

Proof.

A direct computation gives

	$\displaystyle\big\\|B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\big\\|_{\mathrm{HS}}^{2}=\big\\|(B_{\nu}^{\pm}-\nu^{-1}C^{\pm})e_{\nu,0}\big\\|_{L^{2}}^{2}+\sum_{\ell\in\mathbb{Z}\setminus\{0\}}\big\\|(B_{\nu}^{\pm}-\nu^{-1}C^{\pm})e_{\nu,\ell}\big\\|_{L^{2}}^{2}$
	$\displaystyle=2\sum_{\ell\geq 1}\Big(\frac{1}{4\ell(\nu+\ell)}-\frac{1}{4\nu\ell}\Big)^{2}+4\sum_{\ell\geq 1}\Big(\frac{(1+\frac{2\ell}{\nu})^{\frac{1}{2}}}{4\ell(\nu+\ell)}-\frac{1}{4\nu\ell}\Big)^{2}$
	$\displaystyle=\frac{1}{8\nu^{2}}\sum_{\ell\geq 1}\frac{1}{(\nu+\ell)^{2}}+\frac{1}{4\nu^{2}}\sum_{\ell\geq 1}\frac{1}{(\nu+\ell)^{2}}\frac{(\frac{\ell}{\nu})^{2}}{((1+\frac{2\ell}{\nu})^{\frac{1}{2}}+1+\frac{\ell}{\nu})^{2}}$
	$\displaystyle\leq\frac{3}{8\nu^{2}}\sum_{\ell\geq 1}\frac{1}{(\nu+\ell)^{2}},$

which proves the result. ∎

Lemma 2.11.

We have

\displaystyle\big\|C^{\pm}\big\|_{L^{2}\to L^{2}}=1.

Proof.

First, we observe that the definition of the operator $C^{\pm}$ reduces to

\displaystyle c_{0,0}=\pm i\frac{\pi}{4},\qquad c_{k,k^{\prime}}=\frac{\delta_{kk^{\prime}}}{4k},\qquad c_{0,k}=c_{k,0}=\frac{(-1)^{k+1}}{4k},

for $k,k^{\prime}\in\mathbb{Z}\setminus\{0\}$ . We then introduce the unitary operator

\displaystyle W:L^{2}((0,\infty))\ni e_{\nu,k}\mapsto\tilde{e}_{k}\in L^{2}((0,2\pi)),\qquad\text{where}\qquad\forall k\in\mathbb{Z},\quad\tilde{e}_{k}(t):=(-1)^{k}(2\pi)^{-\frac{1}{2}}e^{ikt}.

In the orthonormal basis $(\tilde{e}_{k})_{k\in\mathbb{Z}}$ of $L^{2}((0,2\pi))$ , the operator $\tilde{C}:=WCW^{*}\cong(\tilde{c}_{\ell,\ell^{\prime}})_{\ell,\ell^{\prime}\in\mathbb{Z}}$ decomposes as

\displaystyle\tilde{c}_{0,0}=\pm i\frac{\pi}{4},\qquad\tilde{c}_{k,k^{\prime}}=\frac{\delta_{kk^{\prime}}}{4k},\qquad\tilde{c}_{0,k}=\tilde{c}_{k,0}=-\frac{1}{4k}.

Moreover, its integral kernel is

	$\displaystyle\tilde{C}^{\pm}(t,s)$	$\displaystyle=\frac{1}{2\pi}\sum_{j,k\in\mathbb{Z}}\tilde{c}_{k,j}e^{-ijs}e^{ikt}$
		$\displaystyle=\frac{1}{2\pi}\Big(\tilde{c}_{0,0}-\frac{1}{4}\sum_{k\in\mathbb{Z}\setminus\{0\}}\frac{e^{-iks}}{k}-\frac{1}{4}\sum_{k\in\mathbb{Z}\setminus\{0\}}\frac{e^{ikt}}{k}+\frac{1}{4}\sum_{k\in\mathbb{Z}\setminus\{0\}}\frac{e^{ik(t-s)}}{k}\Big).$

Let us consider, for instance, the $+$ case. Then $\tilde{c}_{0,0}=i\frac{\pi}{4}$ . Using that $\sum_{k\in\mathbb{Z}\setminus\{0\}}\frac{e^{ikt}}{k}=i(\pi-t)$ for $t\in(0,2\pi)$ , we obtain

\displaystyle\tilde{C}^{+}(t,s)

\displaystyle=\frac{1}{2\pi}\Big(i\frac{\pi}{4}+\frac{i}{4}(\pi-s)-\frac{i}{4}(\pi-t)+\frac{i}{4}\big((\pi-t+s)\mathbf{1}_{t\geq s}-(\pi-s+t)\mathbf{1}_{t<s}\big)\Big)=\frac{i}{4}\mathbf{1}_{t\geq s}.

Hence we have shown that, for all $\varphi\in L^{2}((0,2\pi))$ , $\tilde{C}^{+}$ is given by

\displaystyle\forall t\in(0,2\pi),\qquad(\tilde{C}^{+}\varphi)(t)=\frac{i}{4}\int_{0}^{t}\varphi(s)\mathrm{d}s.

The norm of this Volterra operator is standard to compute. Let $\varphi\in L^{2}((0,2\pi))$ be a normalised eigenstate, associated to a real eigenvalue $\mu$ , of $(\tilde{C}^{+})^{*}\tilde{C}^{+}$ . We have

\displaystyle\forall t\in(0,2\pi),\qquad((\tilde{C}^{+})^{*}\tilde{C}^{+}\varphi)(t)=\frac{1}{16}\int_{t}^{2\pi}\int_{0}^{s}\varphi(s^{\prime})\,\mathrm{d}s^{\prime}\mathrm{d}s=\mu f(t).

Hence $\varphi$ is a solution to $-16\mu\varphi^{\prime\prime}=\varphi$ with the boundary conditions $\varphi^{\prime}(0)=0$ , $\varphi(2\pi)=0$ . This implies that the eigenvalues of $(\tilde{C}^{+})^{*}\tilde{C}^{+}$ are $(2n+1)^{-2}$ , with $n\in\mathbb{N}$ , and therefore $\|\tilde{C}^{+}\|_{L^{2}\to L^{2}}=1$ . Since $W$ is unitary, this also shows that $\|C^{+}\|_{L^{2}\to L^{2}}=1$ . One can argue in the same way for $C^{-}$ . ∎

Now we can conclude the proof of Theorem 2.5.

Proof of Theorem 2.5.

It suffices to write, using Lemma 2.11,

	$\displaystyle\Big\|\big\\|A_{\nu}^{\pm}\big\\|_{L^{2}\to L^{2}}-\frac{1}{\nu}\Big\|$	$\displaystyle=\Big\|\big\\|A_{\nu}^{\pm}\big\\|_{L^{2}\to L^{2}}-\big\\|\nu^{-1}C^{\pm}\big\\|_{L^{2}\to L^{2}}\Big\|$
		$\displaystyle\leq\big\\|A_{\nu}^{\pm}-B_{\nu}^{\pm}+B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\big\\|_{L^{2}\to L^{2}}$
		$\displaystyle\leq\big\\|r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}\big\\|_{L^{2}\to L^{2}}+\big\\|B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\big\\|_{\mathrm{HS}},$

and then apply (2.14) and Lemma 2.10. ∎

2.3 The Hankel transform

In the remainder of this section, we establish Theorem 2.1 and Proposition 2.8. The main ingredient in the proof is the usual Hankel transform of order $\nu$ :

\displaystyle\mathcal{H}_{\nu}[\varphi](k):=\int_{0}^{\infty}\sqrt{kr}J_{\nu}(kr)\varphi(r)\,\mathrm{d}r,

which extends to a unitary operator on $L^{2}$ such that $\mathcal{H}_{\nu}^{-1}=\mathcal{H}_{\nu}$ . It can further be extended to a suitable space of distributions (similarly as the Fourier transform).

It is well known that the Hankel transform of order $\nu$ diagonalises the Bessel operator $h_{\nu}$ , in the sense that for any bounded measurable function $f:\mathbb{R}\to\mathbb{R}$ and $\varphi,\psi\in L^{2}$ ,

\displaystyle\langle\varphi,f(h_{\nu})\psi\rangle_{L^{2}}=\langle\mathcal{H}_{\nu}[\varphi],f(k^{2})\mathcal{H}_{\nu}[\psi]\rangle_{L^{2}}.

(2.15)

In particular, for all $\varepsilon>0$ and $\varphi,\psi\in L^{2}$ such that $r^{-1}\varphi,r^{-1}\psi\in L^{2}$ ,

\displaystyle\big\langle\varphi,r^{-1}(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}\psi\big\rangle_{L^{2}}=\int_{0}^{\infty}\frac{1}{k^{2}-(1\pm i\varepsilon)}\,\overline{\mathcal{H}_{\nu}[r^{-1}\varphi](k)}\,\mathcal{H}_{\nu}[r^{-1}\psi](k)\,\mathrm{d}k.

(2.16)

Let

\displaystyle\hat{\mathcal{Q}}_{\nu}:=\big\{\hat{\varphi}\in L^{2}_{\mathrm{loc}}((0,\infty))\,:\,\hat{\varphi}^{\prime}\in L^{2}((0,\infty))\big\}=\dot{H}^{1}_{0}((0,\infty)),

equipped with the scalar product defined by

\displaystyle\langle\hat{\varphi},\hat{\psi}\rangle_{\hat{\mathcal{Q}}_{\nu}}:=\langle\hat{\varphi},h_{\nu}\hat{\psi}\rangle_{L^{2}},

for any $\hat{\varphi},\hat{\psi}\in\hat{\mathcal{Q}}_{\nu}$ . We introduce the unitary operator

\displaystyle\mathcal{U}_{\nu}:L^{2}\ni\psi\mapsto\mathcal{H}_{\nu}[r^{-1}\psi]\in\hat{\mathcal{Q}}_{\nu},

whose inverse is given by $\mathcal{U}_{\nu}^{-1}\hat{\psi}=r\mathcal{H}_{\nu}[\hat{\psi}]$ , for any $\hat{\psi}\in\hat{\mathcal{Q}}_{\nu}$ . The facts that $\mathcal{U}_{\nu}$ is well defined and unitary follow from the properties of the Hankel transform, since

\displaystyle\langle\varphi,\psi\rangle_{L^{2}}=\langle r^{-1}\varphi,r^{2}r^{-1}\psi\rangle_{L^{2}}=\langle\mathcal{H}_{\nu}[r^{-1}\varphi],h_{\nu}\mathcal{H}_{\nu}[r^{-1}\psi]\rangle_{L^{2}}=\langle\mathcal{H}_{\nu}[r^{-1}\varphi],\mathcal{H}_{\nu}[r^{-1}\psi]\rangle_{\hat{\mathcal{Q}}_{\nu}}.

We shall use the classical identity

\mathcal{H}_{\nu}[r^{-3/2}J_{\nu}](k)=\frac{1}{2\nu}\left(k^{\nu+1/2}\mathbf{1}_{(0,1]}(k)+k^{-\nu+1/2}\mathbf{1}_{(1,\infty)}(k)\right),

(2.17)

which follows from [17, 6.573.2].

Let

\displaystyle\hat{e}_{\nu,0}(k):=(2\nu)^{-\frac{1}{2}}\big[k^{\nu+\frac{1}{2}}\mathbf{1}_{(0,1]}(k)+k^{-\nu+\frac{1}{2}}\mathbf{1}_{(1,\infty)}(k)\big].

Note that $\hat{e}_{\nu,0}=(2\nu)^{\frac{1}{2}}\mathcal{H}_{\nu}[r^{-\frac{3}{2}}J_{\nu}]=(2\nu)^{\frac{1}{2}}\mathcal{U}_{\nu}(r^{-\frac{1}{2}}J_{\nu})$ , see (2.17). In particular $\hat{e}_{\nu,0}\in\hat{\mathcal{Q}}_{\nu}$ and $\|\hat{e}_{\nu,0}\|^{2}_{\hat{\mathcal{Q}}_{\nu}}=2\nu\|r^{-\frac{1}{2}}J_{\nu}\|_{L^{2}}^{2}=1$ . We have the following orthogonal decomposition of $\hat{\mathcal{Q}}_{\nu}$ (equipped with $\langle\cdot,\cdot\rangle_{\hat{\mathcal{Q}}_{\nu}}$ ):

\displaystyle\hat{\mathcal{Q}}_{\nu}=\hat{\mathcal{Q}}_{\nu}^{0}\oplus\mathrm{Span}(\hat{e}_{\nu,0})\oplus\hat{\mathcal{Q}}_{\nu}^{\infty},

where

\displaystyle\hat{\mathcal{Q}}_{\nu}^{0}:=\{\hat{\varphi}\in\hat{\mathcal{Q}}_{\nu}\,:\,\mathrm{supp}(\hat{\varphi})\subset(0,1]\text{ and }\hat{\varphi}(1)=0\},\quad\hat{\mathcal{Q}}_{\nu}^{\infty}:=\{\hat{\varphi}\in\hat{\mathcal{Q}}_{\nu}\,:\,\mathrm{supp}(\hat{\varphi})\subset[1,\infty)\text{ and }\hat{\varphi}(1)=0\}.

Observe that any $\hat{\varphi}\in\hat{\mathcal{Q}}_{\nu}$ decomposes in this direct sum as

\displaystyle\hat{\varphi}(k)=\big(\hat{\varphi}(k)-\hat{\varphi}(1)k^{\nu+\frac{1}{2}}\big)\mathbf{1}_{(0,1)}(k)+\hat{\varphi}(1)\hat{e}_{\nu,0}(k)+\big(\hat{\varphi}(k)-\hat{\varphi}(1)k^{-\nu+\frac{1}{2}}\big)\mathbf{1}_{(1,\infty)}(k).

Using the unitary operator $\mathcal{U}_{\nu}$ , the orthogonal decomposition (with respect to the $L^{2}$ scalar product) in $L^{2}$ becomes

L^{2}=E_{\nu}^{0}\oplus\operatorname{Span}(e_{\nu,0})\oplus E_{\nu}^{\infty}

with $E_{\nu}^{\#}:=\mathcal{U}_{\nu}^{-1}\hat{\mathcal{Q}}_{\nu}^{\#}$ and $e_{\nu,0}:=\mathcal{U}_{\nu}^{-1}[\hat{e}_{\nu,0}]=(2\nu)^{\frac{1}{2}}p_{\nu}=(2\nu)^{\frac{1}{2}}r^{-\frac{1}{2}}J_{\nu}.$ From the proof of Proposition 2.2, we know that in this orthogonal decomposition we have

\displaystyle\mathrm{Im}\,A_{\nu}^{\pm}=\pm\frac{\pi}{4\nu}(0\oplus 1\oplus 0).

Moreover, Proposition 2.3 yields the real part of the coefficient $a_{00}$ through the identity (2.5).

2.4 The orthonormal basis at low energies

Our first task is to construct an orthogonal basis of $E_{\nu}^{0}$ where the matrix elements of $\mathrm{Re}\,A_{\nu}^{\pm}$ can be computed explicitly.

For all integer $\ell\geq 1$ , we introduce the notation

\displaystyle\hat{\varphi}_{\nu,\ell}(k):=k^{\nu+\frac{1}{2}}(1-k^{2})P_{\ell-1}^{(\nu,1)}(1-2k^{2})\mathbf{1}_{(0,1)}(k),

(2.18)

where $P_{\ell-1}^{(\nu,1)}$ is the usual Jacobi polynomial. We will use the following lemma.

Lemma 2.12.

Let $\nu>0$ . Then, for all $\ell\geq 1$ , $\hat{\varphi}_{\nu,\ell}\in\hat{\mathcal{Q}}_{\nu}^{0}$ . Moreover, for all $0<k<1$ ,

\displaystyle(h_{\nu}\hat{\varphi}_{\nu,\ell})(k)=\frac{4\ell(\ell+\nu)}{1-k^{2}}\hat{\varphi}_{\nu,\ell}(k).

(2.19)

Proof.

The definition (2.18) shows that $\hat{\varphi}_{\nu,\ell}\in\hat{\mathcal{Q}}_{\nu}^{0}$ . We now prove (2.19). Let $0<k<1$ . To shorten the notation, set $P:=P_{\ell-1}^{(\nu,1)}$ . We compute, setting $x=1-2k^{2}$ ,

	$\displaystyle\Big(-\frac{\mathrm{d}^{2}}{\mathrm{d}k^{2}}+\frac{\nu^{2}-\frac{1}{4}}{k^{2}}\Big)k^{\nu+\frac{1}{2}}(1-k^{2})P(1-2k^{2})$
	$\displaystyle=\Big(\frac{1-x}{2}\Big)^{-\frac{1}{2}}\Big(-16\,\Big(\frac{1-x}{2}\Big)^{2}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}+12\,\Big(\frac{1-x}{2}\Big)\frac{\mathrm{d}}{\mathrm{d}x}+(\nu^{2}-\mbox{$\frac{1}{4}$})\Big)\Big(\frac{1-x}{2}\Big)^{\frac{\nu}{2}-\frac{1}{4}}\Big(\frac{1+x}{2}\Big)P(x)$
	$\displaystyle=2^{-\frac{\nu}{2}+\frac{3}{4}}(1-x)^{-\frac{1}{2}}\Big(-2(1-x)^{2}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}+3(1-x)\frac{\mathrm{d}}{\mathrm{d}x}+\frac{1}{2}(\nu^{2}-\mbox{$\frac{1}{4}$})\Big)(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P(x).$		(2.20)

Now,

\displaystyle\frac{\mathrm{d}}{\mathrm{d}x}\big[(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P(x)\big]=(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P^{\prime}+\big((1-x)^{\frac{\nu}{2}-\frac{1}{4}}-(\mbox{$\frac{\nu}{2}-\frac{1}{4}$})(1-x)^{\frac{\nu}{2}-\frac{5}{4}}(1+x)\big)P

and

	$\displaystyle\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}\big[(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P(x)\big]$	$\displaystyle=(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P^{\prime\prime}+\big(2(1-x)^{\frac{\nu}{2}-\frac{1}{4}}-\mbox{$(\nu-\frac{1}{2})$}(1-x)^{\frac{\nu}{2}-\frac{5}{4}}(1+x)\big)P^{\prime}$
		$\displaystyle\quad+\big(\mbox{$\frac{1}{4}(\nu-\frac{1}{2})(\nu-\frac{5}{2})$}(1-x)^{\frac{\nu}{2}-\frac{9}{4}}(1+x)-\mbox{$(\nu-\frac{1}{2})$}(1-x)^{\frac{\nu}{2}-\frac{5}{4}}\big)P.$

Thus

	$\displaystyle\Big(-2(1-x)^{2}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}+3(1-x)\frac{\mathrm{d}}{\mathrm{d}x}+\mbox{$\frac{1}{2}(\nu^{2}-\frac{1}{4})$}\Big)(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P(x)$
	$\displaystyle=-2(1-x)^{\frac{\nu}{2}+\frac{7}{4}}(1+x)P^{\prime\prime}$
	$\displaystyle\quad+\big[-4(1-x)^{\frac{\nu}{2}+\frac{7}{4}}+2(\nu-\mbox{$\frac{1}{2}$})(1-x)^{\frac{\nu}{2}+\frac{3}{4}}(1+x)+3(1-x)^{\frac{\nu}{2}+\frac{3}{4}}(1+x)\big]P^{\prime}$
	$\displaystyle\quad+\big[-\mbox{$\frac{1}{2}(\nu-\frac{1}{2})$}(\nu-\mbox{$\frac{5}{2}$})(1-x)^{\frac{\nu}{2}+\frac{7}{4}}(1+x)+2(\nu-\mbox{$\frac{1}{2}$})(1-x)^{\frac{\nu}{2}+\frac{3}{4}}$
	$\displaystyle\qquad\quad+3(1-x)^{\frac{\nu}{2}+\frac{3}{4}}-3(\mbox{$\frac{\nu}{2}-\frac{1}{4}$})(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)+\mbox{$\frac{1}{2}(\nu^{2}-\frac{1}{4})$}(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)\big]P$
	$\displaystyle=-2(1-x)^{\frac{\nu}{2}+\frac{7}{4}}(1+x)P^{\prime\prime}$
	$\displaystyle\quad+\big[-4(1-x)^{\frac{\nu}{2}+\frac{7}{4}}+2(\nu+1)(1-x)^{\frac{\nu}{2}+\frac{3}{4}}(1+x)\big]P^{\prime}$
	$\displaystyle\quad+2(\nu+1)(1-x)^{\frac{\nu}{2}+\frac{3}{4}}P$
	$\displaystyle=-2(1-x)^{\frac{\nu}{2}+\frac{3}{4}}\Big[(1-x^{2})P^{\prime\prime}+\big[2(1-x)-(\nu+1)(1+x)\big]P^{\prime}-(\nu+1)P\Big]$
	$\displaystyle=-2(1-x)^{\frac{\nu}{2}+\frac{3}{4}}\Big[(1-x^{2})P^{\prime\prime}+\big[1-\nu-(3+\nu)x\big]P^{\prime}-(\nu+1)P\Big]$
	$\displaystyle=-2(1-x)^{\frac{\nu}{2}+\frac{3}{4}}\big[-(\ell-1)(\ell+\nu+1)-(\nu+1)\big]P,$

where in the last equality we have used that the Jacobi polynomial $P$ satisfies the differential equation $(1-x^{2})y^{\prime\prime}+(1-\nu-(\nu+3)x)y^{\prime}=-(\ell-1)(\ell+\nu+1)y$ . Inserting the last equality into (2.20), and using that $(\ell-1)(\ell+\nu+1)+(\nu+1)=\ell(\ell+\nu)$ we obtain

	$\displaystyle\big[h_{\nu}\hat{\varphi}_{\nu,\ell}\big](k)$	$\displaystyle=2^{-\frac{\nu}{2}+\frac{7}{4}}(1-x)^{\frac{\nu}{2}+\frac{1}{4}}\ell(\ell+\nu)P(x)$
		$\displaystyle=4\ell(\ell+\nu)k^{\nu+\frac{1}{2}}P(1-2k^{2})$
		$\displaystyle=4(1-k^{2})^{-1}\ell(\ell+\nu)\hat{\varphi}_{\nu,\ell}(k).$

This concludes the proof of the lemma. ∎

As above we use the notations $p_{\nu}(r)=r^{-\frac{1}{2}}J_{\nu}(r)$ and we set $\tilde{p}_{\nu}:=(2\nu)^{\frac{1}{2}}\,p_{\nu}=e_{\nu,0}$ .

Lemma 2.13.

Let $\nu>0$ . The sequence $(e_{\nu,-\ell})_{\ell\geq 1}$ defined by

\displaystyle\forall\ell\geq 1,\qquad e_{\nu,-\ell}:=\tilde{p}_{\nu+2\ell},

is an orthonormal basis of $E_{\nu}^{0}$ .

Remark 2.14.

We mention that the fact that $(\tilde{p}_{\nu+2\ell})_{\ell\geq 1}$ is an orthonormal family in $L^{2}$ follows, e.g., from [17, 6.574.2]. Likewise, the fact that $r^{-1}\tilde{p}_{\nu+2\ell}\in L^{2}$ for any $\ell\geq 1$ follows from [17, 6.574.2] (since $\nu>0$ and $\ell\geq 1$ ). These two well-known facts also follow from the following proof. The main point of the lemma is that $(\tilde{p}_{\nu+2\ell})_{\ell\geq 1}$ is a basis of $E_{\nu}^{0}$ .

Proof.

We compute $\mathcal{U}_{\nu}[\tilde{p}_{\nu+2\ell}]=\mathcal{H}_{\nu}[r^{-1}\tilde{p}_{\nu+2\ell}]$ . To this end, we write (see [17, 6.574.1])

	$\displaystyle\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}](k)=k^{\frac{1}{2}}\int_{0}^{\infty}r^{-1}J_{\nu}(kr)J_{\nu+2\ell}(r)\mathrm{d}r$	$\displaystyle=\frac{k^{\nu+\frac{1}{2}}\Gamma(\nu+\ell)}{2\Gamma(\ell+1)\Gamma(\nu+1)}\mathbf{F}(\nu+\ell,-\ell;\nu+1;k^{2})\mathbf{1}_{(0,1)}(k)$
		$\displaystyle=\frac{k^{\nu+\frac{1}{2}}}{2(\nu+\ell)}P_{\ell}^{(\nu,-1)}(1-2k^{2})\mathbf{1}_{(0,1)}(k),$

where $\mathbf{F}$ and $P_{\ell}^{(\nu,-1)}$ stand for the usual Hypergeometric function and Jacobi polynomial, respectively. In particular we see that $\mathrm{supp}\,\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}]\subset(0,1]$ and $\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}](1)=0$ for any $\ell\geq 1$ .

Now, using the relation $P_{\ell}^{(\nu,-1)}(x)=\frac{\nu+\ell}{\ell}\frac{x+1}{2}P_{\ell-1}^{(\nu,1)}(x)$ (see [30, (4.22.2)]), we obtain

\displaystyle\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}](k)=\frac{k^{\nu+\frac{1}{2}}(1-k^{2})}{2\ell}P_{\ell-1}^{(\nu,1)}(1-2k^{2})\mathbf{1}_{(0,1)}(k)=\frac{\hat{\varphi}_{\nu,\ell}(k)}{2\ell}.

(2.21)

For any $\ell,\ell^{\prime}\geq 1$ , we can then write

$\displaystyle\big\langle\mathcal{U}_{\nu}p_{\nu+2\ell},\mathcal{U}_{\nu}p_{\nu+2\ell^{\prime}}\big\rangle_{\mathcal{Q}_{\nu}}$	$\displaystyle=\big\langle\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}],\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}^{\prime}]\big\rangle_{\mathcal{Q}_{\nu}}$
	$\displaystyle=\frac{1}{4\ell\ell^{\prime}}\int_{0}^{1}\hat{\varphi}_{\nu,\ell}(k)[h_{\nu}\hat{\varphi}_{\nu,\ell^{\prime}}](k)\mathrm{d}k$
	$\displaystyle=\frac{\ell^{\prime}+\nu}{\ell}\int_{0}^{1}\hat{\varphi}_{\nu,\ell}(k)(k^{2}-1)^{-1}\hat{\varphi}_{\nu,\ell^{\prime}}(k)\mathrm{d}k$
	$\displaystyle=\frac{\ell^{\prime}+\nu}{\ell}\int_{0}^{1}k^{2\nu+1}(k^{2}-1)P_{\ell-1}^{(\nu,1)}(1-2k^{2})P_{\ell^{\prime}-1}^{(\nu,1)}(1-2k^{2})\mathrm{d}k$
	$\displaystyle=\frac{1}{2^{\nu+3}}\frac{\ell^{\prime}+\nu}{\ell}\int_{-1}^{1}(1-x)^{\nu}(1+x)P_{\ell-1}^{(\nu,1)}(x)P_{\ell^{\prime}-1}^{(\nu,1)}(x)\mathrm{d}x$
	$\displaystyle=\frac{\delta_{\ell\ell^{\prime}}}{2(2\ell+\nu)},$	(2.22)

where the last equality follows from the well-known orthogonality relation for Jacobi polynomials (see, e.g., [17, 7.391.1]). Inserting the normalisation $\tilde{p}_{\nu+2\ell}=\sqrt{2(\nu+2\ell)}p_{\nu+2\ell}$ , we deduce that $(\mathcal{U}_{\nu}\tilde{p}_{\nu+2\ell})_{\ell\geq 1}$ is an orthonormal family in $\hat{\mathcal{Q}}_{\nu}^{0}$ , and therefore that $(\tilde{p}_{\nu+2\ell})_{\ell\geq 1}$ is an orthonormal family in $E_{\nu}^{0}$ .

Since in addition it is well known that the family of Jacobi polynomials $(P_{\ell-1}^{(\nu,1)})_{\ell\geq 1}$ forms an orthogonal basis of $L^{2}((-1,1))$ for the weight $(1-x)^{\nu}(1+x)$ , we can conclude from the previous computation that $(\tilde{p}_{\nu+2\ell})_{\ell\geq 1}$ is indeed an orthonormal basis of $E_{\nu}^{0}$ . ∎

Lemma 2.15.

Let $\nu>0$ . Let $\ell,\ell^{\prime}\geq 1$ . Then

\displaystyle\big\langle e_{\nu,-\ell},(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,-\ell^{\prime}}\big\rangle_{L^{2}}=-\frac{\delta_{\ell\ell^{\prime}}}{4\ell(\nu+\ell)}.

Proof.

For all $\varepsilon>0$ , using (2.15) and (2.21), we have

	$\displaystyle\big\langle p_{\nu+2\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu+2\ell^{\prime}}\big\rangle_{L^{2}}$
	$\displaystyle=\mathrm{Re}\,\int_{0}^{1}\frac{1}{k^{2}-(1\pm i\varepsilon)}\overline{\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}](k)}\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell^{\prime}}](k)\mathrm{d}k$
	$\displaystyle=\frac{1}{4\ell\ell^{\prime}}\mathrm{Re}\,\int_{0}^{1}\frac{k^{2\nu+1}(1-k^{2})^{2}}{k^{2}-(1\pm i\varepsilon)}P_{\ell-1}^{(\nu,1)}(1-2k^{2})P_{\ell^{\prime}-1}^{(\nu,1)}(1-2k^{2})\mathrm{d}k.$

Letting $\varepsilon\to 0$ (using Lebesgue’s dominated convergence theorem) yields

	$\displaystyle\lim_{\varepsilon\to 0}\big\langle p_{\nu+2\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu+2\ell^{\prime}}\big\rangle_{L^{2}}$
	$\displaystyle=-\frac{1}{4\ell\ell^{\prime}}\int_{0}^{1}k^{2\nu+1}(1-k^{2})P_{\ell-1}^{(\nu,1)}(1-2k^{2})P_{\ell^{\prime}-1}^{(\nu,1)}(1-2k^{2})\mathrm{d}k.$

Using (2.22), we then obtain

\displaystyle\lim_{\varepsilon\to 0}\big\langle p_{\nu+2\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu+2\ell^{\prime}}\big\rangle=-\frac{\delta_{\ell\ell^{\prime}}}{8\ell(2\ell+\nu)(\ell+\nu)}.

Inserting the normalisation $e_{\nu,-\ell}=\tilde{p}_{\nu+2\ell}=\sqrt{2(\nu+2\ell)}p_{\nu+2\ell}$ , the result follows. ∎

We also need to compute the matrix elements of $\mathrm{Re}\,A_{\nu}^{\pm}$ between the basis elements $e_{\nu,-\ell}$ , $\ell\geq 1$ , of $E_{\nu}^{0}$ and $e_{\nu,0}$ .

Lemma 2.16.

Let $\nu>0$ and $\ell\geq 1$ . Then

\displaystyle\big\langle e_{\nu,\ell},(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}\big\rangle_{L^{2}}=\frac{(-1)^{\ell}}{4\ell(\nu+\ell)}\Big(1+\frac{2\ell}{\nu}\Big)^{\frac{1}{2}}.

Proof.

We proceed as in the proof of Lemma 2.15. For all $\varepsilon>0$ , using (2.15) and (2.21), we have

	$\displaystyle\big\langle p_{\nu+2\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu}\big\rangle_{L^{2}}$
	$\displaystyle=\mathrm{Re}\,\int_{0}^{1}\frac{1}{k^{2}-(1\pm i\varepsilon)}\overline{\mathcal{H}_{\nu}[r^{-1}p_{\nu+2\ell}](k)}\mathcal{H}_{\nu}[r^{-1}p_{\nu}](k)\mathrm{d}k$
	$\displaystyle=\frac{1}{4\nu\ell}\mathrm{Re}\,\int_{0}^{1}\frac{k^{2\nu+1}(1-k^{2})}{k^{2}-(1\pm i\varepsilon)}P_{\ell-1}^{(\nu,1)}(1-2k^{2})\mathrm{d}k.$

Letting $\varepsilon\to 0$ (using Lebesgue’s dominated convergence theorem) yields

	$\displaystyle\lim_{\varepsilon\to 0}\big\langle p_{\nu+2\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu}\big\rangle_{L^{2}}$
	$\displaystyle=-\frac{1}{4\nu\ell}\int_{0}^{1}k^{2\nu+1}P_{\ell-1}^{(\nu,1)}(1-2k^{2})\mathrm{d}k$
	$\displaystyle=-\frac{2^{-\nu-2}}{4\nu\ell}\int_{-1}^{1}(1-x)^{\nu}P_{\ell-1}^{(\nu,1)}(x)\mathrm{d}x=\frac{(-1)^{\ell}}{8\nu\ell(\nu+\ell)},$

where the last equality follows from the change of variables $x\mapsto-x$ , the identity $P_{\ell-1}^{(\nu,1)}(-x)=(-1)^{\ell-1}P_{\ell-1}^{(1,\nu)}(x)$ and [17, 7.391.4]. Inserting the normalisation $e_{\nu,-\ell}=\tilde{p}_{\nu+2\ell}=\sqrt{2(\nu+2\ell)}p_{\nu+2\ell}$ , the result follows. ∎

2.5 The orthonormal basis at high energies

Now we turn to the construction of an orthonormal basis to decompose $A_{\nu}^{\pm}$ in the high energy region $E_{\nu}^{\infty}$ . In this region, it is not clear to us how to define a suitable orthonormal basis with an explicit expression, as in the case of the basis $(\tilde{p}_{\nu+2\ell})_{\ell\geq 1}$ of $E_{\nu}^{0}$ . Nevertheless the Hankel transforms of the basis elements (multiplied by $r^{-1}$ ) that we construct remain explicit, with a form close to that of the previous section. First, for all integer $\ell\geq 1$ and $\nu>0$ , we set

\displaystyle\hat{\psi}_{\nu,\ell}(k):=k^{-\nu+\frac{1}{2}}(1-k^{-2})P_{\ell-1}^{(\nu,1)}(1-2k^{-2})\mathbf{1}_{(1,\infty)}(k).

(2.23)

Lemma 2.17.

Let $\nu>0$ . Then, for all $\ell\geq 1$ , $\hat{\psi}_{\nu,\ell}\in\hat{\mathcal{Q}}_{\nu}^{\infty}$ . Moreover, for all $k>1$ ,

\displaystyle(h_{\nu}\hat{\psi}_{\nu,\ell})(k)=\frac{4\ell(\ell+\nu)}{k^{2}(k^{2}-1)}\hat{\psi}_{\nu,\ell}(k).

(2.24)

Proof.

Let $\ell\geq 0$ . Since $\nu>0$ , the fact that $\hat{\psi}_{\nu,\ell}\in\hat{\mathcal{Q}}_{\nu}^{\infty}$ directly follows from the definition (2.23).

To prove (2.24), we argue as in the proof of Lemma 2.12. Let $k>1$ . To shorten the notations, set $P:=P_{\ell-1}^{(\nu,1)}$ . We compute, setting $x=1-2k^{-2}$ ,

	$\displaystyle\big[h_{\nu}\hat{\psi}_{\nu,\ell}\big](k)=\Big(-\frac{\mathrm{d}^{2}}{\mathrm{d}k^{2}}+\frac{\nu^{2}-\frac{1}{4}}{k^{2}}\Big)k^{-\nu+\frac{1}{2}}(1-k^{-2})P(1-2k^{-2})$
	$\displaystyle=\Big(-16\,\Big(\frac{1-x}{2}\Big)^{3}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}+12\,\Big(\frac{1-x}{2}\Big)^{2}\frac{d}{dx}+\Big(\frac{1-x}{2}\Big)(\nu^{2}-\mbox{$\frac{1}{4}$})\Big)\Big(\frac{1-x}{2}\Big)^{\frac{\nu}{2}-\frac{1}{4}}\Big(\frac{1+x}{2}\Big)P(x)$
	$\displaystyle=2^{-\frac{\nu}{2}-\frac{3}{4}}(1-x)\Big(-2(1-x)^{2}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}+3(1-x)\frac{d}{dx}+\frac{1}{2}(\nu^{2}-\mbox{$\frac{1}{4}$})\Big)(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P(x).$		(2.25)

Now, from the computation in the proof of Lemma 2.12, we have

	$\displaystyle\Big(-2(1-x)^{2}\frac{\mathrm{d}^{2}}{\mathrm{d}x^{2}}+3(1-x)\frac{\mathrm{d}}{\mathrm{d}x}+\frac{1}{2}(\nu^{2}-\mbox{$\frac{1}{4}$})\Big)(1-x)^{\frac{\nu}{2}-\frac{1}{4}}(1+x)P(x)$
	$\displaystyle=2\ell(\ell+\nu)(1-x)^{\frac{\nu}{2}+\frac{3}{4}}P(x).$

Inserting this equality into (2.25), we obtain

	$\displaystyle\big[h_{\nu}\hat{\psi}_{\nu,\ell}\big](k)$	$\displaystyle=2^{-\frac{\nu}{2}+\frac{1}{4}}\ell(\ell+\nu)(1-x)^{\frac{\nu}{2}+\frac{7}{4}}P(x)$
		$\displaystyle=4\ell(\ell+\nu)k^{-\nu-\frac{7}{2}}P(1-2k^{-2})$
		$\displaystyle=4\ell(\ell+\nu)k^{-4}(1-k^{-2})^{-1}k^{-\nu+\frac{1}{2}}(1-k^{-2})P(1-2k^{-2})$
		$\displaystyle=4\ell(\ell+\nu)k^{-2}(k^{2}-1)^{-1}\hat{\psi}_{\nu,\ell}(k).$

This proves the lemma. ∎

Now we set

\displaystyle\psi_{\nu,\ell}:=\mathcal{U}_{\nu}^{-1}\hat{\psi}_{\nu,\ell}.

Lemma 2.18.

Let $\nu>0$ . Then for all $\ell\geq 1$ , $\psi_{\nu,\ell}\in E_{\nu}^{\infty}$ and for all $\ell,\ell^{\prime}\geq 1$ ,

\displaystyle\langle\psi_{\nu,\ell^{\prime}},\psi_{\nu,\ell}\rangle_{L^{2}}=\frac{2\ell^{2}}{2\ell+\nu}\delta_{\ell\ell^{\prime}}.

Proof.

By the definition of $\mathcal{U}_{\nu}$ and $\psi_{\nu,\ell}$ , we have $\psi_{\nu,\ell}\in E_{\nu}^{\infty}$ for all $\ell\geq 1$ , and

\displaystyle\langle\psi_{\nu,\ell^{\prime}},\psi_{\nu,\ell}\rangle_{L^{2}}

\displaystyle=\big\langle\hat{\psi}_{\nu,\ell^{\prime}},\hat{\psi}_{\nu,\ell}\big\rangle_{\hat{\mathcal{Q}}_{\nu}}=\int_{1}^{\infty}\hat{\psi}_{\nu,\ell^{\prime}}(k)\big[h_{\nu}\hat{\psi}_{\nu,\ell}\big](k)\,\mathrm{d}k.

(2.26)

Therefore, using (2.24) and the change of variables $x=1-2k^{-2}$ , we get

$\displaystyle\langle\psi_{\nu,\ell^{\prime}},\psi_{\nu,\ell}\rangle_{L^{2}}$	$\displaystyle=4\ell(\ell+\nu)\int_{1}^{\infty}k^{-2}(k^{2}-1)^{-1}\hat{\psi}_{\nu,\ell^{\prime}}(k)\hat{\psi}_{\nu,\ell}(k)\mathrm{d}k$
	$\displaystyle=4\ell(\ell+\nu)2^{-\nu-3}\int_{-1}^{1}(1-x)^{\nu}(1+x)P^{(\nu,1)}_{\ell^{\prime}-1}(x)P^{(\nu,1)}_{\ell-1}(x)\mathrm{d}x$
	$\displaystyle=2\ell(\ell+\nu)\frac{\ell}{(2\ell+\nu)(\ell+\nu)}\delta_{\ell\ell^{\prime}}.$	(2.27)

Here the last equality follows from the well-known identity for Jacobi polynomials (see, e.g., [17, 7.391.1]).∎

Setting for $\ell\geq 1$ ,

\displaystyle e_{\nu,\ell}:=\Big(\frac{2\ell+\nu}{2\ell^{2}}\Big)^{\frac{1}{2}}\psi_{\nu,\ell},

the previous lemma provides an orthonormal basis for $E_{\nu}^{\infty}$ .

Lemma 2.19.

Let $\nu>0$ . The sequence $(e_{\nu,\ell})_{\ell\geq 1}$ is an orthonormal basis of $E_{\nu}^{\infty}$ .

Proof.

Lemma 2.18 shows that $(e_{\nu,\ell})_{\ell\geq 1}$ is an orthonormal family in $E_{\nu}^{\infty}$ . The fact that it is an orthonormal basis follows from the proof of Lemma 2.18 (see (2.27)) and the well-known property that the family of Jacobi polynomials $(P_{\ell-1}^{(\nu,1)})_{\ell\geq 1}$ is an orthonormal basis of $L^{2}((-1,1))$ for the weight $(1-x)^{\nu}(1+x)$ . ∎

Next we have to compute the matrix elements of $\mathrm{Re}\,A_{\nu}^{\pm}$ in the orthonormal basis of $E_{\nu}^{\infty}$ we just constructed. Unfortunately, it does not yield a diagonal matrix as in the case of $E_{\nu}^{0}$ .

Lemma 2.20.

Let $\nu>0$ . Then for all $\ell,\ell^{\prime}\geq 1$ ,

\displaystyle\big\langle e_{\nu,\ell^{\prime}},(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,\ell}\big\rangle_{L^{2}}=\frac{1}{4\ell\ell^{\prime}}\frac{n!\,\Gamma(\nu+m)}{(m-1)!\,\Gamma(\nu+n+1)}\Big(1+\frac{2\ell}{\nu}\Big)^{\frac{1}{2}}\Big(1+\frac{2\ell^{\prime}}{\nu}\Big)^{\frac{1}{2}},

where $n:=\max(\ell,\ell^{\prime})$ and $m:=\min(\ell,\ell^{\prime})$ .

Proof.

Using Lebesgue’s dominated convergence theorem and the change of variables $x=1-2k^{-2}$ as before, we have

	$\displaystyle\big\langle\psi_{\nu,\ell^{\prime}},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i0))^{-1}r^{-1}\psi_{\nu,\ell}\big\rangle_{L^{2}}$
	$\displaystyle=\lim_{\varepsilon\to 0}\int_{1}^{\infty}\frac{\hat{\psi}_{\nu,\ell^{\prime}}(k)\hat{\psi}_{\nu,\ell}(k)}{k^{2}-(1\pm i\varepsilon)}\mathrm{d}k$
	$\displaystyle=\int_{1}^{\infty}k^{-2\nu-3}(k^{2}-1)P_{\ell^{\prime}-1}^{(\nu,1)}(1-2k^{-2})P_{\ell-1}^{(\nu,1)}(1-2k^{-2})\,\mathrm{d}k$
	$\displaystyle=2^{-\nu-2}\int_{-1}^{1}(1-x)^{\nu-1}(1+x)P^{(\nu,1)}_{\ell^{\prime}-1}(x)P^{(\nu,1)}_{\ell-1}(x)\,\mathrm{d}x.$

To evaluate the integral, let $n:=\max(\ell,\ell^{\prime})$ , $m=\min(\ell,\ell^{\prime})$ . If $m=1$ , the result follows from [17, 7.391.4]. If $m\geq 2$ , write $P^{(\nu,1)}_{m}(x)=P^{(\nu,1)}_{m-1}(1)+(1-x)R_{m-2}(x)$ , where $R_{m-2}$ is some polynomial of degree $m-2$ . Since $m-2<n-1$ , the orthogonality of Jacobi polynomials implies that

	$\displaystyle\int_{-1}^{1}(1-x)^{\nu-1}(1+x)P^{(\nu,1)}_{n-1}(x)P^{(\nu,1)}_{m}(x)\,\mathrm{d}x$	$\displaystyle=\int_{-1}^{1}(1-x)^{\nu-1}(1+x)P^{(\nu,1)}_{n-1}(x)P^{(\nu,1)}_{m-1}(1)\,\mathrm{d}x$
		$\displaystyle=\binom{m-1+\nu}{m-1}\int_{-1}^{1}(1-x)^{\nu-1}(1+x)P^{(\nu,1)}_{n-1}(x)\,\mathrm{d}x.$

Applying [17, 7.391.4], we obtain

\displaystyle\big\langle\psi_{\nu,\ell^{\prime}},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i0))^{-1}r^{-1}\psi_{\nu,\ell}\big\rangle_{L^{2}}

\displaystyle=\frac{1}{2}\binom{m-1+\nu}{m-1}\frac{\Gamma(\nu)\Gamma(n+1)}{\Gamma(\nu+n+1)}=\frac{1}{2\nu}\frac{n!\,\Gamma(\nu+m)}{(m-1)!\,\Gamma(\nu+n+1)}.

We conclude from the normalisation $e_{\nu,\ell}=\frac{\sqrt{2\ell+\nu}}{\ell\sqrt{2}}\psi_{\nu,\ell}$ . ∎

As in the case of $E_{\nu}^{0}$ , we finally need to compute the matrix elements of $\mathrm{Re}\,A_{\nu}^{\pm}$ between the basis elements $e_{\nu,\ell}$ of $E_{\nu}^{\infty}$ and $e_{\nu,0}=\tilde{p}_{\nu}$ .

Lemma 2.21.

Let $\nu>0$ and $\ell\geq 1$ . Then

\displaystyle\big\langle e_{\nu,\ell},(\mathrm{Re}\,A_{\nu}^{\pm})e_{\nu,0}\big\rangle_{L^{2}}=\frac{1}{4\ell(\nu+\ell)}\Big((-1)^{\ell}-\frac{\ell!}{(\nu)_{\ell}}\Big)\Big(1+\frac{2\ell}{\nu}\Big)^{\frac{1}{2}}.

Proof.

We proceed again as in the proof of Lemma 2.15. For all $\varepsilon>0$ , using (2.15) and (2.21), we have

	$\displaystyle\big\langle\psi_{\nu,\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu}\big\rangle_{L^{2}}$
	$\displaystyle=\mathrm{Re}\,\int_{1}^{\infty}\frac{1}{k^{2}-(1\pm i\varepsilon)}\hat{\psi}_{\nu}^{(\ell)}(k)\mathcal{H}_{\nu}[r^{-1}p_{\nu}](k)\,\mathrm{d}k$
	$\displaystyle=\frac{1}{2\nu}\mathrm{Re}\,\int_{1}^{\infty}\frac{k^{-2\nu+1}(1-k^{-2})}{k^{2}-(1\pm i\varepsilon)}P_{\ell-1}^{(\nu,1)}(1-2k^{-2})\,\mathrm{d}k.$

Letting $\varepsilon\to 0$ (using Lebesgue’s dominated convergence theorem) yields

	$\displaystyle\lim_{\varepsilon\to 0}\big\langle\psi_{\nu,\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}r^{-1}p_{\nu}\big\rangle_{L^{2}}$
	$\displaystyle=\frac{1}{2\nu}\int_{1}^{\infty}k^{-2\nu-1}P_{\ell-1}^{(\nu,1)}(1-2k^{-2})\,\mathrm{d}k$
	$\displaystyle=\frac{2^{-\nu-2}}{\nu}\int_{-1}^{1}(1-x)^{\nu-1}P_{\ell-1}^{(\nu,1)}(x)\,\mathrm{d}x=-\frac{1}{4\nu(\nu+\ell)}\Big((-1)^{\ell}-\frac{\ell!}{(\nu)_{\ell}}\Big),$

where the last equality follows from [17, 7.391.2]. From the normalisations $e_{\nu,0}=\tilde{p}_{\nu}=\sqrt{2\nu}p_{\nu}$ and $e_{\nu,\ell}=\frac{\sqrt{2\ell+\nu}}{\ell\sqrt{2}}\psi_{\nu,\ell}$ , the result follows. ∎

2.6 Proof of Theorem 2.1 and Proposition 2.8

Theorem 2.1 is a consequence of the results obtained in the previous sections.

Proof of Theorem 2.1.

It suffices to combine Propositions 2.2 and 2.3, together with Lemmas 2.15, 2.16, 2.20 and 2.21. ∎

Now we turn to the proof of Proposition 2.8 for the operator $B_{\nu}^{\pm}$ defined in (2.12). Clearly, if $\ell\leq-1$ or $\ell^{\prime}\leq-1$ , then

\displaystyle\big\langle e_{\nu,\ell},B_{\nu}^{\pm}e_{\nu,\ell}\big\rangle_{L^{2}}=\big\langle e_{\nu,\ell},A_{\nu}^{\pm}e_{\nu,\ell}\big\rangle_{L^{2}}.

In what follows we compute $\langle e_{\nu,\ell},B_{\nu}^{\pm}e_{\nu,\ell}\rangle_{L^{2}}$ for $\ell\geq 0$ and $\ell^{\prime}\geq 0$ .

Lemma 2.22.

Let $\nu>0$ . Then for all $\ell,\ell^{\prime}\geq 1$ ,

\displaystyle\big\langle e_{\nu,\ell},(\mathrm{Re}\,B_{\nu}^{\pm})e_{\nu,\ell}\big\rangle_{L^{2}}=\frac{\delta_{\ell\ell^{\prime}}}{4\ell(\nu+\ell)}.

Proof.

Arguing as in the proof of the Lemma 2.21 and using (2.24) from Lemma 2.17, we have

	$\displaystyle\big\langle\psi_{\nu,\ell^{\prime}},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i0))^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}\psi_{\nu,\ell}\big\rangle_{L^{2}}$
	$\displaystyle=\lim_{\varepsilon\to 0}\int_{1}^{\infty}\frac{\hat{\psi}_{\nu,\ell^{\prime}}(k)\hat{\psi}_{\nu,\ell}(k)}{k^{2}(k^{2}-(1\pm i\varepsilon))}\,\mathrm{d}k$
	$\displaystyle=\frac{1}{4\ell(\nu+\ell)}\int_{1}^{\infty}\hat{\psi}_{\nu,\ell^{\prime}}(k)\big[h_{\nu}\hat{\psi}_{\nu,\ell}\big](k)\,\mathrm{d}k=\frac{1}{4\ell(\nu+\ell)}\big\langle\psi_{\nu,\ell^{\prime}},\psi_{\nu,\ell}\big\rangle_{L^{2}}.$

The result then follows from Lemma 2.18 and the fact that $e_{\nu,\ell}=\frac{\sqrt{2\ell+\nu}}{\ell\sqrt{2}}\psi_{\nu,\ell}$ . ∎

Lemma 2.23.

Let $\nu>0$ and $\ell\geq 1$ . Then

\displaystyle\big\langle e_{\nu,\ell},(\mathrm{Re}\,B_{\nu}^{\pm})e_{\nu,0}\big\rangle_{L^{2}}=\frac{(-1)^{\ell+1}}{4\ell(\nu+\ell)}\Big(1+\frac{2\ell}{\nu}\Big)^{\frac{1}{2}}.

Proof.

The same computations as in the Lemma 2.21 give

	$\displaystyle\lim_{\varepsilon\to 0}\big\langle\psi_{\nu,\ell},r^{-1}\mathrm{Re}\,(h_{\nu}-(1\pm i\varepsilon))^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}p_{\nu}\big\rangle$
	$\displaystyle=\frac{1}{2\nu}\int_{1}^{\infty}k^{-2\nu-3}P_{\ell-1}^{(\nu,1)}(1-2k^{-2})\,\mathrm{d}k$
	$\displaystyle=\frac{2^{-\nu-3}}{\nu}\int_{-1}^{1}(1-x)^{\nu}P_{\ell-1}^{(\nu,1)}(x)\,\mathrm{d}x=\frac{(-1)^{\ell+1}}{4\nu(\nu+\ell)},$

where the last equality follows from [17, 7.391.2]. From the normalisations $e_{\nu,0}=\sqrt{2\nu}p_{\nu}$ and $e_{\nu,\ell}=\frac{\sqrt{2\ell+\nu}}{\ell\sqrt{2}}\psi_{\nu,\ell}$ , the result follows. ∎

Lemma 2.24.

Let $\nu>0$ . Then

\displaystyle\big\langle e_{\nu,0},r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}e_{\nu,0}\big\rangle_{L^{2}}=\frac{1}{4\nu^{2}}.

Proof.

Recall that

\displaystyle\mathcal{U}_{\nu}[e_{\nu,0}]=\mathcal{H}_{\nu}[r^{-1}e_{\nu,0}](k)=(2\nu)^{-\frac{1}{2}}\big[k^{\nu+\frac{1}{2}}\mathbf{1}_{(0,1]}(k)+k^{-\nu+\frac{1}{2}}\mathbf{1}_{(1,\infty)}(k)\big].

This implies

\displaystyle\big\langle e_{\nu,0},r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}e_{\nu,0}\big\rangle=\frac{1}{2\nu}\int_{1}^{\infty}k^{-2\nu-1}\mathrm{d}k=\frac{1}{4\nu^{2}}.

∎

We can then conclude the proof of Proposition 2.8:

Proof of Proposition 2.8.

It suffices to use the definition (2.12) of $B_{\nu}^{\pm}$ together with Theorem 2.1 and Lemmata 2.22, 2.23 and 2.24. ∎

3 One-channel estimates

Based on the matrix representation of $A_{\nu}^{\pm}$ given by Theorem 2.1, we have already proved parts (i) and (iii) of Theorem 1.1, as well as Theorem 2.4 giving a quantitative lower bound to the boundary value of the weighted resolvent of $h_{\nu}$ . Now we establish additional uniform bounds, in particular parts (ii) and (iv) of Theorem 1.1. The main result of this section is the latter, which we restate here as follows.

Theorem 3.1.

Let $\nu>0$ . One has

\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}=\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}\leq\mathcal{C}(\nu),

(3.1)

where $\mathcal{C}(\nu)$ is defined in (1.8).

The theorem is complemented by uniform lower bounds established in Section 3.4.

3.1 Basic facts

Recall the definition of $h_{\nu}$ with $\nu>0$ given in (1.4). It is clear from the Hardy inequality

\forall\psi\in C_{0}^{\infty}((0,\infty)),\quad\Big\|\frac{\psi}{r}\Big\|_{L^{2}}^{2}\leq 4\,\|\partial_{r}\psi\|_{L^{2}}^{2}.

(3.2)

that $h_{\nu}$ is non-negative. Moreover, its form domain is given by $\mathcal{Q}(h_{\nu})=H^{1}_{0}((0,\infty))$ .

If $\nu\geq 1$ , then the operator domain of $h_{\nu}$ is given by

\mathcal{D}(h_{\nu})=\Big\{\psi\in L^{2}\,:\,-\partial^{2}_{r}\psi+\frac{\nu^{2}-\frac{1}{4}}{r^{2}}\psi\in L^{2}\Big\},

where $-\partial^{2}_{r}\psi+(\nu^{2}-\frac{1}{4})r^{-2}\psi$ should be understood in the sense of distributions. Actually, $h_{\nu}$ is essentially selfadjoint on $C_{0}^{\infty}((0,\infty))$ provided that $\nu\geq 1$ .

If $0<\nu<1$ , then

\mathcal{D}(h_{\nu})=\Big\{\psi\in L^{2}\,:\,-\partial^{2}_{r}\psi+\frac{\nu^{2}-\frac{1}{4}}{r^{2}}\psi\in L^{2}\text{ and }\lim_{r\to 0}\mathcal{W}(\psi,\psi_{0})(r)=0\Big\},

where $\mathcal{W}$ stands for the Wronskian and $\psi_{0}(r):=r^{\frac{1}{2}+\nu}$ . In this case ( $0<\nu<1$ ), the operator initially defined on $C_{0}^{\infty}((0,\infty))$ has several self-adjoint extensions.

In any case, i.e. for all $\nu>0$ , an element of $\mathcal{D}(h_{\nu})$ is a continuous function that behaves like $cr^{\frac{1}{2}+\nu}$ near $0$ , for some $c\in\mathbb{R}$ . Moreover, for all $\nu>0$ , the spectrum of $h_{\nu}$ is $\sigma(h_{\nu})=[0,\infty)$ and $h_{\nu}$ is homogeneous, in the sense that defining, for $\tau>0$ , the unitary operator $u_{\tau}$ in $L^{2}$ , by

(u_{\tau}\psi)(r):=\tau^{1/2}\,\psi(\tau r)\,,

(3.3)

we have

u_{\tau}h_{\nu}u_{\tau}^{-1}=\tau^{-2}h_{\nu}.

(3.4)

For more details on the family of operators $h_{\nu}$ , we refer to [11] and references therein.

3.2 Preliminaries

The Hardy inequality (3.2) has the following direct consequences.

Lemma 3.2.

Let $\nu>0$ . Then, in the sense of quadratic forms in $L^{2}$ defined on $H^{1}((0,\infty))$ ,

r^{-2}\leq\nu^{-2}\,h_{\nu}\qquad\mbox{and}\qquad-\partial_{r}^{2}\leq\max(1,(4\nu^{2})^{-1})\,h_{\nu}.

(3.5)

Proof.

Note that the Hardy inequality (3.2) gives $-\partial_{r}^{2}-\frac{1}{4r^{2}}\geq 0$ . To prove the first inequality in (3.5), we then write

\frac{1}{r^{2}}\leq\frac{1}{\nu^{2}}\Big(\frac{\nu^{2}}{r^{2}}-\partial_{r}^{2}-\frac{1}{4r^{2}}\Big)=\frac{h_{\nu}}{\nu^{2}}.

To prove the second inequality in (3.5), we distinguish two cases. If $\nu^{2}-\frac{1}{4}\geq 0$ , the inequality $-\partial_{r}^{2}\leq h_{\nu}$ is obvious. If $\nu^{2}-\frac{1}{4}<0$ , it suffices to write, using $r^{-2}\leq\nu^{-2}h_{\nu}$ ,

-\partial_{r}^{2}=h_{\nu}+\frac{\frac{1}{4}-\nu^{2}}{r^{2}}\leq h_{\nu}+\frac{\frac{1}{4}-\nu^{2}}{\nu^{2}}h_{\nu}=(4\nu^{2})^{-1}h_{\nu}.

This proves the lemma. ∎

The Hardy inequality (3.2) gives the following easy bound for the weighted Bessel resolvent away from the spectrum $[0,\infty)$ .

Lemma 3.3.

Let $\nu>0$ and $z\in\mathbb{C}\setminus[0,\infty)$ . Then

\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}\leq\nu^{-2}\Big(1+\frac{|z|}{\mathrm{dist}(z,[0,\infty))}\Big).

Proof.

Estimate (3.5) implies that

\big\|r^{-1}h_{\nu}^{-\frac{1}{2}}\big\|_{L^{2}\to L^{2}}=\big\|h_{\nu}^{-\frac{1}{2}}r^{-1}\big\|_{L^{2}\to L^{2}}\leq\nu^{-1}.

(3.6)

Together with the estimate

\left\|h_{\nu}(h_{\nu}-z)^{-1}\right\|_{L^{2}\to L^{2}}\leq 1+|z|\left\|(h_{\nu}-z)^{-1}\right\|_{L^{2}\to L^{2}}\leq 1+\frac{|z|}{\mathrm{dist}(z,[0,\infty))},

this proves the lemma. ∎

For $\mathrm{Re}\,z\leq 0$ , the same Hardy argument gives the following sharper uniform estimate.

Lemma 3.4.

Let $\nu>0$ and $z\in\mathbb{C}$ be such that $\mathrm{Re}\,z\leq 0$ , $z\neq 0$ . Then

\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}\leq\nu^{-2}.

(3.7)

Proof.

Using (3.5) and $\mathrm{Re}\,z\leq 0$ , we have, for all $\psi\in\mathcal{D}(h_{\nu})$ ,

\displaystyle\nu^{2}\Big\|\frac{\psi}{r}\Big\|_{L^{2}}^{2}\leq\langle\psi,h_{\nu}\psi\rangle_{L^{2}}-\mathrm{Re}\,z\,\|\psi\|_{L^{2}}^{2}=\mathrm{Re}\,\langle\psi,(h_{\nu}-z)\psi\rangle_{L^{2}}\leq\Big\|\frac{\psi}{r}\Big\|_{L^{2}}\,\big\|r(h_{\nu}-z)\psi\big\|_{L^{2}}\,,

which yields

\Big\|\frac{\psi}{r}\Big\|_{L^{2}}\leq\nu^{-2}\big\|r(h_{\nu}-z)\psi\big\|_{L^{2}}.

Let $\varphi\in C_{0}^{\infty}((0,\infty))$ . Applying the previous inequality to $\psi=(h_{\nu}-z)^{-1}r^{-1}\varphi$ (which belongs to $\mathcal{D}(h_{\nu})$ ), we obtain

\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\varphi\right\|_{L^{2}}\leq\nu^{-2}\|\varphi\|_{L^{2}}.

Since $C_{0}^{\infty}((0,\infty))$ is dense in $L^{2}$ , this concludes the proof. ∎

The next proposition is the standard maximum principle reduction from the resolvent set to the boundary rays of a sector around the positive half–line. To keep the one channel argument self-contained, we spell out the proof after recording the two scalar input lemmata.

Proposition 3.5.

Let $\nu>0$ . The map

(0,\infty)\ni\varepsilon\mapsto\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}

is non-increasing and we have

\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}=\lim_{\varepsilon\to 0}\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}.

For $\varepsilon>0$ , we write

S_{\varepsilon}:=\mathbb{C}\setminus\{z\in\mathbb{C}\,:\,|\mathrm{Im}\,z|\leq\varepsilon\mathrm{Re}\,z\}.

The following two lemmata employ an explicit expression of the integral kernel of $(h_{\nu}-z)^{-1}$ .

Lemma 3.6.

Let $\nu>0$ and let $\phi,\psi\in C_{0}^{\infty}((0,\infty))$ . Then

\lim_{\mathrel{{\mathop{|z|\to\infty}\limits_{z\not\in[0,\infty)}}}}\left\langle\phi,(h_{\nu}-z)^{-1}\psi\right\rangle_{L^{2}}=0.

Proof.

For $z\in\mathbb{C}\setminus[0,\infty)$ , let $(r,R)\mapsto k_{\nu,z}(r,R)$ denote the integral kernel of $(h_{\nu}-z)^{-1}$ . Then $k_{\nu,z}$ satisfies, for all $r,R>0$ ,

|k_{\nu,z}(r,R)|\leq\frac{C_{\nu}}{|z|^{1/2}}\min(1,|z|^{\frac{1}{4}}r^{\frac{1}{2}})\min(1,|z|^{\frac{1}{4}}R^{\frac{1}{2}}),

(3.8)

for some constant $C_{\nu}>0$ , see, e.g., [11, Prop. 4.3]. One deduces from this that, for any $\phi,\psi\in\mathrm{C}_{0}^{\infty}((0,\infty))$ and $|z|$ large enough,

\left|\left\langle\phi,(h_{\nu}-z)^{-1}\psi\right\rangle\right|\leq\frac{C_{\nu}}{|z|^{1/2}}\|\phi\|_{L^{\infty}}\|\psi\|_{L^{\infty}},

which proves the lemma. ∎

Lemma 3.7.

Let $\varepsilon>0$ , $\nu>0$ , and let $\phi,\psi\in C_{0}^{\infty}((0,\infty))$ . Then

\lim_{\mathrel{{\mathop{z\to 0}\limits_{z\in S_{\varepsilon}}}}}\left\langle\phi,(h_{\nu}-z)^{-1}\psi\right\rangle_{L^{2}}=\left\langle\phi,h_{\nu}^{-1}\psi\right\rangle_{L^{2}}.

Proof.

Estimate (3.8) supplies the required domination, so Lebesgue’s dominated convergence theorem applies. ∎

Proof of Proposition 3.5.

Since $S_{\varepsilon_{2}}\subset S_{\varepsilon_{1}}$ whenever $0<\varepsilon_{1}<\varepsilon_{2}$ , the map

\varepsilon\mapsto\sup_{z\in S_{\varepsilon}}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}

is non-increasing. Also every fixed $z_{0}\in\mathbb{C}\setminus[0,\infty)$ belongs to $S_{\varepsilon}$ for all sufficiently small $\varepsilon>0$ . Hence

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}=\lim_{\varepsilon\to 0}\sup_{z\in S_{\varepsilon}}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}.

It remains to replace $S_{\varepsilon}$ by its boundary. Let $\phi,\psi\in C_{0}^{\infty}((0,\infty))$ . The function

z\mapsto\left\langle\phi,r^{-1}(h_{\nu}-z)^{-1}r^{-1}\psi\right\rangle_{L^{2}}

is analytic on $S_{\varepsilon}$ and continuous up to the boundary; continuity at the vertex $z=0$ follows from Lemma 3.7, applied to $r^{-1}\phi$ and $r^{-1}\psi$ . By Lemma 3.6, applied to the same two test functions, this scalar function tends to zero as $|z|\to\infty$ . The maximum modulus principle on truncated sectors, followed by the limits at zero and infinity, gives

\sup_{z\in S_{\varepsilon}}\left|\left\langle\phi,r^{-1}(h_{\nu}-z)^{-1}r^{-1}\psi\right\rangle_{L^{2}}\right|=\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left|\left\langle\phi,r^{-1}(h_{\nu}-z)^{-1}r^{-1}\psi\right\rangle_{L^{2}}\right|.

The density of $C_{0}^{\infty}((0,\infty))$ in $L^{2}$ then gives the corresponding equality for the operator norms, and the proposition follows. ∎

Remark 3.8.

It follows from scaling invariance that the norm $\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\|_{L^{2}\to L^{2}}$ is constant along the ray $e^{i\arg(z)}(0,\infty)$ . More precisely, recalling the definition (3.3) of the unitary transform $u_{\tau}$ on $L^{2}$ , we have, by (3.4),

	$\displaystyle\left\\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\\|_{L^{2}\to L^{2}}$	$\displaystyle=\left\\|U_{\tau}r^{-1}(h_{\nu}-z)^{-1}r^{-1}U_{\tau}^{-1}\right\\|_{L^{2}\to L^{2}}$
		$\displaystyle=\left\\|r^{-1}(h_{\nu}-\tau^{2}z)^{-1}r^{-1}\right\\|_{L^{2}\to L^{2}}\,,$

for any $\tau>0$ . Together with Proposition 3.5, this justifies the equality in the statement of Theorem 3.1, namely

\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}=\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}.

3.3 Proof of Theorem 3.1

Now we are in a position to prove Theorem 3.1. The proof is based on the method of multipliers. An identity obtained from this method in the context of the Bessel operators reads as follows.

Lemma 3.9.

Let $\nu>0$ and $z\in\mathbb{C}\setminus[0,\infty)$ be such that $\mathrm{Re}\,z>0$ . Let $f\in C_{0}^{\infty}((0,\infty))$ and let $\psi:=(h_{\nu}-z)^{-1}f$ . Then

	$\displaystyle\\|\partial_{r}\psi^{-}\\|_{L^{2}}^{2}+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\int_{0}^{\infty}r\big\|\partial_{r}\psi^{-}\big\|^{2}$
	$\displaystyle=\mathrm{Re}\,\int_{0}^{\infty}\Big(\frac{\|z_{2}\|}{\sqrt{z_{1}}}rf^{-}\overline{\psi^{-}}-\frac{\|z_{2}\|}{\sqrt{z_{1}}}(\nu^{2}-\mbox{$\frac{1}{4}$})\frac{\|\psi^{-}\|^{2}}{r}+2rf^{-}\overline{\partial_{r}\psi^{-}}-2\frac{\nu^{2}-\frac{1}{4}}{r}\psi^{-}\overline{\partial_{r}\psi^{-}}\Big),$		(3.9)

where

z_{1}:=\mathrm{Re}\,z,\quad z_{2}:=\mathrm{Im}\,z,\quad\psi^{-}(r):=e^{-i\sqrt{z_{1}}\mathrm{sgn}(z_{2})r}\psi(r),

and likewise for $f^{-}$ .

Proof.

The identity can be deduced from [8, Rem. 3.1]. For completeness, and because the half–line setting involves both a singular Bessel potential and the boundary point $r=0$ , we give the short derivation in the present context.

The weak formulation of $(h_{\nu}-z)\psi=f$ reads

\int_{0}^{\infty}\bar{\phi}^{\prime}\psi^{\prime}=z\int_{0}^{\infty}\bar{\phi}\psi+\int_{0}^{\infty}\bar{\phi}\tilde{f}\qquad\mbox{with}\qquad\tilde{f}:=f-\frac{\nu^{2}-\frac{1}{4}}{r^{2}}\,\psi\,,

(3.10)

where $\phi$ is any function from the form domain of $h_{\nu}$ . The test functions used below are admissible because $\psi$ is smooth and exponentially decaying at infinity (this follows, e.g., from the explicit expression of the integral kernel of $(h_{\nu}-z)^{-1}$ ), while $\psi(r)\sim cr^{\frac{1}{2}+\nu}$ as $r\to 0$ for some constant $c$ . This asymptotic also eliminates all boundary terms at $r=0$ .

Choosing $\phi:=r\psi^{\prime}+\frac{1}{2}\psi$ in (3.10), taking twice the real part and integrating by parts gives

2\int_{0}^{\infty}|\psi^{\prime}|^{2}=-2z_{2}\,\mathrm{Im}\,\int_{0}^{\infty}\psi\,r\overline{\psi^{\prime}}+2\mathrm{Re}\,\int_{0}^{\infty}\tilde{f}\,r\overline{\psi^{\prime}}+\mathrm{Re}\,\int_{0}^{\infty}\tilde{f}\,\overline{\psi}.

Choosing $\phi:=\psi$ and taking the real part gives

\int_{0}^{\infty}|\psi^{\prime}|^{2}=z_{1}\int_{0}^{\infty}|\psi|^{2}+\mathrm{Re}\,\int_{0}^{\infty}\tilde{f}\,\overline{\psi}.

Finally, choosing $\phi:=r\psi$ and taking respectively the real and imaginary parts gives

\int_{0}^{\infty}r|\psi^{\prime}|^{2}=z_{1}\int_{0}^{\infty}r|\psi|^{2}+\mathrm{Re}\,\int_{0}^{\infty}r\tilde{f}\,\overline{\psi}

and

\mathrm{Im}\,\int_{0}^{\infty}\overline{\psi}\,\psi^{\prime}=z_{2}\int_{0}^{\infty}r|\psi|^{2}+\mathrm{Im}\,\int_{0}^{\infty}r\overline{\psi}\,\tilde{f}.

No singular term appears on the left-hand side of the third identity.

We now take the first identity minus the second one, add $|z_{2}|z_{1}^{-1/2}$ times the third one, and subtract $2z_{1}^{1/2}\mathop{\mathrm{sgn}}\nolimits(z_{2})$ times the fourth one. Substituting $\tilde{f}=f-(\nu^{2}-\frac{1}{4})r^{-2}\psi$ and using

\partial_{r}\psi^{-}=e^{-i\sqrt{z_{1}}\mathop{\mathrm{sgn}}\nolimits(z_{2})r}\big(\psi^{\prime}-i\sqrt{z_{1}}\mathop{\mathrm{sgn}}\nolimits(z_{2})\psi\big),\qquad f^{-}=e^{-i\sqrt{z_{1}}\mathop{\mathrm{sgn}}\nolimits(z_{2})r}f,

we arrive at the desired identity (3.9). ∎

Proof of Theorem 3.1.

Let $z\in\mathbb{C}\setminus[0,\infty)$ . If $z_{1}=\mathrm{Re}\,z\leq 0$ (and $z\neq 0$ ), we know from Lemma 3.4 that

\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}\leq\nu^{-2}.

Now we assume that $z_{1}>0$ . Let $f$ and $\psi$ be as in Lemma 3.9. We see that the second term on the left-hand side of (3.9) is non-negative. The first two terms on the right-hand side are controlled by rough Cauchy bounds, since they only affect the auxiliary constants. When $\nu^{2}-\frac{1}{4}\geq 0$ , the second term on the right-hand side of (3.9) is non-positive and can be dropped. More precisely, using the Cauchy–Schwarz inequality and

|z_{2}|\,\|\psi\|_{L^{2}}^{2}=\Big|\mathrm{Im}\,\int_{0}^{\infty}\overline{\psi}f\Big|\leq\Big\|\frac{\psi}{r}\Big\|_{L^{2}}\|rf\|_{L^{2}},

we obtain

	$\displaystyle\Big\|\int_{0}^{\infty}\Big(\frac{\|z_{2}\|}{\sqrt{z_{1}}}rf^{-}\overline{\psi^{-}}-\frac{\|z_{2}\|}{\sqrt{z_{1}}}(\nu^{2}-\mbox{$\frac{1}{4}$})\frac{\|\psi^{-}\|^{2}}{r}\Big)\Big\|$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\\|rf\\|_{L^{2}}^{\frac{3}{2}}\Big\\|\frac{\psi}{r}\Big\\|^{\frac{1}{2}}_{L^{2}}+\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big\\|\frac{\psi}{r}\Big\\|^{\frac{3}{2}}_{L^{2}}\\|rf\\|_{L^{2}}^{\frac{1}{2}}\Big)$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\mbox{$\frac{3}{4}$}\\|rf\\|^{2}_{L^{2}}+\mbox{$\frac{1}{4}$}\Big\\|\frac{\psi}{r}\Big\\|^{2}_{L^{2}}+\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big(\mbox{$\frac{3}{4}$}\Big\\|\frac{\psi}{r}\Big\\|^{2}_{L^{2}}+\frac{1}{4}\\|rf\\|^{2}_{L^{2}}\Big)\Big)$
	$\displaystyle=\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\Big(\mbox{$\frac{3}{4}$}+\mbox{$\frac{1}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\\|rf\\|^{2}_{L^{2}}+\Big(\mbox{$\frac{1}{4}$}+\mbox{$\frac{3}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\Big\\|\frac{\psi}{r}\Big\\|^{2}_{L^{2}}\Big)$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\Big(\mbox{$\frac{3}{4}$}+\mbox{$\frac{1}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\\|rf\\|^{2}_{L^{2}}+\nu^{-2}\Big(\mbox{$\frac{1}{4}$}+\mbox{$\frac{3}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|^{2}_{L^{2}}\Big),$		(3.11)

where we used $\|r^{-1}\psi\|_{L^{2}}=\|r^{-1}\psi^{-}\|_{L^{2}}$ and Lemma 3.2 in the last inequality. For the last term on the right-hand side of (3.9), we integrate by parts, obtaining

2\mathrm{Re}\,\int_{0}^{\infty}\frac{\nu^{2}-\frac{1}{4}}{r}\psi^{-}\overline{\partial_{r}\psi^{-}}=\int_{0}^{\infty}\frac{\nu^{2}-\frac{1}{4}}{r}\partial_{r}|\psi|^{2}=\Big(\nu^{2}-\frac{1}{4}\Big)\int_{0}^{\infty}\frac{|\psi|^{2}}{r^{2}}.

(3.12)

Inserting (3.11) and (3.12) into (3.9), and estimating the third term on the right-hand side of (3.9) by the Cauchy–Schwarz inequality, we obtain

	$\displaystyle\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}\nu^{-2}\Big(\mbox{$\frac{1}{4}$}+\mbox{$\frac{3}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\Big)\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|_{L^{2}}^{2}$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\Big(\mbox{$\frac{3}{4}$}+\mbox{$\frac{1}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\\|rf\\|^{2}_{L^{2}}+2\\|rf\\|_{L^{2}}\\|\partial_{r}\psi^{-}\\|_{L^{2}}\Big).$

To simplify the notations, we rewrite the previous inequality as

\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)\big\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\|_{L^{2}}^{2}-2\|rf\|_{L^{2}}\big\|\partial_{r}\psi^{-}\big\|_{L^{2}}-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{2,\nu}\|rf\|^{2}_{L^{2}}\leq 0,

(3.13)

where $C_{j,\nu}$ are the positive explicit constants depending on $\nu$ following from the equations above.

Assume first that $\nu\leq\frac{1}{\sqrt{3}}$ . From (3.13) and Lemma 3.2, we obtain

\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)\big\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\|_{L^{2}}^{2}-2\max(1,(2\nu)^{-1})\|rf\|_{L^{2}}\big\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\|_{L^{2}}-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{2,\nu}\|rf\|^{2}_{L^{2}}\leq 0.

Computing the roots of the previous quadratic expression in $\|h_{\nu}^{\frac{1}{2}}\psi^{-}\|_{L^{2}}$ , this yields

\displaystyle\big\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\|_{L^{2}}\leq\frac{2\max(1,(2\nu)^{-1})+\Big[4\max(1,(2\nu)^{-2})+4\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu})\Big)\sqrt{\frac{|z_{2}|}{z_{1}}}C_{2,\nu}\Big]^{\frac{1}{2}}}{2\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)}\|rf\|_{L^{2}}.

Using again Lemma 3.2, we arrive at

\displaystyle\Big\|\frac{\psi}{r}\Big\|_{L^{2}}\leq\frac{2\max(1,(2\nu)^{-1})+\Big[2\max(1,(4\nu)^{-2})+4\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu})\Big)\sqrt{\frac{|z_{2}|}{z_{1}}}C_{2,\nu}\Big]^{\frac{1}{2}}}{2\nu\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)}\|rf\|_{L^{2}}.

(3.14)

This inequality holds for any $\psi$ of the form $\psi=(h_{\nu}-z)^{-1}f$ with $f\in C_{0}^{\infty}((0,\infty))$ . Applying it to $\psi=(h_{\nu}-z)^{-1}r^{-1}\varphi$ with $\varphi\in C_{0}^{\infty}((0,\infty))$ (and recalling from Lemma 3.3 that $r^{-1}(h_{\nu}-z)^{-1}r^{-1}$ extends to a bounded operator on $L^{2}((0,\infty))$ ), we obtain

	$\displaystyle\big\\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\varphi\big\\|_{L^{2}}$
	$\displaystyle\leq\frac{2\max(1,(2\nu)^{-1})+\Big[4\max(1,(2\nu)^{-2})+4\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{1,\nu})\Big)\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{2,\nu}\Big]^{\frac{1}{2}}}{2\nu\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{1,\nu}\Big)}\\|\varphi\\|_{L^{2}}$		(3.15)

for any $\varphi\in C_{0}^{\infty}((0,\infty))$ . Since $C_{0}^{\infty}((0,\infty))$ is dense in $L^{2}$ , estimate (3.15) leads to

		$\displaystyle\big\\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\big\\|_{L^{2}\to L^{2}}$
		$\displaystyle\leq\frac{2\max(1,(2\nu)^{-1})+\Big[4\max(1,(2\nu)^{-2})+4\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{1,\nu})\Big)\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{2,\nu}\Big]^{\frac{1}{2}}}{2\nu\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{1,\nu}\Big)}.$		(3.16)

Applying now Proposition 3.5, we deduce that

		$\displaystyle\sup_{z_{1}>0,\,z_{2}\neq 0}\big\\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\big\\|$
		$\displaystyle\leq\lim_{\varepsilon\to 0}\sup_{z_{2}=\varepsilon z_{1}}\frac{2\max(1,(2\nu)^{-1})+\Big[4\max(1,(2\nu)^{-2})+4\Big(1-\sqrt{\varepsilon}\,C_{1,\nu})\Big)\sqrt{\varepsilon}\,C_{2,\nu}\Big]^{\frac{1}{2}}}{2\nu\Big(1-\sqrt{\varepsilon}\,C_{1,\nu}\Big)}$
		$\displaystyle=\max\Big(\frac{1}{\nu^{2}},\frac{2}{\nu}\Big).$		(3.17)

This proves (3.1) in the case where $\nu\leq\frac{1}{\sqrt{3}}$ .

Suppose now that $\nu>\frac{1}{\sqrt{3}}$ . Going back to (3.13), we have

	$\displaystyle\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{1,\nu}\Big)^{2}\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|_{L^{2}}^{4}$
	$\displaystyle\leq\Big(2\\|rf\\|_{L^{2}}\big\\|\partial_{r}\psi^{-}\big\\|_{L^{2}}-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{2,\nu}\\|rf\\|^{2}_{L^{2}}\Big)^{2}$
	$\displaystyle\leq 4\\|rf\\|^{2}_{L^{2}}\big\\|\partial_{r}\psi^{-}\big\\|^{2}_{L^{2}}+\frac{\|z_{2}\|}{z_{1}}C_{2,\nu}^{2}\\|rf\\|^{4}_{L^{2}}$
	$\displaystyle=4\\|rf\\|^{2}_{L^{2}}\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|^{2}_{L^{2}}-4\Big(\nu^{2}-\mbox{$\frac{1}{4}$}\Big)\Big\\|\frac{\psi}{r}\Big\\|_{L^{2}}^{2}\\|rf\\|_{L^{2}}^{2}+\frac{\|z_{2}\|}{z_{1}}C_{2,\nu}^{2}\\|rf\\|^{4}_{L^{2}}.$

Computing the discriminant of the previous quadratic expression in $\|h_{\nu}^{\frac{1}{2}}\psi^{-}\|_{L^{2}}^{2}$ , we obtain

\displaystyle 16\|rf\|^{4}_{L^{2}}-4\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)^{2}\Big(4\Big(\nu^{2}-\mbox{$\frac{1}{4}$}\Big)\Big\|\frac{\psi}{r}\Big\|_{L^{2}}^{2}\|rf\|_{L^{2}}^{2}-\frac{|z_{2}|}{z_{1}}C_{2,\nu}^{2}\|rf\|^{4}_{L^{2}}\Big)\geq 0,

which implies

\displaystyle 16\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)^{2}\Big(\nu^{2}-\mbox{$\frac{1}{4}$}\Big)\Big\|\frac{\psi}{r}\Big\|_{L^{2}}^{2}\leq\Big(16+4\Big(1-\sqrt{\frac{|z_{2}|}{z_{1}}}C_{1,\nu}\Big)^{2}\frac{|z_{2}|}{z_{1}}C_{2,\nu}^{2}\Big)\|rf\|^{2}_{L^{2}}.

(3.18)

We can then proceed exactly as before, using (3.18) instead of (3.14), which leads to (3.1) in the case where $\nu>\frac{1}{\sqrt{3}}$ . This concludes the proof of the theorem. ∎

3.4 One-channel lower bounds

We record the channel lower bounds that complement the upper estimate in Theorem 3.1 and the lower bound in Theorem 2.4. Using a logarithmic trial state near the origin, the latter can be improved for small values of $\nu$ , giving, for any $z\in\mathbb{C}\setminus[0,\infty)$ , the lower bound $\nu^{-2}$ for the norm of the weighted resolvent of $h_{\nu}$ . This in turn yields part (ii) of Theorem 1.1.

Proposition 3.10.

Let $\nu>0$ . Then, for all $z\in\mathbb{C}\setminus[0,\infty)$ ,

\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}\geq\frac{1}{\nu^{2}}.

(3.19)

Proof.

Let $z\in\mathbb{C}\setminus[0,\infty)$ . Let $\chi\in C_{c}^{\infty}([0,\infty))$ be such that $\mathbf{1}_{[0,1]}\leq\chi$ , and let $\eta\in C^{\infty}((0,\infty))$ be such that $\eta=0$ on $(0,\frac{1}{2}]$ and $\eta=1$ on $[1,\infty)$ . We set $\eta_{\sigma}(r):=\eta(\frac{r}{\sigma})$ . For $0<\sigma<1$ , define

u_{\sigma}(r):=r^{\frac{1}{2}}\chi(r)\eta_{\sigma}(r).

Since $\chi\eta_{\sigma}\in C^{\infty}_{c}((0,\infty))$ , we have $u_{\sigma}\in H^{2}((0,\infty))$ and

\partial_{r}^{2}u_{\sigma}=-\frac{u_{\sigma}}{4r^{2}}+r^{-\frac{1}{2}}(\chi\eta_{\sigma})^{\prime}+r^{\frac{1}{2}}(\chi\eta_{\sigma})^{\prime\prime}.

This implies that $u_{\sigma}\in\mathcal{D}(h_{\nu})$ and that $(h_{\nu}-z)u_{\sigma}\in\mathcal{D}(r)$ with

r(h_{\nu}-z)u_{\sigma}=\nu^{2}\frac{u_{\sigma}}{r}-zru_{\sigma}+r^{\frac{1}{2}}(\chi\eta_{\sigma})^{\prime}+r^{\frac{3}{2}}(\chi\eta_{\sigma})^{\prime\prime}.

(3.20)

For $0<\sigma<1$ and $z\in\mathbb{C}\setminus[0,\infty)$ , we set

v_{\sigma}(z):=r(h_{\nu}-z)u_{\sigma},\qquad\text{so that}\qquad u_{\sigma}=(h_{\nu}-z)^{-1}r^{-1}v_{\sigma}(z).

We claim that

\frac{\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}v_{\sigma}(z)\|_{L^{2}}}{\|v_{\sigma}(z)\|_{L^{2}}}=\frac{\|r^{-1}u_{\sigma}\|_{L^{2}}}{\|r(h_{\nu}-z)u_{\sigma}\|_{L^{2}}}\geq\frac{1}{\nu^{2}-C_{\chi,\eta}(1+|z|)(\ln\sigma)^{-1}},

(3.21)

for some positive constant $C_{\chi,\eta}$ depending on $\chi$ and $\eta$ but independent of $\sigma$ and $z$ . To prove (3.21), we first observe that

\|r^{-1}u_{\sigma}\|_{L^{2}}^{2}=\int_{0}^{\infty}r^{-1}(\chi\eta_{\sigma})^{2}\geq\int_{\sigma}^{1}r^{-1}=-\ln\sigma>0.

(3.22)

Next we estimate each term from the right-hand side of (3.20). First,

\|zru_{\sigma}\|^{2}_{L^{2}}\leq|z|^{2}\int_{0}^{\infty}r^{3}(\chi\eta_{\sigma})^{2}\leq C_{\chi,\eta}|z|^{2}.

Moreover,

\big\|r^{\frac{1}{2}}(\chi\eta_{\sigma})^{\prime}\big\|_{L^{2}}^{2}\leq C_{\chi,\eta}+C_{\chi,\eta}\sigma^{-2}\int_{\frac{\sigma}{2}}^{\sigma}r\leq C_{\chi,\eta},

and a similar computation gives

\big\|r^{\frac{3}{2}}(\chi\eta_{\sigma})^{\prime\prime}\big\|_{L^{2}}^{2}\leq C_{\chi,\eta}.

Letting $A_{\sigma}(z):=zru_{\sigma}+r^{\frac{1}{2}}(\chi\eta_{\sigma})^{\prime}+r^{\frac{3}{2}}(\chi\eta_{\sigma})^{\prime\prime}$ , we have shown that $\|A_{\sigma}(z)\|_{L^{2}}\leq C_{\chi,\eta}(1+|z|)$ . Using (3.20) and (3.22), we can write

	$\displaystyle\frac{\\|r^{-1}u_{\sigma}\\|_{L^{2}}}{\\|r(h_{\nu}-z)u_{\sigma}\\|_{L^{2}}}$	$\displaystyle\geq\frac{\\|r^{-1}u_{\sigma}\\|_{L^{2}}}{\nu^{2}\\|r^{-1}u_{\sigma}\\|_{L^{2}}+\\|A_{\sigma}(z)\\|_{L^{2}}}$
		$\displaystyle=\frac{1}{\nu^{2}+\\|A_{\sigma}(z)\\|_{L^{2}}\\|r^{-1}u_{\sigma}\\|_{L^{2}}^{-1}}$
		$\displaystyle\geq\frac{1}{\nu^{2}-C_{\chi,\eta}(1+\|z\|)(\ln\sigma)^{-1}}.$

This establishes (3.21), which in turn yields that, for all $0<\sigma<1$ and $z\in\mathbb{C}\setminus[0,\infty)$ ,

\big\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\big\|_{L^{2}\to L^{2}}\geq\frac{1}{\nu^{2}-C_{\chi,\eta}(1+|z|)(\ln\sigma)^{-1}}.

Letting $\sigma\to 0$ concludes the proof. ∎

Remark 3.11.

Together with Lemma 3.4, Proposition 3.10 implies that, for all $\nu>0$ and $z\in\mathbb{C}$ such that $\mathrm{Re}\,z\leq 0$ , $z\neq 0$ ,

\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}=\frac{1}{\nu^{2}}.

(3.23)

Combining Proposition 3.10 and Theorem 3.1, we deduce part (ii) of Theorem 1.1:

Corollary 3.12.

Let $0<\nu\leq\frac{1}{2}$ . One has

\|A_{\nu}^{\pm}\|_{L^{2}\to L^{2}}=\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}(h_{\nu}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}=\frac{1}{\nu^{2}}.

(3.24)

4 Application to the multidimensional operators

4.1 The Laplacian

We now pass from the one-channel Bessel operators in $L^{2}((0,\infty))$ to the Laplacian in $L^{2}(\mathbb{R}^{d})$ by the usual partial wave decomposition. As a consequence of Theorem 3.1, we obtain Theorem 1.2 for the $d$ -dimensional Laplacian with $d\geq 3$ . Let us abbreviate

N_{d}(z):=\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}.

The main ingredient is the following proposition, which is a direct consequence of the well-known expression of the Laplacian in spherical coordinates $(r,\sigma)\in\mathbb{R}_{+}\times\mathbb{S}^{d-1}$ , with $\mathbb{R}_{+}:=(0,\infty)$ . Note also that (4.2) establishes Proposition 1.3 stated in the introduction.

Proposition 4.1.

Let $d\geq 3$ . For all $z\in\mathbb{C}\setminus[0,\infty)$ , we have

\displaystyle N_{d}(z)=\sup_{\ell\geq 0}\big\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})},

(4.1)

where $\nu_{\ell,d}:=\ell+\frac{d}{2}-1$ . In particular, the maps

d\mapsto N_{2d+1}(z)\quad\text{ and }\quad d\mapsto N_{2d+2}(z)

(4.2)

are non-increasing on $\mathbb{N}^{*}:=\{1,2,3,\dots\}$ .

Proof.

By using the spherical coordinates $(r,\sigma)$ , let

U_{1}:L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}_{+},r^{d-1}\mathrm{d}r)\otimes L^{2}(\mathbb{S}^{d-1},\mathrm{d}\sigma)

be the unitary operator defined by $(U_{1}f)(r,\sigma):=f(r\sigma)$ and recall that

U_{1}\Delta U_{1}^{*}=\partial_{r}^{2}+\frac{d-1}{r}\partial_{r}+\frac{1}{r^{2}}\Delta_{\mathbb{S}^{d-1}},

(4.3)

in the sense of quadratic forms in $L^{2}(\mathbb{R}^{d})$ defined on $C_{0}^{\infty}(\mathbb{R}^{d})$ , where $\Delta_{\mathbb{S}^{d-1}}$ stands for the Laplace–Beltrami operator on the sphere $\mathbb{S}^{d-1}$ . Denote by $\lambda_{\ell}:=-\ell^{2}-\ell(d-2)$ , $\ell\geq 0$ , the eigenvalues of $\Delta_{\mathbb{S}^{d-1}}$ and let $K_{\ell}$ be the associated eigenspaces.

Let

U_{2}:L^{2}(\mathbb{R}_{+},r^{d-1}\mathrm{d}r)\otimes L^{2}(\mathbb{S}^{d-1},\mathrm{d}\sigma)\to\bigoplus_{\ell=0}^{\infty}L^{2}(\mathbb{R}_{+},r^{d-1}\mathrm{d}r)\otimes K_{\ell}

(4.4)

be defined by $U_{2}(f\otimes\varphi)=\bigoplus_{\ell=0}^{\infty}(f\otimes\pi_{K_{\ell}}\varphi)$ and extended by linearity, where $\pi_{K_{\ell}}$ stands for the orthogonal projection onto $K_{\ell}$ in $L^{2}(\mathbb{S}^{d-1},\mathrm{d}\sigma)$ . Then $U_{2}$ is a unitary operator and

U_{2}\Big(\partial_{r}^{2}+\frac{d-1}{r}\partial_{r}+\frac{1}{r^{2}}\Delta_{\mathbb{S}^{d-1}}\Big)U_{2}^{*}=\bigoplus_{\ell=0}^{\infty}\Big(\partial_{r}^{2}+\frac{d-1}{r}\partial_{r}+\frac{\lambda_{\ell}}{r^{2}}\Big)\otimes 1_{K_{\ell}}.

Now let

u_{3}:L^{2}(\mathbb{R}_{+},r^{d-1}\mathrm{d}r)\to L^{2}(\mathbb{R}_{+},\mathrm{d}r)=L^{2}(\mathbb{R}_{+})

be the unitary operator defined by $(u_{3}f)(r)=r^{(d-1)/2}f$ . Then

u_{3}\Big(\partial_{r}^{2}+\frac{d-1}{r}\partial_{r}+\frac{\lambda_{\ell}}{r^{2}}\Big)u_{3}^{*}=\partial_{r}^{2}-\frac{\nu_{\ell}^{2}-\frac{1}{4}}{r^{2}},\qquad\nu_{\ell,d}:=\ell+\frac{d}{2}-1,

for all $\ell\geq 0$ .

Finally, let

U:=\Big(\bigoplus_{\ell=0}^{\infty}u_{3}\otimes 1_{K_{\ell}}\Big)U_{2}U_{1}:L^{2}(\mathbb{R}^{d})\to\bigoplus_{\ell=0}^{\infty}L^{2}(\mathbb{R}_{+})\otimes K_{\ell}.

(4.5)

Taking the Friedrichs extensions, we deduce from the previous identities that

U(-\Delta)U^{*}=\bigoplus_{\ell=0}^{\infty}h_{\nu_{\ell,d}}\otimes 1_{K_{\ell}},

and therefore

|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}=U\Big(\bigoplus_{\ell=0}^{\infty}r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\otimes 1_{K_{\ell}}\Big)U^{*}.

(4.6)

This yields

	$\displaystyle N_{d}(z)$	$\displaystyle=\sup_{\ell\geq 0}\big\\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\otimes 1_{K_{\ell}}\big\\|_{L^{2}(\mathbb{R}_{+})\otimes K_{\ell}\to L^{2}(\mathbb{R}_{+})\otimes K_{\ell}}$
		$\displaystyle=\sup_{\ell\geq 0}\big\\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}.$		(4.7)

The fact that the maps in (4.2) are nonincreasing directly follows from the last equality since $\nu_{\ell,d}=\ell+\frac{d}{2}-1$ . ∎

Theorem 1.2 now follows by combining Proposition 4.1 with Theorem 3.1.

Proof of Theorem 1.2.

For any $d\geq 3$ , we deduce from Proposition 4.1 that

\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}\\ =\sup_{\ell\geq 0}\big\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}.

(4.8)

Now if $d=3$ , then $\nu_{0,3}=\frac{1}{2}$ and $\nu_{\ell,3}\geq\frac{3}{2}>\frac{1}{\sqrt{3}}$ for $\ell\geq 1$ . Thus Theorem 3.1 gives

\displaystyle\sup_{\ell\geq 0}\big\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}\leq\max\Big(\frac{2}{\nu_{0,3}},\sup_{\ell\geq 1}\frac{1}{\sqrt{\nu_{\ell,3}^{2}-\frac{1}{4}}}\Big)=4.

Likewise, if $d\geq 4$ , then Theorem 3.1 gives

\displaystyle\sup_{\ell\geq 0}\big\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}\leq\sup_{\ell\geq 0}\frac{1}{\sqrt{\nu_{\ell,d}^{2}-\frac{1}{4}}}=\frac{1}{\sqrt{\nu_{0,d}^{2}-\frac{1}{4}}}=\frac{2}{\sqrt{(d-1)(d-3)}}.

This concludes the proof of Theorem 1.2. ∎

Similarly, Proposition 3.10 gives the lower bound in the lowest channel for any $z\in\mathbb{C}\setminus[0,\infty)$ .

Proposition 4.2.

Let $d\geq 3$ and let $z\in\mathbb{C}\setminus[0,\infty)$ . Then

\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}\geq\left(\frac{d}{2}-1\right)^{-2}.

(4.9)

Proof.

Let $z\in\mathbb{C}\setminus[0,\infty)$ . Using (4.8) and Proposition 3.10, we obtain

	$\displaystyle\left\\|\|x\|^{-1}(-\Delta-z)^{-1}\|x\|^{-1}\right\\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}$	$\displaystyle=\sup_{\ell\geq 0}\big\\|r^{-1}\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}$
		$\displaystyle\geq\sup_{\ell\geq 0}\frac{1}{\nu_{\ell,d}^{2}}=\frac{1}{\nu_{0,d}^{2}}=\left(\frac{d}{2}-1\right)^{-2}.$

This concludes the proof of Proposition 4.2. ∎

By combining Proposition 4.2 and Theorem 1.2, we deduce the sharp result (1.10) in the three-dimensional case. In dimension $4$ we have the bounds

1\leq\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{4})\to L^{2}(\mathbb{R}^{4})}\leq\frac{2}{\sqrt{3}}\approx 1.15.

(4.10)

A similar consequence of the argument used in the proof of Proposition 4.1 is the following remark regarding the two dimensional case.

Remark 4.3.

Suppose that $d=2$ . Restricting the weighted resolvent to the orthogonal of the vectors with zero angular momentum, we have the following uniform resolvent estimate:

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\Pi^{\perp}\right\|_{L^{2}(\mathbb{R}^{2})\to L^{2}(\mathbb{R}^{2})}\leq\frac{2}{\sqrt{3}},

where $\Pi$ is the orthogonal projection onto $U^{*}(L^{2}(\mathbb{R}_{+})\otimes K_{0})U$ (see the proof of Proposition 4.1 for the notations, in particular (4.5) for the definition of $U$ ) and $\Pi^{\perp}=1-\Pi$ .

The lower bound given by Proposition 4.2 can be improved for large $d$ , by the same argument, using the lower bound for $\|A_{\nu}^{\pm}\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}$ obtained in Theorem 2.4.

Proposition 4.4.

Let $d\geq 3$ . One has

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}\geq\left(\frac{d}{2}-1\right)^{-1}\left(\frac{\pi^{2}}{12}+\frac{1}{8}\left(\frac{d}{2}-1\right)^{-2}\right)^{\frac{1}{2}}.

The estimates of Propositions 4.2 or 4.4 and Theorem 1.2 do not match in the limit of large dimensions. To obtain the correct asymptotics, one can use the asymptotics of $\|A_{\nu}^{\pm}\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}$ obtained in Theorem 2.5.

Proposition 4.5.

One has

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}=\frac{2}{d}+O(d^{-3/2}),\qquad d\to\infty.

(4.11)

Proof.

By (3.1) and (4.8), we have

\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}

\displaystyle=\sup_{\ell\geq 0}\big\|A_{\nu_{\ell,d}}^{\pm}\big\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}.

Using Theorem 1.2, we deduce that

\displaystyle\big\|A_{\nu_{0,d}}^{\pm}\big\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}\leq\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}

\displaystyle\leq\frac{2}{\sqrt{(d-1)(d-3)}},

and therefore, by (2.11),

\displaystyle\frac{1}{\nu_{0,d}}-\frac{C}{\nu_{0,d}^{3/2}}\leq\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}

\displaystyle\leq\frac{2}{\sqrt{(d-1)(d-3)}},

for some positive constant $C$ . Since $\nu_{0,d}=\frac{d-1}{2}$ , this implies (4.11). ∎

Finally, the same argument as in the proof of Proposition 4.2 combined with Proposition 2.2 also allows us to recover Simon’s identity (1.11) proved originally in [31].

Proof of (1.11).

One has

	$\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\\|\mathrm{Im}\,\|x\|^{-1}(-\Delta-z)^{-1}\|x\|^{-1}\right\\|_{L^{2}\to L^{2}}$	$\displaystyle=\sup_{z\in\mathbb{C}\setminus[0,\infty)}\sup_{\ell\geq 0}\big\\|r^{-1}\mathrm{Im}\,\big(h_{\nu_{\ell,d}}-z\big)^{-1}r^{-1}\big\\|_{L^{2}\to L^{2}}$
		$\displaystyle=\frac{\pi}{4\nu_{0,d}}=\frac{\pi}{2(d-2)},$		(4.12)

which establishes the desired result. ∎

4.2 Inverse square potentials

The same partial wave computation applies without change to the operator $H_{d,c}$ introduced in (1.6). The restriction $c>-c_{d}$ is exactly the condition to guarantee that the lowest effective Bessel order be positive.

Proof of Theorem 1.5.

We use the notation and the unitary transform $U$ from the proof of Proposition 4.1. Since the multiplication operator $|x|^{-2}$ is radial, it preserves the spherical harmonic channels. After the radial conjugation by $u_{3}$ , the $\ell$ -th channel becomes

-\partial_{r}^{2}+\frac{\nu_{\ell,d,c}^{2}-\frac{1}{4}}{r^{2}}=h_{\nu_{\ell,d,c}},

where $\nu_{\ell,d,c}^{2}=\nu_{\ell,d}^{2}+c.$ Thus, for all $z\in\mathbb{C}\setminus[0,\infty)$ ,

	$\displaystyle\left\\|\|x\|^{-1}(H_{d,c}-z)^{-1}\|x\|^{-1}\right\\|_{L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})}$
	$\displaystyle\qquad=\sup_{\ell\geq 0}\left\\|r^{-1}(h_{\nu_{\ell,d,c}}-z)^{-1}r^{-1}\right\\|_{L^{2}(\mathbb{R}_{+})\to L^{2}(\mathbb{R}_{+})}.$

The sequence $\ell\mapsto\nu_{\ell,d,c}$ is increasing, and the function $\mathcal{C}$ defined in (1.8) is decreasing. Theorem 3.1 therefore gives the desired bound (1.14).

Finally, if $0<\nu_{0,d,c}\leq 1/2$ , the lower bound in Proposition 3.10, applied to the lowest channel, gives

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(H_{d,c}-z)^{-1}|x|^{-1}\right\|\geq\frac{1}{\nu_{0,d,c}^{2}}.

This is exactly $\mathcal{C}(\nu_{0,d,c})$ in the range $0<\nu_{0,d,c}\leq 1/2$ . ∎

Finally, we establish identity (1.12) for the imaginary part of the weighted resolvent from Theorem 1.4.

Proof of (1.12).

The same channel decomposition gives, for every $z\in\mathbb{C}\setminus[0,\infty)$ ,

\displaystyle|x|^{-1}\mathrm{Im}\,(H_{d,c}-z)^{-1}|x|^{-1}=U\Big(\bigoplus_{\ell=0}^{\infty}r^{-1}\mathrm{Im}\,(h_{\nu_{\ell,d,c}}-z)^{-1}r^{-1}\otimes 1_{K_{\ell}}\Big)U^{*}.

Therefore Proposition 2.2 yields

\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}\mathrm{Im}\,(H_{d,c}-z)^{-1}|x|^{-1}\right\|=\sup_{\ell\geq 0}\frac{\pi}{4\nu_{\ell,d,c}}=\frac{\pi}{4\nu_{0,d,c}},

because $\ell\mapsto\nu_{\ell,d,c}$ is increasing. ∎

4.3 The Aharonov–Bohm operators

Now we establish Theorems 1.6 and Theorem 1.4 for the two-dimensional magnetic operator $-\Delta_{\alpha}$ introduced in (1.7).

Proof of Theorem 1.6.

We proceed as in the proof of Theorem 1.2 in Section 4.1. In polar coordinates $(r,\sigma)\in\mathbb{R}_{+}\times\mathbb{S}^{1}$ , the Aharonov–Bohm operator reads (cf. (4.3))

U_{1}(-\Delta_{\alpha})U_{1}^{*}=-\partial_{r}^{2}-\frac{1}{r}\partial_{r}+\frac{1}{r^{2}}(-i\partial_{\sigma}-\alpha)^{2}.

The spectrum of $(-i\partial_{\sigma}-\alpha)^{2}$ in $L^{2}(\mathbb{S}^{1})$ with domain $H^{1}(\mathbb{S}^{1})$ is purely discrete. The eigenvalues and normalised eigenfunctions read $\{(m-\alpha)^{2}\}_{m\in\mathbb{Z}}$ and $\{\phi_{m}(\sigma):=(2\pi)^{-1/2}e^{im\sigma}\}_{m\in\mathbb{Z}}$ . In the definition (4.4) of the unitary transform $U_{2}$ , let us replace the eigenspaces $K_{\ell}$ by the one dimensional subspaces spanned by $\phi_{m}$ and let sum over $m\in\mathbb{Z}$ instead of $\ell\in\mathbb{N}$ . Similarly, let us perform the same replacements in the definition (4.5) of the unitary transform $U$ . Then (cf. (4.6))

|x|^{-1}(-\Delta_{\alpha}-z)^{-1}|x|^{-1}=U\Big(\bigoplus_{m\in\mathbb{Z}}r^{-1}\big(h_{\nu_{m}(\alpha)}-z\big)^{-1}r^{-1}\otimes 1_{\{\phi_{m}\}}\Big)U^{*},

where

\nu_{m}(\alpha):=|m-\alpha|\,.

Consequently (cf. (4.8)),

\left\||x|^{-1}(-\Delta_{\alpha}-z)^{-1}|x|^{-1}\right\|_{L^{2}(\mathbb{R}^{2})\to L^{2}(\mathbb{R}^{2})}=\sup_{m\in\mathbb{Z}}\big\|r^{-1}\big(h_{\nu_{m}(\alpha)}-z\big)^{-1}r^{-1}\big\|_{L^{2}\to L^{2}}.

(4.13)

This identity is an analogue of Proposition 4.1.

To apply Theorem 3.1, notice that, for any $\alpha\not\in\mathbb{Z}$ , there exists a unique $m_{\alpha}\in\mathbb{Z}$ such that $\alpha\in(m_{\alpha},m_{\alpha+1})$ . For all integer $m<m_{\alpha}$ or $m>m_{\alpha+1}$ , one has $\nu_{m}(\alpha)>1>\frac{1}{\sqrt{3}}$ , so the upper bound of (3.1) with $\nu:=\nu_{m}(\alpha)$ is less than $\frac{1}{\sqrt{1-\frac{1}{4}}}<4$ . If $\mathrm{dist}(\alpha,m_{\alpha})<\mathrm{dist}(\alpha,m_{\alpha+1})$ , then $\nu_{m_{\alpha+1}}(\alpha)\in(\frac{1}{2},1)$ , so the upper bound of (3.1) with $\nu:=\nu_{m_{\alpha+1}}(\alpha)$ (respectively, with $\nu:=\nu_{m_{\alpha}}(\alpha)$ ) is less than $\frac{2}{\frac{1}{2}}=4$ or $\frac{1}{\sqrt{\frac{1}{3}-\frac{1}{4}}}<4$ (respectively, equals $\frac{1}{\mathrm{dist}(m_{\alpha},\alpha)^{2}}>4$ ). If $\mathrm{dist}(\alpha,m_{\alpha})>\mathrm{dist}(\alpha,m_{\alpha+1})$ , the role of $m_{\alpha}$ and $m_{\alpha+1}$ is exchanged in the previous argument. Finally, if $\mathrm{dist}(\alpha,m_{\alpha})=\mathrm{dist}(\alpha,m_{\alpha+1})$ , then $\nu_{m_{\alpha}}(\alpha)=\frac{1}{2}=\nu_{m_{\alpha+1}}(\alpha)$ , so the upper bound of (3.1) with $\nu:=\nu_{m_{\alpha+1}}(\alpha)=\nu_{m_{\alpha+1}}(\alpha)$ is given by $\frac{1}{\mathrm{dist}(m_{\alpha},\alpha)^{2}}=4$ . In summary, applying Theorem 3.1, we get

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\big\|r^{-1}\big(h_{\nu_{m}(\alpha)}-z\big)^{-1}r^{-1}\big\|_{L^{2}\to L^{2}}\leq\frac{1}{\mathrm{dist}(\alpha,\mathbb{Z})^{2}}

Since this bound is independent of $m$ , together with equality (4.13) it implies the upper bound of Theorem 1.6.

To prove the lower bound of Theorem 1.6, we combine (4.13) and Proposition 3.10 as follows:

	$\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\sup_{m\in\mathbb{Z}}\big\\|r^{-1}\big(h_{\nu_{m}(\alpha)}-z\big)^{-1}r^{-1}\big\\|_{L^{2}\to L^{2}}$	$\displaystyle\geq\sup_{z\in\mathbb{C}\setminus[0,\infty)}\big\\|r^{-1}\big(h_{\nu_{m_{*}}(\alpha)}-z\big)^{-1}r^{-1}\big\\|_{L^{2}\to L^{2}}$
		$\displaystyle\geq\frac{1}{\nu_{m_{*}(\alpha)}^{2}}\,,$

where $m_{*}\in\mathbb{Z}$ is such that $\mathrm{dist}(\alpha,\mathbb{Z})=|m_{*}-\alpha|$ . This together with equality (4.13) implies the lower bound of Theorem 1.6. ∎

Finally, we establish identity (1.13) for the imaginary part of the weighted resolvent from Theorem 1.4.

Proof of (1.13).

Using the decomposition in the proof of Theorem 1.6, we have

\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}\mathrm{Im}\,(-\Delta_{\alpha}-z)^{-1}|x|^{-1}\right\|=\sup_{m\in\mathbb{Z}}\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\|r^{-1}\mathrm{Im}\,(h_{\nu_{m}(\alpha)}-z)^{-1}r^{-1}\right\|_{L^{2}\to L^{2}}.

Proposition 2.2 gives

\sup_{m\in\mathbb{Z}}\frac{\pi}{4\nu_{m}(\alpha)}=\frac{\pi}{4\,\mathrm{dist}(\alpha,\mathbb{Z})},

which proves (1.13). ∎

5 The multiplier method for the Laplacian

This final section records a direct route for proving Theorem 1.2. It is kept because it represents an alternative way to establish the quantitative uniform resolvent estimate, which admits a generalisation to potential perturbations which are not necessarily radial.

5.1 Preliminaries

We begin by justifying that for any $z\in\mathbb{C}\setminus[0,\infty)$ , $|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}$ indeed extends to a bounded operator on $L^{2}:=L^{2}(\mathbb{R}^{d})$ .

Lemma 5.1.

Let $d\geq 3$ and $z\in\mathbb{C}\setminus[0,\infty)$ . Then

\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}\leq\Big(\frac{2}{d-2}\Big)^{2}\Big(1+\frac{|z|}{\mathrm{dist}(z,[0,\infty))}\Big).

Proof.

It suffices to combine Hardy’s inequality (valid for every $d\geq 3$ )

\forall\psi\in H^{1}(\mathbb{R}^{d})\,,\qquad\|\nabla\psi\|_{L^{2}}^{2}\geq\left(\frac{d-2}{2}\right)^{2}\big\||x|^{-1}\psi\big\|_{L^{2}}\,,

(5.1)

which implies

\big\||x|^{-1}(-\Delta)^{-\frac{1}{2}}\big\|_{L^{2}\to L^{2}}=\big\|(-\Delta)^{-\frac{1}{2}}|x|^{-1}\big\|_{L^{2}\to L^{2}}\leq\frac{2}{d-2},

(5.2)

and the estimate

\left\|(-\Delta)(-\Delta-z)^{-1}\right\|_{L^{2}\to L^{2}}\leq 1+|z|\left\|(-\Delta-z)^{-1}\right\|_{L^{2}\to L^{2}}\leq 1+\frac{|z|}{\mathrm{dist}(z,[0,\infty))}.

∎

Next we establish another easy result showing that if $\mathrm{Re}\,z\leq 0$ , $z\neq 0$ , the norm of the weighted resolvent $|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}$ is bounded by the inverse of the constant in Hardy’s inequality (5.1).

Lemma 5.2.

Let $d\geq 3$ and $z\in\mathbb{C}$ be such that $\mathrm{Re}\,z\leq 0$ , $z\neq 0$ . Then

\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}\leq\Big(\frac{2}{d-2}\Big)^{2}.

(5.3)

Proof.

Writing, for all $\psi\in H^{2}(\mathbb{R}^{d})$ ,

\displaystyle\|\nabla\psi\|_{L^{2}}^{2}\leq\|\nabla\psi\|_{L^{2}}^{2}-\mathrm{Re}\,z\,\|\psi\|_{L^{2}}^{2}=\mathrm{Re}\,\langle\psi,(-\Delta-z)\psi\rangle_{L^{2}}\leq\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\,\big\||x|(-\Delta-z)\psi\big\|_{L^{2}}\,,

we deduce from Hardy’s inequality (5.1) that

\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\leq\Big(\frac{2}{d-2}\Big)^{2}\big\||x|(-\Delta-z)\psi\big\|_{L^{2}}.

Let $\varphi\in C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ . Applying the previous inequality to $\psi=(-\Delta-z)^{-1}|x|^{-1}\varphi$ (which belongs to $H^{2}(\mathbb{R}^{d})$ ), we obtain

\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\varphi\right\|_{L^{2}\to L^{2}}\leq\Big(\frac{2}{d-2}\Big)^{2}\|\varphi\|_{L^{2}}.

Since $C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ is dense in $L^{2}$ , this concludes the proof. ∎

In order to estimate the norm of $|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}$ for $\mathrm{Re}\,z>0$ , $\mathrm{Im}\,z\neq 0$ , we will use the following proposition which will be proven thanks to the maximum modulus principle.

Proposition 5.3.

Let $d\geq 3$ . The map

(0,\infty)\ni\varepsilon\mapsto\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}

(5.4)

is non-increasing and we have

\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}=\lim_{\varepsilon\to 0}\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}.

(5.5)

The proof of Proposition 5.3 relies on the following well-known lemmata.

Lemma 5.4.

Let $\varphi,\psi\in C_{0}^{\infty}(\mathbb{R}^{d})$ . Then

\lim_{\mathrel{{\mathop{|z|\to\infty}\limits_{z\not\in[0,\infty)}}}}\left\langle\varphi,(-\Delta-z)^{-1}\psi\right\rangle_{L^{2}}=0.

Proof.

If $\mathrm{dist}\big(z,\sigma(-\Delta)\big)\to\infty$ , then the result immediately follows from (1.1). To cover the other situations at the same time, let us consider $z=k^{2}$ with $\mathrm{Im}\,k>0$ . Inspired by [12, Sec. 3.1.1], we employ the representation formula

\frac{1}{a^{2}-k^{2}}=\int_{0}^{\infty}\frac{\sin(at)}{a}\,e^{ikt}\,\mathrm{d}t

valid with every $a>0$ . We apply the formula with the help of the Fourier transform as follows:

	$\displaystyle\|k\|\,\left\|\left\langle\varphi,(-\Delta-z)^{-1}\psi\right\rangle_{L^{2}}\right\|$	$\displaystyle=\|k\|\,\left\|\left\langle\hat{\varphi},(\|\cdot\|^{2}-k^{2})^{-1}\hat{\psi}\right\rangle\right\|$
		$\displaystyle=\left\|\int_{\mathbb{R}^{d}}\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\int_{0}^{\infty}\frac{\sin(\|\xi\|t)}{\|\xi\|}\,k\,e^{ikt}\,\mathrm{d}t\,\mathrm{d}\xi\right\|$
		$\displaystyle=\left\|\int_{0}^{\infty}\int_{\mathbb{R}^{d}}\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\,\cos(\|\xi\|t)\,e^{ikt}\,\mathrm{d}\xi\,\mathrm{d}t\right\|\,,$

where the second line follows by an integration by parts in $t$ and Fubini’s theorem. Note that $\hat{\varphi},\hat{\psi}$ belong to the Schwartz space $\mathscr{S}(\mathbb{R}^{d})$ . Splitting the integration in $t$ over $(0,1)$ and $(1,\infty)$ , the first integral is bounded by a constant independent of $k$ :

\displaystyle\left|\int_{0}^{1}\int_{\mathbb{R}^{d}}\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\,\cos(|\xi|t)\,e^{ikt}\,\mathrm{d}\xi\,\mathrm{d}t\right|

\displaystyle\leq\|\hat{\varphi}\|_{L^{2}(\mathbb{R}^{d})}\,\|\hat{\psi}\|_{L^{2}(\mathbb{R}^{d})}\int_{0}^{1}\mathrm{d}t

\displaystyle=\|\varphi\|_{L^{2}(\mathbb{R}^{d})}\,\|\psi\|_{L^{2}(\mathbb{R}^{d})}\,.

For the second integral, we integrate by parts (twice) with respect to $|\xi|$ as follows:

	$\displaystyle\left\|\int_{1}^{\infty}\int_{\mathbb{R}^{d}}\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\,\cos(\|\xi\|t)\,e^{ikt}\,\mathrm{d}\xi\,\mathrm{d}t\right\|$	$\displaystyle=\left\|\int_{1}^{\infty}\int_{\mathbb{R}^{d}}\left(\frac{\xi}{\|\xi\|}\cdot\nabla\right)^{2}\left[\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\right]\frac{\sin(\|\xi\|t)}{t^{2}}\,e^{ikt}\,\mathrm{d}\xi\,\mathrm{d}t\right\|$
		$\displaystyle\leq\\|\hat{\varphi}\\|_{H^{2}(\mathbb{R}^{d})}\,\\|\hat{\psi}\\|_{H^{2}(\mathbb{R}^{d})}\int_{1}^{\infty}\frac{\mathrm{d}t}{t^{2}}=\\|\hat{\varphi}\\|_{H^{2}(\mathbb{R}^{d})}\,\\|\hat{\psi}\\|_{H^{2}(\mathbb{R}^{d})}\,.$

In summary, there exists a positive constant $C$ independent of $z$ such that

\left|\left\langle\varphi,(-\Delta-z)^{-1}\psi\right\rangle_{L^{2}}\right|\leq\frac{C}{\sqrt{|z|}}

for every $z\in\mathbb{C}\setminus[0,\infty)$ . (In fact, $C\leq\|\hat{\varphi}\|_{H^{2}(\mathbb{R}^{d})}\,\|\hat{\psi}\|_{H^{2}(\mathbb{R}^{d})}$ .) This, in particular, implies the claim of the lemma. ∎

For any $0<\varepsilon<1$ , we set

S_{\varepsilon}:=\mathbb{C}\setminus\{z\in\mathbb{C}\,|\,|\mathrm{Im}\,z|\leq\varepsilon\,\mathrm{Re}\,z\}.

(5.6)

Lemma 5.5.

Let $\varphi,\psi\in C_{0}^{\infty}(\mathbb{R}^{d})$ . Then

\lim_{\mathrel{{\mathop{z\to 0}\limits_{z\in S_{\varepsilon}}}}}\left\langle\varphi,(-\Delta-z)^{-1}\psi\right\rangle_{L^{2}}=\left\langle\phi,(-\Delta)^{-1}\psi\right\rangle_{L^{2}}.

(5.7)

Proof.

Let $z\in S_{\varepsilon}$ and write $z_{1}:=\mathrm{Re}\,z$ , $z_{2}:=\mathrm{Im}\,z$ . We have

\displaystyle\left\langle\varphi,(-\Delta-z)^{-1}\psi\right\rangle_{L^{2}}=\int_{\mathbb{R}^{d}}\big(|\xi|^{2}-z\big)^{-1}\overline{\hat{\varphi}(\xi)}\hat{\psi}(\xi)\,\mathrm{d}\xi.

(5.8)

In order to apply Lebesgue’s dominated convergence theorem, we bound $||\xi|^{2}-z|^{-1}$ as follows: If $z_{1}\leq 0$ , then

\big||\xi|^{2}-z\big|^{-1}=\big(\big(|\xi|^{2}-z_{1}\big)^{2}+z_{2}^{2}\big)^{-\frac{1}{2}}\leq|\xi|^{-2}.

(5.9)

If $0<z_{1}<\varepsilon^{-1}z_{2}$ , we have

	$\displaystyle\big\|\|\xi\|^{2}-z\big\|^{2}$	$\displaystyle=\|\xi\|^{4}-2z_{1}\|\xi\|^{2}+z_{1}^{2}+z_{2}^{2}$
		$\displaystyle\geq\|\xi\|^{4}-2z_{1}\|\xi\|^{2}+(1+\varepsilon^{2})z_{1}^{2}$
		$\displaystyle=\Big[\big(1+\varepsilon^{2}\big)^{-\frac{1}{2}}\|\xi\|^{2}-\big(1+\varepsilon^{2}\big)^{\frac{1}{2}}z_{1}\Big]^{2}+\varepsilon^{2}\big(1+\varepsilon^{2}\big)^{-1}\|\xi\|^{4}$
		$\displaystyle\geq\varepsilon^{2}\big(1+\varepsilon^{2}\big)^{-1}\|\xi\|^{4},$

and therefore

\big||\xi|^{2}-z\big|^{-1}\leq\big(1+\varepsilon^{2}\big)^{-\frac{1}{2}}\varepsilon^{-1}|\xi|^{-2}.

(5.10)

Using (5.9), (5.10) and the fact that $\xi\mapsto|\xi|^{-1}\hat{\varphi}(\xi)$ , $\xi\mapsto|\xi|^{-1}\hat{\psi}(\xi)$ are square integrable (by Hardy’s inequality (5.1), since $\hat{\varphi}$ , $\hat{\psi}$ belong to the Schwartz space $\mathscr{S}(\mathbb{R}^{d})$ ), we can apply Lebesgue’s dominated convergence theorem in the integral in (5.8). This gives (5.7). ∎

Now we are ready to prove Proposition 5.3.

Proof of Proposition 5.3.

First we show that

\displaystyle\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}=\lim_{\varepsilon\to 0}\sup_{z\in S_{\varepsilon}}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}.

(5.11)

Since $S_{\varepsilon_{2}}\subset S_{\varepsilon_{1}}$ if $0<\varepsilon_{1}<\varepsilon_{2}<1$ , the map

(0,\infty)\ni\varepsilon\mapsto\sup_{z\in S_{\varepsilon}}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}

(5.12)

is non-increasing. This implies

\displaystyle\lim_{\varepsilon\to 0}\sup_{z\in S_{\varepsilon}}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}\leq\sup_{z\in\mathbb{C}\setminus[0,\infty)}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}.

To prove the converse inequality, we note that, for all $z_{0}\in\mathbb{C}\setminus[0,\infty)$ we have $z_{0}\in S_{\varepsilon}$ for $\varepsilon>0$ small enough and therefore

\displaystyle\left\||x|^{-1}(-\Delta-z_{0})^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}\leq\lim_{\varepsilon\to 0}\sup_{z\in S_{\varepsilon}}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}.

This implies (5.11).

It remains to establish that, for all $\varepsilon>0$ ,

\displaystyle\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}=\sup_{z\in S_{\varepsilon}}\left\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\right\|_{L^{2}\to L^{2}}.

(5.13)

Let $\varepsilon>0$ and $\varphi,\psi\in C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ . The map

z\mapsto\left\langle\varphi,|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\psi\right\rangle_{L^{2}}

(5.14)

is analytic on $S_{\varepsilon}$ and continuous up to the boundary (the continuity at $z=0$ follows from Lemma 5.5). Moreover since $\varphi,\psi\in C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ , it follows from Lemma 5.4 that

\left|\left\langle\varphi,|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\psi\right\rangle_{L^{2}}\right|\xrightarrow[|z|\to\infty]{}0.

(5.15)

Therefore, by the maximum modulus principle,

\displaystyle\sup_{z\in S_{\varepsilon}}\left|\left\langle\varphi,|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\psi\right\rangle_{L^{2}}\right|=\sup_{|\mathrm{Im}\,z|=\varepsilon\mathrm{Re}\,z}\left|\left\langle\varphi,|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\psi\right\rangle_{L^{2}}\right|.

Since $C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ is dense in $L^{2}$ , this implies (5.13). ∎

5.2 The alternative proof of Theorem 1.2

Now we are in a position to give an alternative multiplier proof of the free upper bound (1.9). The argument is based on the following identity.

Lemma 5.6.

Let $d\geq 2$ . Let $z\in\mathbb{C}\setminus[0,\infty)$ be such that $\mathrm{Re}\,z>0$ . Let $\psi\in\mathscr{S}(\mathbb{R}^{d})$ and let $f:=(-\Delta-z)\psi$ . Then

	$\displaystyle\\|\nabla\psi^{-}\\|^{2}+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\int_{\mathbb{R}^{d}}\|x\|\Big(\big\|\nabla\psi^{-}\big\|^{2}-\frac{d-1}{2}\frac{\|\psi\|^{2}}{\|x\|^{2}}\Big)$
	$\displaystyle=\mathrm{Re}\,\int_{\mathbb{R}^{d}}\Big((d-1)\bar{\psi}f+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\|x\|\bar{\psi}f+2f^{-}x\cdot\overline{\nabla\psi^{-}}\Big),$		(5.16)

where

z_{1}:=\mathrm{Re}\,z,\quad z_{2}:=\mathrm{Im}\,z,\quad\psi^{-}(x):=e^{-i\sqrt{z_{1}}\mathrm{sgn}(z_{2})|x|}\psi(x),

(5.17)

and likewise for $f^{-}$ .

Proof.

The identity (5.16) in a more general context of electromagnetic perturbations of the Laplacian and $\psi$ being merely $H^{1}(\mathbb{R}^{d})$ is due to [7, Lem. 4.1] (see also [9, Eq. (37)] for an abridged proof). We give a proof in the present context to make the paper self-contained. The painful regularisation procedure of [7], originally due to [8], is not needed here because $\psi\in\mathscr{S}(\mathbb{R}^{d})$ .

Consider the equation

-\Delta\psi=z\psi+f\,,

(5.18)

where $\psi,f\in\mathscr{S}(\mathbb{R}^{d})$ . First, multiplying (5.18) by $\bar{\phi}:=x\cdot\nabla\bar{\psi}+\frac{d}{2}\bar{\psi}$ , integrating over $\mathbb{R}^{d}$ , taking twice the real part of the obtained identity and integrating by parts, we obtain

2\int_{\mathbb{R}^{d}}|\nabla\psi|^{2}=-2\,z_{2}\,\mathrm{Im}\,\int_{\mathbb{R}^{d}}\psi\,x\cdot\nabla\bar{\psi}+2\,\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,x\cdot\nabla\bar{\psi}+d\,\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,\bar{\psi}\,.

(5.19)

Noticing that $\bar{\phi}=iA\bar{\psi}$ , where $A:=\frac{1}{2}(x\cdot(-i\nabla)+(-i\nabla)\cdot x)$ is the dilation operator, (5.19) can be considered as a consequence of the virial identity $2(-\Delta)=i[-\Delta,A]$ . Second, multiplying (5.18) by $\bar{\psi}$ , integrating over $\mathbb{R}^{d}$ , taking the real part of the obtained identity and integrating by parts, we get

\int_{\mathbb{R}^{d}}|\nabla\psi|^{2}=z_{1}\int_{\mathbb{R}^{d}}|\psi|^{2}+\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,\bar{\psi}\,.

(5.20)

Third, multiplying (5.18) by $|x|\bar{\psi}$ , integrating over $\mathbb{R}^{d}$ , taking the real and imaginary part of the obtained identity and integrating by parts, we respectively get

\int_{\mathbb{R}^{d}}|x|\,|\nabla\psi|^{2}-\frac{d-1}{2}\,\int_{\mathbb{R}^{d}}\frac{|\psi|^{2}}{|x|}=z_{1}\int_{\mathbb{R}^{d}}|x|\,|\psi|^{2}+\mathrm{Re}\,\int_{\mathbb{R}^{d}}|x|\,\bar{\psi}\,f

(5.21)

and

\mathrm{Im}\,\int_{\mathbb{R}^{d}}\bar{\psi}\,\frac{x}{|x|}\cdot\nabla\psi=z_{2}\int_{\mathbb{R}^{d}}|x|\,|\psi|^{2}+\mathrm{Im}\,\int_{\mathbb{R}^{d}}|x|\,\bar{\psi}\,f\,.

(5.22)

Note that the singular term in (5.21) is well defined due to the weighted Hardy inequality

\int_{\mathbb{R}^{d}}|x|\,|\nabla\psi|^{2}\geq\left(\frac{d-1}{2}\right)^{2}\int_{\mathbb{R}^{d}}\frac{|\psi|^{2}}{|x|}

(5.23)

valid for $d\geq 2$ . Taking the clever sum

\eqref{i1}-\eqref{i2}+\frac{|z_{2}|}{\sqrt{z_{1}}}\,\eqref{i3}-2\,\sqrt{z_{1}}\,\mathop{\mathrm{sgn}}\nolimits(z_{2})\,\eqref{i4}\,,

(5.24)

we arrive at the desired identity (5.16). ∎

We now give the alternative multiplier proof of Theorem 1.2.

Alternative proof of Theorem 1.2.

Estimate (1.9) for $d=3$ follows from the pointwise bound $|G_{z}(x,y)|\leq G_{0}(x,y)$ , where $G_{z}$ denotes the integral kernel of the resolvent $(-\Delta-z)^{-1}$ (see, e.g., [14, Sec. 2]) and the consequence of Hardy’s inequality recalled in (5.2).

Now we prove (1.9) for $d\geq 4$ . By Lemma 5.2 we know that, if $z_{1}=\mathrm{Re}\,z\leq 0$ (and $z\neq 0$ ), then the explicit bound (5.3) holds. Now we consider $z\in\mathbb{C}\setminus[0,\infty)$ with $z_{1}>0$ , $\psi\in\mathscr{S}(\mathbb{R}^{d})$ and we set $f=(-\Delta-z)\psi$ . If $f=0$ then $\psi=0$ . In what follows we assume that $f\neq 0$ . With the notations from Lemma 5.6, we deduce from (5.16) (and the weighted Hardy inequality (5.23) showing that the second term on the left-hand side of (5.16) is non-negative) that

$\displaystyle\\|\nabla\psi^{-}\\|$	$\displaystyle\leq\mathrm{Re}\,\int_{\mathbb{R}^{d}}\Big((d-1)\bar{\psi}f+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\|x\|\bar{\psi}f+2f^{-}x\cdot\overline{\nabla\psi^{-}}\Big)$
	$\displaystyle=\mathrm{Re}\,\int_{\mathbb{R}^{d}}\|x\|f^{-}\Big((d-1)\frac{\overline{\psi^{-}}}{\|x\|}+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\overline{\psi^{-}}+2\frac{x}{\|x\|}\cdot\overline{\nabla\psi^{-}}\Big)$
	$\displaystyle\leq\\|\|x\|f\\|_{L^{2}}\Big\\|(d-1)\frac{\overline{\psi^{-}}}{\|x\|}+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\overline{\psi^{-}}+2\frac{x}{\|x\|}\cdot\overline{\nabla\psi^{-}}\Big\\|_{L^{2}}.$	(5.25)

Next we compute

	$\displaystyle\Big\\|(d-1)\frac{\overline{\psi^{-}}}{\|x\|}+\frac{\|z_{2}\|}{\sqrt{z_{1}}}\overline{\psi^{-}}+2\frac{x}{\|x\|}\cdot\overline{\nabla\psi^{-}}\Big\\|_{L^{2}}^{2}$
	$\displaystyle=(d-1)^{2}\Big\\|\frac{\psi}{\|x\|}\Big\\|_{L^{2}}^{2}+\frac{z_{2}^{2}}{z_{1}}\\|\psi\\|_{L^{2}}^{2}+4\Big\\|\frac{x}{\|x\|}\cdot\overline{\nabla\psi^{-}}\Big\\|_{L^{2}}^{2}$
	$\displaystyle\quad+2(d-1)\frac{\|z_{2}\|}{\sqrt{z_{1}}}\int_{\mathbb{R}^{d}}\frac{\|\psi\|^{2}}{\|x\|}+2(d-1)\int_{\mathbb{R}^{d}}\frac{x}{\|x\|^{2}}\cdot\nabla\|\psi\|^{2}+2\frac{\|z_{2}\|}{\sqrt{z_{1}}}\int_{\mathbb{R}^{d}}\frac{x}{\|x\|}\cdot\nabla\|\psi\|^{2}.$		(5.26)

Integrating by parts in the last two integrals and using $(d-1)^{2}-2(d-1)(d-2)=-(d-1)(d-3)$ , we obtain

\displaystyle\Big\|(d-1)\frac{\overline{\psi^{-}}}{|x|}+\frac{|z_{2}|}{\sqrt{z_{1}}}\overline{\psi^{-}}+2\frac{x}{|x|}\cdot\overline{\nabla\psi^{-}}\Big\|_{L^{2}}^{2}=-(d-1)(d-3)\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}^{2}+\frac{z_{2}^{2}}{z_{1}}\|\psi\|_{L^{2}}^{2}+4\Big\|\frac{x}{|x|}\cdot\overline{\nabla\psi^{-}}\Big\|_{L^{2}}^{2}.

Inserting this into (5.25) yields

\displaystyle\|\nabla\psi^{-}\|_{L^{2}}^{4}

\displaystyle\leq\||x|f\|_{L^{2}}^{2}\Big(-(d-1)(d-3)\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}^{2}+\frac{z_{2}^{2}}{z_{1}}\|\psi\|_{L^{2}}^{2}+4\Big\|\frac{x}{|x|}\cdot\overline{\nabla\psi^{-}}\Big\|_{L^{2}}^{2}\Big).

Now, we have

|z_{2}|\,\|\psi\|_{L^{2}}^{2}=\Big|\mathrm{Im}\,\int_{\mathbb{R}^{d}}\bar{\psi}f\Big|\leq\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\||x|f\|_{L^{2}},

and therefore we arrive at

\displaystyle\|\nabla\psi^{-}\|_{L^{2}}^{4}

\displaystyle\leq\||x|f\|_{L^{2}}^{2}\Big(-(d-1)(d-3)\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}^{2}+\frac{|z_{2}|}{z_{1}}\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\||x|f\|_{L^{2}}+4\|\nabla\psi^{-}\|_{L^{2}}^{2}\Big).

Computing the discriminant of the previous quadratic expression in $\|\nabla\psi^{-}\|_{L^{2}}^{2}$ , it follows that

\displaystyle 16\||x|f\|_{L^{2}}^{4}+4\Big(-(d-1)(d-3)\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}^{2}+\frac{|z_{2}|}{z_{1}}\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\Big)\||x|f\|_{L^{2}}^{2}\geq 0.

Since we have assumed that $f\neq 0$ , this yields

\displaystyle(d-1)(d-3)\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}^{2}-\frac{|z_{2}|}{z_{1}}\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\||x|f\|_{L^{2}}-4\||x|f\|_{L^{2}}^{2}\leq 0.

(5.27)

The discriminant of the previous quadratic expression in $\||x|^{-1}\psi\|_{L^{2}}$ reads

(z_{1}^{-2}z_{2}^{2}+16(d-1)(d-3))\||x|f\|_{L^{2}}^{2},

and therefore (5.27) implies

\displaystyle\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\leq\frac{\frac{|z_{2}|}{z_{1}}+\sqrt{\frac{z_{2}^{2}}{z_{1}^{2}}+16(d-1)(d-3)}}{2(d-1)(d-3)}\||x|f\|_{L^{2}}.

(5.28)

This inequality holds for any $\psi\in\mathscr{S}(\mathbb{R}^{d})$ , with $f=(-\Delta-z)\psi$ . Let now $\varphi\in\mathscr{S}(\mathbb{R}^{d})$ and let $\eta\in C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ . Observe that $(-\Delta-z)^{-1}|x|^{-1}\eta\varphi\in\mathscr{S}(\mathbb{R}^{d})$ . Indeed, with $\mathcal{F}$ the usual Fourier transform, and using that $\eta$ is supported away from $0$ and that $z\notin[0,\infty)$ , we have

\varphi\in\mathscr{S}(\mathbb{R}^{d})\,\Rightarrow\,|x|^{-1}\eta\varphi\in\mathscr{S}(\mathbb{R}^{d})\,\Rightarrow\,\mathcal{F}(|x|^{-1}\eta\varphi)\in\mathscr{S}(\mathbb{R}^{d})\,\Rightarrow\,(|\xi|^{2}-z)^{-1}\mathcal{F}(|x|^{-1}\eta\varphi)\in\mathscr{S}(\mathbb{R}^{d}).

Applying then (5.28) to $\psi=(-\Delta-z)^{-1}|x|^{-1}\eta\varphi$ (and recalling from Lemma 5.1 that $|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}$ extends to a bounded operator on $L^{2}$ ) gives

\displaystyle\big\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\eta\varphi\big\|_{L^{2}}\leq\frac{\frac{|z_{2}|}{z_{1}}+\sqrt{\frac{z_{2}^{2}}{z_{1}^{2}}+16(d-1)(d-3)}}{2(d-1)(d-3)}\|\eta\varphi\|_{L^{2}},

(5.29)

for any $\varphi\in\mathscr{S}(\mathbb{R}^{d})$ and $\eta\in C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ . Since the set $\{\varphi\in\mathscr{S}(\mathbb{R}^{d})\,|\,\exists\eta\in C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\}),\,\eta\varphi=\varphi\}$ is dense in $L^{2}(\mathbb{R}^{d})$ (as it contains $C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ , inequality (5.29) leads to

\displaystyle\big\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\big\|_{L^{2}\to L^{2}}\leq\frac{\frac{|z_{2}|}{z_{1}}+\sqrt{\frac{z_{2}^{2}}{z_{1}^{2}}+16(d-1)(d-3)}}{2(d-1)(d-3)}.

(5.30)

Applying now Proposition 5.3, we deduce that

	$\displaystyle\sup_{z_{1}>0\,,z_{2}\neq 0}\big\\|\|x\|^{-1}(-\Delta-z)^{-1}\|x\|^{-1}\big\\|_{L^{2}\to L^{2}}$	$\displaystyle\leq\lim_{\varepsilon\to 0}\sup_{z_{2}=\varepsilon z_{1}}\frac{\varepsilon+\sqrt{\varepsilon^{2}+16(d-1)(d-3)}}{2(d-1)(d-3)}$
		$\displaystyle=\frac{2}{\sqrt{(d-1)(d-3)}}.$		(5.31)

This concludes the alternative proof of (1.9). ∎

Remark 5.7.

Similarly as in Remark 3.8, by scaling invariance, the norm $\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\|_{L^{2}\to L^{2}}$ is constant along the ray $e^{i\arg(z)}(0,\infty)$ . More precisely, given any positive number $\tau$ , and defining the unitary transform $U_{\tau}$ on $L^{2}$ by

(U_{\tau}\psi)(x):=\tau^{d/2}\,\psi(\tau x)\,,

we have

	$\displaystyle\left\\|\|x\|^{-1}(-\Delta-z)^{-1}\|x\|^{-1}\right\\|_{L^{2}\to L^{2}}$	$\displaystyle=\left\\|U_{\tau}\|x\|^{-1}(-\Delta-z)^{-1}\|x\|^{-1}U_{\tau}^{-1}\right\\|_{L^{2}\to L^{2}}$
		$\displaystyle=\left\\|\|x\|^{-1}(-\Delta-\tau^{2}z)^{-1}\|x\|^{-1}\right\\|_{L^{2}\to L^{2}}\,,$

for any $\tau>0$ . We did not use this property in our proof of Theorem 1.2.

Remark 5.8.

The map

\mathbb{C}\setminus[0,\infty)\ni z\mapsto|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\in\mathcal{B}(L^{2})

is analytic and, as proved in [20], the limits

|x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}:=\lim_{\varepsilon\to 0^{+}}|x|^{-1}(-\Delta-\lambda\pm i\varepsilon)^{-1}|x|^{-1}

exist in the operator norm, for any $\lambda>0$ . Moreover, the maps

(0,\infty)\ni\lambda\mapsto|x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}\in\mathcal{B}(L^{2})

are continuous. However, $|x|^{-1}(-\Delta-z)^{-1}|x|^{-1}$ does not converge to $|x|^{-1}(-\Delta)^{-1}|x|^{-1}$ in $\mathcal{B}(L^{2})$ when $z\in\mathbb{C}\setminus[0,\infty)$ , $z\to 0$ . Indeed, if $z$ is any fixed complex number in $\mathbb{C}\setminus[0,\infty)$ , then by Remark 5.7,

		$\displaystyle\left\\|\|x\|^{-1}(-\Delta-\tau^{2}z)^{-1}\|x\|^{-1}-\|x\|^{-1}(-\Delta)^{-1}\|x\|^{-1}\right\\|_{L^{2}\to L^{2}}$
		$\displaystyle=\left\\|\|x\|^{-1}(-\Delta-z)^{-1}\|x\|^{-1}-\|x\|^{-1}(-\Delta)^{-1}\|x\|^{-1}\right\\|_{L^{2}\to L^{2}}.$

Hence if the left-hand side tends to $0$ as $\tau\to 0$ , this implies that the right-hand side vanishes for any $z\in\mathbb{C}\setminus[0,\infty)$ , which is certainly not true. Similarly, $|x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}$ does not converge to $|x|^{-1}(-\Delta)^{-1}|x|^{-1}$ in $\mathcal{B}(L^{2})$ .

Applying Theorem 1.2 to a vector of the form $|x|(-\Delta-z)\psi$ with $\psi$ regular enough (e.g., $\psi\in C_{0}^{\infty}(\mathbb{R}^{d})$ ) and then letting $\mathrm{Im}\,z\to 0$ , we obtain, for instance in the case $d\geq 4$ , that

\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}}\leq\frac{2}{\sqrt{(d-1)(d-3)}}\,\big\||x|(-\Delta-\lambda)\psi\big\|_{L^{2}},\quad\lambda>0,\quad\psi\in C_{0}^{\infty}(\mathbb{R}^{d}).

(5.32)

For $\lambda\leq 0$ , Hardy’s inequality gives a better estimate. The following proposition shows that even for positive $\lambda$ , the constant on the right-hand side of (5.32) can be improved.

Proposition 5.9.

Let $d\geq 3$ . For every $\lambda>0$ and $\psi\in C_{0}^{\infty}(\mathbb{R}^{d})$ , one has

\Big\|\frac{\psi}{|x|}\Big\|_{L^{2}(\mathbb{R}^{d})}\leq C_{d}\,\big\||x|(-\Delta-\lambda)\psi\big\|_{L^{2}(\mathbb{R}^{d})}\qquad\mbox{with}\qquad C_{d}:=\begin{cases}4&\mbox{if}\quad d=3\,,\\ 1&\mbox{if}\quad d=4\,,\\ \displaystyle\frac{1}{\sqrt{d(d-4)}}&\mbox{if}\quad d\geq 5\,.\end{cases}

(5.33)

The proof of Proposition 5.9 is based on the following well-known virial identity (see, e.g., [7, Lem. 4.6]).

Lemma 5.10.

Let $\lambda\in\mathbb{R}$ , $\psi\in C_{0}^{\infty}(\mathbb{R}^{d})$ and $f:=(-\Delta-\lambda)\psi$ . Then

2\int_{\mathbb{R}^{d}}|\nabla\psi|^{2}=d\,\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,\overline{\psi}+2\,\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,x\cdot\overline{\nabla\psi}\,.

(5.34)

Proof.

(5.34) is a special case of the previously established (5.19) (for whose validity the hypothesis $z\in\mathbb{C}\setminus[0,\infty)$ of Lemma 5.6 is not needed). ∎

Now we are ready to prove Proposition 5.9.

Proof of Proposition 5.9.

Using the notations of Lemma 5.10, and starting with the virial identity (5.34), we estimate the right-hand side as follows:

	$\displaystyle\left\|d\,\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,\overline{\psi}+2\,\mathrm{Re}\,\int_{\mathbb{R}^{d}}f\,x\cdot\overline{\nabla\psi}\right\|$	$\displaystyle=\left\|\mathrm{Re}\,\int_{\mathbb{R}^{d}}\|x\|\,f\left(d\,\frac{\overline{\psi}}{\|x\|}+2\,\frac{x}{\|x\|}\cdot\overline{\nabla\psi}\right)\right\|$
		$\displaystyle\leq\big\\|\|x\|f\big\\|_{L^{2}(\mathbb{R}^{d})}\,\left\\|d\,\frac{\overline{\psi}}{\|x\|}+2\,\frac{x}{\|x\|}\cdot\overline{\nabla\psi}\right\\|_{L^{2}(\mathbb{R}^{d})}\,.$

Integrating by parts, we have

	$\displaystyle\left\\|d\,\frac{\overline{\psi}}{\|x\|}+2\,\frac{x}{\|x\|}\cdot\overline{\nabla\psi}\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}$	$\displaystyle=d^{2}\,\left\\|\frac{\psi}{\|x\|}\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}+4\,\left\\|\nabla\psi\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}+2\,d\int_{\mathbb{R}^{d}}\frac{x}{\|x\|^{2}}\cdot\nabla\|\psi\|^{2}$
		$\displaystyle=d^{2}\,\left\\|\frac{\psi}{\|x\|}\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}+4\,\left\\|\nabla\psi\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}-2\,d\,(d-2)\left\\|\frac{\psi}{\|x\|}\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}$
		$\displaystyle=4\,\left\\|\nabla\psi\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}-d\,(d-4)\left\\|\frac{\psi}{\|x\|}\right\\|_{L^{2}(\mathbb{R}^{d})}^{2}\,.$

Consequently, (5.34) implies

4\,\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{d})}^{4}\leq\big\||x|f\big\|_{L^{2}(\mathbb{R}^{d})}^{2}\left(4\,\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{d})}^{2}-d\,(d-4)\left\|\frac{\psi}{|x|}\right\|_{L^{2}(\mathbb{R}^{d})}^{2}\right)\,.

(5.35)

$d=3$ Estimate (5.33) follows from the pointwise bound $|G_{z}(x,y)|\leq G_{0}(x,y)$ , where $G_{z}$ denotes the integral kernel of the resolvent $(-\Delta-z)^{-1}$ (see, e.g., [14, Sec. 2]) and the consequence of Hardy’s inequality recalled in (5.2). Anyway, let us show how to deduce it from (5.35) with the help of the Hardy inequality (5.1):

4\,\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{3})}^{4}\leq\big\||x|f\big\|_{L^{2}(\mathbb{R}^{3})}^{2}\left(4\,\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{3})}^{2}+3\left\|\frac{\psi}{|x|}\right\|_{L^{2}(\mathbb{R}^{3})}^{2}\right)\leq 16\,\big\||x|f\big\|_{L^{2}(\mathbb{R}^{3})}^{2}\,\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{3})}^{2}\,.

Dividing by $\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{3})}^{2}$ and applying the Hardy inequality (5.1) once more, we conclude with

\left\|\frac{\psi}{|x|}\right\|_{L^{2}(\mathbb{R}^{3})}^{2}\leq 16\,\big\||x|f\big\|_{L^{2}(\mathbb{R}^{3})}^{2}\,,

which proves (5.33).

$d=4$ In four dimensions, (5.35) and the Hardy inequality (5.1) directly imply

\left\|\frac{\psi}{|x|}\right\|_{L^{2}(\mathbb{R}^{4})}^{2}\leq\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{4})}^{2}\leq\big\||x|f\big\|_{L^{2}(\mathbb{R}^{4})}^{2}\,,

which gives (5.33).

$d\geq 5$ In these higher dimensions, we regard (5.35) as the quadratic inequality

4\,A^{2}-4\,\big\||x|f\big\|_{L^{2}(\mathbb{R}^{3})}^{2}\,A-d\,(d-4)\,\big\||x|f\big\|_{L^{2}(\mathbb{R}^{3})}^{2}\left\|\frac{\psi}{|x|}\right\|_{L^{2}(\mathbb{R}^{d})}^{2}\leq 0

for the unknown $A:=\left\|\nabla\psi\right\|_{L^{2}(\mathbb{R}^{d})}^{2}$ Since the discriminant reads

D:=16\,\big\||x|f\big\|_{L^{2}(\mathbb{R}^{d})}^{2}\left(\big\||x|f\big\|_{L^{2}(\mathbb{R}^{d})}^{2}-d\,(d-4)\left\|\frac{\psi}{|x|}\right\|_{L^{2}(\mathbb{R}^{d})}^{2}\right)\,,

estimate (5.33) follows by requiring $D\geq 0$ . ∎

Remark 5.11.

Using in particular the maximum modulus (or the Phragmén–Lindelöf) principle, one obtains

\sup_{z\in\mathbb{C}\setminus[0,\infty)}\||x|^{-1}(-\Delta-z)^{-1}|x|^{-1}\|_{L^{2}\to L^{2}}=\sup_{\lambda>0}\||x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}\|_{L^{2}\to L^{2}}.

See the proof of Theorem 1.2 above for a similar argument. It is therefore tempting to think that Proposition 5.9 should imply that $\||x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}\|_{L^{2}\to L^{2}}\leq C_{d}$ , with $C_{d}$ the constant given by Proposition 5.9. This is however wrong since, as mentioned in the introduction, we have the lower bound $\||x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}\|_{L^{2}\to L^{2}}\geq\frac{\pi}{2(d-2)}$ and $\frac{\pi}{2(d-2)}>C_{d}$ for $d\geq 5$ . The fact that Proposition 5.9 does not imply a bound on $\||x|^{-1}(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}\|_{L^{2}\to L^{2}}$ is due to the lack of integrability of vectors of the form $(-\Delta-\lambda\pm i0)^{-1}|x|^{-1}\varphi$ , even if $\varphi$ belongs to $C_{0}^{\infty}(\mathbb{R}^{d}\setminus\{0\})$ .

Acknowledgements

The first author was partially supported by the Progetto Ricerca Scientifica 2024 “Wave dynamics in heterogeneous media” of Sapienza University, and by the Gruppo Nazionale per l’Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) The last author was supported by the grant no. 26-21940S of the Czech Science Foundation.

References

[1] J. A. Barceló, L. Vega, and M. Zubeldia, The forward problem for the electromagnetic Helmholtz equation with critical singularities, Adv. Math. 240 (2013), 636–671.
[2] J.-M. Bouclet and H. Mizutani, Uniform resolvent and Strichartz estimates for Schrödinger equations with critical singularities, Trans. Amer. Math. Soc. 370 (2018), no. 10, 7293–7333.
[3] T. A. Bui, P. D’Ancona, X. T. Duong, J. Li, and F. K. Ly, Weighted estimates for powers and smoothing estimates of Schrödinger operators with inverse square potentials, J. Differential Equations 262 (2017), no. 3, 2771–2807.
[4] N. Burq, F. Planchon, J. G. Stalker, and A. S. Tahvildar-Zadeh, Strichartz estimates for the wave and Schrödinger equations with the inverse-square potential, J. Funct. Anal. 203 (2003), 519–549.
[5] N. Burq, F. Planchon, J. G. Stalker, and A. S. Tahvildar-Zadeh, Strichartz estimates for the wave and Schrödinger equations with potentials of critical decay, Indiana Univ. Math. J. 53 (2004), no. 6, 1665–1680.
[6] C. Cazacu and D. Krejčiřík, The Hardy inequality and the heat equation with magnetic field in any dimension, Comm. Partial Differential Equations 41 (2016), 1056–1088.
[7] L. Cossetti, L. Fanelli, and D. Krejčiřík, Absence of eigenvalues of Dirac and Pauli Hamiltonians via the method of multipliers, Comm. Math. Phys. 379 (2020), 633–691.
[8] L. Cossetti and D. Krejčiřík, Absence of eigenvalues of non-selfadjoint Robin Laplacians on the half-space, Proc. London. Math. Soc. 121 (2020), 584–616.
[9] , The virial theorem and the method of multipliers in spectral theory, Mathematical Models for Interacting Dynamics on Networks, Trends in Mathematics, Birkhäuser Cham, 2025.
[10] P. D’Ancona, Kato smoothing and Strichartz estimates for wave equations with magnetic potentials, Comm. Math. Phys. 335 (2015), 1–16.
[11] J. Dereziński and S. Richard, On Schrödinger operators with inverse square potentials on the half–line, Ann. Henri Poincaré 18 (2017), 869–928.
[12] S. Dyatlov and M. Zworski, Mathematical theory of scattering resonances, Amer. Math. Soc., 2019.
[13] L. Fanelli, D. Krejčiřík, and L. Vega, Absence of eigenvalues of two dimensional magnetic Schrödinger operators, J. Funct. Anal. 275 (2018), 2453–2472.
[14] , Spectral stability of Schrödinger operators with subordinated complex potentials, J. Spectr. Theory 8 (2018), 575–604.
[15] L. Fanelli, J. Zhang, and J. Zheng, Uniform resolvent estimates for critical magnetic Schrödinger operators in 2D, Int. Math. Res. Not. IMRN 2023, no. 20, 17656–17703.
[16] R. L. Frank, Eigenvalue bounds for Schrödinger operators with complex potentials, Bull. Lond. Math. Soc. 43 (2011), 745–750.
[17] I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series, and products., Amsterdam: Elsevier/Academic Press, 8th edition, 2015.
[18] M. Hansmann and D. Krejčiřík, The abstract Birman-Schwinger principle and spectral stability, J. Anal. Math. 148 (2022), 361–398.
[19] T. Kato, Wave operators and similarity for some non-selfadjoint operators, Math. Ann. 162 (1966), 258–279.
[20] T. Kato and K. Yajima, Some examples of smooth operators and the associated smoothing effect, Rev. Math. Phys. 1 (1989), 481–496.
[21] C. E. Kenig, A. Ruiz, and C. D. Sogge, Uniform sobolev inequalities and unique continuation for second order constant coefficient differential operators, Duke Math. J. 55 (1987), 329–347.
[22] D. Krejčiřík, Complex magnetic fields: An improved Hardy-Laptev-Weidl inequality and quasi-selfadjointness, SIAM J. Math. Anal. 51 (2019), 790–807.
[23] A. Laptev and T. Weidl, Hardy inequalities for magnetic Dirichlet forms, Oper. Theory Adv. Appl. 108 (1999), 299–305.
[24] H. Mizutani, J. Zhang, and J. Zheng, Uniform resolvent estimates for Schrödinger operator with an inverse square potential, J. Funct. Anal. 278 (2020), no. 4, 108350.
[25] K. Mochizuki, Resolvent estimates for magnetic Schrödinger operators and their applications to related evolution equations, Rend. Istit. Mat. Univ. Trieste 42 (2010), 143–164.
[26] K. Mochizuki, Uniform resolvent estimates for magnetic Schrödinger operators and smoothing effects for related evolution equations Publ. RIMS Kyoto Univ. 46 (2010), 741–754.
[27] K. Pankrashkin and S. Richard, Spectral and scattering theory for the Aharonov-Bohm operators, Rev. Math. Phys. 23 (2011), 53–81.
[28] A. P. Prudnikov, Yu. A. Brychkov O. I. Marichev, Integrals and Series, Vol. 2, Special Functions, Gordon & Breach Sci. Publ., New York, 1986.
[29] M. Reed, B. Simon, Methods of modern mathematical physics. IV Functional Analysis, Academic Press, New York-London, 1980.
[30] G. Szegö, Orthogobal Polynomials, 4th edition, Amer. Math. Soc., Colloquium Publications, 1975.
[31] B. Simon, Best constants in some operator smoothness estimates, J. Funct. Anal. 107 (1992), no. 3, 66–71.

	$\displaystyle\Big\|\big\\|A_{\nu}^{\pm}\big\\|_{L^{2}\to L^{2}}-\frac{1}{\nu}\Big\|$	$\displaystyle=\Big\|\big\\|A_{\nu}^{\pm}\big\\|_{L^{2}\to L^{2}}-\big\\|\nu^{-1}C^{\pm}\big\\|_{L^{2}\to L^{2}}\Big\|$
		$\displaystyle\leq\big\\|A_{\nu}^{\pm}-B_{\nu}^{\pm}+B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\big\\|_{L^{2}\to L^{2}}$
		$\displaystyle\leq\big\\|r^{-1}h_{\nu}^{-1}\mathbf{1}_{h_{\nu}\geq 1}r^{-1}\big\\|_{L^{2}\to L^{2}}+\big\\|B_{\nu}^{\pm}-\nu^{-1}C^{\pm}\big\\|_{\mathrm{HS}},$

	$\displaystyle\Big\|\int_{0}^{\infty}\Big(\frac{\|z_{2}\|}{\sqrt{z_{1}}}rf^{-}\overline{\psi^{-}}-\frac{\|z_{2}\|}{\sqrt{z_{1}}}(\nu^{2}-\mbox{$\frac{1}{4}$})\frac{\|\psi^{-}\|^{2}}{r}\Big)\Big\|$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\\|rf\\|_{L^{2}}^{\frac{3}{2}}\Big\\|\frac{\psi}{r}\Big\\|^{\frac{1}{2}}_{L^{2}}+\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big\\|\frac{\psi}{r}\Big\\|^{\frac{3}{2}}_{L^{2}}\\|rf\\|_{L^{2}}^{\frac{1}{2}}\Big)$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\mbox{$\frac{3}{4}$}\\|rf\\|^{2}_{L^{2}}+\mbox{$\frac{1}{4}$}\Big\\|\frac{\psi}{r}\Big\\|^{2}_{L^{2}}+\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big(\mbox{$\frac{3}{4}$}\Big\\|\frac{\psi}{r}\Big\\|^{2}_{L^{2}}+\frac{1}{4}\\|rf\\|^{2}_{L^{2}}\Big)\Big)$
	$\displaystyle=\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\Big(\mbox{$\frac{3}{4}$}+\mbox{$\frac{1}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\\|rf\\|^{2}_{L^{2}}+\Big(\mbox{$\frac{1}{4}$}+\mbox{$\frac{3}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\Big\\|\frac{\psi}{r}\Big\\|^{2}_{L^{2}}\Big)$
	$\displaystyle\leq\sqrt{\frac{\|z_{2}\|}{z_{1}}}\Big(\Big(\mbox{$\frac{3}{4}$}+\mbox{$\frac{1}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\\|rf\\|^{2}_{L^{2}}+\nu^{-2}\Big(\mbox{$\frac{1}{4}$}+\mbox{$\frac{3}{4}$}\Big\|\nu^{2}-\mbox{$\frac{1}{4}$}\Big\|\Big)\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|^{2}_{L^{2}}\Big),$		(3.11)

	$\displaystyle\Big(1-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{1,\nu}\Big)^{2}\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|_{L^{2}}^{4}$
	$\displaystyle\leq\Big(2\\|rf\\|_{L^{2}}\big\\|\partial_{r}\psi^{-}\big\\|_{L^{2}}-\sqrt{\frac{\|z_{2}\|}{z_{1}}}C_{2,\nu}\\|rf\\|^{2}_{L^{2}}\Big)^{2}$
	$\displaystyle\leq 4\\|rf\\|^{2}_{L^{2}}\big\\|\partial_{r}\psi^{-}\big\\|^{2}_{L^{2}}+\frac{\|z_{2}\|}{z_{1}}C_{2,\nu}^{2}\\|rf\\|^{4}_{L^{2}}$
	$\displaystyle=4\\|rf\\|^{2}_{L^{2}}\big\\|h_{\nu}^{\frac{1}{2}}\psi^{-}\big\\|^{2}_{L^{2}}-4\Big(\nu^{2}-\mbox{$\frac{1}{4}$}\Big)\Big\\|\frac{\psi}{r}\Big\\|_{L^{2}}^{2}\\|rf\\|_{L^{2}}^{2}+\frac{\|z_{2}\|}{z_{1}}C_{2,\nu}^{2}\\|rf\\|^{4}_{L^{2}}.$

	$\displaystyle\frac{\\|r^{-1}u_{\sigma}\\|_{L^{2}}}{\\|r(h_{\nu}-z)u_{\sigma}\\|_{L^{2}}}$	$\displaystyle\geq\frac{\\|r^{-1}u_{\sigma}\\|_{L^{2}}}{\nu^{2}\\|r^{-1}u_{\sigma}\\|_{L^{2}}+\\|A_{\sigma}(z)\\|_{L^{2}}}$
		$\displaystyle=\frac{1}{\nu^{2}+\\|A_{\sigma}(z)\\|_{L^{2}}\\|r^{-1}u_{\sigma}\\|_{L^{2}}^{-1}}$
		$\displaystyle\geq\frac{1}{\nu^{2}-C_{\chi,\eta}(1+\|z\|)(\ln\sigma)^{-1}}.$

	$\displaystyle\|k\|\,\left\|\left\langle\varphi,(-\Delta-z)^{-1}\psi\right\rangle_{L^{2}}\right\|$	$\displaystyle=\|k\|\,\left\|\left\langle\hat{\varphi},(\|\cdot\|^{2}-k^{2})^{-1}\hat{\psi}\right\rangle\right\|$
		$\displaystyle=\left\|\int_{\mathbb{R}^{d}}\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\int_{0}^{\infty}\frac{\sin(\|\xi\|t)}{\|\xi\|}\,k\,e^{ikt}\,\mathrm{d}t\,\mathrm{d}\xi\right\|$
		$\displaystyle=\left\|\int_{0}^{\infty}\int_{\mathbb{R}^{d}}\overline{\hat{\varphi}(\xi)}\,\hat{\psi}(\xi)\,\cos(\|\xi\|t)\,e^{ikt}\,\mathrm{d}\xi\,\mathrm{d}t\right\|\,,$

Quantitative uniform resolvent estimates

Abstract

1 Introduction

1.1 Reduction to one-dimensional Bessel operators

1.2 Multidimensional examples beyond the Laplacian

1.3 Main results for the one-dimensional Bessel operators

Theorem 1.1 (Bessel operators).

1.4 Main results for the multidimensional examples

Theorem 1.2 (Laplacian).

Proposition 1.3 (Dimensional monotonicity).

Theorem 1.4 (Exact imaginary parts).

Theorem 1.5 (Inverse square potentials).

Theorem 1.6 (Aharonov–Bohm potentials).

1.5 Applications

Kato smoothing and evolution control

Spectral enclosures and stability

1.6 Structure of the paper

2 The matrix computation for one Bessel channel

Theorem 2.1.

2.1 A lower bound to the norm of Aν±A_{\nu}^{\pm}

Proposition 2.2.

Proof.

Proposition 2.3.

Proof.

Theorem 2.4.

Proof.

2.2 Asymptotics of the norm of Aν±A_{\nu}^{\pm} as ν→∞\nu\to\infty

Theorem 2.5.

Remark 2.6.

Remark 2.7.

Proposition 2.8.

Remark 2.9.

Lemma 2.10.

Proof.

Lemma 2.11.

Proof.

Proof of Theorem 2.5.

2.3 The Hankel transform

2.4 The orthonormal basis at low energies

Lemma 2.12.

Proof.

Lemma 2.13.

Remark 2.14.

Proof.

Lemma 2.15.

Proof.

Lemma 2.16.

Proof.

2.5 The orthonormal basis at high energies

Lemma 2.17.

Proof.

Lemma 2.18.

Proof.

Lemma 2.19.

Proof.

Lemma 2.20.

Proof.

Lemma 2.21.

Proof.

2.6 Proof of Theorem 2.1 and Proposition 2.8

Proof of Theorem 2.1.

Lemma 2.22.

Proof.

Lemma 2.23.

Proof.

Lemma 2.24.

Proof.

Proof of Proposition 2.8.

3 One-channel estimates

Theorem 3.1.

3.1 Basic facts

3.2 Preliminaries

Lemma 3.2.

Proof.

Lemma 3.3.

Proof.

Lemma 3.4.

Proof.

Proposition 3.5.

Lemma 3.6.

2.1 A lower bound to the norm of $A_{\nu}^{\pm}$

2.2 Asymptotics of the norm of $A_{\nu}^{\pm}$ as $\nu\to\infty$