On real resonances for three-dimensional Schrödinger operators with point interactions

1.   Introduction and background: Dispersive properties of the point-singular perturbed Schrödinger equation

A typical obstacle to the dispersive and scattering properties of the time evolution group associated with the Schrödinger equation

in the unknown $u\equiv u(t, x)$, where $t\in\mathbb{R}$, $x\in\mathbb{R}^d$ ($d\in\mathbb{N}$), and $V:\mathbb{R}^d\to\mathbb{R}$ is a given measurable potential, is the existence of non-trivial solutions to

for some $\mu\in\mathbb{R}$.

In those cases, relevant in a variety of contexts, where $V$ is sufficiently localised ('with short range') and/or is a suitably small perturbation of the Laplacian, non-trivial $L^2(\mathbb{R}^d)$-solutions to (1.2) are interpreted as bound states of the associated Schrödinger operator, and if $\mu > 0$ one refers to it as an eigenvalue embedded in the continuum. Solutions to (1.2) in weaker $L^2$-weighted spaces are generally known instead as resonances (a notion that we shall explicitly define in due time for the purposes of the present analysis), and they too affect the dispersive and scattering behaviour of the propagator defined by (1.1).

When $V\in L^2_{\mathrm{loc}}(\mathbb{R}^d)$, and in (1.2) $u\in H^2_{\mathrm{loc}}(\mathbb{R}^d)$, it was first proved by Kato ^[30] that positive eigenvalues are absent, and by Agmon ^[1] and by Alshom and Smith ^[10] that positive resonances are absent too. For rougher (non-$L^2_{\mathrm{loc}}$) potentials, positive eigenvalues were excluded by Ionescu and Jerison ^[27] and by Koch and Tataru ^[32] by means of suitable Carleman-type estimates which imply, owing to a unique continuation principle ^[29,31], that the corresponding eigenfunctions must be compactly supported and hence vanish. Absence of positive resonances whose associated resonant state $u$ (solution to (1.2)) satisfies appropriate radiation conditions at infinity, was proved by Georgiev and Visciglia ^[20] for $L^{d/2}_{\mathrm{loc}}$-potentials decaying as $|x|^{-(1+\varepsilon)}$ or faster.

A closely related and equally challenging context, which this work is part of, is the counterpart problem of existence or non-existence of spectral obstructions (eigenvalues or resonances) when the potential $V$ in (1.1) and (1.2) is formally replaced by a finite number of delta-like bumps localised at certain given points in space. This corresponds to a well-established rigorous construction of point-like perturbations of the free Laplacian usually referred to as 'Schrödinger operator with point interaction' (we refer to the monograph ^[4] as the standard reference in this field).

We recall that among the various equivalent, yet conceptually alternative ways of defining on $L^2(\mathbb{R}^d)$ the formal operator

obtained by adding to the free Laplacian $N$ singular perturbations centred at the points $y_1, \dots, y_N\in\mathbb{R}^d$ and of magnitude, respectively, $\nu_1, \dots, \nu_N\in\mathbb{R}$, one is to obtain (1.3) as the limit of Schrödinger operators with actual potentials $V^{(j)}_\varepsilon(x-y_j)$ each of which, as $\varepsilon\downarrow 0$, spikes up to a delta-like profile, the support shrinking to the point $\{y_j\}$, and another way is to define (1.3) as a self-adjoint extension of the restriction of $-\Delta$ to smooth functions supported away from the $y_j$'s. Either approach reproduces the free Laplacian unless when $d = 1, 2, 3$, in which case one obtains a non-trivial perturbation of $-\Delta$.

In this work we shall focus on $d = 3$ spatial dimensions and therefore fix a collection

of $N$ distinct points where the perturbation is supported at. It turns out [4,Section Ⅱ.1.1] that the operator

is densely defined, real symmetric, and non-negative on $L^2(\mathbb{R}^d)$, and admits a $N^2$-real-parameter family of self-adjoint extensions, each of which acts as the free negative Laplacian on functions that are supported away from the interaction centres. The most relevant extensions constitute the $N$-parameter sub-family

of so-called 'local' extensions, that is, the rigorous version of (1.3): for them the functions $u$ in the domain of self-adjointness are only qualified by certain local boundary conditions at each singularity centre, with no pairing between distinct centres, which take the explicit form

Physically, each $\alpha_j$ is proportional to the inverse scattering length of the interaction supported at $y_j$. In particular, if for some $j\in\{1, \ldots, N\}$ one has $\alpha_j = \infty$, then no actual interaction is present at the point $y_j$, and in practice things are as if one discards it. When $\alpha = \infty$, one recovers the the Friedrichs extension of $-\Delta\upharpoonright C_0^\infty(\mathbb{R}^3 \!\setminus\! Y)$, namely the self-adjoint negative Laplacian with domain $H^2(\mathbb R^3)$. We may henceforth assume, without loss of generality, that $\alpha$ runs over $\mathbb R^N$.

The perturbations $-\Delta_{\alpha, Y}$ of $-\Delta$ have a long history of investigation and in Section 2 we shall list a number of properties that rigorously qualify them and are relevant for our subsequent analysis. In the mathematical literature they were introduced and characterised for the case $N = 1$ by Berezin and Faddeev ^[12], Albeverio, Høegh-Krohn, and Streit ^[6], Nelson ^[35], Albeverio, Fenstad, and Høegh-Krohn ^[3], and Albeverio and Høegh-Krohn ^[5]. For generic $N\geqslant 1$ centres, $-\Delta_{\alpha, Y}$ was rigorously studied first by Albeverio, Fenstad, and Høegh-Krohn ^[3], and subsequently characterised by Zorbas ^[45], Grossmann, Høegh-Krohn, and Mebkhout ^[24,25], Dąbrowski and Grosse ^[14], and more recently by Arlinskiĭ and Tsekanovskiĭ ^[11], and by Goloshchapova, Malamud, and Zastavnyi ^[22,23].

The analysis of the dispersive and scattering properties of the Schrödinger propagator $e^{ \mathrm{i} t\Delta_{\alpha, Y}}$, $t\in\mathbb{R}$, has been an active subject as well. A class of $L^p\to L^q$ dispersive estimates were established by D'Ancona, Pierfelice, and Teta ^[15] (in weighted form), and by Iandoli and Scandone ^[26] (removing the weights used in ^[15] in the largest regime possible of the $(p, q)$-indices). The $L^p$-boundedness of the wave operators for the pair $(-\Delta_{\alpha, Y}, -\Delta)$ in the regime $p\in(1, 3)$ (from which dispersive and Strichartz estimates can be derived by intertwining $-\Delta$ and $-\Delta_{\alpha, Y}$), as well as the $L^p$-unboundedness of the wave operators when $p = 1$ or $p\in[3, +\infty]$, was proved by Dell'Antonio et al. ^[16] (with counterpart results by Duchêne, Marzuola, and Weinstein ^[17] in $d = 1$ and Cornean, Michelangeli, and Yajima ^[13] in $d = 2$ dimensions).

In analogy with the ordinary Schrödinger equation (1.1), also the dispersive features of the singular point-perturbed Schrödinger equation

strictly depend on the possible presence of eigenvalues or resonances for $-\Delta_{\alpha, Y}$, and indeed in the above-mentioned works ^[15,16,26] special assumptions on the choice of $\alpha$ and $Y$ are often made so as to ensure that no spectral obstruction occurs.

In fact (see Theorem 2.2 below for the complete summary and references), the spectrum $\sigma(-\Delta_{\alpha, Y})$ only consists of an absolutely continuous component $[0, +\infty)$ which is also the whole essential spectrum, plus possibly a number of non-positive eigenvalues. Thus, as usual, for the purposes of the dispersive analysis, one considers $P_{\mathrm{ac}}\, e^{ \mathrm{i} t\Delta_{\alpha, Y}}$, namely the action of the singular Schrödinger propagator on the sole absolutely continuous subspace of $L^2(\mathbb{R}^3)$, and additionally one has to decide whether possible resonances are present.

When $N = 1$ the picture is completely controlled: $-\Delta_{\alpha, Y}$ has only one negative eigenvalue if $\alpha < 0$, and has only a resonance, at zero, if $\alpha = 0$; correspondingly the integral kernel of the propagator $e^{ \mathrm{i} t\Delta_{\alpha, Y}}$ is explicitly known, as found by Scarlatti and Teta ^[39] and Albeverio, Brzeźniak, and Dąbrowski ^[2], from which $L^p\to L^q$ dispersive estimates are derived directly, as found in ^[15]. In certain regimes of $p, q$ slower decay estimates do emerge in the resonant case $\alpha = 0$, as opposed to the non-resonant one.

For generic $N$ perturbation centres, it is again well understood (see Theorem 2.2 below) that at most $N$ non-positive eigenvalues can add up to the absolutely continuous spectrum $[0, +\infty)$ of $-\Delta_{\alpha, Y}$. In particular, as discussed by one of us in [38,Sect. 3], a zero-energy eigenvalue may occur (see also [4,page 485]).

The study of resonances for generic $N$ has been attracting a considerable amount of attention. As explained in Section 2, it is known since the already mentioned work ^[25] by Grossmann, Høegh-Krohn, and Mebkhout (see also [4,Sect. Ⅱ.1.1]), that resonances and eigenvalues $z^2$ of $-\Delta_{\alpha, Y}$ are detected, on an equal footing, by the singularity of an auxiliary $N\times N$ square matrix $\Gamma_{\alpha, Y}(z)$ depending on $z\in\mathbb{C}$. Real negative resonances (thus $z = \mathrm{i}\lambda$ with $\lambda > 0$) are excluded by the arguments of ^[25]. A zero resonance may occur, and one of us ^[38] qualified this possibility in terms of a convenient low-energy resolvent expansion which is at the basis of our definition 2.5 below. Complex resonances ($\mathfrak{Im}z < 0)$ have been investigated by Albeverio and Karabash ^[7,8,9] and Lipovský and Lotoreichik ^[34], using techniques on the localisation of zeroes of exponential polynomials, and turn out to lie mostly within certain logarithmic strips in the complex $z$-plane. Real positive resonances (thus, $z\in\mathbb{R}\setminus\{0\}$) have been recently excluded by Galtbayar and Yajima ^[19], and implicitly also by Sjöstrand ^[40], and by Goloshchapova, Malamud, and Zastavnyi ^[22,23].

In this work we supplement this picture by demonstrating the absence of positive resonances for $-\Delta_{\alpha, Y}$ with an argument that has the two-fold virtue of being particularly compact as compared to the general settings of ^[40] and ^[22,23], and exploiting the explicit structure of the matrix $\Gamma_{\alpha, Y}(z)$, unlike the abstract reasoning of ^[19] (further comments in this respect are cast at the end of Section 3). As such, the approach that we present here has its own autonomous interest.

Moreover, we have already mentioned that the absence of positive resonances for an ordinary Schrödinger operator $-\Delta+V$ is typically proved with Carleman's estimate, whereas for the singular version $-\Delta_{\alpha, Y}$ it appears to be very hard to use those classical techniques – and indeed our proof relies on a direct analysis based on the explicit formula for the resolvent: this makes any proof of absence of resonances surely valuable.

In Section 2 we present the rigorous context within which our main result is formulated. In particular, we survey the definition and the basic properties of the singular point-perturbed Schrödinger operator $-\Delta_{\alpha, Y}$ and we formulate the precise definition of resonance.

The proof of our main theorem is then discussed in Section 3, together with a few additional comments for comparison with the previous literature.

We conclude our presentation in Section 4 with some final remarks that connect our main theorem with recent dispersive and scattering results for the Schrödinger evolution of the singular point-perturbed Laplacian, and highlight interesting open questions.

Notation. For vectors in $x, y\in\mathbb{R}^d$ the Euclidean norm and scalar products shall be denoted, respectively, by $|x|$ and $x\cdot y$, whereas for the action of an operator (or a matrix, in particular) $A$ on the vector $v$ we shall simply write $Av$. The expression $\delta_{j, k}$ denotes the Kronecker delta. By $\mathbb{1}$ and $\mathbb{O}$ we shall denote, respectively, the identity and the zero operator, irrespectively of which vector space they act on, which will be clear from the context. By $\overline{z}$ and $Z^*$ we shall denote, respectively, the complex conjugate of a scalar $z\in\mathbb{Z}$ and the transpose conjugate of a square matrix $Z$ with complex entries. We shall use the shorthand $\langle x\rangle: = \sqrt{1+x^2}$ for $x\in\mathbb{R}$. By $\mathcal{B}(X, Y)$ we shall denote the space of bounded linear operators from the Banach space $X$ to the Banach space $Y$. For a vector $\psi$ in a Hilbert space $\mathcal{H}$ the $\mathcal{H}\to \mathcal{H}$ rank-one orthogonal projection onto the span of $\psi$ shall be indicated with $|\psi\rangle\langle\psi|$. The rest of the notation is standard or will be declared in due time.

2.   Set-up and main result

Let us start by collecting an amount of well known facts concerning the three-dimensional singular point-perturbed Schrödinger operator that we informally referred to in the course of the previous Section.

Let us fix $N\in\mathbb{N}$, a collection $Y = \{y_1, \dots, y_N\}$ of distinct points in $\mathbb{R}^3$, and a multi-index parameter $\alpha\equiv(\alpha_1, \dots, \alpha_N)\in\mathbb{R}^N$.

For $z\in \mathbb{C}$ and $x, y, y'\in \mathbb{R}^3$, let us set

and

Clearly, the map $z\mapsto \Gamma_{\alpha, Y}(z)$ has values in the space of $N\times N$ symmetric, complex valued matrices, and is entire. Therefore, $z\mapsto \Gamma_{\alpha, Y}(z)^{-1}$ is meromorphic on $\mathbb{C}$ and hence the subset $\mathcal{E}_{\alpha, Y}\subset\mathbb{C}$ of poles of $\Gamma_{\alpha, Y}(z)^{-1}$ is discrete. Let us further define

where $\mathbb{C}^+$ (resp., $\mathbb{C}^-$) denotes as usual the open complex upper (resp., lower) half-plane.

Definition 2.1. Let $z\in\mathbb{C}^+\setminus\mathcal{E}_{\alpha, Y}^+$. The operator $-\Delta_{\alpha, Y}$ is defined on the domain

by the action

It is straightforward to check that at fixed $z$ the decomposition (2.4) of a generic element in $\mathcal{D}(-\Delta_{\alpha, Y})$ is unique, and that the space $\mathcal{D}(-\Delta_{\alpha, Y})$, as well as the action of $-\Delta_{\alpha, Y}$ on a generic function of its domain, are actually independent of the choice of $z$. Moreover, on $H^2$-functions $F$ vanishing at all points of $Y$ one has $-\Delta_{\alpha, Y}F = -\Delta F$.

Theorem 2.2 (Basic properties of the point-perturbed Schrödinger operator).

(i) The operator $-\Delta_{\alpha, Y}$ is self-adjoint on $L^2(\mathbb{R}^3)$ and extends the operator $-\Delta\upharpoonright C_0^\infty(\mathbb{R}^3 \!\setminus\! Y)$. The Friedrichs extension of the latter, namely $-\Delta$ with domain $H^2(\mathbb{R}^3)$, corresponds to the formal choice $\alpha = \infty$ in Definition 2.1.

(ii) If $u\in\mathcal{D}(-\Delta_{\alpha, Y})$ and $u|_{\mathcal{U}} = 0$ for some open subset $\mathcal{U}\subset\mathbb{R}^3$, then $(-\Delta_{\alpha, Y}u)|_{\mathcal{U}} = 0$.

(iii) The set $\mathcal{E}_{\alpha, Y}^+$ of poles of $\Gamma_{\alpha, Y}(z)^{-1}$ in the open complex half-plane consists of at most $N$ points that are all located along the positive imaginary semi-axis, and for $z\in\mathbb{C}^+\setminus\mathcal{E}^+$ one has the resolvent identity

(iv) The essential spectrum of $-\Delta_{\alpha, Y}$ is purely absolutely continuous and coincides with the non-negative half-line, the singular continuous spectrum is absent, and there are no positive eigenvalues:

(v) There is a one-to-one correspondence between the poles $z = \mathrm{i}\lambda\in\mathcal{E}_{\alpha, Y}^+$ of $\Gamma_{\alpha, Y}(z)^{-1}$ and the negative eigenvalues $-\lambda^2$ of $-\Delta_{\alpha, Y}$, counting the multiplicity. The eigenfunctions associated with the eigenvalue $-\lambda^2 < 0$ have the form

where $(c_1, \dots, c_N)\in\ker\Gamma_{\alpha, Y}(\mathrm{i}\lambda)$. In the special case $N = 1$ $(Y = \{y\})$,

and the unique negative eigenvalue, when it exists, is non-degenerate and with eigenfunction $\mathcal{G}_{-4\pi \mathrm{i}\alpha}^{y}$.

Theorem 2.2 is a collection of classical results from ^[24,25,45], which are discussed in detail, e.g., in [4,Sect. Ⅱ.1.1].

In addition to Theorem 2.2, the spectral behavior of $-\Delta_{\alpha, Y}$ on the real line, and in particular the nature of the spectral point $z^2 = 0$, was discussed by one of us in ^[38], and we shall now review those results.

Let us first remark, as emerges from Theorem 2.2, that the eigenvalue zero is absent when $N = 1$, but may occur when $N\geqslant 2$: examples of configurations of the $y_j$'s that produce a null eigenvalue are shown in [38,Sect. 3].

In ^[38] a limiting absorption principle for $-\Delta_{\alpha}$ was established, in the spirit of the classical Agmon-Kuroda theory for the free Laplacian ^[1,33], and a low-energy resolvent expansion was produced, analogously to the case of regular Schrödinger operators with scalar potential ^[1,28].

Theorem 2.3 (^[38]). Let $\sigma > 0$ and let ${\bf{B}}_{\sigma}$ be the Banach space

(i) For every $z\in \mathbb C^+\setminus\mathcal{E}_{\alpha, Y}^+$ one has $(-\Delta_{\alpha, Y}-z^2\mathbb{1})^{-1}\in{\bf{B}}_{\sigma}$, and the map $\mathbb C^+\setminus\mathcal{E}_{\alpha, Y}^+\ni z\mapsto (-\Delta_{\alpha, Y} -z^2\mathbb{1})^{-1}\in {\bf{B}}_{\sigma}$ can be continuously extended to $\mathbb R\setminus\mathcal{E}_{\alpha, Y}^0$.

(ii) In a real neighborhood of $z = 0$, one has the expansion

for some $R_{-2}, R_{-1}\in {\bf{B}}_{\sigma}$ and some continuous ${\bf{B}}_{\sigma}$-valued map $z\mapsto R_0(z)$. Moreover, $R_{-2}\neq \mathbb{O}$ if and only if zero is an eigenvalue for $-\Delta_{\alpha, Y}$.

In view of Theorem 2.3(ⅱ) and of the heuristic idea of a resonance as the existence of a non-$L^2$ solution $u$ to $-\Delta_{\alpha, Y}u = z^2u$, it is natural to say that $-\Delta_{\alpha, Y}$ has a zero resonance when $(-\Delta_{\alpha, Y}-z^2\mathbb{1})^{-1} = O(z^{-1})$ as $z\to 0$, that is, with respect to the low-energy asymptotics (2.7), when $R_{-2} = \mathbb{O}$ and $R_{-1}\neq \mathbb{O}$.

There is an equivalent way to formulate such an occurrence (the proof of which is also deferred to Section 3).

Lemma 2.4. The following facts are equivalent:

(i) in the asymptotics (2.7), $R_{-2} = \mathbb{O}$ and $R_{-1}\neq \mathbb{O}$;

(ii) the matrix $\Gamma_{\alpha, Y}(0)$ is singular, but zero is not an eigenvalue of $-\Delta_{\alpha, Y}$.

The resolvent identity (2.6) and Lemma 2.4 above finally motivate the following precise notion of resonance.

Definition 2.5. Let $z\in\mathbb{C}$. The operator $-\Delta_{\alpha, Y}$ has a resonance at $z^2$ if the matrix $\Gamma_{\alpha, Y}(z)$ is singular, but $z^2$ is not an eigenvalue of $-\Delta_{\alpha, Y}$.

As announced in previous Section's introduction, in this work we focus on the exclusion of positive resonances – informally speaking, to draw a parallel to (1.2), we shall conclude that the equation

admits no non-trivial solutions, be they in $L^2(\mathbb{R}^3)$ (non-existence of embedded eigenvalues, as seen already in Theorem 2.2(iv)) or outside of $L^2(\mathbb{R}^3)$ (non-existence of embedded resonances). More precisely, in this work we prove that for any $z\in\mathbb{R}\setminus\{ 0\}$, the spectral point $\mu = z^2$ is not a resonance.

In view of Definition 2.5 and Theorem 2.2(iv), the absence of positive resonances is tantamount as the non-singularity of $\Gamma_{\alpha, Y}(z)$ for any $z\in\mathbb{R}\setminus\{ 0\}$. This is precisely the form of our main result.

Theorem 2.6. For every $\alpha\in \mathbb R^N$, every collection $Y = \{y_1, \ldots y_N\}$ of $N$ distinct points in $\mathbb{R}^3$, and every $z\in\mathbb{R}\setminus\{ 0\}$, the matrix $\Gamma_{\alpha, Y}(z)$ is non-singular. Equivalently, the self-adjoint operator $-\Delta_{\alpha, Y}$ has no real positive resonances.

Theorem 2.6, combined with Theorem 2.3, essentially completes the picture of the spectral theory for three-dimensional Schrödinger operators with finitely many point interactions.

For completeness of presentation, at the end of the proof in the following Section, we shall comment on the comparison with the previous literature on the positive resonances of $-\Delta_{\alpha, Y}$, and in the subsequent Section 4 we shall connect our main theorem with recent dispersive and scattering results for the propagator $e^{ \mathrm{i} t \Delta_{\alpha, Y}}$, together with some relevant open problems.

3.   Proof of the main Theorem and additional remarks

This Section is mainly devoted to the proof of Theorem 2.6. In terms of the notation (2.3), one has to prove that the set $\mathcal{E}_{\alpha, Y}^0\setminus\{0\}$ of non-zero poles of $\Gamma_{\alpha, Y}(z)^{-1}$ on the real line is empty.

In practice it suffices to only consider $z > 0$, for $\Gamma_{\alpha, Y}(-z) = \Gamma_{\alpha, Y}(z)^*$ for any $z\in\mathbb{R}$, and hence $\mathcal{E}_{\alpha, Y}^0$ is symmetric with respect to $z = 0$.

In fact, $\mathcal{E}_{\alpha, Y}^0$ is also finite. Indeed, $\mathcal{E}_{\alpha, Y}^0\subset\mathcal{E}_{\alpha, Y}$ is a discrete set, and for $z\in\mathbb{R}$ one has $\Gamma_{\alpha, Y}(z) = -\frac{ \mathrm{i} z}{4\pi}\mathbb{1}+\Lambda_{\alpha, Y}(z)$ where the matrix norm of $\Lambda_{\alpha, Y}(z)$ is uniformly bounded in $z$, therefore $\Gamma_{\alpha, Y}(z)$ is invertible for large enough $z$.

Let us first present the proof of Lemma 2.4, which was at the basis of the definition of resonance for $-\Delta_{\alpha, Y}$.

Proof of Lemma 2.4. We recall [28,Section 2] that $(-\Delta-z^2\mathbb{1})^{-1}\in{\bf{B}}_{\sigma}$ for every $z\in \mathbb C^+$, and the map $\mathbb C^+\ni z\mapsto (-\Delta -z^2\mathbb{1})^{-1}\in {\bf{B}}_{\sigma}$ can be continuously extended to the real line. Moreover, we observe that the map $\mathbb R\ni z\to |\mathcal{G}^{y_1}_z\rangle\langle\overline{\mathcal{G}^{y_2}_z}|\in{\bf{B}}_{\sigma}$ is continuous for any $y_1, y_2\in \mathbb R^3$. Owing to these facts, one compares the limits $z\to 0$ in the resolvent identity (2.6), in the resolvent expansion (2.7), and in the low-energy expansion

established in [38,Proposition 5], and concludes

We can now prove the desired equivalence.

(ⅰ)$\, \Rightarrow\,$(ⅱ). Since $R_{-2} = \mathbb{O}$, Theorem 2.3(ⅱ) guarantees that $z = 0$ is not an eigenvalue for $-\Delta_{\alpha, Y}$. Moreover, since $R_{-1}\neq \mathbb{O}$, by (3.1) also $A_{-1}\neq \mathbb{O}$, which implies in particular that $\Gamma_{\alpha, Y}(0)$ is singular.

(ⅱ)$\, \Rightarrow\,$(ⅰ). Since $z = 0$ is not an eigenvalue for $-\Delta_{\alpha, Y}$, Theorem 2.3(ⅱ) guarantees that $R_{-2} = \mathbb{O}$, whence also $A_{-2} = \mathbb{O}$ owing to (3.1). Since $\Gamma_{\alpha, Y}(0)$ is singular, necessarily $A_{-1}\neq \mathbb{O}$, whence also $R_{-1}\neq \mathbb{O}$ again owing to (3.1).

Our argument for Theorem 2.6 is based upon the following useful result in linear algebra. We denote by $Sym_{N}(\mathbb R)$ the space of $N\times N$ symmetric real matrices.

Lemma 3.1. Let $A, B \in Sym_{N}(\mathbb R)$, and assume furthermore that $B$ is positive definite. Then $A- \mathrm{i} B$ is non-singular.

Proof. Suppose for contradiction that $A- \mathrm{i} B$ is singular. Then there exist $v, w\in \mathbb R^n$, at least one of which is non-zero, such that

We can exclude for sure that $v = 0$, for in this case $Bw+ \mathrm{i} Aw = 0$, whence in particular $Bw = 0$ and therefore $w = 0$, against the assumption that $v+ \mathrm{i} w\neq 0$. Applying $B^{-1}$ to the identity above, separating real and imaginary parts, and setting $C: = B^{-1}A$, one gets

which implies $C^2v = -v$. As $v\neq 0$, the conclusion is that $-1$ is an eigenvalue for $C^2$. However, the matrix $C = B^{-1}A$ can be diagonalised over $\mathbb R$, since it is similar to the symmetric matrix $B^{-1/2}AB^{-1/2}$, and therefore $C^2$ is similar to a positive semi-definite matrix. Hence $-1$ cannot be in the spectrum of $C^2$.

We shall also make use of the following property.

Lemma 3.2. Let $N, d\in\mathbb{N}$ and let $y_1, \dots, y_N$ be distinct vectors in $\mathbb{R}^d$. Then there exists a unit vector $a\in\mathbb{R}^d$ such that the numbers $a\cdot v_j$ are all distinct.

Proof. Let $\sigma_d$ be the area measure on the unit sphere $\mathbb{S}^{d-1}: = \{x\in \mathbb R^d\; |\; |x| = 1\}$. For every pair $(j, k)$, with $j, k\in\{1, \ldots, N\}$, $j\neq k$, let us consider the set

Since $y_j\neq y_k$, $P_{jk}$ is an hyperplane in $\mathbb R^d$, whence $\sigma_d(\mathbb{S}^{d-1}\cap P_{jk}) = 0$. It follows that the set

satisfies $\sigma_{d}(Q) = \sigma_{d}(\mathbb{S}^{d-1})$, and in particular $Q$ is non-empty. If we choose $a\in Q$, then by construction the numbers $a\cdot y_j$ are all distinct.

We are ready to prove the main Theorem.

Proof of Theorem 2.6. Let $z > 0$: owing to the scaling property

valid for every $\lambda > 0$, it is enough to prove that $\Gamma_{\alpha, Y}: = \Gamma_{\alpha, Y}(1)$ is non-singular for any choice of $\alpha$ and $Y$.

Using (2.1) and (2.2), we re-write

with $A_{\alpha, Y}, B_Y\in Sym_N(\mathbb R)$ given explicitly by

where the real function ${\rm{sinc}}(x)$ is defined by

Owing to Lemma 3.1, the thesis follows when one proves that $B_Y$ is positive definite.

Based on an immediate integration in polar coordinates, it is convenient to express

where $\sigma_2(p)$ is the area measure on the unit sphere $\mathbb{S}^2\subset \mathbb R^3$.

For generic $v\equiv(v_1, \dots, v_N)$ let $B_Y[v]: = v\cdot B_Y v$ be the quadratic form associated to $B_Y$. Since

then the matrix $B_Y$ is positive semi-definite.

To demonstrate that $B_Y$ is actually positive definite, we specialise to the present context the clever argument by Castel, Filbir, and Szwarc ^[46] (in ^[46] the general question of linear independence of exponential maps over subsets of $\mathbb R^d$ is addressed).

Assume that for some $v\in\mathbb{R}^3\setminus\{0\}$ one has $B_Y v = 0$. From the above computation of $B_Y[v]$, one deduces

We show, by induction on $N$, that the latter identity implies $v = 0$. The case $N = 1$ is obvious. Let $N\geqslant 2$, and in (*) let us consider all possible $p\in\mathbb{S}^2$ of the form

for two fixed vectors $a, b\in\mathbb{R}^3$ such that $|a| = |b| = 1$, $a\perp b$, and the scalars $\alpha_j: = y_j\cdot a$ with $j\in\{1, \dots, N\}$ are all distinct. Lemma 3.2 ensures that this choice of $a$ and $b$ is possible. It is non-restrictive to assume $\alpha_1 < \cdots < \alpha_N$, and let us also set $\beta_j: = y_j\cdot b$, $j\in\{1, \dots, N\}$. Then (*) reads

In fact, since l.h.s. above depends analytically on $t$, such an identity holds true for every $t\in \mathbb C$. Specialising it for $t = - \mathrm{i}\tau$, $\tau\in\mathbb{R}$, it takes the form

and also, upon dividing by $\exp(\alpha_N\sinh\tau+ \mathrm{i}\beta_N\cosh\tau)\neq 0$,

Since $\alpha_j < \alpha_N$ for $j < N$, taking in the latter expression $\tau$ arbitrarily large and positive implies necessarily $v_N = 0$. By the inductive assumption that (*) implies $v = 0$ when it is considered with $N-1$ instead of $N$, one concludes that also $v_1 = \cdots = v_{N-1} = 0$.

In the remaining part of this Section we comment on how the absence of positive resonances for $-\Delta_{\alpha, Y}$ could be also read out from the already mentioned works ^{[19,22,23,40]}.

In ^[40] a wide class of compactly supported perturbations of the Laplacian is considered, originally introduced by Sjöstrand and Zworski in ^[41], which include the case of finitely many point interactions. For these operators, it is proved a version of Rellich's uniqueness theorem [40,Theorem 2.4], which combined with the classical unique continuation principle for the Laplacian allows to deduce the absence of real resonances.

In ^[19] the reasoning is based upon the computation of the residue of $\Gamma_{\alpha, Y}(z)^{-1}$ at a generic pole $z\in\mathcal{E}_{\alpha, Y}^0$. From the resolvent identity (2.6) and the information that $(-\Delta_{\alpha, Y}-z^2\mathbb{1})^{-1}$ is bounded for $z\in\mathbb{R}\setminus\{0\}$ it is shown that one obtains instead a zero value, thus contradicting the fact that $z$ is a pole. This approach requires the additional knowledge (already available, as seen in Theorem 2.2(iv)) that $-\Delta_{\alpha, Y}$ has no positive eigenvalues, and by-passes the explicit structure of the matrix $\Gamma_{\alpha, Y}(z)$.

In ^[22,23] strictly speaking no reference to (positive) resonances of $-\Delta_{\alpha, Y}$ is made. The operator $-\Delta_{\alpha, Y}$ and its main properties are recovered by means of the alternative framework of boundary triplets and Weyl function. Then the positive definiteness of what we denoted here by $B_Y$ is indirectly alluded to by considering another matrix, with similar structure, and proving for the latter the positive definiteness by means of general properties of positive definite functions like our $\mathrm{sinc}(x)$.

4.   Connection with the dispersive properties of $e^{ \mathrm{i} t\Delta_{\alpha, Y}}$ and open problems

In this short, concluding section we return to the general subject of the dispersive properties of the singular point-perturbed Schrödinger equation (1.4), in order to emphasize the connection of our Theorem 2.6 with recent dispersive and scattering results for the propagator $e^{ \mathrm{i} t\Delta_{\alpha, Y}}$.

We already mentioned in Section 1 that $L^1\to L^\infty$ (and hence by interpolation general $L^p\to L^q$) dispersive estimates for $e^{ \mathrm{i} t\Delta_{\alpha, Y}}$ were first proved in ^[15], with a suitable weight that accounts for the singularity $|x-y_j|^{-1}$ of $e^{ \mathrm{i} t\Delta_{\alpha, Y}}f$, whereas the subsequent work ^[26] reproduced such estimates without weight in the regime $q\in[2, 3)$. In ^[15] the authors assumed that $\Gamma_{\alpha, Y}(z)$ be non-singular for every $z\in\mathbb{R}$. In view of Theorem 2.6, it is sufficient to impose that zero is regular (it is neither an eigenvalue, nor a resonance).

On a related note, we mentioned that in ^[16] the $L^p$-boundedness for the wave operator for the pair $(-\Delta_{\alpha, Y}, -\Delta)$ was proved for all possible $p$'s, namely $p\in(1, 3)$, under the implicit assumption of the absence of a zero eigenvalue as well as of positive resonances. The latter condition is now rigorously confirmed by Theorem 2.6.

This rises up two interesting open questions. First, one would like to investigate whether, in the spirit of ^[15,26], an obstruction at zero in the form of a zero-energy eigenvalue or resonance would still allow one to derive certain (possibly slower) dispersive estimates. Let us recall that the counterpart problem for ordinary Schrödinger operators with spectral obstruction at zero has been intensively studied, significantly by Jensen and Kato ^[28], Rauch ^[36], Rodnianski and Schlag ^[37], Erdoğan and Schlag ^[18] and Yajima ^[42].

Analogously, it would be of interest to understand whether some $L^p$-boundedness of the wave operator for $(-\Delta_{\alpha, Y}, -\Delta)$ is still valid in the presence of a zero-energy eigenvalue, thus extending the recent results by Goldberg and Green ^[21] and by Yajima ^[43,44] on the $L^p$-boundedness of wave operators for the ordinary Schrödinger operators with threshold singularities.

These are questions that surely deserve to be investigated.

Acknowledgments

Partially supported by the Alexander von Humboldt foundation.

Conflict of interest

All authors declare no conflicts of interest in this paper.