Regularity of stable solutions of a Lane-Emden type system

Similar documents
Liouville theorems for stable Lane-Emden systems and biharmonic problems

The critical dimension for a fourth order elliptic problem with singular nonlinearity

Perturbations of singular solutions to Gelfand s problem p.1

Heat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control

Extremal Solutions and Instantaneous Complete Blow-up for Elliptic and Parabolic Problems

On the Hénon-Lane-Emden conjecture

Non-homogeneous semilinear elliptic equations involving critical Sobolev exponent

Minimization problems on the Hardy-Sobolev inequality

Non-radial solutions to a bi-harmonic equation with negative exponent

New exact multiplicity results with an application to a population model

SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction

MULTIPLE SOLUTIONS FOR CRITICAL ELLIPTIC PROBLEMS WITH FRACTIONAL LAPLACIAN

Global unbounded solutions of the Fujita equation in the intermediate range

Stability properties of positive solutions to partial differential equations with delay

A simplified proof of a conjecture for the perturbed Gelfand equation from combustion theory

Note on the Chen-Lin Result with the Li-Zhang Method

A DEGREE THEORY FRAMEWORK FOR SEMILINEAR ELLIPTIC SYSTEMS

Saddle Solutions of the Balanced Bistable Diffusion Equation

Integro-differential equations: Regularity theory and Pohozaev identities

Ciudad Universitaria C.P., México D.F., México b Departamento de Ingeniería Matemática and CMM, Universidad de Chile, Santiago, Chile

Hardy Rellich inequalities with boundary remainder terms and applications

ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.

LIFE SPAN OF BLOW-UP SOLUTIONS FOR HIGHER-ORDER SEMILINEAR PARABOLIC EQUATIONS

Revisit the biharmonic equation with negative exponent in lower dimensions

EXISTENCE OF NONTRIVIAL SOLUTIONS FOR A QUASILINEAR SCHRÖDINGER EQUATIONS WITH SIGN-CHANGING POTENTIAL

A LIOUVILLE THEOREM FOR p-harmonic

Symmetry and monotonicity of least energy solutions

Liouville Theorems for Integral Systems Related to Fractional Lane-Emden Systems in R N +

HARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx.

Rigidity Results for Elliptic PDEs

Kato s inequality when u is a measure. L inégalité de Kato lorsque u est une mesure

ON NONHOMOGENEOUS BIHARMONIC EQUATIONS INVOLVING CRITICAL SOBOLEV EXPONENT

Nonlinear elliptic systems with exponential nonlinearities

EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction

Presenter: Noriyoshi Fukaya

UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS

Applications of the compensated compactness method on hyperbolic conservation systems

INTERNATIONAL PUBLICATIONS (USA)

Hardy inequalities, heat kernels and wave propagation

EXISTENCE AND REGULARITY OF WEAK SOLUTIONS FOR SINGULAR ELLIPTIC PROBLEMS. 1. Introduction In this article we study the quasilinear elliptic problem

NONLINEAR SCHRÖDINGER ELLIPTIC SYSTEMS INVOLVING EXPONENTIAL CRITICAL GROWTH IN R Introduction

Blow-up for a Nonlocal Nonlinear Diffusion Equation with Source

ON A CONJECTURE OF P. PUCCI AND J. SERRIN

SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction

A NOTE ON ZERO SETS OF FRACTIONAL SOBOLEV FUNCTIONS WITH NEGATIVE POWER OF INTEGRABILITY. 1. Introduction

UNIQUENESS RESULTS ON SURFACES WITH BOUNDARY

PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION

BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED

Existence of a ground state and blow-up problem for a nonlinear Schrödinger equation with critical growth

EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY

I Results in Mathematics

EXISTENCE AND NONEXISTENCE OF ENTIRE SOLUTIONS TO THE LOGISTIC DIFFERENTIAL EQUATION

Homogenization and error estimates of free boundary velocities in periodic media

Existence of Multiple Positive Solutions of Quasilinear Elliptic Problems in R N

Sharp blow-up criteria for the Davey-Stewartson system in R 3

EXISTENCE OF SOLUTIONS FOR CROSS CRITICAL EXPONENTIAL N-LAPLACIAN SYSTEMS

Some aspects of vanishing properties of solutions to nonlinear elliptic equations

Blow up points of solution curves for a semilinear problem

NONHOMOGENEOUS ELLIPTIC EQUATIONS INVOLVING CRITICAL SOBOLEV EXPONENT AND WEIGHT

AN APPLICATION OF THE LERAY-SCHAUDER DEGREE THEORY TO THE VARIABLE COEFFICIENT SEMILINEAR BIHARMONIC PROBLEM. Q-Heung Choi and Tacksun Jung

MINIMAL SURFACES AND MINIMIZERS OF THE GINZBURG-LANDAU ENERGY

COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS

SIMULTANEOUS AND NON-SIMULTANEOUS BLOW-UP AND UNIFORM BLOW-UP PROFILES FOR REACTION-DIFFUSION SYSTEM

Existence of Positive Solutions to Semilinear Elliptic Systems Involving Concave and Convex Nonlinearities

EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem

Elliptic equations with one-sided critical growth

The De Giorgi-Nash-Moser Estimates

arxiv: v1 [math.ap] 31 May 2013

Journal of Inequalities in Pure and Applied Mathematics

AN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS. To the memory of our friend and colleague Fuensanta Andreu

A Semilinear Elliptic Problem with Neumann Condition on the Boundary

The Existence of Nodal Solutions for the Half-Quasilinear p-laplacian Problems

Nonexistence of solutions for quasilinear elliptic equations with p-growth in the gradient

THE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS

EXISTENCE AND UNIQUENESS OF POSITIVE SOLUTION OF A LOGISTIC EQUATION WITH NONLINEAR GRADIENT TERM

BIHARMONIC WAVE MAPS INTO SPHERES

Liouville-type theorems and decay estimates for solutions to higher order elliptic equations

EXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS

SOME PROPERTIES ON THE CLOSED SUBSETS IN BANACH SPACES

Global bifurcations of concave semipositone problems

OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS

Para el cumpleaños del egregio profesor Ireneo Peral

On uniqueness in the inverse conductivity problem with local data

A global solution curve for a class of free boundary value problems arising in plasma physics

Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition

EQUAZIONI A DERIVATE PARZIALI. STEKLOV EIGENVALUES FOR THE -LAPLACIAN

A non-local critical problem involving the fractional Laplacian operator

Uniqueness of Positive Solutions for a Class of p-laplacian Systems with Multiple Parameters

Publication IV. c 2011 arxiv.org. Reprinted with permission.

arxiv: v1 [math.ap] 28 Aug 2018

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT

Existence of positive solution for system of quasilinear elliptic systems on a bounded domain

Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator

Exact multiplicity of boundary blow-up solutions for a bistable problem

arxiv: v1 [math.ap] 4 Nov 2013

VANISHING-CONCENTRATION-COMPACTNESS ALTERNATIVE FOR THE TRUDINGER-MOSER INEQUALITY IN R N

Xiyou Cheng Zhitao Zhang. 1. Introduction

Estimations universelles pour les solutions d EDP elliptiques non linéaires

Symmetry and non-uniformly elliptic operators

ẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.

Transcription:

Regularity of stable solutions of a Lane-Emden type system Craig Cowan Department of Mathematics University of Manitoba Winnipeg, MB R3T N Craig.Cowan@umanitoba.ca February 5, 015 Abstract We examine the system given by u = λ(v + 1) p Ω v = γ(u + 1) θ Ω, u = v = 0 Ω, where λ, γ are positive parameters and where 1 < p θ and where Ω is a smooth bounded domain in R N. We show the extremal solutions associated with the above system are bounded provided N < 1 + (θ + 1) pθ 1 pθ(p + 1) + θ + 1 pθ(p + 1) θ + 1 pθ(p + 1). θ + 1 In particular this shows that the extremal solutions are bounded for any 1 < p θ provided N 10. 010 Mathematics Subject Classification: 35 J 57 Key words: Extremal solution, Stable solution, Regularity of solutions, Systems. 1

1 Introduction In this work we examine the following system: u = λ(v + 1) p Ω (N) λ,γ v = γ(u + 1) θ Ω, u = v = 0 Ω, where Ω is a bounded domain in R N, λ, γ > 0 are positive parameters and where p, θ > 1. Our interest is in the regularity of the extremal solutions associated with (N) λ,γ. In particular we are interested when the extremal solutions of (N) λ,γ are bounded, since one can then apply elliptic regularity theory to show the extremal solutions are classical solutions. The main approach will be to utilize the stability of the extremal solutions to obtain added regularity. We remark that for this system, when p θ, it is not clear how to utilize the stability in a meaningful way and this work represents a first work in this direction. The nonlinearities we examine naturally fit into the following class: (R): f is smooth, increasing, convex on R with f(0) = 1 and f is f(u) superlinear at (i.e. lim u u = ). 1.1 Second order scalar case For a nonlinearity f of type (R) consider the following second order scalar analog of the above system given by { u = λf(u) in Ω (Q) λ u = 0 on Ω. This scalar equation is now quite well understood whenever Ω is a bounded smooth domain in R N. See, for instance, [1,, 3, 4, 14, 15, 17, 18, 19, 1]. We now list the properties one comes to expect when studying (Q) λ. There exists a finite positive critical parameter λ, called the extremal parameter, such that for all 0 < λ < λ there exists a a smooth minimal solution u λ of (Q) λ. By minimal solution, we mean here that if v is another solution of (Q) λ then v u λ a.e. in Ω. For each 0 < λ < λ the minimal solution u λ is semi-stable in the sense that λf (u λ )ψ dx Ω Ω ψ dx, ψ H 1 0 (Ω),

and is unique among all the weak semi-stable solutions. The map λ u λ (x) is increasing on (0, λ ) for each x Ω. This allows one to define u (x) := lim λ λ u λ (x), the so-called extremal solution, which can be shown to be a weak solution of (Q) λ. In addition one can show that u is the unique weak solution of (Q) λ. See [18]. There are no solutions of (Q) λ (even in a very weak sense) for λ > λ. A question which has attracted a lot of attention is whether the extremal function u is a classical solution of (Q) λ. This is of interest since one can then apply the results from [11] to start a second branch of solutions emanating from (λ, u ). The answer typically depends on the nonlinearity f, the dimension N and the geometry of the domain Ω. We now list some known results. [11] Suppose f(u) = e u. If N < 10 then u is bounded. For N 10 and Ω the unit ball u (x) = log( x ). [4] Suppose f satisfies (R) but without the convexity assumption and Ω is the unit ball. Then u is bounded for N < 10. In view of the above result this is optimal. On general domains, and if f satisfies (R), then u is bounded for N 3 [1]. Recently this has been improved to N 4 provided the domain is convex (again one can drop the convexity assumption on f), see [3]. We now examine the generalization of (N) λ,γ given by (P ) λ,γ u = λf(v) Ω v = γg(u) Ω, u = v = 0 Ω, where f and g satisfy (R). Define Q = {(λ, γ) : λ, γ > 0}, U := {(λ, γ) Q : there exists a smooth solution (u, v) of (P ) λ,γ }, and set Υ := U Q. Note that Υ is the analog of λ for the above system. A generalization of (P ) λ,γ was examined in [0] and many results were obtained, including 3

Theorem. (Montenegro, [0]) Suppose f and g satisfy (R). Then 1. U is nonempty.. For all (λ, γ) U there exists a smooth, minimal solution of (P ) λ,γ. 3. For each 0 < σ < there is some 0 < λ σ < such that U {(λ, σλ) : 0 < λ} is given by {(λ, σλ) : 0 < λ < λ σ} H where H is either the empty set or {(λ σ, σλ σ)}. The map σ λ σ is bounded on compact subsets of (0, ). Fix 0 < σ < and let (u λ, v λ ) denote the smooth minimal solution of (P ) λ,σλ for 0 < λ < λ σ. Then u λ (x), v λ (x) are increasing in λ and hence u (x) := lim u λ (x), λ λ σ v (x) := lim u λ (x), λ λ σ are well defined and can be shown to be a weak solution of (P ) λ σ,σλ σ. Our notation will vary slightly from above. Let (λ, γ ) Υ and set σ := γ λ. Define Γ σ := {(λ, σλ) : λ < λ < λ } and we let (u, v ), called the extremal solution associated with (P ) λ,γ, be the pointwise limit of the minimal solutions along the ray Γ σ as λ λ. As mentioned above (u, v ) is a weak solution of (N) λ,γ in a suitable sense. The following result shows that the minimal solutions are stable in some suitable sense and this will be crucial in obtaining regularity of the extremal solutions associated with (N) λ,γ. Theorem. (Montenegro [0]) Let (λ, γ) U and let (u, v) denote the minimal solution of (P ) λ,γ. Then (u, v) is semi-stable in the sense that there is some smooth 0 < ζ, χ H0 1 (Ω) and 0 η such that ζ = λf (v)χ + ηζ, χ = γg (u)ζ + ηχ, in Ω. (1) We give an alternate proof of a result which is slightly different than the above one, but which is sufficient for our purposes. Fix (λ, γ ) Υ, σ := γ λ and let (u λ, v λ ) denote minimal solution of (P ) λ,γ on the ray Γ σ. Taking a derivative in λ of (P ) λ,σλ shows that ζ = λf (v λ ) χ + f(v λ ), χ = λσg (u λ ) ζ + σg(u λ ) in Ω, where ζ := λ u λ and χ := λ v λ. Using the monotonicity of u λ, v λ and the maximum principle shows that ζ, χ > 0. We now recall some known results regarding the regularity of the extremal solutions associated with various systems. In what follows Ω is a bounded domain in R N. 4

In [7] the following system u = λe v Ω (E) λ,γ v = γe u Ω, u = v = 0 Ω, was examined. It was shown that if 3 N 9 and N 8 < γ λ < 8 N, then the extremal solution (u, v ) is bounded. Note that not only does the dimension N play a role but how close (λ, γ ) are to the diagonal γ = λ plays a role. When γ = λ one can show that the above system reduces to the scalar equation u = λe u. We remark that we were unable to extend the methods used in [7] to handle (N) λ,γ except in the case where p = θ. In [9] the system (P ) λ,γ u = λf (u, v), v = γg(u, v) in Ω, with u = v = 0 on Ω was examined examined in the cases where F (u, v) = f (u)g(v), G(u, v) = f(u)g (v) (resp. F (u, v) = f(u)g (v), G(u, v) = f (u)g(v)) and were denoted by (G) λ,γ (resp. (H) λ,γ ). It was shown that the extremal solutions associated with (G) λ,γ were bounded provided Ω was a convex domain in R N where N 3 and f and g satisfied conditions similar to (R). Regularity results regarding (H) λ,γ we also obtained in the case where at least one of of f and g were explicit nonlinearities given by (u + 1) p or e u. Remark 1. 1. After the completion of this work we were made aware of a recent improvement of the results from [7]. In [13] they examined (E) λ,γ and showed the extremal solutions are bounded provided N 9. They make no assumptions on the parameters λ, γ. In their approach they obtain, independent of our work, the stability inequality (4). We also mention a related work [10] which examined fourth order scalar problems but used similar techniques.. Since there was considerable delay from the completion of this paper to publication, we now mention various works that have appeared in the meantime. The ideas of this work were extended in [6] to prove various Liouville theorems regarding solutions of u = v p and v = u θ in 5

and R N +. The methods developed here and then extended to Liouville theorems in [6] were then extended to handle some Liouville theorems related to some Hénon systems in R N in [16] (this work also contains many other results). Regarding Lane-Emden system we finally mention the work [5] which obtained optimal Liouville theorems related to the Lane-Emden system u = v p, v = u θ in R N, at least in the case of radial solutions. We finally mention a very powerful new monotonicity method which was developed in [1] to prove Liouville theorems related to stable solutions of u = u p 1 u in R N. Since that work [1] work there has been many extensions of that method. R N Main Results We now state our main results. Theorem 1. Suppose that 1 < p θ, (λ, γ ) Υ and let (u, v ) denote the extremal solution associated with (N) λ,γ. Suppose that N (θ + 1) < 1 + pθ(p + 1) pθ(p + 1) pθ(p + 1) +. () pθ 1 θ + 1 θ + 1 θ + 1 Then u, v are bounded. Remark. We are interested in obtaining lower bounds on the right hand side of (), in the case where 1 < p θ, and so we set f(p, θ) := 4(θ+1) pθ 1 t 0, where t 0 := pθ(p + 1) pθ(p + 1) pθ(p + 1) +. θ + 1 θ + 1 θ + 1 Note that () is satisfied exactly when N < f(p, θ). We now rewrite f using the change of variables z = pθ(p+1) θ+1 to arrive at ) f(p, z) = 4p z p and the transformed domain is given by ( z + z z D = {(p, z) : p 1, p z p + p}. A computer algebra system easily shows that f > 8 on D. Also note that p f > 0 on D and so we see that f > 8 on D and hence we see that the extremal solutions are bounded for any 1 < p θ provided N 10., 6

There are two main steps in proving the above theorems. We first show that minimal solutions of (P ) λ,γ, which are semi-stable in the sense of (1), satisfy a stability inequality which is reminiscent of semi-stability in the sense of the second order scalar equations. This is given by Lemma 1. The second ingredient will be a pointwise comparison result between u and v, given in Lemma. We remark that this was motivated by [] and a similar result was used in [8]. Lemma 1. Let (u, v) denote a semi-stable solution of (P ) λ,γ in the sense of (1). Then f λγ (v)g (u)φψ φ + ψ, (3) for all φ, ψ H0 1 (Ω). Taking φ = ψ gives f λγ (v)g (u)φ φ (4) for all φ H 1 0 (Ω). Proof. Since (u, v) is a semi-stable solution of (P ) λ,γ there is some 0 < ζ, χ H0 1 (Ω) smooth such that ζ ζ λf (v) χ ζ, χ χ γg (u) ζ, in Ω. χ Let φ, ψ Cc (Ω) and multiply the first equation by φ and the second by ψ and integrate over Ω to arrive at λf (v) χ ζ φ φ, γg (u) ζ χ ψ ψ, where we have utilized the result that for any sufficiently smooth E > 0 we have E E φ φ, for all φ Cc (Ω). We now add the inequalities to obtain (λf (v)φ ) χ ζ + (γg (u)ψ ) ζ χ φ + ψ. (5) Now note that λγf (v)g (u)φψ tλf (v)φ + 1 t γg (u)ψ, 7

for any t > 0. Taking t = χ(x) ζ(x) gives λγf (v)g(u)φψ (λf (v)φ ) χ ζ + (γg (u)ψ ) ζ χ, and putting this back into (5) gives the desired result. Lemma. Let (u, v) denote a smooth solution of (N) λ,γ and suppose that θ p > 1. Define { ( ) 1 } γ(p + 1) p+1 α := max 0, 1. λ(θ + 1) Then λ(θ + 1)(v + 1 + α) p+1 γ(p + 1)(u + 1) θ+1 in Ω. (6) Proof. Let (u, v) denote a smooth solution of (N) λ,γ and define w := v + 1 + α C(u + 1) t where C := ( ) 1 γ(p + 1) p+1 θ + 1 and t := λ(θ + 1) p + 1 1. Note that w 0 on Ω and define Ω 0 := {x Ω : w(x) < 0}. If Ω 0 is empty then we are done so we suppose that Ω 0 is nonempty. Note that since w 0 on Ω we have w = 0 on Ω 0. A computation shows that w = γ(u+1) θ Ct(u+1) t 1 λ(v+1) p +Ct(t 1)(u+1) t u, in Ω, and since t 1 we have w γ(u + 1) θ Ct(u + 1) t 1 λ(v + 1) p in Ω. Note that we have, by definition, and so we have v + 1 v + 1 + α < C(u + 1) t in Ω 0, w γ(u + 1) θ C p+1 tλ(u + 1) tp+t 1 in Ω 0, but the right hand side of this is zero and hence we have w 0 in Ω 0 with w = 0 on Ω 0 and hence w 0 in Ω 0, which is a contradiction. So Ω 0 is empty. 8

Proof of Theorem 1. Let (λ, γ ) Υ and let σ := γ λ and suppose that (u, v) denotes a minimal solution of (N) λ,γ on the ray Γ σ. Put φ := (v + 1) t 1, where 1 < t, into (4) to obtain λγpθ (v + 1) p 1 θ 1 (u + 1) ((v + 1) t 1) t (v + 1) t v, and multiply (N) λ,γ by (v + 1) t 1 1 and integrate by parts to obtain t (v + 1) t v = t γ (u + 1) θ ((v + 1) t 1 1). t 1 Equating these and expanding the squares and dropping some positive terms gives λγpθ (v + 1) p 1 θ 1 (u + 1) (v + 1) t t γ (u + 1) θ (v + 1) t 1 t 1 + λγpθ (v + 1) p 1 θ 1 (u + 1) (v + 1) t. (7) We now use Lemma to get a lower bound for I := (v + 1) p 1 θ 1 (u + 1) (v + 1) t, but we need to rework the pointwise estimate (6) first. From (6) we have γ(p + 1) θ+1 (u + 1) (v + 1 + α) p+1, λ(θ + 1) and for all δ > 0 there is some C(δ) > 0 such that (v + 1 + α) p+1 (1 + δ)(v + 1) p+1 + C(δ)α p+1. From this we see that there is some C 1 = C 1 (δ, p, α) such that (v + 1) p+1 γ(p + 1) (u + 1) θ+1 C 1. λ(θ + 1) 1 + δ We now rewrite I as I = (u + 1) θ 1 (v + 1) t 1 (v + 1) p+1, 9

and use the above estimate to show that γ(p + 1) 1 I (u + 1) θ (v + 1) t 1 C 1 λ(θ + 1) δ + 1 (u + 1) θ 1 (v + 1) t 1. (8) We now return to (7) and write the left hand side, where ε > 0 is small, as ε λγpθi + (1 ε) λγpθi, and we leave the first term alone and we use the above lower estimate for I on the second term. Putting this back into (7) and after some rearranging one arrives at ε λγpθi + γk (u + 1) θ (v + 1) t 1 λγpθ I 1 + (1 ε) λγpθc 1 I (9) where K := I 1 := (1 ε) pθ(p + 1) 1 + δ θ + 1 t t 1, (v + 1) p 1 (u + 1) θ 1 (v + 1) t, and I := (u + 1) θ 1 (v + 1) t 1. For the moment we assume the following claims: for all T > 1 and k > 1 1 I (u + 1) θ (v + 1) t 1 + Ω T θ 1 T θ+1 k t 1 + 1 (u + 1) θ 1 p 1 k p+1 (v + 1) +t, (10) and I 1 1 T t (u + 1) θ 1 p 1 (v + 1) +t + Ω T p 1 +t k θ 1 + T p+1 t k θ+1 (u + 1) θ (v + 1) t 1. (11) Putting (10) and (11) back into (9) one arrives at an estimate of the form K 1 I + K (u + 1) θ (v + 1) t 1 C(ε, p, θ, T, k, δ) (1) 10

where K 1 := ε λγpθ λγpθ (1 ε) λγpθ C1 T t k p+1, and K := γk λγpθ T p+1 t (1 ε) λγpθ C 1, k θ+1 T θ+1 and where C(ε, p, θ, T, k, δ) is a positive finite constant which is uniform on the ray Γ σ. Define pθ(p + 1) pθ(p + 1) pθ(p + 1) t 0 := + θ + 1 θ + 1 θ + 1 and note that t 0 > 1 for all p, θ > 1. Fix 1 < t < t 0 and hence pθ(p + 1) t θ + 1 t 1 > 0. We now fix ε > 0 and δ > 0 sufficiently small such that K > 0. We now fix T > 1 sufficiently large such that ε λγpθ λγpθ T t (this is the first two terms from K 1 ) and γk (1 ε) λγpθ C 1 T θ+1 (the first and third terms from K ) are positive and bounded away from zero on the the ray Γ σ. We now take k > 1 sufficiently big such that K 1, K are positive and bounded away from zero on the ray Γ σ and hence we have estimates of the form: for all 1 < t < t 0 there is some C t > 0 such that (u + 1) θ (v + 1) t 1 C t, (13) where C t is some finite uniform constant on the ray Γ σ. Using the pointwise lower estimate (6) for v + 1 gives: for all 1 < t < t 0 there is some C t <, uniform along the ray Γ σ, such that (θ+1)(t 1) θ+( (u + 1) p+1 C t, (14) 11

and hence this estimate also holds if one replaces u with u. We now let 1 < t < t 0 and note that 1 v = (u + 1) θ v (u + 1) θ (v + 1) (u + 1) θ (v + 1) t 1 C t, γ by (13) and hence we can pass to the limit and see that v H0 1 (Ω). We now proceed to show that v is bounded in low dimensions. First note that v γ = (u + 1) θ = (u + 1) θ v + 1 v + (u + 1) θ v + 1 in Ω. To show that v is bounded it is sufficient, since v H0 1 (Ω), to show that (u +1) θ v +1 L T (Ω) for some T > N. Using (6) and passing to the limit one sees there is some C > 0 such that (u + 1) θ v + 1 C(u + 1) pθ 1 p+1 + C in Ω, and so (u +1) θ v +1 L T (Ω) for some T > N provided (pθ 1) N p + 1 < θ + (θ + 1)(t 0 1), p + 1 after considering (14). This rearranges into N < 1 + (θ + 1) p + 1 t 0, which is the desired result. We now use (N) λ,γ and elliptic regularity to see that u is also bounded. Proof of Claims (10) and (11). We first prove (10). We write I as I = + +, u+1 T u+1<t,v+1 k u+1<t,v+1>k where the integrands are the same as in I. Note that the first integral is less than or equal u+1>t (u + 1) θ 1 ( u + 1 T ) θ+1 (v + 1) t 1 1 T θ+1 (u + 1) θ (v + 1) t 1. 1

The second integral is trivial to get upper estimate on. One estimates the third integral in the same way as the first to see that 1 (u + 1) θ 1 p 1 u+1<t,v+1>k k p+1 (v + 1) +t. Combining these estimates gives (10). We now prove (11). We write I 1 = + +, v+1 T v+1<t,u+1<k v+1<t,u+1 k where the integrands are the same as I 1. Note that the first integral is less than or equal ( ) (v + 1) p 1 v + 1 t +t (u + 1) θ 1, T v+1>t and this is less than or equal 1 T t (v + 1) p 1 +t (u + 1) θ 1. The second integral is easily estimated. We rewrite the third integral as ( ) (v + 1) t 1 (u + 1) θ (v + 1) p+1 t (u + 1) θ 1, v+1<t,u+1 k and we now estimate the terms inside the bracket in the obvious manner. Combining these gives (11). References [1] H. Brezis, T. Cazenave, Y. Martel, A. Ramiandrisoa; Blow up for u t u = g(u) revisited, Adv. Diff. Eq., 1 (1996) 73-90. [] H. Brezis and L. Vazquez, Blow-up solutions of some nonlinear elliptic problems, Rev. Mat. Univ. Complut. Madrid 10 (1997), no., 443 469. [3] X. Cabré, Regularity of minimizers of semilinear elliptic problems up to dimension four, Comm. Pure Appl. Math. 63 (010), no. 10, 136-1380. [4] X. Cabré and A. Capella, Regularity of radial minimizers and extremal solutions of semilinear elliptic equations, J. Funct. Anal. 38 (006), no., 709 733. 13

[5] W. Chen, L. Dupaigne and Marius Ghergu, A new critical curve for the Lane-Emden system, Discrete and Continuous Dynamical Systems 013. [6] C. Cowan, Liouville theorems for stable Lane-Emden systems and biharmonic problems, Nonlinearity 6 (013) 357-371 [7] C. Cowan, Regularity of the extremal solutions in a Gelfand system problem, Advanced Nonlinear Studies, Vol. 11, No. 3, p. 695 Aug., 011. [8] C. Cowan, P. Esposito and N. Ghoussoub, Regularity of extremal solutions in fourth order nonlinear eigenvalue problems on general domains. Discrete Contin. Dyn. Syst. 8 (010), no. 3, 10331050. [9] C. Cowan and M. Fazly, Regularity of the extremal solutions associated to elliptic systems, J. Differential Equations 57(014)4087-4107 [10] C. Cowan and N. Ghoussoub, Regularity of semi-stable solutions to fourth order nonlinear eigenvalue problems on general domains, Calc. Var. Partial Differential Equations 49 (014), no. 1-, 91-305 [11] 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.07-18. [1] Juan Davila, Louis Dupaigne, Kelei Wang and J. Wei, A Monotonicity Formula and a Liouville-type Theorem for a Fourth Order Supercritical Problem, Advances in Mathematics 58(014), 40-85. [13] L. Dupaigne, A. Farina and B. Sirakov, Regularity of the extremal solutions for the Liouville system, Geometric Partial Differential Equations proceedings CRM Series Volume 15, 013, pp 139-144. [14] P. Esposito, N. Ghoussoub, Y. Guo; Mathematical Analysis of Partial Differential Equations Modeling Electrostatic MEMS, Research Monograph, Courant Lecture Notes in Mathematics, 0. Courant Institute of Mathematical Sciences, New York; American Mathematical Society, Providence, RI, 010. xiv+318 pp. [15] P. Esposito, N. Ghoussoub and Y. Guo, Compactness along the branch of semi-stable and unstable solutions for an elliptic problem with a singular nonlinearity, Comm. Pure Appl. Math. 60 (007), 1731 1768. 14

[16] M. Fazly and N. Ghoussoub, ON THE HÉNON-LANE-EMDEN CON- JECTURE, DISCRETE AND CONTINUOUS DYNAMICAL SYS- TEMS, Volume 34, Number 6, June 014 pp. 513-533 [17] N. Ghoussoub and Y. Guo, On the partial differential equations of electro MEMS devices: stationary case, SIAM J. Math. Anal. 38 (007), 143-1449. [18] Y. Martel, Uniqueness of weak extremal solutions of nonlinear elliptic problems, Houston J. Math. 3 (1997), 161-168. [19] F. Mignot, J-P. Puel; Sur une classe de problemes non lineaires avec non linearite positive, croissante, convexe, Comm. Partial Differential Equations 5 (1980) 791-836. [0] M. Montenegro, Minimal solutions for a class of elliptic systems, Bull. London Math. Soc.37 (005) 405-416. [1] G. Nedev, Regularity of the extremal solution of semilinear elliptic equations, C. R. Acad. Sci. Paris Sr. I Math. 330 (000), no. 11, 997 100. [] P. Souplet, The proof of the Lane-Emden conjecture in four space dimensions. Adv. Math. 1 (009), no. 5, 1409147. 15