Show Summary Details
More options …

Editor-in-Chief: Radulescu, Vicentiu / Squassina, Marco

IMPACT FACTOR 2018: 6.636

CiteScore 2018: 5.03

SCImago Journal Rank (SJR) 2018: 3.215
Source Normalized Impact per Paper (SNIP) 2018: 3.225

Mathematical Citation Quotient (MCQ) 2017: 1.89

Open Access
Online
ISSN
2191-950X
See all formats and pricing
More options …
Volume 9, Issue 1

# Monotonicity formulas for coupled elliptic gradient systems with applications

Mostafa Fazly
/ Henrik Shahgholian
Published Online: 2019-06-06 | DOI: https://doi.org/10.1515/anona-2020-0010

## Abstract

Consider the following coupled elliptic system of equations

$(−Δ)sui=(u12+⋯+um2)p−12uiin Rn,$

where 0 < s ≤ 2, p > 1, m ≥ 1, u = $\begin{array}{}\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ and ui : ℝn → ℝ. The qualitative behavior of solutions of the above system has been studied from various perspectives in the literature including the free boundary problems and the classification of solutions. For the case of local scalar equation, that is when m = 1 and s = 1, Gidas and Spruck in [26] and later Caffarelli, Gidas and Spruck in [6] provided the classification of solutions for Sobolev sub-critical and critical exponents. More recently, for the case of local system of equations that is when m ≥ 1 and s = 1 a similar classification result is given by Druet, Hebey and Vétois in [17] and references therein. In this paper, we derive monotonicity formulae for entire solutions of the above local, when s = 1, 2, and nonlocal, when 0 < s < 1 and 1 < s < 2, system. These monotonicity formulae are of great interests due to the fact that a counterpart of the celebrated monotonicity formula of Alt-Caffarelli-Friedman [1] seems to be challenging to derive for system of equations. Then, we apply these formulae to give a classification of finite Morse index solutions. In the end, we provide an open problem in regards to monotonicity formulae for Lane-Emden systems.

MSC 2010: 35J60; 35J50; 35B35; 35B45

## 1.1 Background

Our objective in this paper is to establish monotonicity formulae for solutions of the following coupled elliptic system of equations

$(−Δ)sui=|u|p−1uiin Rn$(1.1)

and classify finite Morse index solutions when 0 < s ≤ 2, p > 1, m ≥ 1, u = $\begin{array}{}\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ and ui : ℝn → ℝ. Note that when m = 1 the above is known as the Lane-Emden equation. An important feature of system (1.1) is the variational structure of the coupled equations of the form

$−Δu=∇H(u)in Rn.$(1.2)

Note that (1.1) is a particular case of (4.33) for H(u) = $\begin{array}{}\frac{1}{p+1}\end{array}$ |u|p+1. Due to this variational structure, qualitative behaviour of solutions of system (4.33) has been studied extensively in the context of elliptic partial differential equations from both pure and applied mathematics perspectives. Let us mention that Andersson et al. in [4] considered system (1.1), up to a negative sign, when p = 0 and s = 1 which reads

$Δui=ui|u|χ{|u|>0}.$(1.3)

This is closely related to minimizers of the energy

$∫∑i=1m|∇ui|2+2|u|dx,$(1.4)

and they studied the regularity of free boundaries. To prove regularity results in [4], they established a monotonicity formula that is inspired by the one given by Weiss in [44, 45]. We also refer interested readers to [3, 7] for regularity results on cooperative systems and to the book of Petrosyan, Shahgholian and Uraltseva in [36] for more information. It should be remarked that as a particular case, one can consider the case of two components, that is m = 2, and p = 3 namely

$(−Δ)su=(u2+v2)u in Rn,(−Δ)sv=(u2+v2)v in Rn.$(1.5)

Note that the above system is a special case of the nonlinear Schrödinger system that is

$(−Δ)su=(μ1u2+βv2)u in Rn,(−Δ)sv=(μ2u2+βv2)v in Rn,$(1.6)

where μ1, μ2, β are constants. We would like to mention that most of the results in the paper are valid for (1.6) as well. We refer interested readers to [34, 41, 42] and references therein for more information regarding classification of solutions of system (1.6) for various parameters μ1, μ2 and β.

## 1.2 Tools and methodology (Historical development)

Semilinear elliptic equations (with almost four decades of history) are probably the most widely and intensely studied equations in PDE. The simplest equation in this class is expressed as

$Δu=f(u),$(1.7)

and exhibits as many possible features, as the right hand side nonlinearity f(u) may do. This equation has also been studied from so many different perspectives, that any attempt to try to list them here would probably fail. There are however two different types of questions, of interest to the current authors, that seem to be interconnected and developed in parallel, but almost separately:

1. Analysis of the singular set: {x : ∇u = 0}.

2. Stability and Liouville’s type questions.

One can actually add a third less studied, and probably more complicated problem to the above, which is:

3. Structure of the blow-up set: {x : |u(x)| = ∞}.

In studying the above problems, experts have (in many cases) developed parallel tools and ideas to handle technical passages of the analysis of questions in each cases. Two of these tools, also utilized in this paper, are Monotonicity functional, and Blow-up/down Analysis.

The particular monotonicity functional used here combines the energy functional along with a balanced term (see below for explicit form). The use of monotonicity functional (in the way we present here) can be traced back1 to the work of Fleming in [25] for area minimizing currents which has subsequently been proved by Allard in [2] for stationary rectifiable n-varifolds, that was later developed further by others. The monotonicity functional that we shall use in this paper originates in the study of harmonic maps by Price in [37] and Schoen and Uhlenbeck in [40], see also [39]. And in connection to free boundary problems, it was first used by Ou [5], and developed later by Weiss in [45] and others in various forms. In the context of semilinear heat equations, a similar monotonicity formula is given by Giga and Kohn in [28] and for the corresponding elliptic equation by Pacard [35].

The second main tool, blow-up/down analysis, has its origin in local regularity theory and the so-called linearization technique (also called harmonic blow-up). To study local structure of level surfaces of a solution to the above semilinear problem, one considers scaled version of the problem and classifies the limit manifold of such scalings. The limit manifold naturally carries information of the local behavior of the solutions and hence one can then with some devices link back this behavior to the local problem and deduce the expected properties. It is now that the role of monotonicity functional become crucial and indispensable in that one can use the fact that when blowing up a solution, the monotonicity functional, call it E(r, u), being monotone will have a limit. At the same time this functional has a nice scaling properties E(rs, u) = E(s, ur) for r, s > 0 and ur a correctly scaled version of u. Hence one obtains

$E(0+,u)=limr→0E(rs,u)=limr→0E(s,ur)=E(s,u0),$(1.8)

so that E(s, u0) is constant. Now a strong version of the monotonicity functional asserts that the only time E is constant is when the function u0 is homogeneous, where the order of homogeneity is dictated by E. Since the same can be done by blowing-down, we will have that u0, u are homogeneous of same order. Next if one can prove that the homogenous solutions are unique we must then have u0 = u and hence E(0+, u) = E(∞, u), and we arrive at E(r, u) = constant. Therefore u is itself homogeneous, and already classified.

In the regularity theory of level surfaces the classification of homogeneous global solutions is one of the key elements of the theory, and many time a very hard nut to crack, in higher dimensions; in two dimensions homogenous solutions can be computed directly. In the stability theory and Liouville type problems this part is based on computations of eigenvalues of the corresponding Laplace-Beltrami on the sphere, as also done in this paper. As shown here below (and in many preceding papers by experts in this area) the only homogenous solutions in the appropriate space are the trivial solutions. In particular, the energy functional becomes zero and hence E(r, u) = 0. From here one deduces that u ≡ 0 is the only solution satisfying the given condition.

It is worthwhile remarking the fundamental feature of these problems, that in some cases depend on the dimension and order of homogeneity are directly a consequence of the eigenvalues of Laplace-Beltrami, where in the case of index-theory results in the computation of eigenvalues for nonlinear Laplace-Beltrami, that in turn implies that non-trivial solutions do not exists in certain range of values. These ideas have their origin in the regularity theory of minimal surfaces, that boils down to proving minimal cones do not exists in dimensions less than eight.

## 1.3 Problem setting

Our main technique is to derive monotonicity formulae for solutions of (1.1) for various values of parameter 0 < s ≤ 2. To provide such monotonicity formulae we consider various cases.

Case s = 1. Consider the following functional for every λ > 0 and x0 ∈ ℝn

$E1(u,λ,x0):=λ−n+2p+1p−1∫Bλ(x0)12∑i=1m|∇ui|2−1p+1|u|p+1+1p−1λ−n+2p+1p−1−1∫∂Bλ(x0)|u|2,$(1.9)

then the following monotonicity formula holds for classical solutions of (1.1).

#### Theorem 1.1

Suppose that u is a solution of (1.1) for s = 1. Then, E1 is a nondecreasing function of λ and in fact

$dE1dλ=λ−n+1+2p+1p−1∫∂Bλ(x0)∑i=1m∂ui∂r+2p−1uir2,$(1.10)

where $\begin{array}{}\frac{\mathrm{\partial }}{\mathrm{\partial }r}\end{array}$ is polar derivative.

Note that for the case of single equations, that is when m = 1, similar monotonicity formulae are given by Pacard in [35] and Weiss in [44, 45]. For the case of systems that is when m ≥ 1, very recently, Andersson, Shahgholian, Uraltseva and Weiss in [4] provided a monotonicity formula for solutions of system (1.1) when p = 0 and applied it to study free boundary problems. Let us mention that very similar monotonicity formulae appear in the field of harmonic maps that is

$−Δu=|∇u|2u,$(1.11)

where u : ℝn → 𝕊m–1, see Evans in [18, 19] and reference therein.

Case s = 2. For every λ > 0 and x0 ∈ ℝn define

$E2(u,λ,x0):=λ4p+1p−1−n∫Bλ(x0)12∑i=1m|Δui|2−1p+1|u|p+1−4p−1p+3p−1−nλ1+8p−1−n∫∂Bλ(x0)|u|2−4p−1p+3p−1−nddλλ8p−1+2−n∫∂Bλ(x0)|u|2+12λ3ddλλ8p−1+1−n∫∂Bλ(x0)∑i=1m4p−1λ−1ui+∂ui∂r2+12ddλλ4p+1p−1−n∫∂Bλ(x0)∑i=1m|∇ui|2−∂ui∂r2+12λ4p+1p−1−n−1∫∂Bλ(x0)∑i=1m|∇ui|2−∂ui∂r2.$(1.12)

Then the following monotonicity formula holds for the fourth order Lane-Emden system.

#### Theorem 1.2

Suppose that n ≥ 5, p > $\begin{array}{}\frac{n+4}{n-4}\end{array}$ and u is a solution of (1.1) when s = 2. For any λ > 0 and x0 ∈ ℝn

$dE2(u,λ,x0)dλ≥Cλ8p−1+2−n∫∂Bλ(x0)∑i=1m4p−1r−1ui+∂ui∂r2,$(1.13)

where E2 is defined by (1.12) and C is independent from λ.

Note also that the above monotonicity formula is provided by Davila, Dupaigne, Wang and Wei in [15] for the case of single equations that is when m = 1. For the case of fractional Laplacian, we provide monotonicity formulae for the extension problems. Assume that uiC2σ(ℝn), σ > s > 0 and

$∫Rn|ui(y)|(1+|y|)n+2sdy<∞,$(1.14)

for each 1 ≤ im. The fractional Laplacian of ui when 0 < s < 1 denoted by

$(−Δ)sui(x):=p.v.∫Rnui(x)−ui(y)|x−y|n+2sdy,$(1.15)

is well-defined for every x ∈ ℝn. Here p.v. stands for the principle value. It is by now standard that the fractional Laplacian can be seen as a Dirichlet-to-Neumann operator for a degenerate but local diffusion operator in the higher-dimensional half-space ℝn+1, see Caffarelli and Silvestre in [8]. In other words, for uiC2σL1(ℝn, (1+|y|n+2s)dy) when σ > s and 0 < s < 1, there exists v = $\begin{array}{}\left({v}_{i}{\right)}_{i=1}^{m}\end{array}$ such that $\begin{array}{}{v}_{i}\in {C}^{2}\left({\mathbb{R}}_{+}^{n+1}\right)\cap C\left(\overline{{\mathbb{R}}_{+}^{n+1}}\right),{y}^{1-2s}{\mathrm{\partial }}_{y}{v}_{i}\in C\left(\overline{{\mathbb{R}}_{+}^{n+1}}\right)\end{array}$ and

$∇⋅(y1−2s∇vi)=0in R+n+1,vi=uion ∂R+n+1,−limy→0y1−2s∂yvi=κs|v|p−1vion ∂R+n+1,$(1.16)

for the following constant κs,

$κs:=Γ(1−s)22s−1Γ(s).$(1.17)

For the case of 1 < s < 2, there are various definitions for the fractional operator (–Δ)s, see [9, 12, 24, 47]. From the Fourier transform one can define the fractional Laplacian by

$(−Δ)s^ui(ζ)=|ζ|2su^i(ζ),$(1.18)

or inductively by (–Δ)s = (–Δ)s–1 o(–Δ). Note that Yang in [47] gave a counterpart of the extension problem by Caffarelli and Silvestre in [8] for the fractional Laplacian (–Δ)s, where s is any positive, noninteger number. In other words, he showed that the higher order fractional Laplacian operator can also be regarded as the Dirichlet-to-Neumann map for an extension function satisfying a higher order elliptic equation in the upper half space with one extra spatial dimension. More precisely, there exists an extension function viW2,2($\begin{array}{}{\mathbb{R}}_{+}^{n+1}\end{array}$, yb) such that

$Δb2vi=0 in R+n+1,limy→0yb∂yvi=0 in ∂R+n+1,limy→0yb∂yΔbvi=Cn,s|v|p−1vi in Rn,$(1.19)

where b := 3 – 2s and Δb vi := yb div(ybvi). We refer interested readers to Case and Chang [9] and Chang and Gonzales [12] as well. We are now ready to provide monotonicity formulae for various parameters s > 0. Case 0 < s < 1. Let v = $\begin{array}{}\left({v}_{i}{\right)}_{i=1}^{m}\end{array}$ be a solution of the extension problem (1.16). Now define the energy functional for any λ > 0 and x0$\begin{array}{}\mathrm{\partial }{\mathbb{R}}_{+}^{n+1}\end{array}$ as

$Es(v,λ,x0):=λ2s(p+1)p−1−n12∫R+n+1∩Bλy1−2s∑i=1m|∇vi|2dxdy−κsp+1∫∂R+n+1∩Bλ|v|p+1dx+λ2s(p+1)p−1−n−1sp+1∫∂Bλ∩R+n+1y1−2s|v|2dσ.$(1.20)

We are now ready to provide a monotonicity formula for solutions of (1.16) when 0 < s < 1.

#### Theorem 1.3

Suppose that 0 < s < 1. Let v = $\begin{array}{}\left({v}_{i}{\right)}_{i=1}^{m}\end{array}$ where each $\begin{array}{}{v}_{i}\in {C}^{2}\left({\mathbb{R}}_{+}^{n+1}\right)\cap C\left(\overline{{\mathbb{R}}_{+}^{n+1}}\right)\end{array}$ be a solution of (1.16) such that y1–2s y viC($\begin{array}{}\overline{{\mathbb{R}}_{+}^{n+1}}\end{array}$). Then, Es is a nondecreasing function of λ. Furthermore,

$dEsdλ=λ2s(p+1)p−1−n+1∫∂B(x0,λ)∩R+n+1y1−2s∑i=1m∂vi∂r+2sp−1vir2dσ,$(1.21)

where Es provided in (1.20).

Case 1 < s < 2. Suppose that v = $\begin{array}{}\left({v}_{i}{\right)}_{i=1}^{m}\end{array}$ is a solution of the extension problem (1.19). Similarly, now define the energy functional for any λ > 0 and x0$\begin{array}{}\mathrm{\partial }{\mathbb{R}}_{+}^{n+1}\end{array}$ as

$Es(v,λ,x0):=λ2sp+1p−1−n∫R+n+1∩Bλ(x0)12y3−2s∑i=1m|Δbvi|2−Cn,sp+1∫∂R+n+1∩Bλ(x0)|v|p+1−sp−1p+2s−1p−1−nλ−3+2s+4sp−1−n∫R+n+1∩∂Bλ(x0)y3−2s|v|2−sp−1p+2s−1p−1−nddλλ4sp−1+2s−2−n∫R+n+1∩∂Bλ(x0)y3−2s|v|2+12r3ddλλ4sp−1+2s−3−n∫R+n+1∩∂Bλ(x0)y3−2s∑i=1m2sp−1λ−1vi+∂vi∂r2+12ddλλ2sp+1p−1−n∫R+n+1∩∂Bλ(x0)y3−2s∑i=1m|∇vi|2−∂vi∂r2+12λ2sp+1p−1−n−1∫R+n+1∩∂Bλ(x0)y3−2s∑i=1m|∇vi|2−∂vi∂r2.$(1.22)

Considering the above energy functional, we now provide a monotonicity formula for solutions of (1.19) when 1 < s < 2.

#### Theorem 1.4

Assume that n > 2s and p > $\begin{array}{}\frac{n+2s}{n-2s}\end{array}$. Let v = $\begin{array}{}\left({v}_{i}{\right)}_{i=1}^{m}\end{array}$ be a solution of (1.19) then E(v, λ, x0) is a nondecreasing function of λ > 0. In addition,

$dEs(v,λ,x0)dλ≥C(n,s,p)λ4sp−1+2s−2−n∫R+n+1∩∂Bλ(x0)y3−2s∑i=1m2sp−1r−1vi+∂vi∂r2,$(1.23)

where Es is given by (1.22) and C(n, s, p) is independent from λ.

We also refer interested readers to Davila, Dupaigne and Wei in [14] and to Wei and the first author in [24] for a similar monotonicity formula for the case of scalar equations and 0 < s < 1 and 1 < s < 2, respectively. Before we state our main results let us present some backgrounds regarding classification of solutions of (1.1) in the absence of stability. We provide such classifications for scalar equations and systems separately.

#### Remark 1.1

Note that monotonicity formulae provided as Theorems 1.1-1.4 hold for the following system with a slightly more general right-hand side,

$(−Δ)sui=|u|p−1αiui++βiui−in Rn,$(1.24)

where αi and βi are positive constants for 1 ≤ im.

Suppose that m = 1 when (1.1) turns into a single equation. We first consider the local operator cases meaning s = 1 and s = 2. Suppose that s = 1 and the parameter p is in the subcritical case that is when 1 < p < pS(n, 1) where

$pS(n,1):=∞ if n≤2,n+2n−2 if n>2.$(1.25)

For this case, there is a very well-known classification result of Gidas and Spruck in [26] stating that the only nonnegative solution of the Lane-Emden equation (1.1) with s = 1 is the trivial solution. For the critical case, that is when p = pS(n, 1), Caffarelli, Gidas and Spruck in [6] that there is a unique (up to translation and rescaling) positive solution for the Lane-Emden equation. For the fourth order Lane-Emden equation, that is when s = 2, Wei and Xu in [43] provided a similar classification result for the subcritical 1 < p < pS(n, 2) and the critical cases p = pS(n, 2) when

$pS(n,2):=∞ if n≤4,n+4n−4 if n>4,$(1.26)

see also [33]. Note that for the fractional Laplacian operator 0 < s < 1, such a classification result is given by Li [32] and Chen-Li-Ou [11] where the critical exponent is

$pS(n,s)=∞ if n≤2s,n+2sn−2s if n>2s.$(1.27)

For the case of system of equations, that is when m ≥ 2, Druet, Hebey and Vétois in [17] provided a classification result for solutions of (1.1) where s = 1 via the moving sphere method. Suppose that p = $\begin{array}{}\frac{n+2}{n-2}\end{array}$ and u = $\begin{array}{}\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ is a nonnegative classical solution of (1.1) where s = 1. Then they proved that there exist x0 ∈ ℝn, λ > 0, Λ$\begin{array}{}{\mathcal{S}}_{+}^{m-1}\end{array}$ such that

$u(x)=λλ2+1n(n−2)|x−x0|2n−22Λ.$(1.28)

We also refer interested readers to [16, 29, 30] where the authors studied various counterparts of system (1.1). Note that the following singular function

$us(x)=A|x|−2sp−1for A∈Rm with |A|p−1=22sΓ(n2−sp−1)Γ(s+sp−1)Γ(sp−1)Γ(n−2s2−sp−1),$(1.29)

solves (1.1) in ℝn ∖ {0} for a supercritical parameter p that is p > pS(n, s). Before we state our main results, let us define the notion of stable solutions.

#### Definition 1.1

We say a solution u of (1.1) is stable outside a compact set if there exists R0 > 0 such that

$∑i=1m∫|u|p−1ϕi2+(p−1)∑i,j=1m∫|u|p−3uiujϕiϕj≤∑i=1m||ϕi||H˙s(Rn)2,$(1.30)

for any ϕ = $\begin{array}{}\left({\varphi }_{i}{\right)}_{i=1}^{m}\end{array}$ where ϕi$\begin{array}{}{C}_{c}^{\mathrm{\infty }}\end{array}$(ℝnBR0) for 1 ≤ im.

Here is our main result.

#### Theorem 1.5

Suppose 0 < s ≤ 2 and n > 2s. Let u be a solution of (1.1) that is stable outside a compact set. Then either for 1 < p < pS(n, s) or for p > pS(n, s) and

$pΓ(n2−sp−1)Γ(s+sp−1)Γ(sp−1)Γ(n−2s2−sp−1)>Γ(n+2s4)2Γ(n−2s4)2,$(1.31)

each component ui must be identically zero. For the case of Sobolev critical exponent, that is when p = pS(n, s), a solution u has finite energy that is

$||u||Lp+1(Rn)p+1=||u||H˙s(Rn)2<∞.$(1.32)

In this case, if in addition u is stable, then each component ui must be identically zero.

As a direct consequence, the above theorem implies that the only nonnegative solution for system (1.5) when s = 1 is the trivial solution for dimensions n < 12 and n ≠ 4. Here, is how this article is structured. For the rest of this article, we provide a proof for Theorem 1.5 considering various cases for parameter s > 0. In Section 2, we consider the case s = 1 that is when the operator in (1.1) is the local Laplacian operator. In Section 3, we let s = 2 that refers to the bi-Laplacian operator. Lastly, in Section 4, we consider nonlocal cases 0 < s < 1 and 1 < s < 2. For these non-integer parameters, the operator in (1.1) is a fractional Laplacian operator. For all cases 0 < s ≤ 2, we apply blow-down analysis arguments and monotonicity formulae.

## 2 Local Case: Laplacian operator

In this section we assume that s = 1. Therefore system (1.1) turns into the following form

$−Δui=|u|p−1uiin Rn.$(2.1)

Note that this is a particular case of system (4.33). For a general nonlinearity H : ℝm → ℝ such that ∇H ≥ 0, it is proved in [22] that bounded stable solutions of (4.33) are constant when n ≤ 4. For radial solutions, this is known to hold in more dimensions that is when n ≤ 9 without any sign assumptions on the nonlinearity, see [21]. In addition it is known that at least for the case of radial solutions the dimension n = 9 is the optimal dimension. For the rest of this section we prove the following classification of finite Morse index solutions of (2.1).

#### Theorem 2.1

Suppose that u = $\begin{array}{}\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ is a finite Morse index solution of (2.1) when m ≥ 1 and n ≥ 3. Let 1 < p < $\begin{array}{}\frac{n+2}{n-2}\end{array}$ and $\begin{array}{}\frac{n+2}{n-2}\end{array}$ < p < pc(n) where

$pc(n)=∞ if n≤10,(n−2)2−4n+8n−1(n−2)(n−10) if n≥11.$(2.2)

Then each ui must be identically zero. For the Sobolev critical exponent p = $\begin{array}{}\frac{n+2}{n-2}\end{array}$, a solution u has finite energy that is

$∫Rn|u|p+1=∑i=1m∫Rn|∇ui|2<∞.$(2.3)

In this case, if in addition u is stable, then each component ui must be identically zero.

Note that for the case of scalar equations, that is when m = 1, the above theorem is given by Farina in [20]. He used a Moser iteration type argument for the proof. We refer interested reader to Crandall and Rabinowitz in [13] for a similar approach. To prove the above theorem, we apply a blow-down analysis argument as well as the monotonicity formula given as Theorem 1.1. We now derive a few elliptic estimates.

#### Lemma 2.1

Suppose that u is a stable solution of (1.1). Then, for any R > 1

$∫BR|u|p+1≤CRn−2p+1p−1 and ∑i=1m∫BR|∇ui|2≤CRn−2p+1p−1,$(2.4)

where C is a positive constant that is independent from R.

#### Proof

Test the stability inequality on ϕi = ui ζR where ζR : $\begin{array}{}{C}_{c}^{\mathrm{\infty }}\end{array}$(ℝn) → ℝ and ζR ≡ 1 on BR and ζR ≡ 0 on ℝnB2R with ∥ζRL(B2RBR)CR–1. Now multiply both sides of (1.1) with $\begin{array}{}{u}_{i}{\zeta }_{R}^{2}\end{array}$ and integrate by parts. Equating the inequalities that one gets from this and from stability finishes the proof.□

Applying the Hölder’s inequality, we get the following L2 estimate.

#### Corollary 2.1

Suppose that u is a stable solution of (1.1). Then, for any R > 1

$∫BR|u|2≤CRn−4p−1,$(2.5)

where C is a positive constant that is independent from R.

In this part, we present a classification of stable homogeneous solutions. This is a key point in our proof of Theorem 2.1.

#### Theorem 2.2

Suppose that $\begin{array}{}u=\left({u}_{i}{\right)}_{i=1}^{m}\text{\hspace{0.17em}}\text{\hspace{0.17em}}for\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{i}={r}^{-\frac{2}{p-1}}{\psi }_{i}\left(\theta \right)\end{array}$ is a stable solution of (1.1). Then, each ψi is identically zero provided $\begin{array}{}\frac{n+2}{n-2} where pc(n) is given by (2.2)

#### Proof

We omit the proof here, since a similar argument will be provided in the proof of Theorem 3.2 for the fourth order case. □

#### Proof of Theorem 2.1

The proof is based on a blow-down analysis and it relies on the monotonicity formula provided as Theorem 1.1. We omit the details, since similar arguments will be provided for the poof of Theorem 3.1. For the case of Sobolev critical exponent, one can conclude the result via applying the Pohozaev indentity. □

## 3 Local Case: Bi-Laplacian Operator

In this section, we consider the following fourth order system

$Δ2ui=|u|p−1uiinRn.$(3.1)

This section is devoted to the proof of the following theorem.

#### Theorem 3.1

Suppose that $\begin{array}{}u=\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ is a finite Morse index solution of (3.1) when m ≥ 1 and n ≥ 5. Let 1 < p < $\begin{array}{}\frac{n+4}{n-4}\text{\hspace{0.17em}}\text{\hspace{0.17em}}and\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{n+4}{n-4}\end{array}$ < p < c(n) where

$p¯c(n)= ∞if n≤12,n+2−n2+4−nn2−8n+32n−6−n2+4−nn2−8n+32if n≥13.$(3.2)

Then, each ui must be identically zero. For the Sobolev critical exponent $\begin{array}{}p=\frac{n+4}{n-4}\end{array}$, a solution u has finite energy that is

$∫Rn|u|p+1=∑i=1m∫Rn|Δui|2<∞.$(3.3)

In this case, if in addition u is stable, then each component ui must be identically zero.

In order to prove the above theorem, we are required to establish some a priori estimates on solutions.

#### Lemma 3.1

Suppose that $\begin{array}{}u=\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ be a smooth stable solution of (3.1) and set $\begin{array}{}w=\left({w}_{i}{\right)}_{i=1}^{m}\end{array}$ where wi = Δ ui. Then the following estimate holds,

$∫Rn(|w|2+|u|p+1)ζ2≤C∫Rn|u|2|∇Δζ||∇ζ|+|Δζ|2+|Δ|∇ζ|2|+|u||w||∇ζ|2,$(3.4)

for a test function $\begin{array}{}\zeta :{C}_{c}^{\mathrm{\infty }}\left({\mathbb{R}}^{n}\right)\to \mathbb{R}\end{array}$.

#### Proof

Test the stability inequality on uiζ where ζ is a test function then multiply both sides of the ith equation of (1.1) with u–i ζ2. Equating these completes the proof. □

Applying an appropriate test function yields the following estimate.

#### Corollary 3.1

Under the same assumptions as in (3.1), there exists a constant C such that

$∫BR(x)|w|2+|u|p+1≤CR−4∫B2R(x)∖BR(x)|u|2+CR−2∫B2R(x)∖BR(x)|u||w|,$(3.5)

and therefore

$∫BR(x)|w|2+|u|p+1≤CRn−4p+1p−1.$(3.6)

#### Proof

Set the test function $\begin{array}{}{\zeta }_{R}\in {C}_{c}^{1}\left({\mathbb{R}}^{n}\right)\end{array}$ in (3.1) where 0 ≤ ζR ≤ 1 being the following test function

$ζR(x)=1,if |x|2R,$

satisfying ∥∇ζRR−1 and ∥ΔζRR−2. □

We now provide classification of stable homogeneous solutions. Note first that the following Hardy-Rellich inequality with the best constant holds. Suppose that $\begin{array}{}h:{C}_{c}^{2}\left({\mathbb{R}}^{n}\right)\to \mathbb{R}\end{array}$ then

$∫Rn|Δh|2dx≥n2(n−4)216∫Rnh2|x|4dx.$(3.7)

This inequality implies that the singular solution given by (1.29) is stable if and only if

$p|A|p−1=pα(α+2)(n−α−2)(n−α−4)≤n2(n−4)216,$(3.8)

where $\begin{array}{}\alpha :=\frac{4}{p-1}\end{array}$ and 𝓐 given in (1.29).

#### Theorem 3.2

Suppose that $\begin{array}{}u=\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ for $\begin{array}{}{u}_{i}={r}^{-\frac{4}{p-1}}{\psi }_{i}\left(\theta \right)\end{array}$ is a stable solution of (3.1). Then, each ψi ≡ 0 provided $\begin{array}{}\frac{n+4}{n-4} where p̄c(n) given by (3.2).

#### Proof

It is straightforward to see that $\begin{array}{}\psi =\left({\psi }_{i}{\right)}_{i=1}^{m}\end{array}$ satisfies

$Δθ2ψi−αΔθψi+βψi=|ψ|p−1ψi,$(3.9)

for

$α:=(q+2)(n−4−q)+q(n−2−q)andβ:=q(q+2)(n−4−q)(n−2−q),$(3.10)

where $\begin{array}{}q:=\frac{4}{p-1}\end{array}$. Multiplying both sides of (3.9) with ψi and integrating over 𝓢n−1, we conclude

$∑i=1m∫Sn−1|Δψi|2+α|∇ψi|2+β∫Sn−1|ψ|2=∫Sn−1|ψ|p+1.$(3.11)

We now test the stability inequality (1.30) on $\begin{array}{}\varphi =\left({\varphi }_{i}{\right)}_{i=1}^{m}\end{array}$ for $\begin{array}{}{\varphi }_{i}:={r}^{-\frac{n-4}{2}}{\psi }_{i}\left(\theta \right){\eta }_{ϵ}\left(r\right)\end{array}$. Here, ηϵ is a standard cut-off function $\begin{array}{}{\eta }_{ϵ}\in {C}_{c}^{1}\left({\mathbb{R}}_{+}\right)\end{array}$ at the origin and at infinity that is ηϵ = 1 for ϵ < r < ϵ−1 and ηϵ = 0 for either r < ϵ/2 or r > 2/ϵ. Applying similar ideas provided in [24], we get

$p∫Sn−1|ψ|p+1≤∑i=1m∫Sn−1|Δψi|2+n(n−4)2|∇ψi|2+n2(n−4)216∫Sn−1|ψ|2.$(3.12)

Combining (3.11) and (3.12) we get

$(p−1)∑i=1m∫Sn−1|Δψi|2+pα−n(n−4)2|∇ψi|2+pβ−n2(n−4)216|ψ|2≤0.$(3.13)

Note the coefficients p − 1, $\begin{array}{}p\alpha -\frac{n\left(n-4\right)}{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}p\beta -\frac{{n}^{2}\left(n-4{\right)}^{2}}{16}\end{array}$ are positive when $\begin{array}{}\frac{n+4}{n-4} where pc(n) is given by (3.2). □

#### Proof of Theorem 3.1

The proof of the Sobolev critical case relies on applying the Pohozaev identity and we omit it here. We now provide a sketch of the proof when $\begin{array}{}p>\frac{n+4}{n-4}\end{array}$ in a few steps.

• Step 1

limr→∞ E(u, r, 0) < ∞. Note that the energy E(u, r, 0) is nondecreasing in r, as given in Theorem 4.14. Therefore,

$E(u,r,0)≤r−1∫r2rE(u,t,0)dt≤r−2∫r2r∫tt+rE(u,λ,0)dλdt.$(3.14)

Applying estimates given in Corollary 3.1 imply that the right-hand side of (3.14) is bounded.

• Step 2

Define $\begin{array}{}{u}_{i}^{\lambda }\left(x\right)={\lambda }^{\frac{4}{p-1}}{u}_{i}\left(\lambda x\right)\end{array}$ for each 1 ≤ im where $\begin{array}{}u=\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ is a stable solution of (1.1). Then $\begin{array}{}{u}_{i}^{\lambda }\to {u}_{i}^{\mathrm{\infty }}\end{array}$ where $\begin{array}{}{u}_{i}^{\mathrm{\infty }}\in {W}_{loc}^{2,2}\left({\mathbb{R}}^{n}\right)\cap {L}_{loc}^{p+1}\left({\mathbb{R}}^{n}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}^{\mathrm{\infty }}=\left({u}_{i}^{\mathrm{\infty }}{\right)}_{i=1}^{m}\end{array}$ is a stable solution of (1.1). Set $\begin{array}{}{w}_{i}^{\lambda }\left(x\right):={\lambda }^{\frac{4}{p-1}+2}{w}_{i}\left(\lambda x\right)\end{array}$. From Corollary 3.1, we have

$∫Br(x)|wλ|2+|uλ|p+1≤Crn−4p+1p−1.$(3.15)

From elliptic estimates, up to a subsequence, $\begin{array}{}{u}_{i}^{\lambda }\to {u}_{i}^{\mathrm{\infty }}\end{array}$ for each 1 ≤ im weakly in $\begin{array}{}{W}_{loc}^{2,2}\left({\mathbb{R}}^{n}\right)\cap {L}_{loc}^{p+1}\left({\mathbb{R}}^{n}\right)\end{array}$. From compactness embeddings and applying interpolation we arrive at $\begin{array}{}{u}_{i}^{\lambda }\to {u}_{i}^{\mathrm{\infty }}\end{array}$ in $\begin{array}{}{L}_{loc}^{q}\left({\mathbb{R}}^{n}\right)\end{array}$ for any q ∈ [1, p + 1). Note also that u is a stable solution, since uλ is a stable solution and we can send λ to infinity.

• Step 3

u is a homogeneous solution. This is a direct consequence of the monotonicity formula and the following fact

$limλ→∞[E(uλ,R,0)−E(uλ,r,0)]=0.$(3.16)

The left-hand side of (3.16) is bounded from below by

$E(uλ,R,0)−E(uλ,r,0)≥C(n,p)∑i=1m∫BR∖Br4p−1|x|−1uiλ+∂uiλ∂r2|x|8p−1+2−ndx$(3.17)

$=C(n,p)∑i=1m∫BR∖Br4p−1|x|−1ui∞+∂ui∞∂r2|x|8p−1+2−ndx.$(3.18)

This implies that for each 1 ≤ im, we have

$ui∞(x)=|x|−4p−1ui∞x|x|.$(3.19)

This completes the proof of this step.

• Step 4

limr→∞ E(u, r, 0) = 0. Since each $\begin{array}{}{u}_{i}^{\mathrm{\infty }}\end{array}$ is a homogenous function, Theorem 3.2 implies that u = 0. Therefore, $\begin{array}{}\underset{\lambda \to \mathrm{\infty }}{lim}{u}_{i}^{\lambda }=0\end{array}$ strongly in L2(B4) for each i = 1, ⋯, m that is

$limλ→∞∫B4|uiλ|2=0and limλ→∞∫B4|uiλwiλ|=0.$(3.20)

Applying Corollary 3.1, we conclude

$limλ→∞∫B4|wλ|2+|uλ|p+1=0.$(3.21)

On the other hand, there exists r0 > 0 such that

$limi→∞||uλi||W2,2(∂Br0)=0.$(3.22)

This implies that

$limi→∞E(u,λir0,0)=limi→∞E(uλi,r0,0).$(3.23)

The fact that E is nondecreasing, that is given as a monotonicity formula in Theorem 1.2, completes the proof. □

## 4 Nonlocal Case: Lower and higher order fractional Laplacian

In this section, we consider system (1.1) with the fractional Laplacian operator (−Δ)s where s ∈ (0, 2) for s ≠ 1 and we establish Theorem 1.5. We first note that the following Hardy inequality holds for n > 2s

$∫Rn|ξ|2s|h^|2dξ>Λn,s∫Rn|x|−2sh2dx,$(4.1)

for any $\begin{array}{}h\in {C}_{c}^{\mathrm{\infty }}\left({\mathbb{R}}^{n}\right)\end{array}$ where the optimal constant is given by

$Λn,s=22sΓ(n+2s4)2Γ(n−2s4)2.$(4.2)

Note that $\begin{array}{}{\mathit{\Lambda }}_{n,2}=\frac{{n}^{2}\left(n-4{\right)}^{2}}{16}\end{array}$ meaning (4.1) recovers (3.7) for s = 2. For more information interested readers are encouraged to see [29] by Herbst (and also [46]). We now provide a classification result for homogeneous solutions. Note that for the case of scalar equations, that is when m = 1, this classification was given in [24]. Note also that the proof is valid regardless of magnitude of the parameter s.

#### Theorem 4.1

Suppose $\begin{array}{}{u}_{i}={r}^{-\frac{2s}{p-1}}{\psi }_{i}\left(\theta \right)\end{array}$ is a stable solution of (1.1) for s ∈ (0, 2) and s ≠ 1. Then, each ψi vanishes identically, provided $\begin{array}{}p>\frac{n+2s}{n-2s}\end{array}$ and

$pΓ(n2−sp−1)Γ(s+sp−1)Γ(sp−1)Γ(n−2s2−sp−1)>Γ(n+2s4)2Γ(n−2s4)2.$(4.3)

#### Proof

Since $\begin{array}{}u=\left({u}_{i}{\right)}_{i=1}^{m}\end{array}$ is a solution of (1.1), each ψi satisfies

$|x|−2psp−1|ψ|p−1(θ)ψi(θ)=∫|x|−2sp−1ψi(θ)−|y|−2sp−1ψi(σ)|x−y|n+2sdy=|x|−2psp−1[∫ψi(θ)−t−2sp−1ψi(θ)(t2+1−2t<θ,σ>)n+2s2tn−1dtdσ+∫t−2sp−1(ψi(θ)−ψi(σ)(t2+1−2t<θ,σ>)n+2s2tn−1dtdσ],$(4.4)

where we have used the change of variable |y| = t|x|. Simplifying the above, for each i, we obtain

$ψi(θ)An,s+∫Sn−1K2sp−1(<θ,σ>)(ψi(θ)−ψi(σ))dσ=|ψ|p−1(θ)ψi(θ),$(4.5)

where

$An,s:=∫0∞∫Sn−11−t−2sp−1(t2+1−2t<θ,σ>)n+2s2tn−1dσdt,$(4.6)

and

$K2sp−1(<θ,σ>):=∫0∞tn−1−2sp−1(t2+1−2t<θ,σ>)n+2s2dt.$(4.7)

Multiplying (4.5) with ψi and integrating we get

$∫Sn−1|ψ|2(θ)An,s+∫Sn−1K2sp−1(<θ,σ>)|ψ(θ)−ψ(σ)|2dθdσ=∫Sn−1|ψ|p+1(θ)dθ.$(4.8)

We now test the stability inequality (1.30) for $\begin{array}{}{\varphi }_{i}\left(x\right)={r}^{-\frac{n-2s}{2}}{\psi }_{i}\left(\theta \right){\eta }_{ϵ}\left(r\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{i}={r}^{-\frac{2s}{p-1}}{\psi }_{i}\left(\theta \right)\end{array}$ with the same ηϵ(r) as the one given in the proof of Theorem 3.2. Applying similar ideas provided in [24], we conclude

$Λn,s∫Sn−1|ψ|2+∫Sn−1Kn−2s2(<θ,σ>)|ψ(θ)−ψ(σ)|2dσ≥p∫Sn−1|ψ|p+1,$(4.9)

when

$Λn,s:=∫0∞∫Sn−11−tn−2s2(t2+1−2t<θ,σ>)n+2s2tn−1dσdt.$(4.10)

Combining (4.9) and (4.8), we end up with

$(Λn,s−pAn,s)∫Sn−1|ψ|2+∫Sn−1(Kn−2s2−pK2sp−1)(<θ,σ>)|ψ(θ)−ψ(σ)|2dσ≥0.$(4.11)

The fact that Kα is decreasing in α implies $\begin{array}{}{K}_{\frac{n-2s}{2}}<{K}_{\frac{2s}{p-1}}\end{array}$ for $\begin{array}{}p>\frac{n+2s}{n-2s}\end{array}$. Therefore, $\begin{array}{}{K}_{\frac{n-2s}{2}}-p{K}_{\frac{2s}{p-1}}<0\end{array}$. On the other hand the assumption of the theorem implies that Λn,sp An,s < 0. Therefore, each ψi vanishes identically. This completes the proof. □

## 4.1 Lower Order Fractional Laplacian Operator

In this part, we show that Theorem 1.5 holds when 0 < s < 1. To do so we provide the following estimate first.

#### Lemma 4.1

Suppose that $\begin{array}{}p\ne \frac{n+2s}{n-2s}\end{array}$. Let u be a solution of (1.1) that is stable outside a ball BR0 and v satisfies (1.16). Then there exists a constant C > 0 such that

$∫BRy1−2s|v|2≤CRn+2−2sp+1p−1,$(4.12)

for any R > 3R0.

#### Lemma 4.2

Let u be a solution of (1.1) that is stable outside a ball BR0 and v satisfies (1.16). Then there exists a positive constant C such that

$∫BR∩∂R+n+1|v|p+1+∑i=1m∫BR∩R+n+1y1−2s|∇vi|2≤CRn−2sp+1p−1.$(4.13)

#### Proof of Theorem 1.5

when 0 < s < 1. We omit the proof here since arguments are very similar to the ones which will be provided for the case of 1 < s < 2. □

## 4.2 Higher Order Fractional Laplacian Operator

As the last past of this section, we shall restrict ourselves to the case 1 < s < 2. Let us start with the following integral estimate on stable solutions.

#### Lemma 4.3

Let u be a solution of (1.1) that is stable outside a ball BR0 and v satisfies (1.19). Then there exists a positive constant C such that

$∫∂R+n+1|v|p+1η2+∑i=1m∫R+n+1yb|Δbvi|2η2≤C∫R+n+1yb|v|2|Δbη|2+|Δb|∇η|2|+|∇η⋅∇Δbη|$(4.14)

$+C∑i=1m∫R+n+1yb|vi||Δbvi||∇η|2.$(4.15)

#### Proof

Here is a sketch of the proof. Multiply the ith equation of (1.19) with yb uiη2 where η is a test function to get

$Cn,s∫∂R+n+1|v|p+1η2=∑i=1m∫R+n+1ybΔbviΔb(viη2).$(4.16)

Note also that for each i these identities hold for a test function η

$ΔbviΔb(viη2)−|Δb(viη)|2=−vi2|Δbη|2+2viΔbvi|∇η|2−4|∇vi⋅∇η|2−4viΔbη∇vi⋅∇η,$(4.17)

$Δb(viη)=ηΔbvi+viΔbη+2∇vi⋅∇η.$(4.18)

Applying these identifies together with (4.16) one can see that

$Cn,s∫∂R+n+1|v|p+1η2=∑i=1m∫R+n+1yb|Δb(viη)|2+2∑i=1m∫R+n+1ybviΔbvi|∇η|2−4∑i=1m∫R+n+1yb|∇vi⋅∇η|2+∫R+n+1yb|v|2(|Δbη|2+2∇η⋅∇Δbη).$(4.19)

Testing the stability inequality (1.30) on ϕi = uiη and applying (4.19) one can finish the proof. □

One can set η to be the standard test function to get the following estimate.

#### Corollary 4.1

With the same assumption as Lemma 4.3. Then there exists a positive constant C such that

$∫BR∩∂R+n+1|v|p+1+∑i=1m∫BR∩R+n+1yb|Δbvi|2≤CR−4∫BR∩R+n+1yb|v|2.$(4.20)

Here we provide more decay estimates of solutions. These Lemmata are main tools in our proof of Theorem 1.5. We omit the proofs here and we refer interested readers to see the proof of Lemma 4.5-4.6 in [24] where similar arguments are applied.

#### Lemma 4.4

Suppose that $\begin{array}{}p\ne \frac{n+2s}{n-2s}\end{array}$. Let u be a solution of (1.1) that is stable outside a ball BR0 and v satisfies (1.19). Then there exists a constant C > 0 such that

$∫BRyb|v|2≤CRn+4−2sp+1p−1,$(4.21)

for any R > 3R0.

#### Lemma 4.5

Let u be a solution of (1.1) that is stable outside a ball BR0 and v satisfies (1.19). Then there exists a positive constant C such that

$∫BR∩∂R+n+1|v|p+1+∑i=1m∫BR∩R+n+1yb|Δbvi|2≤CRn−2sp+1p−1.$(4.22)

#### Proof of Theorem 1.5

when 1 < s < 2. The proof of the case p = pS(n, s) is based on the Pohozaev identity provided in [24, 38], and we omit it here. The proof is based on the monotonicity formula that is Theorem 1.4 and a blow-down analysis.

• Step 1

The energy is bounded, that is limλ→∞ E(v, 0, λ) < ∞. This is a direct consequence of the monotonicity formula and Lemma 4.4 and Lemma 4.5. Similar to the proof of Theorem 3.1 we have

$E(v,λ,0)≤λ−2∫λ2λ∫tt+λE(v,y,0)dydt.$(4.23)

Lemma 4.5 and Lemma 4.4 imply that the right-hand side of (4.23) is bounded.

• Step 2

The sequence $\begin{array}{}{v}_{i}^{\lambda }\end{array}$ converges weakly in $\begin{array}{}{H}_{loc}^{1}\left({\mathbb{R}}^{n},{y}^{3-2s}dxdy\right)\end{array}$ to a function $\begin{array}{}{v}_{i}^{\mathrm{\infty }}\end{array}$ where each $\begin{array}{}{v}_{i}^{\mathrm{\infty }}\end{array}$ is homogeneous for 1 ≤ im and therefore they are zero. Note that the convergence part is a direct consequence of the elliptic estimates. We now show that each $\begin{array}{}{v}_{i}^{\mathrm{\infty }}\end{array}$ is homogeneous. From the boundedness of the energy we have

$limk→∞E(v,R2λk,0)−E(v,R1λk,0)=0.$(4.24)

From this and applying the scaling invariance of the energy and also the monotonicity formula we get

$0=lim infk→∞∑i=1m∫(BR2∖BR1)∩R+n+1y3−2sr4sp−1+2s−2−n2sp−1r−1viλk+∂viλk∂r2dydx$(4.25)

$≥∑i=1m∫(BR2∖BR1)∩R+n+1y3−2sr4sp−1+2s−2−n2sp−1r−1vi∞+∂vi∞∂r2dydx,$(4.26)

since we have the weak convergence of $\begin{array}{}\left({v}_{i}^{{\lambda }_{k}}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{to}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{v}_{i}^{\mathrm{\infty }}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{H}_{loc}^{1}\left({\mathbb{R}}^{n},{y}^{3-2s}dydx\right)\end{array}$. Therefore,

$2sp−1r−1vi∞+∂vi∞∂r=0a.e. inR+n+1,$(4.27)

for each 1 ≤ im.

• Step 3

limλ→∞ E(v, λ, 0) = 0. Note that the monotonicity formula implies that

$E(v,λ,0)≤λ−1∫λ2λE(t)dt≤sup[λ,2λ]I+Cλ−n−1+2s(p+1)p−1∫B2λ∖Bλ|v|2,$(4.28)

where

$I(v,λ)=I(vλ,1)=12∑i=1m∫R+n+1∩B1y3−2s|Δbviλ|2dxdy−κsp+1∫∂R+n+1∩B1|vλ|p+1dx.$(4.29)

Note that lim λ→∞ I(v, λ) = 0. On the other hand, from the fact that u is a homogenous solution we have

$limλ→∞uiλ=0,$(4.30)

strongly in L2(B4). Therefore,

$limλ→∞∫B4|uλ|2=0.$(4.31)

This implies that limλ→∞ E(v, λ, 0) = 0 and completes the proof. □

We end this section with the following open problem in regards to monotonicity formulae for the Lane-Eden system.

#### Open Problem 1

Consider the Lane-Emden system for any parameters s > 0 and 1 < q < p

$(−Δ)su=vpinRn,(−Δ)sv=uqinRn.$(4.32)

Proving a monotonicity formula for the above system, similar to the ones given in Theorem 1.1-1.4, seems more challenging to derive. Note that for the case of 1 < q = p and 1 = q < p such monotonicity formulae are known, see the introduction. Needless to mention that the Lane-Emden system (4.32) is not a gradient system, meaning that it is not of the following form

$(−Δ)su=∇H(u)inRn.$(4.33)

where u : ℝn → ℝm.

Lastly, let us mention that in [23] we apply monotonicity formulae derived in this article to prove regularity of free boundaries and partial regularity of weak solutions for certain coupled elliptic systems.

## Acknowledgement

The first author appreciates H. Shahgholian and J. Wei for their hospitality during his visits to KTH and UBC, and he is thankful to E. Hebey for his talk at the conference on the occasion of Michael Struwe’s 60th birthday at ETH.

The first author is partially supported by University of Texas at San Antonio Start-up Grant. The second author was partially supported by Swedish Research Council.

## References

• [1]

H.W. Alt, L.A. Caffarelli and A. Friedman, Variational problems with two phases and their free bound-aries, Trans. Amer. Math. Soc., 282(2) (1984) pp. 431- 461.

• [2]

W. K. Allard. On the first variation of a varifold. Ann. of Math., 95 (1972) pp. 417-491.

• [3]

J. Andersson, Optimal regularity for the Signorini problem and its free boundary, Inventiones Math., 204 (2016) pp. 1-82.

• [4]

J. Andersson, H. Shahgholian, N. N. Uraltseva. and G. S. Weiss, Equilibrium points of a singular cooperative system with free boundary, Advances in Mathematics, 280 (2015) pp. 743-771.

• [5]

B. Ou, Global solutions to a free boundary problem. Comm. Partial Differential Equations, 19 (1994) pp. 369-397.

• [6]

L. Caffarelli, B. Gidas and J. Spruck, Asymptotic symmetry and local behavior of semilinear elliptic equations with critical Sobolev growth, Comm. Pure Appl. Math., 42 (1989) pp. 271-297

• [7]

L. Caffarelli, H. Shahgholian and K. Yeressian, A minimization problem with free boundary related to a cooperative system, Duke Math. J., 167 (2018) pp. 1825-1882.

• [8]

L. Caffarelli and L. Silvestre, An extension problem related to the fractional Laplacian, Comm. Partial Differential Equations, 32 (2007) pp. 1245-1260.

• [9]

J. Case and S. A. Chang, On fractional GJMS operators, Comm. Pure Appl. Math., 69 (2016) pp. 1017-1061.

• [10]

W. X. Chen and C. Li, Classification of solutions of some nonlinear elliptic equations, Duke Math. J., 63 (1991) pp. 615-622.Google Scholar

• [11]

W. Chen and C. Li, B. Ou, Classification of solutions for an integral equation, Comm. Pure Appl. Math., 59 (2006) pp. 330-343.

• [12]

Sun-Yung A. Chang and M. Gonzalez, Fractional Laplacian in conformal geometry, Advances in Mathematics, 226 (2011) pp. 1410-1432.

• [13]

M.G. Crandall and P.H. Rabinowitz, Some continuation and variation methods for positive solutions of nonlinear elliptic eigenvalue problems, Arch. Rat. Mech. Anal., 58 (1975) pp. 207-218.

• [14]

J. Davila, L. Dupaigne and J. Wei, On the fractional Lane-Emden equation, Trans. Amer. Math. Soc., 369 (2017) pp. 6087-6104.

• [15]

J. Davila, L. Dupaigne, K. Wang and J. Wei, A monotonicity formula and a Liouville-type theorem for a fourth order supercritical problem, Advances in Mathematics, 258 (2014) pp. 240-285.

• [16]

O. Druet and E. Hebey, Stability for strongly coupled critical elliptic systems in a fully inhomogeneous medium, Analysis and PDE, (2009) pp. 305-359. Google Scholar

• [17]

O. Druet, E. Hebey and J. Vétois, Bounded stability for strongly coupled critical elliptic systems below the geometric threshold of the conformal Laplacian, Journal of Functional Analysis, 258 (2010) pp. 999-1059.

• [18]

E. C. Evans, Partial regularity for stationary harmonic maps into spheres. Arch. Rational Mech. Anal., 116 (1991) pp. 101-113.

• [19]

E. C. Evans, Monotonicity formulae for variational problems, Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 371 (2013), no. 2005, 9 pp. Google Scholar

• [20]

A. Farina; On the classification of solutions of the Lane-Emden equation on unbounded domains ofn, J. Math. Pures Appl., (9) 87 (2007) pp. 537-561.

• [21]

M. Fazly, Rigidity results for stable solutions of symmetric systems, Proc. Amer. Math. Soc., 143 (2015) pp. 5307-5321.

• [22]

M. Fazly and N. Ghoussoub, De Giorgi type results for elliptic systems, Calc. Var. Partial Differential Equations, 47 (2013) pp. 809-823.

• [23]

M. Fazly and H. Shahgholian, Regularity of free boundaries and partial regularity of weak solutions for coupled elliptic systems, In Progress. Google Scholar

• [24]

M. Fazly and J. Wei, On finite Morse index solutions of higher order fractional Lane-Emden equations, American Journal of Mathematics, 139 (2017) pp. 433-460.

• [25]

W. H. Fleming, On the oriented Plateau problem. Rend. Circ. Mat. Palermo, (2) 11 (1962) pp. 69-90 Google Scholar

• [26]

B. Gidas and J. Spruck, A priori bounds for positive solutions of nonlinear elliptic equations, Comm. Partial Differential Equations, 6 (1981) pp. 883-901.

• [27]

B. Gidas, W. M. Ni and L. Nirenberg, Symmetry of positive solutions of nonlinear elliptic equations inn, in Mathematical Analysis and Applications, Part A, Adv. in Math. Suppl. Stud., 7a, Academic Press, New York-London, 1981, 369-402.Google Scholar

• [28]

Y. Giga and R. V. Kohn, Asymptotically self-similar blow-up of semilinear heat equations, Comm. Pure Appl. Math., 38 (1985) pp. 297-319.

• [29]

E. Hebey, Critical elliptic systems in potential form, Adv. Differential Equations, 11 (2006) pp. 511-600. Google Scholar

• [30]

E. Hebey and P. D. Thizy, Stationary Kirchhoff systems in closed high dimensional manifolds, Commun. Contemp. Math., 18, 1550028 (2016) 53 pp.

• [31]

D. D. Joseph and T. S. Lundgren, Quasilinear Dirichlet problems driven by positive sources, Arch. Rational Mech. Anal., 49 (1972/73) pp. 241- 269.

• [32]

Y. Li, Remark on some conformally invariant integral equations: the method of moving spheres, J. Eur. Math. Soc. (JEMS), 6 (2004) pp. 153-180. Google Scholar

• [33]

C. S. Lin, A classification of solutions of a conformally invariant fourth order equation inN, Comment. Math. Helv., 73 (1998) pp. 206-231.

• [34]

B. Noris, H. Tavares, S. Terracini and G. Verzini, Uniform Hölder bounds for nonlinear Schrödinger systems with strong competition, Comm. Pure Appl. Math., 63 (2010) pp. 267-302.

• [35]

F. Pacard, Partial regularity for weak solutions of a nonlinear elliptic equation, manuscripta mathematica, 79 (1993) pp. 161-172.

• [36]

A. Petrosyan, H. Shahgholian. and N. Uraltseva, Regularity of free boundaries in obstacle-type problems, Graduate Studies in Mathematics, 136. American Mathematical Society, Providence, RI, 2012. x+221. Google Scholar

• [37]

P. Price. A monotonicity formula for Yang-Mills fields, Manuscripta Math., 43 (1983) pp. 131-166.

• [38]

X. Ros-Oton and J. Serra, Local integration by parts and Pohozaev identities for higher order fractional Laplacians, Discrete Contin. Dyn. Syst. A, 35 (2015) pp. 2131-2150.

• [39]

R. Schoen, Analytic aspect of the harmonic map problem, Seminar on nonlinear partial differential equations (Berkeley, Calif., 1983), 321-358, Math. Sci. Res. Inst. Publ., 2, Springer, New York, 1984. Google Scholar

• [40]

R. Schoen and K. Uhlenbeck, A regularity theory for harmonic maps, J. Differential Geometry, 17 (1982) pp. 307-335.

• [41]

S. Terracini, G. Verzini and A. Zilio, Uniform Hölder regularity with small exponent in competition-fractional diffusion systems, Discrete Contin. Dyn. Syst., 34 (2014) pp. 2669-2691. Google Scholar

• [42]

J. Wei and T. Weth, Radial solutions and phase separation in a system of two coupled schrödinger equations, Arch. Rational Mech. Anal. 190 (2008) pp. 83-106.

• [43]

J. Wei and X. Xu; Classification of solutions of higher order conformally invariant equations, Math. Ann., 313 (1999) pp. 207-228.

• [44]

G. S. Weiss, A homogeneity improvement approach to the obstacle problem, Invent. Math., 138 (1999) pp. 23-50.

• [45]

G. S. Weiss, Partial regularity for weak solutions of an elliptic free boundary problem, Communications in Partial Differential Equations 23 (1998) pp. 439-455.

• [46]

D. Yafaev, Sharp constants in the Hardy-Rellich inequalities, Journal of Functional Analysis, 168 (1999) pp. 121-144.

• [47]

R. Yang, On higher order extensions for the fractional Laplacian, preprint. http://arxiv.org/pdf/1302.4413v1.pdf

## Footnotes

• 1

It is worth mentioning that a different type of almost monotonicity functional was used by Arne Beurling in his thesis: A. Beurling, Etudes sur un probleme de majoration, thesis, (Uppsala), 1933. The monotonicity functional of A. Beurling states that for a jordan curve y, the product of the harmonic measures for both sides of the curve can be controlled as follows: ω1(Br(z) ∩ y) ⋅ ω2(Br(z) ∩ y) ≤ Ar2, where zy, and ω1, ω2 denote the harmonic measures on each side of the curve.

Accepted: 2019-01-06

Published Online: 2019-06-06

Published in Print: 2019-03-01

Citation Information: Advances in Nonlinear Analysis, Volume 9, Issue 1, Pages 479–495, ISSN (Online) 2191-950X, ISSN (Print) 2191-9496,

Export Citation