arxiv: v1 [math.ap] 20 Dec 2018

Similar documents
Travelling waves. Chapter 8. 1 Introduction

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

Half of Final Exam Name: Practice Problems October 28, 2014

B5.6 Nonlinear Systems

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

A plane autonomous system is a pair of simultaneous first-order differential equations,

Nonlinear Autonomous Systems of Differential

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

Stability of Dynamical systems

ẋ = 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.

Nonlinear dynamics & chaos BECS

MATH 215/255 Solutions to Additional Practice Problems April dy dt

2.10 Saddles, Nodes, Foci and Centers

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

B5.6 Nonlinear Systems

Problem Set Number 5, j/2.036j MIT (Fall 2014)

STABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS

7 Two-dimensional bifurcations

Convergence Rate of Nonlinear Switched Systems

7 Planar systems of linear ODE

MCE693/793: Analysis and Control of Nonlinear Systems

Problem set 7 Math 207A, Fall 2011 Solutions

1 Lyapunov theory of stability

Lyapunov Stability Theory

Now I switch to nonlinear systems. In this chapter the main object of study will be

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL

Asymptotic Convergence of the Steepest Descent Method for the Exponential Penalty in Linear Programming

Part II. Dynamical Systems. Year

Section 9.3 Phase Plane Portraits (for Planar Systems)

ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)

The Liapunov Method for Determining Stability (DRAFT)

Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains

LMI Methods in Optimal and Robust Control

Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

Modeling and Analysis of Dynamic Systems

Calculus and Differential Equations II

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Mathematical Economics. Lecture Notes (in extracts)

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);

Discretized Fast-Slow Systems near Pitchfork Singularities

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

STABILITY. Phase portraits and local stability

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

Nonlinear Control Lecture 5: Stability Analysis II

Mathematical Modeling I

Stability Implications of Bendixson s Criterion

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait

Homoclinic and Heteroclinic Connections. for a Semilinear Parabolic Equation. Marek Fila (Comenius University)

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

arxiv: v2 [math.ap] 28 Nov 2016

2D-Volterra-Lotka Modeling For 2 Species

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

Chapter III. Stability of Linear Systems

MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation

Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia. 1. Introduction

Nonlinear Control Lecture 2:Phase Plane Analysis

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)

Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

Global unbounded solutions of the Fujita equation in the intermediate range

The Higgins-Selkov oscillator

ODE Final exam - Solutions

Introduction to the Phase Plane

2 A Model, Harmonic Map, Problem

Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Periodic Sinks and Observable Chaos

Threshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations

Hybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

A CONVEX-CONCAVE ELLIPTIC PROBLEM WITH A PARAMETER ON THE BOUNDARY CONDITION

THE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF

Course 212: Academic Year Section 1: Metric Spaces

One Dimensional Dynamical Systems

Canards at Folded Nodes

Practice Problems for Final Exam

Differential equations, comprehensive exam topics and sample questions

1 2 predators competing for 1 prey

Stability of Stochastic Differential Equations

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)

REUNotes08-ODEs May 30, Chapter Two. Differential Equations

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang

Example of a Blue Sky Catastrophe

Department of Mathematics IIT Guwahati

Perturbation Theory of Dynamical Systems

8.1 Bifurcations of Equilibria

In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation

Lecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Math 273 (51) - Final

BIFURCATION FOR NON LINEAR ORDINARY DIFFERENTIAL EQUATIONS WITH SINGULAR PERTURBATION

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

Transcription:

arxiv:181.8418v1 [math.ap] Dec 18 Radial solutions of scaling invariant nonlinear elliptic equations with mixed reaction terms arie-françoise Bidaut-Véron, arta Garcia-Huidobro Laurent Véron Abstract We study global properties of positive radial solutions of u = u p + u p in R where p > 1 and is a real number. We prove the existence or the non-existence of ground states and of solutions with singularity at. 1 athematics Subject Classification. 35J6, 35B8, 684. Key words. elliptic equations; ground states; linearization; homoclinic and heteroclinic orbits; Floquet integral. Contents 1 Introduction General properties of the system 8.1 Reduction to autonomous equation and system................... 8.1.1 The standard reduction............................ 8.1. The geometry of the vector field H...................... 8.1.3 Graphic representation of the vector field H................. 9.1.4 Other reduction................................ 11. Regular solutions and ground states......................... 1.3 Explicit singular solutions............................... 13.3.1 Upper estimate of the regular solutions................... 16.4 Linearization of the system 1.17) near equilibria.................. 17.4.1 Linearization at, )............................. 17.4. Linearization at a fixed point P := X, Y )................4.3 Energy and Lyapunov functionals for system 1.17)..............4.4 Comparison results............................... 4 Laboratoire de athématiques et Physique Théorique, Université de Tours, 37 Tours, France. E-mail: veronmf@univ-tours.fr Departamento de atematicas, Pontifica Universidad Catolica de Chile Casilla 37, Correo, Santiago de Chile. E-mail: mgarcia@mat.puc.cl Laboratoire de athématiques et Physique Théorique, Université de Tours, 37 Tours, France. E-mail: veronl@univ-tours.fr 1

Quasilinear elliptic equations with source reaction 3 Study of ground states of 1.17) when > 4 3.1 Behaviour near equilibrium.............................. 4 3. Existence or nonexistence of ground states...................... 8 4 Study of ground states of 1.17) when < 31 4.1 Ground states and large solutions when <................... 39 4. The case <, 3 and p +........................ 43 4.3 The case <, 3 and < p < +.................... 43 4.4 The case < and 1 < p <.......................... 45 4.5 The case <, = and p > 1.......................... 48 4.6 The case <, = 1 and p > 1.......................... 49 4.7 Appendix: on-existence of cycle in the case =................ 5 1 Introduction The aim of this article is to study local and global properties of positive radial solutions of the equation u = u u + u p, 1.1) in R or R \ {} where p > 1 and is a real parameter. This is a particular case of the following class of equations u = u u + u q, 1.) where q > 1 which has been the subject or many works in the radial case when <, where a basic observation is that the two terms u u and u q are in competition. The first work in that case is due to Chipot and Weissler [14] who, in particular, solved completely the case = 1, then Serrin and Zou [19] performed a very detailed analysis. uch less is known in the case >. Under the scaling transformation T k defined for k > by 1.) becomes Therefore, if q p, 1.) can be reduced to u k := T k [u]x) = k ukx), 1.3) u k = u k u k + k p q) u k q, 1.4) u = u u ± u p. 1.5) oreover, when q < p, the limit equation of 1.4) when k is the Lane-Emden equation the exponent p is dominant. The other scaling transformation u = u u, 1.6) v k := S k [u]x) = k q q 1 ukx), 1.7)

Quasilinear elliptic equations with source reaction 3 transforms 1.) into v k = k q p q) q 1 v k v k + v k q, 1.8) and if q > p, the limit equation of 1.8) when k is the Riccati equation v = v q, 1.9) the exponent q is dominant. In [14] and [19] most of the study deals with the case q p. In the critical case i.e. when q = p p + 1, 1.1) then not only the sign of but also its value plays a fundamental role, with a delicate interaction with the exponent p. otice that an equivalent form of 1.1) is v = λ v v ± v q 1.11) with λ >. In the critical case first studies in the case < are due to Chipot and Weissler [14] for = 1. The case was left open by Serrin and Zou [19] and the first partial results are due to Fila and Quittner [16] and Voirol [, 3]. The case > was not considered. The equation 1.1) is the stationary part of the associated parabolic equation t u u u u u q =. 1.1) which is studied in [14] and [15], one of the aim being to find conditions for the blow-up of positive solutions. A general survey with several open problems can be found in []. In the non radial case an important contribution dealing with a priori estimates of local positive solutions of 1.) and existence or non-existence of entire positive solution in R is due to the authors [7]. In this paper we complete the results of [7] in presenting a quite exhaustive study of the radial solutions of 1.1) for any real number. The radial solutions of 1.1) are functions r ur) defined in, ) where they satisfy u rr 1 u r = u u + u r p. 1.13) r Because of the invariance of 1.13) under the transformation T k there exists an autonomous variant of 1.1) obtained by putting ur) = r xt) with t = ln r, 1.14) then with K = x tt + Lx t K x + x x + )p and L = p K x x t = 1.15) )p + ) = K. 1.16)

Quasilinear elliptic equations with source reaction 4 Putting yt) = r u r r), then xt), yt)) satisfies the system x t = x y y t = Ky + x x + y p. 1.17) We are mainly interested in the trajectories of the system which remain in the first quarter Q={x, y) R : x >, y > }. Among these trajectories, the ones corresponding to ground states, i.e. positive C solutions u of 1.13) are defined on [, ). They verify u r ) = and actually they are C on, )). Using the invariance of the equation under T k all the ground states can be derived by scaling from a unique one which satisfies u) = 1. Since it is easy to prove that such a solution u is decreasing, in the variables x, y), a ground state is a trajectory of 1.17) in Q, defined on R and satisfying denoted by T reg. lim t yt) xt) =. The corresponding trajectory is Contrarily to the case of the Lane-Emden equation 1.6), there exists no natural Lyapunov function when. This makes the study much more delicate and it is based upon a phase plane analysis. The solutions of 1.13) invariant under T k for any k > correspond to constant solutions of 1.15) and have the form where X is a positive root of X + Ur) = Xr for all r >, 1.18) ) p X + K =. 1.19) This equation plays a fundamental role in the description of the set of solutions of 1.13). The following constant, defined for = 1, and p > 1 or 3 and 1 < p has an important role in the description of the set roots of 1.19), ) p )p µ := µ ) = p + 1). 1.) p This set is described in the following proposition. Proposition 1 1- If, equation 1.19) admits a positive root, necessarily unique, if and only if 3 and p >. - If < and p, equation 1.19) admits a unique positive root X. 3- If < and either = 1, and p > 1 or 3 and 1 < p, there exists no positive root of 1.19) if µ < <, a unique positive root if = µ < and two positive roots X 1, < X, if < µ. We also set Y = X and P = X, Y ) resp. Y j, = X j, and P j, = X j,, Y j, ), for j=1,) and define the corresponding singular solutions U r) = X r resp. U j, r) = X j, r ).

Quasilinear elliptic equations with source reaction 5 Recall briefly the description of the positive solutions of the Lane-Emden equation 1.6), i.e. = : there exists radial ground states if and only if 3 and p + +. If p = these ground states are explicit and they satisfy lim r r ur) = c >. There exist infinitely many singular solutions u ondulating around U. ote that a ground state corresponds to a homoclinic orbit at for system 1.17) and these singular solutions are cycles surrounding P. We recall that an orbit of 1.17) which connects two different equilibria resp. the same equilibrium) when t R is called heteroclinic resp. homoclinic). Without the radiality assumptions, and using a delicate combination of refined Bernstein techniques and Keller-sserman estimate we have obtained in [7, Th. C, Th. D] a series of general a priori estimates for any positive solution of 1.1), in an arbitrary domain of R in the case 1, p > 1 and q = > := no ground state. ) p ) 4p and >. In particular we proved there that if p > 1 and ) p, or if, p < +3 1 and > equation 1.1) admits In the sequel we describe the ground states and the singular global solutions of 1.13) in R \ {}. Concerning the ground states, we discuss according to the sign of and the value of p. The next value of appears when we linearize the system 1.17) at the equilibrium P, =, p) = p + 1) )p ) 4p) p )) + 4) 1. 1.1) Then is positive resp. negative) if p > + + resp. p < and we set µ = ). It is easy to see that if = then the characteristic values of the linearized operator at P are purely imaginary. otice that is positive resp. negative) according p > + resp. p < + ). Theorem A Let 1, p > 1 and >. 1- If 1 < p + if 3, and any p > 1 if = 1,, then equation 1.13) admits no ground state. - If 3 and p > +, there exist constants min, max verifying < < min max, such that - if < < min there exist ground states u satisfying ur) U r) when r. - if = min or = max there exists a ground state u minimal at infinity, that is satisfying lim r r ur) = c >. - for > max there exists no radial ground state. The values of min and max appear as transition values for which the ground state still exists but it is smaller than the others at infinity; it is of order r instead of r. They are not explicit but they can be estimated in function of and p. It is a numerical evidence that min = max in the phase plane analysis of system 1.17) and we conjecture that this is true. When = min or max, the system 1.17) admits homoclinic trajectories. We prove that the system 1.17) admits a Hopf bifurcation when =. When p > + we also prove the existence of different types of positive singular solutions

Quasilinear elliptic equations with source reaction 6 Theorem A Let 3. 1- If < p + for any > there exist a unique always up to scaling transformation) positive singular solutions u of 1.13) satisfying lim r ur) = X and lim r ur) = c > r r - If p > +, then i) If > max, there exists a unique singular solution u of 1.13) with the same behaviour as in 1. ii) If < < min there exist positive singular solutions u ondulating around U on R. In terms of the system 1.17) the 1) and -i) correspond to the existence of a heteroclinic orbit in Q connecting P to, ) and ii) to the existence of a cycle in Q surrounding P. When is negative, the precise description of the trajectories of 1.17) depends also on the value of p with respect to +. It is proved in [7, Th. B, E] that for 3 and 1 < p < there exists ɛ > such that if ɛ equation 1.1) admits no positive solution in R. The same conclusion holds if 3, 1 < p or = and p > 1) and > µ. We first consider the case p for which there exists a unique explicit singular solution U, and the results present some similarity with the ones of Theorem A. Theorem B Let 3, p and <. Then 1- If p +, then equation 1.13) admits ground states u. oreover they satisfy ur) U r) as r. - If p < +, there exist numbers µ min and µ max verifying < µ < µ min µ max < µ 1), such that i) for < µ max there exist ground states u such that ur) U r) when r. ii) for = µ min or for = µ max there exist ground states minimal at infinity in the sense that ur) cr when r, c >. iii) for µ min < < there exists no radial ground state. Here also the value of µ min, µ max are not explicit and we conjecture that they coincide. The next result presents some similarity with Theorem A. Theorem B Let 3 and < p < +. i) If < < there exists a unique up to scaling) positive singular solution u of 1.13), such that ur) U r) when r and ur) cr when r for some c >. ii) If µ min < < there exist positive singular solutions u ondulating around U on [, ) and singular solution ondulating around U in a neighbourhood of and satisfying ur) cr for some c > when r. In terms of the system 1.17), i) corresponds to a heteroclinic orbit connecting P and, ), while ii) to the existence of a periodic solution in Q around P, and the existence of a solution in Q converging to, ) at and having a limit cycle at t = which is a periodic orbit around P.

Quasilinear elliptic equations with source reaction 7 The situation is more complicated when 1 < p < and < µ because there exist two explicit singular solutions U 1, and U, which coincide when = µ. Theorem C Let <, 3 and 1 < p < constants µ min and µ max verifying such that, or = and p > 1. Then there exist two µ µ < µ min µ max < µ 1), 1- If < µ max then equation 1.13) admits ground states u either ondulating around U, or such that ur) U, r) as r. - If = µ min or = µ max there exists a ground state u such that ur) U 1, r) as r. 3- If µ min < < there exists no radial ground state. Here again µ min and µ max appear as transition values for which the ground state still exists but it is smaller than the others at infinity: it behaves like U 1 instead of U. The proof of this theorem is very elaborate in particular in the case =. In the case = 1 the result is already proved in [14]. The nonexistence of ground state, not necessarily radial for > µ is proved in [1] and independently in [7] with a different method. In the radial case it was obtained much before in the case = 1 in [14] and then by Fila and Quittner [16] who raised the question whether the condition µ min < < is optimal for the non-existence of radial ground state. This question received a negative answer in the work of Voirol [] who extended the domain of non-existence to µ ɛ < µ. The next result is the counterpart of Theorem C when dealing with singular solutions. Theorem C Let <, 3 and 1 < p <. i) If < µ there exist positive singular solutions u such that ur) U 1, r) as r and ur) cr with c > when r. ii) If < µ there exists a unique up to scaling positive singular solution u, such that ur) = U, r) as r and ur) = U 1, r) as r. Furthermore ur) > U 1, r) for all r >. iii) If ˆµ min < < µ there exist positive singular solutions u ondulating around U, at and such that ur) U 1, r) as r, and positive singular solutions u ondulating around U, on R. iv) If = µ min or = ˆµ max there exists a positive singular solutions u different from U 1, such that ur) U 1, r) when r and r. v) If µ min < < ˆµ max there exists a positive singular solution u such that lim r ur) = r c > and either ondulating around U, or such that ur) U, r) when r. vi) If 3 and = µ, there exist positive singular solutions u satisfying lim r ur) = r c > and ur) U µ r) as r. In terms of the system 1.17) i) corresponds to a heteroclinic orbit connecting P 1, to, ); ii) to a heteroclinic orbit connecting P, to P 1, ; iii) to a trajectory having a periodic orbit around P, for limit cycle at and converging to P, at and to a periodic orbit around

Quasilinear elliptic equations with source reaction 8 P, ; iv) corresponds to homoclinic orbit at P 1, ; v) corresponds to a trajectory connecting, ) at and either converging to P, at or having a periodic orbit around P, for limit set at ; vi) corresponds to a heteroclinic orbit connecting from, ) to P µ. Acknowledgements This article has been prepared with the support of the collaboration programs ECOS C14E8 and FODECYT grant 11654 for the three authors. General properties of the system.1 Reduction to autonomous equation and system.1.1 The standard reduction We recall that if u is a C 3 function defined on some interval I [, ) verifying 1.13) and if then x satisfies the autonomous equation ur) = r xt) with t = ln r, x tt + Lx t K x + x x + p x x t =,.1) on lni) where K and L are defined in 1.16). Setting u r = r yt), then xt), yt)) satisfies { xt = H 1 x, y).) y t = H x, y), where H 1 x, y) = x y H x, y) = Ky + x x + y p. and we denote by H the vector field of R with components H 1 and H..3).1. The geometry of the vector field H Let us denote by Q := {x, y) : x >, y > } the first quadrant. The vector field is inward in resp. outward of) Q on the axis {x, y) : x >, y = } resp. {x, y) : x =, y > }). We set { L := x, y) Q : y = x } { ) and C := x, y) Q : x = Ky y p 1 } p.4) ) and ψy) = Ky y p 1 p. Then x t = on L and y t = on C. The curves L and C have [ zero, one or two intersections in Q according the value of K and. If, K >, then ) 1 ) C, p p K ] [ ) p) p, ] ) K. The points, ), P and, ) K ) belong to C. The function ψ is increasing on, K q ) ) ) and decreasing on K q, ) K ).

Quasilinear elliptic equations with source reaction 9 If < and K >, ψ is concave and increasing on, ) with ψy) = ) 1 p y 1 + o1)) as y. If < and K <, ψ is still concave and increasing on ) K, ) with the same asymptotic as above. We quote below the possible connected components of Q \ L C). { } } x A = x, y) : y < {x, y) : Ky + x p + y p < = {x, y) : x t <, y t < }. { } } x B = x, y) : y > {x, y) : Ky + x p + y p < = {x, y) : x t >, y t < }. { } } x C = x, y) : y > {x, y) : Ky + x p + y p > = {x, y) : x t >, y t > }. { } } } x D = x, y) : y < {x, y) : Ky + x p + y p > {x, y) : x > X }, = {x, y) : x t <, y t > } {x, y) : x > X,. { } } } x E = x, y) : y < {x, y) : Ky + x p + y p > {x, y) : x < X } 1, = {x, y) : x t <, y t > } {x, y) : x < X 1,. These connected components are A, B, C, D if K, < or K, >. A, C, D if K < and µ < <. A, C, D, E if K < and = µ. A, B, C, D, E if K < and < µ..1.3 Graphic representation of the vector field H We present below some graphics of the vector field H associated to system 1.17). 3 y C D) x t < y t > L A) x t < y t < P C) x t > y t > B) x t > y t < o x Figure 1: >, K >.

C) x t > y t > Quasilinear elliptic equations with source reaction 1 o x y C D) x t < y t > A) x t < y t < L P x t > y t < C) x t > y t > B) o x Figure : <, K. y C D) x t < y t > A) x t < y t < L C) x t > y t > o x y A) x t < y t < D) Figure 3: µ < <, K <. C x t < y t > L P x t > y t < C) x t > y t > B) o x

x t > y t < B) C) x t > y t > x E) t < y t > P 1, Quasilinear elliptic equations with source reaction 11 o x y C D) x t < y t > A) x t < y t < L E) x t < y t > P µ C) x t > y t > o 1 x Figure 4: = µ, K <. y A) x t < y t < C D) x t < y t > L P, x t > y t < B) C) x t > y t > x t < y t > E) P 1, o x.1.4 Other reduction y Figure 5: < µ, K <. A) x t < y t < The following change of unknowns, already used in [9] when =, C D) x t < y t > L σt) = yt) xt) = r u rr) ur) E) x t < y t > and zt) P µ = x xt) = r u ur),.5) yt) u r r) C) x t > y t > o 1 x

Quasilinear elliptic equations with source reaction 1 and valid if u r, transforms 1.17) into a Kolmogorov system with vector field V = V 1, V ) σ t = σ σ + + z + σ 1 ) σz = V 1 σ; z) z t = z pσ z σ 1 ).6) σz = V σ; z). Since σ and z are in factor the two axis {σ = } and {z = } are trajectories, actually not admissible for.) in view of.5). The system is singular on these two axis however it can be desingularized in setting σ = σ k+1 and z = z k+1 for some integer k > p + 1, which transforms.6) into a new nonsingular Kolmogorov system, σ t = 1 k + 1 σ z t = 1 k + 1 z σ k+1 + + z k+1 + p σ k+1 z k+1 σ σ) k+1 z 1 ) k+1 σ σ) k+1 z 1 ) k+1. Therefore no other trajectory can intersect them in finite time and the quadrant Q := {σ, z) : σ >, z > } is invariant. Furthermore σz = r u. It is noticeable that if = the initial system is quadratic and regular. The system.6) will be used in the most delicate cases. It corresponds to the differentiation of the initial equation.1).. Regular solutions and ground states Definition.1 A regular solution of 1.13) is a C solution defined on some maximal interval [, r ) satisfying u) = u > and u r ) =. A ground state is a nonnegative C solution defined on [, ). The existence and uniqueness of a regular solution is standard by the Cauchy-Lipschitz integral method. If u is a C solution it satisfies u r < on, r ). Indeed r 1 u r r) is decreasing near, hence u r < on some maximal interval, r 1 ), r ) and u r r 1 ) = if r 1 < r. If ur 1 ) = then u by uniqueness. If ur 1 ) > then u rr < near r 1 which would imply that ur) < ur 1 ) for r 1 ɛ r < r 1 which contradict the negativity of u r on, r 1 ). Hence ur 1 ) < which implies that u r r) < on the maximal interval, r ) where u >. Thus, if u is a ground state u r < on, ). Hence the trajectory of a ground state expressed in the system 1.17) lies in Q and expressed in the system.6) it lies in the quadrant Q. It is easy to check that the regular solution u := u u, such that u) = u satisfies ) ur) = u 1 u r p + 1) u r ) 1 + o1)) as r..8) 4p + ) + )p + ) q Under the scaling transformation T k, u can be transformed into the regular solution 1.13) u := u 1 satisfying u) = 1. If one considers the system 1.17) the transformation T k becomes the time shift which transforms t xt), yt)) into t xt + ln k), yt + ln k)), and the.7)

Quasilinear elliptic equations with source reaction 13 trajectory xt), yt)) of 1.17) corresponding to a ground state is therefore uniquely determined and denoted by T reg and satisfies yt) lim xt), yt)) =, ), lim t t xt) = yt) and lim t x p t) = 1..9) Hence in the system.5) there holds on the corresponding trajectory lim σt) =, lim zt) =..1) t t.3 Explicit singular solutions Explicit self-similar solutions of 1.13), necessarily under the form u = Ar, play a fundamental role in the study, whenever they exist. The following result covers Proposition 1. Proposition. 1- Let. Then there exists a unique self-similar solution of 1.13) if and only if 3 and p >. We denote it by U r) = X r, where X > depends also on, and p. To this solution is associated the equilibrium X, X ) of the system 1.17). Furthermore the mappings X is continuous and decreasing on [, ), X increasing and there holds i) X = ii) ) 1 K, ) K 1 1 ) p ) K p ) + X is ).11) K. - Let <. If 3 and p there exists a unique self-similar solution of 1.13) U r). The mapping X is continuous and decreasing, X holds { ) 1 ) } K max, p) X is decreasing and there ) 1 ) ) K + p)..1) 3- Let <. If = 1, and p > 1, or 3 and 1 < p <, there exists no self-similar solution of 1.13) if µ < < where µ = µ ) > is defined in 1.). If = µ there exists a unique self-similar solution U µ r). If < µ < there exist two-self-similar solutions U 1, r) = X 1, r and U, r) = X, r with X 1, < X,. Furthermore the mappings X 1, and X, are continuous, respectively increasing and decreasing on

Quasilinear elliptic equations with source reaction 14, µ ), while X1, and X, Furthermore there holds for large enough; are respectively decreasing and increasing. i) X j, µ = ii) iii) ) p K p ) K < X1, < ) ) p) ) 1, ) K K 1 1 p ) 1 K < X, < )K ) ), ) ) p)..13) Proof. The function U = X r is a self-similar solution of 1.13) if and only if X is a positive root of ) p f x) := x + x K =..14) Equivalently P = X, Y ) = Y, Y ) is a fixed point of system 1.17), where Y is the positive root of ) p f y) = y + y K =..15) The use of the variable y is a little easer than x. Since X is explicit if =, we shall study the cases. 1- Case >. If K >, equivalently p >, and, f is an increasing function tending to at and negative at y =. Hence Y is the unique positive root of.15). If K <, f is positive on [, ), hence no such solution exists. Since ) p Y + Y K =, by the implicit function theorem, Y is C 1. For > >, f y) > f y) for all y >. Hence Y is decreasing on [, ). Actually the expression of f shows more, namely that Y and ) p ) K + Y is increasing on [, ). Furthermore Y < K = Y < > K = Y > ) K ) K, 1 ) p ) 1 K p ), +

Quasilinear elliptic equations with source reaction 15 and we get.11). - Case <. Clearly f has a minimum at y = y, where y, = p + 1 ) 1 ) p and f y, ) = p ) p K..16) p + 1 -i) If 3 and p, then f ) <, f is decreasing on, y, ) and increasing on y,, ), hence f y) = has a unique positive root Y > y,. Since for < < and y >, f y) > f y), the mapping Y is continuous and decreasing and Y is increasing. Then the left-hand side of.1) is clear. ext we put A p = ) p, B = K, ξ = X, a = B A and η = ξ A. Then φη) = η η a =. Since φ1 + a) = 1 + a) 1 a a 1 + a) 1 a a p a 1 + a) 1 + a) + a p )) we derive η 1 + a, which implies the right-hand side of.1). -ii) If = 1, or 3 and 1 < p <, then K <. Hence, if f y, ) >, or equivalently µ <, the equation f y) = has no positive root, if = µ, it has a double root Y µ where Y µ = ) p K p ) 1 = µ Y µ = p + 1 p K,.17) and if < µ the equation.15) has two positive roots < Y 1, < Y, and since f does not vanish at Y j,, they are C 1 functions of, ), respectively increasing and decreasing. Since, K <, we obtain from f Y ) =, and ) p Y Y 1, > ) K,, < Y, = Y, < For a sharper estimate, we have for large enough, Ỹ = ) K 1 K )K ) ) p) ) ) = f Ỹ ) <.

Quasilinear elliptic equations with source reaction 16 Hence ) ) K K < Y1, < 1 K )K ) ). Similarly ) ) p) 1 K 1 p ) < Y, < ) ) p). The estimates.13) follow..3.1 Upper estimate of the regular solutions We first recall the following estimate in the case, consequence of the fact that the positive solutions of 1.1) are superharmonic and proved in a more general setting in [7, Prop..1 ]. Proposition.3 1- There exists no positive solution of 1.13) in R, ), R if and either = 1, and p > 1 or 3 and 1 < p. In particular there exists no ground state. - If 3, p >, and u is a positive solution of 1.13) in R, ), R. Then there exists ρ R such that and ur) min u r r) min { ) 1, { ) ) 1 p ) ), p ) ) Furthermore, if R =, inequalities.18), and.19) hold with ρ =. ) } r, for all r > ρ,.18) ) } r for all r > ρ..19) The next estimate is verified by any ground state, independently of the sign of. Proposition.4 Let p > 1 and 1. Then the ground state u of 1.13) with u) = 1 satisfies } ur) min {1, c,p, r and u r r) c,p, r r >..) Proof. The trajectory T reg starts from, ) and enters the region C. If it stays in C, then x is increasing on R. Since y < x, if xt) tends to some finite limit as t, it implies that the limit set of the trajectory exists. It cannot be a cycle since x is monotone, thus it is one of the equilibrium of the system. Hence xt) X j, for j = 1 or if K < and < µ or xt) X if K and either > or µ. This implies that.) holds. If xt) tends to, so y does because yt) x p t) K xt) xt) ) p as t.

Quasilinear elliptic equations with source reaction 17 Therefore y t Cy p for some C > which would imply that t y 1 p t)+c)t is increasing, which is impossible. ext we suppose that the trajectory leaves C by crossing the line L. Since it cannot enters B through C in the case K <, µ ), it leaves C by intersecting L and we denote by t 1 the first time where T reg intersects L. Then x t t 1 ) = and x tt t 1 ) = y t t 1 ) <. Therefore t 1 is a local maximum. ow the trajectory cannot cross again the half line { x, y) : x = xt 1 ), y > xt 1) for any t t 1. } because on it there holds x t = xt 1) y <. Hence xt) xt 1) In the same way, either y is increasing on R and since x t = x y and x is bounded, y cannot tend to infinity when t, thus yt) y >, or y is not monotone and T reg crosses C at a first value t, necessarily larger than t 1 and where x t t ) <. Then y tt t ) = p xt ) x t t ) < and t is a local maximum of y. Therefore t t xt), yt)) remains in the subset of Q bordered by the portion of trajectory of T reg for t t and {x, y) : < x xt ), y = yt )}. This implies that yt) yt ) for all t R. oticing that ur) 1 since u is decreasing, we get the conclusion. Remark. The above method does not give an explicit estimate of the upper bounds of x and y and such a bound can be estimated in some cases. If it follows from Proposition.3 that for any p > 1 there holds ) 1 ur) r, thus this estimate is independent of. If µ ) < <, it is easy to check that the function v = ln u satisfies v ae )v.1) with a = 1 µ )) >. We derive that the function wr) = r 1 v r r)) is increasing, with limit l, ]. Hence This implies e 1 p)vr) e 1 p)v) + a)r 4 } ur) min {1, c a 1 r..).4 Linearization of the system 1.17) near equilibria.4.1 Linearization at, ) The linearized system at, ) is x t = x y y t = Ky,.3)

Quasilinear elliptic equations with source reaction 18 which admits the eigenvalues ote that λ λ 1 =. λ 1 = K λ =..4) a) Assume that 3 and p >, equivalently K >, Then, ) is a saddle point. There exists a unique unstable trajectory T unst such that yt) lim xt), yt)) =, ) and lim =,.5) t t xt) and more precisely, from.8), xt) = e t 1 et p + 1) 4p+ ) e t 1 + o1)) p + 1) + p) p + 1) ) yt) = e pt 1 p + 1) + p + 1) + p e )t 1 + o1)) as t..6) From Definition. and the lines which follow, the unstable trajectory T unst coincides with the regular trajectory T reg. This is included in the region C for t < T for some < T. There exists also a unique stable trajectory T st such that yt) lim xt), yt)) =, ) and lim t t xt) = + K =..7) Since >, T st belongs to the region A for t > T 1 for some T 1 <. Remark. If T reg Q satisfies lim xt), yt)) =, ), then T reg = T st and the corresponding t solution is a ground state. The same conclusion holds, if T st Q satisfies xt), yt)) = lim t, ). Such a solution is called a homoclinic orbit at, ). Because of the uniqueness of the stable and unstable trajectories of a saddle point, it is unique in the class of solutions satisfying.7). Equivalently the class of ground states u of 1.13) satisfying ur) cr for some c > is then a one parameter family characterized by u) = u. b) Assume that 3 and 1 < p <. Then K < and < λ 1 < λ. Hence, ) is a source and all the trajectories of 1.17) in some neighbourhood of, ) converge to, ) when t. Among those trajectories there exists one fast trajectory which satisfies yt) lim t xt) =. It is actually the regular trajectory T reg. There exist also infinitely many slow trajectories which satisfy yt) lim t xt) =. c) If p =, then K = and λ 1 = < λ =. We still find the regular trajectory T reg associated to λ and the corresponding eigenvector 1, ). By the central manifold theorem

Quasilinear elliptic equations with source reaction 19 corresponding to λ 1 there exists an invariant curve passing through, ) with slope. Using the matched asymptotic expansion method, one finds that if < there exists a solution x of.1) such that xt) C 1 1 t 1 when t, i.e. ur) C 1 1 r ln r)) 1 when r, and if > there exists a solution x of.1) such that xt) C 1 1 t) 1 when t, equivalently ur) C 1 1 r ln 1 r ))1 when r. d) If =,, ) is a source with λ 1 = λ =, with corresponding eigenspace 1, ). The linearized problem is equivalent to equation x tt 4 q 1 x 4 t + q 1) x = with general solutions xt) = ae t q 1 + bte t q 1 for some real parameters a, b. Hence there exists infinitely trajectories of.1) tending to when t and they are tangent to 1, ) at, ). The regular trajectory T reg corresponds to b = and the other trajectories correspond to singular solutions u of 1.1). They satisfy ur) b ln r as r and there holds u = u p + u p πbδ, in the sense of distributions in B ɛ for some ɛ >. ext we give a general result in case the system admits only one equilibrium in Q. and R. If u is a regular solution the following tetra- Lemma.5 Let 3, p > chotomy occurs: i) either lim r r ur) = c for some c >, ii) or ur) U r) as r, iii) or ur) has an ω-limit cycle surrounding P, iv) or ur) changes sign for some r >. Proof. By assumption P is the unique equilibrium. The trajectory T reg starts from, ) and remains in the region C where x t, y t > for t t. If t =, u is a ground state, hence it is bounded from Proposition.4. Its ω-limit set is non-empty. Because x and y are monotone, it converges when t to some point which is necessarily P. If t <, then at t = t the trajectory leaves C through L since it cannot enter in B, and it enters the region D where x t <, y t >. oreover x t t ) = and xt ) > X. Then three possibilities occur: α) either xt) X monotonicaly when t ; thus the trajectory converges to P. β) either xt) monotonicaly when t. Since, ) is a saddle point, then T reg =T st. This implies that T reg is a homoclinic trajectory at, ). γ) or there exists t 1 > t such that x t t 1 ) =. Then xt 1 ) < X. Hence T reg enters the region B and by continuity there exists t < t such that xt ) = xt 1 ) and yt ) < yt 1 ). Therefore the bounded region of R bordered by the segment I = {x, y) : x = xt 1 ), yt ) < y < yt 1 )} and the portion of T reg defined by {xt), yt)) T reg : t t t 1 } is positively invariant notice that x t > on I) and it contains P and no other equilibrium. Therefore either the trajectory converges to P or it admits an ω-limit cycle which is a closed orbit surrounding P.

Quasilinear elliptic equations with source reaction.4. Linearization at a fixed point P := X, Y ) Suppose that P or P j, ) exists. Then setting x, y) = X, Y )+x, y), the linearized system at this point is x t = x y y t = px x + p Y K ) y,.8) The eigenvalues of its matrix are the roots of ) p T λ) = λ p + 1 Y L λ + K p p + 1 Y,.9) using equation.15). If is such that p p = L, then K = K L = K + =, and we denote by such a value of which is characterized by Y Y p p + 1 Y = L..3) Since Y is a positive root of.15), we get ) p Y = A := K p + 1 p L = )) + 4, p) hence, well defined for, is given by = p + 1)L p 1 A ) p ) 1 p + 1) )p ) =..31) 4p) p )) + 4) 1 That is 1.1). Remark. We see that > resp. < ) if and only if 3 and p > + resp. 1 < p < + if 3 and any p > 1 if = ). If = and 3, the eigenvalues of the linearized system are purely imaginary. If p = +, then =. It is known that in that case the point P = ), ) ) is a center for the system 1.17) associated to u = u + and that there exist infinitely many cycles turning around P with equation 1 y + x xy = E..3) It can be verified that E <, in particular using the function F defined in.35). For and 1 < p <, there always holds µ and more precisely Lemma.6 If 3 and 1 < p <, then < µ. If = and p > 1, then = µ.

Quasilinear elliptic equations with source reaction 1 Proof. From Proposition. and identities.17) and.9), there holds, µ < µ µ Y µ < µy K < L. But L K = series of equivalences. and the conclusion follows. If = we just replace < by = in the above Lemma.7 Assume 3. If p > + and = > then P is a weak sink and a Andronov-Hopf bifurcation point. If p < + and = < then P is a weak source. Proof. We recall that a weak sink is an asymptotically stable equilibrium which attracts the nearby points as t at a rate slower than the usual exponential rate. A weak source is a weak sink of the system obtained by changing t into t see [17, Chap. 9]). We write x = x X, ȳ = y Y and obtain the new nonlinear system where with c 1 = x t = x ȳ ȳ t = px x + p Y h x, ȳ) = c 1 x + c ȳ + c 3 x 3 + c 4 ȳ 3 +..., p) X p, c p) = p + 1) Y ȳ Kȳ + h x, ȳ),.33) p)p ) c 3 = X p 3 6, c p) 4 = 3p + 1) 3 Y Setting α = 4 +, we have px ) = α, since X = Y and satisfies.3), thus c 1 = α )L )L, c =, c 3 = Y p + 1)Y 3p + 1) Y..34) In order to compute the Lyapunov coefficients we transform the system by setting The new system is γ =, s = γt, γw = x ȳ. x s = w ) w s = x 1 h x, x γ γw. By integrating the first line and using the expansion of h, we obtain [ w s = x 1 ) x ) x 3 γ c 1 x + c γw + c 3 x 3 + c 4 γw +... p+3.

Quasilinear elliptic equations with source reaction This can be written in the following way w s = x + ν, x + ν 1,1 xw + ν, w + ν 3, x 3 + ν,1 x w + ν 1, xw + ν,3 w 3 +... By [17, Th. 9..3], the Lyapunov coefficient is given by Λ = ν, + 3ν,3 + ν 1,1 ν, + ν, ) which yields by computation γλ = α 3c 4 4c ) c ) + 1) α + c 1 ) = α p + 1) Y L. If L > resp. L < ) P is a weak sink resp. a weak source)..4.3 Energy and Lyapunov functionals for system 1.17) If xt) is a solution of.1) we set Then, if xt), yt)) Q, we have x F t t) = Lx t x t ) ) x x = L y + y q Hence, if L >, F is monotone in the region F t) := x t + x p + 1 Kx,.35) ) p xt = Lx t + y q ) x t ) y. { x, y) Q : x y > }, located under L. Remark. a) This function was introduced classically in the case =, leading in particular to.3) when p = +. b) With this function we can deduce an upper estimate of regular solution, completing, Proposition.3 and Proposition.4, namely if > and p + any ground state satisfies ) 1 p + 1 ur) K r for all r >..36) ext we construct a Lyapunov functional adapting the method initiated by [] and already used in [3] and [4]. Lemma.8 We define on R J x, y) = Kx x p + 1 ) p 3 p + 1) x 1 3p + 1 ) x y,.37) and, if xt), yt)) is a solution of.38), Vt) = J xt), yt)). If xt), yt)) = x, y) Q we have ) [ ) ) p )] x x x V t t) = y L y p y ) x ) p ) ).38) x = L y p y x y.

Quasilinear elliptic equations with source reaction 3 Consequently, the function t Vt) is decreasing if > and 1 < p +, and increasing if < and p +. Proof. We recall the ansatz introduced in [] for finding a Lyapunov function for a system of the form { xt = fx, y).39) y t = gx, y). If fx, y) = y = hx), then consider the function Lx, y) = y hx) fx, t)dt x gt, ht))dt..4) In the case of system 1.17), hx) = x, and we find Lx, y) = J x, y). Then.38) and the conclusion follow. Remark. If L > we set H = x, y) Q : ) p x y p x y = L ) 1, or x L = y = q q 1. Then R := θ 1 V t t) = xt), yt)) H { x, y) Q : oreover H is starshapped with respect to and we set θh = x, y) Q : } x = y. ) p x y p x y L ) 1, or x L = y q q 1. If > we have V t t) if xt), yt)) R and V t t) if xt), yt)) R c Q. If <, the signs of V t t) in the same regions are opposite. to In [7] we also used a function introduced in [19] for equation 1.). When q = p it reduces p + 1 Zr) = r a u r + u + a uu ) r r + u u r p,.41) where a = ) 1) p+3. Since r = e t, we find that in Q, Zr) = e )L p+3 t The function satisfies the relation ) p + 1 y + x axy + xy p..4) Z r p p + 1 u r Z = r a 1 U,.43)

Quasilinear elliptic equations with source reaction 4 where U = 1)) uu r p + 3) r = 1)) p + 3) xy pp + 3) L + L + ote that U has a constant sign in Q if L <. p + 1) ru r pp + 3) y p + 1) ) ) r p+3.44).4.4 Comparison results Lemma.9 Let 1, p > 1 and, R such that <. Then, as long as they lie below the line L, i.e. x t >, the regular trajectories T reg and T reg associated to and respectively do not intersect. Furthermore T reg is below T reg. Proof. We use the expansion.6) and deduce, for i = or, y i t) = e t 1 x i t) + i e t + p + o1)) )1 ).45) as t. Hence T reg is below T reg for t t for some t R. Suppose that the trajectories intersect for a first time at some point x, y ) below L. Since the system is autonomous, there will exist two solutions of the systems relative to and satisfying at the same time t, x t ) = x t ) = x and y t ) = y t ) = y. From.), x t t ) = x t ) > and p y t t ) = y t t ) + )y < y t t ). Hence the intersection of the two trajectories is transverse and the slope of T reg at this intersection point is smaller that the one of T reg which is impossible. The second assertion follows immediately. 3 Study of ground states of 1.17) when > When and either 3 and 1 < p, or = 1, and p > 1, there exists no ground state by Proposition.3; in this section we assume 3 and p > 3.1 Behaviour near equilibrium. Since and p > there exists a unique equilibrium P Q. The sign of the real part of the characteristic roots of the linearization of the system.) at P depends on the position of with respect to defined in.31). Lemma 3.1 Assume > and 3. 1) If < p +, then P is a source. ) If p > + ; then P is a sink when < and it is a source when >. 3) If p > +, then P a weak sink. Furthermore, if < < ɛ, for ɛ small enough, there exists a periodic trajectory in Q surrounding P. t

Quasilinear elliptic equations with source reaction 5 Proof. Step 1: Assume < p +. The linearized system is given in.8). Because >, the product of the characteristic roots given by equation.9) is positive since it is given by K p p + 1 Y = p + 1 K + p ) p Y p + 1. The sum or the real part) of the characteristic roots is equal to positive, as p + implies L <. Hence P is a source. Step : Assume p > + p p + 1 Y L which is, hence L >. As in Step 1, the product of the characteristic roots is positive. By Proposition.-1), Y is an increasing function of, then the sum or p the real part of) of the characteristic roots, given by p + 1 Y L is also increasing and vanishes if =. It is negative if < < and positive if >. It implies assertion. Step 3: If =, and p > +, then P is a weak sink by Lemma.7. The appearance of the limit cycle, which is the called the Andronov-Hopf bifurcation, occurs for > when is small enough see [17, Chap. 9]). Remark. The product of the characteristic roots is also expressed by K p p + 1 Y = p + 1 Y + ) Y. Hence it is positive for any, p > 1. ext we give some sufficient conditions for nonexistence of a periodic solution or a homoclinic orbit at, ). Lemma 3. Assume 3. If > and 1 < p + the system 1.17) admits no closed orbit in Q. If < and p > +, the system 1.17) admits no cycle in Q surrounding P. Proof. If γ is a non-trivial closed orbit it corresponds either to a T -periodic solution or a solution such that xt), yt)) = lim xt), yt)). lim t The function V defined in.37) is monotone and it satisfies either V) = VT ) or lim Vt) = t lim Vt). Hence it is constant and by.38) it implies yt) = t xt) for all t, a contradiction. ext we suppose p > +, < and that there exists a T -periodic solution xt), yt) with trajectory γ Q surrounding P, hence P belongs to the bounded connected component Γ of R \ {γ} bordered by γ. Since P is a sink or a weak sink by Lemma 3.1, there exists a neighbourhood O of P such that all the trajectories issued from O converge to P as t. Hence any trajectory issued from O, necessarily contained in Γ, has an α-limit set in Γ which is either a stationary point different from P, which is excluded, or a limit cycle t

Quasilinear elliptic equations with source reaction 6 {xt), yt)} t [,τ) := γ Γ τ is its period). This limit cycle is not stable, hence, by Floquet s theory τ H1 xt), yt)) + H ) xt), yt)) τ ) p dt = x y p + 1 y t) L dt 3.1) We perform the change of unknowns x = x X, ȳ = y Y used in Lemma.7 which leads to the system.33). The explicit value of the remaining term is h x, ȳ) = x + X ) p X p px x + ȳ + Y ) p Y = Φ x) + Ψȳ), p p Y where Φ and Ψ are defined accordingly. It is positive by convexity because >. Since from.33), τ τ xt)dt ȳt)dt =, and τ = px xt)dt + p τ τ τ p + 1 Y ȳt)dt K ȳt)dt + h x, ȳ)dt we derive p) p p Y + p ) τ p + 1 Y K ȳt)dt <. Using the equation.15) satisfied by Y it yields ) K ) τ p + 1 Y ȳt)dt <. By.11)-ii), K τ τ p + 1 Y >, hence ȳt)dt <, therefore τ 1 yt)dt < Y and by concavity, 1 τ yt)) 1 τ ) dt yt)dt < Y τ τ. 3.3) Combining 3.1) and 3.3), we obtain Since Y which contradicts 3.4). ȳ 3.) L < p p + 1 Y. 3.4) is decreasing by Proposition.-1-ii) we have for = p p + 1 Y L p p + 1 Y L, Remark. Up to changing the sense of variation of Vt), the proof of the first assertion shows that there exists no closed orbit in Q if < and p +. However the proof of the second assertion is not valid when <. The nonexistence of any periodic solution can also be proved when the equilibrium is a node i.e. the two characteristic values are real with the same sign).

Quasilinear elliptic equations with source reaction 7 Lemma 3.3 1- Let 3 and p > +. There exists a unique and explicit > such that for any, P is a repelling node, degenerate if =. If < <, P is a repelling spiraling point. If 3 1 and < <, P is an attracting spiraling point. - If 11 and + < p < 1 4 1 there exists a unique and explicit < 1 < such that if P is an attracting node if < < 1 and an attracting spiraling point if 1 < <. 3- If > or if < < 1, there exist no periodic trajectory in Q around P, neither no homoclinic trajectory in Q at, ) surrounding P. Proof. The characteristic values are real if and only if ) p D = p + 1 Y L 8 p =. K p + 1 Y L + + 1 p ) p + 1 Y ) p p + 1 Y L + 1 ) 3.5) Hence 3.5) is satisfied if one of the following conditions holds: i) ii) Y p + 1 ) L 1, p Y p + 1 ) 3.6) L + 1. p It is easy to check that for p > +, one has < L + 1 < qk. Since >, by Proposition.-1), the mapping Y [, ). Therefore there exists a unique such that Using.15) we find = p + 1 p is continuous, increasing and from [, ) into Y = p + 1 ) L + 1. 3.7) p Concerning the lower bound in i), there holds ) p K p + 1 p K + ) 1 p ) 1. 3.8) L 1 > p > 1 1 4 and 11. 3.9) Then we can define 1 by 1 Y 1 = p + 1 ) L 1. 3.1) p

Quasilinear elliptic equations with source reaction 8 which leads to ote that p Hence 1Y 1 1 = p + 1 p ) p = L 1 < L = p K p 1 p K + + ) 1 p Y ) 1. 3.11) L + 1 = p Y. 1 < <. 3.1) ext we prove 3) by adapting an argument introduced in [1] for quadratic systems. We return to system.33) that we write under the form { xt = a x ȳ 3.13) ȳ i = c x + dȳ + Φ x) + Ψȳ) which defines a, c, d with Φ and Ψ given by 3.13), and the trinomial.9) for characteristic values endows the form T λ) = λ d + a)λ + ad + c. In the range of values of, the discriminant D = d + a) 4ad + c) = d a) 4c is positive. We consider the intersection of a straight line l passing through P with equation ȳ = A x with a trajectory xt), ȳt)). Then { U = ȳt A x t = A + d a)a + c) x + Φ x) + Ψȳ) = A + d a)a + c) x + Φ x) + ΨA x). We can choose A such that A + d a)a + c = since D >. Since Φ and Ψ achieves positive values, we derive from the expression of h that U > for x. This proves that any closed orbit around P or passing by P can intersect l only one time which is a contradiction. 3. Existence or nonexistence of ground states Proposition 3.4 Let 3, < p + and >. Then there exists no radial ground state, and there exists a singular solutions ur) such that lim r ur) = X r lim r r ur) = c for some c >. Proof. Assume that T reg remains in Q, by Lemma.5 we have three possibilities: α) either T reg converges to P when t, which is impossible since P is a source, β) or T reg has a limit cycle at, and this is impossible by Lemma 3., γ) or T reg converges to, ) when t, hence it is a homoclinic orbit. The function V defined in.37) is decreasing by Lemma.8. Since V ) = V ) =, it is identically and and

Quasilinear elliptic equations with source reaction 9 so is V t. This implies that xt) yt) = for all t which is a contradiction. Hence T reg does not remain in Q. We denote by Ξ the connected region of Q bordered by the semi-axis {x, y) : x =, y > } and T reg. Since H is outward on the semi-axis, Ξ is negatively invariant. By Section 3.4.1, yt), ) is a saddle point, hence the stable trajectory T st = {xt), yt)} satisfies lim t xt) = which implies ur) cr for some c >. Its α-limit set αt st ) cannot be a limit cycle as we have seen it above. If it contains, ) it implies again that V, which is monotone, is equal to, hence V t and xt) yt), which is impossible. Hence αt st) contains P. Since P is a source it implies that T st converges to P when t. Proposition 3.5 Let 3 and p > +. Then for any < there exists a ground state u which satisfies u U at. Proof. If < < resp. = ), P is a sink resp. a weak sink). Suppose first that the trajectory T reg does not stay in Q, then it leaves Q at some point, y s ) with y s >. As a consequence, the stable trajectory T st at, ) remains in the negatively invariant region Ξ defined in the proof of Proposition 3.4. Since it cannot converge to P when t it admits a limit cycle surrounding P which contradicts Lemma 3.. Therefore T reg Q and, again using Lemma 3., either it converges to P when t and the proof is complete, or to, ) and T reg = T st is a homoclinic trajectory. The trace of the linearized system.3) at, ) is equal to K = L <. Therefore, from [17, Th. 9.3.3] the connection is attracting and the trajectories inside the bounded region T bordered by the homoclinic trajectory T reg spiral towards it when t. Hence any such trajectory inside T either has a limit cycle when t which is impossible by Lemma 3. or converges to P which is also impossible. Consequently there exists no homoclinic trajectory at, ) which ends the proof. ext we study the case where is large enough. We have already proved in [7] that for any p > 1, there exists =, p) > see introduction) such that if > there exists no ground state, radial or non-radial. In the radial case we have a more precise result. Proposition 3.6 Let 3 and p > +. Then for any there exists no ground state, but there exist singular solutions u which satisfy lim r ur) = X and lim r ur) = c >. r r Proof. Since >, P is a source by Lemma 3.1, and there exist no periodic orbit neither a homoclinic trajectory at, ) by Lemma 3.3. Thus T reg leaves Q through the semiaxis {, y) : y > } by Lemma.5. As in the proof of Proposition 3.5 the stable trajectory T st at, ) remains in the negatively invariant region Ξ already defined. Then it converges necessarily to P when t. ext we study the case < <. Theorem 3.7 Let 3, p > +. There exist two positive real numbers min and max such that < min max < such that, 1- For < < min there exist ground states ondulating around U when r and positive singular solutions ondulating around U on [, ).

Quasilinear elliptic equations with source reaction 3 - For = min and for = max there exist ground states u such that lim r r ur) = c >. 3- For max < < there exists no ground state and there exist singular solutions such that lim r ur) = X or turning around P when r, and lim r ur) = c >. r r Proof. Since is subject to vary, we put it in exponent in the different specific trajectories of the system. For < <, P is a source and the trajectories converging toward this point when t are spiralling. We have three possibilities: i) either T reg leaves Q at some point, y s ), y s >, ii) or T reg has a ω-limit set which is a periodic orbit surrounding P, iii) or T reg converges to, ) when t. If i) holds, then T st remains in the region Ξ := Ξ of Q bordered by T reg and the semi-axis {, y) : y > }. Then, either it converges to P when t, or it admits an unstable from outside) α-limit cycle. We denote by x st, y st ) the first backward intersection of T st with the straight line L. Furthermore T reg intersects L at some point necessarily unique) x reg, y reg). Because of the relative position of T st and T reg, there holds g) = x reg x st >. If ii) holds, then we claim that T st leaves Q at some point x s, ) with x s >. Indeed, if T st Q, it is bounded by Proposition.3. Hence the α-limit set is non-empty. It cannot be, ) since T st T reg and T st cannot converge to P or have a limit cycle around P since it would imply that T st T reg { }. Hence T st intersects L at some point x st, y st ) and T reg intersect the first time L and x reg, y reg). Because T reg lies in the region of Q bordered by T st and the semi-axis {x, ) : x > }, there holds If iii) holds then g) =. g) = x reg x st <. The function g defined for any, ) is continuous. We know that there exists ɛ > such that for any, + ɛ), T reg has a ω-limit cycle surrounding P, hence g) <. If =, T reg leaves Q at some point, y s ) with y s > from Proposition 3.6 and the intersection of T reg with the semi-axis {, y) : y > } is transverse at, y s ). By continuity with respect to the parameter, for any ɛ, ], for ɛ > small enough, T reg leaves Q at some point of the semi-axis. So we are in situation i) and g) >. Then, by continuity of g there exist min and max such that < min max < such that g min ) = g max ) =. If = min or = max, T st = T reg and the trajectory T reg is a homoclinic at, ). For < < min we are in situation ii), which proves 1. For > max we are in situation i) and either T st converges to P or has a α-limit cycle surrounding P when t. Theorems A and A follow from the previous results.