On the hardy-Hénon heat equation with an inverse square potential
Abstract.
We study Cauchy problem for the Hardy-Hénon parabolic equation with an inverse square potential, namely,
where , and or , . We establish sharp fixed time-time decay estimates for heat semigroups in weighted Lebesgue spaces, which is of independent interest. As an application, we establish:
-
•
Local well-posedness (LWP) in scale subcritical and critical weighted Lebesgue spaces.
-
•
Small data global existence in critical weighted Lebesgue spaces.
-
•
Under certain conditions on and we show that local solution cannot be extended to global one for certain initial data in the subcritical regime. Thus, finite time blow-up in the subcritical Lebesgue space norm is exhibited.
-
•
We also demonstrate nonexistence of local positive weak solution (and hence failure of LWP) in supercritical case for the Fujita exponent. This indicates that subcriticality or criticality are necessary in the first point above.
In summary, we establish a sharp dissipative estimate and addresses short and long time behaviors of solutions. In particular, we complement several classical results and shed new light on the dynamics of the considered equation.
Key words and phrases:
Hardy-Hénon equation, Inverse square potential, Dissipative estimate, well-posesness, finite time blow up2000 Mathematics Subject Classification:
35K05, 35K67, 35B30, 35K571. Introduction
1.1. Fixed-time estimates for heat semigroup
Consider the linear heat equation associated with the inverse square potential, namely
(1.1) |
where In this paper, we assume that unless it is explicitly specified. The Schrödinger operator with inverse square potentials
is initially defined with domain . Then it is extended as an unbounded operator in weighted Lebesgue space that generates a positive semigroup provided and where , defined by
(1.2) |
are the roots of , see [24, Theorems 3.2, 3.3]. Here, the weighted Lebesgue space is defined by the norm
The study of is motivated from physics and mathematics spanning areas such as combustion theory, the Dirac equation with Coulomb potential, quantum mechanics and the study of perturbations of classic space-time metrics. See e.g. [30, 21, 8] and the references therein.
The aim of this article is to understand the dynamics of solutions of Hardy-Hénon heat equations (1.1) and (1.6) when a singular potential is present, in light of the research programme initiated by Zhang [32], Pinsky [27, 28], Ioku et al. in [17, 18], Ishige [19] and Ishige-Kawakami in [20], and Bhimani-Haque [3] (cf. [5, 6, 4]). We also note that there is a extensive literature on Hardy-Hénon heat equation without potential, i.e. (1.6) with we refer to recent work of Chikami et al. in [10, 9] and the references therein, see also Remark 1.1.
We begin by stating our dissipative estimates in weighted Lebesgue spaces in the following theorem.
Theorem 1.1.
1.2. Hardy-Hénon equations (HHE) with inverse-square potential
We consider (1.1) with an inhomogeneous power type nonlinearity:
(1.6) |
where , and and or We assume that the non-linearity function satisfies the following conditions:
(1.7) |
The typical examples of would be
The potential is called Hénon type if and is called Sobolev type if . The equation (1.6) with is known as a Hardy parabolic equation, while that with is known as a Hénon parabolic equation. Equation (1.6) is called Hardy-Hénon parabolic equation with an inverse square potential. The elliptic part of (1.6) when , i.e.
was proposed by Hénon [14] as a model to study the rotating stellar systems and has been extensively studied in scientific community, see e.g. [12].
The equation (1.6) is invariant under the following scale transformation:
More precisely, if is a solution to (1.6), then so is with the rescaled initial data . Then the following identity holds
Hence, if and satisfy
then the identity holds for any , i.e., the norm is invariant with respect to . Denote
(1.8) |
We say Cauchy problem (1.6) scale
(1.9) |
Remark 1.2.
For we get (often denoted by for shorthand). In particular, when , we have So is the critical Lebesgue space without weight.
We recall the notion of well-posedness in the sense of Hadamard.
Definition 1.1 (well-posedness).
Let and
-
–
We say that is an -integral solution on to (1.6) if and satisfies
(1.10) for any . Maximum of such is denoted by .
-
–
Let be Banach spaces. Then (1.6) is called locally well-posed (in short LWP) from to if, for each bounded , there exist and a Banach space so that
-
(a)
for all , (1.6) has a unique integral solution
-
(b)
is continuous from to
If we say (1.6) is locally well-posed in . If , then we say (1.6) is globally well-posed in .
-
(a)
Remark 1.3.
We briefly mention some history on several facets of (1.6). We define Fujita exponent by
which is often known to divide the existence and nonexistence of positive global solutions.
- (1)
-
(2)
Taking in (1.6), we get classical Hardy-Hénon heat equation
(1.12) In this case, Chikami et al. in [9] introduced weighted Lebesuge space to treat potential and establish well-posedness results. Later, Chikami et al. in [10] generalize these results in weighted Lorentz spaces. In this paper, we could establish analogue of these results in the presence of potential, i.e. for (1.6) with and relaxed conditions on other parameters . See Remarks 1.4 and 1.7 below.
-
(3)
Several authors considered (1.6) with some mild restriction on external potential:
(1.13) and showed sharp contrast between existence of classical global solution and finite time blow-up in norm by finding appropriate Fujita exponent. We recall some of them here:
- (a)
- (b)
1.3. Dynamics of HHE with inverse square potential
We are now ready to state our well-posedness result in the following theorem.
Theorem 1.2 (Well-posedness: subcritical and critical case).
Let and be as defined in (1.2). Let
(1.14) |
and satisfies
(1.15) |
Let , and be as in (1.8) and satisfy
(1.16) |
Then Cauchy problem (1.6) is locally well-posed in and for the critical case we also have small data global existence. In the subcritical case, if we impose further restriction
(1.17) |
then one has uniqueness in .
Theorem 1.2 is new for and Up to now, we could not know the well-posedness of (1.6) with in the mere spaces but in weighted spaces. See Remark 1.1(6). We prove Theorem 1.2 via fixed point argument. To this end, the main new ingredient required is our fixed-time estimate established in Theorem 1.1.
Remark 1.4.
We have several comments on Theorem 1.2.
- –
- –
- –
- –
- –
- –
- –
We now strengthen and complement Theorem 1.2 by establishing following result.
Theorem 1.3 (Finite time blow-up for large data in the subcritical case).
Assume that Let be as in Theorem 1.2 (so local wellposedness for (1.6) holds). Let satisfies for 111for example or . Further assume
(1.18) |
Then there exists initial data such that . Moreover if , one has a unique blow-up solution to (1.6) with initial data in the following sense: there exist a unique solution of (1.6) defined on such that
Remark 1.5.
We have several comments for Theorem 1.3.
-
–
For the critical case , similar blowup happens in a Kato norm: If , one would have for certain choice of . See Section 4 for definition of Kato norm.
- –
- –
-
–
Assume ,
and the hypothesis on from Theorem 1.2. Let or or . Then Theorem 1.3 reveals that, there exists data in such that the local solution established in Theorem 1.2 cannot be extend to global in time. In the critical case, it also says that small data assumption in Theorem 1.2 is essentially optimal to establish global existence.
Definition 1.2 (weak solution).
Remark 1.6.
Proceeding as [16, Proposition 3.1] it follows that -integral solutions are weak solution. In that case .
We shall now turn our attention to supercritical case. In this case, we show that there exists positive initial data in that do not generate a (weak) local solution to (1.6). Specifically, we have the following theorem.
Theorem 1.4 (Nonexistence of local positive weak solution in supercritical case).
Remark 1.7.
- –
- –
- –
2. Preliminaries
Notations: The symbol means whereas mean . By we denote . The notation means for some universal constant By we mean . By we mean and .
We shortly denote unweighed Lebesgue space norm by The Schwartz space is denoted by , and the space of tempered distributions is denoted by For and , we introduce the weighted local Lebesgue space given by
where is the set of measurable functions on .
2.1. Lorentz space
The Lorentz space is the space of all complex-valued measurable functions such that where is defined by
with , and denotes the Lebesgue measure on . Therefore
Let us gather some useful results on Lorentz spaces relevant to subsequent our proofs.
Lemma 2.1 (Lemmata 2.2, 2.5 in [26]).
Let , . Then
-
(1)
, the usual Lebesgue -norm.
-
(2)
if .
-
(3)
for .
Lemma 2.2 (Theorems 2.6, 3.4 in [26]).
We have the following inequalities in Lorentz spaces:
-
(1)
(Hölder’s inequality) Let and is such that . Then .
-
(2)
(Young’s inequality) Let and is such that . Then .
2.2. Heat kernel estimate
Let be the symmetric (in variable) heat kernel associated with the operator , i.e.
see [24, Proposition 3.6.]. Then we have the following bounds for :
3. Dissipative estimates in weighted Lebesgue spaces
In order to prove Theorem 1.1, we first show it is enough to prove for (Lemma 3.1), then using a duality (Lemma 3.2) we show it is enough to prove for . Then we crucially use a known heat kernel estimate (Theorem A) to achieve the desired result.
Lemma 3.1.
Let , and Then is bounded from into if and only if is bounded from into with
(3.1) |
for any .
Proof.
Lemma 3.2.
Let and and . Let for and for set . Then
if and only if
Proof.
Note that
This completes the proof. ∎
Proof of Theorem 1.1 (Sufficiency part).
Assume that (1.4) and (1.5) hold. In view of Lemma 3.1 it is enough to prove the case i.e.
(3.2) |
For and applying Theorem A we achieve
(3.3) |
where with as in Theorem A. Set
Then using and (3.3) we have
Splitting the integrations in variable we obtain
Now we show that each of these terms is dominated by which would prove (3.2) to conclude the proof.
Estimate for IV: Using boundedness of and changing the order of integration and Hölder’s inequality we obtain
IV | ||||
where in the last step we have used the hypothesis
Estimate for III: Note that for , we have , and as is radially decreasing we have therefore
III | ||||
where in the last step we have used the hypothesis as in the estimate for IV and the fact that is Schwartz class function.
Estimate for II: We claim that uniformly for all . In fact when we have for all . On the other hand when , using , for we have
This proves the claim. Then
II | ||||
using the claim above and the hypothesis as in the estimate for IV.
Estimate for : Let us treat case by case.
Case By hypothesis . Then using Young’s inequality
Case If , then using Young’s and Holder’s inequalities in Lorentz spaces i.e. Lemma 2.2 we have
If , then by using Lemma 3.3 (1) we choose so that
(3.4) |
and then using Young’s and Holder’s inequalities we achieve
Case Using and Young’s inequality
Note that we have as and (1.4) is assumed. Then by choosing satisfying (3.4) and using Young’s and Holder’s inequalities we achieve
By hypothesis then
Case Since we proceed as in above case and prove the estimate for . Now with the assumption using Lemma 3.3 (2) we choose satisfying
(3.5) |
and obtain
Case If , by Lemma 3.3 (3) we choose satisfying
(3.6) |
so that
If , then we claim . We need to show i.e. . Infact if , then , a contradiction as .
Next we claim . This is because using .
Above claim shows , then we have . Therefore
Case Follows from duality Lemma 3.2 and the above cases. This completes the proof. ∎
Proof of Theorem 1.1 (Necessity part).
Necessity of , : Let be supported in and equal to in with . Then and hence by hypothesis (1.3), we have . On the other hand for
where we have used Theorem A in the second step and in the last step. Since , we must have . Using symmetry of heat kernel see (4) in Remark 1.1. it follows that . This proof is a major modification made to [22, Section 4] where , was treated.
Necessity of : Let If , then using (1.3), we have in (and hence pointwise a.e.) as . Since , using semigroup property, we have in as . Thus which is a contraction.
Necessity of : We prove this by modifying the proof in case in [29, Remark 10]. Let be a smooth non-negative function with support in and take with . Then for and
where we have used Theorem A in the second step, the fact in the third step. Now and . Therefore for we have from (1.3) that
but and as . Therefore we must have . ∎
Lemma 3.3.
4. Local and small data global well-posedness
In this section we prove the well-posedness in critical and subcritical case i.e. when (recall that and ). In order to prove Theorem 1.2, we introduce the Kato space depending on four parameters .
Definition 4.1 (Kato space).
Let and , set . Then the Kato space is defined by
endowed with the norm
Remark 4.1.
By Theorem 1.1, we immediately get the following result (in fact these results are equivalent):
Lemma 4.1.
Let and . Then
if and only if
(4.1) |
Recall that by solution we meant integral solution and therefore, we introduce a nonlinear mapping given by
A fixed point of this map would essentially be a solution to (1.6). Next using Lemma 4.2, we establish the nonlinear estimates in Kato spaces with appropriate conditions on the parameters.
Proposition 4.1 (Nonlinear estimate, sub-critical & critical case).
Remark 4.2.
Note that .
Remark 4.3.
Lemma 4.2.
Assume and let , , and
(4.8) |
then for and we have
Proof.
Proof of Proposition 4.1.
Let us first establish two claims:
Claim I: Let be as in Definition 4.1. Then
(4.9) |
Proof of Claim I: Note that implies . First inequality in (4.6) says . Thus
Claim II:
(4.10) |
Proof of Claim II: For the subcritical case we have Proof of claim:
where in the first and third inequalities we used and in the last step we used . Proof for the ease , we only need to make the first, third nonstrict inequalities by strict inequalities (using ) and last strict inequality by equality (using ). This proves Claim II.
Now note that (4.3), (4.5) implies (4.8) for . By Lemma 4.2 with and (1.7) we have
(4.11) | |||||
where the last inequality is due to the fact . Recall and is convergent if . Taking (4.2), (4.6), (4.9) into account, note that the last time-integral in (4.11) is bounded by
This together with (4.11) implies the first part of the result.
Remark 4.5 (Hypotheses of Proposition 4.1).
In the next result, we prove that there exists parameter such that (4.1) in Lemma 4.1 and (4.3), (4.5), (4.6) in Proposition 4.1 are satisfied.
Lemma 4.3.
Proof.
We need
(4.14) |
and
(4.15) |
Now (4.14) follows if we chosse so that
Choose such that
so that (4.15) is satisfied. Then choose so that
which is possible as as a consequence of (1.15). This completes the proof.
The furthermore more part is clear. ∎
As we are done with linear estimate Lemma 4.1 and nonlinear estimate 4.1 and existence of parameter we are in a position to prove the following well-posedness result which implies Theorem 1.2.
Theorem 4.1 (Local well-posedness in the subcritical weighted Lebesgue space).
Let , satisfy (1.14) , (1.15). Let , satisfy the subcriticality condition defined in (1.9) and
(4.16) |
Let and satisfy hypothesis (4.1) of Lemma 4.1, and hypotheses (4.3), (4.5), (4.6) of Proposition 4.1. Then the Cauchy problem (1.6) is locally well-posed in for arbitrary data . More precisely, the following assertions hold.
-
(1)
(Existence) For any there exist a positive number and an -integral solution to (1.6) satisfying
Moreover, the solution can be extended to the maximal interval .
-
(2)
(Uniqueness in ) Let If satisfy (1.10) with , then on
-
(3)
(Continuous dependence on initial data) For any initial data and in let and be the corresponding existence time given by part (1). Then there exists a constant depending on and such that the corresponding solutions and satisfy
for
-
(4)
(Blow-up criterion in subcritical case )) If then Moreover, the following lower bound of blow-up rate holds: there exists a positive constant independent of such that
(4.17) for .
-
(5)
(Blow-up criterion in critical case ) If is an -integral solution constructed in the assertion (1) and then
-
(6)
(Small data global existence in critical case ) There exists depending only on and such that if satisfies (or in view of Lemma 4.1), then and .
Proof of Theorem 4.1.
Existence in Kato space : Define
with the metric
Then by Lemma 4.1, Proposition 4.1, for we have
(4.18) |
and
(4.19) | |||||
Subcritical case :
Using (4.18), (4.19) and choosing and small enough so that , we find is a contraction in sub-critical case (and hence we have existence of unique solution ). This proves (1), (2).
Critical case : Note that using a density argument we have
Thus we choose so that and where as in (4.18). Then by using (4.18), (4.19) for we have
and
Thus is a contraction in .
This proves (1).
Solution is in :
Using Lemma 4.1, Proposition 4.1
and
Since , solution is indeed in rest of the results follows as in classical case.
Uniqueness, continuous dependency, blow-up, small data global existence are usual as in classical case. ∎
Remark 4.6.
Proof of Theorem 1.2.
Remark 4.7.
In [9], Kato space was not used and hence the they did not achieve uniqueness in mere .
5. Finite time blow-up and nonexistence results
In this section we establish that in the sub-critical and critical case there exists initial data for which solution established by Theorem 1.2 cannot be extended globally in time. Then blow-up alternative (see Theorem 4.1) implies solution must blow-up in finite time. On the other hand for super-critical case, we shall prove that there exists data such that no local weak (hence integral) solution exists.
Before proving the above two we establish the following important lemma which will be used in both he proofs.
Lemma 5.1.
Proof.
Let
where is such that on and supported in . We note that for we have and the estimate
(5.1) | |||||
by choosing
We define a function given by
We note that , since . By using the weak form (1.2), non-negativity of , the above estimate (5.1), Hölder’s inequality and Young’s inequality, the estimates hold:
(5.2) | |||||
where , i.e., and is defined by
The last equality holds only when (1.18) holds. Now from (5.2), we have
which completes the proof. ∎
5.1. Finite time blow-up in critical and subcritical case
The proof of this theorem is based on the arguments of [15, Proposition 2.2, Theorem 2.3] where lifespan of solution for nonlinear Schrödinger equation is studied.
Proof of Theorem 1.3.
Let be a parameter. We take an initial data as , where is given by
(5.3) |
with satisfying
(5.4) |
Then we see and hence by Theorem 1.2, we can define the maximal existence time . Moreover the solution with initial data would be nonnegative as heat kernel is so. Since , it follows from a change of variable and then Lemma 5.1 that for any
which implies
(5.5) |
where .
Claim: There exists such that if , then .
Indeed, on the contrary we assume that for a sequence . Since , we have . The following estimates hold:
which a contradiction and hence the claim is established.
5.2. Nonexistece of weak solution in the supercritical case
In this subsection we give a proof of Theorem 1.4. We only give a sketch of the proof. For the details, we refer to [15, Proposition 2.4, Theorem 2.5] where nonlinear Schrödinger equation is studied.
Proof of Theorem 1.4.
Let . Suppose that the conclusion of Theorem 1.4 does not hold. Then there exists a positive weak solution on to (1.6) (See Definition 1.2) with any initial data in particular for given by (5.3) with satisfying
(5.6) |
Note that such choice is possible as and (1.20) i.e. . Now (5.6) implies . For we have using Lemma 5.1
(5.7) |
Combining Lemma 5.1 and (5.7), we obtain
which leads to a contradiction, as satisfies
This completes the proof. ∎
Acknowledgement: S Haque is thankful to DST–INSPIRE (DST/INSPIRE/04/2022/001457) & USIEF–Fulbright-Nehru fellowship for financial support. S Haque is also thankful to Harish-Chandra Research Institute & University of California, Los Angeles for their excellent research facilities.
References
- [1]
- [2] B. Ben Slimene, S. Tayachi, and F. B. Weissler, Well-posedness, global existence and large time behavior for Hardy–Hénon parabolic equations, Nonlinear Analysis: Theory, Methods & Applications, 152 (2017), pp. 116–148.
- [3] D. G. Bhimani and S. Haque, On inhomogeneous heat equation with inverse square potential, arXiv:2210.09910.
- [4] D. G. Bhimani, M. Majdoub, and R. Manna, Heat equations associated to harmonic oscillator with exponential nonlinearity, arXiv:2306.02828, (2023).
- [5] D. G. Bhimani, R. Manna, F. Nicola, S. Thangavelu, and S. I. Trapasso, Phase space analysis of the Hermite semigroup and applications to nonlinear global well-posedness, Adv. Math., 392 (2021), pp. Paper No. 107995, 18.
- [6] , On heat equations associated with fractional harmonic oscillators, Fract. Calc. Appl. Anal., 26 (2023), pp. 2470–2492.
- [7] H. Brezis and T. Cazenave, A nonlinear heat equation with singular initial data, Journal D’Analyse Mathématique, 68 (1996), pp. 277–304.
- [8] N. Burq, F. Planchon, J. G. Stalker, and A. S. Tahvildar-Zadeh, Strichartz estimates for the wave and Schrödinger equations with the inverse-square potential, J. Funct. Anal., 203 (2003), pp. 519–549.
- [9] N. Chikami, M. Ikeda, and K. Taniguchi, Optimal well-posedness and forward self-similar solution for the Hardy-Hénon parabolic equation in critical weighted Lebesgue spaces, Nonlinear Anal., 222 (2022), pp. Paper No. 112931, 28.
- [10] N. Chikami, M. Ikeda, K. Taniguchi, and S. Tayachi, Unconditional uniqueness and non-uniqueness for Hardy–Hénon parabolic equations, Mathematische Annalen, (2024), pp. 1–61.
- [11] H. Fujita, On the blowing up of solutions of the Cauchy problem for , J. Fac. Sci. Univ. Tokyo Sect. I, 13 (1966), pp. 109–124 (1966).
- [12] N. Ghoussoub and A. Moradifam, Functional inequalities: new perspectives and new applications, vol. 187 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, 2013.
- [13] A. Haraux and F. B. Weissler, Non-uniqueness for a semilinear initial value problem, Indiana University Mathematics Journal, 31 (1982), pp. 167–189.
- [14] M. Hénon, Numerical experiments on the stability of spherical stellar systems, Astronomy and Astrophysics, Vol. 24, p. 229 (1973), 24 (1973), p. 229.
- [15] M. Ikeda and T. Inui, Some non-existence results for the semilinear Schrödinger equation without gauge invariance, J. Math. Anal. Appl., 425 (2015), pp. 758–773.
- [16] M. Ikeda and Y. Wakasugi, Small data blow-up of -solution for the nonlinear Schrödinger equation without gauge invariance, Differential and Integral Equations, 26 (2013), p. 1275.
- [17] N. Ioku, G. Metafune, M. Sobajima, and C. Spina, estimates for homogeneous operators, Commun. Contemp. Math., 18 (2016), pp. 1550037, 14.
- [18] N. Ioku and T. Ogawa, Critical dissipative estimate for a heat semigroup with a quadratic singular potential and critical exponent for nonlinear heat equations, J. Differential Equations, 266 (2019), pp. 2274–2293.
- [19] K. Ishige, On the Fujita exponent for a semilinear heat equation with a potential term, J. Math. Anal. Appl., 344 (2008), pp. 231–237.
- [20] K. Ishige and T. Kawakami, Critical Fujita exponents for semilinear heat equations with quadratically decaying potential, Indiana Univ. Math. J., 69 (2020), pp. 2171–2207.
- [21] H. Kalf, U.-W. Schmincke, J. Walter, and R. Wüst, On the spectral theory of Schrödinger and Dirac operators with strongly singular potentials, in Spectral theory and differential equations (Proc. Sympos., Dundee, 1974; dedicated to Konrad Jörgens), Lecture Notes in Math., Vol. 448, Springer, Berlin-New York, 1975, pp. 182–226.
- [22] R. Killip, C. Miao, M. Visan, J. Zhang, and J. Zheng, Sobolev spaces adapted to the Schrödinger operator with inverse-square potential, Mathematische Zeitschrift, 288 (2018), pp. 1273–1298.
- [23] G. Metafune, L. Negro, and C. Spina, Sharp kernel estimates for elliptic operators with second-order discontinuous coefficients, Journal of Evolution Equations, 18 (2018), pp. 467–514.
- [24] , Maximal regularity for elliptic operators with second-order discontinuous coefficients, J. Evol. Equ., 21 (2021), pp. 3613–3637.
- [25] C. Miao, B. Yuan, and B. Zhang, Well-posedness of the Cauchy problem for the fractional power dissipative equations, Nonlinear Anal., 68 (2008), pp. 461–484.
- [26] R. O’Neil, Convolution operators and spaces, Duke Mathematical Journal, 30 (1963), pp. 129–142.
- [27] R. Pinsky, The Fujita exponent for semilinear heat equations with quadratically decaying potential or in an exterior domain, J. Differential Equations, 246 (2009), pp. 2561–2576.
- [28] R. G. Pinsky, Existence and nonexistence of global solutions for in , J. Differential Equations, 133 (1997), pp. 152–177.
- [29] S. Tayachi and F. B. Weissler, New life-span results for the nonlinear heat equation, Journal of Differential Equations, 373 (2023), pp. 564–625.
- [30] J. L. Vazquez and E. Zuazua, The Hardy inequality and the asymptotic behaviour of the heat equation with an inverse-square potential, J. Funct. Anal., 173 (2000), pp. 103–153.
- [31] F. B. Weissler, Existence and nonexistence of global solutions for a semilinear heat equation, Israel J. Math., 38 (1981), pp. 29–40.
- [32] Q. S. Zhang, The quantizing effect of potentials on the critical number of reaction-diffusion equations, J. Differential Equations, 170 (2001), pp. 188–214.