Saïd Benachour
email: benachou@iecn.u-nancy.fr

Grzegorz Karch
email: karch@math.uni.wroc.pl

Philippe Laurençot
email: laurenco@mip.ups-tlse.fr

Asymptotic profiles of solutions to viscous Hamilton-Jacobi equations

Keywords: 2000 Mathematics Subject Classification: 35K15, 35B40. 1 Diffusive Hamilton-Jacobi equation, self-similar large time behavior, Laplacian unilateral estimates Equation de Hamilton-Jacobi diffusive, comportement asymptotique autosimilaire, estimations unilatérales du Laplacien

The large time behavior of solutions to the Cauchy problem for the viscous Hamilton-Jacobi equation u t -∆u + |∇u| q = 0 is classified. If q > q c := (N + 2)/(N + 1), it is shown that non-negative solutions corresponding to integrable initial data converge in W 1,p (R N ) as t → ∞ toward a multiple of the fundamental solution for the heat equation for every p ∈ [1, ∞] (diffusion-dominated case). On the other hand, if 1 < q < q c , the large time asymptotics is given by the very singular self-similar solutions of the viscous Hamilton-Jacobi equation.

For non-positive and integrable solutions, the large time behavior of solutions is more complex. The case q ≥ 2 corresponds to the diffusion-dominated case. The diffusion profiles in the large time asymptotics appear also for q c < q < 2 provided suitable smallness assumptions are imposed on the initial data. Here, however, the most important result asserts that under some conditions on initial conditions and for 1 < q < 2, the large time behavior of solutions is given by the self-similar viscosity solutions to the non-viscous Hamilton-Jacobi equation z t + |∇z| q = 0 supplemented with the initial datum z(x, 0) = 0 if x = 0 and z(0, 0) < 0.

Résumé

Nous classifions le comportement asymptotique des solutions du problème de Cauchy pour l'équation de Hamilton-Jacobi avec diffusion u t -∆u+|∇u| q = 0. Si q > q c := (N +2)/(N +1), nous montrons que, lorsque t → ∞, les solutions intégrables et positives convergent dans W 1,p (R N ) vers un multiple de la solution fondamentale de l'équation de la chaleur pour tout p ∈ [1, ∞] (diffusion dominante). Ensuite, si 1 < q < q c , le comportement asymptotique est décrit par la solution très singulière auto-similaire de l'équation de Hamilton-Jacobi avec diffusion.

En ce qui concerne les solutions intégrables et négatives, la situation est plus complexe. Le terme de diffusion est de nouveau dominant si q ≥ 2, ainsi que lorsque q c < q < 2 pourvu que la donnée initiale soit suffisamment petite. Ensuite, pour 1 < q < 2, nous identifions une classe de données initiales pour laquelle le comportement asymptotique des solutions est donné par une solution de viscosité auto-similaire de l'équation de Hamilton-Jacobi z t + |∇z| q = 0 avec la condition initiale (non continue) z(x, 0) = 0 si x = 0 et z(0, 0) < 0.

Introduction

We investigate the large time behavior of integrable solutions to the Cauchy problem for the viscous Hamilton-Jacobi equation

u t -∆u + |∇u| q = 0 , x ∈ R N , t > 0 , (1.1) u(x, 0) = u 0 (x) , x ∈ R N , (1.2)
where q > 1. The dynamics of the solutions to (1.1)-(1.2) is governed by two competing effects, namely those resulting from the diffusive term -∆u and those corresponding to the "hyperbolic" nonlinearity |∇u| q . Our aim here is to figure out whether one of these two effects rules the large time behavior, according to the values of q and the initial data u 0 . Since the nonlinear term |∇u| q is non-negative, it acts as an absorption term for non-negative solutions and as a source term for non-positive solutions. We thus consider separately non-negative and non-positive solutions. Let us outline our main results now.

For non-negative initial data, it is already known that diffusion dominates the large time behavior for q > q c := (N + 2)/(N + 1) and that the nonlinear term only becomes effective for q < q c [START_REF] Amour | Global existence and decay for viscous Hamilton-Jacobi equations[END_REF][START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF][START_REF] Ben-Artzi | Decay of mass for a semilinear parabolic equation[END_REF][START_REF] Biler | Asymptotic properties of solutions of the viscous Hamilton-Jacobi equation[END_REF]. We obtain more precise information in Theorems 2.1 and 2.2 below. In particular, if q ∈ (1, q c ) and the initial datum decays sufficiently rapidly at infinity, there is a balance between the diffusive and hyperbolic effects: the solution u(t) behaves for large t like the very singular solution to (1.1), the existence and uniqueness of which have been established in [START_REF] Benachour | Very singular solutions to a nonlinear parabolic equation with absorption. I. Existence[END_REF][START_REF] Benachour | Very singular solutions to a nonlinear parabolic equation with absorption. II -Uniqueness[END_REF][START_REF] Qi | The self-similar profiles of generalized KPZ equation[END_REF].

For non-positive initial data, there are two critical exponents q = q c and q = 2, as already noticed in [START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF], and the picture is more complicated. More precisely, the diffusion governs the large time dynamics for any initial data if q ≥ 2 and for sufficiently small initial data if q ∈ (q c , 2), and we extend the result from [START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF]Proposition 2.2] in that case (cf. Theorem 2.3, below). On the other hand, when q ∈ (1, 2), we prove that, for sufficiently large initial data, the large time behavior is governed by the nonlinear reaction term. This fact is also true for any initial datum u 0 ≡ 0 if N ≤ 3 and q is sufficiently close to 1. We actually conjecture that the nonlinear reaction term always dominates in the large time for any non-zero initial datum as soon as q ∈ (1, q c ).

Let us finally mention that, when q ∈ (q c , 2), there is at least one (self-similar) solution for which there is a balance between the diffusive and hyperbolic effects for large times [START_REF] Ben-Artzi | The local theory for viscous Hamilton-Jacobi equations in Lebesgue spaces[END_REF].

Before stating more precisely our results, let us recall that for every initial datum u 0 ∈ W 1,∞ (R N ) the Cauchy problem (1.1)-(1.2) has a unique global-in-time solution which is classical for positive times, that is

u ∈ C(R N × [0, ∞)) ∩ C 2,1 (R N × (0, ∞)) .
In addition, this solution satisfies the estimates

u(t) ∞ ≤ u 0 ∞ and ∇u(t) ∞ ≤ ∇u 0 ∞ for all t > 0. (1.3)
Moreover, by the maximum principle, u 0 ≥ 0 implies that u ≥ 0 and u 0 ≤ 0 ensures that u ≤ 0. We refer the reader to [START_REF] Amour | Global existence and decay for viscous Hamilton-Jacobi equations[END_REF][START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF][START_REF] Gilding | The Cauchy problem for u t = ∆u+|∇u| q[END_REF] for the proofs of all these preliminary results. In addition, a detailed analysis of the well-posedness of (1.1)-(1.2) in the Lebesgue spaces L p (R N ) may be found in the recent paper [START_REF] Ben-Artzi | The local theory for viscous Hamilton-Jacobi equations in Lebesgue spaces[END_REF].

Notations. The notation to be used is mostly standard. For 1 ≤ p ≤ ∞, the L p -norm of a Lebesgue measurable real-valued function v defined on R N is denoted by v p . We will always denote by • X the norm of any other Banach space X used in this paper. Also, W 1,∞ (R N ) denotes the Sobolev space consisting of functions in L ∞ (R N ) whose first order generalized derivatives belong to L ∞ (R N ). The space of compactly supported and

C ∞ -smooth functions in R N is denoted by C ∞ c (R N ), and C 0 (R N ) is the set of continuous functions u such that lim R→∞ sup |x|≥R {|u(x)|} = 0 .
For a real number r, we denote by r + := max {r, 0} its positive part and by r -:= max {-r, 0} its negative part. The letter C will denote generic positive constants, which do not depend on t and may vary from line to line during computations. Throughout the paper, we use the critical exponent

q c := N + 2 N + 1 .

Results and comments

As already outlined, the large time behavior of solutions to (1.1)-(1.2) is determined not only by the exponent q of the nonlinear term |∇u| q but also by the sign, size, and shape of the initial conditions. In the present paper, we attempt to describe this variety of different asymptotics of solutions, imposing particular assumptions on initial data. In order to present our results in the most transparent form, we divide this section into subsections.

Non-negative initial conditions

In Theorems 2.1 and 2.2 below, we always assume that

u 0 is a non-negative function in L 1 (R N ) ∩ W 1,∞ (R N ) , u 0 ≡ 0 , (2.1)
and we denote by u = u(x, t) the corresponding non-negative solution of the Cauchy problem (1.1)-(1.2). In that case, we recall that t -→ u(t) 1 is a non-increasing function and that |∇u| belongs to L q (R N × (0, ∞)). In addition, [START_REF] Amour | Global existence and decay for viscous Hamilton-Jacobi equations[END_REF][START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF][START_REF] Ben-Artzi | Decay of mass for a semilinear parabolic equation[END_REF], for details). Since we would have I ∞ = u 0 1 > 0 for the linear heat equation, we thus say that diffusion dominates the large time behavior when I ∞ > 0, that is, when q > q c . We first consider the diffusion-dominated case.

I ∞ := lim t→∞ R N u(x, t) dx = R N u 0 (x) dx - ∞ 0 R N |∇u(x, s)| q dx ds (2.2) satisfies I ∞ > 0 if q > q c and I ∞ = 0 if q ≤ q c (cf.
Theorem 2.1 Suppose (2.1) and that q > q c . For every p ∈ [1, ∞],

lim t→∞ t (N/2)(1-1/p) u(t) -I ∞ G(t) p = 0 (2.3) and lim t→∞ t (N/2)(1-1/p)+1/2 ∇u(t) -I ∞ ∇G(t) p = 0. (2.4)
Here, G(x, t) = (4πt) -N/2 exp(-|x| 2 /(4t)) is the fundamental solution of the heat equation.

When p = 1, the relation (2.3) is proved in [START_REF] Biler | Asymptotic properties of solutions of the viscous Hamilton-Jacobi equation[END_REF] and Theorem 2.1 extends the convergence of u(t) towards a multiple of

G(t) to W 1,p (R N ), p ∈ [1, ∞].
Remark 2.1 Theorem 2.1 holds true when I ∞ = 0 (i.e. for q ≤ q c ) as well, but in that case, the relation (2.3) says only that u(t) p tends to 0 as t → ∞ faster than t -(N/2)(1-1/p) .

Our next theorem is devoted to the balance case 1 < q < q c when a particular self-similar solution of (1.1) appears in the large time asymptotics. Theorem 2.2 Suppose (2.1). Assume that q ∈ (1, q c ) and, moreover, that

ess lim |x|→∞ |x| a u 0 (x) = 0 with a = 2 -q q -1 . (2.5) For every p ∈ [1, ∞], lim t→∞ t (N/2)(1-1/p)+(a-N )/2 u(t) -W (t) p = 0 (2.6)
and lim t→∞ t (N/2)(1-1/p)+(a-N )/2+1/2 ∇u(t) -∇W (t) p = 0, (2.7)

where W (x, t) = t -a/2 W (xt -1/2 , 1) is the very singular self-similar solution to (1.1).

For the existence and uniqueness of the very singular solution to (1.1), we refer the reader to [START_REF] Benachour | Very singular solutions to a nonlinear parabolic equation with absorption. I. Existence[END_REF][START_REF] Benachour | Very singular solutions to a nonlinear parabolic equation with absorption. II -Uniqueness[END_REF][START_REF] Qi | The self-similar profiles of generalized KPZ equation[END_REF]. Notice also that the initial datum u 0 is integrable by assumption (2.5) since a > N for 1 < q < q c . Remark 2.2 In the critical case q = q c , it is also expected that u(t) converges towards a multiple of G(t) with a correction in the form of an extra logarithmic factor resulting from the absorption term. This conjecture is supported by what is already known for non-negative solutions to the Cauchy problem w t -∆w + w (N +2)/N = 0 (see, e.g., [START_REF] Vázquez | Asymptotic behaviour of nonlinear parabolic equations. Anomalous exponents[END_REF] and the references therein).

Non-positive initial conditions

We now turn to non-positive solutions and assume that

u 0 is a non-positive function in L 1 (R N ) ∩ W 1,∞ (R N ) , u 0 ≡ 0 . (2.8)
We denote by u = u(x, t) the corresponding non-positive solution of the Cauchy problem (1.1)-(1.2). In that case, we recall that t -→ u(t) 1 is a non-decreasing function and put

I ∞ := inf t≥0 R N u(x, t) dx = -sup t≥0 u(t) 1 ∈ [-∞, -u 0 1 ] . (2.9) Substituting u = -v in (1.1)-(1.2) we obtain that v = v(x, t) is a non-negative solution to v t -∆v -|∇v| q = 0, v(x, 0) = -u 0 (x), (2.10)
which has been studied in [START_REF] Ben-Artzi | The local theory for viscous Hamilton-Jacobi equations in Lebesgue spaces[END_REF][START_REF] Gilding | The Cauchy problem for the KPZ equation, prépublication LAMFA[END_REF][START_REF] Gilding | The Cauchy problem for u t = ∆u+|∇u| q[END_REF][START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF].

We start again with the diffusion-dominated case.

Theorem 2.3 Suppose (2.8).

a) Assume that q ≥ 2. Then I ∞ > -∞ and |∇u| belongs to L q (R N × (0, ∞)). In addition, I ∞ is given by (2.2) and the relations (2.3) and (2.4) hold true for every

p ∈ [1, ∞].
b) Assume that q ∈ (q c , 2). There exists ε = ε(N, q) such that, if

u 0 1 ∇u 0 (N +1)q-(N +2) ∞ < ε , (2.11)
then the conclusions of part a) are still valid.

The fact that I ∞ > -∞ under the assumptions of Theorem 2.3 is established in [START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF], together with the relation (2.3) for p = 1. We extend here this convergence to

W 1,p (R N ), p ∈ [1, ∞].
The smallness assumption imposed in (2.11) is necessary to obtain the heat kernel as the first term of the asymptotic expansion of solutions. This is an immediate consequence of the following theorem and the subsequent discussion.

Theorem 2.4 Suppose (2.8) and that q ∈ (q c , 2).

a) There exists a non-positive self-similar solution

V = V (x, t) = t -(2-q)/(2(q-1)) V (xt -1/2 , 1)
to (1.1) such that

lim t→∞ t (N/2)(1-1/p) V (t) p = ∞ and lim t→∞ t (N/2)(1-1/p)+1/2 ∇V (t) p = ∞ for all p ∈ [1, ∞]. b) There is a constant K = K(q) ≥ 0 such that, if u 0 ∈ W 2,∞ (R N ) satisfies u 0 ∞ (∆u 0 ) + 1-2/q ∞ > K (2.12) then lim t→∞ u(t) ∞ > 0. (2.13)
The first part of Theorem 2.4 is proved in [START_REF] Ben-Artzi | The local theory for viscous Hamilton-Jacobi equations in Lebesgue spaces[END_REF] while the second assertion is new. Let us point out here that, for the Hamilton-Jacobi equation w t + |∇w| q = 0, the L ∞ -norm of solutions remains constant throughout time evolution, while it decays to zero for the linear heat equation. We thus realize that, under the assumptions of Theorem 2.4 b), the diffusive term is not strong enough to drive the solution to zero in L ∞ as t → ∞ and the large time dynamics is therefore ruled by the Hamilton-Jacobi term |∇u| q .

Unfortunately, the conditions (2.11) and (2.12) do not involve the same quantities. Still, we can prove that if u 0 fulfils

u 0 ∞ D 2 u 0 1-2/q ∞ > K
(which clearly implies (2.12) since q < 2), the quantity u 0 1 ∇u 0 (N +1)q-(N +2) ∞ cannot be small. Indeed, there is a constant C depending only on q and N such that

u 0 ∞ D 2 u 0 1-2/q ∞ q(N +1)/2 ≤ C u 0 1 ∇u 0 q(N +1)-(N +2) ∞ . (2.14) For the proof of (2.14), put B = u 0 ∞ D 2 u 0 1-2/q ∞
and note that the Gagliardo-Nirenberg inequalities

u 0 ∞ ≤ C ∇u 0 N/(N +1) ∞ u 0 1/(N +1) 1 , ∇u 0 ∞ ≤ C D 2 u 0 (N +1)/(N +2) ∞ u 0 1/(N +2) 1 , imply that ∇u 0 (2-q)(N +2) ∞ ≤ C D 2 u 0 (2-q)(N +1) ∞ u 0 2-q 1 = C B -q(N +1) u 0 q(N +1) ∞ u 0 2-q 1 ≤ C B -q(N +1) ∇u 0 qN ∞ u 0 2 1 ,
whence the above claim.

We next show that the second assertion of Theorem 2.4 is also true when q ∈ (1, q c ).

Theorem 2.5 Suppose (2.8) and that q ∈ (1, q c ]. There is a constant

K = K(q) ≥ 0 such that, if u 0 ∈ W 2,∞ (R N ) fulfils (2.
12), then (2.13) holds true. Furthermore, if N ≤ 3 and 1 < q < 4/ 1 + √ 1 + 2N , then K(q) = 0.

We actually conjecture that K(q) = 0 for any q ∈ (1, q c ), but we have yet been unable to prove it.

The last result confirms the domination of the Hamilton-Jacobi term for large times when (2.13) holds true and provides precise information on the large time behavior.

Theorem 2.6 Let q ∈ (1, 2). Assume that u 0 ∈ C 0 (R N ) fulfils (2.8
) and is such that

M ∞ := lim t→∞ u(t) ∞ > 0 . (2.15) Then lim t→∞ u(t) -Z M∞ (t) ∞ = 0 , (2.16)
where Z M∞ is given by

Z M∞ (x, t) := -M ∞ -(q -1) q -q/(q-1)
|x| t 1/q q/(q-1) +

(2.17)

for (x, t) ∈ R N × (0, ∞). In fact, Z M∞ is the unique viscosity solution in BUC(R N × (0, ∞)) to z t + |∇z| q = 0 in R N × (0, ∞) (2.18)
with the bounded and lower semicontinuous initial datum z(x, 0) = 0 if x = 0 and z(0, 0) = -M ∞ .

The last assertion of Theorem 2.6 follows from [START_REF] Strömberg | The Hopf-Lax formula gives the unique viscosity solution[END_REF]. Moreover, Z M∞ is actually given by the Hopf-Lax formula

Z M∞ (x, t) = inf y∈R N -M ∞ 1 {0} (y) + (q -1) q -q/(q-1) |x -y| q/(q-1) t -1/(q-1) for (x, t) ∈ R N × (0, ∞), where 1 {0} denotes the characteristic function of the set {0}. Observe that Z M∞ is a self-similar solution to (2.18) since Z M∞ (x, t) = Z M∞ (xt -1/q , 1).
If N = 1, the convergence stated in Theorem 2.6 extends to the gradient of u. Proposition 2.1 Assume that N = 1 and consider a non-positive function u 0 in W 1,1 (R) ∩ W 1,∞ (R). Under the assumptions and notations of Theorem 2.6, we have also

lim t→∞ t (1-1/p)/q u x (t) -Z M∞,x (t) p = 0 for p ∈ [1, ∞).
In fact, if N = 1 and

u 0 ∈ W 1,1 (R) ∩ W 1,∞ (R), the function U := u x is a solution to the convection-diffusion equation U t -U xx + (|U| q ) x = 0 , x ∈ R , t > 0 , (2.19) with initial datum U(0) = u 0x and satisfies R U(x, t) dx = R u 0x (x) dx = 0 , t ≥ 0 . (2.20)
The large time behavior of non-negative or non-positive integrable solutions to (2. [START_REF] Karch | On zero mass solutions of viscous conservation laws[END_REF]) is now well-identified [START_REF] Escobedo | Large time behavior for convection-diffusion equations in R N[END_REF][START_REF] Escobedo | Asymptotic behavior and source-type solutions for a diffusion-convection equation[END_REF] but this is far from being the case for solutions satisfying (2.20). In this situation, some sufficient conditions on U(0) are given in [START_REF] Karch | On zero mass solutions of viscous conservation laws[END_REF] for the solution to (2.19) to exhibit a diffusion-dominated large time behavior. Also, convergence to N-waves is studied in [START_REF] Kim | An Oleinik type estimate for a convection-diffusion equation and convergence to N-waves[END_REF] but, for solutions satisfying (2.20), no condition is given in that paper which guarantees that U(t) really behaves as an N-wave for large times. As a consequence of our analysis, we specify such a condition and also provide several new information on the large time behavior of solutions to (2. [START_REF] Karch | On zero mass solutions of viscous conservation laws[END_REF]) satisfying (2.20). Results on the large time behavior of solutions to equation (2.19) satisfying the condition (2.20) are reviewed in the companion paper [START_REF] Benachour | Asymptotic profiles of solutions to convection-diffusion equations[END_REF].

We finally outline the contents of the paper: the next section is devoted to some preliminary estimates. Theorems 2.1 and 2.3 (diffusion-dominated case) are proved in Section 4 and Theorem 2.2 in Section 5. The remaining sections are devoted to the "hyperbolic"-dominated case: Theorems 2.4 and 2.5 are proved in Section 5 and Theorem 2.6 and Proposition 2.1 in Section 6.

Preliminary estimates

Let us first state a gradient estimate for solutions to (1.1) which is a consequence of [4, Theorem 1] (see also [START_REF] Gilding | The Cauchy problem for u t = ∆u+|∇u| q[END_REF]Theorem 2]). Note that, in this section, we do not impose a sign condition on the solution u to (1.1).

Proposition 3.1 Assume that u = u(x, t) is the solution to (1.1)-(1.2) corresponding to the initial datum u 0 ∈ W 1,∞ (R N ). For every q > 1, there is a constant C 1 > 0 depending only on q such that ∇u(t) ∞ ≤ C 1 u 0 1/q ∞ t -1/q , for all t > 0 . (3.1)
Proof. Setting v = u+ u 0 ∞ , it readily follows from (1.1) and the maximum principle that v is a non-negative solution to (1.1). By [4, Theorem 1], there is a constant C depending only on q such that ∇v (q-1)/q (t) ∞ ≤ C t -1/q , t > 0 .

Since ∇v = (q/(q -1)) v 1/q ∇v (q-1)/q and |u(x, t)| ≤ u 0 ∞ , we further deduce that

∇u(t) ∞ = ∇v(t) ∞ ≤ C v(t) 1/q ∞ ∇v (q-1)/q (t) ∞ ≤ C u 0 1/q ∞ t -1/q , whence (3.1).
Next, we derive estimates for the second derivatives of solutions to (1.1)-(1.2) when q ∈ (1, 2]. Proposition 3.2 Under the assumptions of Proposition 3.1, if q ∈ (1, 2], the Hessian matrix D 2 u = u x i x j 1≤i,j≤N of u satisfies

D 2 u(x, t) ≤ ∇u 0 2-q ∞ q (q -1) t Id , (3.2) D 2 u(x, t) ≤ C 2 u 0 (2-q)/q ∞ t 2/q Id , (3.3) for (x, t) ∈ R N × (0, ∞), where C 2 is a positive constant depending only on q. Furthermore, if u 0 ∈ W 2,∞ (R N ), D 2 u(x, t) ≤ D 2 u 0 ∞ Id . (3.4)
In Proposition 3.2, Id denotes the identity matrix of M N (R). Given two matrices

A and B in M N (R), we write A ≤ B if Aξ • ξ ≤ Bξ • ξ for every vector ξ ∈ R N .
For q = 2, the estimates (3.2) and (3.3) follow from the analysis of Hamilton [START_REF] Hamilton | A matrix Harnack estimate for the heat equation[END_REF] (since, if f is a non-negative solution to the linear heat equation f t = ∆f , the function -ln f solves (1.1) with q = 2). In Proposition 3.2 above, we extend that result to any q ∈ (1, 2].

Remark 3.1 The estimates (3.2) and (3.3) may also be seen as an extension to a multidimensional setting of a weak form of the Oleinik type gradient estimate for scalar conservation laws. Indeed, if N = 1 and

U = u x , then U is a solution to U t -U xx + (|U| q ) x = 0 in R × (0, ∞). The estimates (3.2) and (3.3) then read U x ≤ C U(0) 2-q ∞ t -1 and U x ≤ C u 0 (2-q)/q ∞ t -2/q
for t > 0, respectively, and we thus recover the results of [START_REF] Feireisl | The L 1 -stability of constant states of degenerate convection-diffusion equations[END_REF][START_REF] Kim | An Oleinik type estimate for a convection-diffusion equation and convergence to N-waves[END_REF] in that case.

Proof of Proposition 3.2. For 1 ≤ i, j ≤ N, we put w ij = u x i x j . It follows from equation (1.1) that

w ij,t -∆w ij = -q |∇u| q-2 N k=1 u x k w jk x i = -q |∇u| q-2 N k=1 w ik w jk -q |∇u| q-2 N k=1 u x k w jk,x i (3.5) -q (q -2) |∇u| q-4 N k=1 u x k w ik N k=1 u x k w jk . Consider now ξ ∈ R N \ {0} and set h = N i=1 N j=1 w ij ξ i ξ j .
Multiplying (3.5) by ξ i ξ j and summing up the resulting identities yield

h t -∆h = -q |∇u| q-2 N k=1 N i=1 w ik ξ i 2 -q |∇u| q-2 ∇u • ∇h -q (q -2) |∇u| q-4 N i=1 N j=1 u x j w ij ξ i 2 . (3.6)
Thanks to the following inequalities

|∇u| q-4 N i=1 N j=1 u x j w ij ξ i 2 ≤ |∇u| q-4 N j=1 |u x j | 2 N j=1 N i=1 w ij ξ i 2 ≤ |∇u| q-2 N k=1 N i=1 w ik ξ i 2 ,
and

h 2 ≤ |ξ| 2 N k=1 N i=1 w ik ξ i 2 ,
and since q ≤ 2, the right-hand side of identity (3.6) can be bounded from above. We thus obtain

h t -∆h ≤ -q (q -1) |∇u| q-2 N k=1 N i=1 w ik ξ i 2 -q |∇u| q-2 ∇u • ∇h ≤ -q |∇u| q-2 ∇u • ∇h - q (q -1) |∇u| q-2 |ξ| 2 h 2 .
Consequently,

Lh ≤ 0 in R N × (0, ∞) , (3.7)
where the parabolic differential operator L is given by

Lz := z t -∆z + q |∇u| q-2 ∇u • ∇z + q (q -1) |∇u| q-2 |ξ| 2 z 2 .
On the one hand, since q ∈ (1, 2] and |∇u(x, t)| ≤ ∇u 0 ∞ , it is straightforward to check that

H 1 (t) := 1 h(0) ∞ + q (q -1) t |ξ| 2 ∇u 0 2-q ∞ -1 , t > 0 , satisfies LH 1 ≥ 0 with H 1 (0) ≥ h(x, 0) for all x ∈ R N . The comparison principle then entails that h(x, t) ≤ H 1 (t) for (x, t) ∈ R N × (0, ∞), from which we conclude that h(x, t) ≤ h(0) ∞ ≤ D 2 u 0 ∞ |ξ| 2 ,
and

h(x, t) ≤ |ξ| 2 ∇u 0 2-q ∞ q (q -1) t .
In other words, (3.2) and (3.4) hold true.

On the other hand, we infer from (3.1) that

H 2 (t) := 2 C 2-q 1 |ξ| 2 q 2 (q -1) u 0 (2-q)/q ∞ t 2/q , t > 0 , satisfies LH 2 ≥ 0 with H 2 (0) = +∞ ≥ h(x, 0) for all x ∈ R N .
We then use again the comparison principle as above and obtain (3.3).

Remark 3.2 Since q ∈ (1, 2] and ∇u may vanish, the proof of Proposition 3.2 is somehow formal because of the negative powers of |∇u| in (3.6). It can be made rigorous by first considering the regularised equation

u ε t -∆u ε + |∇u ε | 2 + ε 2 p/2 =
0 for ε ∈ (0, 1), and then letting ε → 0 as in [START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF].

In fact, we need a particular case of Proposition 3.2.

Corollary 3.1 Under the assumptions of Proposition 3.2 ∆u(x, t) ≤ C 3 ∇u 0 2-q ∞ t , (3.8) ∆u(x, t) ≤ C 4 u 0 (2-q)/q ∞ t 2/q , (3.9)
for (x, t) ∈ R N × (0, +∞), where C 3 and C 4 are positive constants depending only on q and N.

Furthermore, if

u 0 ∈ W 2,∞ (R N ), sup x∈R N ∆u(x, t) ≤ sup x∈R N ∆u 0 (x) , t ≥ 0 . (3.10)
Proof. Consider i ∈ {1, . . . , N} and define ξ i = (ξ i j ) ∈ R N by ξ i i = 1 and ξ i j = 0 if j = i. We take ξ = ξ i in (3.7) and obtain that Lu x i x i ≤ 0, that is, (u x i x i ) t -∆u x i x i + q |∇u| q-2 ∇u.∇u x i x i + q (q -1) |∇u| q-2 u 2

x i x i ≤ 0 in R N × (0, ∞). Summing the above inequality over i ∈ {1, . . . , N} and recalling that

|∆u| 2 ≤ N N i=1 u 2 x i x i ,
we end up with (∆u) t -∆ (∆u) + q |∇u| q-2 ∇u.∇ (∆u) + q (q -1)

|∇u| q-2 N |∆u| 2 ≤ 0 in R N × (0, ∞).
We next proceed as in the proof of Proposition 3.2 to complete the proof of Corollary 3.1.

Diffusion-dominated case

The proofs of Theorems 2.1 and 2.3 rely on some properties of the non-homogeneous heat equation which we state now. Similar results have already been used in [START_REF] Biler | Asymptotic properties of solutions of the viscous Hamilton-Jacobi equation[END_REF][START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF].

Theorem 4.1 Assume that u = u(x, t) is the solution of the Cauchy problem to the linear non-homogeneous heat equation

u t = ∆u + f (x, t), x ∈ R N , t > 0, (4.1) u(x, 0) = u 0 (x), x ∈ R N , (4.2) with u 0 ∈ L 1 (R N ) and f ∈ L 1 (R N × (0, ∞)). Then lim t→∞ u(t) -I ∞ G(t) 1 = 0 , (4.3)
where

I ∞ := lim t→∞ R N u(x, t) dx = R N u 0 (x) dx + ∞ 0 R N f (x, t) dx dt .
Assume further that there is

p ∈ [1, ∞] such that f (t) ∈ L p (R N ) for every t > 0 and lim t→∞ t 1+(N/2)(1-1/p) f (t) p = 0 . (4.4) Then lim t→∞ t (N/2)(1-1/p) u(t) -I ∞ G(t) p = 0 , (4.5)
and lim t→∞ t (N/2)(1-1/p)+1/2 ∇u(t) -I ∞ ∇G(t) p = 0 . (4.6)

Proof. We first observe that the assumptions on u 0 and f warrant that I ∞ is finite, and we refer to [START_REF] Biler | Asymptotic properties of solutions of the viscous Hamilton-Jacobi equation[END_REF] for the proof of (4.3). We next assume (4.4) and prove (4.6). Let T > 0 and t ∈ (T, ∞). By the Duhamel formula,

∇u(t) = ∇G(t -T ) * u(T ) + t T ∇G(t -τ ) * f (τ ) dτ.
It follows from the Young inequality that

t (N/2)(1-1/p)+1/2 ∇u(t) -∇G(t -T ) * u(T ) p ≤ C t (N/2)(1-1/p)+1/2 (T +t)/2 T (t -τ ) -(N/2)(1-1/p)-1/2 f (τ ) 1 dτ + C t (N/2)(1-1/p)+1/2 t (T +t)/2 (t -τ ) -1/2 f (τ ) p dτ ≤ C t t -T (N/2)(1-1/p)+1/2 ∞ T f (τ ) 1 dτ + C sup τ ≥T τ (N/2)(1-1/p)+1 f (τ ) p t (T +t)/2 (t -τ ) -1/2 τ -1/2 dτ ≤ C t t -T (N/2)(1-1/p)+1/2 ∞ T f (τ ) 1 dτ + C sup τ ≥T τ (N/2)(1-1/p)+1 f (τ ) p .
Also, classical properties of the heat semigroup (see, e.g., [START_REF] Duoandikoetxea | Moments, masses de Dirac et décomposition de fonctions[END_REF]) ensure that

lim t→∞ t (N/2)(1-1/p)+1/2 ∇G(t -T ) * u(T ) - R N u(x, T ) dx ∇G(t -T ) p = 0 ,
and lim t→∞ t (N/2)(1-1/p)+1/2 ∇G(t -T ) -∇G(t) p = 0 for every p ∈ [1, ∞]. Since, by elementary calculations, we have

∇u(t) -I ∞ ∇G(t) p ≤ ∇u(t) -∇G(t -T ) * u(T ) p + ∇G(t -T ) * u(T ) - R N u(x, T ) dx ∇G(t -T ) p + R N u(x, T ) dx -I ∞ ∇G(t -T ) p + |I ∞ | ∇G(t -T ) -∇G(t) p ,
the previous relations imply that lim sup t→∞ t (N/2)(1-1/p)+1/2 ∇u(t) -

I ∞ ∇G(t) p ≤ C ∞ T f (τ ) 1 dτ + sup τ ≥T τ (N/2)(1-1/p)+1 f (τ ) p + R N u(x, T ) dx -I ∞ .
The above inequality being valid for any T > 0, we may let T → ∞ and conclude that (4.6) holds true. The assertion (4.5) then follows from (4.3) and (4.6) by the Gagliardo-Nirenberg inequality.

Proof of Theorem 2.1. Since u is non-negative, we infer from [4, Eq. ( 17)] that there is a constant C = C(q) such that ∇u (q-1)/q (t) ∞ ≤ C u(t/2) (q-1)/q ∞ t -1/2 , t > 0 .

Also, u is a subsolution to the linear heat equation and therefore satisfies

u(t) p ≤ G(t) * u 0 p ≤ C t -(N/2)(1-1/p) u 0 1 , t > 0 ,
for every p ∈ [1, ∞] by the comparison principle. Since ∇u = (q/(q -1)) u 1/q ∇u (q-1)/q , we obtain that

t (N/2)(1-1/p)+1 |∇u(t)| q p ≤ C t (N +2-q(N +1))/2 -→ t→∞ 0 for p ∈ [1, ∞],
because q > (N + 2)/(N + 1). Theorem 2.1 then readily follows by Theorem 4.1 with f (x, t) = -|∇u(x, t)| q .

Proof of Theorem 2.3, part a). Since q ≥ 2, we infer from [START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF] that I ∞ is finite and negative and that

∇u ∈ L q (R N × (0, ∞)) . (4.7) Setting b := ∇u 0 q-2 ∞ , it follows from (1.3) that u t -∆u ≥ -b |∇u| 2 in R N × (0, ∞).
The comparison principle then entails that u ≥ w, where w is the solution to

w t -∆w = -b |∇w| 2 , w(0) = u 0 .
The Hopf-Cole transformation h := e -bw -1 then implies that h solves

h t -∆h = 0 , h(0) = e -bu 0 -1 . Therefore, for t > 0, 0 ≤ -bw(x, t) ≤ h(x, t) ≤ h(t) ∞ ≤ C t -N/2 h(0) 1 ≤ C t -N/2 , since u 0 ∈ L 1 (R N ) ∩ L ∞ (R N ). Recalling that 0 ≥ u ≥ w, we end up with u(t) ∞ ≤ C t -N/2 , t > 0 . (4.8) It next follows from [17, Theorem 2] that ∇u(t) ∞ ≤ C u(t/2) ∞ t -1/2 , t > 0 ,
which, together with (4.8), yields

∇u(t) ∞ ≤ C t -(N +1)/2 , t > 0 . (4.9)
Recalling (1.3), we also have

∇u(t) ∞ ≤ C (1 + t) -(N +1)/2 , t ≥ 0 . (4.10) We next put A 1 (t) := sup τ ∈(0,t) τ 1/2 ∇u(τ ) 1 ,
which is finite by [START_REF] Ben-Artzi | The local theory for viscous Hamilton-Jacobi equations in Lebesgue spaces[END_REF]. Since q ≥ 2 and N ≥ 1, we infer from the Duhamel formula and (4.10) that, for α ∈ (0, 1/2),

t 1/2 ∇u(t) 1 ≤ C u 0 1 + C t 1/2 t 0 (t -τ ) -1/2 ∇u(τ ) q q dτ ≤ C + C t 1/2 t 0 (t -τ ) -1/2 (1 + τ ) -(q-1)(N +1)/2 ∇u(τ ) 1 dτ ≤ C + C t 1/2 t 0 (t -τ ) -1/2 (1 + τ ) -1 τ -1/2 A 1 (τ ) dτ ≤ C + C α -1/2 (1-α)t 0 (1 + τ ) -1 τ -1/2 A 1 (τ ) dτ + C t 1/2 A 1 (t) t (1-α)t (t -τ ) -1/2 2 2 + t τ -1/2 dτ ≤ C + C α -1/2 t 0 (1 + τ ) -1 τ -1/2 A 1 (τ ) dτ + C A 1 (t) 1 1-α (1 -τ ) -1/2 τ -1/2 dτ , whence 1 -C α 1/2 A 1 (t) ≤ C(α) 1 + t 0 (1 + τ ) -1 τ -1/2 A 1 (τ ) dτ .
Consequently, there is α 0 ∈ (0, 1/2) sufficiently small such that

A 1 (t) ≤ B 1 (t) := C(α 0 ) 1 + t 0 (1 + τ ) -1 τ -1/2 A 1 (τ ) dτ for t ≥ 0. Now, for t ≥ 0, dB 1 dt (t) = C(α 0 ) (1 + t) -1 t -1/2 A 1 (t) ≤ C(α 0 ) (1 + t) -1 t -1/2 B 1 (t) ,
from which we deduce that

A 1 (t) ≤ B 1 (t) ≤ B 1 (0) exp C(α 0 ) t 0 (1 + τ ) -1 τ -1/2 dτ ≤ C(α 0 ) .
We have thus proved that

∇u(t) 1 ≤ C t -1/2 , t > 0 . (4.11)
We finally infer from (4.9), (4.11) and the Hölder inequality that

t (N/2)(1-1/p)+1 |∇u(t)| q p ≤ C t (N +2-q(N +1))/2 -→ t→∞ for p ∈ [1, ∞],
and we conclude as in the proof of Theorem 2.1.

Proof of Theorem 2.3, part b). Since q ∈ (q c , 2), we obtain from [START_REF] Ph | On the growth of mass for a viscous Hamilton-Jacobi equation[END_REF] that there is ε > 0 such that, if u 0 fulfils (2.11), then I ∞ is finite and negative and there are C > 0 and δ > 0 such that ∇u(t) q q ≤ C t -1 (1 + t) -δ , t > 0 . (4.12)

In particular,

|∇u| q ∈ L 1 (R N × (0, ∞)) and lim t→∞ t |∇u(t)| q 1 = 0 . (4.13) We next claim that ∇u(t) ∞ ≤ C t -(N +1)/2 , t > 0 . (4.14)
Indeed, we fix r ∈ (q c , q) such that r < N/(N -1) and define s = r/(r -1) and a sequence (r i ) i≥0 by

r 0 = 1 q and r i+1 = (N + 1) r -(N + 2) 2r + q r r i , i ≥ 0 .
We now proceed by induction to show that, for each i ≥ 0, there is

K i ≥ 0 such that ∇u(t) ∞ ≤ K i t -(N +1)/2 + t -r i , t > 0 . (4.15)
Thanks to (3.1), the assertion (4.15) is true for i = 0. Assume next that (4.15) holds true for some i ≥ 0. We infer from (4.12), (4.15) and the Duhamel formula that

∇u(t) ∞ ≤ C u 0 1 t -(N +1)/2 + C t/2 0 (t -τ ) -(N +1)/2 ∇u(τ ) q q dτ + C t t/2 (t -τ ) -(N/2)(1-1/r)-1/2 ∇u(τ ) q sq dτ ≤ C t -(N +1)/2 u 0 1 + t/2 0 ∇u(τ ) q q dτ + C t t/2 (t -τ ) -(N/2)(1-1/r)-1/2 ∇u(τ ) q/r ∞ ∇u(τ ) q/s q dτ ≤ C t -(N +1)/2 + C I(t) ,
where

I(t) := t t/2 (t -τ ) -(N/2)(1-1/r)-1/2 τ -(N +1)/2 + τ -r i q/r τ -1/s dτ .
Since r < N/(N -1) and q > q c , we have

I(t) ≤ C t t/2 (t -τ ) -(N/2)(1-1/r)-1/2 τ -(q(N +1))/2r + τ -(qr i )/r τ -1/s dτ ≤ C t -((N +1)r-(N +2))/2r t -(q(N +1))/2r + t -(qr i )/r ≤ C t -(N +1)/2 t -((N +1)q-(N +2))/2r + t -r i+1 ≤ C t -(N +1)/2 + t -r i+1 for t ≥ 1. Consequently, for t ≥ 1, ∇u(t) ∞ ≤ K i+1 t -(N +1)/2 + t -r i+1 , while (1.
3) implies that the same inequality is valid for t ∈ [0, 1] for a possibly larger constant K i+1 . Thus (4.15) is true for i + 1, which completes the proof of (4.15). To obtain (4.14), it suffices to note that r i → ∞ since q > r. Now, owing to (4.13) and (4.14), we are in a position to apply Theorem 4.1 and conclude that (2.3) and (2.4) holds true for p = 1 and p = ∞. The general case p ∈ (1, ∞) then follows by the Hölder inequality.

Convergence towards very singular solutions

The goal of this section is to prove Theorem 2.2. Recall that we assume that 1 < q < q c and that u 0 is a non-negative and integrable function satisfying in addition ess lim |x|→∞ |x| a u 0 (x) = 0 , (5.1) with a = (2 -q)/(q -1) ∈ (N, ∞). We define R(u 0 ) := inf {R > 0 , |x| a u 0 (x) ≤ γ q a.e. in {|x| ≥ R}} , where γ q := (q -1) (q-2)/(q-1) (2 -q) -1 and observe that R(u 0 ) is finite by (5.1).

Denoting by u the corresponding solution to (1.1) and introducing

τ (u 0 ) := (N + 2) -q(N + 1) (N + 1)q -N 1-q R(u 0 ) 2 ,
we infer from [5, Lemma 2.2 & Proposition 2.4] that there is a constant C 1 depending only on N and q such that

t (a-N )/2 u(t) 1 + t a/2 u(t) ∞ + t (a+1)/2 ∇u(t) ∞ ≤ C 1 (5.2)
for each t > τ (u 0 ) and

u(x, t) ≤ Γ q (|x| -R(u 0 )) , t > 0 , |x| > R(u 0 ) . (5.3)
Here, Γ q is given by Γ q (r) = γ q r -a , r ∈ (0, ∞).

Let us observe at this point that decay estimates for ∇u(t) in L p can be deduced from (5.2) and the Duhamel formula. Lemma 5.1 For p ∈ [1, ∞], there is a constant C(p) depending only on N, q and p such that t ((a+1)p-N )/2p ∇u(t) L p ≤ C(p) for t > τ (u 0 ) . (5.4) Proof. Indeed, since u is non-negative, it follows from [4, Theorem 1] that ∇u (q-1)/q (t) ∞ ≤ C(q) t -1/q for t > 0, which, together with (5.2) and the Duhamel formula entails that, for t > τ (u 0 ),

∇u(t) 1 ≤ ∇G(t/2) * u(t/2) 1 + t t/2 ∇G(t -s) * |∇u| q 1 ds ≤ C t -1/2 u(t/2) 1 + C t t/2 (t -s) -1/2 ∇u (q-1)/q (s) q ∞ u(s) 1 ds ≤ C t -(a+1-N )/2 + C t t/2 (t -s) -1/2 s -(a+2-N )/2 ds ≤ C t -(a+1-N )/2 .
Interpolating between (5.2) and the above estimate yields (5.4).

In order to investigate the large time behavior of u, we use a rescaling method and introduce the sequence of rescaled solutions (u k ) k≥1 defined by

u k (x, t) = k a u(kx, k 2 t) , (x, t) ∈ R N × [0, ∞) , k ≥ 1 .
Lemma 5.2 For k ≥ 1, we have

t (a-N )/2 u k (t) 1 + t a/2 u k (t) ∞ + t (a+1)/2 ∇u k (t) ∞ ≤ C 1 (5.5) for t > τ k := τ (u 0 ) k -2 and u k (x, t) ≤ Γ q |x| - R(u 0 ) k for |x| > R(u 0 ) k and t > 0 . (5.6)
Proof. It is straightforward to check that, for each k ≥ 1, u k is the solution to (1.1) with initial datum u k (0) and satisfies estimates (5.5) and (5.6) as a consequence of (5.2) and (5.3).

We next use (1.1) and the non-negativity of u k to control the behavior of u k (x, t) for large x uniformly with respect to k. For k ≥ 1, t > 0 and R ≥ 0, we put

I k (R, t) := {|x|≥R} u k (x, t) dx + t 0 {|x|≥R} |∇u k (x, t)| q dxdt . (5.
I k (R, t) = 0 . (5.8) Proof. Let ̺ be a non-negative function in C ∞ (R N ) such that 0 ≤ ̺ ≤ 1 and ̺(x) = 0 if |x| ≤ 1 2 and ̺(x) = 1 if |x| ≥ 1 .
For R > 0 and x ∈ R N , we set ̺ R (x) = ̺(x/R). As u k is a non-negative solution to (1.1), we have

I k (R, t) ≤ u k (x, t) ̺ R (x) dx + t 0 |∇u k (x, s)| q ̺ R (x) dxds ≤ u k (x, 0) ̺ R (x) dx + t 0 u k (x, s) |∆̺ R (x)| dxds ≤ k a-N {|x|≥kR/2} u 0 (x) dx + |∆̺| ∞ R 2 t 0 {R/2≤|x|≤R}
u k (x, s) dxds . (5.9)

Owing to (5.1) and (5.6), we further obtain that, for R ≥ 1 + 4 R(u 0 ),

I k (R, t) ≤ k a-N {|x|≥kR/2} Γ q |x| 2 dx + |∆̺| ∞ R 2 t 0 {R/2≤|x|≤R} Γ q |x| - R(u 0 ) k dxds ≤ C R -(a-N ) + T |∆̺| ∞ R 2 {R/2≤|x|≤R} Γ q R 4 dx ≤ C(T, ̺) R -(a-N ) .
Lemma 5.3 then readily follows since a > N.

We finally study the behavior of u k for small times.

Lemma 5.4 Let r > 0. There is a positive constant C(r) depending only on q, N and r such that

{|x|≥r} u k (x, t) dx ≤ C(r) sup |x|≥kr/2
{|x| a u 0 (x)} + t (5.10)

for t > τ k and k ≥ 4 R(u 0 )/r.
Proof. We fix r > 0 and use the same notations as in the proof of Lemma 5.3. Thanks to the properties of ̺, we infer from (5.9) with R = r that, for t > τ k and

k ≥ 4 R(u 0 )/r, {|x|≥r} u k (x, t) dx ≤ u k (x, t) ̺ r (x) dx ≤ k a-N {|x|≥kr/2} u 0 (x) dx + |∆̺| ∞ r 2 t 0 {r/2≤|x|≤r} u k (x, s) dxds ≤ C(̺, r) sup |x|≥kr/2 {|x| a u 0 (x)} + t ,
where we have used (5.6) to obtain the last inequality.

Proof of Theorem 2.2. Owing to Lemma 5.2 and Lemma 5.3 we may proceed as in [START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF]Theorem 3] to prove that there are a subsequence of (u k ) (not relabeled) and a non-negative function

u ∞ ∈ C((0, ∞); L 1 (R N )) ∩ L q ((s, ∞) × R N )) ∩ L ∞ (s, ∞; W 1,∞ (R N )) satisfying u ∞ (t) = G(t -s) * u ∞ (s) - t s G(t -τ ) * |∇u ∞ (τ )| q dτ and lim k→∞ sup τ ∈[s,t] u k (τ ) -u ∞ (τ ) 1 = 0 (5.11)
for every s > 0 and t > s.

It remains to identify the behavior of u ∞ as t → 0. On the one hand, consider r > 0 and t > 0. Since τ k → 0 as k → ∞, we have t > τ k for k large enough and it follows from Lemma 5.4, (5.1) and (5.11)

that 0 ≤ {|x|≥r} u ∞ (x, t) dx ≤ C(r) t . Consequently, lim t→0 {|x|≥r} u ∞ (x, t) dx = 0 . (5.12)
On the other hand, consider M > 0 and set k M := M 1/(a-N ) . For k ≥ k M , we denote by v k the solution to (1.1) with initial datum v k (0) given by v k (x, 0

) := M k N u 0 (kx), x ∈ R N . Since a > N, we have v k (0) ≤ u k (0) for k ≥ k M and the comparison principle warrants that v k (x, t) ≤ u k (x, t) , (x, t) ∈ R N × [0, ∞) , k ≥ k M . (5.13)
We next observe that (v k (0)) converges narrowly towards (M u 0 1 ) δ as k → ∞ (δ denoting the Dirac mass at x = 0). We then proceed as in [START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF] to conclude that lim

k→∞ sup τ ∈[s,t] v k (τ ) -S M (τ ) 1 = 0
for every s > 0 and t > s, where S M denotes the unique non-negative solution to (1.1) with initial datum (M u 0 1 ) δ [START_REF] Benachour | Global solutions to viscous Hamilton-Jacobi equations with irregular initial data[END_REF]. Recalling (5.11) and (5.13), we realize that

S M (x, t) ≤ u ∞ (x, t) , (x, t) ∈ R N × (0, ∞) .
The above inequality being valid for any M > 0, it is then straightforward to deduce that lim

t→0 {|x|≤r} u ∞ (x, t) dx = ∞ . (5.14)
In other words, u ∞ is a very singular solution to (1.1) and the uniqueness of the very singular solution to (1.1) (cf. [START_REF] Benachour | Very singular solutions to a nonlinear parabolic equation with absorption. II -Uniqueness[END_REF][START_REF] Qi | The self-similar profiles of generalized KPZ equation[END_REF]) implies that u ∞ = W , where W is the very singular solution to (1.1), see Theorem 2.2. The uniqueness of the limit actually entails that the whole sequence (u k ) k≥1 converges towards W in C([s, t]; L 1 (R N )) for s > 0 and t > s. Expressed in terms of u, we have thus shown that

lim t→∞ t (a-N )/2 u(t) -W (t) 1 = 0 . (5.15)
Finally, it follows from (5.2), (5.15) and the Gagliardo-Nirenberg inequality that (2.6) holds true.

The last step of the proof is to obtain the convergence (2.7) for the gradients. Consider p ∈ [1, ∞], t > 0 and α ∈ (0, 1). By the Duhamel formula, we have

A p (t) := t ((a+1)p-N )/2p ∇(u -W )(t) L p ≤ t ((a+1)p-N )/2p ∇G((1 -α)t) * (u -W )(αt) L p + t ((a+1)p-N )/2p t αt ∇G(t -s) * (|∇u(s)| q -|∇W (s)| q ) L p ds ≤ C(α) t (a-N )/2 (u -W )(αt) 1 + C t ((a+1)p-N )/2p t αt (t -s) -1/2 s -1/2 ∇(u -W )(s) L p ds ,
where we have used the fact that

max { ∇u(s) ∞ , ∇W (s) ∞ } ≤ C s -(a+1)/2
by (5.2) and the properties of W in order to obtain the last inequality. Consequently, by the definition of A p (t) and the change of variables s → ts, we obtain

A p (t) ≤ C(α) t (a-N )/2 (u -W )(αt) 1 + C t ((a+1)p-N )/2p t αt (t -s) -1/2 s -1/2 s -((a+1)p-N )/2p A p (s) ds ≤ C(α) t (a-N )/2 (u -W )(αt) 1 + C 1 α (1 -s) -1/2 s -1/2 s -((a+1)p-N )/2p A p (st) ds . Now, introducing A p (∞) := lim sup t→+∞ A p (t) ≥ 0 ,
which is finite by (5.4), we may let t → +∞ in the above inequality and use (5.15) to conclude that

A p (∞) ≤ C 1 α (1 -s) -1/2 s -1/2 s -((a+1)p-N )/2p ds A p (∞) .
Finally, the choice of α < 1 sufficiently close to 1 readily yields that A p (∞) = 0, from which (2.7) follows.

6 Proofs of Theorems 2.4 and 2.5

Proof of Theorem 2.4, part a). The required non-positive self-similar solution

V = V (x, t) = t -(2-q)/(2(q-1)) V x t -1/2 , 1
is constructed and studied in [START_REF] Ben-Artzi | The local theory for viscous Hamilton-Jacobi equations in Lebesgue spaces[END_REF]Theorem 3.5]. In particular, it is shown that the selfsimilar profile V(x) := V (x, 1) is a radially symmetric bounded C 2 function. Moreover, the profile V and its first derivative V ′ both decay exponentially as |x| → ∞ (see [

m(t) ≤ -u 0 ∞ + T (∆u 0 ) + ∞ + C u 0 (2-q)/q ∞ t T τ -2/q dτ
for T > 0 and t > T . Since q < 2, we may let t → ∞ in the above inequality and obtain with the choice T = u 0 (2-q)/2 ∞ (∆u 0 ) + -q/2 ∞ that there is a constant K depending only on q such that

m ∞ ≤ -u 0 ∞ + K q/2 (∆u 0 ) + (2-q)/2 ∞ u 0 (2-q)/2 ∞ . (6.1)
Therefore, if u 0 ∞ > K (∆u 0 ) + (2-q)/q ∞ , we readily conclude from (6.1) that m ∞ < 0, whence (2.13).

Proof of Theorem 2.5. The proof of the first assertion of Theorem 2.5 is the same as that of Theorem 2.4, part b), hence we skip it. We next assume that N ≤ 3 and that 1 < q < 4/(1 + √ 1 + 2N). For t > 0, we put

ℓ(t) := u(t) ∞ (∆u(t)) + 1-2/q ∞ .
Since u is a non-positive subsolution to the linear heat equation, we infer from classical properties of the heat semigroup that

u(t) ∞ ≥ G(t) * u 0 ∞ ≥ C t -N/2
for t large enough. As q < 2, this estimate and (3.9) entail that, for t large enough,

ℓ(t) ≥ C t (4(2-q)-N q 2 )/2q 2 -→ t→∞ ∞ , since q < 4/(1 + √ 1 + 2N ).
Consequently, there exists t 0 large enough such that ℓ(t 0 ) > K(q) and we may apply the first assertion of Theorem 2.5 to t -→ u(t 0 + t) to complete the proof.

Under the assumptions of Theorem 2.4, part b) or Theorem 2.5, we may actually bound the L 1 -norm of u(t) from below and improve significantly [21, Proposition 2.1]. Proposition 6.1 Assume that u 0 satisfies (2.8) and that

M ∞ := lim t→∞ u(t) ∞ > 0 .
Then there is a constant C = C(N, q, u 0 ) such that

u(t) 1 ≥ C t N/q , t ≥ 0 . (6.2) Proof. We fix t > 0. For k ≥ 1, let x k ∈ R N be such that u(t) ∞ -1/k ≤ -u(x k , t).
For R > 0, it follows from (3.1) and the time monotonicity of u(t) ∞ that

u(t) 1 ≥ - {|x-x k |≤R} u(x, t) dx ≥ {|x-x k |≤R} (-u(x k , t) -|x -x k | ∇u(t) ∞ ) dx ≥ C R N N u(t) ∞ - 1 k - C 1 R N +1 N + 1 u 0 1/q ∞ t -1/q ≥ C R N M ∞ - 1 k -C ′ R t -1/q .
Letting k → ∞ and choosing R = M ∞ t 1/q /(2 C ′ ) yields the claim (6.2).

7 Proof of Theorem 2.6 and Proposition 2.1

Proof of Theorem 2.6. STEP 1. Recall that, by (2.8), u 0 is a non-positive function. We assume further that u 0 is compactly supported in a ball B(0, R 0 ) of R N for some R 0 > 0.

For λ ≥ 1, we introduce

u λ (x, t) := u(λx, λ q t) , (x, t) ∈ R N × (0, ∞) , which solves u λ,t + |∇u λ | q = λ q-2 ∆u λ in R N × (0, ∞) (7.1)
with initial datum u λ (0). Lemma 7.1 There is a constant C = C(N, q, u 0 ∞ ) such that, for t ≥ 0 and λ ≥ 1,

u λ (t) ∞ + t 1/q ∇u λ (t) ∞ + t u λ,t (t) ∞ ≤ C . (7.2) Proof. It first follows from (1.3) that u λ (t) ∞ = u(λ q t) ∞ ≤ u 0 ∞ , while Proposition 3.1 yields ∇u λ (t) ∞ = λ ∇u(λ q t) ∞ ≤ C 1 u 0 1/q ∞ t -1/q .
We next infer from [START_REF] Gilding | The Cauchy problem for the KPZ equation, prépublication LAMFA[END_REF]Theorem 5] that

u λ,t (t) ∞ = λ q u t (λ q t) ∞ ≤ λ q C(N, q) u 0 ∞ (λ q t) -1 = C(N, q) u 0 ∞ t -1 ,
which completes the proof.

Owing to Lemma 7.1, we may apply the Arzelà-Ascoli theorem and deduce that there are a subsequence of (u λ ) (not relabeled) and a non-positive function u

∞ ∈ C(R N × (0, ∞)) such that u λ -→ u ∞ in C(B(0, R) × (t 1 , t 2 )) (7.3)
for any R > 0 and 0 < t 1 < t 2 . It also follows from (7.3) and Lemma 7.1 that u ∞ (t) ∈ BUC(R N ) and satisfies

u ∞ (t) ∞ + t 1/q ∇u ∞ (t) ∞ + t u ∞,t (t) ∞ ≤ C (7.4)
for each t > 0. We next introduce the function H λ : R × R N × S N (R) → R defined by H λ (ξ 0 , ξ, S) := ξ 0 + |ξ| q -λ q-2 tr(S) , where S N (R) denotes the subset of symmetric matrices of M N (R) and tr(S) denotes the trace of the matrix S. On the one hand, we notice that (7.1) reads

H λ (u λ,t , ∇u λ , D 2 u λ ) = 0 in R N × (0, ∞)
and that H λ is elliptic. On the other hand, H λ converges uniformly on every compact subset

of R × R N × S N (R) towards H ∞ : R × R N → R given by H ∞ (ξ 0 , ξ) := ξ 0 + |ξ| q .
Therefore, for every τ > 0, u ∞ (.+τ ) is the unique viscosity solution to (2.18) with initial datum u ∞ (τ ) ( see, e.g., [10, Proposition IV.1] and [9, Theorem 4.1]). In addition, since u ∞ (τ ) is bounded and Lipschitz continuous by (7.4), we infer from [14, Section 10.3, Theorem 3] that u ∞ (. + τ ) is given by the Hopf-Lax formula u ∞ (x, t + τ ) = inf y∈R N u ∞ (y, τ ) + (q -1) q -q/(q-1) |x -y| q/(q-1) t -1/(q-1) (7.5)

for (x, t) ∈ R N × [0, ∞).
It remains to identify the behavior of u ∞ (t) as t → 0. Consider first x ∈ R N , t ∈ (0, ∞) and s ∈ (0, t). We infer from (3.9) and (7.1) that

u λ (x, t) ≤ u λ (x, s) + λ q-2 t s ∆u λ (x, σ) dσ ≤ u λ (x, s) + λ q-2 t s λ 2 C (λ q σ) -2/q dσ ≤ u λ (x, s) -C λ q-2
t (q-2)/q -s (q-2)/q . Since q ∈ (1, 2), we may pass to the limit as λ → ∞ in the previous inequality and use (7.3) to deduce that t → u ∞ (x, t) is non-increasing for every x ∈ R N . Since u ∞ is bounded by (7.4), we may thus define u ∞ (0) by for n large enough. Letting n → ∞ then completes the proof of Theorem 2.6.

Proof of Lemma 7.2. Let h 0 be a non-positive function in C ∞ c (R) such that h 0 (y) = -u 0 ∞ if y ∈ (-R 0 , R 0 ) (recall that u 0 is compactly supported in B(0, R 0 )). We denote by h the solution to the one-dimensional viscous Hamilton-Jacobi equation h t -h yy + |h y | q = 0 in R × (0, ∞) , h(0) = h 0 in R .

For i ∈ {1, . . . , N} and (x, t) ∈ R N × (0, ∞), we put h i (x, t) := h(x i , t) and notice that h i is the solution to (1.1) with initial datum h i (0) ≤ u 0 . The comparison principle then entails that h(x i , t) = h i (x, t) ≤ u(x, t) ≤ 0 , (x, t) ∈ R N × (0, ∞) . (7.8)

We next introduce w := h y and notice that w is the solution to the one-dimensional convection-diffusion equation w t -w yy + (|w| q ) y = 0 in R × (0, ∞) , (7.9) w(0) = w 0 := h 0,y in R .

The comparison principle then entails that b(y, t) ≤ w(y, t) ≤ a(y, t) , (y, t) ∈ R × (0, ∞) , (7.10) where b ≤ 0 and a ≥ 0 denote the solutions to (7.9) with initial data b(0) = -w - 0 ≤ 0 and a(0) = w + 0 ≥ 0. Since w 0 ∈ L 1 (R), it follows from [START_REF] Escobedo | Asymptotic behavior and source-type solutions for a diffusion-convection equation[END_REF] that lim t→∞ b(t) -Σ -B (t) 1 = lim t→∞ a(t) -Σ A (t) 1 = 0 , (7.11) where B := b(0) 1 , A := a(0) 1 , and, for M ∈ R, Σ M is the source solution to the one-dimensional conservation law Σ M,t + (|Σ M | q ) y = 0 in R × (0, ∞) , Σ(0) = M δ 0 in R .

Here, δ 0 denotes the Dirac mass in R centered at y = 0. The source solution Σ M is actually given by Σ M (y, t) := y 1/(q-1) (qt) -1/(q-1) 1 [0,ξ M (t)] (y) , ξ M (t) := q M q -1 (q-1)/q t 1/q , if M ≥ 0, and Σ M (y, t) := -|y| 1/(q-1) (qt) -1/(q-1) 1 [-η M (t),0] (y) , η M (t) := q -M q -1 (q-1)/q t 1/q , if M ≤ 0 (see, e.g., [START_REF] Liu | Source-solutions and asymptotic behavior in conservation laws[END_REF]). In particular, Σ M satisfies λ Σ M (λy, λ q t) = Σ M (y, t) for (λ, y, t) ∈ (0, ∞) × R × (0, ∞) . (7.12) Now, let t > 0 and set ̺(t) := N 1/2 max {ξ A (t), η -B (t)} ≤ C t 1/q . If x ∈ R N is such that |x| > ̺(t), there is i ∈ {1, . . . , N} such that |x i | > max {ξ A (t), η -B (t)}, whence either x i > ξ A (t) or x i < -η -B (t). In the latter case, we infer from (7.8), (7.10) and (7.12) that 0 ≥ u λ (x, t) ≥ h(λx i , λ q t) = λx i -∞ w(y ′ , λ q t) dy

′ ≥ λ x i -∞ b(λy ′ , λ q t) dy ′ ≥ λ x i -∞
(b(λy ′ , λ q t) -Σ -B (λy ′ , λ q t)) dy ′ ≥ -(b -Σ -B )(λ q t) 1 .

Similarly, if x i > ξ A (t), (7.8), (7.10) and (7.12) yield 0 ≥ u λ (x, t) ≥ -(a -Σ A )(λ q t) 1 .

Therefore, if x ∈ R N is such that |x| > ̺(t), then |u λ (x, t)| ≤ max { (a -Σ A )(λ q t) 1 , (b -Σ -B )(λ q t) 1 } . (7.13) Passing to the limit as λ → ∞ in (7.13) and using (7.3) and (7.11) provide the first assertion of Lemma 7.2. We next use once more (7.3) and (7.13) to conclude that (7.7) holds true.

Proof of Proposition 2.1. We keep the notations of the proof of Theorem 2.6 and introduce U λ (x, t) := u λ,x (x, t) = λ u x (λx, λ q t) , (x, t) ∈ R × (0, ∞) .

It follows from (7.1) and Lemma 7.1 that U λ,t + (|U λ | q ) x = λ q-2 U λ,xx , (x, t) ∈ R × (0, ∞) , and U λ (t) 1 ≤ u 0,x 1 and t 1/q U λ (t) ∞ ≤ C (7.14) for t > 0. We recall that, by Theorem 2.6, the family (u λ ) converges towards Z M∞ in C(R N × [t 1 , t 2 ]) for any t 2 > t 1 > 0. Owing to (7.14), we readily conclude that (U λ ) converges weakly-⋆ towards Z M∞,x in L ∞ (R N × (t 1 , t 2 )) for any t 2 > t 1 > 0. We may then proceed along the lines of [START_REF] Escobedo | Asymptotic behavior and source-type solutions for a diffusion-convection equation[END_REF]Section 3] to show that (U λ ) converges towards Z M∞,x in L 1 (R) as λ → ∞. Expressing this convergence result in terms of U = u x and using (3.1) yield Proposition 2.1 by interpolation.

7 ) 5 . 3

 753 Lemma For every T > 0, we have lim

u

  ∞ (x, 0) := sup t>0 {u ∞ (x, t)} ∈ (-∞, 0] for x ∈ R N . ) -Z M∞ (t) ∞ ≤ 2 u n 0 -u 0 ∞

Acknowledgements. We thank Professor Herbert Koch for pointing out to us Ref. [18] and Professor Brian Gilding for useful comments on Proposition 3.2. The preparation of this paper was partially supported by the KBN grant 2 P03A 002 24, the POLONIUM project ÉGIDE-KBN No. 05643SE, and the EU contract HYKE No. HPRN-CT-2002-00282.

In particular, u ∞ (0) is a lower semicontinuous function as the supremum of continuous functions.

More information on u ∞ (0) are consequences of the next result.

Lemma 7.2 For each t > 0, there is ̺(t) > 0 such that u ∞ (x, t) = 0 if |x| > ̺(t) and

Moreover, ̺(t) → 0 as t → 0. Taking Lemma 7.2 for granted, we see that (7.6) and Lemma 7.2 imply that u ∞ (x, 0) = 0 for x = 0 since ̺(t) → 0 as t → 0. We set ℓ := -u ∞ (0, 0), so that

and fix (x, t) ∈ R N × (0, ∞). We will now proceed along the lines of [START_REF] Strömberg | The Hopf-Lax formula gives the unique viscosity solution[END_REF] to show that u ∞ (x, t) = Z ℓ (x, t) (recall that Z ℓ is defined in (2.17)). Introducing the notation µ := (q -1) q -q/(q-1) t -1/(q-1) , it follows from (7.6) and Lemma 7.2 that, for 0 < σ < τ

|x -y| q/(q-1) -|x| q/(q-1) , while, for 0 < σ < τ and |y| ≥ ̺(σ), u ∞ (y, σ) + µ |x -y| q/(q-1) ≥ 0 .

The previous bounds from below and (7.5) entail that

, we may pass to the limit as σ → 0 in the above inequality and obtain

for τ > 0. Letting τ → 0 yields u ∞ (x, t) ≥ min 0, -ℓ + µ |x| q/(q-1) = Z ℓ (x, t) .

On the other hand, (7.5) and (7.6) ensure that

whence u ∞ (x, t) ≤ Z ℓ (x, t) by the continuity of u ∞ in R N × (0, ∞). We have thus shown that u ∞ = Z ℓ . In particular, u ∞ (t) ∞ = ℓ for t ≥ 0. But (2.15) and (7.7)

whence ℓ = M ∞ and u ∞ = Z M∞ . For t > 0, the sequence (u λ (t)) has thus only one possible cluster point in L ∞ (R N ) as λ → ∞, from which we conclude that the whole family (u λ (t)) converges to

Setting λ = t 1/q and using the self-similarity of Z M∞ , we are finally led to (2.16).

STEP 2. We now consider an arbitrary function u 0 ∈ C 0 (R N ) fulfilling (2.8) and such that (2.15) holds true. There is a sequence (u n 0 ) of non-positive functions in

For n ≥ 1, we denote by u n the solution to (1.1) with initial datum u n 0 and put

By [START_REF] Gilding | The Cauchy problem for u t = ∆u+|∇u| q[END_REF]Corollary 4.3], we have

from which we readily deduce that

Consequently, M n ∞ -→ M ∞ as n → ∞ and (2.15) guarantees that M n ∞ > 0 for n large enough. The analysis performed in the previous step then implies that

for n large enough. Therefore,