Recent rigidity results for graphs with prescribed mean curvature

This survey describes some recent rigidity results obtained by the authors for the prescribed mean curvature problem on graphs $u : M \rightarrow \mathbb{R}$. Emphasis is put on minimal, CMC and capillary graphs, as well as on graphical solitons for the mean curvature flow, in warped product ambient spaces. A detailed analysis of the mean curvature operator is given, focusing on maximum principles at infinity, Liouville properties, gradient estimates. Among the geometric applications, we mention the Bernstein theorem for positive entire minimal graphs on manifolds with non-negative Ricci curvature, and a splitting theorem for capillary graphs over an unbounded domain $\Omega \subset M$, namely, for CMC graphs satisfying an overdetermined boundary condition.


Introduction
The purpose of this survey is to illustrate some recent rigidity results, obtained by the authors, about the prescribed mean curvature problem on a complete Riemannian manifold . It is a pleasure to dedicate our work to Alberto Farina on the occasion of his birthday, and to have the opportunity to discuss some of his significant contributions to the theory. Hereafter, we shall always restrict to real valued functions ∶ → ℝ, and their associated entire graphical hypersurfaces Σ ↪ ℝ × described by the map (throughout the paper, the word entire will mean defined on all of ). Taking into account the huge literature on hypersurfaces with prescribed mean curvature, notably on minimal and constant mean curvature (CMC) ones, some choices have to be made, and we decided to exclusively focus on references considering the graphical case in codimension 1. Thus, unfortunately, important works concerning embedded (immersed) hypersurfaces have been omitted, unless they provide new results in the graphical case too. Concerning the prescribed mean curvature problem in Euclidean space, that will be briefly recalled in the next section, we recommend the survey [55] for an excellent account, that also includes Liouville theorems for more general operators.

The prescribed curvature problem in ℝ ×
To put our results into perspective, we briefly recall the main achievements concerning the problem on ℝ . We begin with entire minimal graphs in ℝ +1 , described by functions ∶ ℝ → ℝ solving the minimal surface equation that is the Euler-Lagrange equation of the area functional with respect to compactly supported variations. In 1915, Bernstein [8] proved that the only solutions of (1) on ℝ 2 are affine functions (his proof, highly nontrivial, has been perfected in [79,96], where a missing detail was corrected). Since then, various other proofs appeared in the literature, each one exploiting some different peculiarities of the two dimensional setting, especially: the possibility to use holomorphic function theory and the uniformization theorem [101,77,65,48], or the quadratic volume growth of ℝ 2 coupled with the specific form of the Gauss-Bonnet theorem [108,123]. A presentation of some of them can be found in Farina's survey [55]. The problem whether Bernstein's rigidity theorem holds for minimal graphs over ℝ if ≥ 3 challenged researchers for long, and the efforts to answer this question gave rise to many new and deep techniques in Geometric Analysis. The first proof of Bernstein theorem that allowed for a possible higher dimensional generalization was given by Fleming [67], and since then, in few years, the problem has been completely solved affirmatively in dimension ≤ 7 by works of De Giorgi ( [42,43], = 3), Almgren ([1], = 4) and Simons ([125], ≤ 7), while Bombieri, De Giorgi and Giusti [12] constructed the first counterexample if ≥ 8 (a large class of further examples was later given by Simon [124]). Summarizing, the following theorem holds: If satisfies some a-priori bounds, then more rigidity is expected. Indeed, as a consequence of the local gradient estimate for minimal graphs ∶ ( ) ⊂ ℝ → ℝ due to Bombieri, De Giorgi and Miranda [11]: for some constants ( ), one deduces both that positive minimal graphs on ℝ are constant, for each ≥ 2 (ℬ2) and that minimal graphs on ℝ with at most linear growth have bounded gradient, thus, by work of Moser [98] and De Giorgi [42,43], they are affine.
Further splitting results in the spirit of (ℬ3) can be obtained by just assuming the boundedness of some of the directional derivatives of . In particular, improving on work of Bombieri and Giusti [13], Farina in [53,54] showed that is affine provided that − 7 partial derivatives of are bounded on one side. The possible existence of a similar results for entire minimal graphs defined on more general complete Riemannian manifolds ( , ) heavily depends on the geometry of . If ℝ× is endowed with the product metric ⟨ , ⟩ = d 2 + then the (normalized) mean curvature of the image ( ) in the upward pointing direction is given by formally the same as for ℝ , where now , div are the gradient and divergence in ( , ). In particular, entire minimal graphs solve (1) on . However, for instance if is the hyperbolic space ℍ , (1) admits uncountably many bounded solutions. Indeed, the Plateau problem at infinity is solvable for every continuous on the boundary at infinity ∞ ℍ , that is, there exists solving (1) on ℍ and approaching at infinity. The result was shown in [100] for = 2, while the higher dimensional case, to our knowledge, first appeared in 2010 as a special case of [52]. We note in passing that, as a consequence of Theorem 3.2 below, any solution of Plateau's problem on ℍ has a bounded gradient.
Recently, various interesting works investigated the Plateau's problem at infinity on Cartan-Hadamard (CH) manifolds , that is, on complete, simply connected manifolds with nonpositive sectional curvature Sec, elucidating the sharp thresholds of Sec to guarantee its solvability or not (cf. [111,18,19,20,78,22]). The picture is quite subtle and an updated account can be found in [75]. In particular, having denoted with the distance from a fixed origin ∈ , we mention that the Plateau problem at infinity on such is solvable whenever any of the following pinching conditions hold outside of a compact set ( ): ≤ Sec ≤ − 2 (by [111], see also [18,Cor 1.7]), for constants , > 0, ′ > 1. Hereafter, an inequality Sec ≤ ( ) will shortly mean Sec( ) ≤ ( ) for every ∈ and 2-plane ≤ .
In general, a two-sided bound on Sec is necessary, although a (suitable) upper bound suffices on rotationally symmetric CH manifolds (cf. [22]). Indeed, in [78] the authors constructed a 3-dimensional CH manifold with Sec ≤ − 2 < 0, that admits bounded solutions of the minimal surface equation, but for which the Dirichlet problem at infinity is not solvable for any nonconstant boundary datum . The dimension restriction is not removable, since for = 2 the one-sided condition Sec ≤ − 2 < 0 suffices for the solvability of the Plateau problem at infinity, see [71]. We also mention the surprising existence of parabolic, entire minimal graphs in 2 × ℝ, with 2 a CH surface satisfying Sec ≤ − 2 < 0, obtained in [34,71]. It is therefore tempting to ask whether results similar to those for ℝ may hold on manifolds with non-negative curvature, and in view of the structure theory initiated by Cheeger and Colding in [24,25,26,27], a particularly intriguing class to investigate in this respect is that of manifolds with Ric ≥ 0. This is enforced by the analogy with the behaviour of harmonic functions on . In fact, positive harmonic functions on a manifold with Ric ≥ 0 are constant [131,30], and harmonic functions with linear growth provide splitting directions for the tangent cones at infinity of , and for itself provided that Ric ≥ 0 is strengthened to Sec ≥ 0 [28]. The parallelism is even more evident by considering that (1) rewrites as the harmonicity of on Σ with its induced metric.
While we are aware of no references considering the full Bernstein theorem (ℬ1) on spaces different from ℝ , and in fact the problem seems quite challenging even on manifolds with nonnegative sectional curvature, recently some progresses were made on (ℬ2), (ℬ3). Below, we comment on the following recent result, independently due to the authors and to Ding [44]: Under the further condition Sec ≥ − 2 for some constant > 0, Theorem 1.1 was proved by Rosenberg, Schulze and Spruck in [117]. If the Ricci and sectional curvatures are allowed to be negative somewhere, then the situation is different: in this respect, we quote the recent [22] by Casteras, Heinonen and Holopainen, ensuring the constancy of positive minimal graphs with at most linear growth on manifolds with asymptotically non-negative sectional curvature and only one end. Regarding property (ℬ3), the rigidity of manifolds supporting linearly growing, nonconstant minimal graphs was studied in [45,46] on manifolds satisfying the curvature tensor decays quadratically.
Together with equation (1), that we will investigate in Section 3, in the present survey we will also consider energies containing a potential term, of the type with 0 < ∈ ( ) and ∈ 1 (ℝ). The corresponding Euler-Lagrange equation becomes with = ′ . This family of equations includes the capillarity equation ( ) ( ) = , for some constant > 0, that has been thoroughly studied, for instance, in [64].

Notation 1.2.
Hereafter, given two positive functions , , with ≳ we mean ≥ for some positive constant , and similarly for ≲ . We write ∈ 1 (+∞) to denote and similarly for ∈ 1 (0 + ). When there is no risk of confusion, we write ∞ instead of +∞, −∞.

Remark 1.4.
Condition ≢ 0, together with ( ) ≡ 0, guarantees that is one sided bounded and thus the applicability of (ℬ2). It is therefore essential, as the example of affine functions show.
Remark 1.5. The result establishes that the only graphs with constant mean curvature (CMC) defined on the entire ℝ must be minimal, a result first proved by Heinz [76] for = 2, and later extended independently by Chern [31] and Flanders [66] to all ≥ 3. As a matter of fact, Salavessa in [118] observed that the result is true on each manifold with zero Cheeger constant, in particular, on complete manifolds with subexponential volume growth.
In particular, in the assumptions of Theorem 1.3 the only solution of the capillarity equation on ℝ is ≡ 0. In Theorem 1.1 below, we shall prove that Theorem 1.3 holds on each manifold with Ric ≥ 0; precisely, the implication ( ) ≡ 0 only requires that has polynomial volume growth of geodesic balls, a condition implied (via Bishop-Gromov comparison), but not equivalent, to Ric ≥ 0. On the contrary, as we shall see below, Theorem 1.3 fails on manifolds with faster than polynomial volume growth, cf. Remark 2.22. We note that the restriction < 1 can be relaxed for the capillarity case, the sharp threshold being < 2 by recent work of Farina and Serrin [60], and we will come back on this later, see Theorem 2.28 below.
Remark 1.6. Regarding the case of non-constant weights , let > 0 on ℝ + , and assume that is increasing for large enough . If then (7) is necessary and sufficient to guarantee that solutions of This result can be obtained, for instance, by combining Theorems 10.3 and Corollary 10.23 in [10]. The quadratic lower bound for is sharp, since positive entire solutions were constructed in [49,9] when (7) holds and Remark 1.7. We mention that (7) is sufficient for the validity of the Liouville property (9) on complete manifolds satisfying for some constant > 0, where is the distance from a fixed origin , and cut( ) is the cut-locus of . This result follows, for example, as a particular case of [10,Cor 10.23], and condition (10) is sharp for its validity. Furthermore, (7) is also a necessary condition for (9) for a large class of manifolds, including each Cartan-Hadamard space, see [10,Thm. 2.31]. By Bishop-Gromov comparison theorem, note that (10) implies the following quadratic exponential bound for the volume of balls centered at : The absence of a growth requirement on in Theorem 1.3 might be partially justified by considering the Laplacian and the mean curvature operator as members of the the family of quasilinear operators for ∈ (ℝ + 0 ) ∩ 1 (ℝ + ) with ′ > 0 on ℝ + , with respectively the choices ( ) = and ( ) = ∕ denote the pre-Legendre transform of , note that realizes a diffeomorphism from ℝ + 0 to [0, ∞ ). If ∞ = ∞, then the integrability qualifying as the right Keller-Osserman condition associated to and Δ writes as see the detailed account in [10,Sec. 10]. Note that (12) rewrites as (7) for the Laplace-Beltrami operator. By [99], in the assumptions of Theorem 1.3 on , (12) is necessary and sufficient for the vanishing of a solution of Δ = ( ), at least in the case (∞) = ∞ (that implies ∞ = ∞, cf. also [10,Remark 10.2]). However, the mean curvature operator does not satisfy ∞ = ∞, and Theorem 1.3 might suggest that a Liouville theorem for solutions of (4), in the spirit of (9), could hold on reasonably well-behaved manifolds , arguably those satisfying (10) or (11), without the need to control the growth of at infinity. However, we shall see below that a suitably defined Keller-Osserman condition will appear again as soon as has a volume growth faster than polynomial. The guiding example is again = ℍ , where the existence problem for (4) was recently studied in [14].

Further ambient spaces, and MCF solitons
Ambient manifolds of the type ℝ × leave aside various cases of interest, notably interesting representations of the hyperbolic space ℍ +1 . This and other examples motivated us to consider graphs in more general ambient manifolds̄ , that foliate topologically as products ℝ × along the flow lines of relevant, in general not parallel, vector fields. A sufficiently large class for our purposes is that of warped products for some positive ℎ ∈ ∞ (ℝ). Note that = ℎ( ) is a conformal field with geodesic flow lines. Given ∶ → ℝ, we define the graph that we call a geodesic graph. For instance, ℍ +1 admits the following representations: (i) ℍ +1 = ℝ × ℝ , with the slices of constant , called horospheres, being flat Euclidean spaces with normalized mean curvature 1 in the direction of − . This follows by changing variables = − log 0 in the upper half-space model (ii) ℍ +1 = ℝ× cosh ℍ , with the slices { = 0 } for constant 0 , called hyperspheres, being totally umbilical hyperbolic spaces of normalized mean curvature = tanh( 0 ) in the direction of − , and sectional curvature −(cosh 0 ) −2 .
In cases (i) and (ii) we say that Σ is, respectively, a geodesic graph over horospheres and hyperspheres. Do Carmo and Lawson in [47] investigated geodesic graphs in ℍ +1 with constant mean curvature ∈ [1, 1], both over horospheres and hyperspheres, and proved the following Bernstein type theorem: if Σ is over a horosphere, then Σ is a horosphere (in particular, = ±1) if Σ is over a hypersphere, then Σ is a hypersphere.
Their result is part of a more general statement, that draws the rigidity of a properly embedded CMC hypersurface Σ from the rigidity of its trace on the boundary at infinity ∞ ℍ +1 . The proof, relying on the moving plane method, seems difficult to adapt to variable mean curvature, and the search for an alternative approach was one among the geometric motivations that lead us to write the monograph [10].
The equations for geodesic graphs with prescribed mean curvature is more transparent when written in terms of the parameter of the flow Φ of , related to by the equation Set so ∶ → and Σ writes as { = ( )}. Clearly, = when̄ = ℝ × with the Riemannian product metric. Computing the normalized mean curvature of Σ with respect to the upward-pointing unit normal we obtain the following equation for ∶ → : where is the derivative of with respect to . Clearly, if = 0 this is an equation of type for suitable , , . Equation (19) will be our focus for the first part of this survey. Another prescription of that leads to an equation of the same type arises from the study of mean curvature flow solitons. Briefly, consider a mean curvature flow (MCF) evolving from Σ, that is, a solution ∶ [0, ] × →̄ of Differentiating at = 0, we see that a MCF soliton satisfies the identitȳ ( , ) = . Interesting solitons arise when is a conformal or Killing field, and in our case the choice = ±ℎ( ) allows to rewrite (18) as the following equation: with = −1 ( ). In particular, a self-translator in ℝ +1 , that is, a soliton for the vertical direction , satisfies We mention that the existence problem for equations tightly related to (20), describing certain -minimal graphs in ℝ × , has recently been considered in [21].
Remark 1.8. Another relevant class of ambient spaces is that of products × ℎ ℝ with metric =̄ + ℎ( ) 2 d 2 , for some 0 < ℎ ∈ ∞ ( ). In this case = is Killing but not parallel, unless ℎ is constant. For instance, ℍ +1 admits two such warped product decompositions, according to whether has one or two fixed points at infinity: the first can be obtained by isolating the coordinate = in the upper half-space model (14), leading to the second can be written as with ∶ ℍ → ℝ the distance from a fixed origin in ℍ , and corresponds, in the upper halfspace model, to the fibration of ℍ +1 via Euclidean lines orthogonal to the totally geodesic hypersphere { 2 0 + | | 2 = 1}. There is no rigidity for graphs in ℍ +1 along these decompositions: the Plateau's problem at infinity for graphs with constant ∈ (0, 1) is always solvable, by work of Guan and Spruck [74] for (23), and Ripoll and Telichevesky [110] for (22). Existence for the prescribed mean curvature equation on more general products × ℎ ℝ is studied in [17].

Analytic behaviour of the mean curvature operator
Hereafter, all solutions of differential equations and inequalities will be meant in the weak sense. In our study of the rigidity problem for equations of the type on a complete, connected manifold , we assume for suitable ∈ [0, 1] and ∈ ℝ, reducing the problem to the investigation of solutions of Note that, for ∈ (0, 1), the RHS is allowed to vanish both as | | → 0 and as | | → ∞.
As we shall see, the arbitrariness in the choice of will allow to reach the optimal thresholds for the rigidity of solutions to (27). Our goal in this section is to identify sharp conditions that force any non-constant solution of (24), possibly already matching some a-priori bound, to satisfy ( ) ≡ 0 on . If has just isolated zeroes, we would directly deduce that must be constant, otherwise the rigidity issue for solutions of requires further tools and assumptions, described in the next section. The argument to deduce ( ) ≡ 0 depends on the fact that any non-constant solution of (26), for suitable or provided that grows slowly enough, is automatically bounded from above and satisfies where, hereafter in this survey, * ∶= sup , * ∶= inf .

Remark 2.1 (Constant solutions). We observe that a constant function solves (26) if and only if either
Therefore, in what follows we concentrate on non-constant solutions.
The conclusion ( * ) ≤ 0 should be interpreted as a maximum principle at infinity, in the spirit of [106]. Indeed, if is compact and we assume ∈ 2 ( ) and = 1, ( * ) ≤ 0 follows readily by evaluating (26) at a point 0 realizing * , and using that | ( 0 )| = 0, 2 ( 0 ) ≤ 0. We shall be more detailed in the next section. To constrain solutions of (24), the idea is to apply the maximum principle at infinity both to and tō ∶= − , which is possible if, for instance, ℱ satisfies conditions like In this case,̄ solves and̄ (̄ * ) ≤ 0 rewrites as ( * ) ≥ 0. If, for instance, is increasing, we would eventually conclude the identity ( ) ≡ 0 on . Taking into account these facts, in the present section we focus on the maximum principle at infinity for solutions of inequality (26), and on related Liouville properties. Our purpose is to guide the reader through the results in [10], discuss the sharpness of the assumptions and their interplay, and illustrate some geometric applications. For reasons of room, no proofs of the analytic results will be given or sketched, but their applications to geometric problems will be commented in some details.

Weak and strong maximum principles at infinity
For bounded from above, the necessity to deduce ( * ) ≤ 0 even if * is not attained motivated the introduction of the fundamental Omori-Yau maximum principles, [103,131,30], that generated an active area of research related to potential theory on manifolds [93,94], with applications to a variety of geometric problems. We refer the reader to [3, 10, 106] for a detailed account and an extensive set of references. In its original formulation, the Omori-Yau principle for the Laplace-Beltrami operator is the property that, for every ∈ 2 ( ) bounded from above, there exists a sequence { } ⊂ satisfying Clearly, if * is attained at some then the sequence can be chosen to be constantly equal to . The Omori-Yau principle is also called, in [3,10], the strong maximum principle at infinity. The fact that the gradient condition | ( )| < 1∕ is unnecessary in many applications motivated the introduction, in [106], of the related weak maximum principle at infinity, asking for the validity, along a suitable sequence { }, just of bounds Examples show that the two are indeed different, that is, there exist manifolds satisfying the weak but not the strong maximum principle, cf. [15,10]. The weak principle turns out to have important links with the theory of stochastic processes, explored in [106] and recalled later in this paper in more detail. In our setting, we shall look at suitable formulations of weak and strong principles, adapted to solutions with low regularity and to mean curvature type operators. Extension to more general operators can be found in [10] (quasilinear ones) and in [93,94] (fully nonlinear ones).

Remark 2.3.
If is complete, by Ekeland's principle the set { > , | | < } is non-empty for each choice of < * and > 0. The restriction ∈ 1 ( ) is made necessary for such set to be open. Suitable notions of (SMP ∞ ) for viscosity solutions that are less regular than 1 have been studied in [93,94] in the setting of fully nonlinear potential theory.
Remark 2.4. For = 1 and ∈ 2 ( ), then (SMP ∞ ) can equivalently be formulated as the existence of a sequence { } ⊂ such that making contact with the original formulation of Omori and Yau.
The definition of (WMP ∞ ) is localized on upper level sets { > } of , for ∈ ℝ. However, with the aid of the Pasting Lemma, it can easily be seen that the following statements are equivalent (cf. [10,Prop. 7.4]): on that is bounded from above satisfies ( * ) ≤ 0.
The applications described below will motivate the study of the above maximum principles at infinity. In particular, we mention that although the (WMP ∞ ) can be successfully used to study the rigidity of entire minimal graphs and MCF solitons, to obtain sharp results for hypersurfaces satisfying (18) with more general choices of we had to appeal to the full strentgh of the (SMP ∞ ). Examples show that some restriction on the geometry of has to be required in order to deduce the validity of (WMP ∞ ) and (SMP ∞ ). Concerning the weak maximum principle, we prove in [10, Thm. 7.5] the following sharp Theorem 2.5 (Thm. 7.5 in [10]). Let be a complete Riemannian manifold, and let If then (WMP ∞ ) holds. In particular, for ∈ (ℝ), any non-constant solution ∈ Lip loc ( ) of that is bounded from above satisfies ( * ) ≤ 0.
Concerning the strong maximum principle at infinity, it is still an open problem to decide whether the volume growth (32) is sufficient. In [10], we obtained sharp conditions by imposing the next lower bound on the Ricci curvature: for constants ≥ 0, ≥ −2, where is the distance from a fixed origin . Namely, we compare to the radially symmetric model being ℝ endowed, in polar coordinates ( , ), with the metric where d 2 is the metric on the unit sphere −1 and ∈ ∞ (ℝ + 0 ) is a solution of The case = 0 corresponds to choosing the Euclidean space as model , described by the function ( ) = , while the case = 0 corresponds to choosing the Hyperbolic space of curvature − 2 , for which By Proposition 2.11 in [107] and the Bishop-Gromov comparison theorem, the volume of the geodesic balls and in and , respectively, relate as follows: if > −2, as → ∞, wherē = (1 + √ 1 + 4 2 )∕2. Therefore, the inequalities imply any of the conditions (69) with > 0: respectively if either > −2 or = −2. Thus, the growth thresholds in the following (sharp) sufficient condition for (SMP ∞ ) in [10, Thm. 8.5] turn out to be the same as for (WMP ∞ ).

Application: minimal and prescribed mean curvature graphs
We present the following applications of the weak and strong maximum principles at infinity, that generalize Do Carmo-Lawson's theorem stated in the introduction. We begin with entire minimal graphs, proving under ( ), there exists no entire minimal graph over .
Remark 2.8. The integrability conditions in ( ), ( ) mean, loosely speaking, that the warped product̄ opens up faster than linearly for large . The case of cones ℎ( ) = is thus borderline. To recover the results in [47], the choice = ℍ and ℎ( ) = cosh leads to the constancy of entire minimal graphs over hyperspheres in ℍ +1 , while the choice = ℝ and ℎ( ) = allows us to recover the non-existence of entire minimal graphs over horospheres.
Concerning graphs with non-constant mean curvature, (WMP ∞ ) seems to be not sufficient to conclude rigidity for the more involved equation (18), and we use the full strength of (SMP ∞ ) to obtain the following Theorem. For simplicity, we state the result for the warping function ℎ( ) = cosh , but the proof can easily be extended to each ℎ ∈ 2 (ℝ) satisfying for some constant > 0. Fix a constant 0 ∈ (−1, 1), and consider an entire geodesic graph of ∶ → ℝ with prescribed curvature ( ) ≥ − 0 in the upward direction. Then, is bounded from above and satisfies * ≤ arctanh( 0 ).

Remark 2.10.
In particular, the choice = ℍ enables us to recover the statement of Do Carmo-Lawson's result [47] for CMC graphs over hyperspheres.
Proof: sketch. Defining , ( ) and ( ) as in (15) and (16) If, by contradiction, the following upper level set of (hence, of ) is non-empty for some ∈ ℝ + : then, for > 0 is small enough, on the (non-empty) set Since is bounded above, applying Theorem 2.6 with = 1, = 0, = 1 we reach the desired contradiction.
To prove ( ), suppose that | ( )| ≥ 1 on . Since cosh is even, the graph of − has curvature − ( ) in the upward direction. Thus, up to replacing with − we can suppose that ( ) ≥ 1. Applying the first part of the theorem to any 0 > −1 we obtain * ≤ arctanh( 0 ), and the non-existence of follows by letting 0 → −1.
Our results also apply to Schwarzschild, ADS-Schwarzschild and Reissner-Nordström-Tangherlini spaces, considered in [37], that we now briefly recall. Having fixed a mass parameter > 0 and a compact Einstein manifold ( , ) with Ric = ( − 1) , the Schwarzschild space is the product where and 0 ( ) is the unique positive root of ( ) = 0;̄ generates a Lorentzian manifold ℝ ×̄ with metric − ( )d 2 + ⟨ , ⟩ that solves the Einstein field equation in vacuum with zero cosmological constant. Similarly, given̄ ∈ {1, 0, −1} the ADS-Schwarzschild space (with, respectively, spherical, flat or hyperbolic topology according to whether̄ = 1, 0, −1) is the manifold (39) with respectively and 0 ( ), as before, the unique positive root of ( ) = 0. They generate (static) solutions of the vacuum Einstein field equation with negative cosmological constant, normalized to be − ( + 1)∕2. The change of variables allows to write 1̄ = ℝ + × ℎ . Similarly, the Reissner-Nordström-Tangherlini space is a charged black-hole solution of Einstein equation, described by (39) with the choice with ∈ [− , ] being the charge. We consider solitons with respect to the conformal field A direct application of Theorem 2.5 shows the following: Theorem 2.11 (Thm. E in [37]). There exists no entire graph in the Schwarzschild and ADS-Schwarzschild space (with spherical, flat or hyperbolic topology), over a complete , that is a soliton with respect to the field √ ( ) . Differentiating, If is compact, this is sufficient to guarantee that (43) does not admit solutions, by the classical maximum principle, and settles the case of Schwarzschild and ADS-Schwarzschild spaces with spherical topology. In the remaining cases, observe that ( ) ∼ for large , thus in view of (42) ℎ( ) = ( ) ∼ as → ∞.
As a consequence, the flow parameter ( ), hence , is bounded from above. The Einstein condition and Bishop-Gromov comparison theorem imply To conclude, we apply Theorem 2.5 to deduce ( * ) ≤ 0 and to reach the required contradiction.

Maximum principles at infinity for solutions with controlled growth
Theorem 7.5 in [10] also applies to solutions with a controlled growth at infinity.
for some constant > 0. Let ∈ Lip loc ( ) be a non-constant solution of such that for some > 0 satisfying If one of the following properties hold: then, is bounded above on and ( * ) ≤ 0.
for some explicit constant ≥ 0 that vanishes in various cases of interest. Note that, for instance, (49) readily imply ≤ 0 whenever (46) and > 0 hold.
Remark 2.14 (The limiting case = 0). This corresponds to the most "singular" case, and sometimes it requires an ad-hoc treatment. Conditions (47) and (48) rephrase for = 0 in the following equivalent form: for somē > 1. Furthermore, because of [10, Rem. 7.6], the conclusion of the theorem holds under the weaker assumption We underline that ( ) guarantees the absence of solutions with polynomial growth to Remark 2.15 (Sharpness). In Subsection 7.4 of [10], we show that the radially symmetric model described in Subsection 2.1 admits an unbounded, radial solution ∈ Lip loc ( ), increasing as a function of the distance from the origin and solving that coincides with for ≥ 2 whenever ∈ [0, 1], ≤ + 1 and one of the following conditions hold: while it coincides with ( ) ∕ log for ≥ 2 when either The existence of under any of 1), … , 6) above shows the sharpness of the parameter ranges (44) and (47), and of the growth conditions (46) for and (48) for | |. In particular, -in (2), all the assumptions are satisfied but the second or third in (48), where the liminf is ∞, while in (3), the liminf in the third in (48) is finite but bigger than the threshold 2 − for ≤ 1; -in (4) the requirements in the first of (48) are all met, but + = ( ) instead of + = ( ).

Application: Solitons for the Mean Curvature Flow
Our first result considers solitons for the MCF in ℝ × that are realized as the graph of ∶ → ℝ, and move in the direction of the parallel field (called graphical self-translators). In Euclidean space ℝ +1 = ℝ × ℝ , the bowl soliton (cf. [2] and [32,Lem. 2.2]) and the nonrotational entire graphs in [134] for ≥ 3, provide examples of entire, convex graphs that are self-translators and grow of order 2 at infinity. In [10, Thm.2.28], we prove that the growth 2 is sharp for a much larger class of ambient spaces.
Then, there exists no graphical self-translator ∶ → ℝ with respect to the vertical direction satisfying Proof: sketch. It is a direct application of Theorem 2.12 ( > 0) or Theorem 2.5 ( = 0), The generality of Theorem 2.12, and of its proof, allow for applications to self-translators in Euclidean space whose translation vector field differs from the vertical field . The next result relates to [4], where the authors proved that hyperplanes are the only complete selftranslators whose Gauss image lies in a spherical cap of of radius < ∕2: it is of interest also because it shows the sharpness of Theorem 2.12, in its more general form described in Remark 2.13, even with respect to the constant in (49).
Then, Σ cannot be a self-translator with respect to any parallel field whose angle ∈ (0, ∕2) with the horizontal hyperplane ℝ satisfies tan >̂ .
In particular, if̂ = 0, then Σ cannot be a self-translator with respect to a vector which is not tangent to the horizontal ℝ .
Proof: sketch. Up to a rotation of coordinates on ℝ and a reflection in the ℝ-axis, we can assume that ⟨ , ⟩ > 0 and = cos 1 + sin , where 1 is the gradient of the first coordinate function 1 . By (17) and (18), the soliton equation = ⟨ , ⟩ satisfied by ( ) = ( ( )) = ( ) can be written as The operator in the left-hand side is self-adjoint on ℝ endowed with the weighted mea- Assume by contradiction that (57) holds, in particular, is non-constant. Applying Theorem 2.12 to (58) (in the strengthened form described in Remark 2.13) with the choices = sin , = 1, = 1, = 0, we conclude that sin ≤ ̂ , for some explicit constant ≥ 0. The value of , that also depends on (59), becomes = cos . Therefore, we deduce that̂ ≥ tan , which is a contradiction.
With a very similar proof, we can also consider graphical self-expanders for the MCF, i.e. graphs that move by MCF along the integral curves of the position vector field Proof. Set ( ) = | | for ∈ ℝ , and write the position field along Σ as follows: The soliton equation = ( , ) satisfied by ( ) = ( ( )) = ( ) becomes The only constant solution is ≡ 0. By contradiction, assume the existence of a non-constant solution . An explicit computation shows that lim inf →∞ log vol 2 ∕2 ( ) 2 < ∞.
Applying Theorem 2.12, adapted to weighted volumes, to both and − , with the choices = 1, = 0, we get 0 ≤ ≤ 0, which is the final contradiction.

Keller-Osserman condition and Liouville theorems
To obtain the rigidity of solutions of the prescribed mean curvature equation without assuming an a-priori control on , on general manifolds one needs to grow sufficiently fast at infinity. A notable exception is Theorem 1.3 on ℝ , of which we now comment an improved version; in [10, Thm. 10.30], we modified the original capacity argument in [121,128] to apply to solutions of provided that has a polynomial volume growth, and we prove the next result. Note that the inequality corresponds to the limiting case = 0 in (45). We mention that the method is successfully used for more general equations, in particular we refer to works of Mitidieri and Pohozaev [97], D'Ambrosio and Mitidieri [38,39], Farina and Serrin [60,61].
admitting unbounded solutions to (61), for a positive constant, with arbitrarily small polynomial growth. The positivity of the limit in (62) is also sharp, because otherwise, by Theorem 2.12 and Remark 2.14, solutions with polynomial growth to (61) with a positive constant do not exist.
As a consequence of Theorem 2.20 and of Bernstein Theorem 1.1, we directly deduce the following Hence, gives rise to a minimal graph in ℝ × , that is one sided bounded because ≢ 0. Theorem 3.2 guarantees that ≡ is constant, and plugging in the equation we get ( ) = 0.
As Remark 2.22 shows, when the growth of | | is faster than polynomial one needs some growth conditions on to guarantee that solutions of on are bounded from above and satisfy ( * ) ≤ 0. In [10], for the class of inequalities and when ′ ( ) considering the homeomorphism ∶ [0, ∞) → [0, ∞) given by we identify the following Keller-Osserman condition: with as in (5). Observe that (65) does not depend on the growth/decay of , but the latter provides the thresholds to guarantee that (65) is sufficient to ensure that solutions of (63) are bounded from above and satisfy ( * ) ≤ 0. The meaning of the bounds on is described in detail in [10]. However, for the mean curvature operator and in our case of interest still (∞) < ∞. The proposal in [91] to modify (64) for the mean curvature operator by setting turns out to lead to sharp results, and in our prototype case (66) for some > , and assume that ∈ Lip loc ( ) is a non-constant solution of If either  [130]. (Sharpness of > ). We consider the radially symmetric model in Subsection 2.1. Because of the asymptotic growth of log | | in (35), given ∈ [0, 1], the condition

Remark 2.26
while, for ∈ (0, 1], condition The case = 0 in the second of (69) will not be considered in the present counterexample.
Keeping this in mind, we prove the optimality of condition > . The computations in Subsection 10.4 of [10] (cf. also Remark 10.29 therein) show that the radially symmetric function constructed in Remark 2.15, that is positive and coincides with for large , solves provided that ∈ [0, 1] and either Hence, > is necessary for the conclusion in Theorem 2.24. The last restriction in the second of (71) is also optimal: in the "Euclidean setting" = −2, if = and < + 1 then entire solutions of (70), with the equality sign, are constant if they have polynomial growth, see [61,Thm. 12] and also Example 4 at p. 4402 therein.
To extend the above theorem to more general satisfying the Keller-Osserman condition (67), we had to use a different approach inspired by the classical Phrágmen-Lindelöff method, see [10, sec. 10.3]: under the validity of (67), the key point is to construct a suitable family of radial blowing-up solutions that locally converge to a constant function as → 0. The family { } is used to constrain via the use of a comparison principle. The presence of the gradient term in the RHS of (68) requires to control the gradient of and in a neighbourhood of the set of maxima of − , and therefore we need to restrict ourselves to ∈ 1 ( ).
The construction of is rather delicate, especially in the effort to avoid unnecessary technical conditions. It requires to estimate Δ from above, and therefore, by comparison theory, it needs a control on the Ricci curvature from below. We require the same bounds as those in (2.6), and prove the following result that, because of the discussion at the end of Subsection 2.1, well fits with Theorem 2.24.
for constants ≥ 0, ≥ −2, with the distance from a fixed origin . Choose and let ∈ (ℝ) be non-decreasing. Then, under the validity of the Keller-Osserman condition with as in (5), any non-constant solution ∈ 1 ( ) of is bounded above and satisfies ( * ) ≤ 0.

Application: the capillarity equation
Applying Theorem 2.24 to the capillarity problem, we deduce the following corollary. To our knowledge, this is the first result on the capillarity equation that allows to have geodesic balls with growth of their volume faster than polynomial.
for some < 2. If there exists > 0 such that then the only solution of the capillarity equation Proof: sketch. Evidently, ≡ 0 is the only constant solution. Assume the existence of a nonconstant solution . Up replacing with − , we can assume that { > 0} ≠ ∅. By assumption, there exists < 1 sufficiently close to 1 in such a way that

The capillarity equation implies
Therefore, + ∶= max{ , 0} is a non-constant solution of div where we used that 1− ∕ √ 1 + 2 is bounded on ℝ. The inclusion of the "artificial" gradient term enables the function ( + ) = + to satisfy the Keller-Osserman condition 1 > , and therefore we can apply Theorem 2.24 to obtain * + ≤ 0, contradiction.

Gradient estimates and Bernstein theorems in ℝ ×
As we saw above, under general assumptions on , , solutions of However, to infer the constancy of solving (75), one needs rather different arguments (and more binding assumptions) than those leading to the results in the previous section; indeed, while the above theorems apply to differential inequalities, those in this section are very specific to the equality case, unless in the special situation when has slow volume growth. This consideration is not surprising, as it parallels the case of harmonic functions: positive solutions of Δ = 0 are constant on each complete manifold with Ric ≥ 0 by [131,30], while the constancy of every positive solution of Δ ≤ 0 is equivalent to the parabolicity of (cf. [119,72]) that, for complete manifolds with Ric ≥ 0, is equivalent to the slow volume growth condition As a matter of fact, (76) is sufficient for the parabolicity of a complete manifold , regardless to any curvature condition (cf. [72,Thm. 7.5]). The next example is instructive:  (1). Denoting with the distance from an origin (0, ), observe that for some constant > 0. Furthermore, is bounded on if and only if ≥ 3. Therefore, if ≥ 3, assumptions (77) are not sufficient to deduce the constancy of positive solutions of (75), not even of bounded ones. As shown in [22], if the first in (77) is strengthened to having asymptotically non-negative Sectional curvature, that is, if and if has only one end, then positive solutions of (75) that are bounded on one side are constant.
The situation for manifolds with slow volume growth, including the above example when = 2, is quite rigid. Suppose that is a positive solution of (75) on , and assume that satisfies (76). Let be the metric on the graph Σ. By a calibration argument that can be found in [88], the volume of the portion of (Σ, ) inside an ambient geodesic ball ⊂ × ℝ centered at ( , ( )) satisfies with the geodesic ball of radius in ( , ) centered at . In particular, since the geodesic ball in (Σ, ) is contained in Σ ∩ , its volume satisfies (76). Hence, Σ is parabolic. The constancy of follows since > 0 solves Δ = 0 on Σ, with Δ the Laplace-Beltrami operator of . Summarizing, A complete manifold with slow volume growth, in the sense of (76), satisfies (ℬ2): entire positive minimal graphs over are constant.
Inequality (79) allows to "transplant" a volume condition on to a corresponding one on Σ, in particular condition (76) that suffices for parabolicity, without any a-priori bounds on . A natural question that arises is whether there is a way to transplant the parabolicity itself from to Σ. Currently, we are not aware of any example of entire, minimal graph over a complete, parabolic manifold that is not parabolic.
If the volume growth of ( , ) does not satisfy (76), the main tool to prove the constancy of solutions of (75) is a local gradient bound for positive solutions ∶ ⊂ → ℝ + , as in [11]. However, the possible lack of a Sobolev inequality on manifolds satisfying Ric ≥ −( − 1) 2 (or Sec ≥ − 2 ) forced to devise a method different from the one in [11]. Typically, in a manifold setting the estimate is achieved by refining the classical idea of Korevaar [82], as for instance in [117]. There, the authors inferred the following bound for positive solutions ∶ ⊂ → ℝ under the requirements that Ric ≥ 0 and Sec ≥ − 2 on : for some = ( ). Further improvements can be found in [40,41,22]. The proof of the constancy of proceeds by first showing the global bound | | ≤ on , for some constant > 0, that allows to interpret the mean curvature operator as a uniformly elliptic linear operator on ( , ) acting on . At this stage, for instance, the quadratic lower bound (77) on Sec suffices to show | | ≤ for every positive minimal graph with at most linear growth, [22]. Next, assumptions like Ric ≥ 0, or Sec satisfying (78), and possessing only one end, as in [22] serve to guarantee the validity of the global Harnack inequality sup ≤ ′ inf , from which the constancy of follows. We mention that, recently, in [44] the author proved Theorem 1.1 on each manifold supporting global doubling and a weak Neumann-Poincaré inequalities, a class that includes manifolds with Ric ≥ 0. His idea is to avoid the need of gradient estimates by transplanting the doubling and weak Neumann-Poincaré inequalities (the latter, only for suitable functions) from to the graph Σ, and then proving the Harnack inequality for the harmonic function on Σ.

Ricci curvature and gradient estimate for minimal graphs
Our approach in [35] to obtain Theorem 1.1 still relies on a gradient estimate. Precisely, we prove Then In the particular case Ω = , If equality holds in (82) at some point, then = 0 and is constant.
Despite we found no explicit example, we feel likely that the bound (82) be sharp also for > 0, in the sense that the constant √ − 1 cannot be improved. Our estimate should be compared to the one for positive harmonic functions on manifolds satisfying (80), obtained by P. Li and J. Wang [87] by refining Cheng-Yau's argument: and its version for sets with boundary in [92,Thm. 2.24]: Indeed, to reach our goal we reworked Korevaar's technique, that parallels Cheng-Yau's estimate though with some notable difference. We also tried to simplify the method and highlight its key features, in such a way to obtain neat constants in (81). Writing ∇, ‖⋅‖ for the connection and norm in the graph metric on Σ, and letting denote the upward pointing unit normal of Σ, the Jacobi equation for −1 = ⟨ , ⟩ leads to the identity where ‖II Σ ‖ 2 is the second fundamental form of Σ, Ric is the Ricci curvature of × ℝ, and ℒ is the weighted Laplacian The lower bound on the Ricci curvature implies so the negativity of Ric can be absorbed by defining and noting that satisfies However, from (86) it is not clear that satisfies a maximum principle on Ω if Ω is unbounded. Note that (86) somehow corresponds, in the case of harmonic functions, to the equation for = | |∕ = | log | (cf. [87]): for a suitable weighted Laplacian ℒ ′ . The term with̄ 3 appears because of the refined Kato inequality for solutions of Δ = 0, and allows to satisfy the Keller-Osserman condition (7). Although Cheng-Yau-Li-Wang's method proceeds by localizing (87) with suitable cut-offs, the validity of (83), equivalently, (sup ̄ ) ≤ 0, is expectable also in view of Theorem 2.24 above. Indeed, integral estimates close to those in Theorem 2.24 were used, in place of localization tecniques, to prove the generalization of (83) and (84) for -harmonic functions without the need to bound the sectional curvature of , see [127,92]. On the contrary, in (86) it seems difficult to relate, if possible, the refined Kato inequality for ‖∇ 2 ‖ 2 = −2 ‖II Σ ‖ 2 to in order that the resulting matches (7). Therefore, to overcome the problem we localize as in Korevaar's method:

Sketch of the proof of Theorem 3.2. Fix
We claim that the following set is empty for every choice of constant > 0: Once the claim is shown, (81) follows by letting → 0 and then ↓ √ − 1. By contradiction, assume that Ω ′′ ≠ ∅, and let , ′ , be positive small numbers to be chosen later. Consider a non-negative exhaustion on Ω ′′ (that is, { ≤ } is compact in Ω ′′ for every constant ), and define is non-empty. Computing ℒ 0 we get and moreover where we used Young inequality with free parameter ′ > 0 to be specified later. In [117,22,40,41], at this stage is chosen to be the distance function on ( , ) from a fixed origin.
Since the components and of −1 and −1 are related by then we deduce ‖∇ ‖ 2 ≤ | | 2 = 1, while because of minimality with the components of 2 . Since the eigenvalues of with respect to are 1 and −2 , that might be very different, to estimate Δ an information on the full Hessian of the distance function (hence, an assumption involving Sec) seems unavoidable. Indeed, one can weaken the bound on Sec to a bound involving an intermediate Ricci curvature, but a control just on Ric does not suffice. Our first observation is that, in (89), the quantities Δ and ‖∇ ‖ 2 need not be estimated separately, but the last term in the right-hand side will be small for suitably small << ′ provided that We recognize that (90) is satisfied by the logarithm of a solution of which makes contact with the equations studied in the previous section, and in fact suggests to produce by mean of potential theory. We use the information that Σ is calibrated, so one can control the volume of balls in Σ in terms of the measures of Ω and Ω in : by modifying the method leading to (79), we prove in [35,Lem. 1,2] that Ω ′′ can be isometrically embedded as an open subset of a complete manifold ( , ℎ) whose geodesic balls ℎ satisfy, under either of assumptions ( ) and ( ) in Theorem 3.2, Therefore, by [106,Thm.4.1] the weak maximum principle holds for Δ ℎ on ( , ℎ), namely, for ∈ (ℝ) every solution of Δ ℎ = ( ) that is bounded from above satisfies ( * ) ≤ 0.
In particular, the only non-negative solution of By [106], the validity of the weak maximum principle for Δ ℎ turns out to be equivalent to the stochastic completeness of ( , ℎ), that is, to the property that the minimal Brownian motion on is non-explosive. Since (92) implies then the stochastic completeness of also follows by [72,Thm. 9.1].
Next, we use the AK-duality principle (Ahlfors-Khas'minskii duality) recently established in [95,93,94]: in our case of interest, the AK-duality guarantees that the weak maximum principle is equivalent to the possibility to produce solutions of (91) on ( , ℎ). We mention that the AK-duality works in greater generality for subsolutions of fully nonlinear PDEs of the type, under mild assumptions on ℱ, and states the equivalence between a maximum principle on unbounded open sets for solutions of (94), and the existence of suitable families of exhaustions that solve ℱ( , , d , While, for general ℱ, the construction of is subtle and proceeds by stacking solutions of obstacle problems, in our case the linearity of Δ ℎ allows to produce in a much simpler way: fix a small relatively compact ball and an exhaustion Ω ↑ by means of relatively compact, smooth open sets, and for each let solve The maximum principle implies that 0 < < 1 on Ω ∖ , thus setting = 1 on ∖Ω produces a locally Lipschitz, viscosity solution of Δ ℎ ≤ on ∖ . By elliptic estimates, ↓ locally smoothly on ∖ , for some 0 ≤ ∈ ∞ ( ∖ ) solving Since ( , ℎ) is stochastically complete, necessarily ≡ 0. Thus, given ∈ ℕ there exists ( ) such that ( ) < 2 − on Ω . The series converges locally uniformly to a locally Lipschitz exhaustion satisfying Δ ℎ ≤ on ∖ .
The full set of properties in (91) can be matched by extending in , then suitably translating and rescaling , and eventually smoothing by using Greene-Wu approximation techniques [73].
Having shown the existence of in (90), inserting into (89) and using that We infer the existence of ′ sufficiently small, and of ( ′ , ) << ′ small enough that and thus, up to further reducing ′ , we conclude from (88) a contradiction that concludes the proof.

Ricci curvature and gradient estimate for CMC graphs
The above method allows for generalizations to the CMC case, obtained in the very recent [36], that apply to the rigidity of capillary graphs over unbounded regions. We detail the application in the next section. Although the guiding idea is the same as the one for Theorem 3.2, the difficulty to deal with nonzero makes the statement of the next result, and its proof, more involved. Then, In particular, if Ω = , Remark 3.5. Because of the non-existence result of Heinz [76], Chern-Flanders [31,66] and Salavessa [118] for entire CMC graphs with nonzero on manifolds with subexponential volume growth, if ≠ 0 the bound for solutions on the entire is meaningful only if > 0.
3. if > 0 and ≥ 0, we can choose Note that, in the limit → 0, the above choices allow us to recover inequality (81) for minimal graphs.
It is interesting to compare our estimates with the known restrictions on the mean curvature of entire CMC graphs. By [118], integrating on an open set Ω ⋐ the equation for the prescribed mean curvature, one has in particular, up to replacing with − , | | does not exceed the Cheeger constant of : By Cheeger's inequality [23], the bottom of the spectrum of Δ on satisfies 1 ( ) ≥ ( ) 2 ∕4. On the other hand, Brooks upper bound [16] for 1 ( ) on manifolds satisfying Ric ≥ −( − 1) 2 , guarantees that Summarizing, if Ric ≥ −( − 1) 2 then any entire cmc graph must satisfy In particular, the second case in (99) cannot happen for entire graphs. The bound in (100) is optimal since, in the hyperbolic space ℍ of curvature − 2 , Salavessa in [118] constructed a positive, radial, entire graph with mean curvature for each choice of ∈ (0, −1 ]. Furthermore, the graph has bounded gradient if < −1 .

Splitting for capillary graphs
In the previous sections, we gave conditions to ensure that solutions of (19) satisfy ( ) ≡ 0, and in some cases that they are constant. In this section, we give conditions to guarantee that the existence of a non-constant ∶ → ℝ solving the prescribed mean curvature problem forces to split as a product × ℝ, with only depending on the ℝ-variable. The use of solutions of PDEs to prove splitting type results is classical in Geometric Analysis, and dates back at least to the theorem of Cheeger-Gromoll [29] for manifolds with Ric ≥ 0 and more than one end. In the same period, Schoen and Yau [119] investigated the topology of complete, stable minimal hypersurfaces without boundary Σ ↪̄ +1 in a complete manifold with nonnegative sectional curvature, and discovered the following integral inequality for harmonic functions ∶ Σ → ℝ: (see also [89] for another proof). The formula was later used by Li and Wang [89] to prove that any such hypersurface must split as a Riemannian product × ℝ, provided that it is properly immersed and has more than one end. The elegant techniques were based on the theory of harmonic functions developed in [85], and relate to their previous splitting theorems for manifolds with Ric ≥ −( − 1) 2 and positive spectrum [86,87].
In [56], motivated by De Giorgi and Gibbons' conjectures for the Allen-Cahn equation Farina independently discovered an integral formula, tightly related to (102), valid for each quasilinear operator. He used the integral inequality to establish the 1D-symmetry of certain Here, stability means the the linearized operator at is non-negative, and is said to have 1Dsymmetry if there exists ∈ −1 such that ( ) =̄ ( ⋅ ). However, he didn't publish the results up until 2008, in a joint paper with Sciunzi and Valdinoci [59]. For the mean curvature operator, the formula reads as follows [59,Thm. 2.5]: where, at a point such that ( ) ≠ 0, ⊤ and are the gradient and second fundamental form of the level set { = ( )}. They deduced the next theorems: if ∈ Lip loc (ℝ), stable solutions ∈ 2 (ℝ ) of (101)  The study of capillary graphs in ℝ × motivates the introduction of boundary conditions in (101), and provide a natural setting where the results in the previous sections apply. A capillary hypersurface Σ in a Riemannian manifold with boundarȳ is a CMC hypersurface that meets ̄ at a constant angle. If the angle is ∕2, Σ is said to have free boundary in̄ . Capillarity hypersurfaces naturally have a variational interpretation, see [64,115]. Suppose that Σ is embedded, connected and transversal to ̄ along Σ. Fix a connected, relatively compact open set ⋐̄ in such a way that Σ divides into two pieces and , and note that ∩ is the union of a portion of Σ and of a relatively open subset Ω ⊂ ̄ . Then, Σ is a stationary point for the functional |Σ ∩ | − cos | |, ∈ 0, 2 with respect to variations that are compactly supported in and preserve the volume of , if and only if Σ is a CMC hypersurface satisfying where and̄ are, respectively, the outward normals to Σ ↪ Ω and Σ ↪ Σ.
In [64], Finn studied in depth various aspects of the capillarity problem. However, to the best of our knowledge, rigidity issues have been investigated only in compact ambient spaces̄ with a large amount of symmetries. Especially, capillary hypersurfaces in the unit ball̄ = +1 attracted the attention of researchers, also in view of the link between free boundary minimal hypersurfaces in +1 and the Steklov eigenvalue problem (cf. [68,69] and the references therein). We mention that stable capillary hypersurfaces in +1 have been classified in [116,102,5] (free boundary case), and [115,132], to which we refer for a detailed account. We here address the noncompact, graphical case, so in our setting, we let̄ = Ω×ℝ + 0 for some smooth open domain Ω ⊂ , and we consider Σ to be the graph of ∶ Ω → ℝ + 0 satisfying = 0 on Ω. Having fixed ⋐̄ , we let denote the subgraph of Σ in . Then, Σ is a capillary hypersurface if and only if satisfies the overdetermined boundary value problem for some constant ∈ ℝ, with the unit exterior normal to Ω. Problems like (104) stimulated a vast research over the past 30 years, originally in the semilinear case with ∈ Lip loc (ℝ), a major issue being to characterize domains supporting a bounded solution of (105) (that, agreeing with the literature, we call -extremal domains). Rigidity for relatively compact Ω follows from pioneering works of Serrin and Weinberger [120,135], while the noncompact case has been investigated in the influential [7,6] by Berestycki, Caffarelli and Nirenberg. In [7], the authors conjectured that the only connected, smooth -extremal domains of ℝ with connected complementary are either the ball, the half-space, the cylinder × ℝ − or their complements. The conjecture turns out to be false in full generality, by counterexamples in [122] ( ≥ 3) and [113] ( = 2 and Ω an exterior region). Surprisingly, in ℝ 2 the conjecture is true if both Ω and its complementary are unbounded [112]. Corresponding rigidity results for -extremal domains in ℍ 2 and 2 have been obtained in [50,51].
Typically, the characterization of domains Ω ⊂ ℝ supporting a solution of (105) proceeds by first showing that is monotone in one variable, namely > 0, which is the case for suitable (for instance, derivatives of bistable nonlinearities, [7]). Since satisfies the linearized equation Δ + ′ ( ) = 0, its positivity implies that is stable. In [62], Farina and Valdinoci made the remarkable discovery that the identity corresponding to (103) for Δ + ( ) = 0, meaning remains true even if the support of contains a portion of Ω, provided that satisfies (105) and > 0. In other words, the contributions of boundary terms to (106) on Ω∩spt vanish identically if both and are constant on Ω, a fact that enables to use (106) to characterize domains supporting a solution of (105). This point of view has been further developed in [58] on manifolds with Ric ≥ 0, to obtain splitting results for domains supporting a solution of (105) that is monotone in the direction of a Killing field.
In the literature, overdetermined problems for more general operators have mainly been investigated in the -Laplacian case, leaving the corresponding questions for the mean curvature operator, to the best of our knowledge, mostly unexplored. Taking into account that (104) alone imposes restriction on the geometry of Ω even without overdetermined boundary conditions (cf. [90,33,80], and the survey in [133]), we could expect more rigidity than in the case of the Laplacian. However, as far as we know there is still no attempt to study the equivalent of the Berestycki-Caffarelli-Nirenberg conjecture for (101). We here comment on the very recent [36], where we address this conjecture for constant , proving the following splitting theorem for capillary graphs on manifolds with non-negative Ricci curvature. The result also relates to the work in progress [57], where the authors study the problem for general , with emphasis on bistable nonlinearities. To state the theorem, we recall that a manifold with boundary ( , ℎ) is said to be parabolic if, for each compact subset ⊂ , the capacity cap( ) ∶= inf ∫ | | 2 d ∶ ∈ Lip ( ), ≥ 1 on vanishes. It can be checked that the condition corresponds to the recurrency of the minimal Brownian motion on that is normally reflected on . By [80], the parabolicity of ( , ℎ) is equivalent to the validity of the following principle: with the outward pointing normal to . If = Ω and Ω is an open subset of a complete manifold , by [80] a sufficient condition for Ω to be parabolic is that with the metric ball of radius in centered at a fixed origin of . By an application of Hölder inequality (cf.  for some constant > 0. Split Ω into its connected components { Ω}, 1 ≤ ≤ 0 , possibly 0 = ∞. Let ∈ 2 (Ω) be a non-constant solution of the capillarity problem for some set of constants , , ∈ ℝ, with the exterior normal to Ω, { } a bounded sequence if 0 = ∞, and with the agreement 1 ≤ 2 ≤ … ≤ 0 . Assume that either is constant on Ω, or that If there exists a Killing vector field on Ω with the following properties: sup Ω | | < ∞, ( , ) ≥ 0 on Ω, for every , then: (i) Ω = (0, ) × with the product metric, for some ≤ ∞ and some complete, boundaryless with Ric ≥ 0.
-if > 0, then 1 ≤ 0 and in particular, Furthermore, if is constant on the entire Ω, then the requirement (109) can be omitted and the case = 0 in (iii) occurs.
We stress that is just required to satisfy a mild relation, to make it coherent with the Neumann data: precisely, shall point outwards on components Ω for which the outward pointing normal derivative of is positive (i.e., > 0), leaving the possibility that be tangent to Ω, while shall point inwards if < 0. No condition is imposed on components of Ω where = 0. However, to avoid a degenerate case, we require that is not tangent to the entire portion of Ω where ≠ 0. To conclude, we briefly describe the proof of Theorem 4.2. The argument follows the lines of [58,Thm.5], with some notable improvements depending on the gradient estimate in Theorem 3.4, and can be divided into the following steps.
The conclusion̄ > 0 follows from the strong maximum principle. In particular, d ≠ 0 on Ω.
-We show the validity of the following integral formula on Σ: for every ∈ Lip (Σ), where, at a point such that d ( ) ≠ 0, ∇ ⊤ and are the gradient and second fundamental form of the level set { = ( )} in Σ, and integration is with respect to the volume measure of Σ. Inequality (113) is tightly related to (102) and (106); in particular, the "magic cancellation" of the boundary terms in (106) still holds for (108). As shown in [57], in fact, overdetermined problems associated to a large class of quasilinear operators lead to integral formulas like (113).
-We use again the gradient bound | | ≤ to deduce that for some constant ′ . Therefore, the parabolicity of (Ω, ) enables to deduce the existence of a sequence ↑ 1 such that so taking limits we deduce By adapting [58], the first two conditions guarantee that Ω splits as a Riemannian product along the flow lines of ∇ (everywhere nonvanishing on Ω), and that only depends on the split direction. Also, the constancy of ‖∇ ‖ ∕̄ rewrites as the constancy of ( , ) on Ω.