An overview on the standing waves of nonlinear Schr\"odinger and Dirac equations on metric graphs with localized nonlinearity

TL;DR

This paper reviews the existence and nonexistence of standing waves in nonlinear Schrödinger and Dirac equations on metric graphs with localized nonlinearity, highlighting their similarities, differences, and the nonrelativistic limit convergence.

Contribution

It provides a comparative overview of standing wave solutions for NLSE and NLDE on metric graphs, including the nonrelativistic limit analysis.

Findings

Existence and nonexistence results for standing waves in NLSE and NLDE.

Analysis of subcritical and critical cases for NLSE.

Convergence of bound states from NLDE to NLSE in the nonrelativistic limit.

Abstract

We present a brief overview on the existence/nonexistence of standing waves for the NonLinear Schr\"odinger and the NonLinear Dirac Equations (NLSE/NLDE) on metric graphs with localized nonlinearity. We first focus on the NLSE, both in the subcritical and the critical case, and then on the NLDE, highlighting similarities and differences with the NLSE. Finally, we show how the two equations are related in the nonrelativistic limit, proving the convergence of bound states.

Figures

Equations

$$L^p(\mathcal{G}):=\underset{e\in \mathrm{E}}{⨁}{L}^p(I_e),p\in [1,\infty ],\text{and}{H}^m(\mathcal{G}):=\underset{e\in \mathrm{E}}{⨁}{H}^m(I_e),m\in \mathbb{N},$$

$$ı{\partial }_{t}w=-{\Delta }_{\mathcal{G}}w-|w{|}^{p-2}w,p\ge 2,$$

$$-{\Delta }_{|{\oplus }_{e\in \mathrm{E}}{C}_0^{\infty }(\mathring{I_e})},$$

$$w(t,x)=e^{-i\lambda t}u(x),$$

$$-{\Delta }_{\mathcal{G}}u-|u{|}^{p-2}u=\lambda u.$$

$$(-{\Delta }_{\mathcal{G}}^{\delta ,\alpha }u{)}_{|I_e}:=-u_e^{\prime \prime }\forall e\in \mathrm{E},\forall u\in \mathrm{dom}(-{\Delta }_{\mathcal{G}}^{\delta ,\alpha }),$$

$$\mathrm{dom}(-{\Delta }_{\mathcal{G}}^{\delta ,\alpha }):=\left\{u\in H^2(\mathcal{G}):u\text{ satisfies eqrefeq-cont and eqrefeq-delta}\right\},$$

$$u_{e_1}(\mathrm{v})=u_{e_2}(\mathrm{v}),\forall e_1,e_2\succ \mathrm{v},\forall \mathrm{v}\in \mathcal{K},$$

$$\sum _{e\succ v}\frac{d{u}_e}{d{x}_e}(\mathrm{v})=\alpha u(\mathrm{v}),\forall \mathrm{v}\in \mathcal{K}$$

$$-{\Delta }_{\mathcal{G}}u-{\chi }_{\mathcal{K}}|u{|}^{p-2}u=\lambda u.$$

$$ı{\partial }_t\Psi ={\mathcal{D}}_{\mathcal{G}}\Psi -|\Psi {|}^{p-2}\Psi ,p\ge 2,$$

$$\mathcal{D}:=-ıc\frac{d}{dx}\otimes {\sigma }_1+m{c}^2\otimes {\sigma }_3,$$

$${\sigma }_1:=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right)\text{and}{\sigma }_3:=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right).$$

$${\mathcal{D}}_{\mathcal{G}}\psi -|\psi {|}^{p-2}\psi =\omega \psi .$$

$${\mathcal{D}}_{\mathcal{G}}\psi -{\chi }_{\mathcal{K}}|\psi {|}^{p-2}\psi =\omega \psi .$$

$${\mathcal{D}}_{\mathcal{G}}\psi -|⟨\psi ,{\sigma }_3\psi ⟩{|}^{(p-2)/2}\psi =\omega \psi$$

$$E(u,\mathcal{K},p):=\frac{1}{2}{\int }_{\mathcal{G}}|u^{\prime }{|}^{2}dx-\frac{1}{p}{\int }_{\mathcal{K}}|u{|}^{p}dx,p\ge 2,$$

$${\mathcal{H}}_{\mu }(\mathcal{G}):=\left\{u\in \mathcal{H}(\mathcal{G}):\Vert u{\Vert }_{2,\mathcal{G}}^2=\mu \right\}$$

$$\mathcal{H}(\mathcal{G}):=\left\{u\in H^1(\mathcal{G}):u\text{ saisfies eqrefeq-cont}\right\}.$$

$$\lambda=\lambda(u):=\frac{1}{\mu}\bigg{(}\int_{\mathcal{G}}|u^{\prime}|^{2}\,dx-\int_{\mathcal{K}}|u|^{p}\,dx\bigg{)}.$$

$$u(x)=e^{ı\varphi }v(x)$$

$$\Vert u{\Vert }_{p,\mathcal{G}}^p\le \mathcal{C}(\mathcal{G},p)\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^{\frac{p}{2}-1}\Vert u{\Vert }_{2,\mathcal{G}}^{\frac{p}{2}+1},\forall u\in \mathcal{H}(\mathcal{G}),$$

$$\Vert u{\Vert }_{\infty ,\mathcal{G}}\le \mathcal{C}(\mathcal{G},\infty )\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^{\frac{1}{2}}\Vert u{\Vert }_{2,\mathcal{G}}^{\frac{1}{2}},\forall u\in \mathcal{H}(\mathcal{G})$$

$$E(u,\mathcal{K},p)\ge \frac{1}{2}\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^2-\frac{\mathcal{C}(\mathcal{G},p)}{p}\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^{\frac{p}{2}-1}\Vert u{\Vert }_{2,\mathcal{G}}^{\frac{p}{2}+1},\forall u\in \mathcal{H}(\mathcal{G})$$

$$E(u,\mathcal{K},p)\ge \frac{1}{2}\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^{\frac{p}{2}-1}\left(\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^{\frac{6-p}{2}}-\frac{2{\mu }^{\frac{p+2}{4}}\mathcal{C}(\mathcal{G},p)}{p}\right),\forall u\in {\mathcal{H}}_{\mu }(\mathcal{G})$$

$$\underset{u\in {\mathcal{H}}_{\mu }(\mathcal{G})}{\inf }E(u,\mathcal{K},p)\le 0$$

$$\underset{u\in {\mathcal{H}}_{\mu }(\mathcal{G})}{\inf }E(u,\mathcal{K},p)<0⟹\text{a ground state does exist}.$$

$${\mu }^{\frac{p-2}{6-p}}|\mathcal{K}|>N^{\frac{4}{6-p}}{c}_p,$$

$$c_{p}=\bigg{[}\bigg{(}\frac{p(p-4)}{16}\bigg{)}^{\frac{2}{p-2}}+\frac{p}{8}\bigg{(}\frac{p(p-4)}{16}\bigg{)}^{\frac{4-p}{p-2}}\bigg{]}^{\frac{p-2}{6-p}},$$

$$\mu^{\frac{p-2}{6-p}}|\mathcal{K}|<\bigg{(}\frac{p}{2}\bigg{)}^{\frac{2}{6-p}}\frac{\mathcal{C}(\mathcal{G},p)^{\frac{4-p}{6-p}}}{\mathcal{C}(\mathcal{G},\infty)^{p}},$$

$$\mu \mapsto \sigma \mu \mathcal{G}\mapsto {\sigma }^{\frac{2-p}{6-p}}\mathcal{G}.$$

$$E(u,\mathcal{K},6)\ge \frac{1}{2}\Vert u^{\prime }{\Vert }_{2,\mathcal{G}}^2\left(1-\frac{{\mu }^2\mathcal{C}(\mathcal{G},6)}{3}\right),\forall u\in {\mathcal{H}}_{\mu }(\mathcal{G})$$

Full text

AMS Subject Classification: 35R02, 35Q55, 81Q35, 35Q40, 49J40, 49J35, 58E05, 46T05.

Keywords: metric graphs, NLS, NLD, ground states, bound states, localized nonlinearity, nonrelativistic limit.

1. Introduction

The aim of this paper is to present the state of the art on the study of the standing waves of the NonLinear Schrödinger and the NonLinear Dirac Equations (NLSE/NLDE) on metric graphs with localized nonlinearities.

In the following, by metric graph we mean the locally compact metric space which one obtains endowing a multigraph $\mathcal{G}=(\mathrm{V},\mathrm{E})$ with a parametrization that associates each bounded edge $e\in\mathrm{E}$ with a closed and bounded interval $I_{e}=[0,\ell_{e}]$ of the real line, and each unbounded edge $e\in\mathrm{E}$ with a (copy of the) half-line $I_{e}=\mathbb{R}^{+}$ (an extensive description can be found in [7, 14] and references therein). As a consequence, functions on metric graphs $u=(u_{e})_{e\in\mathrm{E}}:\mathcal{G}\to\mathbb{R},\,\mathbb{C}$ have to be seen as bunches of functions $u_{e}:I_{e}\to\mathbb{R},\,\mathbb{C}$ such that $u_{|_{e}}=u_{e}$. Consistently, Lebesgue and Sobolev spaces are defined as

[TABLE]

and are equipped with the natural norms. Moreover, throughout $L^{p}$-norms are denoted by $\|u\|_{p,\mathcal{G}}$ and the $H^{1}$-norm is denoted by $\|u\|$, for the sake of simplicity.

It is also worth recalling that, in the present paper we limit ourselves to focus on the case of metric graphs $\mathcal{G}$ satisfying the following hypothesis:

• (H1)

$\mathcal{G}$ is connected and noncompact;

• (H2)

$\mathcal{G}$ has a finite number of edges;

• (H3)

$\mathcal{G}$ has a non-empty compact core $\mathcal{K}$ (which is the subgraph of $\mathcal{G}$ consisting of all its bounded edges)

(see, e.g., Figure 1).

The growing interest in the study of evolution equations on metric graphs and networks is due to the fact that they are regarded as effective models for the dynamics of physical systems constrained in branched spatial structures (see [10] and references therein).

In particular, motivated by the study of qualitative properties of Bose-Einstein Condensates in ramified traps (see, e.g., [34]), in the last years a considerable attention has been devoted to the focusing NLSE, i.e.

[TABLE]

where $-\Delta_{\mathcal{G}}$ is a suitable self-adjoint realization of the operator

[TABLE]

and, precisely, on the existence of standing waves of (1). Those are solutions of the form

[TABLE]

with $\lambda\in\mathbb{R}$ and $u\in L^{2}(\mathcal{G})$ solving the stationary version of (1), namely,

[TABLE]

The first results in this direction (e.g., [1, 2, 3, 4]) considered the so-called infinite $N$-star graph (see, e.g., Figure 2), in the case where $-\Delta_{\mathcal{G}}$ is the Laplacian with $\delta$-type vertex conditions, that is

[TABLE]

where

[TABLE]

for some $\alpha\in\mathbb{R}$ ($e\succ\mathrm{v}$ meaning that the edge $e$ is incident at the vertex $\mathrm{v}$, while $\tfrac{du_{e}}{dx_{e}}(\mathrm{v})$ stands for $u_{e}^{\prime}(0)$ or $-u_{e}^{\prime}(\ell_{e})$ according to whether $x_{e}$ is equal to [math] or $\ell_{e}$ at $\mathrm{v}$).

On the other hand, in the case $\alpha=0$ which is usually called Kirchhoff Laplacian (and which will be denoted simply by $-\Delta_{\mathcal{G}}$ in place of $-\Delta_{\mathcal{G}}^{\delta,0}$ in the sequel, for the sake of simplicity) more general graphs have been studied (precisely, any graph satisfying (H1)-(H2)). We mention, in this regard, [7, 8, 9] for a discussion of the existence of ground states (i.e., those standing waves that minimize the energy functional associated with (2)) and [11, 21, 36, 41, 42] concerning, more generally, excited states. We also mention [22] where the same problems are studied in the presence of an external potential.

A modification of this model, introduced in [33, 40], consists of assuming that the nonlinearity affects only the compact core of the graph (which then must be supposed non empty as in (H3)), so that (2) reads

[TABLE]

where $\chi_{\mathcal{K}}$ is the characteristic function of $\mathcal{K}$. The existence of stationary solutions to (5) has been discussed in [47, 48, 50] in the $L^{2}$-subcritical case $p\in(2,6)$ and, more recently, in [26, 27] in the critical case $p=6$. A detailed description of these results will be presented in Section 2.

For the sake of completeness, we also remark that the NLSE on compact graphs (which, in particular, do not fulfill (H1)) has been studied, e.g., in [20, 24, 28, 38]; while the case of one or higher-dimensional periodic graphs (which, in particular, do not fulfill (H2) as, for instance, in Figure 3) has been addressed, e.g., by [6, 25, 32, 43].

Besides the NLS, other dispersive PDEs on metric graphs have been explored in the last years. We mention, for instance, the case of [39] which deals with the dynamics for the Airy equation (motivated by the study of the KDV equation) on star graphs.

Recently, we started a new research project concerning the NonLinear Dirac Equation (NLDE) on metric graphs in [18]. The physical motivations for such a model mainly come from solid state physics and nonlinear optics. The existence of Dirac solitons in Bose-Einstein condensates and optical lattices and their concrete realization in discrete waveguide arrays have been recently investigated in [35, 51]. In that case one may expect to recover the metric graph model in an appropriate scaling regime. We also mention that a rigorous mathematical study of the dynamics and the existence of Dirac solitons on lattices has been recently treated in [12, 15, 17, 16, 31].

The investigation of stationary solution of the NLDE in this context has been first proposed by [46] for the case of the infinite 3-star graph (see Figure 4). However, in [46] the authors considered the case of an extended nonlinearity, i.e.

[TABLE]

where $\mathcal{D}_{\mathcal{G}}$ is a suitable self-adjoint realization of the one-dimensional Dirac operator $\mathcal{D}$, i.e.

[TABLE]

$m>0$ and $c>0$ representing the mass of the generic particle of the system and a relativistic parameter (respectively), and $\sigma_{1}$ and $\sigma_{3}$ representin the Pauli matrices, i.e.

[TABLE]

Here the equation of the standing waves $\Psi(t,x)=e^{-\imath\omega t}\psi(x)$ reads

[TABLE]

In fact, in order to deal with more complex graph topologies, in [18] we restricted ourselves to the case of a Kirchhoff-type extension of the Dirac operator (for details see Section 3.1) and, most importantly, we considered the case of a localized nonlinearity, that is

[TABLE]

In addition, in order to test the consistency of this model we also proved (see Section 4) the convergence (of the bound states) of (7) to (the bound states of) (5) in the nonrelativistic limit, namely as “$c\to\infty$”.

Note that in (7) (as in the equation studied in [46]) the nonlinearity is a pure power, thus being manifestly non covariant. This kind of nonlinearities typically arises in nonlinear optics. We also stress that from a theoretical point of view there is no conceptual contradiction as the NLDE should be intended only as an effective model.

Nevertheless, the equation obtained replacing the nonlinearity $|\psi|^{p-2}\,\psi$ with a covariant one, i.e.

[TABLE]

(the nonlinearity now being defined in terms of a Lorentz-scalar) is great interest. This equation might, indeed, be interpreted as the model for a truly relativistic particle (a Dirac fermion), whose dynamics is constrained to a one-dimensional structure. In this case the model should be Lorentz-covariant (strictly speaking, in this case the equation should be rewritten in covariant form). We expect that this will require some technical adjustments in order to deal with the indefiniteness of the nonlinearity. This might be quite relevant in the study of the nonrelativistic limit, since in [18] we used the fact that a positive nonlinearity allows to get a priori estimates uniform in $c$.

Finally, we mention that the discussion of extended nonlinearities presents some technical problems that we were not able to overcome thus far. However, we plan to investigate this question in a forthcoming paper.

2. Nonlinear Schrödinger equation

The study of the existence/nonexistence and multiplicity of bound states of (5) carried on in [26, 27, 47, 48, 50] is completely based on variational methods.

This is due to the simple observation that $L^{2}$-solutions of (5), with $-\Delta_{\mathcal{G}}$ denoting the Kirchhoff Laplacian, arise also as constrained critical points of the energy functional

[TABLE]

on the manifold

[TABLE]

for some fixed (but arbitrary) $\mu>0$, where

[TABLE]

Moreover, if $u$ is a constrained critical point, also known as bound state, then an easy computation shows that the Lagrange multiplier $\lambda$ can be found as

[TABLE]

Remark 2.1*.*

Note that $\mathcal{H}(\mathcal{G})$ is the form domain of $-\Delta_{\mathcal{G}}$. In addition, it is worth mentioning that the choice of the manifold (8) is also related to the fact that the time dependent counterpart of (5) is mass-preserving.

2.1. Ground states

The first step in the study of the bound states consists in looking for ground states, namely the constrained minimizers of $E(\cdot,\mathcal{K},p)$ on the manifold $\mathcal{H}_{\mu}(\mathcal{G})$.

Ground states can be shown (see [7]) to be of the form

[TABLE]

where $\phi$ is a fixed phase factor and $v$ is real valued. Hence, in the following, in minimization problems we will always limit ourselves to consider real-valued functions.

For more general bound states we cannot prove at the moment that such a property holds. Therefore, while in existence and multiplicity results we will limit to real-valued functions, nonexistence results concern also complex-valued functions.

2.1.1. The subcritical case: $p\in(2,6)$

We start reporting on the existence of constrained minimizers in the so-called $L^{2}$-subcritical case, namely, when $p\in(2,6)$.

For $p>2$ and for every graph satisfying (H1)-(H2), it is possible to prove (see, e.g., [8, Proposition 2.1] and [50, Proposition 4.1] for real-valued functions and [48, Proposition 2.6] for complex-valued functions) that the following Gagliardo-Nirenberg inequalities hold:

[TABLE]

with $\mathcal{C}(\mathcal{G},p)$ and $\mathcal{C}(\mathcal{G},\infty)$ denoting the optimal constants of the respective inequalities. Then

[TABLE]

so that

[TABLE]

and, if $p\in(2,6)$, this immediately entails that the functional $E(\cdot,\mathcal{K},p)$ is bounded from below on the manifold $\mathcal{H}_{\mu}(\mathcal{G})$.

Remark 2.2*.*

Note that, in fact, (10) and (11) hold for every $u\in H^{1}(\mathcal{G})$ under the sole assumption (H2), up to a redefinition of the optimal constants. However, we chose to mention the $\mathcal{H}(\mathcal{G})$-version, which holds even if (H2) is not fullfilled. We remark that the best constants of this version play a crucial role in the analysis of the ground states.

In addition to (12), as a general fact, one can prove that

[TABLE]

and that

[TABLE]

As a consequence, the following results can be established

Theorem 2.1** ([50, Theorems 3.3$\&$3.4],[48, Corollary 3.4]).**

Let $\mathcal{G}$ satisfy (H1)-(H3) and let $p\in(2,6)$. Therefore:

if $p<4$, then there exists a ground state for every $\mu>0$;

**

• 2)

if $p\geq 4$, then:

**

• (i)

whenever

[TABLE]

where $N$ is the number of half-lines of $\mathcal{G}$ and

[TABLE]

there exists a ground state of mass $\mu$;

• (ii)

whenever

[TABLE]

there does not exist any ground state of mass $\mu$.

The proof of the previous theorem is clearly based on (13) and (14). Precisely, from (14) there results that the existence part can be obtained exhibiting a function with negative energy, usually called competitor. This can be achieved for every mass $\mu$ when $p\in(2,4)$, while when $p\in[4,6)$ this is possible only under (15). Such competitors are constant on the compact core and exponentially decreasing on half-lines. In this way they minimize the kinetic energy on the compact core and have the exact qualitative behavior of the bound states on the half-lines (where the problem is linear). Clearly, such competitors cannot be minimizers since they do not fulfill (4) with $\alpha=0$.

On the other hand, (13) shows that, in order to prove nonexistence, it is sufficient to prove that every function possesses a strictly positive energy level. This is actually, the idea behind the first proof of the nonexistence result for ground states. However, there is a different strategy, which is similar to that of the proof of Theorem 2.5, which allows to prove item 2)(ii) in a more straight way and with the sharper threshold given by (16).

Finally, it is worth mentioning that the negativity of the energy levels of the minimizers entails that the associated Lagrange multiplier, given by (9), is negative and, especially, that the results of Theorem 2.1 are invariant under the homotetic transformation

[TABLE]

In other words, if one sets for instance $p=4$, a problem on a graph $\mathcal{G}$ with a given mass constraint $\mu$ is completely equivalent to a problem with half the mass and a “doubled” graph, and viceversa.

2.1.2. The critical case: $p=6$

The study of the so-called $L^{2}$-critical case $p=6$ has been successfully managed only very recently. Here the main difficulty comes from the fact that the boundedness from below of the energy functional strongly depends on the the mass constraint.

Indeed, in this case (12) reads

[TABLE]

and hence the boundedness of functional clearly depends on the value of $\mu$. In particular, it turns out that the discriminating value is related to the best constant of (10), that is, the critical mass $\mu_{\mathcal{G}}$ is defined by

[TABLE]

and is well known when $\mathcal{G}=\mathbb{R},\mathbb{R}^{+}$, i.e.

[TABLE]

In addition, there holds

[TABLE]

and $\mu_{\mathbb{R}},\mu_{\mathbb{R}}^{+}$ are the sole values of the mass constraint at which the problem with the extended nonlinearity on $\mathbb{R},\mathbb{R}^{+}$ (respectively) admits a ground state (see, e.g., [23]).

However, one can easily see that, while for the problem with extended nonlinearity (treated in [9]) $\mu_{\mathcal{G}}$ is the proper parameter to be investigated, in the localized case one should better consider a reduced critical mass

[TABLE]

where $\mathcal{C}(\mathcal{K})$ denotes the optimal constant of the following modified Gagliardo-Nirenberg inequality

[TABLE]

namely

[TABLE]

In fact, using this new parameter one obtains

[TABLE]

which is clearly sharper than (17), as one can easily see that $\mathcal{C}(\mathcal{K})\leq\mathcal{C}(\mathcal{G},6)$.

Existence of ground states has been first investigated in [26] in the simple case of the tadpole graph (see Figure 5). Very recently, more refined results have been obtained without any assumption on the topology of the graph, thus providing a (almost) complete classification of the phenomenology.

In order to state such results, let

[TABLE]

and recall the following definitions for a metric graph:

a graph $\mathcal{G}$ is said to admit a cycle covering if and only if every edge of $\mathcal{G}$ belongs to a cycle, namely either a loop (i.e., a closed path of consecutive bounded edges) or an unbounded path joining the endpoints of two distinct half-lines (which are then identified as a single vertex at infinity);

• -

a graph $\mathcal{G}$ is said to possesses a terminal edge if and only if it contains an edge ending with a vertex of degree one.

Theorem 2.2** ([27, Theorem 1.1]).**

Let $\mathcal{G}$ satisfy (H1)-(H3). Then,

[TABLE]

In addition,

• (i)

if $\mathcal{G}$ has at least a terminal edge (as, for instance, in Figure 6(A)**), then

[TABLE]

and there is no ground state of mass $\mu$ for any $\mu>0$;

• (ii)

if $\mathcal{G}$ admits a cycle-covering (as, for instance, in Figure 6(B)**), then

[TABLE]

and there is no ground state of mass $\mu$ for any $\mu>0$;

• (iii)

if $\mathcal{G}$ has only one half-line and no terminal edges (as, for instance, in Figure 6(C)**), then

[TABLE]

and there is a ground states of mass $\mu$ if and only if $\mu\in[\mu_{\mathcal{K}},\mu_{\mathbb{R}}]$.

• (iv)

if $\mathcal{G}$ has no terminal edges, does not admit a cycle-covering, but presents at least two half-lines (as, for instance, in Figure 6(D)**), then

[TABLE]

and there is a ground states of mass $\mu$ if and only if $\mu\in[\mu_{\mathcal{K}},\mu_{\mathbb{R}}]$, provided that $\mu_{\mathcal{K}}\neq\mu_{\mathbb{R}}$.

Theorem 2.2 presents several differences with respect to its analogous for the extended case established in [9] (in addition to the fact that $\mu_{\mathcal{G}}$ is replaced with $\mu_{\mathcal{K}}$). For details on such differences we refer the reader to [27]. Here, we aim at highlighting a more crucial feature of the previous result. Indeed, while in the subcritical case the topology of the graph plays no role in establishing existence/nonexistence of ground states, here the topological classification of the graphs is the central point, as it occurs in the context of extended nonlinearities. Even though there is no rigorous explanation of such a phenomenon, this seems to be related to the fact that the critical power of the nonlinearity “breaks” the characteristic lenght of the graph, i.e. $\mu^{\frac{p-2}{6-p}}|\mathcal{K}|$, thus making metric properties less relevant in the discussion.

Nevertheless, they preserve a role also in the critical case, at least in cases (iii) and (iv), as shown by the following

Theorem 2.3** ([27, Theorems 1.2$\&$1.3]).**

Estimates (18) and (19) are sharp in general; i.e., for every $\varepsilon>0$ there exist two non-compact metric graphs $\mathcal{G}_{\varepsilon}^{1},\mathcal{G}_{\varepsilon}^{2}$ (with compact cores $\mathcal{K}^{1}_{\varepsilon},\mathcal{K}^{2}_{\varepsilon}$), with one half-line and without terminal edges, such that

[TABLE]

and two non-compact metric graphs $\mathcal{G}_{\varepsilon}^{3},\mathcal{G}_{\varepsilon}^{4}$ (with compact cores $\mathcal{K}^{3}_{\varepsilon},\mathcal{K}^{4}_{\varepsilon}$), without terminal edges and cycle coverings and with at least two half-lines, such that

[TABLE]

In other words Theorem 2.3 shows that estimates (18) and (19) are sharp by exhibiting four suitable sequences of graphs. They can be constructed as follows:

the sequence $\mathcal{G}_{\varepsilon}^{1}$ can be constructed by considering a graph whose compact core does not admit a cycle covering (see, e.g., Figure 7(A)) and letting the lenght of one of its cut-edges, the edges whose removal disconnects the graph (e.g., $\widehat{e}$ in Figure 7(A)), go to infinity;

• 2)

the sequence $\mathcal{G}_{\varepsilon}^{2}$ can be constructed by considering a graph as in Figure 7(B)) and letting the lenght of the compact core go to infinity, keeping the total diameter of the compact core bounded (namely, thickening the compact core);

• 3)

the sequences $\mathcal{G}_{\varepsilon}^{3},\mathcal{G}_{\varepsilon}^{4}$ can be constructed by considering a signpost graph (see, e.g., Figure 7(C)) and letting the lenght of its cut-edge $\widetilde{e}$ go to infinity and to zero, respectively.

The proofs of Theorems 2.2 and 2.3 are quite technical and we refer the reader to [27] for a complete presentation. They are based on a careful analysis of the behavior of the constant $\mathcal{C}(\mathcal{K})$ and on the “graph surgery” and rearrangement techniques developed by [7, 8, 9]. In particular, given a nonnegative function $u\in\mathcal{H}(\mathcal{G})$, one uses the decreasing rearrangement

[TABLE]

and the symmetric rearrangement

[TABLE]

where

[TABLE]

and precisely the facts that

[TABLE]

that

[TABLE]

and that, for every function $u$ which possesses at least two preimages at each level,

[TABLE]

in order to construct suitable minimizing sequences or suitable competitors to exploit level arguments analogous to (14) to get compactness in the limit.

2.2. Bound states

Besides the investigation of ground states, one can also study standing waves which are not necessarily minimizers of the constrained energy, namely bound/excited states. They may arise as constrained critical points of the energy functional on $\mathcal{H}_{\mu}(\mathcal{G})$.

To the best of our knowledge all the known results on existence (and multiplicity), as well as on nonexistence, of general bound states strongly exploit the assumption on the subcriticality of the power nonlinearity (the critical case is open thus far).

2.2.1. Existence results

Existence and multiplicity results in this context have been proved by extending some well known techniques from Critical Point Theory (see, e.g., [5, 13, 45]).

However, the context of metric graphs presents two additional technical issues. First, it is not possible, in general, to gain compactness restricting to symmetric (in some suitable sense) functions if one does not want to restrict the discussion to symmetric graphs. Furthermore, the fact that graphs with non empty compact core are not homothetically invariant entails that multiple solutions of prescribed mass cannot be found by scaling arguments, in general.

The strategy used in [47], in order to solve these problems is, then, the following:

detect the energy levels at which the Palais-Smale condition is satisfied, namely detect the values $c\in\mathbb{R}$ such that any sequence $(u_{n})\in\mathcal{H}_{\mu}(\mathcal{G})$ satisfying

• (i)

$E(u_{n},\mathcal{K},p)\to c$,

• (ii)

$\|dE_{|\mathcal{H}_{\mu}(\mathcal{G})}(u_{n},\mathcal{K},p)\|_{\mathbb{T}^{\prime}_{u_{n}}\mathcal{H}_{\mu}(\mathcal{G})}\to 0$

(with $\mathbb{T}^{\prime}_{u_{n}}\mathcal{H}_{\mu}(\mathcal{G})$ denoting the topological dual of the tangent to the manifold $\mathcal{H}_{\mu}(\mathcal{G})$ at $u_{n}$), admits a subsequence converging in $\mathcal{H}_{\mu}(\mathcal{G})$.

• -

constructing suitable min-max levels.

In this case one can prove that the constrained functional possesses the Palais-Smale property only at negative levels and whence the min-max levels have to be negative. In addition, as the functional is even, it is possible to construct such min-max levels using the Krasnosel’skii genus (see, e.g., [37]); that is, for every $A\subset\mathcal{H}(\mathcal{G})\backslash\{0\}$ closed and symmetric, the natural number defined by

[TABLE]

In fact, one can prove that the suitable levels are

[TABLE]

where

[TABLE]

In this way it is possible to prove the following

Theorem 2.4** ([47, Theorem 1.2]).**

Let $\mathcal{G}$ satisfy (H1)-(H3) and let $p\in(2,6)$. For every $k\in\mathbb{N}$, there exists ${\mu_{k}}>0$ such that for every $\mu\geq{\mu_{k}}$ there exist at least $k$ distinct pairs $(\pm u_{j})$ of bound states of mass $\mu$. Moreover, for every $j=1,\dots,k$,

[TABLE]

where $\varphi_{\mu/j}$ denote the unique positive minimizer of $E(\cdot,\mathbb{R},p)$ constrained on $\mathcal{H}_{\mu/j}(\mathcal{G})$ and $\sigma_{k}(\mu)\to 0$ (exponentially fast) as $\mu\to\infty$. Finally, for each $j$, the Lagrange multiplier $\lambda_{j}$ relative to $u_{j}$ is negative.

The strategy hinted before cannot be easily adapted to the extended problem since in this case the Palais-Smale condition fails also at infinitely many negative energy levels and then it is open how to fix a min-max level in the gap between two of these levels.

On the other hand, it is also unclear how to extend this strategy to the critical case. Indeed the construction of negative energy min-max levels is based on the possibility of putting (on increasing the mass) several suitable truncations of scaled copies of $\varphi_{\mu/j}$ on the compact core of the graph keeping the total energy negative. However, no direct analogous of such a techniques is available for the critical case as the minimizers $\varphi_{\mu}$ do exist only at the critical mass and present a zero energy level.

2.2.2. Nonexistence results

Nonexistence results concern clearly only the regime $p\in[4,6)$ since for $p\in(2,4)$ the existence of bound states is guaranteed by Theorem 2.1.

The interesting part of nonexistence results for bound states in the subcritical regime is that they strongly relies both on metric and on topological features of the graph as one can see by the following

Theorem 2.5** ([48, Theorems 3.2$\&$3.5]).**

Let $\mathcal{G}$ satisfy (H1)-(H3) and let $p\in[4,6)$. Therefore,

• (i)

if the graph $\mathcal{G}$ satisfies

[TABLE]

then, there are no bound states of mass $\mu$ with $\lambda\leq 0$;

**

• (ii)

if $\mathcal{G}$ is a tree (i.e., no loops) with at most one pendant (see, e.g., Figure 8), then there is no bound state of mass $\mu$ with $\lambda\geq 0$, every $\mu>0$.

The most remarkable fact is that the dependence on metric or topological features in the nonexistence result is connected to the sign of the Lagrange multiplier. In particular, the first condition prevents the existence of bound states supported on the whole $\mathcal{G}$, as such functions cannot possess a negative Lagrange multiplier since they are in $L^{2}(\mathcal{G})$. On the other hand, one can check that the second part of Theorem 2.5, holds as well for $p\geq 6$ and even in the extended nonlinearity case.

Finally, it is worth observing that the condition of the second part of Theorem 2.5 is sharp. Indeed, one can easily construct counterexamples whenever the assumption of being a tree with at most one pendant is dropped. Precisely, if $\mathcal{G}$ possesses a loop, then one can define a bound state with a nonnegative Lagrange multiplier supported only on the loop (see, e.g., Figure 9(A)); if $\mathcal{G}$ is a tree with two pendants, instead, then one can define a bound state with nonnegative Lagrange multiplier supported on the path which joins the two pendants (see, e.g., Figure 9(B)).

3. Nonlinear Dirac equation

As already mentioned in Section 1, the study of bound state-solutions of NLDE (7) in the context of metric graphs has been introduced in [46]. To our knowledge, the first rigorous work on this subject for general graphs has been done in [18], for the case of localized nonlinearities.

Such discussion requires a suitable definition of the Dirac operator on graphs and the adaptation of some technques from Critical Point Theory.

Before going more into details, we recall that Dirac-type equations are of spinorial nature and, hence, here we de deal with vector-valued function on graphs ($2$-spinors). More precisely, we consider wave functions $\psi:\mathcal{G}\to\mathbb{C}^{2}$, that can be seen as a family of $2$-spinors on intervals, i.e. $\psi=(\psi_{e})$ with

[TABLE]

or equivalently as a $2$-components vector of functions on graphs, i.e. $\psi=(\phi,\chi)^{T}$, with $\phi=(\phi_{e})$ and $\chi=(\chi_{e})$. Accordingly,

[TABLE]

and

[TABLE]

endowed with the natural norms, which we denote, with a little abuse of notation, by $\|\psi\|_{p,\mathcal{G}}$ and $\|\psi\|$ (in the case $m=1$), respectively.

3.1. Remarks on the Dirac operator on graphs

As we mentioned before, the first step in the study of the NLDE on graphs is to define a self-adjoint realization of the Dirac operator on graphs.

A complete discussion of such a topic can be found, e.g., in [19, 44]. Here, however, we limit ourselves to the extension that we call of Kirchhoff type, since it represents the analogous to the Schrödinger operator with Kirchhoff conditions and, hence, corresponds to the free case, that is, the case where there is no interaction at vertices (for detail on such extension, refer to [18]).

For every fixed $m,\,c>0$, we define the Dirac operator of Kirchhoff-type as the operator

[TABLE]

with action

[TABLE]

$\sigma_{1},\sigma_{3}$ being the Pauli matrices defined in (6), and whose domain is

[TABLE]

where

[TABLE]

$\chi_{e}(\mathrm{v})_{\pm}$ standing for $\chi_{e}(0)$ or $-\chi_{e}(\ell_{e})$ according to whether $x_{e}$ is equal to [math] or $\ell_{e}$ at $\mathrm{v}$.

Such operator can be proved to be self-adjoint and, in addition, possesses spectral properties analogous to the standard Dirac operator on $\mathbb{R}$, that is

[TABLE]

even though, according to the structure of the graph, there might be eigenvalues embedded at any point of the spectrum.

The actual reason for which we call such operator of Kirchhoff type is better explained by Section 4. However, an informal justification is provided by the following formal computation. First, one can see that $\mathcal{D}_{\mathcal{G}}^{2}$ acts as $(-\Delta_{\mathcal{G}})\otimes\mathbb{I}_{\mathbb{C}^{2}}$. In addition, if one considers spinors of the type $\psi=(\phi,0)^{T}$ and assumes that they belong to the domain of $\mathcal{D}_{\mathcal{G}}^{2}$, namely that $\psi\in\mathrm{dom}(\mathcal{D}_{\mathcal{G}})$ and that $\mathcal{D}_{\mathcal{G}}\psi\in\mathrm{dom}(\mathcal{D}_{\mathcal{G}})$, one clearly sees that $\phi\in\mathrm{dom}(-\Delta_{\mathcal{G}})$.

Finally, compared to the Schrödinger case, a big difference is given by the definition of the quadratic form associated with $\mathcal{D}_{\mathcal{G}}$. In particular, here it is not explicitly known since Fourier transform is not available (in a simple way), as we are not in an euclidean space, and since classical duality arguments fail as it is not true in general that $H^{-1/2}(\mathcal{G},\mathbb{C}^{2})$ is the topological dual of $H^{1/2}(\mathcal{G},\mathbb{C}^{2})$.

It can be defined, clearly, using the Spectral Theorem, so that

[TABLE]

where $\mu_{\mathcal{D}_{\mathcal{G}}}^{\psi}$ denotes the spectral measure associated with $\mathcal{D}_{\mathcal{G}}$ and $\psi$. However, such a definition is not useful for computations. A more precise description of the quadratic form and its domain can be obtained arguing as follows (for more details, see [18]).

First, using Real Interpolation Theory, one can prove that

[TABLE]

whence

[TABLE]

On the other hand, according to (22) one can decompose the form domain as the orthogonal sum of the positive and negative spectral subspaces for the operator $\mathcal{D}_{\mathcal{G}}$, i.e.

[TABLE]

As a consequence, if one denotes $\psi^{+}:=P^{+}\psi$ and $\psi^{-}:=P^{-}\psi$ and recalls that it is possible to define a norm for $\mathrm{dom}(\mathcal{Q}_{\mathcal{D}_{\mathcal{G}}})$ as

[TABLE]

(with $\sqrt{|\mathcal{D}_{\mathcal{G}}|}\psi$ given by Borel functional calculus for $\mathcal{D}_{\mathcal{G}}$), then there results that

[TABLE]

Clearly, (23) and (24) are not explicit forms for the domain and the action of the quadratic form associated with $\mathcal{D}_{\mathcal{G}}$. Nevertheless, they present the suitable detalis in order to manage the computations required in the proofs of Theorems 3.1 and 4.1 below.

3.2. Bound states

The study of bound states of (7) presents some relevant differences with respect to the NLS case (5) (a difference which also arises in the extended case).

The main point here is the unboundedness from below of the spectrum of $\mathcal{D}_{\mathcal{G}}$ which makes the associated quadratic form strongly indefinite, even fixing the $L^{2}$-norm. As a consequence, the natural energy functional associated with (7), i.e.

[TABLE]

is unbounded from below even under the mass constraint $\|\psi\|_{2,\mathcal{G}}^{2}=\mu$. Hence, no minimization can be performed and also the search for constrained critical points presents several technical difficulties.

Therefore, the most promising strategy seems to be that of considering the action functional associated with (7) without any constraint. Then for a fixed $\omega\in\mathbb{R}$ one looks for critical points of the functional

[TABLE]

on the form-domain $\mathrm{dom}(\mathcal{Q}_{\mathcal{D}_{\mathcal{G}}})$.

To this aim, one needs to adapt to the metric graphs setting some well known techniques of Critical Point Theory. In particular, a suitable version of the so-called linking technique for even functional is exploited. We refer the reader to [45, 49], for a general presentation of those methods in an abstract setting, together with several applications, and to [29] which deals with nonlinear Dirac equations.

Critical levels for $\mathcal{L}$ are defined as the following min-max levels

[TABLE]

with

[TABLE]

where $h_{t}$ is the usual pseudo-gradient flow associated with $\mathcal{L}(\cdot,\mathcal{K},p)$, $S^{+}_{r}$ is the sphere of radius $r$ of $\mathrm{dom}(\mathcal{Q}_{\mathcal{D}_{\mathcal{G}}})^{+}$ and $\gamma$ denotes again the Krasnosel’skii genus. Indeed, one can check that, if $\mathcal{F}_{N}\neq\emptyset$, then there exists a Palais-Smale sequence $\left(\psi_{n}\right)\subset\mathrm{dom}(\mathcal{Q}_{\mathcal{D}_{\mathcal{G}}})$ at level $\alpha_{N}$, i.e.

[TABLE]

and in addition, there results

[TABLE]

for fixed $\rho,\,\widetilde{\rho}>0$.

As a consequence, in order to prove that $\mathcal{F}_{N}\neq\emptyset$, it suffices to show that the functional possesses a so-called linking geometry; namely, that for every $N\in N$ there exist $R=R(N,p)>0$ and an $N$-dimensional space $Z_{N}\subset\mathrm{dom}(\mathcal{Q}_{\mathcal{D}_{\mathcal{G}}})^{+}$ such that

[TABLE]

where

[TABLE]

and

[TABLE]

and that there exist $r,\rho>0$ such that

[TABLE]

(for a graphic intuition see, e.g., Figure 10).

Finally, checking the validity of the Palais-Smale condition at positive levels for the action functional it is possible to prove the following

Theorem 3.1** ([18, Theorem 2.11]).**

Let $\mathcal{G}$ satisfy (H1)-(H3), $p>2$ and $m,c>0$. Then, for every $\omega\in\mathbb{R}\backslash\sigma(\mathcal{D}_{\mathcal{G}})=(-mc^{2},mc^{2})$ there exists infinitely many (distinct) pairs of bound states of frequency $\omega$ of the NLDE, at mass $m$ and relativistic parameter $c$.

4. Nonrelativistic limit

The NLDE and the NLSE equation are clearly closely related as, physically, the latter should correspond to the nonrelativistic limit of the former. Heuristically, one expect to recover the NLSE from NLDE as the relativistic parameter tends to infinity (namely, “$c\to\infty$”), that is, neglecting relativistic effects. It is, then, particularly interesting to rigorously establish such a connection.

This has been first done in [30] for the Dirac-Fock equations proving the convergence of bound states to those of the (nonrelativistic) Hartree-Fock model.

More in detail, given a sequence $c_{n}\to\infty$, one is interested in the limit behavior of a sequence of bound states with frequencies $0<\omega_{n}<mc^{2}_{n}$. Precisely, one has to choose frequencies such that

[TABLE]

This represents a kind of renormalization, which corresponds to the fact that one has to subtract the rest energy of the particles, as this property is a peculiarity of relativistic theories which is absent in the nonrelativistic setting.

The strategy adopted in [18] is an adaptation of the one developed in [30]. Namely, first one has to establish $H^{1}$ uniform bounds for the sequence of the bound states and then, suitably manipulating (7), one has to prove that the lower component goes to zero while the sequence of upper components is a Palais-Smale sequence for the action functional associated to (5), namely for

[TABLE]

at a fixed $\lambda<0$. Therefore, up to a compactness argument, the existence of a limit satisfying the NLSE is proved.

However, we would like to remark that some relevant differences are present here with respect to [30], that call for some relevant modifications of the proofs. In particular, in [30] bound states were obtained as constrained critical points of the energy. As a consequence uniform boundedness of the $L^{2}$-norms and the non-triviality limit are easily obtained. On the contrary, in the case considered in [18], those properties are not a priori guaranteed and represent two main points of the proofs.

However, for the sake of completeness, we mention that in [30] one of the most delicate parts the proof is to estimate the sequence of the Lagrange multipliers corresponding to the $L^{2}$-constraint, proving that they are in the spectral gap of the operator. This, actually, is one of the reasons for which in [18] we looked for critical points of the action and not for the constrained critical points of the energy, in order to avoid such additional technical difficulties.

In view of all the above remarks, it is possible to state the following

Theorem 4.1** ([18, Theorem 2.12]).**

Let $\mathcal{G}$ satisfy (H1)-(H3), $p\in(2,6)$, $m>0$ and $\lambda<0$. Let also $(c_{n}),(\omega_{n})$ be two real sequences such that

[TABLE]

If $\{\psi_{n}=(\phi_{n},\chi_{n})^{T}\}$ is a bound state of frequency $\omega_{n}$ of the NLDE at relativistic parameter $c_{n}$, then, up to subsequences, there holds

[TABLE]

where $u$ is a bound state of frequency $\lambda$ of the NLSE.

Actually, the bound states of NLDE do not converge exactly to the bound states of the NLSE depicted in (7). In fact, $u$ is a solution of

[TABLE]

The coefficient $2m$ is consistent with the fact that in the nonrelativistic limit the kinetic (free) part of the hamiltonian of a particle is given by

[TABLE]

Remark 4.1*.*

It is also worth observing that Theorem 4.1, in contrast to Theorem 3.1, holds only for a fixed range of exponents, the $L^{2}$-subcritical case $p\in(2,6)$.