This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Oscillating solutions for nonlinear Helmholtz Equations

Rainer Mandel1, Eugenio Montefusco2, Benedetta Pellacci3 Institut für Analysis, Karlsruher Institut für Technologie, Englerstraße 2, 76131 Karlsruhe Rainer.Mandel@kit.edu Dipartimento di Matematica, ”Sapienza” Università di Roma, p.le Aldo Moro 5, 00185, Roma, Italy. eugenio.montefusco@uniroma1.it Dipartimento di Scienze e Tecnologie, Università di Napoli ”Parthenope”, Centro Direzionale, Isola C4 80143 Napoli, Italy. benedetta.pellacci@uniparthenope.it
Abstract.

Existence results for radially symmetric oscillating solutions for a class of nonlinear autonomous Helmholtz equations are given and their exact asymptotic behavior at infinity is established. Some generalizations to nonautonomous radial equations as well as existence results for nonradial solutions are found. Our theorems prove the existence of standing waves solutions of nonlinear Klein-Gordon or Schrödinger equations with large frequencies.

Key words and phrases:
Nonlinear Helmholtz equations, standing waves, oscillating solutions
1991 Mathematics Subject Classification:
35J05, 35J20, 35Q55.
11footnotetext: Research partially supported by the German Research Foundation (DFG) through the grant MA 6290/2-1 and CRC 1173.33footnotetext: Research partially supported by MIUR-PRIN project 2015KB9WPT0062015KB9WPT_{0}06 - PE1, “Gruppo Nazionale per l’Analisi Matematica, la Probabilità e le loro Applicazioni” (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM); University Project ”Sostegno alla ricerca individuale per il trienno 2016-2018”.

1. Introduction

The main aim of this paper is to give existence results for the following class of nonlinear equations

(1.1) Δu=g(u)in N-\Delta u=g(u)\qquad\text{in }\mathbb{R}^{N}

with N1N\geq 1 and assuming that the nonlinearity gg is such that

(1.2) gC1,σ()\displaystyle g\in C^{1,\sigma}(\mathbb{R})  for some σ(0,1),\displaystyle\;\text{ for some $\sigma\in(0,1)$},
(1.3) gg is odd,
(1.4) g(0)>0,\displaystyle\;g^{\prime}(0)>0,
(1.5) α0(0,+]:g is positive\displaystyle\exists\alpha_{0}\in(0,+\infty]\,:\;\text{$g$ is positive} on (0,α0)(0,\alpha_{0}) and negative on (α0,)(\alpha_{0},\infty).

There is a huge literature concerning (1.1) and nonautonomous variants of it under the assumption g(0)<0g^{\prime}(0)<0. Two seminal papers in this context are the contributions by Berestycki-Lions and Strauss [BL83, Str77] who proved the existence of smooth radially symmetric and exponentially decaying solutions for a large class of nonlinearities with this property. We refer to the monographs [AM07, Wil96] for more results in this context. One of the main interests in finding solutions of (1.1) is motivated by the fact that a solution uH1(N)u\in H^{1}(\mathbb{R}^{N}) of (1.1) gives rise to a standing wave, i.e. a solution of the form ψ(x,t)=eiλtu(x)\psi(x,t)=e^{i\lambda t}u(x), of the nonlinear time-dependent Klein-Gordon equation

2ψt2Δψ+V0ψ=f(ψ)(t,x)×N\dfrac{\partial^{2}\psi}{\partial t^{2}}-\Delta\psi+V_{0}\psi=f(\psi)\qquad(t,x)\in\mathbb{R}\times\mathbb{R}^{N}

with f(z)=g(|z|)z|z|f(z)=g(|z|)\tfrac{z}{|z|}. Therefore, the assumption (1.4) amounts to look for standing waves having low frequencies ω<V0\omega<V_{0} and numerous existence results for H1(N)H^{1}(\mathbb{R}^{N})-solutions under this assumption can be found in the references mentioned above. In this paper we deal with nonlinearities satisfying g(0)>0g^{\prime}(0)>0, which gives rise to standing waves with large frequencies ω>V0\omega>V_{0}. Looking at the form of the linearized operator Δg(0)-\Delta-g^{\prime}(0), one realizes that u0=0u_{0}=0 lies in its essential spectrum and we are actually dealing with a class of nonlinear Helmholtz equations. Furthermore, as explained in subsection 2.2 in [BL83], the hypothesis (1.4) has the striking consequence that radially symmetric H1(N)H^{1}(\mathbb{R}^{N}) solutions of (1.1) can not exist, and usual variational methods fail. On the other hand, (1.4) is naturally linked to (1.5); in particular, if g(z)/zg(z)/z decreases in (0,+)(0,+\infty), then (1.4) turns out to be necessary in order to have H1(N)H^{1}(\mathbb{R}^{N}) solutions. Actually, the relevant solutions naturally lie outside this functional space. This fact can also be illustrated by an examining the behaviour of the minimal energy solutions on a sequence of large bounded domains. Namely, in Theorem 4.4 we will show that if one takes a sequence of bounded domains Ωn\Omega_{n} invading N\mathbb{R}^{N}, then (1.4) guarantees the existence of a sequence (un)(u_{n}) of global minimizers of the associated action functional over H01(Ωn)H_{0}^{1}(\Omega_{n}) for sufficiently large nn. But, it results that (un)(u_{n}) converges in Cloc2(N)C^{2}_{{\rm loc}}(\mathbb{R}^{N}) to the constant solution uα0u\equiv\alpha_{0}.

Therefore, under the assumption (1.4), one has to look for solutions in a broader class of functions. Our focus will be on oscillating and localized ones which we define as follows.

Definition 1.1.

A distributional solution uC1,α(N)u\in C^{1,\alpha}(\mathbb{R}^{N}) of (1.1) is called oscillating if it has an unbounded sequence of zeros. It is called localized when it converges to zero at infinity.

Let us notice that, the Strong Maximum Principle implies that oscillating solutions of (1.1) change sign at each of their zeros; so that we are going to find solutions that change sign infinitely many times.

In our study, we will pay particular attention to the following model cases

(1.6) g1(z)=λz+zs+z2\displaystyle g_{1}(z)=-\lambda z+\dfrac{z}{s+z^{2}}\;\; where s>0,λ<1s.\displaystyle\text{ where }s>0,\,\lambda<\frac{1}{s}.
(1.7) g2(z)=k2z|z|p2z,\displaystyle g_{2}(z)=k^{2}z-|z|^{p-2}z,\;\qquad g3(z)=k2z+|z|p2zfor k0.\displaystyle g_{3}(z)=k^{2}z+|z|^{p-2}z\;\text{for $k\neq 0$}.

Our interest in these examples has various motivations. The nonlinearity g1g_{1} is related to the study of the propagation of lights beams in a photorefractive crystals (see [CCMS97, Yan04]) when a saturation effect is taken into account. Differently from the more frequently studied model

g~(z):=λz+z31+sz2,\tilde{g}(z):=-\lambda z+\dfrac{z^{3}}{1+sz^{2}},

see e.g. [CBWM04], g1g_{1} describes a transition from the linear propagation and the saturated one. This difference has important consequences, for instance for g=g~g=\tilde{g} there are H1(N)H^{1}(\mathbb{R}^{N}) solutions of (1.1) (e.g. see Theorem 3.6 in [SZ99]), whereas, as we have already observed, this is not the case if g=g1g=g_{1} due to g1(0)>0g_{1}^{\prime}(0)>0. Notice that, as λ<1/s\lambda<1/s, equation (1.1) for g=g1g=g_{1} can be rewritten in the following form

(1.8) Δuk2u=u3s(s+u2)in Nwith k2=1sλ-\Delta u-k^{2}u=-\frac{u^{3}}{s(s+u^{2})}\quad\text{in }\mathbb{R}^{N}\quad\text{with $k^{2}=\frac{1}{s}-\lambda$}

which allows to settle the problem in H1(N)H^{1}(\mathbb{R}^{N}) in every dimension NN and that shows that also this saturable model is included in the class of the nonlinear Helmholtz equations. The principal difference between (1.7) and (1.6) is that the formers are superlinear and homogeneous nonlinearities, while the latter is not homogeneous and it is asymptotically linear. However, all of them satisfy our general assumptions, with α0(0,+)\alpha_{0}\in(0,+\infty) for g1g_{1} and g2g_{2}, and α0=+\alpha_{0}=+\infty for g3g_{3}.

Up to now nonlinear Helmholtz equations (1.1) have been mainly investigated for the model nonlinearity g3g_{3} or more general superlinear nonlinearities, even not autonomous. In a series of papers [EW14, EW15, Eve15, EW] Evéquoz and Weth proved the existence of radial and nonradial real, localized solutions of this equation under various different assumptions on the nonlinearity. Let us mention that some of the tools used in [EW15] had already appeared in a paper by Gutiérrez [Gut04] where the existence of complex-valued solutions was proved for space dimensions N=3,4N=3,4. Let us first focus our attention on radially symmetric solutions and state our first result, which provides a complete description of the radially symmetric solutions of (1.1).

Theorem 1.2.

Assume (1.2),(1.3),(1.4),(1.5). Then there is a continuum 𝒞={uαC2(N):|α|<α0}\mathcal{C}=\{u_{\alpha}\in C^{2}(\mathbb{R}^{N}):|\alpha|<\alpha_{0}\} in C2(N)C^{2}(\mathbb{R}^{N}) consisting of radially symmetric oscillating solutions of (1.1) having the following properties for all |α|<α0|\alpha|<\alpha_{0}:

  • (i)

    uα(0)=αu_{\alpha}(0)=\alpha,

  • (ii)

    uαL(N)=|α|,uαL(N)2G(α)\|u_{\alpha}\|_{L^{\infty}(\mathbb{R}^{N})}=|\alpha|,\|u_{\alpha}^{\prime}\|_{L^{\infty}(\mathbb{R}^{N})}\leq\sqrt{2G(\alpha)}.

Moreover, for N=1N=1 all these solutions are periodic; whereas, for N2N\geq 2 they are localized and satisfy the following asymptotic behavior:

  • (iii)

    There are positive numbers cα,Cα>0c_{\alpha},C_{\alpha}>0 such that

    cαr(1N)/2|uα(r)|+|uα(r)|+|uα′′(r)|Cαr(1N)/2for all r1.c_{\alpha}r^{(1-N)/2}\leq|u_{\alpha}(r)|+|u_{\alpha}^{\prime}(r)|+|u_{\alpha}^{\prime\prime}(r)|\leq C_{\alpha}r^{(1-N)/2}\quad\text{for all }r\geq 1.

Here a continuum in C2(N)C^{2}(\mathbb{R}^{N}) is a connected subset of C2(N)C^{2}(\mathbb{R}^{N}) with respect to the uniform convergence of the zeroth, first and second derivatives. The continuum 𝒞\mathcal{C} found in Theorem 1.2 is even maximal in the sense that there are no further radially symmetric localized solutions as we will see in section 2. Moreover, conclusion (iii)(iii) states that the property uαLs(N)u_{\alpha}\in L^{s}(\mathbb{R}^{N}) is equivalent to uαW2,s(N)u_{\alpha}\in W^{2,s}(\mathbb{R}^{N}), and this happens if and only if s>2NN1s>\frac{2N}{N-1}. Notice that this implies that uαL2(N)u_{\alpha}\notin L^{2}(\mathbb{R}^{N}), showing again that the solutions, as expected, live outside the commonly used energy space.

Furthermore, let us stress that the behaviour of the nonlinearity beyond α0\alpha_{0} is completely irrelevant, in particular, the negativity of gg on (α0,)(\alpha_{0},\infty) is actually not needed. This is the reason why we do not need to assume any subcritical growth condition on the exponent pp in the model nonlinearities g2,g3g_{2},g_{3}. Let us recall that, in the autonomous setting, Theorem 4 in [EW14] yields nontrivial radially symmetric solutions of (1.1) for superlinear nonlinearities, so that their results hold for the nonlinearity g3g_{3}, but not for g1,g2g_{1},g_{2}.

Theorem 1.2 admits generalizations to some nonautonomous radially symmetric nonlinearities. In particular we can prove a nonautonomous version of this result that applies to the nonlinearities

(1.9) g1(r,z)\displaystyle g_{1}(r,z) =λ(r)z+zs(r)+z2,\displaystyle=-\lambda(r)z+\frac{z}{s(r)+z^{2}},
(1.10) g2(r,z)\displaystyle g_{2}(r,z) =k(r)2z±Q(r)|z|p2z,\displaystyle=k(r)^{2}z\pm Q(r)|z|^{p-2}z,

under suitable assumptions on the coefficients λ(r),s(r),k(r),Q(r)\lambda(r),\,s(r),\,k(r),\,Q(r), see Theorem 2.10 and the Corollaries 2.11 and 2.13. Our results in this context extend Theorem 4 in [EW14] in several directions (see Remark 2.12).

The existence of non-radially symmetric solutions is clearly a more difficult topic and here we can give a partial positive answer in this direction, by exploiting the argument developed in [EW14, EW15, Eve15, EW], where the authors study the equation

(1.11) Δuk2u=f(x,u)in N,-\Delta u-k^{2}u=f(x,u)\quad\text{in }\mathbb{R}^{N},

where f(x,u)f(x,u) is a super-linear nonlinearity satisfying suitable hypotheses, that include, for example, a non-autonomous generalization of our model nonlinearity g3g_{3}, i.e.

f(x,u)=Q(x)|u|p2u, with Q(x)>0.f(x,u)=Q(x)|u|^{p-2}u,\quad\text{ with }\quad Q(x)>0.

Among other results, in [EW15] (Theorem 1.1 and Theorem 1.2) it is shown that if QQ is N\mathbb{Z}^{N}-periodic or vanishing at infinity then there exist nontrivial solutions of (1.11) for pp satisfying 2(N+1)N1<p<2NN2\frac{2(N+1)}{N-1}<p<\frac{2N}{N-2} when N3N\geq 3.

Our contribution to this issue is that the positivity assumption on QQ may be replaced by a negativity assumption in order to make the dual variational approach work, so that, using Fourier transform we show that the main ideas from [EW15] may be modified in such a way that their main results remain true for negative QQ. Our results read as follows.

Theorem 1.3.

Let N3,2(N+1)N1<p<2NN2N\geq 3,\frac{2(N+1)}{N-1}<p<\frac{2N}{N-2} and let QL(N)Q\in L^{\infty}(\mathbb{R}^{N}) be periodic and negative almost everywhere. Then the equation (1.11) has a nontrivial localized oscillating strong solution in W2,q(N)C1,α(N)W^{2,q}(\mathbb{R}^{N})\cap C^{1,\alpha}(\mathbb{R}^{N}) for all q[p,),α(0,1)q\in[p,\infty),\alpha\in(0,1).

Theorem 1.4.

Let N3,2(N+1)N1<p<2NN2N\geq 3,\frac{2(N+1)}{N-1}<p<\frac{2N}{N-2} and let QL(N)Q\in L^{\infty}(\mathbb{R}^{N}) be negative almost everywhere with Q(x)0Q(x)\to 0 as |x||x|\to\infty. Then the equation (1.11) has a sequence of pairs ±um\pm u_{m} of nontrivial localized oscillating strong solutions in W2,q(N)C1,α(N)W^{2,q}(\mathbb{R}^{N})\cap C^{1,\alpha}(\mathbb{R}^{N}) for all q[p,),α(0,1)q\in[p,\infty),\alpha\in(0,1) such that

umLp(N)as m.\|u_{m}\|_{L^{p}(\mathbb{R}^{N})}\to\infty\quad\text{as }m\to\infty.

Since the above results together with those from [EW15] provide some existence results for the Nonlinear Helmholtz equation associated with the nonlinearity g2g_{2} from (1.10), one is lead to wonder whether similar results hold true for asymptotically linear nonlinearities like g1g_{1} in (1.9). Here, the dual variational framework does not seem to be convenient since even the choice of the appropriate function spaces is not clear. A thorough discussion of such nonlinear Helmholtz equations leading to existence results for nonradial solutions still remains to be done.

Let us observe that there is a gap in the admissible range of exponent between Theorem 1.2 and Theorem 1.3. Reading Theorem 1.2 one is naturally lead to the conjecture that nontrivial nonradial solutions in Lp(N)L^{p}(\mathbb{R}^{N}) may be found regardless of any sign condition on QQ and for all exponents p>2NN1p>\frac{2N}{N-1}. On the contrary, Theorem 1.3 only holds for exponents p>2(N+1)N1p>\frac{2(N+1)}{N-1}, so that it is still an open question whether or not nonradial LpL^{p}-solutions exist for p(2NN1,2(N+1)N1]p\in(\frac{2N}{N-1},\frac{2(N+1)}{N-1}].

The paper is organized as follows: In section 2 we present the proof of Theorem 1.2 as well as a generalization to the radial nonautonomous case (see Theorem 2.10 and Corollaries 2.11, 2.13). In section 3 we present the proofs of Theorem 1.3 and Theorem 1.4. In section 4 we will discuss in detail the attempt to obtain a solution by approximating N\mathbb{R}^{N} by bounded domains.

2. Radial solutions

2.1. The autonomous case

Throughout this section we will suppose that (1.2), (1.3), (1.4), (1.5) hold true. We will prove Theorem 1.2 by providing a complete understanding of the initial value problem

(2.1) u′′N1ru=g(u)in (0,),u(0)=α,u(0)=0-u^{\prime\prime}-\dfrac{N-1}{r}u^{\prime}=g(u)\quad\text{in }(0,\infty),\qquad u(0)=\alpha,\;u^{\prime}(0)=0

for α\alpha\in\mathbb{R} and NN\in\mathbb{N}. Notice that our assumptions on gg require that there exists a δ>0\delta>0 such that

g(z)z>0z(δ,δ).g(z)z>0\quad\forall\,z\in(-\delta,\delta).

Such a positivity region is in fact almost necessary as the following result shows.

Proposition 2.1.

Assume that gC()g\in C(\mathbb{R}) satisfies g(z)z<0g(z)z<0 for all zz\in\mathbb{R}. Then there is no nontrivial localized solution and there is no nontrivial oscillating solution uC2(N)u\in C^{2}(\mathbb{R}^{N}) of (1.1).

Proof.

Assume that uC2(N)u\in C^{2}(\mathbb{R}^{N}) is a nontrivial localized or oscillating solution. Then it attains a positive local maximum or a negative local minimum in some point x0Nx_{0}\in\mathbb{R}^{N}. Hence we obtain

Δu(x0)u(x0)=u(x0)g(u(x0))<0,-\Delta u(x_{0})u(x_{0})=u(x_{0})g(u(x_{0}))<0,

a contradiction. \Box

Remark 2.2.

In view of elliptic regularity theory the above result is also true for weak solutions uH1(N)u\in H^{1}(\mathbb{R}^{N}) since these solutions coincide almost everywhere with classical solutions and decay to zero at infinity by Theorem C.3 in [Sim82]. Notice that in case N3N\geq 3 we can deduce the non-existence of H1(N)H^{1}(\mathbb{R}^{N}) solutions from the fact that g(z)z<0g(z)z<0 in \mathbb{R} violates the necessary condition (1.3) in [BL83], see section 2.2 in that paper. In the case N=2N=2 the same follows from Remarque 1 in [BGK83].

First we briefly address the one-dimensional initial value problem

(2.2) u′′=g(u)in (0,),u(0)=α,u(0)=0.-u^{\prime\prime}=g(u)\quad\text{in }(0,\infty),\qquad u(0)=\alpha,\;u^{\prime}(0)=0.

In view of the oddness of gg it suffices to discuss the inital value problem for α0\alpha\geq 0. The uniquely determined solution of the initial value problem will be denoted by uαu_{\alpha} with maximal existence interval (Tα,Tα)(-T_{\alpha},T_{\alpha}) for Tα(0,]T_{\alpha}\in(0,\infty].

Proposition 2.3.

Let N=1N=1. Then the following holds:

  • (i)

    If α=α0\alpha=\alpha_{0}\in\mathbb{R} then uαα0u_{\alpha}\equiv\alpha_{0} and if α=0\alpha=0 then uα0u_{\alpha}\equiv 0.

  • (ii)

    If α>α0\alpha>\alpha_{0} then uαu_{\alpha} strictly increases to ++\infty on (0,Tα)(0,T_{\alpha}).

  • (iii)

    If 0<α<α00<\alpha<\alpha_{0} then uαu_{\alpha} is periodic and oscillating with uα=α\|u_{\alpha}\|_{\infty}=\alpha.

Proof.

Conclusion (i) immediately follows from (1.5). Then we only have to prove (ii) and (iii). For notational convenience we write u,Tu,T instead of uα,Tαu_{\alpha},T_{\alpha}. In the situation of (ii) we set ξ:=sup{s[0,T):u′′(s)>0}\xi:=\sup\{s\in[0,T):u^{\prime\prime}(s)>0\}. From u(0)=α>α0u(0)=\alpha>\alpha_{0} and (2.2) we get u′′(0)=g(u(0))=g(α)>0u^{\prime\prime}(0)=-g(u(0))=-g(\alpha)>0 and thus ξ(0,T]\xi\in(0,T]. We even have ξ=T\xi=T, because otherwise

u(ξ)=α+0ξ0tu′′(s)𝑑s𝑑t>α>α0u(\xi)=\alpha+\int_{0}^{\xi}\int_{0}^{t}u^{\prime\prime}(s)\,ds\,dt>\alpha>\alpha_{0}

and thus u′′(ξ)>0u^{\prime\prime}(\xi)>0 in view of assumption (1.5) and (2.2). This, however, would contradict that ξ\xi is the supremum, hence ξ=T\xi=T. As a consequence, uu is strictly convex on (0,Tα)(0,T_{\alpha}) which implies (ii).

In order to show (iii) we notice that (1.3) implies that solutions are symmetric about critical points and antisymmetric about zeros. Therefore, it suffices to show that uu decreases until it attains a zero. By the choice of α(0,α0)\alpha\in(0,\alpha_{0}) we have u′′(0)<0u^{\prime\prime}(0)<0 so that uu decreases on a right neighbourhood of 0. Exploiting (1.5) and (2.2) we deduce that u′′(s)u^{\prime\prime}(s) is negative whenever 0<u(s)<u(0)<α00<u(s)<u(0)<\alpha_{0}. As a consequence, we obtain that uu decreases as long as it remains positive. Moreover, it cannot be positive on [0,)[0,\infty) since this would imply, thanks 0u(r)α<α00\leq u(r)\leq\alpha<\alpha_{0} and the assumptions (1.4),(1.5),

u′′(r)+c(r)u(r)=0,with c(r):=g(u(r))u(r)c0>0.u^{\prime\prime}(r)+c(r)u(r)=0,\quad\text{with }c(r):=\frac{g(u(r))}{u(r)}\geq c_{0}>0.

Hence, Sturm’s comparison theorem (p.2 in [Swa68]) ensures that uu vanishes somewhere, so that it cannot be positive in [0,+)[0,+\infty), a contradiction. Hence, uu attains a zero and the proof is finished. \Box

Next, we consider the initial value problem (2.1) in the higher dimensional case N2N\geq 2. Again, we may restrict our attention to the case α0\alpha\geq 0 and we will denote with G(z)G(z) the primitive of the function g(s)g(s), such that G(0)=0G(0)=0. The following result furnishes the study of the solution set which are needed in the proof of Again, the uniquely determined solution of the initial value problem (2.1) will be denoted by uαu_{\alpha} with maximal existence interval (Tα,Tα)(-T_{\alpha},T_{\alpha}).

Remark 2.4.

There are many contributions concerning (1.1) in dimension N=1N=1, mainly related to some resonance phenomena. In this context, some “Landesman-Lazer” type conditions, joint with suitable hypotheses on the nonlinearity gg, are assumed in oder to obtain existence of bounded, periodic or oscillating solution, eventually with arbitrarily large LL^{\infty} norm, by taking advantage of the presence on a forcing term in the equation (see [SV14, Ver03] and the references therein). Here the situation is different, as we do not need any monotonicity assumption on gg, nor the knowledge of the asymptotic behavior at infinity of gg is important, as it is in [SV14, Ver03]. Moreover, our solutions satisfy a uniform LL^{\infty} bound, so that the phenomenon we are dealing with is actually different from the resonant one.

Lemma 2.5.

Let N2N\geq 2. Then the following holds:

  • (i)

    If α=α0\alpha=\alpha_{0}\in\mathbb{R} then uαα0u_{\alpha}\equiv\alpha_{0} and if α=0\alpha=0 then uα0u_{\alpha}\equiv 0.

  • (ii)

    If α>α0\alpha>\alpha_{0} then uαu_{\alpha} strictly increases to ++\infty on [0,Tα)[0,T_{\alpha}).

  • (iii)

    If 0<α<α00<\alpha<\alpha_{0} then uαu_{\alpha} is oscillating, localized and satisfies

    (2.3) uαL()=|α|anduαL()2G(α)\|u_{\alpha}\|_{L^{\infty}(\mathbb{R})}=|\alpha|\quad\text{and}\quad\|u_{\alpha}^{\prime}\|_{L^{\infty}(\mathbb{R})}\leq\sqrt{2G(\alpha)}

    as well as

    (2.4) cαr(1N)/2|uα(r)|+|uα(r)|+|uα′′(r)|Cαr(1N)/2for r1c_{\alpha}r^{(1-N)/2}\leq|u_{\alpha}(r)|+|u_{\alpha}^{\prime}(r)|+|u_{\alpha}^{\prime\prime}(r)|\leq C_{\alpha}r^{(1-N)/2}\quad\text{for }r\geq 1

    for some cα,Cα>0c_{\alpha},C_{\alpha}>0 depending on the solution but not on rr.

Proof.

The existence and uniqueness of a twice continuously differentiable solution uα:(Tα,Tα)u_{\alpha}:(-T_{\alpha},T_{\alpha})\to\mathbb{R} can be deduced from Theorem 1 and Theorem 2 in [RW97]. We write again u,Tu,T in place of uα,Tαu_{\alpha},T_{\alpha}. The proof of (i) is direct and assertion (ii) follows similar to the one-dimensional case. Indeed, note that u′′(0)>0u^{\prime\prime}(0)>0 because of

Nu′′(0)=limr0+u′′(r)+N1ru(r)=g(u(0))=g(α)>0.\displaystyle Nu^{\prime\prime}(0)=\lim_{r\to 0^{+}}u^{\prime\prime}(r)+\tfrac{N-1}{r}u^{\prime}(r)=-g(u(0))=-g(\alpha)>0.

Then, letting ξ:=sup{s(0,T):u(s)>0}\xi:=\sup\{s\in(0,T):u^{\prime}(s)>0\}, it results ξ(0,T]\xi\in(0,T]. Assuming by contradiction that ξ<T\xi<T and using that α>α0\alpha>\alpha_{0}, from (1.5) we obtain

ξN1u(ξ)=0ξtN1g(u(t))𝑑t>0\xi^{N-1}u^{\prime}(\xi)=-\int_{0}^{\xi}t^{N-1}g(u(t))\,dt>0

which is impossible, i.e. ξ=T\xi=T. Then, (2.1),(1.5) and the maximality of TT yield (ii). The proof of (iii) is lengthy so that it will be subdivided into four steps.

Step 1: uu decreases to a first zero.  For all r>0r>0 such that 0<u<α00<u<\alpha_{0} on [0,r][0,r] we have

rN1u(r)=0rtN1g(u(t))𝑑t<0,r^{N-1}u^{\prime}(r)=-\int_{0}^{r}t^{N-1}g(u(t))\,dt<0,

showing that uu decreases as long as it remains positive, as in the one-dimensional case. Moreover, the function uu can not remain positive on [0,)[0,\infty) because otherwise v(r):=r(N1)/2u(r)v(r):=r^{(N-1)/2}u(r) would be a positive solution of

(2.5) v′′+c(r)v=0where c(r)=g(u(r))u(r)(N1)(N3)4r2.v^{\prime\prime}+c(r)v=0\qquad\text{where }c(r)=\frac{g(u(r))}{u(r)}-\frac{(N-1)(N-3)}{4r^{2}}.

As in the proof of Proposition 2.3 we observe c(r)c0>0c(r)\geq c_{0}>0 for sufficiently large rr so that Sturm’s comparison theorem tells us that vv vanishes somewhere. This is a contradiction to the positivity of uu and thus uu attains a first zero.

Step 2: uu oscillates and satisfies (2.3).  Let us first show that there are 0=r0<r1<r2<r3<0=r_{0}<r_{1}<r_{2}<r_{3}<\ldots such that all r4jr_{4j} are local maximizers, all r4j+2r_{4j+2} are local minimizers and all r2j+1r_{2j+1} are zeros of uu. Moreover, we will find that all zeros or critical points of uu are elements of this sequence and

(2.6) 2G(u(r0))>u(r1)2>2G(u(r2))>u(r3)2>2G(u(r4))>\displaystyle 2G(u(r_{0}))>u^{\prime}(r_{1})^{2}>2G(u(r_{2}))>u^{\prime}(r_{3})^{2}>2G(u(r_{4}))>\ldots

In order to prove this we consider the function

(2.7) Z(r):=u(r)2+2G(u(r)).Z(r):=u^{\prime}(r)^{2}+2G(u(r)).

and we observe that ZZ decreases as

(2.8) Z(r)=2u(r)(u′′(r)+g(u(r)))=2(N1)ru(r)2<0.\displaystyle Z^{\prime}(r)=2u^{\prime}(r)(u^{\prime\prime}(r)+g(u(r)))=-\frac{2(N-1)}{r}u^{\prime}(r)^{2}<0.

The existence of a first zero r1>0=r0r_{1}>0=r_{0} of uu has been shown in Step 1 and the strict monotonicity of ZZ implies Z(r1)<Z(r0)Z(r_{1})<Z(r_{0}). Concerning the behaviour of uu on [r1,)[r_{1},\infty) there are now three alternatives:

  • (a)

    uu decreases until it attains u(r0)-u(r_{0})

  • (b)

    uu decreases on [r1,)[r_{1},\infty) to some value u[u(r0),0)u_{\infty}\in[-u(r_{0}),0)

  • (c)

    uu decreases until it attains a critical point at some r2>r1r_{2}>r_{1} with u(r0)<u(r2)<0-u(r_{0})<u(r_{2})<0.

Let us show that the cases (a) and (b) do not occur. Indeed, if there exists r>r0r>r_{0} such that u(r)=u(r0)u(r)=-u(r_{0}), then, by (2.7) we deduce that

Z(r)2G(u(r))=2G(u(r0))=Z(r0)Z(r)\geq 2G(u(r))=2G(u(r_{0}))=Z(r_{0})

which is forbidden by (2.8). Then, in particular (2.3) holds. Hence, the case (a) is impossible. Let us now suppose that (b) holds. Then uu_{\infty} has to be a stationary solution of (2.1) and thus u=α0=u(r0)u_{\infty}=-\alpha_{0}=-u(r_{0}). But then

Z(r)2G(u(r))2G(u)=2G(u(r0))=Z(r0)as rZ(r)\geq 2G(u(r))\to 2G(u_{\infty})=2G(-u(r_{0}))=Z(r_{0})\quad\text{as }r\to\infty

which again contradicts (2.8). So the case (c) occurs and there must be a critical point r2r_{2} with

2G(u(r2))=Z(r2)<Z(r1)=u(r1)2<Z(r0)=2G(u(r0)),2G(u(r_{2}))=Z(r_{2})<Z(r_{1})=u^{\prime}(r_{1})^{2}<Z(r_{0})=2G(u(r_{0})),

so that (2.1), (1.5) and (1.3) yield

0>u(r2)>u(r0) and u(r2)=0,u′′(r2)>0.0>u(r_{2})>-u(r_{0})\text{ and }u^{\prime}(r_{2})=0,u^{\prime\prime}(r_{2})>0.

Hence, r2r_{2} is a local minimizer. Using that ZZ is decreasing we can now repeat the argument to get a zero r3>r2r_{3}>r_{2}, a local maximizer r4>r3r_{4}>r_{3}, a zero r5>r4r_{5}>r_{4} and so on. By the strict monotonicity of ZZ one obtains (2.6) and thus (2.3). Notice that this reasoning also shows that there are no further zeros or critical points.

Step 3: uu is localized.   First we show u(r)0u(r)\to 0 as rr\to\infty. Our proof is similar to the one of Lemma 4.1 in [GZ08] and it will be presented for the convenience of the reader. Take the sequence of maximizers {r4j}\{r_{4j}\} and assume by contradiction that u(r4j)z(0,α0)u(r_{4j})\to z\in(0,\alpha_{0}). Then (2.1) and Ascoli-Arzelà Theorem imply that u(+r4j)u(\cdot+r_{4j}) converges locally uniformly to the unique solution ww of (2.2) with w(0)=z,w(0)=0w(0)=z,w^{\prime}(0)=0. Proposition 2.3, (iii) implies that this solution ww is TT-periodic with two zeroes at T/4,3T/4T/4,3T/4. As a consequence, there exists δ>0\delta>0 such that |w|22δ|w^{\prime}|^{2}\geq 2\delta on [T/42δ,T/4+2δ][T/4-2\delta,T/4+2\delta]. Hence, for sufficiently large j0j_{0}\in\mathbb{N} we have for jj0j\geq j_{0}

u(r4j+r)2δfor r[T/4δ,T/4+δ] and r4(j+1)r4j\displaystyle u^{\prime}(r_{4j}+r)^{2}\geq\delta\quad\text{for $r\in[T/4-\delta,T/4+\delta]$ and }\quad r_{4(j+1)}-r_{4j} Tδ.\displaystyle\geq T-\delta.

From this we deduce for jj0j\geq j_{0}

(2.9) u(r)2\displaystyle u^{\prime}(r)^{2} δ,for r[r4j+T/4δ,r4j+T/4+δ],\displaystyle\geq\delta,\qquad\text{for }r\in[r_{4j}+T/4-\delta,r_{4j}+T/4+\delta],
(2.10) r4j\displaystyle r_{4j} r4j0+(jj0)(Tδ)for jj0.\displaystyle\geq r_{4j_{0}}+(j-j_{0})(T-\delta)\qquad\text{for }j\geq j_{0}.

Then, for kj0k\geq j_{0} and r>r4k+T/4+δr>r_{4k}+T/4+\delta we may exploit (2.8) and (2.9) to obtain

Z(r)\displaystyle Z(r) =Z(0)2(N1)0ru(t)2t𝑑t\displaystyle=Z(0)-2(N-1)\int_{0}^{r}\frac{u^{\prime}(t)^{2}}{t}\,dt
Z(0)2(N1)j=j0kr4j+T/4δr4j+T/4+δu(t)2t𝑑t\displaystyle\leq Z(0)-2(N-1)\sum_{j=j_{0}}^{k}\int_{r_{4j}+T/4-\delta}^{r_{4j}+T/4+\delta}\frac{u^{\prime}(t)^{2}}{t}\,dt
Z(0)2(N1)δ2j=j0kr4j+T/4δr4j+T/4+δ1t𝑑t\displaystyle\leq Z(0)-2(N-1)\delta^{2}\sum_{j=j_{0}}^{k}\int_{r_{4j}+T/4-\delta}^{r_{4j}+T/4+\delta}\frac{1}{t}\,dt
=Z(0)2(N1)δ2j=j0kln(r4j+T4+δr4j+T4δ).\displaystyle=Z(0)-2(N-1)\delta^{2}\sum_{j=j_{0}}^{k}\ln\left(\dfrac{r_{4j}+\frac{T}{4}+\delta}{r_{4j}+\frac{T}{4}-\delta}\right).

Let us fix c(δ)>0c(\delta)>0 such that ln(1+x)c(δ)x\ln(1+x)\geq c(\delta)x for 0x2δ/(r4j0+T4δ)0\leq x\leq 2\delta/(r_{4j_{0}}+\frac{T}{4}-\delta). Then (2.10) implies

Z(r)\displaystyle Z(r) Z(0)2(N1)δ2c(δ)j=j0k2δr4j+T4δ\displaystyle\leq Z(0)-2(N-1)\delta^{2}c(\delta)\sum_{j=j_{0}}^{k}\dfrac{2\delta}{r_{4j}+\frac{T}{4}-\delta}
Z(0)2(N1)δ2c(δ)j=j0k2δr4j0+(jj0)(Tδ)+T4δ.\displaystyle\leq Z(0)-2(N-1)\delta^{2}c(\delta)\sum_{j=j_{0}}^{k}\dfrac{2\delta}{r_{4j_{0}}+(j-j_{0})(T-\delta)+\frac{T}{4}-\delta}.

Choosing now k,rk,r sufficiently large we obtain that Z(r)Z(r)\to-\infty because the harmonic series diverges, but (2.3)implies that Z(r)2G(u(r))0Z(r)\geq 2G(u(r))\geq 0, yielding a contradiciton. As a consequence, u(r4j)u(r_{4j}) converges to zero as jj\to\infty and analogously we deduce that also u(r4j+2)0u(r_{4j+2})\to 0. In the end, we obtain u(r)0u(r)\to 0 as r+r\to+\infty.
Since ZZ is decreasing and nonnegative it follows that Z(r)Z[0,Z(0))Z(r)\to Z_{\infty}\in[0,Z(0)) as rr\to\infty. Hence, by (2.7), also |u||u^{\prime}| has a limit at infinity which must be zero because uu converges to 0. Finally, from the differential equation we deduce that u′′(r)0u^{\prime\prime}(r)\to 0 as rr\to\infty, i.e.

(2.11) u(r),u(r),u′′(r)0(r).u(r),u^{\prime}(r),u^{\prime\prime}(r)\to 0\quad(r\to\infty).

As in Lemma 4.2 in [GZ08] we get that for any ε>0\varepsilon>0 there exists Cε>0C_{\varepsilon}>0 such that

(2.12) |u(r)|,|u(r)|,|u′′(r)|Cεr1N2+ε(r1).|u(r)|,|u^{\prime}(r)|,|u^{\prime\prime}(r)|\leq C_{\varepsilon}r^{\frac{1-N}{2}+\varepsilon}\qquad(r\geq 1).

Step 4: Proof of (2.4).  Slightly generalizing the approach from the proof of Theorem 4 in [EW14] we study the function

(2.13) ψ(r):=v(r)2+2rN1G(u(r)),where v(r)=r(N1)/2u.\psi(r):=v^{\prime}(r)^{2}+2r^{N-1}G(u(r)),\quad\text{where }\;v(r)=r^{(N-1)/2}u.

Using the function cc from (2.5) and taking into account (2.1), we obtain that ψ\psi satisfies the following differential equation

ψ(r)\displaystyle\psi^{\prime}(r) =2v(r)[c(r)v(r)]+2r(N1)/2g(u(r))[v(r)N12r(N3)/2u(r)]\displaystyle=2v^{\prime}(r)\left[-c(r)v(r)\right]+2r^{(N-1)/2}g(u(r))\left[v^{\prime}(r)-\frac{N-1}{2}r^{(N-3)/2}u(r)\right]
+2(N1)rN2G(u(r))\displaystyle\quad+2(N-1)r^{N-2}G(u(r))
=(N1)rN2(2G(u(r))u(r)g(u(r)))+(N1)(N3)2r2v(r)v(r).\displaystyle=(N-1)r^{N-2}\left(2G(u(r))-u(r)g(u(r))\right)+\frac{(N-1)(N-3)}{2r^{2}}v(r)v^{\prime}(r).

Taking into account (2.11) and using (1.4) and (1.5) we obtain that there exist C,r0(0,+)C,\,r_{0}\in(0,+\infty) such that

(2.14) u(r)2G(u(r))Crr0.\dfrac{u(r)^{2}}{G(u(r))}\leq C\qquad\forall\,r\geq r_{0}.

Then, exploiting (1.2) and (2.12), we find positive numbers C,C′′,rC^{\prime},C^{\prime\prime},r^{*} such that, for all rrr\geq r^{*}, it results

|(N1)rN2(2G(u(r))u(r)g(u(r)))|(N1)C2r|2G(u(r))u(r)g(u(r))|u(r)22rN1G(u(r))Cr|u(r)|σψ(r)C′′r1+(1N2+ε)σψ(r).\displaystyle\begin{aligned} &\hskip-28.45274pt\big{|}(N-1)r^{N-2}\left(2G(u(r))-u(r)g(u(r))\right)\big{|}\\ &\leq\dfrac{(N-1)C}{2r}\frac{|2G(u(r))-u(r)g(u(r))|}{u(r)^{2}}\cdot 2r^{N-1}G(u(r))\\ &\leq\frac{C^{\prime}}{r}|u(r)|^{\sigma}\psi(r)\\ &\leq C^{\prime\prime}r^{-1+(\frac{1-N}{2}+\varepsilon)\sigma}\psi(r).\end{aligned}

Moreover, using (2.13) and (2.14), we get

|(N1)(N3)2r2v(r)v(r)|\displaystyle\Big{|}\frac{(N-1)(N-3)}{2r^{2}}v(r)v^{\prime}(r)\Big{|} |(N1)(N3)|r2(v(r)2+v(r)2)\displaystyle\leq\frac{|(N-1)(N-3)|}{r^{2}}\cdot(v(r)^{2}+v^{\prime}(r)^{2})
|(N1)(N3)r2(CrN1G(u(r))+v(r)2)\displaystyle\leq\frac{|(N-1)(N-3)}{r^{2}}\cdot(Cr^{N-1}G(u(r))+v^{\prime}(r)^{2})
|(N1)(N3)|(C+1)r2ψ(r).\displaystyle\leq\frac{|(N-1)(N-3)|(C+1)}{r^{2}}\cdot\psi(r).

This yields |ψ(r)|a(r)ψ(r)\left|\psi^{\prime}(r)\right|\leq a(r)\psi(r) for rrr\geq r^{*} and some positive integrable function aa. Dividing this inequality by the positive function ψ(r)\psi(r) and integrating the resulting inequality over [r,)[r^{*},\infty) shows that ψ\psi is bounded from below and from above by a positive number. From this we obtain the lower and upper bounds (2.4) and the proof is finished. \Box

We are now ready to give the proof of Theorem 1.2.

Proof of Theorem 1.2  Let us define the set

𝒞={uα(||)C2(N):|α|<α0}\mathcal{C}=\{u_{\alpha}(|\cdot|)\in C^{2}(\mathbb{R}^{N}):|\alpha|<\alpha_{0}\}

where uαu_{\alpha} denotes the unique solution of the initial value problem (2.1). The set 𝒞\mathcal{C} is a subset of C2(N)C^{2}(\mathbb{R}^{N}), and it is a continuum thanks to the Ascoli-Arzelà Theorem. From Lemma 2.5 we obtain that all elements of 𝒞\mathcal{C} are oscillating localized solutions satisfying (2.3) and (2.4). \Box

Remark 2.6.

Let us mention that an analogous result to Theorem 1.2 in Theorem 1 [GLZ09] and it is applied to a more restrictive class of nonlinearities. Moreover, the above theorem is related to Theorem 4 in [EW14] but we do not need their assumption (g2)(g2). Actually, this hypothesis is not satisfied in our model cases g=g1g=g_{1} or g=g2g=g_{2}.

Remark 2.7.

The arguments from the proof of Theorem 1.2 also show the existence of oscillating localized solutions to initial value problems which are not of nonlinear Helmholtz type. For instance, one can treat concave-convex problems such as

(2.15) Δu=λ|u|q2u+μ|u|p2uin N,-\Delta u=\lambda|u|^{q-2}u+\mu|u|^{p-2}u\quad\text{in }\mathbb{R}^{N},

for 1<q<2<p<1<q<2<p<\infty with λ>0,μ\lambda>0,\mu\in\mathbb{R}, see for instance [ABC94] or [BEP95] for corresponding results on a bounded domain with homogeneous Dirichlet boundary conditions. The existence of solutions is provided by Theorem 1 in [RW97] so that the steps 1,2,3 are proven in the same way as above and we obtain infinitely many radially symmetric, oscillating, localized, solutions of (2.15).

Remark 2.8.

Using nonlinear oscillation theorems instead of Sturm’s comparison theorem we can even extend the above observation towards superlinear nonlinearities gg satisfying λ|z|qg(z)zΛ|z|q\lambda|z|^{q}\leq g(z)z\leq\Lambda|z|^{q} for

2<q2(N+1)N1,N{1,2,3}or2<q2(N1)N2,N4.2<q\leq\frac{2(N+1)}{N-1},N\in\{1,2,3\}\quad\text{or}\quad 2<q\leq\frac{2(N-1)}{N-2},N\geq 4.

Indeed, in the first case the function cc from (2.5) satisfies the estimate c(r)r(1N)(q2)/2|v(r)|q2c(r)\geq r^{(1-N)(q-2)/2}|v(r)|^{q-2} so that Atkinson’s oscillation criterion applies, see the first line and third column of the table on p.153 in [Swa79]. In the second case Noussair’s oscillation criterion result can be used in order finish step 1, see the third line and third column of the table on p.153 in [Swa79].

Remark 2.9.

If zg(z)/zz\mapsto g(z)/z is decreasing, then one can show that the first zero of uαu_{\alpha} is smaller than the first zero of uα~u_{\tilde{\alpha}} whenever 0<α<α~<α00<\alpha<\tilde{\alpha}<\alpha_{0}. Indeed, we set u:=uα,v:=uα~u:=u_{\alpha},v:=u_{\tilde{\alpha}}. Then the interval

I:={t>0:u(s)>v(s)>0 for all s(0,t)}I:=\{t>0:u(s)>v(s)>0\text{ for all }s\in(0,t)\}

is open, connected and nonempty and thus I=(0,r)I=(0,r^{*}) for some r>0r^{*}>0. On its right boundary we either have u(r)=v(r)0u(r^{*})=v(r^{*})\geq 0 or u(r)>v(r)=0u(r^{*})>v(r^{*})=0; so it remains to exclude the first possibility. Using u>v>0u>v>0 on II and (2.1) we have

(2.16) (rN1(uvvu))=rN1uv(g(v)vg(u)u)>0on I.\big{(}r^{N-1}(u^{\prime}v-v^{\prime}u)\big{)}^{\prime}=r^{N-1}uv\Big{(}\frac{g(v)}{v}-\frac{g(u)}{u}\Big{)}>0\quad\text{on }I.

Integrating (2.16) from 0 to rr^{*} the assumption u(r)=v(r)>0u(r^{*})=v(r^{*})>0 leads to

0<(uvvu)(r)=u(r)(uv)(r),hence(uv)(r)>0.0<(u^{\prime}v-v^{\prime}u)(r^{*})=u(r^{*})(u-v)^{\prime}(r^{*}),\quad\text{hence}\quad(u-v)^{\prime}(r^{*})>0.

On the other hand uv>0u-v>0 on I=(0,r)I=(0,r^{*}) and (uv)(r)=0(u-v)(r^{*})=0 implies (uv)(r)0(u-v)^{\prime}(r^{*})\leq 0, a contradiction. Thus u(r)>v(r)=0u(r^{*})>v(r^{*})=0 so that the first zero of vv comes before the first zero of uu.

2.2. The nonautonomous case

In this section we generalize Theorem 1.2 to a nonautonomous setting. Our aim is to identify mild assumptions on a nonautonomous nonlinearity gg that ensure the existence of a continuum of oscillating localized solutions of the initial value problems

(2.17) u′′N1ru=g(r,u),u(0)=α,u(0)=0-u^{\prime\prime}-\dfrac{N-1}{r}u=g(r,u),\qquad u(0)=\alpha,\;u^{\prime}(0)=0

that behave like r(1N)/2r^{(1-N)/2} at infinity in the sense of (2.4). Before formulating such assumptions and stating the corresponding existence result let us mention that our result applies to the nonlinearities (1.9),(1.10) under suitable conditions on the coefficient functions. This will be seen in Corollary 2.11 and Corollary 2.13 at the end of this section. Our existence results for (2.17) will be proven assuming that

(2.18) gC([0,+)×,)g\in C([0,+\infty)\times\mathbb{R},\mathbb{R}) is continuously differentiable w.r.t. rr.

Moreover, we suppose that there exist positive numbers α,α,λ,Λ\alpha_{*},\alpha^{*},\lambda,\Lambda and a locally Lipschitz continuous function g:g_{\infty}:\mathbb{R}\to\mathbb{R} such that

(2.19) limrg(r,)\displaystyle\lim_{r\to\infty}g(r,\cdot) =g()\displaystyle=g_{\infty}(\cdot) uniformly on [α,α][-\alpha_{*},\alpha^{*}]
(2.20) gr(r,z)z\displaystyle g_{r}(r,z)z 0\displaystyle\leq 0 on [0,+)×[α,α],\displaystyle\text{on }[0,+\infty)\times[-\alpha_{*},\alpha^{*}],
(2.21) λz2g(z)z\displaystyle\lambda z^{2}\leq g_{\infty}(z)z g(r,z)zΛz2\displaystyle\leq g(r,z)z\leq\Lambda z^{2} on [0,+)×[α,α].\displaystyle\text{on }[0,+\infty)\times[-\alpha_{*},\alpha^{*}].

These assumptions will allow us to prove the mere existence of an oscillating localized solution. In order to show the desired asymptotic behaviour we need some extra condition ”at infinity” where rr is large and the solution itself is small: We will assume that there exist ε,σ,C>0\varepsilon,\sigma,C>0 and some integrable function kk such that

(2.22) |2G(r,z)zg(r,z)|\displaystyle|2G(r,z)-zg(r,z)| Cz2|ln(z)|1σ,|z|ε,rε1\displaystyle\leq Cz^{2}|\ln(z)|^{-1-\sigma},\qquad|z|\leq\varepsilon,\,r\geq\varepsilon^{-1}
(2.23) gr(r,z)z\displaystyle g_{r}(r,z)z k(r)z2,|z|ε,rε1.\displaystyle\geq-k(r)z^{2},\quad\qquad\qquad|z|\leq\varepsilon,r\geq\varepsilon^{-1}.

These assumptions are rather technical but can be verified easily in concrete situations as we show in the proof of Corollary 2.11. Let us remark that our assumptions (1.2),(1.3),(1.4)(1.5) from the autonomous case (for any choice α=α(|α|,α0)\alpha^{*}=\alpha_{*}\in(|\alpha|,\alpha_{0})) are more restrictive than the assumptions used above. In particular, the following theorem generalizes our autonomous result.

Theorem 2.10.

Let N2N\geq 2. Moreover assume (2.19),(2.20),(2.21) as well as

(2.24) G(0,α)min{G(α),G(α)}for α[α,α].G(0,\alpha)\leq\min\{G_{\infty}(-\alpha_{*}),G_{\infty}(\alpha^{*})\}\quad\text{for }\alpha\in[-\alpha_{*},\alpha^{*}].

Then there is an oscillating, localized solution uu of (2.17) that satisfies u(0)=αu(0)=\alpha as well as

(2.25) uL()max{α,α}anduL()2G(0,α).\|u\|_{L^{\infty}(\mathbb{R})}\leq\max\{-\alpha_{*},\alpha^{*}\}\quad\text{and}\quad\|u^{\prime}\|_{L^{\infty}(\mathbb{R})}\leq\sqrt{2G(0,\alpha)}.

Moreover, if (2.22) and (2.23) hold, then we can find c,C>0c,C>0 such that

(2.26) cr(1N)/2|u(r)|+|u(r)|+|u′′(r)|Cr(1N)/2for r1.cr^{(1-N)/2}\leq|u(r)|+|u^{\prime}(r)|+|u^{\prime\prime}(r)|\leq Cr^{(1-N)/2}\quad\text{for }r\geq 1.
Proof.

The proof of our result follows the same argument of the proof of Theorem 1.2, so we only mention the main differences. For simplicity we only treat the case α>0\alpha>0 with (2.24). The existence of a maximally extended solution of (2.17) follows from a Peano type existence theorem for singular initial value problems, see Theorem 1 in [RW97].

Step 1 is proven as in the autonomous case where the function cc from (2.5) has to be replaced by c(r)=g(r,u(r))/u(r)(N1)(N3)/4r2c(r)=g(r,u(r))/u(r)-(N-1)(N-3)/4r^{2}. Assumption (2.21) ensures that cc is bounded from below by a positive constant as long as 0u(r)<α0\leq u(r)<\alpha so that uu has to attain a first zero. In step 2 one shows that Z(r):=u(r)2+2G(r,u(r))Z(r):=u^{\prime}(r)^{2}+2G(r,u(r)) is nondecreasing due to Gr(r,u(r))0G_{r}(r,u(r))\leq 0 for αu(r)α-\alpha_{*}\leq u(r)\leq\alpha^{*}, see (2.20). Arguing as in the autonomous case we find that uu decreases until it attains a local minimum at some r2>r1r_{2}>r_{1} with α<u(r2)<0-\alpha_{*}<u(r_{2})<0. More precisely one finds a sequence (rj)(r_{j}) such that all r2jr_{2j} are critical points and all r2j+1r_{2j+1} are zeros of uu with the additional property (the counterpart to (2.6))

2G(r0,u(r0))>u(r1)2>2G(r2,u(r2))>u(r3)2>2G(r4,u(r4))>.2G(r_{0},u(r_{0}))>u^{\prime}(r_{1})^{2}>2G(r_{2},u(r_{2}))>u^{\prime}(r_{3})^{2}>2G(r_{4},u(r_{4}))>\ldots.

This and G(r0,u(r0))=G(0,α)G(r_{0},u(r_{0}))=G(0,\alpha) yields the LL^{\infty}-bounds for uu^{\prime} whereas the LL^{\infty}-bounds follow from α<u(r2j)<α-\alpha_{*}<u(r_{2j})<\alpha^{*} for all jj\in\mathbb{N}. Hence, (2.25) is proved so that step 2 is finished. Step 3 is the same as in the proof of Theorem 1.2. Since the reasoning of Lemma 4.2 in [GZ08] may be adapted to our nonautonomous (but asymptotically autonomous) problem we also find (2.12), i.e.

(2.27) |u(r)|,|u(r)|,|u′′(r)|Cεr1N2+ε(r1)|u(r)|,|u^{\prime}(r)|,|u^{\prime\prime}(r)|\leq C_{\varepsilon}r^{\frac{1-N}{2}+\varepsilon}\qquad(r\geq 1)

In step 4 we use (2.22),(2.23) in order to study the asymptotics of the function

ψ(r):=v(r)2+2rN1G(r,u(r))\psi(r):=v^{\prime}(r)^{2}+2r^{N-1}G(r,u(r))

where v(r):=r(N1)/2u(r)v(r):=r^{(N-1)/2}u(r). One shows

ψ(r)\displaystyle\psi^{\prime}(r) =2rN1G(r,u(r))(N1r2G(r,u(r))u(r)g(r,u(r))2G(r,u(r))+Gr(r,u(r))G(r,u(r)))\displaystyle=2r^{N-1}G(r,u(r))\Big{(}\frac{N-1}{r}\frac{2G(r,u(r))-u(r)g(r,u(r))}{2G(r,u(r))}+\frac{G_{r}(r,u(r))}{G(r,u(r))}\Big{)}
+(N1)(N3)2r2v(r)v(r).\displaystyle\quad+\frac{(N-1)(N-3)}{2r^{2}}v(r)v^{\prime}(r).

For sufficiently large rrr\geq r^{*} we get estimates |Gr(r,u(r))|k(r)G(r,u(r))|G_{r}(r,u(r))|\leq k(r)G(r,u(r)). Moreover, with an analogous inequality as in (2.14) as well as (2.22),(2.27) we get

(2.28) 2G(r,u(r))u(r)g(r,u(r))2G(r,u(r))2G(r,u(r))u(r)g(r,u(r))λu(r)2C|ln(u(r))|1σCln(r)1σ\displaystyle\begin{aligned} \frac{2G(r,u(r))-u(r)g(r,u(r))}{2G(r,u(r))}&\leq\frac{2G(r,u(r))-u(r)g(r,u(r))}{\lambda u(r)^{2}}\\ &\leq C|\ln(u(r))|^{-1-\sigma}\\ &\leq C^{\prime}\ln(r)^{-1-\sigma}\end{aligned}

so that we may find as in the autonomous case a positive integrable function aa such that |ψ(r)|a(r)ψ(r)|\psi^{\prime}(r)|\leq a(r)\psi(r). This shows that ψ\psi is bounded from below and from above by a positive number. From this and

(2.29) λz22G(z)2G(r,z)Λz2on [0,+)×[α,α],\lambda z^{2}\leq 2G_{\infty}(z)\leq 2G(r,z)\leq\Lambda z^{2}\quad\text{on }[0,+\infty)\times[-\alpha_{*},\alpha^{*}],

which is a consequence of (2.21), we obtain the lower and upper bounds (2.26) and the proof is finished. \Box

Finally let us apply Theorem 2.10 to the special nonlinearities g1,g2g_{1},g_{2} given in (1.9),(1.10). We obtain the following results.

Corollary 2.11.

Let N2,p>2N\geq 2,p>2 and suppose that k,QC1([0,+),)k,Q\in C^{1}([0,+\infty),\mathbb{R}) are nonincreasing functions with limits k>0k_{\infty}>0 and QQ_{\infty}\in\mathbb{R}, respectively. Then there is a nonempty open interval II containing 0 and a continuum 𝒞={uαC2(N):αI}C2(N)\mathcal{C}=\{u_{\alpha}\in C^{2}(\mathbb{R}^{N}):\alpha\in I\}\subseteq\,C^{2}(\mathbb{R}^{N}) consisting of radially symmetric oscillating classical solutions of the equation

Δuk(|x|)2u=Q(|x|)|u|p2uin N-\Delta u-k(|x|)^{2}u=Q(|x|)|u|^{p-2}u\quad\text{in }\mathbb{R}^{N}

having the properties (2.25),(2.26) stated in Theorem 2.10. In case Q0Q_{\infty}\geq 0 we have I=I=\mathbb{R}.

Proof.

We set

g(r,z)=k(r)2z+Q(r)|z|p2z,g(z)=k2z+Q|z|p2zg(r,z)=k(r)^{2}z+Q(r)|z|^{p-2}z,\quad g_{\infty}(z)=k_{\infty}^{2}z+Q_{\infty}|z|^{p-2}z

and

G(r,z)=k(r)22z2+Q(r)p|z|p,G(z)=k22z2+Qp|z|p.G(r,z)=\frac{k(r)^{2}}{2}z^{2}+\frac{Q(r)}{p}|z|^{p},\quad G_{\infty}(z)=\frac{k_{\infty}^{2}}{2}z^{2}+\frac{Q_{\infty}}{p}|z|^{p}.

By the regularity assumptions on k,Qk,Q we have (2.18). Moreover, k,Q0k^{\prime},Q^{\prime}\leq 0 implies gr(r,z)z0g_{r}(r,z)z\leq 0 for all zz\in\mathbb{R} and thus (2.20). We set

(2.30) I:={α:G(0,α)<supG}={α:G(0,α)<(121p)k2(k2(Q))2p2}.\displaystyle\begin{aligned} I&:=\{\alpha\in\mathbb{R}:G(0,\alpha)<\sup_{\mathbb{R}}G_{\infty}\}\\ &=\Big{\{}\alpha\in\mathbb{R}:G(0,\alpha)<\big{(}\frac{1}{2}-\frac{1}{p}\big{)}k_{\infty}^{2}\Big{(}\frac{{k_{\infty}^{2}}_{~}}{(Q_{\infty})^{-}}\Big{)}^{\frac{2}{p-2}}\Big{\}}.\end{aligned}

Here, (Q)=max{Q,0}(Q_{\infty})^{-}=\max\{-Q_{\infty},0\}. For any given αI\alpha\in I we can choose

0<α=α<(|V|(Q))1p2s.t.G(α)=G(α)=G(0,α).0<\alpha_{*}=\alpha^{*}<\Big{(}\frac{|V_{\infty}|_{~}}{(Q_{\infty})^{-}}\Big{)}^{\frac{1}{p-2}}\;\text{s.t.}\;G_{\infty}(-\alpha_{*})=G_{\infty}(\alpha^{*})=G(0,\alpha).

For this choice of α,α\alpha_{*},\alpha^{*} assumption (2.21) holds. Finally, (2.22) follows from 2G(r,z)zg(r,z)=O(|z|p)2G(r,z)-zg(r,z)=O(|z|^{p}) as z0z\to 0 uniformly with respect to rr and (2.23) holds because gr(r,z)z(2k(r)k(r)+Q(r))z2g_{r}(r,z)z\geq(2k(r)k^{\prime}(r)+Q^{\prime}(r))z^{2} for |z|1|z|\leq 1. Hence, all assumptions of Theorem 2.10 are satisfied and the existence of solutions of (2.17) follows. Due to the unique solvability of these initial value problems and the Theorem of Ascoli-Arzelà they form a continuum in C2()C^{2}(\mathbb{R}) with respect to the C2C^{2}-convergence on compact sets. Finally we remark that (2.30) implies I=I=\mathbb{R} whenever Q0Q_{\infty}\geq 0. \Box

Remark 2.12.

Theorem 2.10 extends Theorem 4 in [EW14] in various directions. First of all, it provides more qualitative information of the solutions such as the W1,W^{1,\infty}-bounds, the oscillating behaviour of the solutions and the lower bounds for their decay at infinity. Additionally, we do not assume any global positivity assumption on ff. Furthermore, our assumption (2.22) is not covered by the hypotheses in [EW14].

The following result can be proved similarly and we state it for completeness.

Corollary 2.13.

Let N2N\geq 2 and suppose that λ,sC1([0,+),)\lambda,s\in C^{1}([0,+\infty),\mathbb{R}) are nondecreasing functions with limits λ,s\lambda_{\infty},s_{\infty}, respectively, such that ss is positive and λ<1/s\lambda_{\infty}<1/s_{\infty}. Then there is a nonempty open interval II containing 0 and a continuum 𝒞={uαC2(N):αI}\mathcal{C}=\{u_{\alpha}\in C^{2}(\mathbb{R}^{N}):\alpha\in I\} in C2(N)C^{2}(\mathbb{R}^{N}) consisting of radially symmetric oscillating classical solutions of the equation

Δu+λ(|x|)u=us(|x|)+u2in N-\Delta u+\lambda(|x|)u=\frac{u}{s(|x|)+u^{2}}\quad\text{in }\mathbb{R}^{N}

having the properties (2.25),(2.26) from Theorem 2.10. In the case λ0\lambda_{\infty}\leq 0 we have I=I=\mathbb{R}.

3. Nonradial solutions

In this section we study equation (1.11) proving Theorem 1.3 and Theorem 1.4. We will follow the argument introduced in [EW15] adapting their methods to our context. First note that, up to rescaling, we may assume k=1k=1 in the following. Let us introduce some notations in order to facilitate the reading. Let Ψ\Psi be the real part of the fundamental solution of the Helmholtz equation Δ1-\Delta-1 on N\mathbb{R}^{N} (see for example (11) in [EW15] ). Performing the transformation v=|Q|1/p|u|p2uv=|Q|^{1/p^{\prime}}|u|^{p-2}u for 1p+1p=1\frac{1}{p}+\frac{1}{p^{\prime}}=1 our problem amounts to solving

(3.1) |v|p2v=|Q|1/p[Ψ(|Q|1/pv)]in N.|v|^{p^{\prime}-2}v=-|Q|^{1/p}[\Psi\ast(|Q|^{1/p}v)]\quad\text{in }\mathbb{R}^{N}.

Notice that the right hand side comes with a negative sign in contrast to [EW15]. This is because we assume QQ to be negative so that Q=|Q|Q=-|Q|. Let us introduce the linear operators ,𝕂p:Lp(N)Lp(N)\mathbb{R},\,\mathbb{K}_{p}:L^{p^{\prime}}(\mathbb{R}^{N})\to L^{p}(\mathbb{R}^{N}) defined by

(3.2) (v)=Ψv,𝕂p(v)=|Q|1/p(|Q|1/pv).\mathbb{R}(v)=\Psi\ast v,\qquad\mathbb{K}_{p}(v)=|Q|^{1/p}\mathbb{R}(|Q|^{1/p}v).

Both \mathbb{R} and 𝕂p\mathbb{K}_{p} are continuous and we have for all f,gLp(N)f,g\in L^{p^{\prime}}(\mathbb{R}^{N})

(3.3) Nf(g)=Nf(Ψg)=limε0N(|ξ|21)f^(ξ)g^(ξ)(|ξ|21)2+ε2𝑑ξ,\int_{\mathbb{R}^{N}}f\mathbb{R}(g)=\int_{\mathbb{R}^{N}}f(\Psi\ast g)=\lim_{\varepsilon\to 0}\int_{\mathbb{R}^{N}}\frac{(|\xi|^{2}-1)\hat{f}(\xi)\hat{g}(\xi)}{(|\xi|^{2}-1)^{2}+\varepsilon^{2}}\,d\xi,

where f^,g^\hat{f},\hat{g} are the Fourier transforms of ff and gg, respectively. In view of the variational structure of (3.1) we define the functionals J,J¯:Lp(N)J,\,\bar{J}:L^{p^{\prime}}(\mathbb{R}^{N})\to\mathbb{R} via the formulas

(3.4) J(v):=1pN|v|p12Nv𝕂p(v)J¯(v):=1pN|v|p+12Nv𝕂p(v)\displaystyle\begin{aligned} J(v)&:=\frac{1}{p^{\prime}}\int_{\mathbb{R}^{N}}|v|^{p^{\prime}}-\frac{1}{2}\int_{\mathbb{R}^{N}}v\mathbb{K}_{p}(v)\\ \bar{J}(v)&:=\frac{1}{p^{\prime}}\int_{\mathbb{R}^{N}}|v|^{p^{\prime}}+\frac{1}{2}\int_{\mathbb{R}^{N}}v\mathbb{K}_{p}(v)\end{aligned}

so that the solutions of (3.1) are precisely the critical points of J¯\bar{J}, see (49) in [EW15]. Notice that the functional JJ is used when QQ is positive. Our main observation is that not only JJ but also J¯\bar{J} has the mountain pass geometry. This follows from the following Lemmas which are the counterparts of Lemma 4.2 and Lemma 5.2 in [EW15]. In the following we will denote with q\|\cdot\|_{q} the standard norm in the Lebesgue space Lq(N)L^{q}(\mathbb{R}^{N}).

Lemma 3.1.

Under the assumptions of Theorem 1.3 there is a function v0Lp(N)v_{0}\in L^{p^{\prime}}(\mathbb{R}^{N}) such that v0p>1,J¯(v0)<0\|v_{0}\|_{p^{\prime}}>1,\bar{J}(v_{0})<0.

Proof.

As in Lemma 4.2(ii) [EW15] it suffices to prove

(3.5) Nz𝕂pz<0\int_{\mathbb{R}^{N}}z\mathbb{K}_{p}z<0

for some zLp(N)z\in L^{p^{\prime}}(\mathbb{R}^{N}) because then one may take v0:=tzv_{0}:=tz for sufficiently large |t||t|. To this end let y𝒮(N)y\in\mathcal{S}(\mathbb{R}^{N}) be a nontrivial Schwartz function satisfying supp(y^)B1(0)\operatorname{supp}(\hat{y})\subset B_{1}(0). For δ>0\delta>0 we set

zδ:=y|Q|1/p1{|Q|>δ},μ:=supp(y^)|y^(ξ)|2|ξ|21𝑑ξ<0,z_{\delta}:=y|Q|^{-1/p}1_{\{|Q|>\delta\}},\qquad\mu:=\int_{\operatorname{supp}(\hat{y})}\frac{|\hat{y}(\xi)|^{2}}{|\xi|^{2}-1}\,d\xi<0,

where 1{|Q|>δ}1_{\{|Q|>\delta\}} is the indicator function of the set {xN:|Q(x)|>δ}\{x\in\mathbb{R}^{N}\,:\,|Q(x)|>\delta\}. Then we have zδLp(N)z_{\delta}\in L^{p^{\prime}}(\mathbb{R}^{N}) and thus 𝕂pzδLp(N)\mathbb{K}_{p}z_{\delta}\in L^{p}(\mathbb{R}^{N}). Hence, by definition of 𝕂p\mathbb{K}_{p}, the function yδ:=|Q|1/pzδ=y1{|Q|>δ}y_{\delta}:=|Q|^{1/p}z_{\delta}=y\cdot 1_{\{|Q|>\delta\}} satisfies

Nzδ(𝕂pzδ)=Nzδ|Q|1/p(|Q|1/pzδ)=Nyδ(yδ).\displaystyle\int_{\mathbb{R}^{N}}z_{\delta}(\mathbb{K}_{p}z_{\delta})=\int_{\mathbb{R}^{N}}z_{\delta}|Q|^{1/p}\mathbb{R}(|Q|^{1/p}z_{\delta})=\int_{\mathbb{R}^{N}}y_{\delta}\mathbb{R}(y_{\delta}).

Since we have |Q|>0|Q|>0 almost everywhere, we get yδyy_{\delta}\to y in Lp(N)L^{p^{\prime}}(\mathbb{R}^{N}) as δ0+\delta\to 0^{+}. Thus the continuity of \mathbb{R} implies that we can choose δ>0\delta>0 so small that the following holds:

Nzδ(𝕂pzδ)<Ny(y)+|μ|2.\displaystyle\int_{\mathbb{R}^{N}}z_{\delta}(\mathbb{K}_{p}z_{\delta})<\int_{\mathbb{R}^{N}}y\mathbb{R}(y)+\frac{|\mu|}{2}.

From this and (3.3) we infer

Nzδ(𝕂pzδ)<limε0N(|ξ|21)|y^(ξ)|2(|ξ|21)2+ε2𝑑ξ+|μ|2=μ+|μ|2<0\displaystyle\int_{\mathbb{R}^{N}}z_{\delta}(\mathbb{K}_{p}z_{\delta})<\lim_{\varepsilon\to 0}\int_{\mathbb{R}^{N}}\frac{(|\xi|^{2}-1)|\hat{y}(\xi)|^{2}}{(|\xi|^{2}-1)^{2}+\varepsilon^{2}}\,d\xi+\frac{|\mu|}{2}=\mu+\frac{|\mu|}{2}<0

which is all we had to show. \Box

Lemma 3.2.

Let the assumptions of Theorem 1.4 hold. Then for every mm\in\mathbb{N} there is an mm-dimensional subspace 𝒲Lp(N)\mathcal{W}\subset L^{p^{\prime}}(\mathbb{R}^{N}) with the following properties:

  • (i)

    Nv𝕂pv<0\displaystyle\int_{\mathbb{R}^{N}}v\mathbb{K}_{p}v<0 for all v𝒲{0}v\in\mathcal{W}\setminus\{0\}.

  • (ii)

    There exists R=R(𝒲)>0R=R(\mathcal{W})>0 such that J¯(v)0\bar{J}(v)\leq 0 for every v𝒲v\in\mathcal{W} with vpR\|v\|_{p^{\prime}}\geq R.

Proof.

Let y1,,ym𝒮(N)y^{1},\ldots,y^{m}\in\mathcal{S}(\mathbb{R}^{N}) be nontrivial Schwartz functions such that

(3.6) j=1msupp(y^j)B1(0),supp(y^j)supp(y^i)=(ij).\bigcup_{j=1}^{m}\operatorname{supp}(\hat{y}^{j})\subset B_{1}(0),\qquad\operatorname{supp}(\hat{y}^{j})\cap\operatorname{supp}(\hat{y}^{i})=\emptyset\quad(i\neq j).

For sufficiently small δ>0\delta>0 we then define

𝒲:=span{zδ1,,zδm}where zδj:=yj|Q|1/p1{|Q|>δ}.\mathcal{W}:=\operatorname{span}\{z_{\delta}^{1},\ldots,z_{\delta}^{m}\}\qquad\text{where }z_{\delta}^{j}:=y^{j}|Q|^{-1/p}1_{\{|Q|>\delta\}}.

Then (3.6) implies that 𝒲\mathcal{W} is mm-dimensional and similar calculations as above show (i) and (ii). \Box

With the aid of the above Lemmas the proofs of our theorems are essentially the same as in [EW14]. We indicate the main steps for the convenience of the reader.

Proof of Theorem 1.3:  Under the given assumptions J¯\bar{J} has the Mountain Pass geometry. Indeed, as in the parts (i),(iii) of Lemma 4.2 in [EW15] one proves that 0 is a strict local minimum and the boundedness of Palais-Smale sequences of J¯\bar{J}. In Lemma 3.1 we proved that there is a v0Lp(N)v_{0}\in L^{p^{\prime}}(\mathbb{R}^{N}) such that v0p>1,J¯(v0)<0\|v_{0}\|_{p^{\prime}}>1,\bar{J}(v_{0})<0. Hence, as in Lemma 6.1 [EW15] the Deformation Lemma implies the existence of a bounded Palais-Smale sequence (vm)(v_{m}) for J¯\bar{J} at its Mountain-Pass level c¯>0\bar{c}>0. Similar to the proof of Theorem 6.2 in [EW15] one has

limmN|Q|1/p(|Q|1/pvm)\displaystyle\lim_{m\to\infty}\int_{\mathbb{R}^{N}}|Q|^{1/p}\mathbb{R}(|Q|^{1/p}v_{m}) =2p2plimm[J¯(vm)+1pJ¯(vm)[vm]]\displaystyle=\frac{2p^{\prime}}{2-p^{\prime}}\lim_{m\to\infty}\Big{[}-\bar{J}(v_{m})+\frac{1}{p^{\prime}}\bar{J}^{\prime}(v_{m})[v_{m}]\Big{]}
=2p2pc¯<0.\displaystyle=-\frac{2p^{\prime}}{2-p^{\prime}}\bar{c}<0.

Then, Theorem 3.1 in [EW15] implies that there are R,ζ>0R,\zeta>0 and points xmNx_{m}\in\mathbb{R}^{N} and a subsequence, still denoted with (vm)(v_{m}), such that

BR(xm)|vm|pζ>0.\int_{B_{R}(x_{m})}|v_{m}|^{p^{\prime}}\geq\zeta>0.

From this point on the reasoning is the same as in [EW15] and we obtain that (vm)(v_{m}) converges weakly to a nontrivial solution vv of (1.11) and the solution uu of the original equation may be found via u=(|Q|1/pv)Lp(N)u=\mathbb{R}(|Q|^{1/p}v)\in L^{p}(\mathbb{R}^{N}). In particular uu satisfies (1.11) and

lim|x||xy|1|u(y)|p𝑑y=0\lim_{|x|\to\infty}\int_{|x-y|\leq 1}|u(y)|^{p}\,dy=0

so that replacing 2 by pp in the proof of Theorem C.3.1 in [Sim82] one proves u(x)0u(x)\to 0 as |x||x|\to\infty, namely that uu is localized. This implies

Δu+(k2+o(1))u=0,as |x|\Delta u+(k^{2}+o(1))u=0,\quad\text{as }|x|\to\infty

so that the PDE version of Sturm’s comparison principle (for instance Theorem 5.1 in [Swa68]) and the Strong Maximum Principle show that uu is oscillating. \Box

Proof of Theorem 1.4:  Lemma 3.2 yields all the required geometrical features of the symmetric Mountain Pass Theorem (see Theorem 6.5 in [Str08]). Moreover, Lemma 5.2 in [EW15] implies that the Palais-Smale condition holds for J¯\bar{J}, giving the existence of pairs of non-trivial localized solutions. The oscillation property follows again from Theorem 5.1 in [Swa68].

\Box

4. On the approximation by bounded domains

In this section we briefly address the question whether localized solutions of (1.1) can be approximated by solutions of the corresponding homogeneous Dirichlet problem on a large bounded domain. We are going to show that this method does not work in general. More precisely, we will prove that the positive minimizers of the Euler functionals associated with the problem on bounded domains diverge in H1(N)H^{1}(\mathbb{R}^{N}) as the domains approach N\mathbb{R}^{N}​. Even though the divergence will only be proved for the sequence of minimizers we believe that the analogous phenomenon occurs for broader classes of finite energy solutions, e.g. constrained minimizers, or solutions with a given upper bound on their nodal domains or on their Morse index. Throughout this section we will assume that the nonlinearity gg satisfies the hypotheses (1.2),(1.3),(1.4) as well as (1.5) with α0(0,+)\alpha_{0}\in(0,+\infty), in order to avoid some sub-critical growth conditions (see Remark 4.3).

Let ΩN\Omega\subset\mathbb{R}^{N} be a bounded domain, and consider the variational problem

(4.1) cΩ:=infH01(Ω)IΩwhere IΩ(u)=12Ω|u|2ΩG(u)c_{\Omega}:=\inf_{H_{0}^{1}(\Omega)}I_{\Omega}\quad\text{where }I_{\Omega}(u)=\frac{1}{2}\int_{\Omega}|\nabla u|^{2}-\int_{\Omega}G(u)

where G(z)G(z) denotes the primitive of gg such that G(0)=0G(0)=0. Notice that (1.2),(1.4),(1.5) implies G(z)C|z|2G(z)\leq C|z|^{2} for some C>0C>0 and for all zz\in\mathbb{R}, so that IΩ:H01(Ω){+}I_{\Omega}:H_{0}^{1}(\Omega)\mapsto\mathbb{R}\cup\{+\infty\} is well-defined. Bounded critical points of IΩI_{\Omega} are classical solutions of the boundary value problem

(4.2) {Δw=g(w)in Ω,wH01(Ω).\begin{cases}-\Delta w=g(w)&\text{in }\Omega,\\ w\in H_{0}^{1}(\Omega).&\end{cases}

In the following proposition we show that IΩI_{\Omega} admits a positive minimizer provided g(0)>λ1(Ω)g^{\prime}(0)>\lambda_{1}(\Omega) holds. More precisely, we have the following result.

Proposition 4.1.

Let ΩN\Omega\subset\mathbb{R}^{N} be a bounded domain and in addition to (1.2),(1.3),(1.4),(1.5) assume g(0)>λ1(Ω)g^{\prime}(0)>\lambda_{1}(\Omega). Then, there exists a global minimizer uΩu_{\Omega} of IΩI_{\Omega} in H01(Ω)H^{1}_{0}(\Omega) which is a solution of (4.2) satisfying 0<uΩ<α00<u_{\Omega}<\alpha_{0} in Ω\Omega.

Proof.

Hypotheses (1.3) and (1.5) imply that G(z)G(α0)G(z)\leq G(\alpha_{0}) holds for every zz\in\mathbb{R}. Hence, for all uH01(Ω)u\in H_{0}^{1}(\Omega) we have

IΩ(u)12Ω|u|2ΩG(α0)=12Ω|u|2|Ω|G(α0),I_{\Omega}(u)\geq\frac{1}{2}\int_{\Omega}|\nabla u|^{2}-\int_{\Omega}G(\alpha_{0})=\frac{1}{2}\int_{\Omega}|\nabla u|^{2}-|\Omega|G(\alpha_{0}),

which shows that IΩI_{\Omega} is coercive and bounded from below. Moreover, if ϕ1\phi_{1} denotes the eigenfunction associated to λ1(Ω)\lambda_{1}(\Omega), then

limt0IΩ(tϕ1)t2=12Ω|ϕ1|2G′′(0)ϕ12=12(λ1(Ω)g(0))Ωϕ12<0,\displaystyle\lim_{t\to 0}\frac{I_{\Omega}(t\phi_{1})}{t^{2}}=\frac{1}{2}\int_{\Omega}|\nabla\phi_{1}|^{2}-G^{\prime\prime}(0)\phi_{1}^{2}=\frac{1}{2}\Big{(}\lambda_{1}(\Omega)-g^{\prime}(0)\Big{)}\int_{\Omega}\phi_{1}^{2}<0,

so that cΩ<0=IΩ(0)c_{\Omega}<0=I_{\Omega}(0). Additionally, IΩI_{\Omega} is weakly sequentially lower semicontinuous so that there exists a minimizer uΩu_{\Omega}, which must be nontrivial because of cΩ<0c_{\Omega}<0. We may assume 0uΩα00\leq u_{\Omega}\leq\alpha_{0} because min{|uΩ|,α0}H01(Ω)\min\{|u_{\Omega}|,\alpha_{0}\}\in H_{0}^{1}(\Omega) is another minimizer of IΩI_{\Omega}. From the strong maximum principle we deduce that uΩu_{\Omega} satisfies 0<uΩ<α00<u_{\Omega}<\alpha_{0} in Ω\Omega as it is nontrivial. \Box

Remarks 4.2.
  • (a)

    If zg(z)/zz\mapsto g(z)/z is decreasing, then the condition g(0)>λ1(Ω)g^{\prime}(0)>\lambda_{1}(\Omega) is even necessary for the existence of a positive solution uH01(Ω)u\in H_{0}^{1}(\Omega). Indeed, testing (4.2) with uu gives

    λ1(Ω)Ωu2Ω|u|2=Ωg(u)u<g(0)Ωu2.\lambda_{1}(\Omega)\,\int_{\Omega}u^{2}\leq\int_{\Omega}|\nabla u|^{2}=\int_{\Omega}g(u)u<g^{\prime}(0)\,\int_{\Omega}u^{2}.

    In particular, note that our model nonlinearities g1,g2g_{1},g_{2} given in (1.6) and (1.7) satisfy this monotonicity property.

  • (b)

    If Ω\Omega is smooth then Theorem 1 in [BO86] shows that the the positive solution of (4.2) is unique provided zg(z)/zz\mapsto g(z)/z is decreasing.

Remark 4.3.

In this section we do not consider the case α0=+\alpha_{0}=+\infty in (1.5) because, without imposing additional growth conditions, the functional IΩI_{\Omega} may not be well-defined in this case and, even if it were, it need not be bounded from below.

Next we study the convergence of the minimizers obtained in Proposition 4.1. To this end, we consider a sequence (Ωn)(\Omega_{n}) of bounded domains satisfying ΩnΩn+1N\Omega_{n}\subset\Omega_{n+1}\subset\mathbb{R}^{N} and nΩn=N\bigcup_{n\in\mathbb{N}}\Omega_{n}=\mathbb{R}^{N}. Since every compact subset of N\mathbb{R}^{N} is covered by finitely many of those bounded domains we observe that λ1(Ωn)0\lambda_{1}(\Omega_{n})\to 0 as nn\to\infty so that, by the above proposition, the existence of positive minimizers is guaranteed for large nn provided that (1.4) holds. We show that the minimizers converge to the constant solution α0\alpha_{0} and therefore do not give any new finite energy solution.

Theorem 4.4.

Assume (1.2),(1.3),(1.4),(1.5) and let (Ωn)(\Omega_{n}) be a sequence of bounded domains such that ΩnΩn+1\Omega_{n}\subseteq\Omega_{n+1} and nΩn=N\cup_{n}\Omega_{n}=\mathbb{R}^{N}. Then, for all sufficiently large nn, there exists a nontrivial minimizer unu_{n} of IΩnI_{\Omega_{n}} on H01(Ωn)H^{1}_{0}(\Omega_{n}) having the following properties:

  • (a)

    0<un<α00<u_{n}<\alpha_{0} in Ωn\Omega_{n},

  • (b)

    IΩn(un)I_{\Omega_{n}}(u_{n})\to-\infty,

  • (c)

    unα0u_{n}\to\alpha_{0} in Cloc(N)C^{\infty}_{loc}(\mathbb{R}^{N}) and unLq(Ωn)\|u_{n}\|_{L^{q}(\Omega_{n})}\to\infty for all q[1,)q\in[1,\infty).

Proof.

Since λ1(Ωn)0\lambda_{1}(\Omega_{n})\to 0, taking into account (1.4) we find n0n_{0} such that, for every nn0n\geq n_{0}, λ1(Ωn)<g(0)\lambda_{1}(\Omega_{n})<g^{\prime}(0). As a consequence, we can apply Proposition 4.1 to deduce that there exists a sequence (un)nn0(u_{n})_{n\geq n_{0}} of positive minimizers of IΩnI_{\Omega_{n}} satisfying conclusion (a). In order to prove conclusion (b) let ϕC0(n)\phi\in C_{0}^{\infty}(\mathbb{R}^{n}) be given with ϕ2=ϕ=1\|\phi\|_{2}=\|\phi\|_{\infty}=1. For every kk\in\mathbb{N} we set

ϕk(x)=1kN/2ϕ(xk+ke1),\phi_{k}(x)=\frac{1}{k^{N/2}}\phi(\tfrac{x}{k}+ke_{1}),

so that ϕk2=1,ϕk1,ϕk20\|\phi_{k}\|_{2}=1,\|\phi_{k}\|_{\infty}\leq 1,\|\nabla\phi_{k}\|_{2}\to 0. Without loss of generality, we may assume σ(0,1)\sigma\in(0,1) from hypothesis (1.2) to be so small that 2+σ[2,2NN2]2+\sigma\in[2,\frac{2N}{N-2}] holds provided N>2N>2. Exploiting (1.2) and (1.5) we obtain positive constants A,CA,\,C such that

g(0)z22G(z)A|z|2+σfor |z|1,ϕk2+σC.g^{\prime}(0)z^{2}-2G(z)\leq A|z|^{2+\sigma}\qquad\text{for }|z|\leq 1,\qquad\|\phi_{k}\|_{2+\sigma}\leq C.

Then, for a fixed positive tmin{(g(0)/(4AC))1/σ,1}t\leq\min\{(g^{\prime}(0)/(4AC))^{1/\sigma},1\} and sufficiently large kk0k\geq k_{0} we have tϕk1\|t\phi_{k}\|_{\infty}\leq 1 so that the following estimate holds

2I(tϕk)=\displaystyle 2I(t\phi_{k})= t2N|ϕk|2g(0)ϕk2+N(g(0)(tϕk)22G(tϕk))\displaystyle t^{2}\int_{\mathbb{R}^{N}}|\nabla\phi_{k}|^{2}-g^{\prime}(0)\phi_{k}^{2}+\int_{\mathbb{R}^{N}}\big{(}g^{\prime}(0)(t\phi_{k})^{2}-2G(t\phi_{k})\big{)}
\displaystyle\leq g(0)2t2+AN|tϕk|2+σt22(g(0)+2tσAC)=:E,\displaystyle-\dfrac{g^{\prime}(0)}{2}t^{2}+A\int_{\mathbb{R}^{N}}|t\phi_{k}|^{2+\sigma}\leq\dfrac{t^{2}}{2}\left(-g^{\prime}(0)+2t^{\sigma}AC\right)=:-E,

where E>0E>0 by the choice of tt. Since the supports of (ϕk)(\phi_{k}) go off to infinity we find some k1k_{1}\in\mathbb{N} such that for all kk1k\geq k_{1} it results

2I(tϕk0+tϕk)E2+2I(tϕk0)+2I(tϕk)32E.\displaystyle 2I(t\phi_{k_{0}}+t\phi_{k})\leq\dfrac{E}{2}+2I(t\phi_{k_{0}})+2I(t\phi_{k})\leq-\frac{3}{2}E.

Inductively, we find k2<k3<k_{2}<k_{3}<\ldots such that for all kkmk\geq k_{m} we have

2I(tϕk0+tϕk1++tϕk)\displaystyle 2I(t\phi_{k_{0}}+t\phi_{k_{1}}+\ldots+t\phi_{k}) E2+2I(tϕk0+tϕk1++tϕkm1)+2I(tϕk)\displaystyle\leq\dfrac{E}{2}+2I(t\phi_{k_{0}}+t\phi_{k_{1}}+\ldots+t\phi_{k_{m-1}})+2I(t\phi_{k})
(1+m/2)E.\displaystyle\leq-(1+m/2)E.

Since for any given mm\in\mathbb{N} supp(tϕk0+tϕk1++tϕk)Ωn(t\phi_{k_{0}}+t\phi_{k_{1}}+\ldots+t\phi_{k})\subset\Omega_{n} for sufficiently large nn, the same estimate holds true for IΩnI_{\Omega_{n}}, yielding conclusion (b).

In order to show (c), note that the sequence (un)(u_{n}) is made of minimizers of IΩnI_{\Omega_{n}} so that

(4.3) Ωng(un)ϕ2Ωn|ϕ|2for all ϕCc1(Ωn).\int_{\Omega_{n}}g^{\prime}(u_{n})\phi^{2}\leq\int_{\Omega_{n}}|\nabla\phi|^{2}\qquad\text{for all }\phi\in C_{c}^{1}(\Omega_{n}).

By the Ascoli-Arzelà Theorem and interior Schauder estimates we find that (un)(u_{n}) converges in Cloc2(N)C^{2}_{loc}(\mathbb{R}^{N}) to some limit function uC2(N)u\in C^{2}(\mathbb{R}^{N}) satisfying 0u(x)α00\leq u(x)\leq\alpha_{0} for all xNx\in\mathbb{R}^{N} as well as (1.1). The Dominated Convergence Theorem allows to pass to the limit in (4.3) and we obtain

Ng(u)ϕ2N|ϕ|2for all ϕCc1(N),\int_{\mathbb{R}^{N}}g^{\prime}(u)\phi^{2}\leq\int_{\mathbb{R}^{N}}|\nabla\phi|^{2}\qquad\text{for all }\phi\in C_{c}^{1}(\mathbb{R}^{N}),

There are now three possibilities: either u0u\equiv 0 or uu is non-constant or uα0u\equiv\alpha_{0}. As a consequence, it is left to show that the first two possibilities do not occur.

First, assume by contraction that u0u\equiv 0, so that for any given compact KK we have unuu_{n}\to u uniformly on KK. Thanks to (1.4) we can find a sufficiently large nn such that g(un)/unδ2:=g(0)/2g(u_{n})/u_{n}\geq\delta^{2}:=g^{\prime}(0)/2. We may choose KK so large such that the fundamental solution ψ\psi of

Δψ+δ2ψ=0\Delta\psi+\delta^{2}\psi=0

changes sign within KK. Then the PDE version of Sturm’s comparison theorem (Theorem 5.1 in [Swa68]) shows that unu_{n} has a zero within KK contradicting the positivity of unu_{n}.
Assume now that uu is non-constant. Arguing as in the proof of Theorem 1.3 in [Far15] and applying Proposition 1.4 in [Far15] one shows that u>0u>0 and u<α0\|u\|_{\infty}<\alpha_{0}. As a consequence, the constant

c0:=min0sug(s)sc_{0}:=\min_{0\leq s\leq\|u\|_{\infty}}\frac{g(s)}{s}

turns out to be positive and, for any compact set KK we can find nn such that g(un)/unc0/2g(u_{n})/u_{n}\geq c_{0}/2. Choosing again KK sufficiently large we get a contradiction as above.

Hence, it turns out that uα0u\equiv\alpha_{0} and in particular we get unLq(Ωn)\|u_{n}\|_{L^{q}(\Omega_{n})}\to\infty for all q[1,)q\in[1,\infty). \Box

References

  • [ABC94] Antonio Ambrosetti, Haïm Brezis, and Giovanna Cerami. Combined effects of concave and convex nonlinearities in some elliptic problems. J. Funct. Anal., 122(2):519–543, 1994.
  • [AM07] Antonio Ambrosetti and Andrea Malchiodi. Nonlinear analysis and semilinear elliptic problems, volume 104 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 2007.
  • [BEP95] Lucio Boccardo, Miguel Escobedo, and Ireneo Peral. A Dirichlet problem involving critical exponents. Nonlinear Anal., 24(11):1639–1648, 1995.
  • [BGK83] Henri Berestycki, Thierry Gallouët, and Otared Kavian. équations de champs scalaires euclidiens non linéaires dans le plan. C. R. Acad. Sci. Paris Sér. I Math., 297(5):307–310, 1983.
  • [BL83] H. Berestycki and P.-L. Lions. Nonlinear scalar field equations. I. Existence of a ground state. Arch. Rational Mech. Anal., 82(4):313–345, 1983.
  • [BO86] Haïm Brezis and Luc Oswald. Remarks on sublinear elliptic equations. Nonlinear Anal., 10(1):55–64, 1986.
  • [CBWM04] Yi-Fan Chen, Kale Beckwitt, Frank W. Wise, and Boris A. Malomed. Criteria for the experimental observation of multidimensional optical solitons in saturable media. Phys. Rev. E, 70:046610, Oct 2004.
  • [CCMS97] Demetrios N. Christodoulides, Tamer H. Coskun, Matthew Mitchell, and Mordechai Segev. Theory of incoherent self-focusing in biased photorefractive media. Phys. Rev. Lett., 78:646–649, Jan 1997.
  • [Eve15] Gilles Evequoz. A dual approach in Orlicz spaces for the nonlinear Helmholtz equation. Z. Angew. Math. Phys., 66(6):2995–3015, 2015.
  • [EW] Gilles Evequoz and Tobias Weth. Branch continuation inside the essential spectrum for the nonlinear schr dinger equation. Preprint, https://arxiv.org/abs/1606.00606.
  • [EW14] Gilles Evequoz and Tobias Weth. Real solutions to the nonlinear Helmholtz equation with local nonlinearity. Arch. Ration. Mech. Anal., 211(2):359–388, 2014.
  • [EW15] Gilles Evequoz and Tobias Weth. Dual variational methods and nonvanishing for the nonlinear Helmholtz equation. Adv. Math., 280:690–728, 2015.
  • [Far15] Alberto Farina. Some symmetry results and Liouville-type theorems for solutions to semilinear equations. Nonlinear Anal., 121:223–229, 2015.
  • [GLZ09] Changfeng Gui, Xue Luo, and Feng Zhou. Asymptotic behavior of oscillating radial solutions to certain nonlinear equations, Part II. Methods Appl. Anal., 16(4):459–468, 2009.
  • [Gut04] Susana Gutiérrez. Non trivial LqL^{q} solutions to the Ginzburg-Landau equation. Math. Ann., 328(1-2):1–25, 2004.
  • [GZ08] Changfeng Gui and Feng Zhou. Asymptotic behavior of oscillating radial solutions to certain nonlinear equations. Methods Appl. Anal., 15(3):285–295, 2008.
  • [RW97] Wolfgang Reichel and Wolfgang Walter. Radial solutions of equations and inequalities involving the pp-Laplacian. J. Inequal. Appl., 1(1):47–71, 1997.
  • [Sim82] Barry Simon. Schrödinger semigroups. Bull. Amer. Math. Soc. (N.S.), 7(3):447–526, 1982.
  • [Str77] Walter A. Strauss. Existence of solitary waves in higher dimensions. Comm. Math. Phys., 55(2):149–162, 1977.
  • [Str08] Michael Struwe. Variational methods, volume 34 of Results in Mathematics and Related Areas. 3rd Series. A Series of Modern Surveys in Mathematics. Springer-Verlag, Berlin, fourth edition, 2008. Applications to nonlinear partial differential equations and Hamiltonian systems.
  • [SV14] Nicola Soave and Gianmaria Verzini. Bounded solutions for a forced bounded oscillator without friction. J. Differential Equations, 256(7):2526–2558, 2014.
  • [Swa68] C. A. Swanson. Comparison and oscillation theory of linear differential equations. Academic Press, New York-London, 1968. Mathematics in Science and Engineering, Vol. 48.
  • [Swa79] C. A. Swanson. Semilinear second-order elliptic oscillation. Canad. Math. Bull., 22(2):139–157, 1979.
  • [SZ99] C. A. Stuart and H. S. Zhou. Applying the mountain pass theorem to an asymptotically linear elliptic equation on 𝐑N{\bf R}^{N}. Comm. Partial Differential Equations, 24(9-10):1731–1758, 1999.
  • [Ver03] Gianmaria Verzini. Bounded solutions to superlinear ODEs: a variational approach. Nonlinearity, 16(6):2013–2028, 2003.
  • [Wil96] Michel Willem. Minimax theorems, volume 24 of Progress in Nonlinear Differential Equations and their Applications. Birkhäuser Boston, Inc., Boston, MA, 1996.
  • [Yan04] Jianke Yang. Stability of vortex solitons in a photorefractive optical lattice. New Journal of Physics, 6(1):47, 2004.