Home p-Laplacian Equations in ℝN with Finite Potential via the Truncation Method
Article Open Access

p-Laplacian Equations in ℝN with Finite Potential via the Truncation Method

  • Xiangqing Liu EMAIL logo and Junfang Zhao
Published/Copyright: November 8, 2016

Abstract

We consider the problem -Δpu+a(x)|u|p-2u=|u|q-2u in N, where 1<p<N, p<q<p*=NpN-p, Δp is the p-Laplacian operator, and the potential function a is positive, bounded and verifies suitable decay assumptions. The existence of infinitely many solutions of the equation is proved via the truncation method.

MSC 2010: 35B05; 35B45

1 Introduction

In this paper, we are concerned with the existence of multiple solutions of the p-Laplacian equation in N,

($P_{0}$) { - Δ p u + a ( x ) | u | p - 2 u = | u | q - 2 u in  N , u W 1 , p ( N ) ,

where 1<p<N, p<q<p*=NpN-p, Δp is the p-Laplacian operator, Δpu=(|u|p-2u), and the potential function a is positive, finite and verifies suitable decay assumptions.

Problem ($P_{0}$) has a variational structure, given by the functional I: W1,p(N) defined by

I ( u ) = 1 p N ( | u | p + a ( x ) | u | p ) 𝑑 x - 1 q N | u | q 𝑑 x .

The embedding W1,p(N)Lq(N) is continuous, but not compact due to the translations, therefore the Palais–Smale condition is not satisfied by I, that is, problem ($P_{0}$) lacks the compactness property. Since the pioneering work [3] of H. Brezis and L. Nirenberg, significant progress has been made in recent decades. In particular, Devillanova and Solimini [7] dealt with the critical growth problem

($P_{a}$) { - Δ u = | u | 2 * - 2 u + λ u in  Ω , u = 0 on  Ω ,

where Ω is an open regular domain of N (N3) and 2*=2NN-2 is the critical exponent, and proved a multiplicity result. The solutions are found as limits of solutions of approximated problems with subcritical growth. The lack of compactness due to the scaling does not allow to deduce that a sequence of approximated solutions must have a converging subsequence but the fact that they solve the approximated problems gives, with the use of a local Pohožaev identity, some extra estimates which lead to a proof of the desired compactness.

This argument was adapted in [5] to problem ($P_{0}$) with p=2. Cerami, Devillanova and Solimini [5] considered the following related problem in a sequence of balls in N:

($P_{b}$) { - Δ u + a ( x ) u = | u | q - 2 u in  B R , u = 0 on  B R ,

where BR is a ball in N, BR={xxN,|x|<R}. Problem ($P_{b}$) possesses infinitely many solutions, which are used as approximation solutions for the problem in the whole space.

Another way to deal with the problem in N is to introduce a coercive potential. Namely we first consider the equation

($P_{W}$) { - Δ u + ( ε W ( x ) + a ( x ) ) u = | u | q - 2 u in  N , u H W ,

where ε>0, W is a coercive function, say W(x)=1+|x|2, and HW is the Hilbert space

H W = { u | u H 1 ( N ) , N W ( x ) u 2 𝑑 x < + }

equipped with the norm

u H W = ( N ( | u | 2 + W ( x ) u 2 ) 𝑑 x ) 1 2 .

The embedding from HW to Lq(N) (2q<2*=2NN-2) is compact and hence problem ($P_{W}$) possesses infinitely many solutions. By taking the limit ε0, we obtain solutions of the original problem with finite potential. This approach was introduced in [8]. In fact, Liu and Wang [8] considered the existence of solutions of some quasilinear elliptic equations, a model problem is the so-called modified nonlinear Schrödinger equation

($P_{m}$) Δ u + 1 2 u Δ u 2 - V ( x ) u + | u | p - 2 u = 0 in N

with p(4,4NN-2).

For the p-Laplacian problems with critical nonlinearities, one may meet with different technical difficulties for different problems, and therefore there are different ways to deal with these kinds of problems. In [4], Cao, Peng and Yan considered the following p-Laplacian problem with critical Sobolev exponent:

{ - Δ p u = | u | p * - 2 u + μ | u | p - 2 u in  Ω , u = 0 on  Ω ,

where 1<p<N, p*=NpN-p, μ>0 and Ω is an open bounded domain in N. They first investigated the equations of subcritical growth and then proved a convergence theorem to obtain solutions of the original problem, following the argument in [7].

As for the existence of multiple solutions of the original problem we need to check that the multiple solutions of the approximated problems do not converge all to the same solution of the limit problem. In both [5] and [7], some estimates on the Morse index are employed, which has been used as one of possible devices to distinguish the limits of the multiple approximated solutions by their original variational characterization. For general p-Laplacian equations with p2 we have no information on the Morse index, therefore the approaches in [7, 5, 8] can not be extended in a straightforward way to problems involving Δp with p2.

In [4], Cao, Peng and Yan proved the existence of multiple solutions, with the help of a claim that if the functional corresponding to the original problem has only finitely many critical values, then the critical set has genus at least two. At this last step the argument seems to be not so easy to follow (see [4, (4.20)]). Actually this is one of the motivations for us to use the truncation method.

Note that every solution of problem ($P_{0}$) decays exponentially. Inspired by this fact, in this paper we are going to use the truncation technique to obtain infinitely many solutions of problem ($P_{0}$). We first consider truncated problems, solutions of which will be used as approximating solutions. By a compactness argument, similar to that in [7, 5], in particular with the use of the Pohožaev formula, theorems of convergence of solutions of approximating problems and existence of multiple solutions of the original problem ($P_{0}$) are proved. The advantage of the truncation method is that the original problems and the approximating problems share some common solutions, more and more as the parameter tends to zero.

Moreover, the truncation method mentioned in this paper may be used to deal with the p-Laplacian problem with critical Sobolev exponent. Although the problem involves critical Sobolev exponent, solutions of the problem are still bounded. So we can use the truncation technique to reduce the critical problem to a subcritical one, see Section 4 for some details.

Let us describe the truncation method in more details. Let φC0() be such that φ(t)=1 for |t|1, φ(t)=0 for |t|2, φ is decreasing in [1,2], and φ is even. For λ>0, (x,s)N×, define

b λ ( x , s ) = φ ( λ e λ 1 + | x | 2 s ) ,
m λ ( x , s ) = 0 s b λ ( x , τ ) 𝑑 τ ,
F λ ( x , s ) = 1 q | s | r | m λ ( x , s ) | q - r , p < r < q ,
f λ ( x , s ) = s F λ ( x , s ) ,

where r(p,q) is a fixed number, for example r=12(p+q). For λ=0 we understand that m0(x,s)s, F0(x,s)1q|s|q and f0(x,s)|s|q-2s. We consider the equation

($P_{\lambda}$) { - Δ p u + a ( x ) | u | p - 2 u = f λ ( x , u ) in  N , u W 1 , p ( N ) .

The embedding from W1,p(N) to Lq(N) (p<q<NpN-p) is not compact, but if unu in W1,p(N), then mλ(,un)mλ(,u) in Lq(N) for λ>0. As a problem with compactness, ($P_{\lambda}$), λ>0, possesses infinitely many solutions. The truncation method has its advantages, that is, if u is a solution of ($P_{\lambda}$), λ>0, and satisfies the estimate

(1.1) | u ( x ) | 1 λ e - λ 1 + | x | 2 for all  x N ,

then automatically u will be a solution of the original problem ($P_{0}$), too. We can prove that as λ0 more and more solutions u of problem ($P_{\lambda}$) satisfy the above estimate (1.1), hence we obtain infinitely many solutions of problem ($P_{0}$).

The truncation technique is not new. People use this technique to deal with various problems, because of various reasons. For example, in equation ($P_{a}$) with 0<λ<λ1, we replace the term |u|2*-2u by u+2*-1 in order to obtain a positive solution. In the quasilinear equation ($P_{m}$) we replace the quasilinear term uΔu2 by uMΔ(uuM), where uM=u if |u|M, uM=±M if ±uM. See also [9, 10] for Hamiltonian systems and for nonlinear wave equations. But as to the authors’ knowledge, the application of the truncation technique to a problem on the unbounded domains by introducing a decay function, such as 1λe-λ1+|x|2, in order to overcome the difficulty of lack of compactness seems to be new.

Now we make the following assumptions on the potential function (see [5]):

  1. a C 1 ( N , ) .

  2. There exist a1,a0>0 such that a0a(x)a1, xN.

  3. For all α>0, we have lim|x|ar(x)eα|x|=+, where ar(x)=x|x|a(x).

  4. There exists c¯>1 such that |a(x)|c¯ar(x) for all xN, |x|c¯.

Our main results are the following two theorems.

Theorem 1.1.

Assume (a1)(a4). Given M>0, there exists μ=μ(M) such that if uW1,p(RN) is a solution of problem ($P_{\lambda}$), λ0 and uM. Then

| u ( x ) | 1 μ e - μ 1 + | x | 2 for all  x N .

Theorem 1.2.

Assume (a1)(a4). Problem ($P_{0}$) has infinitely many solutions.

Throughout the paper, ||p denotes the norm in Lp(N), denotes the norm in W1,p(N), and denote the strong convergence and the weak convergence, respectively. Furthermore, c denotes a (possibly different) positive constant.

2 Uniform Bounds

Let unW1,p(N) be a solution of problem ($P_{\lambda}$) with λ=λn0, n=1,2,. Assume unM. By [12, Theorem 2.1 and Theorem 3.3] (see also [13]), {un} has a profile decomposition

(2.1) u n = u + k Λ U k ( - x n , k ) + r n ,

where uW1,p(N), UkW1,p(N), kΛ, rnW1,p(N) and {xn,k}N such that the following properties hold:

  1. u n u , un(+xn,k)Uk in W1,p(N) as n, kΛ.

  2. | x n , k | , |xn,k-xn,l| as n, k,lΛ, kl.

  3. | u | q q + k Λ | U k | q q lim ¯ n | u n | q q , p<q<p*.

  4. r n 0 in Lq(N) as n, p<q<p*.

Lemma 2.1.

For (x,s)RN×R, the following properties hold:

  1. | s | b λ ( x , s ) | m λ ( x , s ) | , smλ(x,s)0.

  2. | m λ ( x , s ) | min { | s | , 2 λ e - λ 1 + | x | 2 } and m λ ( x , s ) = s if | s | 1 λ e - λ 1 + | x | 2 .

  3. | f λ ( x , s ) | | s | r - 1 | m λ ( x , s ) | q - r .

  4. 1 r s f λ ( x , s ) - F λ ( x , s ) = q - r q r | s | r + 1 | m λ ( x , s ) | q - r - 1 b λ ( x , s ) .

  5. x m λ ( x , s ) = - λ x 1 + | x | 2 ( m λ ( x , s ) - s b λ ( x , s ) ) and

    x F λ ( x , s ) = - ( 1 - r q ) λ x 1 + | x | 2 | s | r | m λ ( x , s ) | q - r - 1 | m λ ( x , s ) - s b λ ( x , s ) | .

Proof.

(1) By the definition of mλ(x,s) and the intermediate value theorem, we have

m λ ( x , s ) = 0 s b λ ( x , τ ) 𝑑 τ = s b λ ( x , ξ ) , 0 < ξ < s .

Then by the definition of bλ(x,s) and the fact that φ is decreasing in [1,2], we get

| s | b λ ( x , s ) = | s | φ ( λ e λ 1 + | x | 2 s ) | s | φ ( λ e λ 1 + | x | 2 ξ ) = | m λ ( x , s ) | .

(2) We consider the case s>0. Since bλ(x,s)1, we have

m λ ( x , s ) = 0 s b λ ( x , τ ) 𝑑 τ s .

For τ>2λe-λ1+|x|2, we have λeλ1+|x|2τ2 and bλ(x,τ)=φ(λeλ1+|x|2τ)=0. Then

m λ ( x , s ) = 0 s b λ ( x , τ ) 𝑑 τ 0 2 λ e - λ 1 + | x | 2 b λ ( x , τ ) 𝑑 τ 2 λ e - λ 1 + | x | 2 .

Hence,

m λ ( x , s ) min { s , 2 λ e - λ 1 + | x | 2 } .

For 0<τs1λe-λ1+|x|2 we have

λ e λ 1 + | x | 2 τ 1 and b λ ( x , τ ) = φ ( λ e λ 1 + | x | 2 τ ) = 1 ,

hence for s1λe-λ1+|x|2 we obtain

m λ ( x , s ) = 0 s b λ ( x , τ ) 𝑑 τ = 0 s 𝑑 τ = s .

(3) Since

(2.2) f λ ( x , s ) = F λ ( x , s ) s = r q | s | r - 2 s | m λ ( x , s ) | q - r + q - r q | s | r | m λ ( x , s ) | q - r - 2 m λ ( x , s ) b λ ( x , s )

and |s|bλ(x,s)|mλ(x,s)|, we have

| f λ ( x , s ) | r q | s | r - 1 | m λ ( x , s ) | q - r + q - r q | s | r - 1 | m λ ( x , s ) | q - r = | s | r - 1 | m λ ( x , s ) | q - r .

(4) By (2.2) and smλ(x,s)0, we get

1 r s f λ ( x , s ) - F λ ( x , s ) = q - r q r | s | r s | m λ ( x , s ) | q - r - 2 m λ ( x , s ) b λ ( x , s )
= q - r q r | s | r + 1 | m λ ( x , s ) | q - r - 1 b λ ( x , s ) .

(5) By the definition of mλ(x,s), we obtain

x m λ ( x , s ) = 0 s x b λ ( x , τ ) 𝑑 τ
= 0 s φ ( λ e λ 1 + | x | 2 τ ) λ e λ 1 + | x | 2 τ λ x 1 + | x | 2 𝑑 τ
= 0 s λ x 1 + | x | 2 τ 𝑑 φ ( λ e λ 1 + | x | 2 τ )
= λ x 1 + | x | 2 ( s φ ( λ e λ 1 + | x | 2 s ) - 0 s φ ( λ e λ 1 + | x | 2 τ ) 𝑑 τ )
= - λ x 1 + | x | 2 ( m λ ( x , s ) - s b λ ( x , s ) ) .

Thus,

x F λ ( x , s ) = q - r q | s | r | m λ ( x , s ) | q - r - 2 m λ ( x , s ) x m λ ( x , s )
= - ( 1 - r q ) λ x 1 + | x | 2 | s | r | m λ ( x , s ) | q - r - 1 | m λ ( x , s ) - s b λ ( x , s ) | .

Lemma 2.2.

Let uW1,p(RN) be a solution of ($P_{\lambda}$), λ>0. Then v=|u| satisfies the following differential inequality:

(2.3) N | v | p - 2 v φ d x + a 0 N v p - 1 φ 𝑑 x N v q - 1 φ 𝑑 x , φ W 1 , p ( N ) , φ 0 .

Proof.

Set vε=u2+ε2-ε, ε>0. Then vεv in W1,p(N) as ε0. By Lemma 2.1 (2) and (3), we have

| f λ ( x , s ) | | s | r - 1 | m λ ( x , s ) | q - r | s | q - 1 .

Then for φC0(N), φ0, we have

N | v | p - 2 v ε φ d x = N | u | p - 2 u u ( u 2 + ε 2 ) 1 2 φ d x
= N | u | p - 2 u ( u ( u 2 + ε 2 ) 1 2 φ ) d x - N | u | p ε 2 ( u 2 + ε 2 ) 3 2 φ 𝑑 x
N | u | p - 2 u ( u ( u 2 + ε 2 ) 1 2 φ ) d x
= - N a ( x ) | u | p - 2 u u ( u 2 + ε 2 ) 1 2 φ 𝑑 x + N f λ ( x , u ) u ( u 2 + ε 2 ) 1 2 φ 𝑑 x
- a 0 N v p - 1 v ( v 2 + ε 2 ) 1 2 φ 𝑑 x + N v q - 1 v ( v 2 + ε 2 ) 1 2 φ 𝑑 x .

Let ε0. By Lebesgue’s dominated convergence theorem we obtain (2.3) for φC0(N), φ0. By approximation, (2.3) holds for φW1,p(N), φ0. ∎

Lemma 2.3.

Let unW1,p(RN) be a solution of ($P_{\lambda}$) with λ=λn0, n=1,2,, and {xn}RN. Suppose u~n=un(+xn)U in W1,p(RN). Then u~nU in Wloc1,p(RN).

Proof.

The sequence u~n satisfies the equation

N | u ~ n | p - 2 u ~ n φ d x + N a ( x + x n ) | u ~ n | p - 2 u ~ n φ 𝑑 x = N f λ n ( x + x n , u ~ n ) φ 𝑑 x

for all φW1,p(N). Let R>0, φC0(N), such that φ(x)=1 for |x|R and φ(x)=0 for |x|2R. The sequence u~n converges in Llocq(N), 1q<p*. We have

N ( | u ~ n | p - 2 u ~ n - | u ~ m | p - 2 u ~ m , u ~ n - u ~ m ) φ 𝑑 x
= - N ( | u ~ n | p - 2 u ~ n - | u ~ m | p - 2 u ~ m , φ ) ( u ~ n - u ~ m ) 𝑑 x
    - N ( a ( x + x n ) | u ~ n | p - 2 u ~ n - a ( x + x m ) | u ~ m | p - 2 u ~ m ) ( u ~ n - u ~ m ) φ 𝑑 x
    + N ( f λ n ( x + x n , u ~ n ) - f λ m ( x + x m , u ~ m ) ) ( u ~ n - u ~ m ) φ 𝑑 x
c B 2 R ( | u ~ n | p - 1 + | u ~ m | p - 1 ) | u ~ n - u ~ m | 𝑑 x + c B 2 R ( | u ~ n | p - 1 + | u ~ m | p - 1 ) | u ~ n - u ~ m | 𝑑 x
    + c B 2 R ( | u ~ n | q - 1 + | u ~ m | q - 1 ) | u ~ n - u ~ m | 𝑑 x
c ( N ( | u ~ n | p + | u ~ m | p + | u ~ n | p + | u ~ m | p ) 𝑑 x ) p - 1 p ( B 2 R | u ~ n - u ~ m | p 𝑑 x ) 1 p
    + ( N ( | u ~ n | q + | u ~ m | q ) 𝑑 x ) q - 1 q ( B 2 R | u ~ n - u ~ m | q 𝑑 x ) 1 q
(2.4) c ( B 2 R | u ~ n - u ~ m | p 𝑑 x ) 1 p + c ( B 2 R | u ~ n - u ~ m | q 𝑑 x ) 1 q 0 as  n , m .

The following elementary inequalities are very useful (see [6]). For p>1 there exists a constant cp such that for ξ,ηN we have

(2.5) ( | ξ | p - 2 ξ - | η | p - 2 η , ξ - η ) c p | ξ - η | p if  p 2 ,
(2.6) ( | ξ | p - 2 ξ - | η | p - 2 η , ξ - η ) c p | ξ - η | 2 ( | ξ | 2 - p + | η | 2 - p ) - 1 if  1 < p < 2 .

For p2, by (2.4) and (2.5) we have

B R | u ~ n - u ~ m | p 𝑑 x c N ( | u ~ n | p - 2 u ~ n - | u ~ m | p - 2 u ~ m , u ~ n - u ~ m ) φ 𝑑 x 0 as  n , m .

For 1<p2, by (2.4) and (2.6) we have

B R | u ~ n - u ~ m | p 𝑑 x c B R | ( | u ~ n | p - 2 u ~ n - | u ~ m | p - 2 u ~ m , u ~ n - u ~ m ) | p 2 ( | u ~ n | 2 - p + | u ~ m | 2 - p ) p 2 𝑑 x
c ( B R ( | u ~ n | p - 2 u ~ n - | u ~ m | p - 2 u ~ m , u ~ n - u ~ m ) 𝑑 x ) p 2 ( B R ( | u ~ n | p + | u ~ m | p ) 𝑑 x ) 2 - p 2
c ( N ( | u ~ n | p - 2 u ~ n - | u ~ m | p - 2 u ~ m , u ~ n - u ~ m ) φ 𝑑 x ) p 2 0 as  n , m ,

since φ=1 in BR and φ0. Hence {u~n} converges in Wloc1,p(N). ∎

Lemma 2.4.

Let the decomposition (2.1) for un hold. Then we have the following results:

  1. v = | u | and V k = | U k | ( k Λ ) satisfy the differential inequality

    (2.7) N | v | p - 2 v φ d x + a 0 N v p - 1 φ 𝑑 x N v q - 1 φ 𝑑 x for all  φ W 1 , p ( N ) , φ 0 .

  2. The index set Λ is finite.

  3. There exists α > 0 such that

    v ( x ) 1 α e - α 1 + | x | 2 , V k ( x ) 1 α e - α 1 + | x | 2 for all  x N , k Λ .

Proof.

(1) Denote vn=|un|. By Lemma 2.2, the sequence vn satisfies the differential inequality

(2.8) N | v n | p - 2 v n φ d x + a 0 N v n p - 1 φ 𝑑 x N v n q - 1 φ 𝑑 x for all  φ W 1 , p ( N ) , φ 0 .

We have unu in W1,p(N). By Lemma 2.3, we know that unu in Wloc1,p(N), consequently we obtain vn=|un||u|=v in Wloc1,p(N). Take the limit n in (2.8). For φC0(N), φ0, we obtain

N | v | p - 2 v φ d x + a 0 N v p - 1 φ 𝑑 x N v q - 1 φ 𝑑 x .

By approximation this inequality holds for φW1,p(N), φ0. Similarly, we can prove that Vk, for kΛ, satisfies (2.7).

(2) By (2.7) and the Sobolev embedding theorem, we have

( N V k q 𝑑 x ) p q S p , q - 1 N ( | V k | p + V k p ) 𝑑 x c 0 N V k q 𝑑 x ,

where Sp,q is the Sobolev constant of the embedding from W1,p(N) to Lq(N):

S p , q = inf u W 1 , p ( N ) { 0 } N ( | u | p + u p ) 𝑑 x ( N u q 𝑑 x ) p q .

Hence N|Uk|q𝑑x=NVkq𝑑xm for some m>0. By property (3) of the decomposition (2.1), Λ is finite.

(3) This fact is well known (see [5, 8]). For the convenience of the readers, we give the proof for the function v.

Step 1. We use Moser’s iteration. For a function v and any T>0 define vT as vT=v if v(x)T and vT=T if v(x)T. For 1ρ<R2 let ψC0(N) such that ψ(x)=1 for |x|ρ, ψ(x)=0 for |x|R and |ψ|2R-ρ for all xN. Choose T>0, r1 and take φ=vvTp(r-1)ψp as the test function in (2.7). We have

v T ( ( p ( r - 1 ) + 1 ) v p ( r - 1 ) ψ p | v | p + p v p ( r - 1 ) + 1 ψ p - 1 | v | p - 2 v ψ ) 𝑑 x
    + v T ( T p ( r - 1 ) ψ p | v | p + p T p ( r - 1 ) v ψ p - 1 | v | p - 2 v ψ ) 𝑑 x + a 0 N v p v T p ( r - 1 ) ψ p 𝑑 x
(2.9) N v q v T p ( r - 1 ) ψ p 𝑑 x .

By (2.9), we get

(2.10) N v T p ( r - 1 ) ψ p | v | p 𝑑 x p N v T p ( r - 1 ) v ψ p - 1 | v | p - 1 | ψ | 𝑑 x + N v q v T p ( r - 1 ) ψ p 𝑑 x .

By (2.10), the Hölder inequality and the Young inequality, we have

(2.11) N v T p ( r - 1 ) ψ p | v | p 𝑑 x c N ( v v T r - 1 ) p | ψ | p 𝑑 x + c N v q v T p ( r - 1 ) ψ p 𝑑 x .

Then by (2.11), we obtain

N | ( v v T r - 1 ψ ) | p 𝑑 x c N | ( v v T r - 1 ) | p ψ p 𝑑 x + c N ( v v T r - 1 ) p | ψ | p 𝑑 x
c r p N | v | p ( v T r - 1 ) p ψ p 𝑑 x + c N ( v v T r - 1 ) p | ψ | p 𝑑 x
(2.12) c r p ( N ( v v T r - 1 ) p | ψ | p 𝑑 x + N v q v T p ( r - 1 ) ψ p 𝑑 x ) .

Thus,

1 r p N | ( v v T r - 1 ψ ) | p 𝑑 x c N ( v v T r - 1 ) p | ψ | p 𝑑 x + c N v q - p ( v v T r - 1 ψ ) p 𝑑 x
c ( R - ρ ) p ( B R ( v v T r - 1 ) p * p p + p * - q 𝑑 x ) p + p * - q p *
+ c ( N v p * 𝑑 x ) q - p p * ( N ( v v T r - 1 ψ ) p * p p + p * - q 𝑑 x ) p + p * - q p *
(2.13) c ( R - ρ ) p ( B R ( v v T r - 1 ) p * d 𝑑 x ) p p * d ,

where d=pp+p*-q<1. Then by the Sobolev embedding theorem and (2.13), we have

( B ρ ( v v T r - 1 ) p * 𝑑 x ) p p * c N | ( v v T r - 1 ψ ) | p 𝑑 x c r p ( R - ρ ) p ( B R ( v v T r - 1 ) p * d 𝑑 x ) p p * d .

So

( B ρ ( v v T r - 1 ) p * 𝑑 x ) 1 p * c r R - ρ ( B R ( v v T r - 1 ) p * d 𝑑 x ) 1 p * d .

If BRvrp*d𝑑x<+, then

(2.14) ( B ρ v r p * 𝑑 x ) 1 r p * ( c r R - ρ ) 1 r ( B R v r p * d 𝑑 x ) 1 r p * d .

Let χ=1d, rj=χj,j=0,1,2,, and ρj=ρ+12j(R-ρ). Then by (2.14) we get

( B ρ j v r j p * 𝑑 x ) 1 r j p * ( c r j ρ j - 1 - ρ j ) 1 r j ( B ρ j - 1 v p * r j - 1 𝑑 x ) 1 r j - 1 p *
= ( c 2 j r j R - ρ ) 1 r j ( B ρ j - 1 v p * r j - 1 𝑑 x ) 1 r j - 1 p * .

We use Moser’s iteration to obtain

( B ρ j v r j p * 𝑑 x ) 1 r j p * j = 1 ( c 2 j χ j R - ρ ) 1 χ j ( B R v p * 𝑑 x ) 1 p * .

Letting j, we have

(2.15) v L ( B ρ ) c ( R - ρ ) d 1 - d v L p * ( B R ) .

For R<2 we prove vLp*(BR)cvLq(B2). Indeed, take r=1 in (2.12), ψ(x)=1 for |x|R, ψ(x)=0 for |x|2, and |ψ|22-R for all xN. Then

( B R v p * 𝑑 x ) p p * c N | ( v ψ ) | p 𝑑 x
c N v p | ψ | p 𝑑 x + c N v q ψ p 𝑑 x
c ( 2 - R ) p ( B 2 v q 𝑑 x ) p q + c ( B 2 v q 𝑑 x ) p q
c ( 2 - R ) p ( B 2 v q 𝑑 x ) p q ,

that is,

(2.16) v L p * ( B R ) c 2 - R v L q ( B 2 ) .

Combining (2.15) and (2.16), we obtain

| v | L ( B 1 ( x ) ) c | v | L q ( B 2 ( x ) ) for all  x N ,

where the constant c is independent of x, hence v(x)0 as |x|.

Step 2. Choose ε and R0 such that vq-p(x)εq-p<12a0 for |x|R0. Then for R>R0 we have

(2.17) N B R | v | p - 2 v φ d x + 1 2 a 0 N B R v p - 1 φ 𝑑 x 0 for all  φ W 0 1 , p ( N B R ) , φ 0 .

For any R>0 define ηR as ηR(x)=0 for |x|R, ηR(x)=1 for |x|R+1 and |ηR|2 for all xN. Taking φ=ηR2v in (2.17), we have

N B R | v | p - 2 v ( η R 2 v + 2 v η R η R ) 𝑑 x + 1 2 a 0 N B R v p η R 2 𝑑 x 0 .

Thus,

N B R + 1 ( | v | p + v p ) η R 2 𝑑 x c | N B R + 1 η R η R v | v | p - 2 v d x | c B R + 1 B R ( | v | p + v p ) 𝑑 x .

For R0>0 let Rn=R0+n. Then for any R>0 we have Rn<R<Rn+1. Define a sequence

a n = N \ B R 0 + n ( | v | p + v p ) 𝑑 x .

Then an+1c(an-an+1), which implies that an+1θan, where θ=cc+1<1. Thus there exists c (independent of n) such that ancθn. Hence,

(2.18) N B R ( | v | p + v p ) 𝑑 x c e - α R ,

and we have

N B R v q 𝑑 x c e - α R .

Step 3. Repeat the argument in Step 1. By (2.18) we have

v ( x ) c e - α | x | for all  x N .

Remark 2.5.

Suppose that vn0 and vn satisfies the differential inequality (2.3), n=1,2,. Suppose vnv in Lq(N). By checking the proof of Lemma 2.4, the sequence vn is uniformly bounded, that is, for some α>0 we have

v n ( x ) 1 α e - α | x | for all  x N .

Lemma 2.6.

Denote

Ω R ( n ) = N { B R k Λ B R ( x n , k ) } .

Then there exists α>0 such that

Ω R ( n ) ( | u n | p + | u n | p ) 𝑑 x 1 α e - α R , | u n ( x ) | 1 α e - α R for all  x Ω R ( n ) .

The proof of Lemma 2.6 is similar to that of Lemma 2.4, so we only sketch it.

Proof.

One can prove the lemma in four steps. Step 1. By Lemma 2.4 and property (4) of the decomposition (2.1), we have

Ω R ( n ) v n q 𝑑 x c N B R v q 𝑑 x + c k Λ N B R V k q 𝑑 x + c N | r n | q 𝑑 x c e - α R + o n ( 1 ) .

Step 2. We use Moser’s iteration to prove the L-estimate

| v n ( x ) | c e - α R + o n ( 1 ) for all  x Ω R ( n ) .

In particular, for any ε>0 there exist n0, R0 such that for nn0 there holds

| v n ( x ) | < ε for all  x Ω R 0 ( n ) .

In the following, we assume nn0, RR0.

Step 3. This is similar to Lemma 2.4. Choose ε, R0 such that vnq-p(x)εq-p<12a0 for xΩR0(n). Then for R>R0 we have

(2.19) Ω R ( n ) | v n | p - 2 v n φ d x + 1 2 a 0 Ω R ( n ) v n p - 1 φ 𝑑 x 0 for all  φ Ω R ( n ) , φ 0 .

For any R>0 define φR as φR(x)=0 for xΩR(n), φR(x)=1 for xΩR+1(n) and |φR|2. Taking φ=φR2vn in (2.19), we have

Ω R ( n ) | v n | p - 2 v n ( φ R 2 v n + 2 v n φ R φ R ) 𝑑 x + 1 2 a 0 Ω R ( n ) v n p φ R 2 𝑑 x 0 .

Thus,

Ω R + 1 ( n ) ( | v n | p + v n p ) φ R 2 𝑑 x c | Ω R + 1 ( n ) φ R φ R v n | v n | p - 2 v n d x |
c Ω R ( n ) Ω R + 1 ( n ) ( | v n | p + v n p ) 𝑑 x ,

where c>0 is independent of n, R. By the definition of φR, we have

Ω R + 1 ( n ) ( | v n | p + v n p ) 𝑑 x c c + 1 Ω R ( n ) ( | v n | p + v n p ) 𝑑 x ,

which implies that there exist positive constants c (independent of n, R) and α such that

(2.20) Ω R ( n ) ( | v n | p + v n p ) 𝑑 x c e - α R .

And we have

Ω R ( n ) v n q 𝑑 x c e - α R .

Step 4. Repeat the argument in Step 2. By (2.20) we have

v n ( x ) c e - α R for all  x Ω R ( n ) .

We follow the idea of [7] to derive a local Pohožaev-type identity with a form as in [5], which is much closer to our case. Choose ψC0(N), tN. Taking tuψ as test function in equation ($P_{\lambda}$) and integrating by parts, we obtain a Pohožaev-type identity

1 p N t a | u | p ψ 𝑑 x - N t x F λ ( x , u ) ψ 𝑑 x = - 1 p N | u | p t ψ d x + N | u | p - 2 t u u ψ d x
(2.21) - 1 p N a ( x ) | u | p t ψ d x + N F λ ( x , u ) t ψ d x .

Assume that unW1,p(N), n=1,2,, is a solution of problem ($P_{\lambda}$) with λ=λn0, unM. Then we have the decomposition (2.1).

Without loss of generality we may assume |xn,1|=min{|xn,k|,kΛ}. Denote xn=xn,1. According to [5] we can construct a sequence of cones Cn, having vertex 12xn and generated by a ball BRn(xn) as follows:

C n = { w N | w = 1 2 x n + λ ( x - 1 2 x n ) , x B R n ( x n ) , λ [ 0 , ) } ,

where Rn satisfies

r ^ k 0 | x n | 2 = r n R n k 0 r n = r ^ | x n | 2 , r ^ = 1 5 ( c ¯ + 1 ) ,

and c¯ is the constant in condition (a4), Λ={1,2,,k0}.

The cone Cn has the following property:

(2.22) C n { B 1 2 r n ( 0 ) k Λ B 1 2 r n ( x n , k ) } = .

We now apply the identity (2.21). Take u=un, t=tn=xn|xn| and ψ=χφR, where χ,φRC0(N) such that χ(x)=0 for xCn, χ(x)=1 for xCn and dist(x,Cn)1, φR(x)=1 for |x|R, and φR(x)=0 for |x|2R. Letting R, we obtain

1 p N t n a | u n | p χ 𝑑 x - N t n x F λ ( x , u n ) χ 𝑑 x
= - 1 p N | u n | p t n χ d x + N | u n | p - 2 t n u n u n χ d x
(2.23)     - 1 p N a ( x ) | u n | p t n χ d x + N F λ ( x , u n ) t n χ d x .

By (2.22) and the definition of χ, the support of χ is contained in the domain Ω=ΩR(n) with R=12rn-1. By Lemma 2.6, the right-hand side of (2.23) decays exponentially, say less than ce-α|xn|. By [5, Lemma 4.2], we have tna12ar for xCn. Similarly by Lemma 2.1 (5), we have tnxFλ(x,un)0 . Hence the left-hand side of (2.23) can be estimated as

1 p N t n a | u n | p χ 𝑑 x - N t n x F λ ( x , u n ) χ 𝑑 x
1 p N t n a | u n | p χ 𝑑 x 1 2 p inf B L ( x n ) a r B L ( x n ) | u n | p 𝑑 x
(2.24) m 4 p inf B L ( x n ) a r ,

where

m = lim n B L ( x n ) | u n | p d x = lim n B L ( 0 ) ( u n ( + x n , 1 ) ) p d x = B L ( 0 ) U 1 p d x > 0

for L large enough. Combining (2.23) and (2.24), we obtain

m 4 p inf B L ( x n ) a r c e - δ | x n | ,

which is impossible by (a3). Thus Λ= and by the decomposition (2.1) we have already proved the following proposition.

Proposition 2.7.

Let unW1,p(RN) be a solution of ($P_{\lambda}$) with λ=λn0, n=1,2,. Assume unM, unu in W1,p(RN). Then {un} converges in Lq(RN).

Proof of Theorem 1.1.

We use an indirect argument. Assume for any μn=1n that there exist unW1,p(N), xnN such that un is a solution of (Pλn), unM. But there holds that

(2.25) | u n ( x n ) | > μ n e - μ n 1 + | x n | 2 .

By Proposition 2.7, the sequence {un} converges in Lq(N). By Remark 2.5, we have

| u n ( x ) | = v n ( x ) 1 α e - α 1 + | x | 2

for some α>0, which contradicts (2.25), provided μn=1n<α. ∎

Corollary 2.8.

Let unW1,p(RN) be a solution of ($P_{\lambda}$) with λ=λn0, n=1,2,. Assume Iλn(un)M, where

(2.26) I λ ( u ) = 1 p N ( | u | p + a ( x ) | u | p ) 𝑑 x - N F λ ( x , u ) 𝑑 x .

Then, up to a subsequence, {un} converges in W1,p(RN) and there exists μ>0 such that

(2.27) | u n ( x ) | 1 μ e - μ 1 + | x | 2 .

Proof.

By Lemma 2.1 (4), we have

M I λ n ( u n ) = I λ n ( u n ) - 1 r D I λ n ( u n ) , u n
= ( 1 p - 1 r ) N ( | u n | p + a ( x ) | u n | p ) 𝑑 x + q - r q r N | u n | r + 1 | m λ ( x , u n ) | q - r - 1 b λ ( x , u n ) 𝑑 x
( 1 p - 1 r ) N ( | u n | p + a ( x ) | u n | p ) 𝑑 x .

Hence {un} is bounded in W1,p(N). By Theorem 1.1, inequality (2.27) holds. By Proposition 2.7, up to a subsequence, {un} converges in Lq(N). Thus,

N ( | u n | p - 2 u n - | u m | p - 2 u m , u n - u m ) 𝑑 x + N a ( x ) ( | u n | p - 2 u n - | u m | p - 2 u m ) ( u n - u m ) 𝑑 x
= N ( f λ n ( x , u n ) - f λ m ( x , u m ) ) ( u n - u m ) 𝑑 x
c N ( | u n | q - 1 + | u m | q - 1 ) | u n - u m | 𝑑 x 0 as  n , m .

By the elementary inequalities (2.5) and (2.6), the sequence {un} converges in W1,p(N). ∎

3 Existence of Multiple Solutions

Problem ($P_{\lambda}$) has a variational structure, given by the functional Iλ defined on W1,p(N) (see (2.26)).

Lemma 3.1.

I λ ( λ 0 ) satisfies the Palais–Smale condition.

Proof.

Let {un}W1,p(N) be a Palais–Smale sequence. By Lemma 2.1 (4), we have

N ( 1 r u n f λ ( x , u n ) - F λ ( x , u n ) ) 𝑑 x 0 .

Then

I λ ( u n ) - 1 r D I λ ( u n ) , u n = ( 1 p - 1 r ) N ( | u n | p + a ( x ) | u n | p ) 𝑑 x + N ( 1 r u n f λ ( x , u n ) - F λ ( x , u n ) ) 𝑑 x
( 1 p - 1 r ) N ( | u n | p + a ( x ) | u n | p ) 𝑑 x .

Hence {un} is bounded in W1,p(N). Assume unu in W1,p(N), unu in Llocq(N), 1q<p*. By Lemma 2.1, we obtain

N ( | u n | p - 2 u n - | u m | p - 2 u m , u n - u m ) 𝑑 x + N a ( x ) ( | u n | p - 2 u n - | u m | p - 2 u m ) ( u n - u m ) 𝑑 x
= N ( f λ ( x , u n ) - f λ ( x , u m ) ) ( u n - u m ) 𝑑 x N ( 2 λ e - λ 1 + | x | 2 ) q - r ( | u n | r - 1 + | u m | r - 1 ) | u n - u m | 𝑑 x
c e - ( q - r ) λ R N B R ( | u n | r + | u m | r ) 𝑑 x + c ( B R ( | u n | r + | u m | r ) 𝑑 x ) r - 1 r ( B R | u n - u m | r 𝑑 x ) 1 r
(3.1) c e - ( q - r ) λ R + c ( B R | u n - u m | r 𝑑 x ) 1 r 0 as  n , m .

By (3.1) and the elementary inequalities (2.5), (2.6), we have that {un} is a Cauchy sequence in W1,p(N). ∎

Now we define a sequence of critical values of Iλ.

c k ( λ ) = inf A Γ k sup u A I λ ( u ) , λ 0 , k = 1 , 2 , ,

where

Γ k = { A A W 1 , p ( N ) , A  compact , - A = A , γ ( A σ - 1 ( S ρ ) ) k  for all  σ G } ,
G = { σ σ C ( W 1 , p ( N ) , W 1 , p ( N ) ) , σ ( - u ) = - σ ( u )  for all  u W 1 , p ( N ) , σ ( u ) = u  if  I 1 ( u ) < 0 } ,
S ρ = { u u W 1 , p ( N ) , u = ρ } for a suitable  ρ > 0 .

For uSρ we have

I 1 ( u ) = 1 p N ( | u | p + a ( x ) | u | p ) 𝑑 x - N F 1 ( x , u ) 𝑑 x
1 p N ( | u | p + a ( x ) | u | p ) 𝑑 x - 1 q N | u | q 𝑑 x
c 0 ( ρ p - ρ q ) .

Choosing ρq-p<12, we have I1(u)m=12c0ρp for uSρ.

The following proposition is known (see [1, 2, 11]).

Proposition 3.2.

Assume 0<λ1. Then we have the following results:

  1. c k ( λ ) > 0 , k=1,2,, are critical values of Iλ.

  2. If c k ( λ ) = c k + 1 ( λ ) = = c k + m - 1 ( λ ) = c , then γ ( K c ( I λ ) ) m , where

    K c ( I λ ) = { u u W 1 , p ( N ) , D I λ ( u ) = 0 , I λ ( u ) = c } .

  3. Assume p = 2 . Then there exists u W 1 , p ( N ) such that I λ ( u ) = c k ( λ ) , DIλ(u)=0 and m*(u)k, where m*() is the augmented Morse index.

Proof of Theorem 1.2.

Given kN, by Corollary 2.8 there exists μk>0 such that if uW1,p(N), DIλ(u)=0 and Iλ(u)dk:=ck(1), then

| u ( x ) | 1 μ k e - μ k 1 + | x | 2 for all  x N .

Choose 0<λk<min{1,μk}. Let u1(λk),,uk(λk) be k solutions of problem ($P_{\lambda}$) with λ=λk, corresponding to the critical values c1(λk)ck(λk). Since Iλ is increasing as λ, we have that c1(λk)ck(λk)dk, u1(λk),,uk(λk) satisfy the estimate

| u ( x ) | 1 μ k e - μ k 1 + | x | 2 1 λ k e - λ k 1 + | x | 2 ,

hence u1(λk),,uk(λk) are solutions of the original problem ($P_{0}$), k is arbitrary and ($P_{0}$) has infinitely many solutions. ∎

Remark 3.3.

We have proved that problem ($P_{0}$) has infinitely many solutions. We can prove a little more.

Claim.

The functional I has an infinite sequence of critical values c1<c2<<ck<. Moreover, if p=2, then limkck=+.

We use an indirect argument. Assume I has only finite critical values c1,,ck. Denote

K = { u u W 1 , p ( N ) , D I ( u ) = 0 } .

Then by Corollary 2.8, the set K is compact. Assume γ(K)=m<+. For 0<λ<1, the functional Iλ has critical values c1(λ)c2(λ)ckm+1(λ). If λ is small enough, they will be critical values of the functional I. We claim c1(λ)<cm+1(λ)<<ckm+1(λ). Otherwise suppose say c=c1(λ)=cm+1(λ). By Proposition 3.2, we have γ(Kc)m+1, where

K c = { u u W 1 , p ( N ) , D I λ ( u ) = 0 , I λ ( u ) = c } K .

This is a contradiction.

Now assume p=2. As in the proof of Theorem 1.2, we have that uk(λk) is a solution of problem ($P_{0}$). But by Proposition 3.2, we can assume m*(uk(λk))k. If the critical values of I are bounded above, then K is compact by Corollary 2.8, and the augmented Morse index of the critical points will be bounded above. We arrive at a contradiction.

Remark 3.4.

In problem ($P_{0}$) the positive nonlinear term has growth p-1. If we replace this term by a(x)|u|l-2u with l(1,q), the main theorems remain true.

4 p-Laplacian Equations with Critical Sobolev Growth

In this section, we briefly discuss how to use the truncation method to deal with the p-Laplacian equations with critical Sobolev growth. We consider the p-Laplacian equation

(4.1) { - Δ p u = | u | p * - 2 u + μ | u | p - 2 u in  Ω , u = 0 on  Ω ,

where 1<p<N, p*=NpN-p, μ>0, and Ω is an open bounded domain in N.

As in the introduction we introduce auxiliary functions. Let φC0() be such that φ(t)=1 for |t|1, φ(t)=0 for |t|2, φ is decreasing in [1,2] and φ is even. For λ>0 define

b λ ( s ) = φ ( λ s ) , m λ ( s ) = 0 s b λ ( τ ) 𝑑 τ ,
F λ ( s ) = 1 p * | s | r | m λ ( s ) | p * - r , f λ ( s ) = d d s F λ ( s ) ,

where r(p,p*) is a fixed number, say r=12(p+p*). For λ=0 we obtain that m0(s)=s, F0(s)=F(s)=1p*|s|p* and f0(s)=f(s)=|s|p*-2s. Besides of (4.1), we consider

(4.2) { - Δ p u = f λ ( u ) + μ | u | p - 2 u in  Ω , u = 0 on  Ω .

As a problem of subcritical growth, (4.2) has infinitely many solutions, thus we have the following convergence theorem.

Theorem 4.1.

Assume N>p2+p. Let λn0, unW01,p(Ω) be a solution of (4.2) with λ=λn, and unc. Then, up to a subsequence, unu in W01,p(Ω), and there exists a constant ν>0 such that

| u n ( x ) | 1 ν .

In particular, for λn<ν, the sequence un will be a solution of the original problem (4.1).

Theorem 4.2.

Assume N>p2+p. Then problem (4.1) has infinitely many solutions.

We are going to discuss the problem of critical growth in detail in a forthcoming paper.

Remark 4.3.

We mentioned problem (4.1) in the introduction. Cao, Peng and Yan [4] considered the equations of subcritical growth

(4.3) { - Δ p u = | u | p * - 2 - ε u + μ | u | p - 2 u in  Ω , u = 0 on  Ω ,

where ε is a small constant. The solutions of (4.3) will be used as the approximated solutions. By following the idea of [7] and overcoming some difficulties due to the p-Laplacian operator, they proved a convergence theorem to obtain solutions of the original problem.


Communicated by Paul Rabinowitz


Funding statement: This work was supported by NSFC 11361077 and Yunnan Province, Young Academic and Technical Leaders Program (2015HB028) and Yunnan Normal University, LianDa Scholar Program.

Acknowledgements

The authors are grateful to the referee for a careful reading of the manuscript, clarifying some details and pointing out some references. In particular, we thank the referee for drawing our attention to Stampacchia lattice properties of W1,p, by which Lemma 2.2 can be proved more directly.

References

[1] A. Ambrosetti and P. H. Rabinowitz, Dual variational methods in critical point theory and applications, J. Funct. Anal. 14 (1973), 349–381. 10.1016/0022-1236(73)90051-7Search in Google Scholar

[2] A. Bahri and P. L. Lions, Morse index of some min-max critical points. I: Application to multiplicity results, Comm. Pure Appl. Math. 61 (1988), 1027–1037. 10.1002/cpa.3160410803Search in Google Scholar

[3] H. Brezis and L. Nirenberg, Positive solutions of nonlinear elliptic equations involving critical Sobolev exponents, Comm. Pure Appl. Math. 36 (1983), 437–477. 10.1002/cpa.3160360405Search in Google Scholar

[4] D. Cao, S. Peng and S. Yan, Infinitely many solutions for p-Laplacian equation involving critical Sobolev growth, J. Funct. Anal. 262 (2012), 2861–2902. 10.1016/j.jfa.2012.01.006Search in Google Scholar

[5] G. Cerami, G. Devillanova and S. Solimini, Infinitely many bound states solutions for some nonlinear scalar field equations, Calc. Var. Partial Differential Equations 23 (2005), 139–168. 10.1007/s00526-004-0293-6Search in Google Scholar

[6] L. Damascelli, Comparision theorems for some quasilinear degenerate elliptic operators and applications to symmetry and monotonity results, Ann. Inst. H. Poincaré Anal. Non Linéaire 15 (1998), 493–516. 10.1016/s0294-1449(98)80032-2Search in Google Scholar

[7] G. Devillanova and S. Solimini, Concentrations estimates and multiple solutions to elliptic problems at critical growth, Adv. Differential Equations 7 (2002), 1257–1280. 10.57262/ade/1356651637Search in Google Scholar

[8] J. Q. Liu and Z. Q. Wang, Multiple solutions for quasilinear elliptic equations with a finite potential well, J. Differential Equation 257 (2014), 2874–2899. 10.1016/j.jde.2014.06.002Search in Google Scholar

[9] P. H. Rabinowitz, Free vibrations for a semilinear wave equation, Comm. Pure Appl. Math. 31 (1978), 31–68. 10.1002/cpa.3160310103Search in Google Scholar

[10] P. H. Rabinowitz, Periodic solutions of Hamiltonian systems, Comm. Pure Appl. Math. 31 (1978), 157–184. 10.21236/ADA096642Search in Google Scholar

[11] S. Solimini, Morse index estimates in min-max theorems, Manuscripta Math. 63 (1989), 421–453. 10.1007/BF01171757Search in Google Scholar

[12] C. Tintarev, Concentration analysis and cocompactness, Concentration Analysis and Applications to PDE (Bangalore 2012), Trends Math., Birkhäuser, Basel (2013), 117–141. 10.1007/978-3-0348-0373-1_7Search in Google Scholar

[13] C. Tintarev and K.-H. Fieseler, Concentration Compactness, Functional-Analytic Grounds and Applications, Imperial College Press, London, 2007. 10.1142/p456Search in Google Scholar

Received: 2016-03-22
Revised: 2016-07-26
Accepted: 2016-10-02
Published Online: 2016-11-08
Published in Print: 2017-07-01

© 2017 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 27.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/ans-2015-5059/html
Scroll to top button