Eigenvalues of one-dimensional non-self-adjoint Dirac operators and applications

TL;DR

This paper studies the eigenvalues of one-dimensional non-self-adjoint Dirac operators with complex potentials, providing new asymptotic results, inequalities, and applications to physics such as wave damping and graphene nanoribbons.

Contribution

It introduces new analysis of eigenvalues for non-self-adjoint Dirac operators, including existence, asymptotics, and inequalities, with applications to physical models.

Findings

Existence and asymptotics of weakly coupled eigenvalues

Lieb-Thirring inequalities for non-self-adjoint operators

Applications to damped wave equations and graphene nanoribbons

Abstract

We analyze eigenvalues emerging from thresholds of the essential spectrum of one-dimensional Dirac operators perturbed by complex and non-symmetric potentials. In the general non-self-adjoint setting we establish the existence and asymptotics of weakly coupled eigenvalues and Lieb-Thirring inequalities. As physical applications we investigate the damped wave equation and armchair graphene nanoribbons.

Equations

$$H=-\mathrm{i}{\partial }_x{\sigma }_1+m{\sigma }_3=\left(\begin{array}{cc}m& -\mathrm{i}{\partial }_x\\ -\mathrm{i}{\partial }_x& -m\end{array}\right)$$

$$\sigma (H)={\sigma }_{e3}(H)=(-\infty ,-m]\cup [m,\infty );$$

$${\sigma }_{e3}(T):=\{z\in \mathbb{C}:T-z\text{ is not Fredholm}\},$$

$${\sigma }_1=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right),{\sigma }_2=\left(\begin{array}{cc}0& -\mathrm{i}\\ \mathrm{i}& 0\end{array}\right),{\sigma }_3=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right).$$

$$V:\mathbb{R}\to {\mathbb{C}}^{2\times 2},\Vert V\Vert \in L^1(\mathbb{R})\cap L^2(\mathbb{R}),$$

$$\Vert V{\Vert }_p^p:={\int }_{\mathbb{R}}\Vert V(x){\Vert }^p\mathrm{d}x.$$

$${\sigma }_{\mathrm{p}}(H+V)\setminus {\sigma }_{e3}(H)\subset {\overline{B}}_{m{r}_0}(m{x}_0)\dot{\cup }{\overline{B}}_{m{r}_0}(-m{x}_0)$$

$$x_0:={\left(\frac{\Vert V{\Vert }_1^4-2\Vert V{\Vert }_1^2+2}{4(1-\Vert V{\Vert }_1^2)}+\frac{1}{2}\right)}^{\frac{1}{2}},r_0:={\left(\frac{\Vert V{\Vert }_1^4-2\Vert V{\Vert }_1^2+2}{4(1-\Vert V{\Vert }_1^2)}-\frac{1}{2}\right)}^{\frac{1}{2}}.$$

$${\int }_{\mathbb{R}}\left(V(x{)}^2+m\{{\sigma }_3,V(x)\}\right)\mathrm{d}x,$$

$$\mathcal{R}(x,y;z)$$

$$\mathcal{N}(x,y;z)$$

$$\zeta (z):=\frac{z+m}{k(z)},k(z):=(z^2-m^2{)}^{\frac{1}{2}}$$

$$\mathcal{Q}(x,y;z)=\mathcal{A}(x)\mathcal{N}(x,y;z)e^{\mathrm{i}k(z)|x-y|}\mathcal{B}(y)$$

$$V=\mathcal{B}\mathcal{A},\mathcal{B}:=U_V|V{|}^{\frac{1}{2}},\mathcal{A}:=|V{|}^{\frac{1}{2}}.$$

$$\Vert \mathcal{A}(x){\Vert }^2=\Vert \mathcal{B}(x){\Vert }^2=\Vert V(x)\Vert .$$

$$Q=L+M;$$

$$\mathcal{L}(x,y)$$

$$\mathcal{M}(x,y)$$

$${\rm Y}={\rm Y}(\zeta (z))=\frac{\mathrm{i}}{2}\left(\begin{array}{cc}\zeta (z)& 0\\ 0& \zeta (z{)}^{-1}\end{array}\right).$$

$$\Vert L\Vert \le \Vert {\rm Y}\Vert \Vert V{\Vert }_1=\frac{1}{2}\max \{|\zeta (z)|,|\zeta (z){|}^{-1}\}\Vert V{\Vert }_1.$$

$$\Vert Q{\Vert }_{{\mathfrak{S}}^2}^2\le \left(\frac{1}{2}+\Vert {\rm Y}{\Vert }_{{\mathfrak{S}}^2}^2\right)\Vert V{\Vert }_1^2=\frac{1}{4}\left(2+|\zeta (z){|}^2+|\zeta (z){|}^{-2}\right)\Vert V{\Vert }_1^2.$$

$$\|M\|=o\big{(}\|\Upsilon\|\big{)},\quad z\to\pm m,\ z\notin\sigma(H).$$

$$\Vert {\rm Y}{\Vert }^{-2}\Vert \mathcal{M}(x,y){\Vert }^2\le C\Vert V(x)\Vert \Vert V(y)\Vert \left(\Vert {\rm Y}{\Vert }^{-2}+1\right),z\notin \sigma (H).$$

$$\underset{k(z)\to 0}{\lim }(e^{\mathrm{i}k(z)|x-y|}-1)=0,$$

$$\underset{k(z)\to 0}{\lim }\Vert {\rm Y}{\Vert }^{-2}{\int }_{{\mathbb{R}}^2}\Vert \mathcal{M}(x,y){\Vert }^2\mathrm{d}x\mathrm{d}y=0.\text{qed}$$

$$U:={\int }_{\mathbb{R}}V(x)\mathrm{d}x.$$

$$z_{+}(\epsilon )=m-\frac{m}{2}{U}_{11}^2{\epsilon }^2+o({\epsilon }^2),\epsilon \to 0+.$$

$$z_{-}(\epsilon )=-m+\frac{m}{2}{U}_{22}^2{\epsilon }^2+o({\epsilon }^2),\epsilon \to 0+.$$

$$I+\epsilon Q(z)=(I+\epsilon M)(I+(I+\epsilon M{)}^{-1}\epsilon L),$$

$$\det (I_2+(I+\epsilon M{)}^{-1}\epsilon L)=0.$$

$$(I+\epsilon M{)}^{-1}L={\hat{A}}_{\epsilon }{\rm Y}\hat{B},$$

Full text

The research of P.S. is supported by the Swiss National Science Foundation, SNF Ambizione grant No. PZ00P2_154786.

1. Introduction

Dirac operators attracted considerable attention in recent years, in particular in the context of non-self-adjoint spectral theory [4, 5, 6, 9, 12, 24, 25], nonlinear Schrödinger equations e.g. [3, 21] or as an effective model for graphene [1, 7, 15, 20]. In this paper we analyze eigenvalues emerging from the thresholds of the essential spectrum of the one-dimensional Dirac operator in $L^{2}({\mathbb{R}})$ perturbed by a general matrix-valued and non-symmetric potential $V$ preserving the essential spectrum.

Our main results include the existence and asymptotics of weakly coupled eigenvalues for the one dimensional Dirac operator (Theorem 2.2) and Lieb-Thirring type inequalities (Theorem 2.4) in the massive as well as the massless case. These results complement the eigenvalue estimates in [6] and also show that the latter are optimal in the weak coupling regime, see Remark 2.3.

As physical applications, we investigate the damped wave equation in $L^{2}({\mathbb{R}})$ and to a two-dimensional model of charge carriers in graphene nanoribbons (or waveguides) with so-called armchair boundary conditions. We emphasize here the inherent non-self-adjoint nature of the former caused by the presence of damping. Moreover, our eigenvalue estimates may be converted to resonance estimates via the well-known method of complex scaling, as in [6].

The application for the damped wave equation (Theorem 3.1) demonstrates a natural effect from the physical point of view: The integrable damping $\varepsilon a_{1}(x)$ cannot affect the essential spectrum; however, for any $\varepsilon>0$, it gives rise to a pair of complex conjugated eigenvalues having the tendency to meet at the real axis. The interpretation of the results for the graphene armchair waveguides is more complicated due to the $4\times 4$ matrix structure and the PDE nature of the problem. Nonetheless, in the simplest setting of a diagonal potential that is constant in the transverse direction, the quantities entering the eigenvalue asymptotics are expressed in terms of the integral of the trace of $V$ only, see Example 3.7 and Theorem 3.6.

The main ingredient in the proofs is the analysis of a Birman-Schwinger operator. Since the problem is not self-adjoint, the existence of eigenvalues in the gap of the essential spectrum does not follow from min-max considerations. Nonetheless, the weak coupling technique, relying on the isolation of a singular part $L$ of the Birman-Schwinger operator, admits a generalization to the non-self-adjoint setting. This is possible since $L$ is of finite rank and so the question of existence and asymptotics of eigenvalues is converted to a matrix problem, which is analyzed with the help of Rouché’s theorem eventually. The proofs of the Lieb-Thirring type bounds are also based on complex analysis techniques, this time on a generalization of Jensen’s identity due to [2]. Inequalities of this type were established in [10] for one- and multidimensional Dirac operators, and improved results in the multidimensional case recently appeared in [5]. The difference of our new estimates compared to the one-dimensional results in [10] is that the weights in the eigenvalue sums are better, which leads to tighter upper bounds for the number of eigenvalues in certain subsets of the complex plane. The price to pay for this improvement is that the eigenvalue sums cannot be controlled by a single $L^{p}$ norm, but only by a combination of two such norms. This phenomenon was already encountered in [5], and the reason for it is a lack of decay of the free resolvent as the spectral parameter tends to infinity.

To avoid technicalities related to domain questions we intentionally require that $V$ is both integrable and square-integrable throughout the entire paper. The technique allowing the omission of the $L^{2}$ assumption is described in [6, Sec. 6]. In the waveguide case, the $L^{2}$ assumption is also convenient, though not essential, when estimating the infinite sums arising from the decomposition of the resolvent, see Remark 3.3. To simplify the presentation of the weak coupling eigenvalue asymptotics (2.25), (2.26) and (3.45), (3.46), we also do not strive for higher order terms in the expansion, although these could in principle be obtained in a similar way as in the Schrödinger case, see e.g. [26, 22]. We have now also all needed ingredients in hand to prove an analogue of the Lieb-Thirring inequalities in Subsection 2.2 for the graphene waveguides. However, the conformal map would be much more involved and we do not pursue this direction.

The paper is organized as follows. In Section 2 we briefly recall the relevant results of [6] for the one dimensional non-self-adjoint Dirac operator and establish the weak coupling eigenvalue asymptotics (Subsection 2.1), and the Lieb-Thirring inequalities (Subsection 2.2). In Section 3, we apply these results to the one-dimensional damped wave equation, (Subsection 3.1), and graphene waveguides (Subsection 3.2).

2. One dimensional Dirac operator

The spectrum of the free operator $H$ with $m\geq 0$ in $L^{2}({\mathbb{R}};{\mathbb{C}}^{2})$

[TABLE]

reads

[TABLE]

here and in the sequel, we use the essential spectrum $\sigma_{\rm e3}$ defined as

[TABLE]

see e.g. [11, Sec. IX] for details, and the Pauli matrices

[TABLE]

The $n\times n$ identity matrix is denoted by $I_{n}$.

As a perturbation, we consider a (possibly complex and non-symmetric) matrix potential

[TABLE]

where $\|V(x)\|$ is the operator norm in ${\mathbb{C}}^{2}$ of the matrix $V(x)$.

Our goal is to analyze the spectrum of $H+V$. Note that the $L^{2}$ condition in (2.5) is imposed for technical reasons only and is not strictly necessary for the results of this section. It offers the advantage that $H+V$ may be defined as an operator sum because $V$ is relatively $H$-compact and hence infinitesimally $H$-bounded. In addition, the relative compactness implies that the essential spectrum is stable, i.e. $\sigma_{\rm e3}(H+V)=\sigma_{\rm e3}(H)$. If only the $L^{1}$ condition is assumed in (2.5), the perturbed operator can be defined by means of a resolvent formula; we refer to [6] for the details.

For any $p\in[1,\infty)$ we set

[TABLE]

It is proved in [6, Thm. 2.1] that if $\|V\|_{1}<1$, then all non-embedded eigenvalues of $H+V$ satisfy

[TABLE]

where

[TABLE]

2.1. Weakly coupled eigenvalues

Here we analyze the point spectrum of $H+\varepsilon V$ as $\varepsilon\to 0+$, i.e. the weak coupling regime. In the self-adjoint setting, a straightforward construction of test functions together with a min-max argument applied to $(H+V)^{2}$ shows that if the matrix

[TABLE]

where the brackets $\{\cdot,\cdot\}$ denote the anticommutator, has a negative eigenvalue, then $H+V$ has an eigenvalue in $(-m,m)$. In the weak coupling regime, (2.9) can be translated (by ignoring the term of $\mathcal{O}(\varepsilon^{2})$) to {\color[rgb]{0,0,0}\int_{{\mathbb{R}}}V_{11}<0} or {\color[rgb]{0,0,0}\int_{{\mathbb{R}}}V_{22}>0}. In Theorem 2.2 below, we prove that the intuition obtained from this simple self-adjoint argument is indeed correct.

The free resolvent $(H-z)^{-1}$, $z\in\rho(H)$, is an integral operator with the kernel (see [6] for details)

[TABLE]

where

[TABLE]

and the square root on ${\mathbb{C}}\setminus[0,\infty)$ is chosen such that $\operatorname{Im}k(z)>0$.

A natural technique is the Birman-Schwinger principle, derived for the Dirac operator e.g. in [6, Thm. 6.1]. The Birman-Schwinger operator $Q(z)$ is an integral operator with the kernel

[TABLE]

where the factorization of $V$ is based on its polar decomposition $V=U_{V}|V|$, namely

[TABLE]

Notice that

[TABLE]

As usual, we split $Q=Q(z)$ into a singular and a regular part $L$ and $M$, respectively, i.e.

[TABLE]

the corresponding ($z$-dependent) kernels read

[TABLE]

where

[TABLE]

Similarly as in [6, Sec. 2], estimating the quadratic form of $L$, we obtain the bound

[TABLE]

For later use, we also notice that

[TABLE]

The following lemma shows that the possible singularities for $z=\pm m$ of the regular part $M$ are weaker than those of $L$.

Lemma 2.1**.**

Let $V$ be as in (2.5), $M$ be the integral operator with the kernel (2.17) and $\Upsilon$ as in (2.18). Then

[TABLE]

Proof.

The proof is inspired by [22, Lem. 1]. Since $\operatorname{Im}k(z)>0$ for $z\notin\sigma(H)$, straightforward estimates and (2.14) show that there exists $C>0$ such that

[TABLE]

Thus, for $\|\Upsilon\|>1$, the function $(x,y)\mapsto\|\Upsilon\|^{-2}\|\mathcal{M}(x,y)\|^{2}$ has an integrable upper bound. Since $\|\Upsilon\|\to\infty$ when $k(z)\to 0$ and

[TABLE]

the dominated convergence theorem yields

[TABLE]

Theorem 2.2**.**

Let $H$ be as in (2.1), $V$ as in (2.5) and let $m>0$. Define the matrix

[TABLE]

If $\operatorname{Re}U_{11}<0$, then, for all sufficiently small $\varepsilon>0$, there exists an eigenvalue $z_{+}(\varepsilon)$ of $H+\varepsilon V$ satisfying

[TABLE]

Similarly, if $\operatorname{Re}U_{22}>0$ then, for all sufficiently small $\varepsilon>0$, there exists an eigenvalue $z_{-}(\varepsilon)$ of $H+\varepsilon V$ satisfying

[TABLE]

Proof.

According to the Birman-Schwinger principle, see [6, Thm. 6.1], $z\in\sigma_{\rm p}(H+\varepsilon V)$ if and only if $-1\in\sigma_{\rm p}(\varepsilon Q(z))$, where $Q$ is as in (2.15). Thus we investigate the invertibility of $1+\varepsilon Q(z)$. Notice that, if $\varepsilon\|M\|<1$, we have

[TABLE]

so the invertibility of $Q(z)$ depends on the invertibility of the second factor in (2.27). We proceed with the analysis of the latter, find for which $z$ it is not invertible and show that, for these $z$, the condition $\varepsilon\|M\|<1$ holds.

The crucial observation is that the kernel of $L$ is separated (in $x$ and $y$), hence $(I+\varepsilon M)^{-1}\varepsilon L$ is of finite rank and so $-1\in\sigma_{\rm p}((I+\varepsilon M)^{-1}\varepsilon L)$ if and only if

[TABLE]

To analyze (2.28), it is convenient to write

[TABLE]

where

[TABLE]

where $\mathcal{A}_{\varepsilon}(x)$ is a matrix satisfying

[TABLE]

Notice also that (with $U$ as in (2.24))

[TABLE]

Employing these observations and Sylvester’s determinant identity, we can rewrite (2.28) as

[TABLE]

Since

[TABLE]

our initial guess (ignoring the smaller terms) for the dependence of $\zeta$ on $\varepsilon$ reads

[TABLE]

Notice that the assumption that $\operatorname{Re}U_{11}<0$ or $\operatorname{Re}U_{22}>0$ is needed since the allowed region in terms of $\zeta$ is $\operatorname{Im}\zeta<0$, see the discussion of [4] after (2.6) there.

Finally, we prove that there is indeed a solution $\zeta_{+}$ of (2.28) in a neighborhood of $\zeta^{0}_{+}$; the reasoning for $\zeta_{-}$ is analogous. To this end, for $|\alpha|<\delta_{0}$, we define $\zeta_{\alpha}:=\zeta^{0}_{+}(1+\alpha)$ and select $\delta_{0}$ so small that, with some $\beta>0$, we have $\operatorname{Im}\zeta_{\alpha}<-\beta<0$ for all $|\alpha|<\delta_{0}$ and $\varepsilon\to 0+$. Observe that with this $\zeta_{\alpha}$, we have (uniformly in $\alpha$) that $\varepsilon\|M\|=o(\varepsilon|\zeta_{+}^{0}|)=o(1)$ as $\varepsilon\to 0+$ and so $r=o(1)$ as $\varepsilon\to 0+$; the former justifies (2.27) in particular.

Take $\zeta=\zeta_{\alpha}$ and define the function

[TABLE]

Since, for all $|\alpha|<\delta_{0}$, we have $\operatorname{Im}\zeta_{\alpha}<-\beta<0$ and $\zeta_{\alpha}\to\infty$ as $\varepsilon\to 0+$, i.e. the corresponding $z_{\alpha}\in\rho(H)$ and $z_{\alpha}\to m$ as $\varepsilon\to 0+$, the function $F$ is holomorphic for $|\alpha|<\delta_{0}$. Moreover, using (2.31), (2.32), we get

[TABLE]

Hence, Rouché’s theorem implies that, for all $0<\varepsilon<\varepsilon_{\delta_{0}}$, functions $F(\alpha)$ and $G(\alpha):=-\alpha$ have the same number of zeros in the ball $B_{\delta_{0}}(0)$. Notice that the same reasoning is valid for any $0<\delta<\delta_{0}$, thus we obtain that the sought solution of (2.28) reads

[TABLE]

The last step, yielding (2.25), is to use the relation (2.11) between $\zeta$ and $z$ and rewrite (2.38) in terms of $z_{+}$. ∎

Remark 2.3**.**

1. \edefnn\selectfonti)

Theorem 2.2 shows that the spectral estimate (2.7) is sharp in the weak coupling regime. Indeed, the latter can be stated as

[TABLE] 2. \edefnn\selectfonti)

Notice that in the proof of Theorem 2.2, we use that $\|M\|=o(\|\Upsilon(\zeta(z))\|)$ as $z\to\pm m$, $z\notin\sigma(H)$ only and not the particular structure of $M$. The latter would be needed to derive more terms in the expansions (2.25), (2.26). 3. \edefnn\selectfonti)

It is known that the weak coupling limit for the Dirac operator is equivalent to the non-relativistic limit. We do not pursue this connection here and refer to [4] for a discussion in this matter.

2.2. Lieb-Thirring inequalities

In the last subsection we have seen that the massive Dirac operator is critical, i.e. an arbitrarily small (non-self-adjoint) perturbation will create an eigenvalue. Here we prove an upper bound for the number of eigenvalues in certain subsets of the complex plane. The upper bound will be a consequence of a Lieb-Thirring type inequality. We prove similar results for the (non-critical) massless Dirac operator.

Theorem 2.4**.**

Let $H$ be as in (2.1) and $V$ as in (2.5). If $m=0$ and $\|V\|_{1}\geq 1$, then

[TABLE]

where each eigenvalue is counted according to its algebraic multiplicity. If $m>0$, then for any $\tau>0$ we have that

[TABLE]

where

[TABLE]

and where $\rho_{0}$ is the unique solution to $\rho e^{\frac{1}{2}(\rho+1)^{2}}=1$ ($\rho_{0}\approx 0.38$).

Remark 2.5**.**

1. \edefnn\selectfonti)

We recall that in the massless case ($m=0$) the spectral inclusion (2.7) states that there are no eigenvalues whenever $\|V\|_{1}<1$. For this reason we are assuming that $\|V\|_{1}\geq 1$ above. 2. \edefnn\selectfonti)

We emphasize the dependence of the bound on the $L^{1}$ norm of $V$. One reason is that this norm is invariant with respect to rescaling of the mass, i.e. the substitution $V\to m^{-1}V(m^{-1}\cdot)$ does not change the $L^{1}$ norm. Secondly, a straightforward adaptation of the proof shows that if instead of $\|V\|\in L^{2}({\mathbb{R}})$ we assume that $\|V\|\in L^{p}({\mathbb{R}})$ for some $p>1$, then (2.40) and (2.41) hold with $\|V\|_{2}^{2}$ replaced by $\|V\|_{p}^{2p/(p-1)}$. 3. \edefnn\selectfonti)

The bounds of Theorem 2.4 should be compared with those of [10]. In the one-dimensional case it is claimed that, for $p>1$ and e.g. for $m=0$, one has

[TABLE]

for any $\tau\in(0,\min(p-1,1))$. In fact, an inspection of the proof shows that the constant $C_{p,\tau}$ still depends on the potential through the parameter $b$ introduced in Section 4.1 there. Our estimate (2.41) yields better weights both for sequences of eigenvalues accumulating to a point in the essential spectrum or tending to infinity. Moreover, the constant is universal; the dependence on (the $L^{2}$ norm of) the potential is exhibited explicitly. The comparison in the massive case is less obvious, and we will not pursue the issue here. 4. iv)

In the self-adjoint case, Lieb-Thirring type inequalities for one-dimensional Dirac operators may be found in [14]. Note that even there the eigenvalue sums cannot be controlled by a single $L^{p}$ norm. Similar inequalities for resonances were established in [18] and [23] in the massless case and in [19] for the massive case.

Theorem 2.4 has the following consequences for the number of eigenvalues $N_{K}(H+V)$ of $H+V$ in a compact subset $K\subset\rho(H)$. For any $\delta,R>0$ and $\epsilon\geq 0$ we set

[TABLE]

Corollary 2.6**.**

If $m=0$ and $\|V\|_{1}\geq 1$, then we have

[TABLE]

If $m>0$, then we have

[TABLE]

Proof.

For $m>0$ the claim follows from the lower bound

[TABLE]

This can be seen by treating the cases $|z|\leq 2m$ and $|z|>2m$ separately and using that $|z^{2}-m^{2}|\geq m\epsilon$ in the first case. The case $m=0$ is even easier. ∎

Proof of Theorem 2.4.

The proof is based on complex analysis. The general approach is explained very well in [8] and we refer the reader to this article for more details. We treat the massive and the massless case separately, starting with the latter. For simplicity we prove (2.40) only for eigenvalues in the upper half plane ${\mathbb{C}}^{+}$; the proof for the lower half plane is analogous. The basic idea in the complex analysis approach is to relate the eigenvalues to the zeros of a holomorphic function. It is convenient to define the following maps (recall that $Q=|V|^{\frac{1}{2}}R(z)V^{\frac{1}{2}}$):

[TABLE]

where ${\rm i}\eta\in\rho(H+V)$ will be chosen momentarily. Note that $\varphi,\nu,\psi$ are conformal maps and the regularized determinant $\det_{2}$ is defined for any $T\in\mathfrak{S}^{2}$ by

[TABLE]

We then have that $\det_{2}(I+T)=0$ if and only if $(I+T)$ is not invertible; see e.g. [27, Thm. 9.2]. In particular, $h(z)=0$ if and only if $z\in\sigma_{\rm p}(H+V)$. Moreover, $h$ is analytic in the Hilbert-Schmidt norm, see e.g. [13]. From [27, Thm. 9.2] and (2.20) we have the uniform bound (notice the $\|\Upsilon\|_{\mathfrak{S}^{2}}^{2}=1/2$ for $m=0$)

[TABLE]

The value of the optimal constant is $\Gamma_{2}=1/2$, see e.g. formula (2.2) in Chaper IV of [17]. This implies the following estimate,

[TABLE]

We will use the $L^{2}$ norm to estimate the second term. First observe that

[TABLE]

where we used the Schwarz inequality in Schatten spaces in the first and the Kato-Seiler-Simon bound [27, Thm. 4.1] in the second estimate; a more precise bound is stated in [6, Theorem 4.3]. We also used that the norm of the kernel of $|R_{0}({\rm i}\eta)|^{1/2}$ is dominated by twice the absolute value of the kernel of $|\sqrt{-\Delta}-{\rm i}\eta|^{-\frac{1}{2}}$ and that the $L^{4}$ norm of the function $||\cdot|-{\rm i}\eta|^{-\frac{1}{2}}$ is bounded from above by $\pi^{\frac{1}{4}}\eta^{-\frac{1}{4}}$. We set

[TABLE]

This guarantees that ${\rm i}\eta\in\rho(H+V)$ since $\gamma^{-1}<1$. By [27, Thm. 9.2] we have the bound

[TABLE]

and so, since $\|V\|_{1}\geq 1$,

[TABLE]

Since $g$ is holomorphic on the unit disk $\mathbb{D}$ and continuous up to the boundary, the following Blaschke condition holds,

[TABLE]

Here and in the following, every zero is counted according to its multiplicity. We have also used the normalization condition $g(0)=1$ in (2.56), which follows from $\psi({\rm i}\eta)=0$. To translate the result back to the $z$-plane we use the distortion bound

[TABLE]

The first inequailty follows from the Koebe distortion theorem, the second from an explicit computation. Combination of (2.55), (2.56) and (2.57) yields

[TABLE]

By the choice of $\eta$, this is equivalent to (2.40) for $z\in{\mathbb{C}}^{+}$.

Next, we consider the massive case. Without loss of generality we may restrict our attention to the case $m=1$; the general case follows by scaling. We use the same maps as in (2.48) except that we replace $\varphi$ by

[TABLE]

Instead of (2.50) we now have (again from (2.20)) the estimate

[TABLE]

Assume first that $\|V\|_{1}<\rho_{0}$. By (2.7) we may choose $\eta>0$ (independent of $V$) such that ${\rm i}\eta\in\rho(H+V)$, i.e. $h({\rm i}\eta)\neq 0$. Since $|\zeta({\rm i}\eta)|=1$, we have

[TABLE]

It follows that

[TABLE]

and thus

[TABLE]

In the first estimate we used the triangle inequality $|h({\rm i}\eta)|\geq 1-|h({\rm i}\eta)-1|$, the monotonicity of the logarithm and the elementary inequality $\log(1-x)\geq-x/(1-x)$ for $x=|h({\rm i}\eta)-1|\in(0,1)$. The second estimate follows from the definition of $\psi$, some elementary inequalities for positive numbers and the fact that $\|V\|_{1}^{2}<\|V\|_{1}$. Since the right hand side is no longer bounded, we cannot use (2.56). Instead, we have to use a more refined result due to [2]. In our case it implies that

[TABLE]

for any $\tau>0$. Using the estimates

[TABLE]

which follow from straightforward albeit tedious computations, directly from the definition of $\psi$ and by Koebe’s distortion theorem we infer from (2.64) that

[TABLE]

This is half of the desired bound (2.41) for $m=1$ and $\|V\|_{1}<\rho_{0}$. Observe that if $\|V\|_{1}<\rho_{0}/2$, then the right hand side of (2.66) is bounded from above by a constant multiple of $\max\{\|V\|_{1},\|V\|_{1}^{2}\}$, while $A_{m}(V)$ is bounded from below by a positive constant. Hence, it remains to prove (2.41) for $\|V\|_{1}\geq\rho_{0}/2$ and with $A_{m}(V)=(1+\|V\|_{2}^{4})^{2+\tau}$. Since we have $|\sqrt{1-\Delta}-{\rm i}\eta|\geq|\sqrt{-\Delta}-{\rm i}\eta|$, (2.52) holds, and we make the same choice of $\eta$ and $\gamma$ as in the massless case. A repetition of the above arguments yields that

[TABLE]

and finally

[TABLE]

For $\|V\|_{1}\geq\rho_{0}/2$, we have that $\gamma^{-1}\leq 4(\gamma\rho_{0}^{2})^{-1}\|V\|_{1}^{2}$. Hence, by the choice of $\eta$ and $\gamma$, this is the desired bound. ∎

3. Applications

3.1. Damped wave equation in $L^{2}({\mathbb{R}})$

Our firs non-self-adjoint application motivated from physics is the damped wave equation

[TABLE]

where the damping $a$ and the potential $q_{0}\in{\mathbb{R}}_{+}$ satisfy

[TABLE]

The second order scalar equation (3.1) can be reformulated as a first order system, suitable for spectral analysis, in the form

[TABLE]

The equivalence to another (perhaps more intuitive) system with $v=u_{t}$ is extensively discussed e.g. in [16]. We view $G$ as an operator in $L^{2}({\mathbb{R}};{\mathbb{C}}^{2})$ and we employ a similarity transformation that brings $G$ to the form studied in Section 2. Let

[TABLE]

Then a straightforward calculation yields

[TABLE]

In the simplest case with $a_{1}=0$, (3.5) immediately gives that

[TABLE]

and thus illustrates the well-known spectral picture exhibiting the effect of damping (the shift of the spectrum to the left). As a corollary of the claims established or recalled in Section 2, we obtain new results in the situation when damping is no longer constant but satisfies (3.2).

Theorem 3.1**.**

Let $G$, $H$, $V$, $a$, $q_{0}$ and $\mu$ be as above. Then

[TABLE]

and the following holds.

1. \edefitn\selectfonti)

For any $\tau>0$, the Lieb-Thirring type inequality (2.41) and the bound on the number of eigenvalues (2.46) hold with

[TABLE]

and with the replacement $H\mapsto{\rm i}H$ and $m\mapsto{\rm i}\mu$. 2. \edefitn\selectfonti)

If $\|V\|_{1}<1$, then

[TABLE]

where $x_{0}$ and $r_{0}$ are as in (2.8) with $\|V\|_{1}$ as in (3.8). 3. \edefitn\selectfonti)

In the weak coupling regime, i.e. $a_{1}$ is replaced everywhere by $\varepsilon a_{1}$ with $\varepsilon>0$: if $\int_{{\mathbb{R}}}a_{1}(x)\,{{\rm d}}x>0$, then, for all sufficiently small $\varepsilon>0$, there are two eigenvalues $z_{\pm}(\varepsilon)$, $z_{-}=\overline{z_{+}}$, of $G$ satisfying

[TABLE]

3.2. Armchair graphene waveguides

As a second application motivated from physics, we consider the two-dimensional Dirac operator of an infinite straight graphene waveguide $\Omega=(-a,a)\times{\mathbb{R}}$,

[TABLE]

where $\tau:=-{\rm i}\partial_{1}+\partial_{2}$ and $\tau^{*}:=-{\rm i}\partial_{1}-\partial_{2}$ is the formal adjoint. The domain of $D$ consists of spinors $\psi\in H^{1}(\Omega;{\mathbb{C}}^{4})$ satisfying so-called armchair boundary conditions

[TABLE]

where $0\leq\Theta<2\pi$ depends on the geometry of the waveguide. It was proved in [15, Prop. 1, 19] that $D$ is self-adjoint and that the spectrum is given by

[TABLE]

where

[TABLE]

To simplify the presentation in the sequel, we restrict ourselves to the case when $\Theta\in(0,\pi)$ and thus $E_{0}=-\xi_{0}>0$; the results can be extended in a straightforward way to the other cases.

Although the algebraic structure of Dirac waveguide operators is more complicated than in the Laplacian (or Schrödinger) case, it might be helpful to remark that the numbers $\{\xi_{n}\}_{n\in{\mathbb{Z}}}$ play the role of spectral thresholds in the essential spectrum of $D$. The corresponding (normalized in $L^{2}((-a,a);{\mathbb{C}}^{4})$) transverse eigenfunctions read

[TABLE]

and the set $\{\Psi_{n}^{\sigma}\}_{n\in{\mathbb{Z}},\sigma\in\{+,-\}}$ forms an orthonormal basis in $L^{2}((-a,a);{\mathbb{C}}^{4})$. To proceed with spectral analysis of perturbations of $D$, we derive a convenient representation of the resolvent of $D$ based on its decomposition into transverse modes. Moreover, we employ a unitary transformation bringing $D$ and its resolvent close to the form of the Dirac operator $H$ investigated in Section 2. Notice that $-\xi_{n}$ plays the role of $m$ in previous formulas.

Lemma 3.2**.**

Let $D$ be as in (3.11). Then, for all $z\notin(-\infty,-E_{0}]\cup[E_{0},\infty)$,

[TABLE]

where

[TABLE]

and, for every $f\in L^{2}(\Omega;{\mathbb{C}}^{4})$,

[TABLE]

Proof.

Notice that, for any $g\in H^{1}({\mathbb{R}};{\mathbb{C}})$, we have

[TABLE]

Thus (3.16) follows by standard arguments. ∎

In the following, we investigate the spectrum of $D+V$ where

[TABLE]

To keep the connection to Section 2, in the proofs we always employ the unitary transformation $\Sigma$ and thus instead of $V$ we in fact work with

[TABLE]

Remark 3.3** (On the assumption $\|V\|\in L^{2}(\Omega)$).**

Similarly to the case of the one-dimensional Dirac operator, the condition $\|V\|\in L^{2}(\Omega)$ is imposed for convenience only. However, the situation is slightly different for the waveguide: We cannot drop the $L^{2}$ norm completely, but merely replace it (at the expense of using a more complicated definition of the sum $D+V$) by an $L^{1+\varepsilon}$ norm, where $\varepsilon>0$ is arbitrary. The $\varepsilon$ loss takes place in an orthogonality argument for estimating an infinite sum in the proof of Lemma 3.5. We do not know if this is just a technical issue.

Lemma 3.4**.**

Let $D$ be as in (3.11) and $V$ as in (3.20). Then

[TABLE]

Proof.

We show $V$ is relatively compact with respect to $D$, so the essential spectrum $\sigma_{\rm e3}$ is preserved, see [11, Thm. IX.2.1]. In the following, we employ the unitary transformation $\Sigma$.

First, using the explicit kernel (2.10) with $z=0$ and $m=-\xi_{n}$ and $\|W\|\in L^{2}(\Omega)$, we verify that $WH_{n}^{-1}P_{n}$ is a Hilbert-Schmidt operator with

[TABLE]

We will show that the series

[TABLE]

having Hilbert-Schmidt operators as summands, is convergent in the operator norm; this will imply that $VD^{-1}$ is compact.

To show the convergence of (3.24), we approximate $W$ by bounded potentials $W_{n}$ defined by

[TABLE]

where $C_{n}$ is chosen such that $\|W\mathbf{1}_{\{x:\|W(x)\|>C_{n}\}}\|_{2}\leq 2^{-n}$; this is possible e.g. by Chebychev’s inequality. In summary, we have chosen $W_{n}$ such that

[TABLE]

By the mutual orthogonality of ${\operatorname{Ran}}P_{n}$, we have for any $N\in\mathbb{N}$ with $N\geq 1$ that

[TABLE]

On the other hand, we have

[TABLE]

The claim is proved. ∎

3.2.1. Weakly coupled eigenvalues in armchair waveguides

We analyze the eigenvalues of $D+\varepsilon V$ emerging from the edges of the essential spectrum and their asymptotics as $\varepsilon\to 0+$. Here $V$ is assumed to satisfy (3.20) and thus $\sigma_{\rm e3}(D+\varepsilon V)=\sigma_{\rm e3}(D)$ for every $\varepsilon\geq 0$ by Lemma 3.2. As in the one-dimensional case, the main tool is the Birman-Schwinger principle. To be able to use the formulas from the one-dimensional case, we employ the unitary transformation $\Sigma$, see (3.17), and representation of the resolvent of $D$ in (3.16).

By standard arguments, it can be verified that $z\in\sigma_{\rm p}(D+\varepsilon V)$ if and only if $-1\in\sigma_{\rm p}(\varepsilon Q)$, where the Birman-Schwinger operator $Q$ has the form

[TABLE]

(the bar denotes the closure) and where $A$, $B$ are multiplication operators by the matrices $\mathcal{A}$ and $\mathcal{B}$ stemming from the polar decomposition of $W$, i.e.

[TABLE]

We further decompose $Q$ into a singular and regular part, namely $Q=L+M$. All these integral operators have explicit kernels; nonetheless, we display here in detail only the formula for $L$. After straightforward manipulations employing the formulas (2.10) and (3.15), we get

[TABLE]

where (with $\Upsilon$ as in (2.18))

[TABLE]

and

[TABLE]

For the regular part $M=M_{1}+M_{2}$, we have

[TABLE]

Lemma 3.5**.**

Let $V$ be as in (3.20) and $M_{1}$, $M_{2}$ be as in (3.34). Then

[TABLE]

Proof.

The estimate of $\|M_{1}\|$ is almost the same as in the proof of Lemma 2.1, we omit the details.

To prove the estimate for $\|M_{2}\|$ let $f,g\in C_{c}^{\infty}(\Omega;{\mathbb{C}}^{4})$, normalized in $L^{2}(\Omega;{\mathbb{C}}^{4})$ and assume that $|z-|\xi_{0}||<1/2$. From formula (3.15) it is straightforward to obtain the estimate

[TABLE]

for all $n\in{\mathbb{Z}}$, all $\sigma\in\{+,-\}$ and for almost all $x_{2}\in{\mathbb{R}}$, where

[TABLE]

An analogue of this inequality holds for $A,f$ replaced by $B,g$. By Schwarz’s and Bessel’s inequality,

[TABLE]

where we used (3.36) in the next-to-last inequality and the Kato-Seiler-Simon inequality, see [27, Thm. 4.1], in the last inequality. The factor $4$ comes from estimating the matrix operator $H_{n}$ in terms of the scalar operator $(-\partial_{2}^{2}+\xi_{n}^{2})^{\frac{1}{2}}$. The supremum over $n$ and $z$ in the final expression is finite. ∎

Having established suitable estimates of the regular part $M$, we prove the main result of this section.

Theorem 3.6**.**

Let $V$ be as in (3.20) and let $W=\Sigma V\Sigma^{-1}$ with $\Sigma$ as in (3.17). Define the matrices

[TABLE]

and

[TABLE]

If $u^{+}\in{\mathbb{C}}$ is a solution of $\det(u^{+}U^{+}+I_{2})=0$ that satisfies $\operatorname{Re}u^{+}>0$, then for any sufficiently small $\varepsilon>0$, there exists an eigenvalue $z_{+}(\varepsilon)$ of $D+\varepsilon V$ satisfying

[TABLE]

Similarly, if $u^{-}\in{\mathbb{C}}$ is a solution of $\det(u^{-}U^{-}+I_{2})=0$ that satisfies $\operatorname{Re}u^{-}<0$, then for any sufficiently small $\varepsilon>0$, there exists an eigenvalue $z_{-}(\varepsilon)$ of $D+\varepsilon V$ satisfying

[TABLE]

Proof.

The strategy and individual steps are the same as in the one-dimensional case (Theorem 2.2) thus we indicate only the differences. Employing a decomposition as in (2.27) and the fact that the kernel of $L$ is separated, we convert the problem $-1\in\sigma_{\rm p}(Q(z))$ into the algebraic equation

[TABLE]

with $U$ as in (3.43) and $U_{1}=\mathcal{O}\left(\frac{\varepsilon\|M\|}{1-\varepsilon\|M\|}\right)$. As an initial guess we consider

[TABLE]

for which we obtain $\det\left(I_{4}+\varepsilon\Upsilon_{2}U\right)=\mathcal{O}(\varepsilon^{2})$ as $\varepsilon\to 0+$; the latter can be verified by the Laplace expansion of the determinant. The rest of the proof follows the lines of the one of Theorem 2.2, employing the estimates on $M$ from Lemma 3.5 and formulas (3.32). ∎

Example 3.7**.**

In particularly simple case where $V=\operatorname{diag}(v_{1},v_{2},v_{3},v_{4})$ with $v_{j}=v_{j}(x_{2})$, $j=1,\dots,4$, and with $a=\Theta^{2}/(8\sin^{2}(\Theta/4))$, straightforward calculations reveal that

[TABLE]

Thus, depending on the sign of $\operatorname{Re}u^{\pm}$, we obtain eigenvalues $z_{\pm}$ obeying (3.45) or (3.46).