Compact anisotropic stable hypersurfaces with free boundary in convex solid cones

We consider a convex solid cone C⊂Rn+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {C}\subset \mathbb {R}^{n+1}$$\end{document} with vertex at the origin and boundary ∂C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\partial \mathcal {C}$$\end{document} smooth away from 0. Our main result shows that a compact two-sided hypersurface Σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma $$\end{document} immersed in C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {C}$$\end{document} with free boundary in ∂C\{0}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\partial \mathcal {C}\setminus \{0\}$$\end{document} and minimizing, up to second order, an anisotropic area functional under a volume constraint is contained in a Wulff-shape. The technique of proof also works for a non-smooth convex cone C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {C}$$\end{document} provided the boundary of Σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma $$\end{document} is away from the singular set of ∂C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\partial \mathcal {C}$$\end{document}.


Introduction
In this paper we study a variational problem for Euclidean hypersurfaces associated to an energy functional of anisotropic character.This means that the energy is computed by integrating an elliptic parametric function (or surface tension) which depends on the normal direction along the hypersurface.As it is explained in the introduction of Taylor [35] these functionals provide a mathematical model to study solid crystals.In the particular case of a constant surface tension we obtain the isotropic case, where the energy is proportional to the Euclidean area.
The surface tension that we consider in this work is the asymmetric norm given by the support function h K of a smooth strictly convex body K ⊂ R n+1 containing the origin in its interior.The corresponding anisotropic area of a two-sided hypersurface Σ with unit normal N has the expression A K (Σ) := Σ h K (N ) dΣ, where dΣ is the Euclidean area element.By using the metric projection onto K it is possible to define also anisotropic counterparts to some of the classical notions in the extrinsic geometry of hypersurfaces, see Reilly [29] and Section 2 for a precise description.When K is the round unit ball about the origin we have h K (N ) = 1, so that we recover the isotropic situation.We point out that the anisotropic geometry of hypersurfaces is typically introduced by means of a function F : S n → R + instead of a convex body K, see Remark 2.3 for more details.
It is well known that the round spheres uniquely solve the isoperimetric problem in R n+1 .In our anisotropic setting, the strictly convex hypersurface ∂K uniquely minimizes, up to translations and dilations centered at the origin, the anisotropic area computed with respect to the outer unit normal among hypersurfaces enclosing the same Euclidean volume.There are different proofs of this statement relying on analytical and geometric techniques, see Taylor [34,35], Fonseca and Müller [10], Brothers and Morgan [3], Milman and Rotem [25], and Cabré, Ros-Oton and Serra [4].The minimizer ∂K is usually referred to as the Wulff shape in honor of the crystallographer Georg Wulff, who first constructed the optimal crystal for a specified integrand [40].
As in the isotropic case, it is interesting to analyze the critical points (stationary hypersurfaces) and the second order minima (stable hypersurfaces) of the anisotropic area under a volume constraint.From the first variational formulas, a hypersurface Σ is anisotropic stationary if and only if its anisotropic mean curvature H K defined in (2.6) is constant, see for instance Palmer [28] or Clarenz [6].This motivated the generalization to the anisotropic context of some classical results for constant mean curvature hypersurfaces in R n+1 .In this direction, an Alexandrov type theorem (resp.a Hopf type theorem) characterizing the Wulff shape among compact embedded hypersurfaces (resp.immersed spheres) with H K constant was obtained by Morgan [26] for n = 1, and by He, Li, Ma and Ge [12] in higher dimension (resp.by He and Li [14], and Koiso and Palmer [22]).On the other hand, Palmer [28] and Winklmann [39] employed the second variation of the anisotropic area to show that, up to translation and homothety, the Wulff shape ∂K is the unique compact, two-sided, anisotropic stable hypersurface immersed in R n+1 .This is an anisotropic extension of a celebrated theorem of Barbosa and do Carmo [2] also discussed by Wente [38].
When we consider a smooth proper domain Ω instead of the whole space R n+1 , we are naturally led to study the anisotropic partitioning problem.Here, we seek minimizers of the functional A K among hypersurfaces inside Ω, possibly with non-empty boundary contained in ∂Ω, and separating a fixed Euclidean volume.The critical points for this problem are called anisotropic stationary hypersurfaces with free boundary in ∂Ω.From the calculus of A ′ K (0), and reasoning as Koiso and Palmer [19], it follows that a hypersurface Σ is anisotropic stationary with free boundary in ∂Ω if and only if H K is constant on Σ and the anisotropic normal N K introduced in (2.2) is tangent to ∂Ω along ∂Σ.The more general case of anisotropic capillary hypersurfaces arises when both, H K and the angle between N K and the inner unit normal to ∂Ω, are constant.These are critical points of an energy functional which involves not only A K (Σ) but also the wetting area of the set bounded by ∂Σ in ∂Ω.We remark that isotropic capillary hypersurfaces in convex domains of R n+1 have been extensively investigated, see for instance Ros and Souam [31], Wang and Xia [36], and the references therein.
An anisotropic stationary hypersurface with free boundary in ∂Ω is stable if A ′′ K (0) is nonnegative under variations preserving the volume separated by Σ and the boundary of Ω.The computation of A ′′ K (0) and the subsequent analysis of anisotropic stable hypersurfaces have been treated in several previous works.In [19,20,21], Koiso and Palmer considered stable capillary surfaces when Ω is a slab of R 3 and the Wulff shape is rotationally symmetric about a line orthogonal to ∂Ω.In [16,18], Koiso classified compact anisotropic stable capillary hypersurfaces disjoint from the edges in wedge-shaped domains of R n+1 .Inside a smooth domain of revolution Ω ⊂ R 3 , and for certain rotationally symmetric surface tensions, Barbosa and Silva [1] established that the totally geodesic disks orthogonal to the revolution axis are the unique compact stationary surfaces with free boundary in ∂Ω and meeting ∂Ω orthogonally.By assuming convexity of Ω they also discussed the stability of these surfaces.More recently, Jia, Wang, Xia and Zhang [15] have employed a Heintze-Karcher inequality and a Minkowski formula to prove that any compact, embedded, anisotropic capillary hypersurface inside an open halfspace of R n+1 is part of a Wulff shape.In [11], besides the computation of a very general second variation formula, Guo and Xia have generalized the argument of Winklmann [39] after Barbosa and do Carmo [2] to show that a compact, two-sided, anisotropic stable capillary hypersurface immersed in a Euclidean half-space of R n+1 is a truncated Wulff shape.
In this paper we are interested in the stability question when the ambient domain is a solid cone C ⊂ R n+1 with vertex at the origin and boundary ∂C smooth away from 0. In the isotropic case it is known that round balls about the vertex provide the unique solutions to the partitioning problem in a (possibly non-smooth) convex cone C, up to translations leaving C invariant, see Lions and Pacella [23], and Figalli and Indrei [9].When ∂C \ {0} is smooth, Ritoré and the author [30] solved the problem by combining an existence result with the classification of compact stable hypersurfaces with free boundary in ∂C.We must also mention the papers of Choe and Park [5], and Pacella and Tralli [27], where an Alexandrov type theorem for compact, embedded, stationary hypersurfaces in convex cones is proven.
In the anisotropic case it is easy to check that the truncated Wulff shape ∂K ∩ C is an anisotropic stationary hypersurface with free boundary in ∂C, see Example 4.2.In [37], Weng found some conditions ensuring that, up to translation and homothety, this is the unique compact and embedded hypersurface with constant anisotropic mean curvature in C. When the cone C is convex, results of Milman and Rotem [25], and Cabré, Ros-Oton and Serra [4] for anisotropic area functionals with homogeneous weights entail that ∂K ∩ C minimizes the anisotropic area (computed with respect to the outer unit normal) for fixed volume.We remark that the techniques in [25] and [4] still hold under weaker regularity assumptions on the cone C and the surface tension.By extending arguments of Figalli and Indrei [9], the uniqueness of ∂K ∩ C as a solution to the partitioning problem in C was recently discussed by Dipierro,Poggesi and Valdinoci [8,Thm. 4.2].The fact that ∂K ∩ C is a minimizer implies, in particular, the weaker property that it is an anisotropic stable hypersurface with free boundary in ∂C.Our main result, Theorem 4.7, is the following uniqueness statement: Up to a translation and a dilation centered at the vertex, any compact, connected, two-sided, anisotropic stable hypersurface immersed in a convex solid cone C with free boundary in ∂C \ {0} is part of the Wulff shape.
We must observe that, in general, we cannot expect ∂K ∩ C to be the only anisotropic stable hypersurface, up to translation and homothety, see Example 4.6.This motivates us to find additional conditions on the cone in order to deduce stronger rigidity consequences.With this idea in mind, in Corollary 4.8 we prove that a compact anisotropic stable hypersurface in a cone C over a smooth strictly convex domain of the sphere S n must be a dilation of ∂K ∩ C. On the other hand, for the case of a Euclidean half-space H, recently treated by Guo and Xia [11], we can conclude that ∂K ∩ H is the only compact anisotropic stable hypersurface, up to translations along the boundary hyperplane.
Our proof of Theorem 4.7 is different from the ones in [30] and [11], which are mainly inspired in Barbosa and do Carmo [2,Thm. (1.3)].Instead, we adapt to our setting the idea employed by Wente [38] to characterize the round spheres as the only compact, two-sided, stable hypersurfaces immersed in R n+1 .Wente applied the stability inequality with the special volume-preserving deformation obtained by parallel hypersurfaces dilated to keep the volume constant.After computing the second variation formulas for this particular variation, he observed that the area decreases strictly unless the hypersurface is a round sphere.By using dilations of anisotropic parallel hypersurfaces, Palmer [28] proved the uniqueness of the Wulff shape as a compact, two-sided, anisotropic stable hypersurface immersed in R n+1 .The variation of a hypersurface Σ by anisotropic parallels is defined as ψ t (p) := p + t N K (p), where p ∈ Σ, t ∈ R and N K is the anisotropic normal introduced in (2.2).With a similar argument, Koiso [16,Thm. 4], [18,Thm. 1] was able to classify compact anisotropic stable hypersurfaces disjoint from the edges in a wedge-shape domain Ω of R n+1 .Moreover, after a suitable translation, the same variation allows to treat the capillary case.We remark that such a deformation is possible because ∂Ω consists of finitely many pieces of hyperplanes.
In a convex cone C different from a half-space, Palmer's deformation does not work in general.Indeed, though N K is tangent to ∂C along ∂Σ, we cannot ensure that ∂C contains a small anisotropic normal segment p + t N K (p) centered at any point p ∈ ∂Σ.To solve this difficulty we replace the anisotropic parallels deformation of Σ with a variation Σ t := ψ t (Σ) associated to the one-parameter group {ψ t } t∈R of a vector field X on R n+1 such that X(0) = 0, X is tangent on ∂C \ {0} and X |Σ = N K .After this, as ∂C is invariant under dilations centered at 0, we can apply to Σ t a suitable dilation in order to construct a deformation that preserves both, ∂C and the volume of Σ.To finish the proof we compute the second derivative of the anisotropic area for this variation, and we discover that it is strictly negative unless Σ is contained in the Wulff shape (up to translation and homothety).
We emphasize that, to prove Theorem 4.7, we only need to calculate A ′′ K (0) for deformations associated to a smooth vector field X on R n+1 which is tangent on ∂C \ {0} and satisfies X |Σ = N K .This is done in Proposition 3.3 for arbitrary compact hypersurfaces with non-empty boundary in R n+1 .As a difference with respect to the anisotropic parallel deformation, a boundary term involving the acceleration vector field Z appears.For anisotropic stationary hypersurfaces in C we see in Proposition 4.4 that this term is related to the extrinsic Euclidean geometry of ∂C.As the formula for A ′′ K (0) in Proposition 4.4 is valid for any smooth domain Ω, this shows that the stability condition is more restrictive when Ω is convex.It is worth mentioning that, in most of the aforementioned works about anisotropic capillary hypersurfaces in a Euclidean domain Ω, this boundary term had no relevance since ∂Ω is totally geodesic.
Finally, we remark that the proof of Theorem 4.7 is still valid in non-smooth convex cones, provided the boundary ∂Σ of an anisotropic stable hypersurface Σ is inside a smooth open part of ∂C.
The paper is organized into three sections besides this introduction.In Section 2 we review some basic facts about the anisotropic geometry of hypersurfaces.In Section 3 we compute the second derivative of the anisotropic area for certain deformations of a compact hypersurface with non-empty boundary.Finally, in Section 4 we prove our classification result of anisotropic stable hypersurfaces in convex solid cones.

Preliminaries
In this section we gather some definitions and results that will be needed throughout this work.Starting with a Euclidean convex body we will use its support function and metric projection to introduce the anisotropic area and notions of anisotropic extrinsic geometry for two-sided hypersurfaces.The relation between this point of view and previous approaches is shown in Remark 2.3.
By a convex body (about the origin) we mean a compact convex set K ⊂ R n+1 containing 0 in its interior.The associated support function h K : R n+1 → R is given by where • , • is the standard scalar product in R n+1 .This defines an asymmetric norm in R n+1 (which is a norm when K is centrally symmetric about 0), see [32,Sect. 1.7.1].For any w = 0, the supporting hyperplane of K with exterior normal w is the set Π w := {p ∈ R n+1 ; p, w = h K (w)}.The corresponding support set is Π w ∩ K.This set is not empty because K is compact.If Π w ∩ K is a single point p, then h K is differentiable at w, and its gradient satisfies (∇h K )(w) = p, see [32,Cor. 1.7.3].
Henceforth, we suppose that K is a strictly convex body with smooth boundary ∂K.Thus, for any w = 0, there is a point π K (w) ∈ ∂K for which Π w ∩ K = {π K (w)}.This provides a map π K : R n+1 * → ∂K, which is called the K-projection, and verifies these identities On the other hand, if η K is the outer unit normal on ∂K, then (η K •π K )(w) = w for any w ∈ S n .This equality and the fact that π K (λ w) = π K (w) for any w = 0 and λ > 0 entail that π K (w) = η −1 K (w/|w|) for any w = 0 (here | • | stands for the Euclidean norm).As η K : ∂K → S n is a diffeomorphism, the support function h K is smooth on R n+1 * .Observe that, when K is the round unit ball about 0, then h K (w) = |w| and π K (w) = w/|w| for any w = 0, whereas η K (p) = p for any p ∈ ∂K = S n .Let ϕ 0 : Σ → R n+1 be a smooth two-sided immersed hypersurface, possibly with smooth boundary ∂Σ.Most of the time we will omit the map ϕ 0 .We will also identify any set S ⊆ Σ with ϕ 0 (S), and the tangent space T p Σ at a point p ∈ Σ with (dϕ 0 ) p (T p Σ).For a fixed smooth unit normal vector field N on Σ, the associated shape operator at p is the endomorphism B p : T p Σ → T p Σ introduced by B p (w) := −D w N .Here D denotes the Levi-Civita connection for the Euclidean metric.It is well known that B p is self-adjoint with respect to the metric induced by the scalar product in R n+1 .
The anisotropic Gauss map or anisotropic normal in Σ is the map N K : Σ → ∂K given by (2.2) By setting ϕ K := h K (N ) we infer from (2.1) that When ∂Σ = ∅ we introduce the anisotropic conormal by equality (2.4) where ν is the inner conormal along ∂Σ.Note that ν K is a normal vector to ∂Σ with N K , ν K = 0.
For any p ∈ Σ we have T NK (p) (∂K) = T p Σ, so that the differential (dN K ) p is an endomorphism of T p Σ.The anisotropic shape operator at p is defined by The anisotropic mean curvature at p is the number where div Σ is the divergence relative to Σ and tr(f ) is the trace of an endomorphism f : It is known, see for instance Palmer [28, p. 3666], that and equality holds if and only p is an anisotropic umbilical point, i.e., (B K ) p is a multiple of the identity map in T p Σ.This fact has a short proof that we reproduce here for the sake of completeness. Since and K is a strictly convex body, then (dπ K ) N (p) can be represented as a positive definite symmetric matrix of order n.Thus, there exist an orthonormal basis {e 1 , . . ., e n } of T p Σ and {λ 1 , . . ., λ n } ⊂ R + such that (dπ K ) N (p) (e j ) = λ j e j for any j = 1, . . ., n. From (2.5) we get (B K ) p (e i ) = n j=1 σ ij λ j e j , where σ ij := B p (e i ), e j .By using the Cauchy-Schwarz inequality in R n with the vectors (σ 11 λ 1 , . . ., σ nn λ n ) and (1, . . ., 1), we obtain Equality holds if and only if σ ij = 0 for any i = j and σ ii λ i = σ jj λ j for any i, j = 1, . . ., n.By denoting α := σ ii λ i , this is equivalent to that (B K ) p (w) = α w for any w ∈ T p Σ, as desired.
Examples 2.1.(i).If K is the round unit ball about 0, then N K = N , B K = B and H K equals the Euclidean mean curvature of Σ.
(ii).For a hyperplane Σ ⊂ R n+1 oriented with a unit normal N , the map N K is constant, so that B K = 0 and H K = 0 in Σ.
(iii).Consider the hypersurface Σ = ∂K with outer unit normal N = η K .Then N K (p) = p, (B K ) p (w) = −w and H K (p) = −1, for any p ∈ Σ and w ∈ T p Σ. Observe that, if K is not centrally symmetric about 0, and we take the inner normal N = −η K on Σ, then similar equalities need not hold and H K may be nonconstant.So, the role of the chosen unit normal N is very important for the computations.
The previous examples show that ∂K and Euclidean hyperplanes are anisotropic umbilical hypersurfaces, i.e., all of its points are anisotropic umbilical.Next, we provide a converse statement after Palmer [28, p. 3666] and Clarenz [7, p. 358].
Proposition 2.2.Let Σ be a two-sided, connected, anisotropic umbilical hypersurface immersed in R n+1 , and such that H K is constant.If H K = 0 then Σ is contained inside a hyperplane.If H K = 0 then, up to a translation and a dilation centered at 0, we have Σ ⊂ ∂K.
Proof.For any p ∈ Σ, there is α(p) ∈ R such that (B K ) p (w) = α(p) w for any w ∈ T p Σ.By definition (2.6) it is clear that α(p) = H K for any p ∈ Σ.By equation (2.5) we deduce that the differential of N K + H K Id vanishes on Σ.In case H K = 0 this implies that N K is a constant vector n, and so Σ is within a hyperplane with unit normal η K (n).When From here we infer that Σ ⊆ (c/H K ) − (1/H K ) (∂K), as we claimed.
We finish this section with the notions of anisotropic area and algebraic volume for a two-sided hypersurface Σ with unit normal N .If dΣ denotes the area element of Σ, then A K (Σ) is given by This coincides with the Euclidean area of Σ when K is the round unit ball about 0. When K is not centrally symmetric about 0, the value of A K (Σ) may depend on the normal N over Σ.On the other hand, for a compact hypersurface Σ, we follow Barbosa and do Carmo [2, Eq. (2.2)] to define For Σ embedded and small enough, the divergence theorem ensures that |V (Σ)| equals the Lebesgue measure of the cone over Σ with vertex at 0. If Σ is immersed then it is possible that V (Σ) = 0.
Consider a dilation δ λ (p) := λ p with λ > 0 and p ∈ R n+1 .A unit normal N λ on the hypersurface δ λ (Σ) is determined by N λ (δ λ (p)) := N (p) for any p ∈ Σ.Therefore, the change of variables formula and the fact that the Jacobian of the diffeomorphism δ λ|Σ : Σ → δ λ (Σ) equals λ n , entail that (2.9)These identities will play a relevant role in the proof of our main result in Section 4.
Remark 2.3.It is usual to introduce the anisotropic area A F associated to an arbitrary function F : S n → R + .This has the expression A F (Σ) := Σ F (N ) dΣ. Convexity assumptions on F are necessary to deduce fine properties for A F and for the F -Laplacian, see for instance Maggi [24,Ch. 20] or Palmer [28].By extending F as a 1-homogeneous function we get an asymmetric norm Φ in R n+1 .The Wulff shape for Φ is the convex body W Φ supported by Φ, see [24,Eq. (20.8)].As we remembered in the Introduction, W Φ minimizes the anisotropic perimeter among sets of the same Euclidean volume.In our context we have F = h K|S n , Φ = h K and W Φ = K.So, the initial convex body K is the optimal shape for the anisotropic area A K defined from its support function h K .In Palmer [28,Sect. 1] the Wulff hypersurface of F is the strictly convex hypersurface defined as φ(S n ), where φ(w) := F (w) w + (∇ S n F )(w).When F = h K|S n it follows from equation (2.1) that φ = π K on S n , so that the corresponding Wulff hypersurface φ(S n ) equals ∂K.

A second variation formula for the anisotropic area
In this section we compute the second derivative of the anisotropic area for certain deformations of an immersed hypersurface with non-empty boundary.The resulting formula will be employed in the proof of our main result in Theorem 4.7.We begin with some preliminary definitions.
A flow (of diffeomorphisms) in R n+1 is a smooth map φ : R n+1 × R → R n+1 such that φ(p, 0) = p for any p ∈ R n+1 , and the map φ t : R n+1 → R n+1 defined by φ t (p) := φ(p, t) is a diffeomorphism for any t ∈ R. Usually we will denote a flow by {φ t } t∈R .The associated velocity vector field is given by When Y is a smooth complete vector field on R n+1 , the corresponding one-parameter group of diffeomorphisms {φ t } t∈R is a flow with velocity Y .
For a hypersurface Σ immersed in R n+1 the variation of Σ induced by a flow {φ t } t∈R is the family {Σ t } t∈R , where Σ t := φ t (Σ) for any t ∈ R. The flow is compactly supported on Σ if there is a compact set C ⊆ Σ such that φ t (p) = p for any p ∈ Σ \ C and t ∈ R. If Σ is a two-sided hypersurface with unit normal N then, along the variation {Σ t } t∈R , we can find a smooth vector field N whose restriction to Σ t provides a unit normal N t with N 0 = N .For a fixed smooth strictly convex body K ⊂ R n+1 with support function h K , the anisotropic area functional for the variation {Σ t } t∈R is the function where Jac φ t is the Jacobian of the diffeomorphism φ t|Σ : Σ → Σ t .
Our objective is to provide a useful expression of A ′′ K (0) for some variations.For that we must compute A ′ K (0) first.We need some preliminary calculations that we gather below.Lemma 3.1.Let K ⊂ R n+1 be a smooth strictly convex body with support function h K .For a twosided hypersurface Σ immersed in R n+1 with unit normal N , the function ϕ K := h K (N ) satisfies , where ∇ Σ is the gradient relative to Σ and N ⊤ K is the tangent projection of the anisotropic normal N K .Moreover, for any smooth vector field U with compact support on Σ and normal component u := U, N , we have where H K is the anisotropic mean curvature and ν K is the anisotropic conormal.
Proof.Take a point p ∈ Σ and a vector w ∈ T p Σ.By using (2.1), (2.2) and that the shape operator B p is a self-adjoint endomorphism of T p Σ, we obtain This implies (3.2).Let us prove (3.3).For any w ∈ R n+1 we denote by w ⊤ and w ⊥ the projections of w with respect to T p Σ and (T p Σ) ⊥ , respectively.For a smooth vector field U on Σ, note that where in the second equality we have employed (2.3).By taking divergences relative to Σ and having in mind (2.6), we get The desired formula follows by using (2.4) and the divergence theorem on Σ.
Next, we compute the first variation of A K .This was previously derived by many authors, see for instance Clarenz [6, Sect.1], or Koiso and Palmer [19,Proof of Prop. 3.1].As in Koiso [17, Lem.9.1], our formula holds for arbitrary deformations of a Euclidean hypersurface with non-empty boundary.We include a short proof that will be helpful in the subsequent calculus of A ′′ K (0).Proposition 3.2.Let K ⊂ R n+1 be a smooth strictly convex body, Σ a two-sided immersed hypersurface with boundary, and {φ t } t∈R a flow in R n+1 with compact support on Σ.Then, we have where H K is the anisotropic mean curvature, u := X, N is the normal component of the velocity vector field X, and ν K is the anisotropic conormal along ∂Σ.
By differentiating under the integral sign in (3.1) and taking into account that φ 0 = Id, we obtain On the one hand it is well known, see Simon [33, §9], that On the other hand, from equations (2.1) and (2.2), it follows that where X ⊤ (p) is the projection of X(p) onto T p Σ.By substituting this information into the previous expression for h ′ p (0) and having in mind the third equality in (3.4), we arrive at (3.9) . Thus, the integrand in (3.6) is the evaluation at p of the function The proof finishes by applying the formula in (3.3) with U = X.
Second variation formulas for the anisotropic area under different hypotheses on Σ and the deformation can be found in Koiso and Palmer [19,Prop. 3.3], [20, Prop.3.3], Barbosa and Silva [1, Prop.3], and Guo and Xia [11,Prop. 3.5].In this work we only need to compute A ′′ K (0) when we move a hypersurface Σ with non-empty boundary by means of some special flows.Proposition 3.3.Let K ⊂ R n+1 be a smooth strictly convex body, Σ a compact two-sided immersed hypersurface with boundary, and X a smooth complete vector field on R n+1 such that X |Σ = N K .Then, for the one-parameter group of diffeomorphisms {φ t } t∈R associated to X, we have where H K is the anisotropic mean curvature, B K is the anisotropic shape operator, ν K is the anisotropic conormal and v := Z, N is the normal component of the vector field Z := D X X.
Proof.For any w ∈ R n+1 the notations w ⊤ and w ⊥ will stand for the projections of w onto T Σ and (T Σ) ⊥ , respectively.For a given point p ∈ Σ we define h p (t) and j p (t) as in (3.5).By differentiating under the integral sign twice in (3.1) and taking into account that φ 0 = Id, we get (3.10) Let us compute all the derivatives in the integrand above.
For the calculus of j ′ p (0) and j ′′ p (0) we refer the reader to Simon [33, §9].From (3.7), (2.6), and the fact that X |Σ = N K , we obtain On the other hand, we have D ei X, e j D ej X, e i , where {e 1 , . . ., e n } is an orthonormal basis of T p Σ. Equation (2.5) and the fact that so that (D ei X) ⊥ = 0 for any i = 1, . . ., n.It is also clear that n i,j=1 D ei X, e j D ej X, e i = tr (B K ) 2 p .
All this together with (3.11) shows that (3.12) Next, we compute h ′ p (0) and h ′′ p (0).Note that X, N = N K , N = ϕ K by (2.3).Hence, equation (3.9) implies that In order to calculate h ′′ p (0) we need an expression for h ′ p (t).From the definition in (3.5) we deduce

1). This entails that
and so , where we have employed the third equality in equation (3.4).Now, by having in mind (3.15), (3.13) and (3.12), we conclude that the integrand in (3.10) is the evaluation at p of the function From here the the proof finishes by applying the formula in (3.3) with U = Z.Lemma 3.4.In the conditions of Proposition 3.3, for any p ∈ Σ, we have Proof.Take an orthonormal basis {e 1 , . . ., e n } of T p Σ. We use the flow {φ t } t∈R associated to X to construct, for any i = 1, . . ., n, a smooth vector field E i around p which is tangent on any Σ t = φ t (Σ) while satisfying E i (p) = e i and [X, We will compute the different terms in the previous equation.
Since |N | 2 = 1 we get D X N , N = 0.By differentiating and applying equation (3.14), we deduce On the other hand, we differentiate twice with respect to X in equality N , E i = 0.By taking into account that D X(p) N = 0 and [X, E i ] = 0, we obtain Since the Riemann curvature tensor vanishes for the standard metric in R n+1 , we infer By substituting (3.17) and (3.18) into (3.16) the proof follows.
Remarks 3.5.(i).If K ⊂ R n+1 is the round unit ball about 0, then the anisotropic area coincides with the Euclidean area and the obtained formulas are well known.
(ii).It is interesting to observe that, unless K is centrally symmetric about 0, the formulas for A ′ K (0) and A ′′ K (0) may depend on the unit normal vector N fixed on Σ. (iii).In Proposition 4.4 we will see that, under some extra conditions, the boundary integrand in the expression of A ′′ K (0) has a geometric interpretation.

Anisotropic stable hypersurfaces in solid cones
In this section we consider a Euclidean solid cone and study compact hypersurfaces immersed in the cone and minimizing the anisotropic area up to second order for deformations preserving the volume of the hypersurface and the boundary of the cone.In this situation we will show that the variational formulas computed in Section 3 can be slightly simplified.After that we will prove our main theorem, where we characterize these second order minima when the cone is convex.It is worth mentioning that, excluding this classification statement, all the results in this section still hold when we replace the cone with any smooth Euclidean open set.
We begin by introducing some notation and definitions.For a domain D ⊂ S n with smooth boundary, the solid cone over D is the set This is a domain of R n+1 with boundary ∂C smooth away from 0. We call ξ to the inner unit normal along ∂C \ {0}.Note that C and ∂C are invariant under the dilations δ λ centered at 0. It is also clear that C coincides with an open half-space when D is an open hemisphere.
Let Σ be a smooth, compact, two-sided hypersurface immersed in C with smooth boundary ∂Σ in ∂C \ {0}.We suppose that Σ ∩ ∂C = ∂Σ, so that Σ \ ∂Σ ⊆ C. We fix a smooth unit normal vector field N on Σ.The inner conormal vector of ∂Σ in Σ is represented by ν.In the planar distribution T (∂Σ) ⊥ we choose the orientation induced by {ν, N }.Thus, for any p ∈ ∂Σ, there is a unique µ(p) ∈ T p (∂Σ) ⊥ such that {ξ(p), µ(p)} is a positively oriented orthonormal basis.Observe that µ is tangent to ∂C and normal to ∂Σ.It is easy to check that these equalities hold along where θ is the oriented angle function between ν and ξ in T (∂Σ) ⊥ .As a consequence, for a fixed smooth strictly convex body K ⊂ R n+1 , the anisotropic normal N K on Σ and the anisotropic conormal ν K along ∂Σ given in (2.2) and (2.4) verify (4.2) A flow of diffeomorphisms {φ t } t∈R in R n+1 is admissible for C if φ t (∂C) = ∂C and φ t (0) = 0, for any t ∈ R. The induced variation {Σ t } t∈R of Σ satisfies ∂Σ t ⊂ ∂C \ {0} and Σ t ∩ ∂C = ∂Σ t , for any t ∈ R.Moreover, the velocity vector field X is tangent on ∂C \ {0} and vanishes at 0. We consider the anisotropic area functional A K (t) introduced in (3.1) and the volume functional V (t), which assigns to any t ∈ R the algebraic volume enclosed by Σ t as defined in (2.8).We computed A ′ K (0) in Proposition 3.2.On the other hand, it is well known ([2, Eq. ( 2. 3)]) that where u := X, N .We say that the flow preserves the volume of Σ if V (t) is constant for any t small enough.This implies that Σ u dΣ = 0. Conversely, for any smooth function u : Σ → R with Σ u dΣ = 0, there is a flow {φ t } t∈R admissible for C, preserving the volume of Σ, and such that X, N = u on Σ, see Barbosa and do Carmo [2, Lem.(2.4)].
A hypersurface Σ in the previous conditions is anisotropic stationary if A ′ K (0) = 0 for any admissible flow for C preserving the volume of Σ.The first variational formulas in Proposition 3.2 and equation (4.3), together with the aforementioned construction of volume-preserving flows, lead to a characterization of anisotropic stationary hypersurfaces.Indeed, we can reason as Koiso and Palmer [19,Prop. 3.1] to deduce the next result, which generalizes the isotropic situation.Example 4.2.Take the hypersurface Σ = ∂K ∩ C with unit normal N = η K (the one pointing outside K).Note that ∂K meets ∂C transversally because K is convex and contains 0 in its interior.Hence, Σ is a compact hypersurface with boundary ∂Σ such that Σ ∩ ∂C = ∂Σ.In Examples 2.1 (iii) we saw that H K (p) = −1 and N K (p) = p, for any p ∈ Σ.As ∂C is invariant under dilations centered at 0 then N K (p) ∈ T p (∂C) for any p ∈ ∂Σ and so, N K , ξ = 0. From the previous proposition we conclude that Σ is anisotropic stationary.
In the isotropic case the orthogonality condition for a stationary hypersurface Σ entails that ν = ξ along ∂Σ.In the next lemma we establish a similar equality for the anisotropic conormal ν K that will be useful in future results.Lemma 4.3.Let Σ be a two-sided hypersurface immersed in C with boundary ∂Σ in ∂C \ {0}.If N K , ξ = 0 in ∂Σ, then cos θ never vanishes and where ϕ K := h K (N ) and θ is the oriented angle function between ν and ξ in T (∂Σ) ⊥ .
Let us compute ν K , ξ .Consider the projection N * K of N K onto T (∂Σ) ⊥ .By using (4.2), equality N K , ξ = 0, and that {ν, N } and {ξ, µ} are orthonormal basis of T (∂Σ) ⊥ , we obtain which is a positive number.Hence, ν K , ξ never vanishes along ∂Σ.On the other hand, by substituting the expression for ξ in (4.1) into equality N K , ξ = 0, we get The definition of ν K and the two previous relations lead to and so ν K , ξ ϕ K cos θ = 1 along ∂Σ.
This implies that cos θ never vanishes and allows to deduce the announced expression for ν K .Now, we can provide new expressions for the derivatives of the anisotropic area when we consider an anisotropic stationary hypersurface Σ in C. On the one hand, the boundary term in the formula for A ′ K (0) obtained in Proposition 3.2 vanishes for any flow {φ t } t∈R admissible for C.This comes from Proposition 4.1 and Lemma 4.3 since the velocity vector field X is tangent to ∂C \ {0}.Hence (4.4) where u := X, N on Σ.Note that, for the flow φ t (p) := e t p, the first equation in (2.9) that A ′ K (0) = nA K (Σ).Thus, by having in mind (4.4) and equation (2.8), it follows that (4.5) This identity is a Minkowski-type formula for compact anisotropic stationary hypersurfaces in C.
On the other hand, the boundary integrand appearing in A ′′ K (0), see Proposition 3.3, can be written in terms of the Euclidean extrinsic geometry of ∂C.We prove this in the next proposition, where we also compute the second derivative of the volume for certain flows.Proposition 4.4.Let Σ be a two-sided anisotropic stationary hypersurface immersed in C with boundary ∂Σ in ∂C \ {0}.Consider a smooth complete vector field X on R n+1 such that X |Σ = N K , X(0) = 0 and X is tangent to ∂C \ {0}.Then, for the admissible flow {φ t } t∈R defined by the oneparameter group of diffeomorphisms associated to X, we have where B K is the anisotropic shape operator of Σ, II is the second fundamental form of ∂C \ {0} with respect to the inner normal, θ is the oriented angle function between ν and ξ in T (∂Σ) ⊥ , and v := Z, N is the normal component of Z := D X X. Remark 4.5.Observe that N K is tangent to ∂C along ∂Σ because Σ is anisotropic stationary.This guarantees the existence of a vector field X in the conditions of the statement and that the term II(N K , N K ) is well defined.Note also that cos θ never vanishes by Lemma 4.3.
Proof of Proposition 4.4.We first check that (4.6) holds.From Proposition 3.3 it suffices to see that By Lemma 4.3 we know that It is clear that X, ξ = 0 along ∂C because X is tangent to ∂C \ {0} and X(0) = 0.By differentiating with respect to X, we obtain 0 = D X X, ξ + X, D X ξ = Z, ξ − II(X, X), so that Z, ξ = II(X, X) = II(N K , N K ) along ∂Σ.
We now turn to the main result of the paper.This is a classification of anisotropic stable hypersurfaces in a Euclidean solid cone.A two-sided hypersurface Σ immersed in C with boundary ∂Σ in ∂C \ {0} is anisotropic stable if A ′ K (0) = 0 and A ′′ K (0) 0 for any flow admissible for C and preserving the volume of Σ.When the cone is convex, the hypersurface ∂K ∩ C is anisotropic stable since it minimizes the anisotropic area (computed with respect to the outer unit normal) among compact hypersurfaces in C separating the same volume, see [25,Cor. 1.2] and [4,Thm. 1.3].So, it is natural to ask if this property characterizes ∂K ∩ C up to translations and dilations centered at 0. The next example shows that, in general, the answer is negative.
Example 4.6.Let C be a convex cone different from a Euclidean half-space and such that ∂C contains a half-hyperplane P (for instance, any solid cone in R 2 satisfies this property).As we observed above, if P + is the open half-space with C ⊂ P + , then ∂K ∩ P + minimizes the anisotropic area in P + for fixed volume.In particular, by applying a suitable translation along P we produce an anisotropic stable hypersurface in C which has not the form p 0 + λ (∂K ∩ C) for some p 0 ∈ R n+1 and λ > 0.
This example illustrates that the optimal conclusion to be deduced for an anisotropic stable hypersurface Σ in C is that Σ ⊂ ∂K, up to translation and homothety.We prove this fact in the next uniqueness statement under our regularity conditions for Σ, C and K. Theorem 4.7.Let C ⊂ R n+1 be a solid convex cone over a smooth domain of S n .Consider a compact, connected, two-sided hypersurface Σ immersed in C with smooth boundary ∂Σ in ∂C \ {0} and such that Σ ∩ ∂C = ∂Σ.If Σ is anisotropic stable for the area A K defined by a smooth strictly convex body K ⊂ R n+1 then, there is p 0 ∈ R n+1 and λ > 0 such that Σ ⊂ p 0 + λ (∂K).
Proof.We will follow the idea explained in the Introduction.As Σ is anisotropic stationary, Proposition 4.1 implies that H K is constant on Σ and N K , ξ = 0 along ∂Σ.So, we can find a smooth complete vector field X on R n such that X |Σ = N K , X(0) = 0 and X is tangent to ∂C \ {0}.It is clear that u = ϕ K on Σ, where u := X, N and ϕ K := h K (N ).Let {ψ t } t∈R be the one-parameter group of diffeomorphisms associated to X.We take the functionals A K (t) := A K (Σ t ) and V (t) := V (Σ t ) associated to the variation Σ t := ψ t (Σ).From the Minkowski formula in equation (4.5) we have H K = 0 and V (Σ) = 0. Hence, there is ε > 0 such that V (t) has the same sign as V (Σ) for any t ∈ (−ε, ε).Next, we apply to Σ t a dilation δ λ(t) (p) := λ(t) p with λ(t) > 0, so that the volume of δ λ(t) (Σ t ) equals the volume of Σ.Since V (δ λ(t) (Σ t )) = λ(t) n+1 V (t), see (2.9), we get , for any t ∈ (−ε, ε).
In particular, λ(0) = 1.We extend λ(t) to a smooth positive function on R. If we define φ t := δ λ(t) • ψ t , for any t ∈ R, then we produce an admissible flow for C that preserves the volume of Σ.Hence, the anisotropic stability of Σ entails that the functional a K (t) := A K (φ t (Σ)) satisfies a ′′ K (0) 0. To prove the theorem we need to compute a ′′ K (0).By equation (2.9) we know that (4.9) Thus, the calculus of a ′′ K (0) relies on the values of λ ′ (0), λ ′′ (0), A ′ K (0) and A ′′ K (0).By using equations (4.3), (4.4) and the fact that u = ϕ K on Σ, we obtain On the other hand, from the expression of λ(t) in (4.8), we have and so, from equation (4.5), it follows that Now we compute A ′′ K (0) and λ ′′ (0).From (4.7) and (4.6), since H K is constant, we get where α := Σ v dΣ and v := D X X, N .Recall that II is the second fundamental form of ∂C \ {0} with respect to the inner normal ξ.On the other hand, from (4.12) we deduce .
When we employ the identity (4.5) in the three summands above, we arrive at Finally, we differentiate into equation (4.9) to infer . By substituting above the expressions in (4.13), (4.16), (4.11), (4.15), and simplifying, we obtain (4.17) The first integrand in the previous formula is nonnegative by (2.7).The convexity of C implies that II(N K , N K ) 0 along ∂Σ.By the first equation in (4.1) we have cos θ = ν, ξ , which is positive along ∂Σ because C is convex and ν points into Σ ⊂ C. Hence, the stability inequality a ′′ K (0) 0 entails that tr(B 2 K ) = nH 2 K on Σ (this means that Σ is anisotropic umbilical) and II(N K , N K ) = 0 along ∂Σ.Since H K = 0, the proof finishes by invoking Proposition 2.2.
The situation described in Example 4.6 leads us to seek additional conditions on the cone in order to deduce stronger uniqueness conclusions.In this direction we can prove the next statement.
Corollary 4.8.Let C ⊂ R n+1 be a solid cone over a smooth strictly convex domain of S n .Consider a compact, connected, two-sided hypersurface Σ immersed in C with smooth boundary ∂Σ in ∂C \ {0} and such that Σ ∩ ∂C = ∂Σ.If Σ is anisotropic stable for the area A K defined by a smooth strictly convex body K ⊂ R n+1 , then Σ = λ (∂K ∩ C) for some λ > 0.
The fact that the base set D ⊂ S n of the cone C is a strictly convex domain means that the second fundamental form of ∂D as a hypersurface of S n is always definite positive with respect to the inner unit normal.This implies that II p (w, w) > 0 for any p ∈ ∂Σ and any w ∈ T p (∂C) \ {0} orthogonal to p. From the identity II(N K , N K ) = 0 at the end of the proof of Theorem 4.7, it follows that (p− p 0 )/λ is proportional to p, for any p ∈ ∂Σ.From here we get p 0 = 0, and this completes the proof.
A special example of convex cone is a half-space H ⊂ R n+1 .Since ∂H is smooth we do not need to assume ∂Σ ⊂ ∂H \ {0}.In this case we can deduce that ∂K ∩ H is the unique compact anisotropic stable hypersurface in H, up to translations along ∂H.A similar result for anisotropic stable capillary hypersurfaces in H has been given by Guo and Xia [11,Thm. 1.1] with a different proof.
Proof.After applying Theorem 4.7 it remains to see that p 0 ∈ ∂H.By reasoning as in the proof of Corollary 4.8 we obtain N K (p) = (p − p 0 )/λ, for any p ∈ Σ.The orthogonality condition N K , ξ = 0 along ∂Σ entails that p − p 0 ∈ T p (∂H) for any p ∈ ∂Σ.As ∂H is a Euclidean hyperplane, the straight line starting from a point p ∈ ∂Σ with p = p 0 and generated by the vector p − p 0 is entirely contained in ∂H.This shows that p 0 ∈ ∂H, as we claimed.
Remark 4.10 (Non-smooth cones).The proof of Theorem 4.7 is still valid when C is an arbitrary open convex cone (the base domain D ⊂ S n need not be smooth), provided the boundary ∂Σ is contained in a smooth open portion of ∂C.For instance, the conclusion holds for a compact anisotropic stable hypersurface Σ disjoint from the edge of a domain C bounded by two transversal hyperplanes.In this direction, Koiso [16,Thm. 4], [18,Thm. 1] has characterized compact anisotropic stable capillary hypersurfaces in wedge-shaped domains of R n+1 .Remark 4.11 (Planar cones).For a solid cone C ⊂ R 2 , the boundary ∂C is the union of two-closed half-lines leaving from 0. Thus, we have II = 0 along ∂C \ {0} and the boundary term in (4.17) disappears.Hence, the conclusion of Theorem 4.7 remains valid even if the cone is not convex.

Proposition 4 . 1 .
A two-sided hypersurface Σ immersed in C with boundary ∂Σ in ∂C \ {0} is anisotropic stationary if and only if the anisotropic mean curvature H K is constant on Σ and N K , ξ = 0 along ∂Σ.