1 Introduction

Surface theory in Riemannian 3-manifolds admitting a Killing vector field has experienced an increasing interest during the last decades. This interest has been specially significant in the case of product spaces \(M\times \mathbb {R}\) as well as in homogeneous 3-manifolds with isometry group of dimension 4, also known as \(\mathbb {E}(\kappa ,\tau )\)-spaces. As a common framework for these spaces, the theory of unit Killing submersions was developed, being the work of Rosenberg, Souam, and Toubiana [23] pioneer in this topic (see also [6, 12, 21, 24]).

The first classification results for unit Killing submersions were obtained by the second author. It is proved in [14] that, given a simply connected surface M and a function \(\tau \in \mathcal {C}^\infty (M)\), there exists a Riemannian submersion \(\pi :\mathbb {E}\rightarrow M\) such that \(\mathbb {E}\) is simply connected and orientable, the fibers of \(\pi \) are the integral curves of a unit Killing vector field \(\xi \) and \(\pi \) has bundle curvature \(\tau \) (i.e., \(\overline{\nabla }_X\xi =\tau X\wedge \xi \) for all vector fields \(X\in \mathfrak {X}(\mathbb {E})\), where \(\overline{\nabla }\) stands for the Levi–Civita connection in \(\mathbb {E}\)). Moreover, \(\pi \) turns out to be unique under these assumptions. Dropping the unitary condition, a Riemannian submersion fulfilling the aforementioned properties will be called a Killing submersion (see Definition 2.1).

This generalization is motivated by the fact that some 3-manifolds admitting non-unitary Killing vector fields have recently gained certain relevance. Among them, we highlight the simply connected homogeneous ones (we refer the reader to the survey by Meeks and Pérez [17] and the references therein), whose spaces of Killing vector fields have dimension at least 3, giving rise to Killing submersion structures in any given Killing direction (see Examples 2.4 and 2.5). Although we will give a complete classification of Killing submersions, it is also important to mention that several authors have already dealt with hypersurfaces in generic n-manifolds admitting a Killing vector field (see, for instance, [2, 5]).

Given an arbitrary Killing submersion \(\pi :\mathbb {E}\rightarrow M\), we can associate to \(\pi \) two geometric functions. One of them is the bundle curvature \(\tau \), a function determining the horizontal part of the curvature 2-form associated with \(\pi \) (see Eq. (2.1)), which generalizes the definition of bundle curvature in the unitary case and measures the non-integrability of the horizontal distribution. The second function is the Killing length \(\mu \), i.e., the length of the Killing vector field. Both \(\tau \) and \(\mu \) are constant along fibers and induce smooth functions on M. This way we can pose the problem of existence and uniqueness of Killing submersions fixing a Riemannian surface M and the functions \(\tau ,\mu \in \mathcal {C}^\infty (M)\), \(\mu >0\).

Section 2 is devoted to give the following complete solution:

  • If M is simply connected, then there exists a Killing submersion \(\pi :\mathbb {E}\rightarrow M\) with bundle curvature \(\tau \) and Killing length \(\mu \), which is unique provided that \(\mathbb {E}\) is simply connected.

    • If M is topologically \(\mathbb {R}^2\), then \(\pi \) is a trivial fibration, and we are able to get an explicit model for \(\pi \) (see Theorem 2.6 and its proof). In particular, \(\mathbb {E}\) is diffeomorphic to \(\mathbb {R}^3\).

    • If M is topologically \(\mathbb {S}^2\), then \(\pi \) admits a global section if and only if \(\int _M\frac{\tau }{\mu }=0\), and in that case \(\mathbb {E}\) is diffeomorphic to \(\mathbb {S}^2\times \mathbb {R}\). Otherwise \(\pi \) is topologically the Hopf fibration and \(\mathbb {E}\) is diffeomorphic to \(\mathbb {S}^3\) (see Theorem 2.9).

  • If \(\pi :\mathbb {E}\rightarrow M\) is a Killing submersion and M and \(\mathbb {E}\) are not simply connected, then there exists a Killing submersion \(\widetilde{\pi }:\widetilde{\mathbb {E}}\rightarrow \widetilde{M}\), being \(\widetilde{M}\) and \(\widetilde{\mathbb {E}}\) the universal Riemannian coverings of M and \(\mathbb {E}\), respectively, and a group G of isometries on \(\widetilde{\mathbb {E}}\) preserving the Killing direction, such that G acts properly discontinuously on \(\widetilde{\mathbb {E}}\) and \(\mathbb {E}=\widetilde{\mathbb {E}}/G\) (see Theorem 2.13).

Existence is therefore guaranteed by the above results, but uniqueness may fail when the base surface is not simply connected (see Example 2.14).

It is also important to remark that Killing submersions are well-understood objects at the level of differential topology (see, for instance, [11, 25]), so our contribution is to give a uniform treatment to all of them in terms of the functions \(\tau \) and \(\mu \), which prove to be useful when analyzing the geometry and topology of surfaces. One example is the duality between mean curvature and bundle curvature shown by the Calabi-type correspondence [13, Theorem 3.7]; another example is the fact that \(\pi :\mathbb {E}\rightarrow M\) admits a global section if and only if either M is non-compact or M is compact and \(\int _M\frac{\tau }{\mu }=0\) (see Proposition 3.3).

In Sect. 3, we will introduce graphs in Killing submersions as sections of the submersion over open sets in the base, and give a divergence-type equation for their mean curvature (see Lemma 3.1). We will prove that any Killing submersion \(\pi :\mathbb {E}\rightarrow M\) over a compact surface M which admits a global section also admits a global minimal section (see Theorem 3.6). The proof is based on minimizing area inside the isotopy class of a global section by means of a result of Meeks, Simon and Yau [20]. This argument is not valid if M is a sphere, but in that case we can employ a Calabi-type duality together with an existence result of Gerhardt [10] for graphs with prescribed mean curvature in Lorentzian warped products. The aforesaid duality is inspired by the work of Albujer and Alías [1] (see also [13, Theorem 3.7]). Let us now give two ways of understanding Theorem 3.6:

  • On the one hand, Theorem 3.6 solves the Bernstein problem in any Killing submersion \(\pi :\mathbb {E}\rightarrow M\) provided that M is compact: Entire graphs with constant mean curvature \(H\in \mathbb {R}\) exist if and only if \(H=0\) and \(\int _M\frac{\tau }{\mu }=0\), and they are unique up to isometries.

  • On the other hand, Theorem 3.6 can be regarded as the solution to a Plateau problem with empty boundary, which from the point of view of calculus of variations is equivalent to minimizing the area functional

    $$\begin{aligned} \mathcal {A}(u)=\int _M\sqrt{\rho ^2+\Vert \nabla u+X\Vert ^2}, \end{aligned}$$
    (1.1)

    where M is a compact surface, \(\rho \in \mathcal {C}^\infty (M)\) is positive, and X is a smooth vector field in M. Theorem 3.6 guarantees the existence of \(u\in \mathcal {C}^\infty (M)\) minimizing (1.1) among all Lipschitz functions on M (note that (1.1) is the area of an entire graph in a Killing submersion over M with Killing length \(\mu =\rho ^{-1}\) and bundle curvature \(\tau =2\rho ^{-1}\,\mathrm {div}(JX)\), see Eq. (3.4)).

Using the work of Simon on the existence of embedded minimal spheres [26], we also deduce that a Killing submersion \(\pi :\mathbb {E}\rightarrow M\), M not necessarily compact, admits an immersed minimal sphere if and only if M is a sphere itself (see Theorem 3.9). These minimal spheres are unique up to ambient isometries if the submersion is trivial by Theorem 3.6, and also if \(\mathbb {E}\) is a homogeneous 3-sphere [16], but uniqueness is not hitherto known to hold in general.

In Sect. 4, we will address the problem of stability of a compact orientable surface with constant mean curvature immersed in the total space of a Killing submersion \(\pi :\mathbb {E}\rightarrow M\). Recall that a two-sided surface \(\Sigma \) immersed in \(\mathbb {E}\) (here two-sided is equivalent to orientable, since \(\mathbb {E}\) is assumed orientable) has constant mean curvature \(H\in \mathbb {R}\) if and only if it is a critical point of \(\mathcal J=\mathrm {Area}-2H\cdot \mathrm {Vol}\), as proved by Barbosa et al. [3]. Then \(\Sigma \) is said stable if

$$\begin{aligned} \mathcal J''(0)=-\int _\Sigma \left( \Delta f+(|A|^2+\mathrm {Ric}(N))f\right) f\ge 0, \end{aligned}$$

for all compactly supported smooth functions \(f\in \mathcal C_0^\infty (\Sigma )\). In other words, \(\Sigma \) is a second-order minimum of \(\mathcal J\) for all compactly supported normal variations of \(\Sigma \) via the trivialization of the normal bundle given by a unit normal vector field N. Here A denotes the shape operator. This is not the notion of stability associated with the isoperimetric problem, where the variations are required to be volume-preserving, but a stronger one (a comprehensive introduction to stability can be found in the survey by Meeks et al. [19]). We will use a characterization due to Fischer-Colbrie [9] stating that \(\Sigma \) is stable if and only if there exists a smooth function \(u>0\) on \(\Sigma \) satisfying \(Lu=0\), where \(L=\Delta +|A|^2+\mathrm {Ric}(N)\) is the so-called stability operator of \(\Sigma \) (a more precise expression for this Schrödinger operator in the Killing submersion setting will be given in Lemma 4.1).

Using the fact that the angle function \(\nu =\langle N,\xi \rangle \) associated with a vertical Killing vector field \(\xi \) lies in the kernel of L, we will show that any compact orientable stable surface immersed in \(\mathbb {E}\) with constant mean curvature must be either a vertical cylinder (i.e., everywhere vertical) or an entire minimal graph (see Theorem 4.2). In particular, if M is not compact and the fibers of \(\pi \) are not compact either, then \(\mathbb {E}\) does not admit compact orientable stable surfaces with constant mean curvature. In other words, the first eigenvalue associated with L is negative for such surfaces in \(\mathbb {E}\), which improves previous results by Meroño and Ortiz [21, Corollaries 4, 6, and 8] and generalizes other results by the second author, Pérez and Rodríguez [15] in \(\mathbb {E}(\kappa ,\tau )\)-spaces. Note also that Theorem 4.2 reveals the topology of a stable surface \(\Sigma \) immersed in \(\mathbb {E}\) with constant mean curvature \(H\in \mathbb {R}\); for instance, if \(\Sigma \) is not a torus, then \(H=0\) and \(\Sigma \) is an entire minimal graph.

Substituting compactness with the weaker hypothesis of completeness and finiteness of the Abresch-Rosenberg total curvature, Espinar and Trejos [7, Theorem 6] provide sharp estimates for the first eigenvalue of L in \(\mathbb {E}(\kappa ,\tau )\), but the problem of finding a sharp estimate in the class of compact surfaces with constant mean curvature H immersed in \(\mathbb {E}(\kappa ,\tau )\) remains open. We also obtain a general (but possibly not sharp) result proving nonexistence of complete stable surfaces with constant mean curvature (see Corollary 4.5) by adapting results of Meeks et al. [19] and Rosenberg [22] to the context of Killing submersions.

As a final remark, although our techniques only produce entire minimal graphs when the base surface is compact, we conjecture that any Killing submersion over a non-compact surface admits a global minimal section. This conjecture is also related to the existence of global spacelike sections with prescribed mean curvature in certain Lorentzian 3-manifolds (see [13]).

2 Killing submersions

Let \(\pi :\mathbb {E}\rightarrow M\) be a Riemannian submersion from a Riemannian 3-manifold \(\mathbb {E}\) to a Riemannian surface M, both of them connected and orientable. A vector \(v\in T\mathbb {E}\) will be called vertical when \(v\in \ker (\mathrm {d}\pi )\) and horizontal when \(v\in \ker (\mathrm {d}\pi )^\bot \). Recall that \(\pi \) is Riemannian if it preserves the length of horizontal vectors.

Definition 2.1

\(\pi :\mathbb {E}\rightarrow M\) is called a Killing submersion if the fibers of \(\pi \) are the integral curves of a complete Killing vector field \(\xi \in \mathfrak {X}(\mathbb {E})\) without zeroes.

It is important to notice that \(\xi \) is not unique under these conditions, since multiplying \(\xi \) by a nonzero real constant also gives a Killing vector field without zeroes generating the same integral curves.

Given a Killing submersion \(\pi :\mathbb {E}\rightarrow M\) and fixing some \(\xi \) satisfying Definition 2.1, we can consider the connection 1-form \(\alpha \in \Omega ^1(\mathbb {E})\), \(\alpha (X)=\langle X,\xi \rangle \), and the curvature 2-form \(\omega =\tfrac{1}{2}\mathrm {d}\alpha \) given by \(\omega (X,Y)=\langle \overline{\nabla }_X\xi ,Y\rangle \) for all \(X,Y\in \mathfrak {X}(\mathbb {E})\), being \(\overline{\nabla }\) the Levi-Civita connection in \(\mathbb {E}\). Since \(\xi \) is Killing, we get that \(\omega \) is skew-symmetric, and the function \(\tau \in \mathcal {C}^\infty (\mathbb {E})\) defined as

$$\begin{aligned} \tau (p)=\frac{-1}{\Vert \xi _p\Vert }\,\omega _p(e_1,e_2), \end{aligned}$$
(2.1)

where \(\{e_1,e_2,\xi _p/\Vert \xi _p\Vert \}\) is a positively oriented orthonormal basis of \(T_p\mathbb {E}\), does not depend on the choice of \(\{e_1,e_2\}\) or \(\xi \). Moreover, if \(\xi \) is unitary, then \(\tau \) satisfies the well-known identity \(\overline{\nabla }_X\xi =\tau X\wedge \xi \) for all \(X\in \mathfrak {X}(\mathbb {E})\), so \(\tau \) will be called the bundle curvature of the Killing submersion, extending previous definitions in the unitary case [6, 14, 23, 24]. The bundle curvature also measures the non-integrability of the horizontal distribution associated with \(\pi \), in the sense that \(\tau =0\) if and only if \(\ker (\mathrm {d}\pi )^\bot \) is integrable (see Eq. (2.5) below).

The 1-parameter group of isometries associated with \(\xi \) will be denoted by \(\{\phi _t\}\), and its elements will be called vertical translations. Note that \(\phi _t\) is defined for all values of \(t\in \mathbb {R}\) since \(\xi \) is assumed complete. Due to the fact that \(\phi _t:\mathbb {E}\rightarrow \mathbb {E}\) is an isometry such that \((\phi _t)_*\xi =\xi \) and \((\phi _t)_*\omega =\omega \), the bundle curvature is constant along the fibers of \(\pi \), and so is the Killing length \(\mu =\Vert \xi \Vert \in \mathcal C^\infty (\mathbb {E})\). It follows that both \(\tau \) and \(\mu \) induce functions in M that will be also denoted by \(\tau ,\mu \in \mathcal C^\infty (M)\).

Definition 2.2

Two Killing submersions \(\pi :\mathbb {E}\rightarrow M\) and \(\pi ':\mathbb {E}'\rightarrow M\), over the same base surface M, are isomorphic if there exists an isometry \(T:\mathbb {E}\rightarrow \mathbb {E}'\) such that \(\pi '\circ T=\pi \).

Given a base surface M and two functions \(\tau ,\mu \in \mathcal {C}^\infty (M)\), \(\mu >0\), our first goal is to classify (up to isomorphism) all Killing submersions over M with bundle curvature \(\tau \) and Killing length \(\mu \). This goal was achieved in the unitary case, provided that M is simply connected in [14]. Nonetheless, it is important to point out some differences in the non-unitary case:

  • The fibers of \(\pi \) are not geodesic in general. As shown in Eq. (4.1) below, if \(\gamma \) is a unit-speed parametrization of a fiber, then \(\overline{\nabla }_{\gamma '}\gamma '=-\frac{1}{\mu }\overline{\nabla }\mu \), so geodesic fibers correspond to the critical points of \(\mu \).

  • If a fiber of \(\pi \) has finite length, then all fibers of \(\pi \) have finite length, but this length may vary from fiber to fiber. In spite of that, there exists a minimal \(\ell >0\) such that \(\phi _\ell :\mathbb {E}\rightarrow \mathbb {E}\) is the identity map. This value of \(\ell \) is the analogue to the constant length of the fiber in the non-unitary case.

Note also that we can change the metric in the vertical direction preserving the submersion structure and making the Killing vector field unitary, so topological results guaranteeing the existence of global sections in the unitary case still hold in the non-unitary case. Hence, if M is simply connected, then a characterization of the existence of global sections follows from [14], but here we are not assuming any topological restriction on M. A concise characterization in terms of \(\tau \) and \(\mu \) will be given in Proposition 3.3 below. At this stage, we only need to guarantee the existence of sections in non-compact open subsets of M, which follows from [11, Section VIII.5] after changing the metric, or when the fibers have infinite length [25, Theorem 12.2].

Next we exhibit some distinguished examples of Killing submersions.

Example 2.3

If \(\tau =0\), then the horizontal distribution is integrable and we obtain the warped product \(M\times _\mu \mathbb {R}\) with one-dimensional fibers, which is the product manifold \(M\times \mathbb {R}\) endowed with the Riemannian metric \(\pi _M^*(\mathrm {d}s_M^2)+\mu ^2\pi _\mathbb {R}^*(\mathrm {d}t^2)\), where \(\mu \) is a function not depending on t, and \(\pi _M\) and \(\pi _\mathbb {R}\) denote the usual projections. If \(\mu \) is constant, we get the Riemannian product space \(M\times \mathbb {R}\).

Example 2.4

Every homogeneous 3-manifold X homeomorphic to \(\mathbb {R}^3\) is isometric to the semi-direct product \(\mathbb {R}^2\ltimes _A\mathbb {R}\) (for some real \(2\times 2\) matrix A) equipped with some left-invariant metric (see [17] for a detailed description of these metrics). This metric can be chosen such that the left-invariant frame

$$\begin{aligned} E_1&=\alpha _{11}(z)\partial _x+\alpha _{21}(z)\partial _y,&E_2&=\alpha _{12}(z)\partial _x+\alpha _{22}(z)\partial _y,&E_3&=\partial _z, \end{aligned}$$

is orthonormal (here \(\alpha _{ij}(z)\) denote the entries of the exponential matrix \(e^{zA}\) and (xyz) represents the usual coordinates in \(\mathbb {R}^3\)). The vector field \(\partial _x\) is right-invariant and Killing, and its integral curves are the fibers of the Killing submersion \((x,y,z)\mapsto (y,z)\). The metric can be expressed as

$$\begin{aligned} \frac{1}{\alpha _{22}^2+\alpha _{21}^2}\mathrm {d}y^2+\mathrm {d}z^2+\frac{\alpha _{22}^2+\alpha _{21}^2}{(\alpha _{11}\alpha _{22}-\alpha _{12}\alpha _{21})^2}\left( \mathrm {d}x-\frac{\alpha _{11}\alpha _{21}+\alpha _{12}\alpha _{22}}{\alpha _{22}^2+\alpha _{21}^2}\mathrm {d}y\right) ^2. \end{aligned}$$

The base surface is \(\mathbb {R}^2\) equipped with the metric \((\alpha _{22}^2+\alpha _{21}^2)^{-1}\mathrm {d}y^2+\mathrm {d}z^2\), and the bundle curvature and Killing length are determined by

$$\begin{aligned} \frac{2\tau }{\mu }&=\frac{\partial }{\partial z}\left( \frac{\alpha _{11}\alpha _{21}+\alpha _{12}\alpha _{22}}{\alpha _{22}^2+\alpha _{21}^2}\right) ,&\mu&=\sqrt{\frac{\alpha _{22}^2+\alpha _{21}^2}{\alpha _{11}\alpha _{22}-\alpha _{12}\alpha _{21}}}. \end{aligned}$$

Any choice of a right-invariant Killing vector field in X gives rise to Killing submersion structures in X. Although they are typically not isomorphic, all of them are topologically the usual projection \(\mathbb {R}^3\rightarrow \mathbb {R}^2\).

Example 2.5

If X is an homogeneous 3-manifold homeomorphic to \(\mathbb {S}^3\), then X is isometric to the three-dimensional Lie group \(\mathrm {SU}(2)\) equipped with some left-invariant metric, and any right-invariant \(\xi \in \mathfrak {X}(X)\) is Killing. If \(\xi \) is not identically zero, then it has no zeros and its integral curves are compact [16]. The space of fibers M is topologically \(\mathbb {S}^2\) and can be easily endowed with the structure of a Riemannian surface such that the natural projection from X to M is a Killing submersion (it suffices to induce in M the metric in the distribution orthogonal to \(\xi \)). This submersion is topologically the Hopf fibration \(\pi _{\mathrm {H}}:\mathbb {S}^3\rightarrow \mathbb {S}^2\).

Examples 2.4 and 2.5 cover all simply connected homogeneous 3-manifolds, except for the Riemannian products \(\mathbb {S}^2(\kappa )\times \mathbb {R}\), \(\kappa >0\), which do not admit any Lie group structures [17, Theorem 2.4]. Nevertheless, the usual projection \(\mathbb {S}^2(\kappa )\times \mathbb {R}\rightarrow \mathbb {S}^2(\kappa )\) is a Killing submersion with \(\tau =0\) and constant \(\mu \).

2.1 Killing submersions over a disk

We obtain the following direct generalization of Theorems 2.8 and 4.2 in [14]. We will give a detailed proof since it provides a constructive and explicit way of producing Killing submersions with prescribed bundle curvature and Killing length.

Theorem 2.6

(Local classification of Killing submersions) Let M be a non-compact simply connected Riemannian surface, and let \(\tau ,\mu \in \mathcal C^\infty (M)\), \(\mu >0\). Then there exists a Killing submersion \(\pi :\mathbb {E}\rightarrow M\) such that

  1. (1)

    the fibers of \(\pi \) have infinite length,

  2. (2)

    \(\tau \) is the bundle curvature of \(\pi \), and

  3. (3)

    \(\mu \) is the length of a Killing field \(\xi \) whose integral curves are the fibers of \(\pi \).

Moreover, such a Killing submersion \(\pi \) is unique up to isomorphism.

Proof

Since M is non-compact and simply connected, there is an isometry \(\varphi :(\Omega ,\mathrm {d}s^2_\lambda )\rightarrow M\), where \(\Omega \subset \mathbb {R}^2\) is the unit disk or the whole \(\mathbb {R}^2\) and \(\mathrm {d}s^2_\lambda =\lambda ^2(\mathrm {d}x^2+\mathrm {d}y^2)\) for some positive \(\lambda \in \mathcal C^\infty (\Omega )\). We will begin by describing how to obtain global coordinates for any Killing submersion over M in terms of \(\lambda \), which will give the key ideas to prove existence and uniqueness afterward.

Let us suppose that \(\pi :\mathbb {E}\rightarrow M\) is a Killing submersion satisfying (1), (2), and (3). The condition (1) yields the existence of a global smooth section \(F_0:M\rightarrow \mathbb {E}\) (see [25, Theorem 12.2]), and it also implies that

$$\begin{aligned} \Psi :\Omega \times \mathbb {R}\rightarrow \mathbb {E},\quad \Psi ((x,y),t)=\phi _t(F_0(\varphi (x,y))) \end{aligned}$$

is a global diffeomorphism, so the following diagram is commutative

where \(\pi _1:\Omega \times \mathbb {R}\rightarrow \Omega \) is the projection over the first factor. Now we can induce in \(\Omega \times \mathbb {R}\) a metric \(\mathrm {d}s^2\) making \(\Psi \) an isometry, so \(\pi _1\) becomes a Killing submersion over \((\Omega ,\mathrm {d}s_\lambda ^2)\). The orthonormal frame \(\{e_1=\frac{1}{\lambda }\partial _x,e_2=\frac{1}{\lambda }\partial _y\}\) in \((\Omega ,\mathrm {d}s_\lambda ^2)\) lifts via \(\pi _1\) to a global frame \(\{E_1,E_2\}\) of the horizontal distribution associated with \(\pi _1\). Hence, \(\{E_1,E_2,E_3=\frac{1}{\mu }\partial _t\}\) forms a global orthonormal frame of \(\mathbb {E}\). Since \(\pi _1(x,y,t)=(x,y)\), there exist \(a,b\in \mathcal {C}^\infty (\Omega )\) such that we can write

$$\begin{aligned} \begin{aligned} (E_1)_{(x,y,t)}&=\tfrac{1}{\lambda (x,y)}\,\partial _x+a(x,y)\partial _t,\\ (E_2)_{(x,y,t)}&=\tfrac{1}{\lambda (x,y)}\,\partial _y+b(x,y)\partial _t,\\ (E_3)_{(x,y,t)}&=\tfrac{1}{\mu (x,y)}\,\partial _t. \end{aligned} \end{aligned}$$
(2.2)

In other words, the metric \(\mathrm {d}s^2\) can be expressed in coordinates as

$$\begin{aligned} \mathrm {d}s^2=\lambda ^2(\mathrm {d}x^2+\mathrm {d}y^2)+\mu ^2\left( \mathrm {d} t-\lambda (a\mathrm {d}x+b\mathrm {d}y)\right) ^2. \end{aligned}$$
(2.3)

Note that \(\partial _t=\mu E_3\) represents the vertical Killing vector field, and there is no loss of generality in supposing that \(\{E_1,E_2,E_3\}\) is positively oriented in \(\mathbb {E}\).

By means of Eq. (2.1), the bundle curvature \(\tau \) of \(\pi _1\) reads

$$\begin{aligned} \tau =\tfrac{-1}{\mu }\omega (E_1,E_2)=-\langle \nabla _{E_1}E_3,E_2\rangle =\langle \nabla _{E_1}E_2,E_3\rangle . \end{aligned}$$
(2.4)

Likewise, using \(E_2\) rather than \(E_1\), we get \(\tau =-\langle \nabla _{E_2}E_1,E_3\rangle \), and adding this to (2.4), we finally obtain

$$\begin{aligned} \tau =\tfrac{1}{2}\langle [E_1,E_2],E_3\rangle =\tfrac{\mu }{2\lambda ^2}\left( (\lambda b)_x-(\lambda a)_y\right) . \end{aligned}$$
(2.5)

Note that the Lie bracket that can be easily computed from (2.2) as

$$\begin{aligned} \begin{aligned}~ [E_1,E_2]&=\tfrac{\lambda _y}{\lambda ^2}E_1-\tfrac{\lambda _x}{\lambda ^2}E_2+\tfrac{\mu }{\lambda ^2}\left( (\lambda b)_x-(\lambda a)_y\right) E_3,\\ [E_1,E_3]&=\tfrac{-\mu _x}{\lambda \mu }E_3,\qquad \qquad [E_2,E_3]=\tfrac{-\mu _y}{\lambda \mu }E_3. \end{aligned} \end{aligned}$$
(2.6)

Using this representation of Killing submersions, the proof follows easily.

Uniqueness. Let us suppose that \(\pi ':\mathbb {E}'\rightarrow M\) is another Killing submersion satisfying the same conditions and construct likewise an isometry \(\Psi ':\Omega \times \mathbb {R}\rightarrow \mathbb {E}'\) which induces functions \(a',b'\in \mathcal C^\infty (\Omega )\) as above. Since \(\tau \) and \(\mu \) coincide for both submersions, Equation (2.5) yields \((\lambda b')_x-(\lambda a')_y=(\lambda b)_x-(\lambda a)_y\). Equivalently, the following identity holds in \(\Omega \):

$$\begin{aligned} (\lambda (b'-b))_x=(\lambda (a'-a))_y. \end{aligned}$$

As \(\Omega \) is simply connected, Poincaré’s lemma guarantees the existence of \(d\in \mathcal C^\infty (\Omega )\) such that \(\lambda (b'-b)=d_y\) and \(\lambda (a'-a)=d_x\). It is not difficult to check that \(R:\Omega \times \mathbb {R}\rightarrow \Omega \times \mathbb {R}\) given by \(R((x,y),t)=((x,y),t-d(x,y))\) satisfies that \(T=\Psi '\circ R\circ \Psi ^{-1}:\mathbb {E}\rightarrow \mathbb {E}'\) is an isometry such that \(\pi '\circ T=\pi \).

Existence. We will show an explicit metric realizing the conditions in the statement, whose expression is inspired by the ideas above. Let us take

$$\begin{aligned} \mathbb {E}=\bigl (\Omega \times \mathbb {R},\mathrm {d}s^2=\lambda ^2(\mathrm {d}x^2+\mathrm {d}y^2)+\mu ^2(\mathrm {d}t+\eta (y\mathrm {d}x-x\mathrm {d}y))^2\bigr ), \end{aligned}$$
(2.7)

where

$$\begin{aligned} \eta =\eta (x,y)=\int _0^1\frac{2s\,\tau (sx,sy)\,\lambda (sx,sy)^2}{\mu (sx,sy)}\mathrm {d}s. \end{aligned}$$
(2.8)

Then \(\pi :\mathbb {E}\rightarrow (\Omega ,\mathrm {d}s_\lambda ^2)\) given by \(\pi (x,y,z)=(x,y)\) becomes a Riemannian submersion with Killing vector field \(\partial _t\) (the coefficients of \(\mathrm {d}s^2\) do not depend upon t) such that \(\Vert \partial _t\Vert =\mu \). Finally (2.7) arises from (2.3) for \(a=\frac{-y\eta }{\lambda }\) and \(b=\frac{x\eta }{\lambda }\). Hence, (2.5) tells us that the bundle curvature associated with (2.7) is \(\frac{\mu }{2\lambda ^2}((x\eta )_x+(y\eta )_y)\), which is equal to \(\tau \) by a simple calculation using (2.8). \(\square \)

Remark 2.7

Condition (1) in the statement is not restrictive at all. If the fibers of \(\pi :\mathbb {E}\rightarrow M\) have finite length, then \(\pi \) can be recovered as a Riemannian quotient of a Killing submersion with fibers of infinite length under an appropriate vertical translation. This is a consequence of Theorem 2.13 below. In particular, the total space of a Killing submersion over a disk is diffeomorphic to \(\mathbb {R}^3\) or to \(\mathbb {R}^2\times \mathbb {S}^1\).

This also implies that the condition (1) in the statement may be replaced by assuming that \(\mathbb {E}\) is simply connected.

2.2 Killing submersions over the sphere

The classification result when the base surface is topologically \(\mathbb {S}^2\) can be carried out in a similar fashion to the unitary case (see [14, Propositions 4.5 and 4.9]). The key ingredients are the horizontal lifts of curves and their holonomy properties, which we formulate next.

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion and \(\alpha :[a,b]\rightarrow M\) a curve of class \(\mathcal {C}^1\). A horizontal lift of \(\alpha \) is a curve \(\widetilde{\alpha }:[a,b]\rightarrow \mathbb {E}\) of class \(\mathcal {C}^1\) such that \(\widetilde{\alpha }'\) is always horizontal and \(\pi \circ \widetilde{\alpha }=\alpha \) in [ab]. The horizontal lift always exists, and it is unique when the point \(\widetilde{\alpha }(a)\) is chosen in the fiber of \(\alpha (a)\).

Proposition 2.8

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion whose fibers have infinite length, and let \(\alpha :[a,b]\rightarrow M\) be a simple \(\mathcal {C}^1\)-curve bounding an orientable relatively compact open set \(O\subset M\). Given a horizontal lift \(\widetilde{\alpha }\) of \(\alpha \), there exists a unique \(d\in \mathbb {R}\) such that \(\phi _d(\widetilde{\alpha }(a))=\widetilde{\alpha }(b)\) and it satisfies

$$\begin{aligned} \left| \int _O\frac{2\tau }{\mu }\right| =|d|. \end{aligned}$$

It is worth highlighting that the proof of Proposition 2.8 relies on the divergence theorem applied to the divergence-type expression for \(\tau \) given by Eq. (2.1) (see also [14, Proposition 3.3]).

Theorem 2.9

(Classification of Killing submersions over the sphere) Let M be \(\mathbb {S}^2\) endowed with some Riemannian metric, and consider \(\tau ,\mu \in \mathcal C^\infty (\mathbb {S}^2)\), with \(\mu >0\). Then there exists a Killing submersion \(\pi :\mathbb {E}\rightarrow M\) such that

  1. (1)

    \(\mathbb {E}\) is simply connected,

  2. (2)

    \(\tau \) is the bundle curvature of \(\pi \), and

  3. (3)

    \(\mu \) is the length of a Killing field \(\xi \) whose integral curves are the fibers of \(\pi \).

Moreover, such a Killing submersion \(\pi \) is unique up to isomorphism.

  1. (a)

    If \(\int _M\frac{\tau }{\mu }=0\), then the length of the fibers of \(\pi \) is infinite and \(\pi \) is isomorphic to

    $$\begin{aligned} \pi _1:(\mathbb {S}^2\times \mathbb {R},\mathrm {d}s^2)\rightarrow \mathbb {S}^2,\qquad \pi _1(p,t)=p, \end{aligned}$$

    for some Riemannian metric \(\mathrm {d}s^2\), with Killing vector field \(\xi _{(p,t)}=\partial _t\).

  2. (b)

    If, on the contrary, \(\int _M\frac{\tau }{\mu }\ne 0\), then the fibers of \(\pi \) have finite length and \(\pi \) is isomorphic to the Hopf fibration

    $$\begin{aligned} \pi _{\mathrm {H}}:(\mathbb {S}^3,\mathrm {d}s^2)\rightarrow \mathbb {S}^2,\qquad \pi _{\mathrm {H}}(z,w)=(2z\bar{w},|z|^2-|w|^2), \end{aligned}$$

    for some Riemannian metric \(\mathrm {d}s^2\) in \(\mathbb {S}^3\) with Killing vector field \(\xi _{(z,w)}=(iz,iw)\). Here \(\mathbb {S}^3\) and \(\mathbb {S}^2\) are the unit spheres in \(\mathbb {C}^2\) and \(\mathbb {R}^3\equiv \mathbb {C}\times \mathbb {R}\), respectively.

Remark 2.10

If condition (1) in the statement of Theorem 2.9 is dropped, then we obtain a quotient under a vertical translation. This follows from Theorem 2.13 below. In the trivial case, the total space of a Killing submersion over a sphere is diffeomorphic to \(\mathbb {S}^2\times \mathbb {R}\) or \(\mathbb {S}^2\times \mathbb {S}^1\). In the non-trivial case, the quotient must be taken under a translation of appropriate length, from where the total space must be diffeomorphic to \(\mathbb {S}^3\) or to the lens space L(n, 1), \(n\ge 2\). Note that these quotients of \(\mathbb {S}^3\) are orientable (L(2, 1) is nothing but the real projective space \(\mathbb {RP}^3\)).

2.3 Killing submersions over a non-simply connected surface

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion over an arbitrary orientable surface M with bundle curvature \(\tau \) and Killing length \(\mu \). Our goal is to show that \(\pi \) can be regarded as the quotient of a Killing submersion over a simply connected surface, which have been classified in the previous sections.

Let \(\rho :\widetilde{M}\rightarrow M\) be the Riemannian universal covering of M, and consider \(\widetilde{\tau }=\tau \circ \rho \) and \(\widetilde{\mu }=\mu \circ \rho \) the lift to \(\widetilde{M}\) of \(\tau \) and \(\mu \), respectively. In view of Theorems 2.6 and 2.9, there exists a unique Killing submersion \(\widetilde{\pi }:\widetilde{\mathbb {E}}\rightarrow \widetilde{M}\) with bundle curvature \(\widetilde{\tau }\), Killing length \(\widetilde{\mu }\), and such that \(\widetilde{\mathbb {E}}\) is simply connected. Let us denote by \(\mathrm {Aut}(\rho )\) the group of automorphisms of \(\rho \) (i.e., \(\mathrm {Aut}(\rho )\) consists of those isometries \(h:\widetilde{M}\rightarrow \widetilde{M}\) such that \(\rho \circ h=\rho \)). It is well known that M is the Riemannian quotient of \(\widetilde{M}\) under the properly discontinuous action of \(\mathrm {Aut}(\rho )\).

We will use \(\mathrm {Aut}(\rho )\) to construct a group G acting properly discontinuously by isometries on \(\widetilde{\mathbb {E}}\), and such that \(\mathbb {E}=\widetilde{\mathbb {E}}/G\), so we will need a way of producing isometries in \(\widetilde{\mathbb {E}}\) by lifting isometries in \(\mathrm {Aut}(\rho )\). We will omit the proof of the following lemma since it is a direct generalization of [14, Theorem 2.8].

Lemma 2.11

Given \(h\in \mathrm {Aut}(\rho )\) and \(p,q\in \widetilde{\mathbb {E}}\) such that \(h(\widetilde{\pi }(p))=\widetilde{\pi }(q)\), there exists a unique isometry \(f:\widetilde{\mathbb {E}}\rightarrow \widetilde{\mathbb {E}}\) such that \(\widetilde{\pi }\circ f=h\circ \widetilde{\pi }\) and \(f(p)=q\).

Remark 2.12

An isometry of the total space of a Killing submersion preserving the vertical direction is called a Killing isometry. It is straightforward that any Killing isometry \(f:\widetilde{\mathbb {E}}\rightarrow \widetilde{\mathbb {E}}\) induces an isometry \(h:\widetilde{M}\rightarrow \widetilde{M}\) such that \(\widetilde{\tau }\circ h=\epsilon \widetilde{\tau }\) and \(\widetilde{\mu }\circ h=a\widetilde{\mu }\) for some \(\epsilon \in \{-1,1\}\) and \(a\in \mathbb {R}\). Lemma 2.11 is the converse when \(\epsilon =a=1\): Given the isometry \(h:\widetilde{M}\rightarrow \widetilde{M}\) preserving \(\widetilde{\tau }\) and \(\widetilde{\mu }\), it can be lifted to some \(f:\widetilde{\mathbb {E}}\rightarrow \widetilde{\mathbb {E}}\), which is unique up to composing with vertical translations. In fact, since \(\epsilon =1\), it can be shown that such an f preserves the orientation in \(\widetilde{\mathbb {E}}\).

Theorem 2.13

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion, and let \(\rho :\widetilde{M}\rightarrow M\) and \(\sigma :\widetilde{\mathbb {E}}\rightarrow \mathbb {E}\) be the universal Riemannian covering maps of M and \(\mathbb {E}\), respectively.

  1. (a)

    There exist a Killing submersion \(\widetilde{\pi }:\widetilde{\mathbb {E}}\rightarrow \widetilde{M}\) such that \(\rho \circ \widetilde{\pi }=\pi \circ \sigma \).

  2. (b)

    There exists a group G of Killing isometries acting properly discontinuously on \(\widetilde{\mathbb {E}}\) such that \(\mathbb {E}=\widetilde{\mathbb {E}}/G\).

Moreover, each \(f\in G\) is associated with some \(h\in \mathrm {Aut}(\rho )\), so the diagram

is commutative. If the fibers of \(\pi \) have infinite length, this correspondence between G and \(\mathrm {Aut}(\rho )\) is bijective. Otherwise, any two isometries associated with the same \(h\in \mathrm {Aut}(\rho )\) differ in a vertical translation of length a multiple of the length of the fiber.

Proof

Let us consider \(\widetilde{\pi }:\widetilde{\mathbb {E}}\rightarrow \widetilde{M}\) as defined above. It suffices to construct the group G acting properly discontinuously on \(\widetilde{\mathbb {E}}\) and such that \(\widetilde{\mathbb {E}}/G=\mathbb {E}\), so the projection to the quotient \(\sigma :\widetilde{\mathbb {E}}\rightarrow \mathbb {E}\) will enjoy the desired properties.

Let us fix \(p_0\in \widetilde{\mathbb {E}}\) and \(q_0\in \mathbb {E}\) with \(\pi (q_0)=\rho (\widetilde{\pi }(p_0))\). Given \(h\in \mathrm {Aut}(\rho )\), we will begin by defining the Killing isometries associated with h. In order to find the initial conditions required by Lemma 2.11, horizontal lifts come in handy.

Let \(\gamma :[0,1]\rightarrow \widetilde{M}\) be a smooth path with \(\gamma (0)=\widetilde{\pi }(p_0)\) and \(\gamma (1)=h(\widetilde{\pi }(p_0))\). Then \(\alpha =\rho \circ \gamma \) is a closed curve in M with \(\alpha (0)=\alpha (1)=\rho (\widetilde{\pi }(p_0))\). Let us consider \(\widehat{\alpha }:[0,1]\rightarrow \mathbb {E}\) the horizontal lift of \(\alpha \) such that \(\widehat{\alpha }(0)=q_0\), and take

$$\begin{aligned} A(h)=\{t\in \mathbb {R}:\phi _t(q_0)=\widehat{\alpha }(1)\}, \end{aligned}$$

where \(\{\phi _t\}\) are the vertical translations in \(\mathbb {E}\). For each \(t\in A(h)\), Lemma 2.11 yields the existence of an isometry \(f:\widetilde{\mathbb {E}}\rightarrow \widetilde{\mathbb {E}}\) such that \(\widetilde{\pi }\circ f=h\circ \widetilde{\pi }\) and \(f(p_0)=\widetilde{\phi }_{-t}(\widehat{\gamma }(1))\), where \(\{\widetilde{\phi }_t\}\) is the group of vertical translations in \(\widetilde{\mathbb {E}}\) and \(\widehat{\gamma }:[0,1]\rightarrow \widetilde{\mathbb {E}}\) is the horizontal lift of \(\gamma \) such that \(\widehat{\gamma }(0)=p_0\). This choice of the initial condition is motivated by the fact that the vertical distance between the endpoints of the horizontal lift of a closed curve must be preserved by these (local) isomorphisms of Killing submersions.

We define G as the set of all the Killing isometries f constructed above when \(t\in A(h)\) and \(h\in \mathrm {Aut}(\rho )\). It is easy to check that G does not depend on the choice of the paths \(\gamma \) as a consequence of Proposition 2.8. Moreover, G is a group acting properly discontinuously on \(\widetilde{\mathbb {E}}\), with quotient space \(\widetilde{\mathbb {E}}/G=\mathbb {E}\). Observe that, if there are two isometries \(f_1\) and \(f_2\) associated with the same h, then \(f_1\circ f_2^{-1}=\widetilde{\phi }_t\) for some \(t\in \mathbb {R}\) such that \(\phi _t(q_0)=q_0\). Hence, \(t=0\) if the fibers have infinite length (so \(f_1=f_2\)) or \(\widetilde{\phi }_t\) is a vertical translation of length a multiple of the length of the fiber, if this length is finite. \(\square \)

Theorem 2.13 yields a constructive way of producing non-simply connected Killing submersions as Riemannian quotients of the simply connected ones by subgroups of Killing isometries. The simplest examples with non-trivial topology arise as quotients of the Heisenberg group, which we analyze next.

Example 2.14

(Non-uniqueness) Let us consider the Heisenberg group \(\mathrm {Nil}_3(\tau )\), which admits a structure of Killing submersion over the Euclidean plane \(\mathbb {R}^2\) with constant \(\tau >0\) and \(\mu =1\), so it is isometric to \(\mathbb {R}^3\) with the metric given by (2.7) for \(\lambda =\mu =1\) and \(\eta =\tau \).

Given \(a, b\in \mathbb {R}\), the isometries of \(\mathrm {Nil}_3(\tau )\) defined by

$$\begin{aligned} f_1(x,y,z)&=(x+1,y,z+\tau y+a),&f_2(x,y,z)&=(x,y+1,z-\tau x+b), \end{aligned}$$

correspond to the translations \(h_1(x,y)=(x+1,y)\) and \(h_2(x,y)=(x,y+1)\), respectively. Observe that the commutator \(f_1\circ f_2\circ f_1^{-1}\circ f_2^{-1}\) maps (xyz) into \((x,y,z+2\tau )\); thus, the group G spanned by \(f_1\) and \(f_2\) has compact quotient \(\mathbb {E}=\mathrm {Nil}_3(\tau )/G\) admitting a Killing submersion structure \(\pi :\mathbb {E}\rightarrow \mathbb {T}=\mathbb {R}^2/H\), where H is the group of isometries of \(\mathbb {R}^2\) spanned by \(h_1\) and \(h_2\). Theorem 2.13 guarantees that these are all Killing submersions over the flat torus \({\mathbb {T}}\) with bundle curvature \(\tau \) and unit Killing vector field.

Let us consider the closed curve \(\alpha :[0,1]\rightarrow {\mathbb {T}}\), \(\alpha (t)=(t,0)\), whose horizontal lift is \(\widetilde{\alpha }(t)=(t,0,0)\). Since (0, 0, 0) and (1, 0, a) are identified by \(f_1\), the vertical distance from \(\widetilde{\alpha }(0)\) to \(\widetilde{\alpha }(1)\) is \(2\tau -a\). As isomorphisms of Killing submersions preserve the horizontal lift of a curve as well as vertical distances, we conclude that different values of a and b give rise to non-isomorphic Killing submersions over \({\mathbb {T}}\) with the same bundle curvature and Killing length.

3 The mean curvature equation

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion and fix a Killing vector field \(\xi \) whose integral curves are the fibers of \(\pi \), which will be supposed to have infinite length throughout this section.

Let \(U\subseteq M\) be an open subset such that there exists a smooth section \(F_0:U\rightarrow \mathbb {E}\). Given \(u:U\rightarrow \mathbb {R}\), we define the Killing graph associated with u with respect to \(F_0\) as the surface in \(\mathbb E\) parametrized by

$$\begin{aligned} F_u:U\rightarrow \mathbb {E},\quad F_u(p)=\phi _{u(p)}(F_0(p)), \end{aligned}$$

where \(\{\phi _t\}\) is the 1-parameter group of isometries associated with \(\xi \). Next lemma gives a divergence-type formula for the mean curvature (generalizing the corresponding formula in the unitary case, see [12, 13]). We will assume that all functions are smooth in the sequel, despite the fact that most of the following arguments work for \(\mathcal C^2\)-graphs.

Lemma 3.1

Given \(u\in \mathcal {C}^\infty (U)\), the mean curvature H of \(F_u\) with respect to a unit normal vector field N along \(F_u\), satisfies (as a function in M)

$$\begin{aligned} 2H\mu =\mathrm {div}(\mu \,\pi _*N), \end{aligned}$$

where \(\mathrm {div}\) denotes the divergence on M and \(\mu =\Vert \xi \Vert \) is the Killing length.

Proof

It is well known that \(2H=\overline{\mathrm {div}}(\overline{N})\), where \(\overline{\mathrm {div}}\) is the ambient divergence and \(\overline{N}\) is an extension to \(\mathbb {E}\) of N. Here we will consider \(\overline{N}\) as the unique extension of N satisfying \((\phi _t)_*\overline{N}=\overline{N}\) for all \(t\in \mathbb R\), i.e., \(\overline{N}\) is constant along the fibers of \(\pi \). Working locally if necessary, we can consider an orthonormal frame \(\{e_1,e_2\}\) in U and its horizontal lift \(\{E_1,E_2\}\), which is completed to an orthonormal frame in \(\pi ^{-1}(U)\) by adding \(E_3=\frac{1}{\mu }\xi \), and

$$\begin{aligned} 2H=\overline{\mathrm {div}}(\overline{N})= \sum _{i=1}^2\langle \overline{\nabla }_{E_i}\overline{N},E_i\rangle +\langle \overline{\nabla }_{E_3}\overline{N},E_3\rangle . \end{aligned}$$
(3.1)

Since \(\pi \) is Riemannian, \(E_1\) and \(E_2\) are horizontal, and \(\overline{N}\) is invariant under vertical translations, we obtain that \(\sum _{i=1}^2\langle \overline{\nabla }_{E_i}\overline{N},E_i\rangle =\mathrm {div}_M(\pi _*N)\). Furthermore,

$$\begin{aligned} \langle \overline{\nabla }_{E_3}\overline{N},E_3\rangle =\langle [E_3,\overline{N}],E_3\rangle =\tfrac{1}{\mu ^2}\left\langle [\xi ,\overline{N}],\xi \right\rangle -\tfrac{1}{\mu }\bigl \langle \overline{N}\bigl (\tfrac{1}{\mu }\bigr )\xi ,\xi \bigr \rangle =\tfrac{1}{\mu }\overline{N}(\mu ). \end{aligned}$$
(3.2)

The first equality follows from Koszul formula and the fact that \(\langle \overline{N},E_3\rangle \) is constant along fibers, the second one follows from the fact that \(E_3=\frac{1}{\mu }\xi \). In the third one, we used that \(\xi \) is Killing to obtain that

$$\begin{aligned} \langle [\overline{N},\xi ],\xi \rangle&=\langle \overline{\nabla }_{\overline{N}}\xi ,\xi \rangle -\langle \overline{\nabla }_\xi \overline{N},\xi \rangle =-\langle \overline{\nabla }_\xi \xi ,\overline{N}\rangle -\langle \overline{\nabla }_\xi \overline{N},\xi \rangle =-\xi (\langle \overline{N},\xi \rangle )=0. \end{aligned}$$

Plugging (3.2) into (3.1), considering H and \(\mu \) as functions on M (recall that \(\mu \) is constant along fibers), and denoting the gradient in M by \(\nabla \), we get

$$\begin{aligned} 2H=\mathrm {div}(\pi _*N)+\tfrac{1}{\mu }\overline{N}(\mu )=\mathrm {div}(\pi _*N)+\tfrac{1}{\mu }\langle \pi _*N,\nabla \mu \rangle _M=\tfrac{1}{\mu }\mathrm {div}\left( \mu \ \pi _*N\right) , \end{aligned}$$

so the statement follows. \(\square \)

Now, in order to work out the term \(\pi _*N\), let us consider \(\overline{u}\in \mathcal C^\infty (\mathbb {E})\) the extension of u by making it constant along fibers, and \(d\in \mathcal C^\infty (\mathbb {E})\) the function that measures the signed vertical Killing distance from a point to the point in \(F_0\) lying on the same fiber. In other words, d is determined by the identity \(\phi _{d(p)}(F_0(\pi (p)))=p\), for all \(p\in \mathbb {E}\). Then the surface parametrized by \(F_u\) is a level surface of the function \(\overline{u}-d\in \mathcal C^\infty (\mathbb {E})\) so we can compute \(N=\overline{\nabla }(\overline{u}-d)/\Vert \overline{\nabla }(\overline{u}-d)\Vert \). If we set \(Z=\pi _*(\overline{\nabla }d)\), it is not difficult to get to

$$\begin{aligned} \pi _*N=\frac{Gu}{\sqrt{\mu ^{-2}+\Vert Gu\Vert ^2}},\qquad Gu=\nabla u-Z, \end{aligned}$$
(3.3)

where the gradient and the norm are computed in M. The denominator in (3.3) is the area element of \(F_u\), i.e., the area of the surface spanned by u is equal to

$$\begin{aligned} \mathcal {A}(u)=\int _U \sqrt{\mu ^{-2}+\Vert Gu\Vert ^2}. \end{aligned}$$
(3.4)

The information about the bundle curvature as well as about the choice of the zero section \(F_0\) is encoded in Z.

Lemma 3.2

The vector field Z satisfies

$$\begin{aligned} \mathrm {div}(JZ)=-\frac{2\tau }{\mu }, \end{aligned}$$

where J denotes a \(\frac{\pi }{2}\)-rotation in TM.

Proof

As the computation is local, we can suppose the Killing submersion is the projection over the first factor \(\pi _1:\Omega \times \mathbb {R}\rightarrow \Omega \), such that \(\{e_1=\frac{1}{\lambda }\partial _x, e_2=\frac{1}{\lambda }\partial _y\}\) and \(\{E_1,E_2,E_3\}\) given by (2.2) are orthonormal frames in \(\Omega \) and \(\Omega \times \mathbb {R}\), respectively, for some \(a,b,\lambda \in \mathcal C^\infty (\Omega )\), \(\lambda >0\). Taking the global initial section \(F_0:\Omega \rightarrow \Omega \times \mathbb {R}\) as \(F_0(x,y)=(x,y,0)\), we get that \(\overline{u}(x,y,z)=u(x,y)\) and \(d(x,y,z)=z\) in this model. Hence,

$$\begin{aligned} \overline{\nabla }d=\sum _{i=1}^3E_i(z)E_i=a E_1+b E_2+\tfrac{1}{\mu }E_3, \end{aligned}$$

which gives \(Z=\pi _*(\overline{\nabla }d)=a e_1+b e_2\). As \(\{e_1,e_2\}\) is orthonormal (and can be assumed positively oriented), we deduce that \(JZ=-be_1+ae_2=-\frac{b}{\lambda }\partial _x+\frac{a}{\lambda }\partial _y\). Using the expression for the divergence in a conformal metric and Eq. (2.5), we finally obtain

$$\begin{aligned} \mathrm {div}(JZ)=\frac{1}{\lambda ^2}\left( -(\lambda b)_x+(\lambda a)_y\right) =-\frac{2\tau }{\mu }. \end{aligned}$$

\(\square \)

Now we are able to characterize the existence of global sections.

Proposition 3.3

(Existence of global sections) Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion with bundle curvature \(\tau \) and Killing length \(\mu \). The submersion \(\pi \) admits a global section if and only if either M is non-compact or M is compact and \(\int _{M}\frac{\tau }{\mu }=0\).

Proof

If M is not compact, the existence follows from [11, Section VIII.5], so let us suppose that M is compact. If M admits a global section, then \(\int _{M}\frac{\tau }{\mu }=0\) as a consequence of Lemma 3.2 and the divergence theorem. Conversely, let us assume that \(\int _{M}\frac{\tau }{\mu }=0\) and take a regular Jordan curve \(\Gamma \subset M\) separating M in two components \(M_1\) and \(M_2\) such that \(\int _{M_1}\frac{\tau }{\mu }=\int _{M_2}\frac{\tau }{\mu }=0\). By Proposition 2.8, \(\Gamma \) lifts to a horizontal closed curve \(\widetilde{\Gamma }\subset \mathbb {E}\). Since each component \(M_i\), \(i\in \{1,2\}\), is not compact, there exist sections \(F_i:\overline{M_i}\rightarrow \mathbb {E}\) such that \(F_i(x)\in \widetilde{\Gamma }\) for all \(x\in \Gamma \). Then \(F: M\rightarrow \mathbb {E}\) given by \(F(x)=F_i(x)\) if \(x\in M_i\) is a global continuous section of \(\pi \) which can be assumed smooth after a perturbation around \(\Gamma \). \(\square \)

Before stating the existence of minimal sections, we need the following lemma, which collects some useful information given by the maximum principle for constant mean curvature surfaces, and will be used repeatedly in the sequel.

Lemma 3.4

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion with bundle curvature \(\tau \) and Killing length \(\mu \). If \(\Sigma \) is a compact surface immersed in \(\mathbb {E}\) with constant mean curvature H and everywhere transversal to the vertical direction, then M is compact, \(\int _M\frac{\tau }{\mu }=0\) and \(\Sigma \) is an entire minimal graph.

Moreover, if such an entire minimal graph exists, then it is unique up to vertical translations and minimizes area among all entire Lipschitz graphs in \(\mathbb {E}\).

Proof

The absolute value of the Jacobian of \(\pi _{|\Sigma }:\Sigma \rightarrow M\) equals \(\frac{1}{\mu }|\langle N,\xi \rangle |\), where N is a unit normal to the surface and \(\xi \) is a vertical Killing vector field. Since \(\Sigma \) is compact and transversal to \(\xi \), we can assume that \(|\langle N,\xi \rangle |\ge \varepsilon \) for some constant \(\varepsilon >0\). As a consequence, \(\pi _{|\Sigma }\) is a covering map with a finite number of sheets. In particular, M is compact and the bundle is trivial (i.e., it admits a global section), so we deduce that \(\int _M\frac{\tau }{\mu }=0\) from Proposition 3.3. Since the submersion is trivial, we do not lose generality if we assume that the fibers of \(\pi \) have infinite length (by considering a certain cover space if necessary).

To prove that \(\Sigma \) is an entire graph, let us assume by contradiction that \(\pi _{|\Sigma }\) has more than one sheet. Then take \(\Sigma '\) a vertical translation of \(\Sigma \) with \(\Sigma \cap \Sigma '=\emptyset \) and move \(\Sigma '\) toward \(\Sigma \) till there is a first contact point p. The normal vector fields to the surfaces coincide at p because \(\langle N,\xi \rangle \) has a global sign. Since there is more than one sheet, the point p comes from two different sheets of \(\Sigma \), and the maximum principle implies that \(\Sigma \) is invariant by a non-trivial vertical translation, contradicting the fact that it is compact and the fibers have infinite length.

To prove that \(\Sigma \) is minimal, Lemma 3.1 applied to the entire graph \(\Sigma \) and the divergence theorem give \(\int _M H\mu =0\), where \(\mu \) is positive, so \(H=0\). Uniqueness up to vertical translations follows from the maximum principle.

Finally, let us show the minimization property. Consider any entire Lipschitz graph \(\Sigma '\), and suppose it does not intersect \(\Sigma \) by applying a vertical translation. If N and \(N'\) stand for unit normal vector fields to \(\Sigma \) and \(\Sigma '\), respectively (\(N'\) is defined almost everywhere), then we can extend N to a unique \(\overline{N}\in \mathfrak {X}(\mathbb {E})\) such that \((\phi _t)_*\overline{N}=\overline{N}\) for all t, where \(\{\phi _t\}\) is the 1-parameter group of vertical translations. Then the ambient divergence \(\overline{\mathrm {div}}(\overline{N})\) identically vanishes since \(\overline{N}\) is normal to a foliation of \(\mathbb {E}\) by minimal surfaces, namely the translations of \(\Sigma \). The divergence theorem for \(\overline{N}\) in the Lipschitz region bounded by \(\Sigma \) and \(\Sigma '\) and Cauchy–Schwarz inequality yield

$$\begin{aligned} \mathrm {Area}(\Sigma )=\left| \int _\Sigma \langle N,N\rangle \right| =\left| \int _{\Sigma '}\langle N',N\rangle \right| \le \mathrm {Area}(\Sigma '), \end{aligned}$$

with equality if and only if \(N'\) is collinear with N almost everywhere, that is to say, if and only if \(\Sigma '\) is a vertical translation of \(\Sigma \). \(\square \)

Remark 3.5

Lemma 3.4 can be extended to the case the mean curvature of \(\Sigma \) does not change sign. Although the maximum principle cannot be applied directly because the mean curvatures of \(\Sigma \) and \(\Sigma '\) may not be properly ordered at the contact point, we could also consider and a Killing submersion over a Riemannian covering of M where \(\Sigma \) lifts to a entire graph. Since the mean curvature of the lifted surface does not change sign either, Lemma 3.1 gives the result.

Theorem 3.6

(Bernstein problem) Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion, and assume that M is compact and \(\int _M\frac{\tau }{\mu }=0\), where \(\tau \) and \(\mu \) denote the bundle curvature and the Killing length, respectively. Then \(\mathbb {E}\) admits an entire minimal graph \(\Sigma \).

  1. (1)

    Up to a vertical translation \(\Sigma \) is the unique entire graph in \(\mathbb {E}\) with constant mean curvature, and minimizes area among all entire graphs in \(\mathbb {E}\).

  2. (2)

    If the fibers of \(\pi \) have infinite length, then \(\Sigma \) is the unique compact minimal surface immersed in \(\mathbb {E}\), up to vertical translations.

Proof

Since M is compact and \(\int _M\frac{\tau }{\mu }=0\), Proposition 3.3 guarantees the existence of global sections, and we can assume that the fibers of \(\pi \) have infinite length without losing generality. The proof is now split into two cases depending of whether or not M is topologically a sphere.

Case I. If M is not homeomorphic to \(\mathbb {S}^2\), then let Q be a quotient of \(\mathbb {E}\) under a vertical translation, and consider \(\Sigma _0\subset Q\) to be a global smooth section. On the one hand, since M is not a sphere, the compact 3-manifold Q is irreducible because its universal cover is topologically \(\mathbb {R}^3\), which is irreducible. On the other hand, the global section \(\Sigma _0\) is incompressible in Q because the inclusion \(\iota : \Sigma _0\rightarrow P\) induces an injective morphism \(\iota _*:\pi _1(\Sigma _0)\rightarrow \pi _1(Q)\) between the homotopy groups. The work of Meeks, Simon, and Yau [20] proves the existence of an embedded minimal surface \(\Sigma '\subset Q\) in the isotopy class of \(\Sigma _0\).

Since \(\Sigma '\) is isotopic to \(\Sigma _0\), it is clear that it lifts to a compact embedded minimal surface \(\Sigma \subset \mathbb {E}\). Since the fibers of \(\pi \) have infinite length, it follows that \(\Sigma \) is an entire graph by an easy application of the maximum principle. The rest of properties in the statement are now a consequence of Lemma 3.4.

Case II. If M is homeomorphic to \(\mathbb {S}^2\), the technique above does not apply (note that \(M\times \mathbb {S}^1\) is not irreducible and a global section is not incompressible). Here we will develop an ad hoc argument based on a Calabi-type duality, and on the existence of an entire spacelike graph with prescribed mean curvature \(\tau \) in \(M\times \mathbb {R}\) when the Lorentzian metric \(\pi _M^*(\mathrm {d}s^2_M)-\mu ^{-2}\pi _\mathbb {R}^*(\mathrm {d}t^2)\) is considered. Such existence follows from the work of Gerhardt [10].

Parametrizing the surface given by Gerhardt as \(\Phi :M\mapsto M\times \mathbb {R}\) with \(\Phi (x)=(x,v(x))\) for some \(v\in \mathcal {C}^\infty (M)\), the fact that \(\Phi \) has prescribed mean curvature \(\tau \) with respect to the aforesaid Lorentzian metric can be expressed as

$$\begin{aligned} \mathrm {div}\left( \frac{\nabla v}{\mu \sqrt{\mu ^2-\Vert \nabla v\Vert ^2}}\right) =\frac{2\tau }{\mu }, \end{aligned}$$
(3.5)

and the spacelike condition reads \(\Vert \nabla v\Vert <\mu \). Given a global smooth section \(F_0:M\rightarrow \mathbb {E}\), define \(d\in \mathcal C^\infty (\mathbb {E})\) as the Killing distance along fibers to \(F_0\). The vector field \(Z=\pi _*(\overline{\nabla }d)\) in M satisfies \(\mathrm {div}(JZ)=\frac{-2\tau }{\mu }\) by Lemma 3.2, and therefore

$$\begin{aligned} \mathrm {div}\left( \frac{\nabla v}{\mu \sqrt{\mu ^2-\Vert \nabla v\Vert ^2}}+JZ\right) =0. \end{aligned}$$
(3.6)

Since M is simply connected, Poincaré’s lemma guarantees the existence of \(u\in \mathcal C^\infty (M)\) such that

$$\begin{aligned} \frac{\nabla v}{\mu \sqrt{\mu ^2-\Vert \nabla v\Vert ^2}}+JZ=J\nabla u. \end{aligned}$$
(3.7)

Defining \(Gu=\nabla u-Z\), and applying J to  (3.7), we get

$$\begin{aligned} Gu=\frac{-J\nabla v}{\mu \sqrt{\mu ^2-\Vert \nabla v\Vert ^2}}. \end{aligned}$$
(3.8)

Taking squared norms in (3.8), we reach

$$\begin{aligned} \sqrt{\mu ^2-\Vert \nabla v\Vert ^2}=\frac{1}{\sqrt{\mu ^{-2}+\Vert Gu\Vert ^2}}. \end{aligned}$$
(3.9)

Now plugging (3.9) into (3.8), leads to

$$\begin{aligned} \frac{\mu \cdot Gu}{\sqrt{\mu ^{-2}+\Vert Gu\Vert ^2}}=-J\nabla v, \end{aligned}$$
(3.10)

and taking the divergence in (3.10), it follows that the graph of u with respect to \(F_0\) is an entire minimal graph (note that \(\mathrm {div}(J\nabla v)=0\)). The remaining properties in the statement follow from Lemma 3.4 and the maximum principle. \(\square \)

Remark 3.7

If M is compact but not simply connected, Gerhardt’s result still holds. Passing to the universal Riemannian covering surface of M, an argument similar to that of Case II in the proof of Theorem 3.6 can be applied giving rise to a complete minimal surface in \(\mathbb {E}\) everywhere transversal to the Killing direction (a so-called vertical multigraph), but it may fail to be an entire graph because of a period problem.

In our opinion, the hypothesis that M is compact in Theorem 3.6 can be dropped. In that case, on the one hand, existence of a global section is guaranteed, so conditions on \(\tau \) or \(\mu \) are not necessary in principle; on the other hand, uniqueness may fail as the solutions to the Bernstein problem in \(\mathbb {R}^3\) or \(\mathrm {Nil}_3(\tau )\) show [4, 8]. Hence, we pose the following conjecture.

Conjecture 3.8

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion such that M is not compact. Then \(\mathbb {E}\) admits an entire minimal graph.

Apart from the trivial Killing submersions over \(\mathbb {S}^2\), Killing submersions with the topology of the Hopf projection also admit minimal spheres. This follows from a result of Simon that \(\mathbb {S}^3\) endowed with any Riemannian metric admits an embedded minimal sphere [26]. These are the only two scenarios where minimal spheres exist.

Theorem 3.9

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion. Then \(\mathbb {E}\) admits a immersed minimal sphere \(\Sigma \) if and only if M is topologically \(\mathbb {S}^2\).

Proof

Since spheres are simply connected, there exists an immersed minimal sphere \(\Sigma \subset \mathbb {E}\) if and only if there exists an immersed minimal sphere \(\widetilde{\Sigma }\subset \widetilde{\mathbb {E}}\), being \(\widetilde{\pi }:\widetilde{\mathbb {E}}\rightarrow \widetilde{M}\) the universal Riemannian covering given by Theorem 2.13. If the fibers of \(\widetilde{\pi }\) have infinite length, then any immersed minimal sphere \(\widetilde{\Sigma }\subset \widetilde{\mathbb {E}}\) must be an entire graph by the maximum principle. Thus, either \(\widetilde{\Sigma }\) is an entire graph or the fibers of \(\pi \) have finite length.

In the first case, \(\widetilde{\pi }\) induces an homeomorphism from \(\widetilde{\Sigma }\) to \(\widetilde{M}\), so \(\widetilde{M}\) is topologically \(\mathbb {S}^2\) and the submersion is trivial. As \(\widetilde{\Sigma }\) is an entire minimal graph and \(\widetilde{\mathbb {E}}\) is homeomorphic to \(\mathbb {S}^2\times \mathbb {R}\), we conclude that \(\widetilde{\Sigma }\) corresponds to one of the examples in Theorem 3.6. In the second case, \(\widetilde{\pi }\) is the Hopf fibration by Theorem 2.9 and the existence is a consequence of Simon’s result [26]. Hence, both cases lead to \(\widetilde{M}\) being topologically \(\mathbb {S}^2\). Theorem 2.13 implies that \(\mathbb {E}\) is a quotient of \(\widetilde{\mathbb {E}}\) under a group G of Killing isometries, and the only orientable quotient of a sphere is the sphere itself, so we deduce that G projects to the trivial group, i.e., G only consists of vertical translations. Therefore, \(\pi :\mathbb {E}\rightarrow M\) is again a Killing submersion over \(\mathbb {S}^2\), and the statement follows (see also Remark 2.10). \(\square \)

4 Compact stable surfaces with constant mean curvature

We begin by introducing the stability operator \(L=\Delta +|A|^2+\mathrm {Ric}(N)\) of a constant mean curvature surface immersed in the total space of a Killing submersion, where A is the shape operator of the immersion and N is a unit normal vector field. Next lemma gives another expression for L in terms of the scalar curvature.

Lemma 4.1

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion with bundle curvature \(\tau \) and Killing length \(\mu \). The stability operator of a constant mean curvature H surface \(\Sigma \) immersed in \(\mathbb {E}\) reads

$$\begin{aligned} L=\Delta -K+4H^2+S-\det (A), \end{aligned}$$

where K denotes the Gaussian curvature, respectively, and the scalar curvature S of \(\mathbb {E}\) (as a function on M) satisfies

$$\begin{aligned} \frac{1}{2}S=K_M-\tau ^2-\frac{\Delta \mu }{\mu }, \end{aligned}$$

being \(K_M\) the Gaussian curvature of M and \(\Delta \) its Laplace operator.

Proof

The computation is local so we shall consider the model metric on \(\mathbb {E}\) given by (2.3) and the metric in the base \(\lambda ^2(\mathrm {d}x^2+\mathrm {d}y^2)\). Using the expression of the Lie bracket (2.6) and Koszul formula, we derive the Levi-Civita connection in \(\mathbb {E}\) applied to the frame given by (2.2)

$$\begin{aligned} \begin{aligned} \overline{\nabla }_{E_1}E_1&=-\frac{\lambda _y}{\lambda ^2}E_2,&\overline{\nabla }_{E_1}E_2&=\frac{\lambda _y}{\lambda ^2}E_1+\tau E_3,&\overline{\nabla }_{E_1}E_3&=-\tau E_2,\\ \overline{\nabla }_{E_2}E_1&=\frac{\lambda _x}{\lambda ^2}E_2-\tau E_3,&\overline{\nabla }_{E_2}E_2&=-\frac{\lambda _x}{\lambda ^2}E_1,&\overline{\nabla }_{E_2}E_3&=\tau E_1,\\ \overline{\nabla }_{E_3}E_1&=-\tau E_2+\frac{\mu _x}{\lambda \mu }E_3,&\overline{\nabla }_{E_3}E_2&=\tau E_1+\frac{\mu _y}{\lambda \mu }E_3,&\overline{\nabla }_{E_3}E_3&=-\frac{1}{\mu }\overline{\nabla }\mu . \end{aligned} \end{aligned}$$
(4.1)

From (4.1), we get the following sectional curvatures of the basic planes spanned by \(\{E_1,E_2,E_3\}\) (R denotes the Riemannian curvature tensor):

$$\begin{aligned} \begin{aligned} R(E_1,E_2,E_2,E_1)&=K_M-3\tau ^2,\\ R(E_1,E_3,E_3,E_1)&=\tau ^2-\tfrac{1}{\mu }E_1(E_1(\mu ))-\tfrac{1}{\lambda \mu }E_2(\lambda )E_2(\mu ),\\ R(E_2,E_3,E_3,E_2)&=\tau ^2-\tfrac{1}{\mu }E_2(E_2(\mu ))-\tfrac{1}{\lambda \mu }E_1(\lambda )E_1(\mu ). \end{aligned} \end{aligned}$$
(4.2)

The scalar curvature S is twice the result of adding up the three quantities in (4.2), which leads to the expression in the statement. Finally, the given expression for L holds in any 3-manifold (see [19, Section 10]). \(\square \)

Note that the angle function of the immersion, defined as \(\nu =\langle N,\xi \rangle \), where \(\xi \) is the vertical Killing vector field, lies in the kernel of the Jacobi operator [3, Proposition 2.12]. Assuming that \(\nu \) is identically zero or never vanishes gives rise to two distinguished families of surfaces in \(\mathbb {E}\):

  • If \(\nu \equiv 0\), then \(\Sigma \) is everywhere vertical, so there exists a curve \(\Gamma \subset M\) such that \(\Sigma =\pi ^{-1}(\Gamma )\) and \(\Sigma \) is called the vertical cylinder over \(\Gamma \).

    Parametrizing \(\Gamma \) as \(\gamma :[a,b]\rightarrow M\) with unit-speed, and taking a horizontal lift T of \(\gamma '\) along \(\Sigma \), then \(\{T,E_3=\frac{1}{\mu }\xi \}\) is a orthonormal frame in \(\Sigma \). Using (4.1) we compute the second fundamental form in this frame as

    $$\begin{aligned} \sigma \equiv \left( \begin{matrix} \langle \overline{\nabla }_TT,N\rangle &{}\langle \overline{\nabla }_TE_3,N\rangle \\ \langle \overline{\nabla }_{E_3}T,N\rangle &{}\langle \overline{\nabla }_{E_3}E_3,N\rangle \end{matrix} \right) =\left( \begin{matrix}\kappa _g&{}\quad \tau \\ \tau &{}\quad -\langle \eta ,\frac{1}{\mu }\nabla \mu \rangle \end{matrix}\right) \end{aligned}$$

    where \(\kappa _g\) is the geodesic curvature of \(\gamma \) in M with respect to the unit conormal \(\eta =\pi _*N\) to \(\gamma \) in M, and \(\nabla \) denotes the gradient in M. Hence, \(\Sigma \) has mean curvature H if and only if \(\kappa _g=2H+\langle \eta ,\frac{1}{\mu }\nabla \mu \rangle \). Using the classical theory of ODEs, if M is complete, then for any point \(x\in M\) and any direction \(v\in T_xM\) there exists a unique constant speed curve \(\gamma :\mathbb {R}\rightarrow M\) with \(\gamma (0)=x\) and \(\gamma '(0)=v\), and such that \(\pi ^{-1}(\Gamma )\) has constant mean curvature H.

    Using the parametrization \((t,s)\mapsto \phi _t(\widetilde{\gamma }(s))\), where \(\{\phi _t\}\) is the group of vertical translations and \(\widetilde{\gamma }\) is a horizontal lift of \(\gamma \), the metric in \(\Sigma \) is expressed as \(\mu (\gamma (s))^2\mathrm {d}t^2+\mathrm {d}s^2\), which is not flat in general as it is in the unitary case (see [6, Proposition 2.10]).

  • If \(\nu \) has no zeroes, then \(\Sigma \) is everywhere transversal to the Killing vector field, and it is called a vertical multigraph. Note that \(\Sigma \) is a graph if and only if additionally \(\pi _{|\Sigma }:\Sigma \rightarrow M\) is injective.

    A vertical multigraph with constant mean curvature is always stable since \(\nu \) is a positive function such that \(L\nu =0\).

Theorem 4.2

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion with bundle curvature \(\tau \) and Killing length \(\mu \), and let \(\Sigma \) be a compact orientable stable surface with constant mean curvature immersed in \(\mathbb {E}\). Then one of the following assertions holds:

  1. (i)

    \(\Sigma \) is an entire minimal graph, M is compact and \(\int _M\frac{\tau }{\mu }=0\).

  2. (ii)

    \(\Sigma \) is a vertical cylinder over a closed curve in M and the fibers of \(\pi \) are compact. Hence, \(\Sigma \) is topologically a torus.

In particular, if both M and the fibers of \(\pi \) are not compact, then \(\mathbb {E}\) does not admit compact orientable stable surfaces with constant mean curvature.

Proof

Stability implies that the angle function \(\nu \), which is bounded on \(\Sigma \) and lies in the kernel of the stability operator, is identically zero or it never vanishes. If \(\nu \equiv 0\), then \(\Sigma =\pi ^{-1}(\Gamma )\) for a certain curve \(\Gamma \subseteq M\). Since \(\Sigma \) is compact, the fibers of \(\pi \) must have finite length and we get to (ii). If \(\nu \) never vanishes, then (i) follows from Lemma 3.4. \(\square \)

On the one hand, if M is compact and \(\int _M\frac{\tau }{\mu }=0\), then there are always surfaces satisfying item (i) by Theorem 3.6. On the other hand, compact vertical cylinders are topologically tori, but a characterization of their stability in terms of the corresponding curve in the base seems to be a tough task. It is important to mention that they may exist; for instance, if \(\Gamma \) is a curve with constant geodesic curvature \(\kappa _g\le 1\) in the hyperbolic plane \(\mathbb {H}^2\), then \(\Gamma \times \mathbb {S}^1\) is a stable torus with constant mean curvature \(H=\frac{1}{2}\kappa _g\) in the Riemannian product space \(\mathbb {H}^2\times \mathbb {S}^1\).

Corollary 4.3

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion, and let \(\Sigma \) be a compact stable surface with constant mean curvature immersed in \(\mathbb {E}\). If \(\Sigma \) is not topologically a torus, then it is an entire minimal graph.

Theorem 4.2 also allows us to give an alternative proof to the classification of compact stable surfaces with constant mean curvature in homogeneous 3-manifolds given by Meeks and Pérez [17, Theorem 4.18].

Corollary 4.4

Let X be a simply connected homogeneous Riemannian 3-manifold. If there exists a compact orientable stable surface \(\Sigma \) with constant mean curvature immersed in X, then X is isometric to \(\mathbb {S}^2\times \mathbb {R}\), being \(\Sigma \) a horizontal slice \(\mathbb {S}^2\times \{t_0\}\).

Proof

If X is homeomorphic to \(\mathbb {R}^3\), then Example 2.4 shows that X is the total space of a Killing submersion with non-compact base or fibers. If X is homeomorphic to \(\mathbb {S}^3\), then there exists a Hopf-like Killing submersion \(\pi :X\rightarrow \mathbb {S}^2\) for which \(\Sigma \) is not vertical by Example 2.4, so \(\Sigma \) should be an entire graph for \(\pi \), but the Hopf fibration does not admit entire sections. We conclude that X isometric \(\mathbb {S}^2\times \mathbb {R}\) (see [17, Theorem 2.4]), and \(\Sigma \) is an entire minimal graph by Theorem 4.2, unique up to vertical translations, so \(\Sigma \) must be a slice \(\mathbb {S}^2\times \{t_0\}\). \(\square \)

Finally, we will give a last characterization of our entire minimal graphs dropping the assumption of compactness. We will apply a result of Meeks, Pérez and Ros [19, Theorem 2.13] based on a previous result of Rosenberg [22].

Corollary 4.5

Let \(\pi :\mathbb {E}\rightarrow M\) be a Killing submersion and suppose that \(\Sigma \) is a complete stable surface with constant mean curvature H immersed in \(\mathbb {E}\).

  1. (1)

    If \(3H^2>\sup _{\pi (\Sigma )}\{\tau ^2-K_M+\frac{1}{\mu }\Delta u\}\), then \(\pi \) is topologically the projection \(\mathbb {S}^2\times \mathbb {R}\rightarrow \mathbb {S}^2\) and \(\Sigma \) is an entire minimal graph.

  2. (2)

    If \(3H^2=\sup _{\pi (\Sigma )}\{\tau ^2-K_M+\frac{1}{\mu }\Delta u\}\) and \(\mu \) is bounded, then \(\Sigma \) is either a vertical multigraph or a vertical cylinder.

Proof

The condition in item (1) is equivalent to \(3H^2+\frac{1}{2}S\ge c>0\) for some \(c>0\) by Lemma 4.1, so [19, Theorem 2.13] guarantees that \(\Sigma \) is homeomorphic to \(\mathbb {S}^2\) or to the real projective plane \(\mathbb {RP}^2\). Taking the universal covering of \(\Sigma \) if necessary, both cases lead to the existence of a constant mean curvature sphere \(\widetilde{\Sigma }\) immersed in \(\mathbb {E}\), which is also stable (see the comments below Assertion 2.2 in [19]). The result now follows from Corollary 4.3.

If the condition in item (2) holds, [19, Theorem 2.13] says that \(\Sigma \) has at most quadratic area growth, so it is in particular parabolic. Hence, [15, Corollary 1] implies that the angle function of \(\Sigma \), which is bounded because \(\mu \) is assumed bounded, is identically zero or never vanishes. This dicotomy reflects the two possible scenarios depicted in the statement. \(\square \)

The hypothesis in the statement of Corollary 4.5 is not expected to be optimal, as shown in [15, Corollary 3] for \(\mathbb {E}(\kappa ,\tau )\)-spaces. Also \(\sup _{\pi (\Sigma )}\{\tau ^2-K_M+\frac{1}{\mu }\Delta u\}\) might be infinite, so in that case no information about the surface is obtained.