Home Boundary Hölder continuity of stable solutions to semilinear elliptic problems in 𝐶1,1 domains
Article
Licensed
Unlicensed Requires Authentication

Boundary Hölder continuity of stable solutions to semilinear elliptic problems in 𝐶1,1 domains

  • Iñigo U. Erneta ORCID logo EMAIL logo
Published/Copyright: November 9, 2024

Abstract

This article establishes the boundary Hölder continuity of stable solutions to semilinear elliptic problems in the optimal range of dimensions n 9 , for C 1 , 1 domains. We consider equations L u = f ( u ) in a bounded C 1 , 1 domain Ω R n , with u = 0 on Ω , where 𝐿 is a linear elliptic operator with variable coefficients and f C 1 is nonnegative, nondecreasing, and convex. The stability of 𝑢 amounts to the nonnegativity of the principal eigenvalue of the linearized equation L f ( u ) . Our result is new even for the Laplacian, for which [X. Cabré, A. Figalli, X. Ros-Oton and J. Serra, Stable solutions to semilinear elliptic equations are smooth up to dimension 9, Acta Math. 224 (2020), 2, 187–252] proved the Hölder continuity in C 3 domains.

Award Identifier / Grant number: MDM-2014-0445-18-1

Award Identifier / Grant number: MTM2017-84214-C2-1-P

Award Identifier / Grant number: PID2021-123903NB-I00

Funding statement: The author acknowledges financial support from MINECO grant MDM-2014-0445-18-1 through the María de Maeztu Program for Units of Excellence in R&D. He is additionally supported by Spanish grants MTM2017-84214-C2-1-P and PID2021-123903NB-I00 funded by MCIN/AEI/10.13039/501100011033 and by ERDF “A way of making Europe”. The author is also supported by Catalan project 2021 SGR 00087.

A Approximating C 1 , 1 domains by smooth ones from the interior

In this appendix, we show that bounded domains of class C 1 , 1 can be approximated from the interior by smooth sets satisfying uniform bounds. This is a well-known result in the literature, and is valid more generally for domains of class C k , α with k 1 and α [ 0 , 1 ] . We include a proof for the sake of completeness. Our proof follows the approach suggested by Gilbarg and Trudinger in [18, Problem 6.9].

First, we recall the definition of C 1 , 1 domains.

Definition A.1

A bounded domain Ω R n is of class C 1 , 1 if, at each point x 0 Ω , there are a ball B = B ρ ( x 0 ) and a one-to-one mapping Ψ of 𝐵 onto U R n such that

  1. Ψ ( B Ω ) R + n ;

  2. Ψ ( B Ω ) R + n ;

  3. Ψ C 1 , 1 ( B ̄ ) and Ψ 1 C 1 , 1 ( U ̄ ) .

Equivalently, Ω is of class C 1 , 1 if each point of Ω has a neighborhood in which Ω is the graph of a C 1 , 1 function of n 1 of the coordinates.

Every such domain can be written as the positive set of a C 1 , 1 function.

Lemma A.2

Let Ω R n be a bounded domain of class C 1 , 1 . Then there is a function Φ C 1 , 1 ( R n ) such that Ω = { Φ > 0 } , Φ = 0 on Ω , and Φ ( x ) 0 for all x Ω .

Proof

By compactness, Ω may be covered by finitely many balls { B j = B ρ j ( x j ) } j = 1 N , with x j Ω , ρ j > 0 , such that there are flattening maps Ψ j C 1 , 1 ( B j ̄ ) as in Definition A.1.

Let ρ > 0 be sufficiently small so that the set B 0 : = { x Ω : dist ( x , Ω ) > ρ } satisfies Ω ̄ j = 0 N B j , and consider a partition of unity { η j } j = 0 N subordinated to the covering { B j } j = 0 N . The function

Φ = η 0 + j = 1 N η j Ψ j n

now satisfies the desired properties. ∎

Remark A.3

Notice that, by construction, the function Φ above is compactly supported and takes negative values in a bounded neighborhood of Ω outside Ω ̄ .

Regularizing Φ and taking appropriate superlevel sets, we obtain the approximation.

Lemma A.4

Suppose that Ω = { Φ > 0 } R n is a bounded domain of class C 1 , 1 , with Φ C 1 , 1 ( R n ) as in Lemma A.2 above. Then there is an exhaustion of Ω by smooth sets Ω k = { Φ k > 0 } ,[5] where the functions Φ k C ( R n ) satisfy

Φ k C 1 ( R n ) + | Φ k | 1 L ( Ω k ) C

for some constant 𝐶 depending only on Φ and Ω. Moreover, we have that Ω k Ω in the sense of the Hausdorff distance.[6]

Proof

Since Ω = { Φ > 0 } and Φ 0 on Ω , by continuity of the gradient, Φ is comparable to the distance function d ( x ) : = dist ( x , Ω ) in Ω ̄ , i.e.,

(A.1) L 1 d ( x ) Φ ( x ) L d ( x ) for x Ω ̄ ,

for some L 1 . For 𝑥 outside Ω, we define d ( x ) : = dist ( x , Ω ) < 0 .

Consider a mollifying sequence ( η ε ) ε > 0 with supp η ε B ε . We will take the functions

(A.2) Φ k : = Φ η ε k 2 L ε k

for an appropriate sequence ε k 0 . Since Φ k C 1 = Φ k L + D 2 Φ k L , recalling that D 2 Φ k L ( R n ) = [ Φ k ] C 0 , 1 ( R n ) (since Φ k is C ) and Φ C 1 , 1 , the uniform bounds for Φ k C 1 stated in the lemma hold by the standard properties of convolutions.

Given δ > 0 , from (A.1), we see that

Φ ( x ) > δ / L for x { x R n : d ( x ) > δ } = { d > δ } ,

and taking the convolution with η ε , it follows that

(A.3) Φ η ε ( x ) > δ / L for x { d > δ + ε } .

Similarly, for δ ̃ > 0 , we have that

Φ ( x ) L δ ̃ for x { d δ ̃ } ,

and regularizing, we obtain

(A.4) Φ η ε ( x ) L δ ̃ for x { d δ ̃ ε } .

Letting δ ̃ = δ / L 2 , since L δ ̃ = δ / L , by (A.3) and (A.4), we have that

{ d > δ + ε } { Φ η ε > δ / L } { d > δ / L 2 ε } ,

and the choice δ = 2 L 2 ε now yields the inclusions

(A.5) { d > ε ( 2 L 2 + 1 ) } { Φ η ε > 2 L ε } { d > ε } .

Next, we construct the sequence ε k in (A.2) inductively. Fix ε 1 > 0 small. For k 1 , we define ε k + 1 : = ε k / { 2 ( 2 L 2 + 1 ) } = ε 1 / { 2 ( 2 L 2 + 1 ) } k . Hence, by (A.5), the sets

Ω k : = { Φ k > 0 } = { Φ η ε k > 2 L ε k }

satisfy

(A.6) Ω k ̄ { d ε k } { d > ε k + 1 ( 2 L 2 + 1 ) } Ω k + 1 Ω .

They clearly exhaust Ω since ε k 0 and thus Ω = k { d ε k 1 } k Ω k . Furthermore, the inclusions (A.6) show that Ω k is at a Hausdorff distance of at most ε k ( 2 L 2 + 1 ) from Ω , and hence Ω k Ω with respect to this distance.

It remains to prove the lower bound for | Φ k | on Ω k , which will also show the boundary Ω k to be smooth. Let x R n with | d ( x ) | = dist ( x , Ω ) = | x x 0 | for some x 0 Ω . Since

Φ ( x ) Φ ( x 0 ) | Φ ( x 0 ) | | Φ ( x 0 ) | [ Φ ] C 0 , 1 ( R n ) | d ( x ) | ,

we have the lower bound

Φ ( x ) Φ ( x 0 ) | Φ ( x 0 ) | 1 2 | Φ | 1 L ( Ω ) 1 for x { ρ < d < ρ } ,

where ρ > 0 depends only on | Φ | 1 L ( Ω ) and [ Φ ] C 0 , 1 ( R n ) . Taking the convolution with η ε k in this last inequality, we obtain

(A.7) Φ k ( x ) Φ ( x 0 ) | Φ ( x 0 ) | 1 2 | Φ | 1 L ( Ω ) 1 for x { ( ρ ε k ) < d < ρ ε k } .

By (A.6), the boundary Ω k is at a distance of at most ε k ( 2 L 2 + 1 ) from Ω . Hence, choosing ε 1 > 0 sufficiently small such that ρ ε 1 > ε 1 ( 2 L 2 + 1 ) , from (A.7) and the definition of ε k = ε 1 / { 2 ( 2 L 2 + 1 ) } k 1 , we deduce

| Φ k ( x ) | Φ k ( x ) Φ ( x 0 ) | Φ ( x 0 ) | 1 2 | Φ | 1 L ( Ω ) 1 on Ω k ,

and therefore | Φ k | 1 L ( Ω k ) 2 | Φ | 1 L ( Ω ) , which concludes the proof. ∎

We conclude this appendix with a description of the flattening procedure and a technical lemma used in the proof of Theorem 1.1 in Section 5.

Let x 0 Ω . Rotating the coordinate axes, we may assume that Φ ( x 0 ) = n Φ ( x 0 ) e n , with | Φ ( x 0 ) | = n Φ ( x 0 ) > 0 . For x = ( x , x n ) R n 1 × R , the map Ψ : R n R n given by

Ψ ( x ) : = ( ( x x 0 ) , Φ ( x ) n Φ ( x 0 ) )

is a local diffeomorphism around x 0 which flattens out the boundary Ω . More precisely, by a quantitative version of the Inverse Function Theorem, we have the following.

Lemma A.5

There are (small) numbers 0 < R 1 < R 2 and ρ > 0 depending only on Φ C 0 , 1 ( R n ) and | Φ | 1 L ( Ω ) such that

(A.8) Ψ ( B R 1 ( x 0 ) Ω ) B ρ / 2 + B ρ + Ψ ( B R 2 ( x 0 ) Ω ) .

Proof

By translation, we may assume that x 0 = 0 Ω . Since the map Ψ satisfies Ψ ( 0 ) = 0 and D Ψ ( 0 ) = I , choosing R 2 > 0 small such that

(A.9) | D Ψ ( x ) D Ψ ( z ) | 1 / 2 for all x , z B R 2 ,

by [20, Chapter XIV, Lemma 1.3], we deduce that, for all y B R 2 / 2 , there is a unique x B R 2 such that Ψ ( x ) = y . Thus, we obtain the second inclusion

(A.10) B R 2 / 2 + Ψ ( B R 2 Ω ) .

Using that [ D Ψ ] C 0 , 1 ( R n ) [ Φ ] C 0 , 1 ( R n ) | Φ | 1 L ( Ω ) , it is easy to check that condition (A.9) is fulfilled if

(A.11) R 2 ( 4 [ Φ ] C 0 , 1 ( R n ) | Φ | 1 L ( Ω ) ) 1 .

To show the first inclusion in (A.8), we proceed as above but considering the inverse map Ψ 1 instead of Ψ. If R 1 > 0 is such that

(A.12) | D Ψ 1 ( x ̃ ) D Ψ 1 ( z ̃ ) | 1 / 2 for all x ̃ , z ̃ B 2 R 1 ,

then, again by [20, Chapter XIV, Lemma 1.3], we deduce

(A.13) Ψ ( B R 1 Ω ) B 2 R 1 + .

For (A.12) to hold, it suffices to take R 1 > 0 sufficiently small such that

(A.14) [ D Ψ 1 ] C 0 , 1 ( B 2 R 1 ) 4 R 1 1 / 2 .

Hence, to quantify the smallness of R 1 , we have to estimate [ D Ψ 1 ] C 0 , 1 ( B 2 R 1 ) .

Let x ̃ = Ψ ( x ) and y ̃ = Ψ ( y ) in B 2 R 1 . If R 1 R 2 , with R 2 satisfying (A.11), then

1 n Φ ( x ) 2 | Φ | 1 L ( Ω ) ,

and hence

| D Ψ 1 ( x ̃ ) D Ψ 1 ( y ̃ ) | 2 2 | Φ ( x ) Φ ( y ) | 2 | n Φ ( x ) | 2 + ( 2 | Φ ( y ) | 2 + | n Φ ( 0 ) | 2 ) | n Φ ( x ) n Φ ( y ) | 2 | n Φ ( x ) | 2 | n Φ ( y ) | 2 8 ( 1 + 6 Φ L ( R n ) 2 | Φ | 1 L ( Ω ) 2 ) | Φ | 1 L ( Ω ) 2 [ Φ ] C 0 , 1 ( R n ) 2 | x y | 2 .

Moreover, since

| x y | [ D Ψ 1 ] L ( B 2 R 1 ) | x ̃ y ̃ | ( 1 + 8 Φ L 2 | Φ | 1 L ( Ω ) 2 ) 1 / 2 | x ̃ y ̃ | ,

by the above, we conclude

[ D Ψ 1 ] C 0 , 1 ( B 2 R 1 ) 2 2 ( 1 + 8 Φ L 2 | Φ | 1 L ( Ω ) 2 ) | Φ | 1 L ( Ω ) [ Φ ] C 0 , 1 ( R n )

for all R 1 R 2 satisfying (A.11). Therefore, if

(A.15) R 1 ( 16 2 ( 1 + 8 Φ L ( R n ) 2 | Φ | 1 L ( Ω ) 2 ) | Φ | 1 L ( Ω ) [ Φ ] C 0 , 1 ( R n ) ) 1 ,

then (A.14) holds and (A.13) follows.

Finally, choosing R 1 , R 2 > 0 satisfying (A.11), (A.15), and 8 R 1 R 2 , we may take ρ = 4 R 1 , and the inclusions (A.10) and (A.13) then yield the claim (A.8). ∎

B On the uniqueness of stable solutions

Here, we prove the uniqueness of stable solutions to non-variational equations involving convex nonlinearities. For this, we employ some fundamental results of Berestycki, Nirenberg, and Varadhan [2] on the principal eigenfunction. Compare the following statement with the one in Dupaigne’s book [14, Proposition 1.3.1].

Proposition B.1

Given Ω R n a bounded domain, let u 1 , u 2 C 0 ( Ω ̄ ) W loc 2 , n ( Ω ) be two stable solutions of the equation L u = f ( u ) in Ω, with u = 0 on Ω . Assume that f C 1 ( R ) is convex. Then either u 1 = u 2 or f ( u ) = μ 1 [ L , Ω ] u on the ranges of u 1 and u 2 .

Remark B.2

Here μ 1 [ L , Ω ] denotes the principal (or smallest) eigenvalue of 𝐿 in Ω, with the sign convention L φ = μ 1 φ . It is characterized by (see [2])

μ 1 [ L , Ω ] = sup { μ : there is a function φ > 0 W loc 2 , n ( Ω ) satisfying L φ + μ φ 0 in Ω } .

Moreover, 𝐿 can be any uniformly elliptic second order operator. In particular, we allow it to have zero order terms.

Remark B.3

In the proof of Theorem 1.1 above, we only need a weaker version of Proposition B.1. Namely, we could assume additionally that u 1 u 2 , which admits a shorter proof. However, the present statement might be more useful in applications.

Proof

Assume that u 1 u 2 and consider the difference w : = u 2 u 1 . Let

Ω + : = { x Ω : u 2 ( x ) > u 1 ( x ) } = { w > 0 } ,

and assume that Ω + 0 (otherwise, we exchange the roles of u 1 and u 2 ). By convexity, we have L w = f ( u 2 ) f ( u 1 ) f ( u 2 ) w , and hence

(B.1) J u 2 w = L w + f ( u 2 ) w 0 in Ω .

By the monotonicity of the principal eigenvalue with respect to the domain, since u 2 is stable, we have that

(B.2) μ 1 [ J u 2 , Ω + ] μ 1 [ J u 2 , Ω ] 0 .

Since w > 0 in Ω + , by (B.1) and (B.2), it follows that

{ J u 2 w + μ 1 [ J u 2 , Ω + ] w 0 in Ω + , w = 0 on Ω + .

Applying [2, Corollary 2.2], we deduce that 𝑤 is a positive principal eigenfunction of J u 2 in Ω + , that is, satisfying

(B.3) J u 2 w + μ 1 [ J u 2 , Ω + ] w = 0 in Ω + .

Using the equation L w = f ( u 2 ) f ( u 1 ) in Ω + , from (B.3), we see that

(B.4) f ( u 2 ) f ( u 1 ) f ( u 2 ) ( u 2 u 1 ) = μ 1 [ J u 2 , Ω + ] ( u 2 u 1 ) 0 in Ω + .

By convexity, we also have the opposite inequality f ( u 2 ) f ( u 1 ) f ( u 2 ) ( u 2 u 1 ) 0 . Hence, by (B.4), we conclude μ 1 [ J u 2 , Ω + ] = 0 , and the nonlinearity 𝑓 is affine in the union of intervals [ u 1 ( x ) , u 2 ( x ) ] with x Ω + .

If 𝑓 is of the form f ( u ) = a u + b in the ranges above, then L w = a w in Ω + , with w > 0 , and therefore a = μ 1 [ L , Ω + ] . Now, to have nontrivial solutions of

L u = μ 1 [ L , Ω + ] u + b ,

the Fredholm alternative forces b = 0 . We conclude that f ( u ) = μ 1 [ L , Ω + ] u in the ranges of u 1 and u 2 in Ω + .

If Ω : = { w < 0 } , then arguing as above with w in place of 𝑤, we deduce that

J u 1 w + μ 1 [ J u 1 , Ω ] w = 0 in Ω ,

with μ 1 [ J u 1 , Ω ] = 0 and f ( u ) = μ 1 [ L , Ω ] u in the ranges of u 2 and u 1 in Ω .

The regularity of 𝑓 and the continuity of the solutions now forces

μ 1 [ L , Ω + ] = μ 1 [ L , Ω ] .

Now, on the one hand, by the stability of u 2 , we have that

0 μ 1 [ J u 2 , Ω ] = μ 1 [ L + μ 1 [ L , Ω + ] , Ω ] = μ 1 [ L , Ω ] μ 1 [ L , Ω + ] .

On the other hand, since Ω + Ω by assumption, the strict monotonicity of μ 1 yields[7]

μ 1 [ L , Ω ] μ 1 [ L , Ω + ] < 0 .

This contradiction forces either Ω + or Ω to be empty, and therefore f ( u ) = μ 1 [ L , Ω ] u in the ranges of u 1 and u 2 , which was the claim. ∎

Acknowledgements

The author wishes to thank Xavier Cabré for useful discussions on the topic of this article, as well as for his encouragement over the years. He also thanks the anonymous referee for their suggestions, which improved the presentation of the paper.

References

[1] H. Berestycki, A. Kiselev, A. Novikov and L. Ryzhik, The explosion problem in a flow, J. Anal. Math. 110 (2010), 31–65. 10.1007/s11854-010-0002-7Search in Google Scholar

[2] H. Berestycki, L. Nirenberg and S. R. S. Varadhan, The principal eigenvalue and maximum principle for second-order elliptic operators in general domains, Comm. Pure Appl. Math. 47 (1994), no. 1, 47–92. 10.1002/cpa.3160470105Search in Google Scholar

[3] H. Brezis, Is there failure of the inverse function theorem?, Morse theory, minimax theory and their applications to nonlinear differential equations, New Stud. Adv. Math. 1, International Press, Somerville (2003), 23–33. Search in Google Scholar

[4] H. Brezis and J. L. Vázquez, Blow-up solutions of some nonlinear elliptic problems, Rev. Mat. Univ. Complut. Madrid 10 (1997), no. 2, 443–469. 10.5209/rev_REMA.1997.v10.n2.17459Search in Google Scholar

[5] X. Cabré, Regularity of minimizers of semilinear elliptic problems up to dimension 4, Comm. Pure Appl. Math. 63 (2010), no. 10, 1362–1380. 10.1002/cpa.20327Search in Google Scholar

[6] X. Cabré, Estimates controlling a function by only its radial derivative and applications to stable solutions of elliptic equations, preprint (2022), https://arxiv.org/abs/2211.13033. Search in Google Scholar

[7] X. Cabré and A. Capella, Regularity of radial minimizers and extremal solutions of semilinear elliptic equations, J. Funct. Anal. 238 (2006), no. 2, 709–733. 10.1016/j.jfa.2005.12.018Search in Google Scholar

[8] X. Cabré, A. Figalli, X. Ros-Oton and J. Serra, Stable solutions to semilinear elliptic equations are smooth up to dimension 9, Acta Math. 224 (2020), no. 2, 187–252. 10.4310/ACTA.2020.v224.n2.a1Search in Google Scholar

[9] X. Cabré, P. Miraglio and M. Sanchón, Optimal regularity of stable solutions to nonlinear equations involving the 𝑝-Laplacian, Adv. Calc. Var. 15 (2022), no. 4, 749–785. 10.1515/acv-2020-0055Search in Google Scholar

[10] X. Cabré and X. Ros-Oton, Regularity of stable solutions up to dimension 7 in domains of double revolution, Comm. Partial Differential Equations 38 (2013), no. 1, 135–154. 10.1080/03605302.2012.697505Search in Google Scholar

[11] F. Costa, G. F. de Souza and M. Montenegro, Extremal solutions of strongly coupled nonlinear elliptic systems and L -boundedness, J. Math. Anal. Appl. 513 (2022), no. 2, Article ID 126225. 10.1016/j.jmaa.2022.126225Search in Google Scholar

[12] C. Cowan and N. Ghoussoub, Regularity of the extremal solution in a MEMS model with advection, Methods Appl. Anal. 15 (2008), no. 3, 355–360. 10.4310/MAA.2008.v15.n3.a7Search in Google Scholar

[13] M. G. Crandall and P. H. Rabinowitz, Some continuation and variational methods for positive solutions of nonlinear elliptic eigenvalue problems, Arch. Ration. Mech. Anal. 58 (1975), no. 3, 207–218. 10.1007/BF00280741Search in Google Scholar

[14] L. Dupaigne, Stable solutions of elliptic partial differential equations, Chapman & Hall/CRC Monogr. Surveys Pure Appl. Math. 143, Chapman & Hall/CRC, Boca Raton 2011. 10.1201/b10802Search in Google Scholar

[15] I. U. Erneta, Stable solutions to semilinear elliptic equations for operators with variable coefficients, Commun. Pure Appl. Anal. 22 (2023), no. 2, 530–571. 10.3934/cpaa.2022162Search in Google Scholar

[16] I. U. Erneta, Energy estimate up to the boundary for stable solutions to semilinear elliptic problems, J. Differential Equations 378 (2024), 204–233. 10.1016/j.jde.2023.09.023Search in Google Scholar

[17] I. M. Gel’fand, Some problems in the theory of quasilinear equations, Amer. Math. Soc. Transl. (2) 29 (1963), 295–381. 10.1090/trans2/029/12Search in Google Scholar

[18] D. Gilbarg and N. S. Trudinger, Elliptic partial differential equations of second order, 2nd ed., Springer, Berlin 2001. 10.1007/978-3-642-61798-0Search in Google Scholar

[19] D. D. Joseph and T. S. Lundgren, Quasilinear Dirichlet problems driven by positive sources, Arch. Ration. Mech. Anal. 49 (1972/73), 241–269. 10.1007/BF00250508Search in Google Scholar

[20] S. Lang, Real and functional analysis, 3rd ed., Grad. Texts in Math. 142, Springer, New York 1993. 10.1007/978-1-4612-0897-6_1Search in Google Scholar

[21] Y. Martel, Uniqueness of weak extremal solutions of nonlinear elliptic problems, Houston J. Math. 23 (1997), no. 1, 161–168. Search in Google Scholar

[22] G. Nedev, Regularity of the extremal solution of semilinear elliptic equations, C. R. Acad. Sci. Paris Sér. I Math. 330 (2000), no. 11, 997–1002. 10.1016/S0764-4442(00)00289-5Search in Google Scholar

[23] P. Sternberg and K. Zumbrun, Connectivity of phase boundaries in strictly convex domains, Arch. Ration. Mech. Anal. 141 (1998), no. 4, 375–400. 10.1007/s002050050081Search in Google Scholar

[24] S. Villegas, Boundedness of extremal solutions in dimension 4, Adv. Math. 235 (2013), 126–133. 10.1016/j.aim.2012.11.015Search in Google Scholar

Received: 2023-07-12
Revised: 2024-09-20
Published Online: 2024-11-09
Published in Print: 2025-02-01

© 2024 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 11.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/crelle-2024-0084/html
Scroll to top button