Implementation of Yang residual power series method to solve fractional non-linear systems

1.   Introduction

This paper is devoted to providing bounds on the location of discrete eigenvalues of operators of the form

acting on $[L^2(\mathbb{R}^d)]^d,$ the Hilbert space of vector fields with components in $L^2(\mathbb{R}^d)$. The operator (1.1) is introduced in a standard way as an $m$-accreative operator obtained as a form sum of the Lamé operator $-\Delta^\ast$ with domain $[H^1(\mathbb{R}^d)]^d$ and a relatively form-bounded potential $V.$ We refer to Appendix A for details. The free Lamé operator $-\Delta^\ast$ is defined as the (self-adjoint) Friedrichs extension of the minimal operator

with $\Delta: = \Delta I_{ \mathbb{C}^{d\times d}}$ the diagonal matrix with the Laplace operator on the diagonal and where $\text{ div }$ is the standard divergence operator (refer to [24,Sec. 2.1] for details). Here we slightly abused the notation using the same symbol $-\Delta^\ast$ both for the symmetric operator (1.2) and its self-adjoint realisation. The two material-dependent constants $\lambda, \mu\in \mathbb{R}$ (tipically called Lamé's coefficients) are assumed to satisfy the standard conditions

which guarantee the strong ellipticity of $- \Delta^\ast$ (refer to [24,Sec. 2.2] for the a brief summary in matter of ellipticity for systems of second order differential operators). $V$ represents the operator of multiplication by $VI_{ \mathbb{C}^{d\times d}},$ with $V\colon \mathbb{R}^d \to \mathbb{C}.$ Being $V$ complex-valued leads to a (possibly) non self-adjoint setting.

The distribution of eigenvalues of self-adjoint operators has been intensively studied for several decades and nowadays the usage of powerful techniques such as spectral theorem and variational principles have become the standard approach for addressing this issue. However, as these tools are no longer available in a non self-adjoint setting, the generalization of spectral bounds to the non self-adjoint framework is not straightforward and requires a diverse strategy.

Nevertheless, a systematic, albeit recent, approach to the study of eigenvalue estimates for perturbed operators with complex-valued potentials, in particular for Schrödinger operators $-\Delta + V$, has been successfully developed and has already a bibliography. It is very well known consequence of Sobolev inequalities that if $V$ is real-valued then the distance from the origin of every eigenvalue $z$ lying in the negative semi-axis (discrete eigenvalue) can be bounded in terms of $L^p$-norm of the potential (see ^[49,57] and ^[13] for a more recent improvement), more precisely the following bound

holds true for every $\gamma \geq 1/2$ if $d = 1,$ $\gamma > 0$ if $d = 2$ and $\gamma\geq 0$ if $d\geq 3,$ with $C_{\gamma, d}$ independent of $V.$

In 2001, Abramov, Aslanyan and Davies ^[1] generalized the previous result to complex-valued potentials showing that every discrete eigenvalue $z$ of the one-dimensional Schrödinger operator $-d^2/dx^2 + V,$ i.e., $z\in \mathbb{C} \setminus [0, \infty),$ lies in the complex plane within a $1/4 ||V||_{L^1(\mathbb{R})}^2$ radius of the origin. In order to overcome the lack of the aforementioned self-adjoint-based tools, the authors introduced a strategy built on the Birman-Schwinger principle (see, the pioneering Birman's paper ^[7] or the recent work ^[46] for a systematic exposition of abstract Birman-Schwinger principles and some interesting explicit applications), which has permeated all the subsequent works on this topic. Ten years later Frank ^[40] went beyond the one dimensional restriction proving the validity of (1.4) for any $d\geq 2$ and $0 < \gamma\leq 1/2,$ covering also the endpoint case $\gamma = 0$ if $d\geq 3.$ His notable accomplishment derives from replacing a pointwise bound for the Green function of $-d^2/dx^2 -z,$ used in ^[1] and highly sensitive of the one dimensional framework, with the much deeper uniform resolvent estimates due to Kenig, Ruiz and Sogge ^[50] (see below for further details). The result in ^[40] partially proved the conjecture in ^[52] according to which (1.4) holds for any $d\geq 2$ and $0 < \gamma\leq d/2$, leaving opened only the case $1/2 < \gamma < d/2.$ This range was covered later by Frank and Simon in ^[45] for radial potential. In the same work the general case was investigated too. In these respects the authors provided the construction of a sequence of real-valued potentials $V_n$ with $||V_n||_p \to 0$ for $d\geq 2$ and for any $p > (d+1)/2$ such that $-\Delta + V_n$ has eigenvalue 1. Even if this does not disprove the conjecture, as it is stated for discrete eigenvalues only, the result is still relevant in view of the recent generalization in ^[45] of Frank's work ^[40] to positive (embedded) eigenvalues. Recent related results can be also found in ^[41].

In passing, let us stress that the extensive bibliography devoted to spectral bounds in a non self-adjoint context that the appearance of the cited work ^[1] stimulated, did not remain long confined to the class of Schrödinger operators and more generally to second order differential operators. Indeed subsequent investigations have shown the robustness of the approach introduced in ^[1] by fruitfully testing it to several other models. Without attempting to be exhaustive we mention ^[48] for an adaptation to the discrete Schrödinger setting, see also ^[51] where matrix-valued damped wave operators are concerned (see also ^[15] for the elasticity counterpart). Lower order operators, such as Dirac or fractional Schrödinger models, are investigated in ^{[18,26,28,30,31,36,58]} (see also ^[27,39]) and in ^[17] respectively in the continuous and discrete scenario; as for higher order operators refer to ^[34,47]. Associated spectral stability results obtained with different techniques and related tools can be found in ^{[3,10,11,16,25,35,44]}.

As a final remark, we mention that since the turn of the millennium, when the interest on non self-adjoint operators started to take hold, this subject has exhibited a rapid development and nowadays the analysis of such Hamiltonians spreads out to cover several diverse branches of spectral theory. Just with the aim of giving a (non comprehensive) overview of the several questions posed in this context, focusing on Schrödinger operators, we cite ^[8,19,43,62] in matter of existence and estimate on the number of eigenvalues, ^{[9,32,33,41,42,52]} as regards with the non self-adjoint counterpart of Lieb-Thirring inequalities and we refer to the aformentioned quoted works for the topic related to eigenvalues bounds for non self-adjoint Schrödinger operators.

The main aim of our paper is to investigate on spectral properties in the elasticity setting, specifically providing bounds on the distribution of discrete (possibly complex) eigenvalues of Lamé operators (1.2) in terms of $L^p$-norm of the potential. In other words we want to show the validity of a suitable analogous of (1.4) in this diverse framework.

A first motivation to our purpose is the trivial observation that in $d = 1$ our operator $-\Delta^\ast$ turns into a scalar differential operator and, even more relevant, it is simply a multiple of the Laplacian, more precisely

Therefore, by virtue of the aforementioned analogous theorem in ^[1] (cfr Theorem 4) for non self-adjoint Schrödinger operators $-d^2/dx^2 + V,$ the following result comes as no surprise and it is indeed its straightforward consequence.

Theorem 1.1. Let $d = 1$ and assume that $V\in L^1(\mathbb{R}).$ Then any eigenvalue $z\in \mathbb{C} \setminus [0, \infty)$ of the perturbed Lamé operator $-\Delta^\ast + V$ satisfies

The validity of Theorem 1.1, compared with Theorem 4 in ^[1], together with the aforementioned extensions to higher dimensional Schrödinger operators, motivate the question of whether an estimate similar to (1.6) holds true in $d\geq 2.$

As already said, starting from the celebrated paper of Abramov, Aslanyan and Davies ^[1], the usage of the Birman-Schwinger principle has been recognized as a crucial tool to get this type of bounds in a non self-adjoint setting of problem. Roughly speaking this permits to re-phrase conveniently the eigenvalues problem associated to a perturbed operator, say $H_0 + V,$ in terms of the eigenvalues problem of an integral operator suitably related to the latter. Precisely, the following proposition holds true.

Proposition 1.1 (Birman-Schwinger principle). Let $z\notin \sigma(H_0).$ Then

with $V_{1/2}: = |V|^{1/2} \text{ sgn }(V).$ (If $z\in \mathbb{C},$ then $\text{ sgn }(z)$ represents the complex signum function defined by $\text{ sgn }(z): = z/|z|$ if $z\neq 0$ and $\text{ sgn }(0): = 0$).

The formal statement of Proposition 1.1 can be made rigorous under suitable hypotheses on $H_0$ and $V.$ We refer again to the recent work [46,Thm. 6,Thm.7] by Hansmann and Krejčiřík for a rigorous exposition of abstract Birman-Schwinger principles.

In the context of elasticity, $H_0$ is replaced by the Lamé operator $-\Delta^\ast,$ therefore the corresponding Birman-Schwinger operator has the form $V_{1/2} (-\Delta^\ast-z)^{-1} |V|^{1/2}$ and makes sense if $z \in \mathbb{C} \setminus [0, \infty)$ as $\sigma(-\Delta^\ast) = [0, \infty).$ This leads to the need for an explicit expression of the resolvent operator associated with $-\Delta^\ast.$ In this regards it turns out that $(-\Delta^\ast -z)^{-1}$ has a favorable form, i.e.,

where $\mathcal{P}$ is the so-called Leray projection operator, customarily used in elasticity to decompose $[L^2(\mathbb{R}^d)]^d$ vector fields into a divergence-free component plus a gradient (refer to the preliminary section for further details).

This expression shows that, as soon as the decomposition, also known as Helmholtz decomposition, $g = \mathcal{P}g + (I-\mathcal{P})g$ is operated, the resolvent operator $(-\Delta^\ast - z)^{-1}$ splits into a sum of two vector-valued resolvent operators associated with the Laplacian acting, respectively, on the components $\mathcal{P}g$ and $(I-\mathcal{P})g$ of $g.$ This fact, together with the validity of the corresponding results for Schrödinger, strongly suggests a positive answer to our question of whether Theorem 1.1 extends to higher dimensions and indeed it is confirmed by the following result proved in this paper.

Theorem 1.2. Let $d\geq 2$ and assume that $V\in L^{\gamma + \frac{d}{2}}(\mathbb{R}^d),$ with $0 < \gamma\leq 1/2$ if $d = 2$ and $0\leq \gamma\leq 1/2$ if $d\geq 3.$ Then any eigenvalue $z \in \mathbb{C} \setminus [0, \infty)$ of the perturbed Lamé operator $-\Delta^\ast + V$ satisfies

with a constant $C_{\gamma, d, \lambda, \mu}$ independent of $V.$

Incidentally, observe that if $d\geq 3$ and $\gamma = 0,$ the previous theorem provides a sufficient condition which guarantees absence of discrete eigenvalues of $-\Delta^\ast + V$ with $V\in L^\frac{d}{2}(\mathbb{R}^d).$ More specifically, if

then the discrete spectrum $\sigma_d(-\Delta^\ast + V)$ is empty. In comparison with this result, we should mention a prior work of the present author ^[24] also related to the problem of establishing sufficient conditions which disprove presence of eigenvalues. In ^[24], with a completely different approach based on the multipliers method as previously applied to Schrödinger operators by Fanelli, Krejčiřík and Vega in ^[37,38], total absence of eigenvalues, i.e., absence of both discrete and embedded eigenvalues, of $-\Delta^\ast + V$ was proved under the following Hardy-type subordination

where $\Lambda$ is a suitable small constant (see condition (4) in ^[24]).

As (1.9), condition (1.10) is intrinsically a smallness condition, it is true, on the other hand it is satisfied by potentials with quite rough local singularities, e.g., $|x|^{-2},$ which, instead, are ruled out by the $L^p$-type condition in Theorem 1.2.

In attempt of including potentials with local stronger singularities, such as inverse-square type, we generalize Theorem 1.2 by measuring the size of the potential in (1.8) in terms of less restrictive norms. Specifically, as first generalization, we consider potentials in the Morrey-Campanato class $\mathcal{L}^{\alpha, p}(\mathbb{R}^d)$ which is defined for $\alpha > 0$ and $1\leq p\leq d/\alpha$ by

In passing, notice that $1/|x|^\alpha\in L^{d/\alpha, \infty}(\mathbb{R}^d)\subset \mathcal{L}^{\alpha, p},$ for $\alpha > 0$ (we emphasize particularly the case $\alpha = 2$) and $1\leq p < d/\alpha,$ however $1/|x|^\alpha\notin L^{d/\alpha} = \mathcal{L}^{\alpha, d/\alpha}.$

More precisely we shall prove the following theorem.

Theorem 1.3. Let $d\geq 2$ and assume that $V\in \mathcal{L}^{\alpha, p}(\mathbb{R}^d)$ with $(d-1)(2\gamma + d)/ 2(d-2\gamma) < p\leq \gamma + {d/2}$ and let $0 < \gamma\leq 1/2$ if $d = 2$ and $0\leq \gamma\leq 1/2$ if $d\geq 3.$ Then any eigenvalue $z \in \mathbb{C} \setminus [0, \infty)$ of the perturbed Lamé operator $-\Delta^\ast + V$ satisfies

with $\alpha = 2d/(2\gamma+d)$ and a constant $C_{\gamma, p, d, \lambda, \mu}$ independent of $V.$

As a byproduct, in higher dimensions, the previous theorem provides a sufficient condition to guarantee absence of discrete eigenvalues. More precisely, the following corollary is immediate consequence of Theorem 1.3.

Corollary 1.1. Let $d\geq 3,$ $(d-1)/2 < p\leq d/2$ and assume $V\in \mathcal{L}^{2, p}(\mathbb{R}^d)$ and

with $C_{0, p, d, \lambda, \mu}$ as in Theorem 1.3 when $\gamma = 0.$ Then the perturbed Lamé operator $-\Delta^\ast + V$ has no eigenvalue in $\mathbb{C}\setminus [0, \infty).$

Observe that Theorem 1.3 does extend Theorem 1.2, indeed from

which holds true for $1\leq p\leq \gamma + d/2,$ in particular it follows that

which immediately gives (1.8) as a consequence of (1.11).

Actually our investigation goes even further providing eigenvalues bounds of the type (1.8) in terms of potentials belonging to the Kerman-Saywer space $\mathcal{KS}_\alpha(\mathbb{R}^d),$ which is defined for $0 < \alpha < d$ by

where the supremum is taken over all dyadic cubes $Q$ in $\mathbb{R}^d.$

As $\mathcal{L}^{\alpha, p}(\mathbb{R}^d)\subset \mathcal{KS}_{\alpha}(\mathbb{R}^d)$ if $p\neq 1$ (see Section 2 in ^[5]) it is true that the Kerman-Sayer class is wider than the Morrey-Campanato class. On the other hand, it turns out that assuming solely $V\in \mathcal{KS}_\alpha$ is not enough to get bound (1.11) with $||V||_{\mathcal{L}^{\alpha, p}(\mathbb{R}^d)}$ replaced by $||V||_{\mathcal{KS}_\alpha(\mathbb{R}^d)}.$ Indeed, additionally, we will ask the potential $V$ to belong to the Muckenhoupt $A_2(\mathbb{R}^d)$ class of weights which is defined, in general, for $1 < p < \infty$ as the set of measurable non-negative function $w$ such that

where $Q$ is any cube in $\mathbb{R}^d$ and $C$ is a constant independent of $Q.$

More precisely, we shall prove the following result.

Theorem 1.4. Let $d\geq 2$ let $1/3\leq \gamma < 1/2$ if $d = 2$ and $0\leq \gamma < 1/2$ if $d\geq 3.$ Let $\alpha = \frac{2d}{2\gamma +d} \beta$ and $\beta = \frac{(d+2\gamma)(d-1)}{2(d-2\gamma)}$ and assume that $V^\beta\in\mathcal{KS}_\alpha(\mathbb{R}^d).$ If moreover $V\in A_2(\mathbb{R}^d),$ then any eigenvalue $z\in \mathbb{C}\setminus [0, \infty)$ of the perturbed Lamé operator $-\Delta^\ast + V$ satisfies

with a constant $C_{\gamma, d, \lambda, \mu}$ independent of $V.$

As a consequence of the previous result, one gets the following corollary on absence of discrete eigenvalues.

Corollary 1.2. Let $d\geq 3$ and assume $V^{\frac{d-1}{2}}\in \mathcal{KS}_{d-1}(\mathbb{R}^d),$ $V\in A_2(\mathbb{R}^d)$ and

with $C_{0, d, \lambda, \mu}$ as in Theorem 1.4 when $\gamma = 0.$ Then the perturbed Lamé operator $-\Delta^\ast + V$ has no eigenvalue in $\mathbb{C}\setminus [0, \infty).$

It is worth comparing Theorem (1.4) with the analogous result in ^[53] (cfr. Theorem 1.1) for Schrödinger operators. Here, the bound (1.13) was obtained without the additional assumption $V\in A_2(\mathbb{R}^d),$ then showing a peculiar feature of the Lamé operator.

Roughly speaking, the philosophy is that thanks to the Helmholtz decomposition which, as shown in (1.7), makes the resolvent operator $(-\Delta^\ast -z)^{-1}$ "behave" like a sum of two resolvent $(-\Delta-z)^{-1},$ at first we can perform our analysis at the level of the much more investigated Schrödinger operators, estimating the two components in (1.7) separately. Then, in order to get bound (1.11) and (1.13), respectively, these two pieces have to be recombined together on weighted $L^2$-spaces and it is at this step that the $A_2$ assumption comes into play. We refer to Section 4 for greater details.

Marginally, observe that, even though the aforementioned strategy underpins the proof of both (1.11) and (1.13), good property of the $\mathcal{L}^{\alpha, p}$ class allowed to drop the $A_2$ assumption in Theorem 1.3 (see Lemma 2.8 below).

It is very well known fact that distinguishing whether $-\Delta + V$ has finite (possibly empty) or infinite discrete spectrum depends on the large $x$ fall-off of the potential. More specifically, it is mainly consequence of the uncertainty principle, quantified by the Hardy inequality

that the borderline is marked by an inverse-square type behavior at infinity.

In scattering theory and in particular in matter of determining existence of wave operators, again the large $x$ behavior of the potential plays a central role, in this context the threshold is given by a Coulomb-type asymptotic decay and positive results require $|x|^{-\alpha}$ with $\alpha > 1.$ Observe that our previous theorems essentially restrict to $|x|^{-\alpha}$ decay with $\alpha > 2d/d+1$ and therefore are not fully satisfactory in the perspective of their possible application to stationary scattering theory. The gap is filled in the following result by using a suitable interpolation argument.

Theorem 1.5. Let $d\geq 2,$ $\gamma > \frac{1}{2}$ and $\alpha > \gamma -\frac{1}{2}.$ Let $q = 2\gamma + (d-1)/2$ and assume $V\in L^q(\langle x \rangle^{2\alpha}dx).$ Then any eigenvalue $z\in \mathbb{C}\setminus [0, \infty)$ of the perturbed Lamé operator $-\Delta^\ast + V$ satisfies

with a constant $C_{\gamma, \alpha, d, \lambda, \mu}$ independent of $V.$

Here we used the notation $\langle x\rangle: = (1+ |x|^2)^{1/2},$ moreover, given a measurable function $w,$ $L^p(w dx)$ stands for the $w$-weighted $L^p$ space on $\mathbb{R}^d$ with measure $w(x) dx.$

The rest of the paper is organized as follows: In Section 2, in attempt of making the paper sufficiently self-contained, we collect some preliminary facts on the Helmholtz decomposition. Here, some properties of Lamé operator are provided too. Among them, particular emphasis will be given to uniform estimates for the resolvent operator $(-\Delta^\ast -z)^{-1}$ which will represent the main ingredient in the proof of our aforementioned results.

Contrarily to the much more investigated Schrödinger operator, up to our knowledge, eigenvalue bounds of the form (1.4) for the perturbed Lamé operator are unknown even in the self-adjoint situation. Although in this case the proof of (1.4) follows almost verbatim the instead-well-known one for Schrödinger, we decided to dedicate Section 3 to prove it anyhow. The advantage of this choice comes out in the possibility of explicitly showing the deep differences and difficulties that arise passing from the self-adjoint to the non self-adjoint framework which, instead, is fully analyzed in Section 4. In particular, Section 4 is devoted to the proof of the main results stated in the introduction. We mention that in the meantime the publication process for this paper was finalized, the results presented here were generalised in ^[14] to cover matrix-valued potential and, more importantly, embedded eigenvalues as well.

Notations.

In this paper both scalar and vector-valued functions are considered. In attempt of produce no confusion, we clarify here that notation like $f, g, u$ are reserved for vector-fields, instead $\phi, \psi$ are set aside for scalar functions.

When the letter adopted for denoting a vector field contains already a subscript, the standard subscript notation $\cdot_j$ for the $j$-th component will be replaced by the superscript $\cdot^{(j)},$ e.g., the $j$-th component of the vector field $u_S$ is indicated by $u_S^{(j)}.$

We use the following definition for the $L^p$- norm of a vector field $u\in [L^p(\mathbb{R}^d)]^d:$

In order to lighten the presentation, in the following we usually abbreviate both $||\cdot||_{L^p(\mathbb{R}^d)}$ and $||\cdot||_{[L^p(\mathbb{R}^d)]^d},$ the $L^p$-norm of scalar and vector-valued functions, respectively, with $||\cdot||_p.$ This, in general, could create some ambiguity, on the other hand, it will be clear from the context and the notation used there if $||\cdot||_p$ stands for one or the other norm.

Again, the notation $\langle \cdot, \cdot \rangle$ will denote both $\langle \cdot, \cdot \rangle_{L^2(\mathbb{R}^d)}$ and $\langle \cdot, \cdot \rangle_{[L^2(\mathbb{R}^d)]^d},$ where the latter extends in an obvious way the usual definition of the former, namely given $f, g \in [L^2(\mathbb{R}^d)]^d,$ one defines

Let $E$ and $F$ be two Banach spaces and let $T\colon E \to F$ be a bounded linear operator from $E$ into $F.$ The notation $||T||_{E\to F}$ will be used to denote the operator norm of $T.$

2.   Preliminaries

This section is concerned with recalling some properties connected with the Helmholtz decomposition together with stating and proving some related consequences on Lamé operators that will strongly enter the proof of our results later. If the first part wants to be just a remainder of very well known results on the Helmholtz decomposition and therefore can be safely skipped by any reader already familiar with this topic, the subsequent subsections, namely Subsection 2.2 and Subsection 2.3, represent an important part of the paper. More specifically, it is there that uniform resolvent estimates for the resolvent operator $(-\Delta^\ast-z)^{-1}$ are provided which are important in their own right.

2.1.   Helmholtz decomposition

Theorem 2.1 (Helmholtz decomposition). Let $d\geq 2$ and let $\Omega\subseteq \mathbb{R}^d$ be either an open, bounded, simply connected, Lipschitz domain or the entire $\mathbb{R}^d.$ Then any square-integrable vector field $f = (f^{(1)}, f^{(2)}, \dots, f^{(d)}) \in [L^2(\Omega)]^d$ can be uniquely decomposed as

where $f_S$ is a divergence-free vector field with null normal derivative and $f_P$ is a gradient. Moreover the two components are orthogonal in a $L^2$- sense. Specifically, the Pythagorean identity

holds true.

Proof. Though the proof simply relies on classical techniques, for sake of completeness we provide a brief sketch of the proof.

Let $f\in [L^2(\Omega)]^d$ and let $\nu$ denote the unit outward normal vector at the boundary $\partial \Omega.$ We consider the trivial decomposition $f = f-\nabla \psi + \nabla \psi.$ In order for $f-\nabla \psi$ to be the divergence-free component $f_S$ of $f$ with null normal derivative, $\psi$ must satisfy the boundary value problem

In passing, observe that (2.2) is a Poisson problem with Neumann boundary conditions, therefore it admits a unique solution $\psi\in H^1(\Omega)$ modulo additive constants.

Let us now prove the uniqueness of the decomposition. Let $f_1, f_2 \in [L^2(\Omega)]^d$ and let $\psi_1, \psi_2\in H^1(\Omega)$ be such that $\text{ div } f_1 = \text{ div } f_2 = 0$ in $\Omega$ and $f_1\cdot \nu = f_2\cdot \nu = 0$ on $\partial \Omega$ and with the property that $f = f_1- \nabla \psi_1 = f_2-\nabla \psi_2$ which gives

Multiplying the latter by $f_1-f_2,$ integrating over $\Omega$ and integrating by parts, one has

It follows that $f_1 = f_2$ and therefore $\psi_1 = \psi_2$ up to an additive constant. To sum up, we have proved that $f$ can be uniquely written as a sum of a divergence free vector field $f_S$ with null normal derivative on the boundary $\partial \Omega$ and a gradient $f_P,$ where an explicit expression for $f_S$ and $f_P$ is provided by the former construction, specifically

with $\psi\in H^1(\Omega)$ unique (up to additive constant) solution of (2.2).

At last, observe that the $L^2$-orthogonality of $f_S$ and $f_P$ and in particular identity (2.1) are immediate consequences of the properties of the two components. This yields the proof.

Remark 2.1. Notice that, as soon as the vector field $f$ is more regular, for instance, say $f\in [H^1(\Omega)]^d$ as it suits our purposes, then the components $f_S$ and $f_P$ of the Helmholtz decomposition are even $H^1$- orthogonal and in particular

Now we are in position to provide the rigorous definition of the so-called Leray projection operator already mentioned in the introduction.

Definition 2.1. Consider the setting of Theorem 2.1. Let $\mathcal{P}$ be the orthogonal projection of $[L^2(\Omega)]^d$ into the subspace of divergence-free vector fields with null normal derivative on $\partial \Omega.$ Then $\mathcal{P}$ is called Leray projection operator. More precisely, for any $f\in [L^2(\Omega)]^d,$ it holds that

where $f = f_S+f_P$ is the Helmholtz decomposition of $f$ as constructed in the previous theorem.

The following result is easily proven.

Lemma 2.1. Consider the setting of Theorem 2.1. The Leray projection operator $\mathcal{P}\colon [L^2(\Omega)]^d\to [L^2(\Omega)]^d$ is bounded.

Proof. By definition, for any $f\in [L^2(\Omega)]^d,$ $\mathcal{P}f: = f-\nabla \psi,$ where $\psi\in H^1(\Omega)$ is the unique solution to the boundary value problem (2.2). In particular, $\psi$ achieves $\int_\Omega |\nabla \psi - f|^2\, dx = \inf_{\phi\in H^1(\Omega)} J(\phi),$ where

Now, from $\inf_{\phi\in H^1(\Omega)} J(\phi)\leq \int_\Omega |f|^2\, dx,$ we immediately get $||\mathcal{P}f||_{[L^2(\Omega)]^d}\leq ||f||_{[L^2(\Omega)]^d},$ which is the thesis.

As showed by the following result, in the specific case $\Omega = \mathbb{R}^d,$ the Leray projection operator has a favorable form in terms of Riesz transform $\mathcal{R} = (\mathcal{R}_1, \mathcal{R}_2, \dots, \mathcal{R}_d)$ defined for any $\phi\in L^2(\mathbb{R}^d),$ in Fourier space, by

Lemma 2.2. Consider the setting of Definition 2.1, fix $\Omega = \mathbb{R}^d$ and let $f\in [L^2(\mathbb{R}^d)]^d$ be any square-integrable vector field on $\mathbb{R}^d.$ Then the $j$-th component of $\mathcal{P}f$ can be written as

for any $j = 1, 2, \dots, d$ and where $\mathcal{R}$ denote the Riesz transform defined in (2.4).

Proof. Let $f\in [L^2(\mathbb{R}^d)]^d.$ It follows from Theorem 2.1 and Definition 2.1 that $\mathcal{P}f = f_S = f-\nabla \psi,$ where $\psi$ satisfies $\Delta \psi = \text{ div } f.$ In Fourier space this yields

In particular, for any $j = 1, 2, \dots, d,$ this gives

and so, by using the Fourier representation of the Riesz transform (2.4),

for any $j = 1, 2, \dots, d.$ From the latter, one immediately gets (2.5) and this concludes the proof.

In passing, observe that since $\psi$ is a solution of (2.2), then it follows easily by elliptic estimates that

(the same would follow from the boundedness of the Leray operator $\mathcal{P}$ (see Lemma 2.1) but in this case we would get a worse bound).

Now, we want to show that as a consequence of the representation (2.6) of $\nabla \psi$ in terms of the Riesz transform, we can get boundedness of type (2.7) replacing the $L^2$-norm with suitable $L^p$-weighted norms.

Notice that, in general, proving boundedness in weighted $L^p$-space does not come as a mere consequence of elliptic estimates as for (2.7) and, in fact, does require a more involved analysis. In our case, it will follow from special well known properties of the Riesz transform that we summarize in the following lemma.

Lemma 2.3 (Boundedness Riesz transform). Let $1 < p < \infty$ and $p'$ such that $1/p+1/p' = 1$ and let $w$ be a weight in the $A_p(\mathbb{R}^d)$-class (see definition (1.12)). Then, for any $j = 1, 2, \dots, d,$ the following bounds on the operator norms of the Riesz transform $\mathcal{R}_j$ hold true:

Moreover, both the bound are sharp, i.e., the best possible bound is established.

Proof. The proof of the sharp bound (2.8) can be found in ^[4] (see also ^[12]), inequality (2.9) can be found in ^[59] (see also ^[23]).

Now we are in position to state and prove the aforementioned boundedness of the operator $\nabla \psi.$ More precisely, we are interested in proving the following result.

Lemma 2.4. Let $u\in \mathcal{S}(\mathbb{R}^d)$ and consider $\psi$ the unique solution to

Then, for any $1 < p < \infty,$ the following estimates hold true:

where $w, c_p, c_{p, d}$ and $r$ are as in Lemma 2.3.

Proof. We will prove only (2.11), the proof of (2.12) is analogous.

Using the representation (2.6), the bound for the Riesz operator (2.8) and the Hölder inequality for discrete measures, we get

This gives (2.11). Bound (2.12) follows in the same way using (2.9).

As a consequence of the previous lemma, we are able to prove the following (almost-) orthogonality result that we will strongly use in the future.

Lemma 2.5. Let $u\in \mathcal{S}(\mathbb{R}^d)$ and consider $u = u_S + u_P$ the Helmholtz decomposition of $u.$ Then, for $1 < p < \infty,$ the following estimates hold true:

where $w, c_p, c_{p, d}$ and $r$ are as in Lemma 2.3.

Proof. The proof is a direct consequence of Lemma 2.4. We know from Theorem 2.1 (see (2.3)) that

with $\psi$ the unique solution to the Poisson problem (2.10). Then, it is easy to see that estimate (2.13) holds true, indeed

where in the last inequality we used (2.11). As it is analogue, we skip the proof of (2.14).

Notice that both Lemma 2.4 and Lemma 2.5 are stated for function in the Schwarz class $\mathcal{S}(\mathbb{R}^d).$ By density, it is straightforward to see that they can be extended to a (weighted) $L^p$ framework.

2.2.   A favorable representation for $(-\Delta^\ast -z)^{-1}$

As already mentioned in the introduction, the usage, as a starting point in our proofs, of an adaptation of the Birman-Schwinger principle to our elasticity context, requires a better understanding of the action of the resolvent operator $(-\Delta^\ast -z)^{-1},$ well-defined for any $z\in \mathbb{C}\setminus [0, \infty),$ associated with the Lamé operator.

The following easy consequence of Helmholtz decomposition will be useful to this end.

Lemma 2.6. Let $d\geq 2$ and let $f$ be a suitably smooth vector field sufficiently rapidly decaying at infinity. Then $-\Delta^\ast$ acts on $f = f_S+f_P$ as

where $f_S$ is a divergence free vector field and $f_P$ a gradient.

Now we are in position to show the validity of identity (1.7), stated in the introduction, which follows as a consequence of (2.15) together with the $H^1$-orthogonality of the components of the Helmholtz decomposition. This is object of the following lemma.

Lemma 2.7. Let $z\in \mathbb{C}\setminus [0, \infty)$ and $g\in [L^2(\mathbb{R}^d)]^d.$ Then the identity

holds true, where $g = g_S+ g_P$ is the Helmholtz decomposition of $g.$

Proof. Given $g\in [L^2(\mathbb{R}^d)]^d,$ we want to obtain an explicit expression of the vector field $f,$ defined by

Observe that, since $z\notin \sigma(-\Delta^\ast) = [0, \infty),$ the previous is equivalent to

Now, writing the Helmholtz decomposition of $f$ and $g,$ namely $f = f_S + f_P$ and $g = g_S + g_P,$ respectively, and using (2.15), the previous can be re-written as

The $H^1$- orthogonality of the two components of the decomposition enables us to split this intertwining equation for both the two components into a system of two decoupled equations, i.e.,

or, equivalently,

Hypothesis (1.3) ensures that $\frac{z}{\mu}, \frac{z}{\lambda + 2\mu}\notin \sigma(-\Delta),$ hence

Plugging these explicit expressions in $f = f_S + f_P,$ from (2.17) we obtain (2.16).

We underline that Lemma 2.7 was already proved in ^[6], its statement and proof were provided also here only for reader's convenience.

2.3.   Uniform resolvent estimates

As already mentioned, in a non self-adjoint framework, the unavailability of a variational characterization of the spectrum causes that Sobolev inequalities no longer suffice to prove spectral bounds. To overcome this lack, uniform resolvent estimates have been recognized as a crucial tool to fruitfully address this problem. For this reason in this subsection we shall prove uniform estimate for the operator $(-\Delta^\ast -z)^{-1}$ that will be the fundamental ingredient in the proof of our main results later on.

Observe that from the representation (2.16), it is reasonable to expect that estimates for $(-\Delta^\ast -z)^{-1}$ should follow as a consequence of the corresponding estimates (if available) for the resolvent operator $(-\Delta-z)^{-1}$ associated to the Laplacian. This is true, indeed, as we will see repeatedly in this paper, the underlying strategy to treat issues concerning the Lamé operator is, actually, to operate at first at the level of the Laplacian taking advantage of the representation (2.16) once the Helmholtz decomposition is operated. On the other hand, this procedure has the main drawback of producing separated outcomes on the single components of the Helmholtz decomposition, which later must be recombined together in order to get meaningful results in application. This need of recombination requires providing suitable (almost-) orthogonality results for the Helmholtz components, roughly speaking some inequality like $||u_S|| + ||u_P||\lesssim ||u||$ (see Lemma 2.5), which turn out to be highly non-trivial to get and require deep result from harmonic analysis and in particular from singular integrals theory (see Lemma 2.3).

The following theorem collects the main estimates for the resolvent operator $(-\Delta-z)^{-1}$ that we will use for our purpose.

Theorem 2.2 (Uniform estimates for $(-\Delta -z)^{-1}$). Let $z\in \mathbb{C}\setminus [0, \infty).$ Then the following estimates for $(-\Delta-z)^{-1}$ hold true.

i) Let $1 < p\leq 6/5$ if $d = 2,$ $2d/(d+2)\leq p \leq 2(d+1)/(d+3)$ if $d\geq 3$ and let $p'$ such that $1/p + 1/p' = 1.$ Then

ii) Let $\alpha > 1/2.$ Then

iii) Let $3/2 < \alpha < 2$ if $d = 2,$ $2d/(d+1) < \alpha \leq 2$ if $d\geq 3$ and let $(d-1)/2(\alpha-1) < p \leq d/\alpha.$ Then for any non-negative function $V$ in $\mathcal{L}^{\alpha, p}(\mathbb{R}^d)$

iv) Let $3/2\leq \alpha < 2$ if $d = 2,$ $d-1\leq \alpha < d$ if $d\geq 3$ and let $\beta = (2\alpha -d+1)/2.$ Then for any non-negative function $V$ such that $|V|^\beta\in \mathcal{KS}_{\alpha}(\mathbb{R}^d)$

Proof. Proof of (2.18) can be found in the work ^[50] by Kenig, Ruiz and Sogge. Estimate (2.19) is proved in the pioneering work by Agmon ^[2] (Lemma 4.1 there), see also ^[63]. Proof of (2.20) was provided by Frank in ^[40] (see also ^[20,22]). Finally, estimate (2.21) is proved by Lee and Seo in ^[53].

Now we are in position to state the corresponding estimate for the resolvent operator $(-\Delta^\ast -z)^{-1}.$

Theorem 2.3. Let $z\in \mathbb{C}\setminus [0, \infty).$ Then, under the same hypotheses of Theorem 2.2, the following estimates for $(-\Delta^\ast-z)^{-1}$ hold true.

If, in addition, $V\in A_2(\mathbb{R}^d),$ then

Remark 2.2. Notice that in order to prove (2.25), that is when the bound on the resolvent operator norm is measured in term of potentials in the Kerman-Saywer space, the additional assumption of $V$ belonging to the $A_2$ class of weights is required, on the contrary this auxiliary hypothesis is not needed when the estimate involves potentials in the Morrey-Campanato class (see (2.24)). This fact is mainly due to a good behavior of functions in the Morrey-Campanato space in relation with the $A_p$ class of weights. This property is clarified in the following result (refer to ^[21], Lemma 1).

Lemma 2.8. Let $V$ be a non-negative function in $\mathcal{L}^{\alpha, p}(\mathbb{R}^d)$ with $0 < \alpha < d$ and $1 < p\leq d/\alpha.$ If $r$ is such that $1 < r < p,$ then $W: = (M V^r)^{1/r}\in A_1(\mathbb{R}^d)\cap \mathcal{L}^{\alpha, p}(\mathbb{R}^d),$ where $M$ denotes the Hardy-Littlewood maximal operator, defined for any $\psi\in L^1_ \rm{loc}(\mathbb{R}^d)$ by $M\psi(x) = \sup_{r > 0}\frac{1}{|B_{r}(0)|}\int_{B_r(0)} |\psi(x-y)|\, dy,$ with $B_r(0)$ the Euclidean ball centered at the origin with radius r. Furthermore there exists a constant $c$ independent of $V$ such that

In passing, notice that $V(x)\leq W(x)$ for almost every $x\in \mathbb{R}^d.$

Now we are in position to prove Theorem 2.3.

Proof of Theorem 2.3. As already mentioned, the proof will basically rely on the interplay between the favorable representation (2.16) of $(-\Delta^\ast -z)^{-1}$ in terms of the resolvent of the Laplace operator together with the uniform estimates for $(-\Delta-z)^{-1}$ summarized in Theorem 2.2 and the orthogonality result stated in Lemma 2.5.

Let us first consider the proof of (2.22).

Let $G$ be any vector-valued function in $[L^p(\mathbb{R}^d)]^d.$ From (2.16) one easily has

where $G = G_S+ G_P$ is the Helmholtz decomposition of $G.$

We shall explicitly estimate the term involving the $S$-component. The analogous term for the $P$- component can be treated similarly. It follows from (2.18) that, for any $j = 1, 2, \dots, d,$

This, along with the sub-additivity of the concave function $|x|^p,$ for $0 < p\leq 1,$ and the Hölder inequality for discrete measures, gives

The same computation for the term involving the $P$- component provides

Plugging the previous two bounds together in (2.26), one has

where $C_{p, d, \lambda, \mu}: = C_{p, d}\, d^\frac{1}{p'} \max\Big\{\mu^{{\frac{d}{2} - \frac{d}{p}}}, (\lambda + 2\mu)^{\frac{d}{2} - \frac{d}{p}} \Big\}.$

Hence, estimate (2.22) follows immediately from the latter as a consequence of (2.13).

Bound (2.23) can be proved with a few modifications to the argument above, namely using (2.19) and (2.14) with $p = 2$ (notice that $\langle x \rangle^{2\alpha}\in A_2(\mathbb{R}^d)$) instead of (2.18) and (2.13), respectively. Similarly bound (2.25) follows as a consequence of (2.21) and (2.14) with $p = 2$.

Now, we turn to the proof of (2.24).

Let $W\in A_1(\mathbb{R}^d)\cap \mathcal{L}^{\alpha, p}(\mathbb{R}^d)$ be as in Lemma 2.8. Observe that, since $W\in \mathcal{L}^{\alpha, p}(\mathbb{R}^d),$ then (2.20) is available. Moreover, as $W\in A_1(\mathbb{R}^d),$ in particular $W\in A_2(\mathbb{R}^d)$ and (2.14) with $p = 2$ is valid too. Therefore, it comes as a slight modification of the argument above proving that the following analogue of (2.25) for $W$ holds true, namely one has

Since $V(x)\leq W(x)$ almost everywhere and as $||W||_{\mathcal{L}^{\alpha, p}(\mathbb{R}^d)}\leq c||V||_{\mathcal{L}^{\alpha, p}(\mathbb{R}^d)},$ using estimate (2.27) one gets

which is (2.24). This concludes the proof of the theorem.

3.   Self-adjoint setting

This section is concerned with the proof of eigenvalue bounds of the form (1.4) for the self-adjoint perturbed Lamé operator. More precisely, we shall prove the following result.

Theorem 3.1. Let $V$ be real-valued and and let $\gamma\geq 1/2$ if $d = 1,$ $\gamma > 0$ if $d = 2$ and $\gamma\geq 0$ if $d = 3.$ Then any negative eigenvalue $z$ of the perturbed Lamé operator $-\Delta^\ast + V$ satisfies

with a constant $C_{\gamma, d, \lambda, \mu}$ independent of $V.$

Here $V_-$ denotes the negative part of $V,$ i.e., $V_-(x): = \max\{-V(x), 0\}.$

Before providing the proof of this theorem, let us comment on the corresponding inequalities of type (3.1) for self-adjoint Schrödinger operators. As the remarks provided below then naturally carry over to Lamé operators, the choice of discussing the case of the Laplacian only, finds its reasons solely in the intent of lightening the discussion.

In the case of $-\Delta + V$ with real-valued potential $V,$ estimate (3.1) was first found by Keller ^[49] in $d = 1$ and, later, generalized to an inequality for the negative eigenvalues power sum known as Lieb-Thirring inequality:

where $\gamma \geq 1/2$ if $d = 1,$ $\gamma > 0$ if $d = 2$ and $\gamma\geq 0$ if $d\geq 3$ (same conditions as in Theorem 3.1) (see ^[57] and ^{[29,54,60,61,64]} for the endpoint cases). In passing, observe that bounds on single eigenvalues, like (3.1), represent a much weaker version of the Lieb-Thirring type inequalities (3.2).

Now, some comment on inequalities (3.1) for $-\Delta + V$ (in fact on the stronger bound (3.2)) are listed below (we refer to ^[56], Chapter 4, for further details).

Remark 3.1. Contrarily to the case of complex-valued potentials, here, as a consequence of the variational principles (no more available in the non self-adjoint context), only the negative part of $V,$ namely $V_-,$ plays a role. Of course, since $-\Delta$ is a non-negative operator, if $V$ is also non-negative then so is $-\Delta + V$ and therefore the variational characterization of the spectrum guarantees that no negative eigenvalues can occur. If $V$ changes its sign, that is if both the positive and negative part of $V = V_+ - V_-$ are non-trivial, it is true that both parts influence the negative eigenvalues, but as $-\Delta + V\geq -\Delta -V_{-},$ it is a consequence of the minimax principle that an upper estimate for the absolute value of negative eigenvalues of $-\Delta - V_-$ provides automatically the same upper estimate for the negative eigenvalues of the complete Hamiltonian $-\Delta + V$ (actually, the same reasoning applies to the eigenvalue power sum). Indeed the effect of $V_+$ on negative eigenvalues is only to increase their size.

Remark 3.2. It is not difficult to see that if $z$ is an eigenvalue of $-\Delta + V$ with eigenfunction $\psi,$ then $\phi_\alpha(\cdot): = \psi(\alpha \, \cdot)$ is an eigenfunction of $-\Delta + V_\alpha(x)$ where $V_{\alpha}(\cdot) = \alpha^2 V(\alpha\, \cdot)$ with eigenvalue $\alpha^2 z.$ By a simple scaling argument, this gives that $p = \gamma + d/2$ is the only possible exponent for which an inequality of the following type

can hold. Notice that here our underlining domain is the entire Euclidean space $\mathbb{R}^d.$

Remark 3.3. Let us underline that there are "natural" constraints on the validity of inequalities of type (3.2) that can be easily justified. We emphasize here the pathological behavior of dimensions $d = 1, 2.$ It is well known that, due to the lack of a Hardy-type inequality, the free Hamiltonian $-\Delta$ is critical in low dimensions, which means that the addition of any arbitrarily small non-trivial negative potential $V$ makes the bottom of the spectrum of the corresponding perturbed operator $-\Delta + V$ negative, thus ensuring existence of negative eigenvalues. On the other hand, if an inequality of the form (3.2) with $\gamma = 0$ holds, then the left-hand side would turn into the counting function of negative eigenvalues and so, as a consequence of the aforementioned criticality, it is an integer greater or equal to one for any such potential. On the contrary, since the right-hand side can be made arbitrarily small, for instance, assuming $L_{0, d} ||V_-||_{L^{d/2}(\mathbb{R}^d)}^{d/2} < 1,$ would give an evident contradiction.

We can now turn to proof of Theorem 3.1. As we will show later, it will come as a consequence of the following lemma which provides dimension-dependent estimates for the expectation value $\langle \psi, V_- \psi\rangle: = \int_{ \mathbb{R}^d} V_-|\psi|^2\, dx$ of the potential energy $V_-$ in the state $\psi\in H^1(\mathbb{R}^d).$

Lemma 3.1. Let $d\geq 1$ and let $V\in L^{\gamma + d/2}(\mathbb{R}^d)$ be a real-valued function. The following estimates for $\langle \psi, V_-\psi \rangle$ hold true.

i) If $d = 1$ and $\gamma \geq \frac{1}{2},$ then

ii) If $d = 2$ and $\gamma > 0,$ then

iii) If $d\geq 3$ and $\gamma \geq 0,$ then

Proof. Let us start with the proof of (3.5).

It is an easy consequence of Hölder inequality that

Being $2\leq 2(2\gamma + d)/(2\gamma + d-2)\leq 2d/(d-2),$ we can use the interpolation inequality to get

Plugging the latter in the former gives (3.5).

Now let us consider $d = 1, 2.$ Estimate (3.4) is immediate consequence of Hölder inequality and the same holds for the case $\gamma = 1/2$ in $d = 1.$ Finally, the remaining case $\gamma > 1/2$ follows, as in the three dimensional framework, from Hölder and interpolation inequality, using that $2\leq 2(2\gamma + 1)/(2\gamma-1) < \infty.$ This proves (3.3) and concludes the proof of the lemma.

In passing, observe that if $\psi\in H^1(\mathbb{R}^d),$ then the norms on the right hand side of (3.3), (3.4) and (3.5) are finite. This is a consequence of the Sobolev embeddings

which are quantified by the inequalities contained in the following lemma (see [55,Ch. 8] or [56,Sec. 2.2.1]).

Lemma 3.2 (Sobolev inequalities). Let $d\geq 1$ and let $\psi\in H^1(\mathbb{R}^d).$

i) If $d = 1,$ then

ii) If $d = 2,$ then

iii) If $d\geq 3,$ then

Here, $S_{2, q}$ and $S_d$ denote the optimal Sobolev constants (^[55,56]).

Now we are in position to prove Theorem 3.1.

3.1.   Proof of Theorem 3.1

In this setting, the variational characterization of the spectrum states that for any $u\in [H^1(\mathbb{R}^d)]^d,$

Therefore, in order to get (3.1) it is sufficient to prove the following lower bound

for any $u\in [H^1(\mathbb{R}^d)]^d$ with $||u||_{[L^2(\mathbb{R}^d)]^d} = 1.$ In order to estimate $\langle u, (-\Delta^\ast - V_-) u \rangle,$ we consider the Helmholtz decomposition $u = u_S + u_P$ of $u.$ It follows from the explicit expression (2.15) of $-\Delta^\ast$ and from the $H^1$-orthogonality of $u_S$ and $u_P$ (see Subsection 2.1) that

Let us start considering the case $d\geq 3.$

Using Sobolev inequality (3.8) on the $j$-th component of the vector-field $u_S$, one gets

where in the last inequality we used the sub-additivity of the concave function $|x|^\frac{d-2}{d}.$

The same computation performed for $u_P$ gives

Now we are in position to estimate $\langle u, V_- u \rangle.$ Using bound (3.5) in Lemma 3.1, two times the Hölder inequality for discrete measures and, finally, the Young inequality, $ab\leq a^p/p + b^q/q$ which holds for all positive $a, b$ and $1/p+1/q = 1,$ we get for some $\varepsilon > 0$

where in the last inequality we simply used the inequality $||u||_{\frac{2d}{d-2}}^2\leq 2 \big[||u_S||_{\frac{2d}{d-2}}^2 + ||u_P||_{\frac{2d}{d-2}}^2 \big].$

Now, plugging (3.11), (3.12) and (3.13) in (3.10), one has

Choosing a suitable small $\varepsilon = : \varepsilon_{\gamma, d, \lambda, \mu},$ one gets (3.9) with

Hence, bound (3.1) is proved if $d\geq 3.$ We skip the proof of the analogous bounds in the lower dimensional cases, namely $d = 1, 2.$ Indeed these follow from the corresponding estimates in Lemma 3.1 and Sobolev inequalities (Lemma 3.2) with minor modifications from the reasoning above.

4.   Non self-adjoint setting: Proof of main results

This section is concerned with the proof of the eigenvalue bounds contained in Theorems 1.2–1.5. As we will see, with the estimates of Theorem 2.3 in hand, the proofs will follow smoothly.

4.1.   Proof of Theorem 1.2

As already mentioned in the introduction, the starting point in our proofs is the Birman-Schwinger principle. In our context it states that if $z\in \mathbb{C}\setminus [0, \infty)$ is an eigenvalue of $-\Delta^\ast + V,$ then $-1$ is an eigenvalue of the Birman-Schwinger operator $V_{1/2} (-\Delta^\ast - z)^{-1} |V|^{1/2}$ on $[L^2(\mathbb{R}^d)]^d.$ This implies that the operator norm of the latter is at least $1.$ Therefore, in order to get the bound (1.8), we are reduced to prove

The same strategy, with the needed modifications, then will be also applied to prove the corresponding bounds in Theorem 1.3, Theorem 1.4 and 1.5.

Providing bounds for $||V_\frac{1}{2} (-\Delta^\ast -z)^{-1} |V|^\frac{1}{2}||_{L^2\to L^2}$ is estimating the quantity $|\langle f, V_{1/2} (-\Delta^\ast - z)^{-1}|V|^\frac{1}{2} g \rangle|,$ for $f, g\in [L^2(\mathbb{R}^d)]^d.$ To simplify the notation we introduce the function $G: = |V|^\frac{1}{2} g.$ Using the Hölder inequality and estimate (2.22) for the resolvent $(-\Delta^\ast -z)^{-1},$ we have

Thus, recalling that $G = |V|^{1/2} g,$ one has

Using the Hölder inequality and its version for discrete measures and the sub-additivity property of the concave function $|x|^p,$ for $0 < p\leq 1,$ one gets

The same estimate for the term involving $g$ gives

Plugging these two bounds in (4.2) we end up with the following inequality

Now, choosing $p = \frac{2(2\gamma + d)}{2\gamma + d + 2}$ (observe that the restriction on $\gamma$ in Theorem 1.2 guarantees that $p$ satisfies the hypotheses in Theorem 2.3) and taking the supremum over all $f, g\in [L^2(\mathbb{R}^d)]^d$ with norm less than or equal to one, we get (4.1). This concludes the the proof of Theorem 1.2.

4.2.   Proof of Theorem 1.3

As in the proof of the previous result, we are reduced to prove the following bound

We pick a strictly positive function $\phi \in \mathcal{L}^{\alpha, p}$ and we define a strictly positive approximation of our potential, that is $V_\varepsilon(x): = \sup\{ |V(x)|, \varepsilon \phi(x)\}.$ Using Cauchy-Schwarz inequality and estimate (2.24) for the resolvent $(-\Delta^\ast -z)^{-1},$ we have

Thus, the theorem is proved once $\varepsilon$ goes to zero, taking the supremum over all $f, g \in [L^2(\mathbb{R}^d)]^d$ with norm less than or equal to one and by choosing $\alpha = \frac{2d}{2\gamma +d}.$

4.3.   Proof of Theorem 1.4

Again we are reduced to prove the following bound

The same strategy used above, with the usage of (2.25) instead of (2.24), gives

Thus, the theorem is proved once $\varepsilon$ goes to zero, taking the supremum over all $f, g \in [L^2(\mathbb{R}^d)]^d$ with norm less than or equal to one and by choosing $\alpha = \frac{d(d-1)}{d-2\gamma}.$

4.4.   Proof of Theorem 1.5

As before, our problem is reduced to proving the following bound

with $q = 2\gamma + (d-1)/2.$ First, observe that from (2.22), with the choice $p = 2(d+1)/(d+3)$ one has

In passing, notice that (1.8) was obtained with the usage of (2.18) with the choice $p = 2(2\gamma + d)/(2\gamma + d+2)$ (see last part of Proof of Theorem 1.2). Then $p = 2(d+1)/(d+3)$ corresponds to the case $\gamma = 1/2$ in (1.8) which gave the aforementioned decay threshold $2d/(d+1)$.

We know from (2.23) that

Using Riesz-Thorin interpolation between estimate (4.4) and (4.5), we get

with $p = 2(d+1)/(d+3),$ $\alpha > 1/2$ and $1/p_\theta + 1/p_\theta' = 1.$

From this fact it is easy to get

Taking the supremum over all $f, g\in [L^2(\mathbb{R}^d)]^d$ with norm less than or equal to one and raising the resulting inequality to the power $p_\theta/(2-p_\theta)$ gives

Here we abuse the notation by using the same symbol $C_{\alpha, d, \lambda, \mu}$ for different constants.

Calling

this clearly gives $\frac{p_\theta}{2-p_\theta} = 2\gamma\, \frac{d+1}{2-\theta + d\theta},$ since we also have $\frac{1}{p_\theta} = \frac{1-\theta}{p} + \frac{\theta}{2},$ with $p = \frac{2(d+1)}{d+3},$ this leads to the constraint $\theta = 1-\frac{d+1}{4\gamma + d-1}.$ With these choices one has

with $q = 2\gamma + (d-1)/2,$ which is the bound (4.3) once renaming $\alpha(2\gamma -1) = \alpha.$

Acknowledgments

The author is deeply grateful to Prof. L. Fanelli for bringing the attention to the setting of the problem and for valuable comments on the draft of the manuscript. The author also thanks Prof. D. Krejčiřík for useful discussions which greatly enrich the paper. The author is greatful to Prof. M. Correggi for helpful references' suggestions. The author thanks the hospitality of Sapienza University of Rome where this work was initiated and Czech Technical University in Prague where it was developed in part. The author gratefully acknowledges financial support by the Deutsche Forschungsgemeinschaft (DFG) through CRC 1173.

Conflict of interest

The author declares no conflict of interest.

A.   Rigorous definition of $-\Delta^\ast + V$

In this appendix we give the rigorous definition of the operator $-\Delta^\ast + V.$ This is introduced as an $m$-accretive operator obtained as a form sum of the free Lamé operator $-\Delta^\ast$ with domain $[H^1(\mathbb{R}^d)]^d$ and a relatively form-bounded potential $VI_{ \mathbb{C}^d\times \mathbb{C}^d}$ with $V\colon \mathbb{R}^d \to \mathbb{C}.$ More precisely, our standing assumption is that there exist numbers $a\in (0, 1)$ and $b\in \mathbb{R}$ such that, for all $u \in [H^1(\mathbb{R}^d)]^d,$

where

Here $\nabla u$ denotes the matrix with entries the first derivatives of $u,$ namely $(\nabla u)_{jk} = \frac{\partial u_j}{\partial x_k},$ whereas $\nabla^T u$ is the transpose matrix of $\nabla u,$ namely $(\nabla^T u)_{jk} = \frac{\partial u_k}{\partial x_j}.$ We refer the reader to [24,Sec. 2.1] for more details. In the following we will show that under the hypotheses of Thm. 1.1–Thm. 1.5, $V$ satisfies (A.1).

As a warm up we consider the easiest case $d = 1$ first. If $u\in H^1(\mathbb{R}),$ then by Sobolev embedding $u\in L^\infty(\mathbb{R}).$ In particular (3.6) holds. Using (3.6), the 1D expression of the Lamé operator (1.5) and the Cauchy Schwarz inequality one easily has that for any $\varepsilon > 0$

Since $\varepsilon$ is arbitrary one gets that if $V\in L^1(\mathbb{R}),$ then (A.1) holds and $-\Delta^\ast + V$ is well defined as an $m$-accretive operator on $H^1(\mathbb{R}).$

In order to tackle the higher dimensional framework $d\geq 2$ we need to state a preliminary lemma. For simplicity, in the following we shorten the notation for the $L^2-L^2$ operator norm $\|\cdot\|_{L^2(\mathbb{R}^d)\to L^2(\mathbb{R}^d)}$ to the more compact $\|\cdot\|.$

Lemma A.1. Let $d\geq 2$ and $z\in \mathbb{C}\setminus [0, \infty).$ For $\gamma$ as in Theorem 1.2, if $V\in L^{\gamma+ \frac{d}{2}}(\mathbb{R}^d),$ then

For $\gamma,$ $\alpha$ and $p$ as in Theorem 1.3, if $V\in \mathcal{L}^{\alpha, p}(\mathbb{R}^d),$ then

For $\gamma,$ $\alpha$ and $\beta$ as in Theorem 1.4, if $V^\beta\in \mathcal{KS}_{\alpha}(\mathbb{R}^d)$ and if $V\in A_2(\mathbb{R}^d),$ then

For $\gamma,$ $\alpha$ and $q$ as in Theorem 1.5, if $V\in L^q(\langle x \rangle^{2\alpha} dx),$ then

The constant $C$ may change from (A.2) to (A.5) but in each estimate it is independent of $V.$

Proof. The proof of estimates (A.2)–(A.5) follows with almost no modification of the proofs of Theorem 1.2–Theorem 1.5, respectively (see Sec. 4).

With Lemma A.1 at hand we see that $V$ satisfies (A.1) whenever the hypotheses of Thm. 1.2–Thm. 1.5 are satisfied. We consider first the case $\gamma\neq 0.$ The special case $\gamma = 0$ will be treated separately.

Let $z\in (-\infty, 0).$ We write estimates (A.2)–(A.5) in a unified way as follows:

where $C$ might depend on $V$ and where $\|V\|$ denotes one of the norms of the right hand side of estimates (A.2)–(A.5). Defining $A: = |V|^{1/2}(-\Delta^\ast -z)^{-1/2}$ one has that $A^\ast = (-\Delta^\ast- z)^{-1/2}|V|^{1/2}$ and consequently $AA^\ast = |V|^{1/2}(-\Delta^\ast-z)^{-1} |V|^{1/2}.$ Using that $\|AA^\ast\| = \|A^\ast\|^2 = \|A\|^2$ we get

Thus, bound (A.1) is obtained in the limit $|z|\to \infty.$

Now it is left the case $\gamma = 0.$ We shall see the explicit computations only in the Lebesgue setting, namely when $V\in L^\frac{d}{2}(\mathbb{R}^d),$ when $V$ belongs to the Morrey-Campanato, Keller-Sawyer and weighted $L^q$ class one can argue similarly.

Let $R > 0,$ we split the left hand side of (A.1) as follows

For $I_{\leq R}$ one easily has

As for $I_{ > R},$ using Hölder and Sobolev inequality (3.8) we have

where in the last step we have used that $\|\nabla u\|_{L^2(\mathbb{R}^d)}^2\leq Ch_0^\ast[u]$ (see [24,Sec. 2.1]).

Plugging the previous estimates for $I_{ > R}$ and $I_{\leq R}$ in (A.6) one eventually has

Since the measure of the set $\{x\in \mathbb{R}^d\colon |V(x)| > R\}$ tends to zero as $R$ goes to infinity, one has that (A.1) follows from (A.7).