Polyharmonic hypersurfaces into pseudo-Riemannian space forms

In this paper, we shall assume that the ambient manifold is a pseudo-Riemannian space form \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N^{m+1}_t(c)$$\end{document}Ntm+1(c) of dimension \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m+1$$\end{document}m+1 and index t (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m\ge 2$$\end{document}m≥2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1 \le t\le m$$\end{document}1≤t≤m). We shall study hypersurfaces \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}^{m}_{t'}$$\end{document}Mt′m which are polyharmonic of order r (briefly, r-harmonic), where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r\ge 3$$\end{document}r≥3 and either \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t'=t$$\end{document}t′=t or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t'=t-1$$\end{document}t′=t-1. Let A denote the shape operator of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}^{m}_{t'}$$\end{document}Mt′m. Under the assumptions that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}^{m}_{t'}$$\end{document}Mt′m is CMC and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{Tr}}}\,{A}^{2}$$\end{document}TrA2 is a constant, we shall obtain the general condition which determines that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}^{m}_{t'}$$\end{document}Mt′m is r-harmonic. As a first application, we shall deduce the existence of several new families of proper r-harmonic hypersurfaces with diagonalizable shape operator, and we shall also obtain some results in the direction that our examples are the only possible ones provided that certain assumptions on the principal curvatures hold. Next, we focus on the study of isoparametric hypersurfaces whose shape operator is non-diagonalizable and also in this context we shall prove the existence of some new examples of proper r-harmonic hypersurfaces (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r \ge 3$$\end{document}r≥3). Finally, we shall obtain the complete classification of proper r-harmonic isoparametric pseudo-Riemannian surfaces into a three-dimensional Lorentz space form.


Introduction
We recall that, in the Riemannian case, harmonic maps are the critical points of the energy functional where ∶ M → N is a smooth map between two Riemannian manifolds (M m , g) and (N n , h) . In particular, is harmonic if and only if it is a solution of the Euler-Lagrange system of equations associated with (1.1), i.e., The left member of (1.2) is a vector field along the map or, equivalently, a section of the pullback bundle −1 TN : it is called tension field and denoted ( ) . In addition, we recall that, if is an isometric immersion, then is a harmonic map if and only if the immersion defines a minimal submanifold of N (see [12,13] for background). Next, in order to define the notion of an r-harmonic map, we consider the following family of functionals which represent a version of order r of the classical energy (1.1). If r = 2s , s ≥ 1: In the case that r = 2s + 1: Here, Δ = d * d represents the Laplacian on the pullback bundle −1 TN . Then a map ∶ (M m , g) → (N n , h) is r-harmonic if, for all variations t , In the case that r = 2 , the functional (1.3) is called bienergy and its critical points are the so-called biharmonic maps. A very ample literature on biharmonic maps is available and we refer to [9,16,31,32] for an introduction to this topic. More generally, the r-energy functionals E r ( ) defined in (1.3), (1.4) have been intensively studied (see [4-6, 20-22, 26-28, 37, 38], for instance).
We say that an r-harmonic map is proper if it is not harmonic (similarly, an r-harmonic submanifold, i.e., an r-harmonic isometric immersion, is proper if it is not minimal). We point out that, as observed by Maeta in his series of papers [20][21][22], in general r-harmonic does not imply r ′ -harmonic for r ′ > r unless the target manifold is flat. (1.4) In our recent work [26], we proved some general results for r-harmonic hypersurfaces into space forms and deduced that the value of the integer r plays a crucial role to generate geometric phenomena which differ substantially from the classical situation corresponding to the biharmonic and triharmonic cases. For instance, if ≥ 3 , there exists no isoparametric hypersurface of m+1 of degree which is proper biharmonic or triharmonic. By contrast, when r ≥ 5 , there are several examples of such hypersurfaces which are proper r-harmonic (see [26]). From the point of view of the differential geometry of submanifolds, the difficulties which one encounters in studying the equations which define a general r-harmonic submanifold are huge. Therefore, a reasonable starting point is to focus on the case that the ambient is a space form N m+1 (c) (here and below, c denotes the sectional curvature) and study CMC hypersurfaces with constant squared norm ‖A‖ 2 of the shape operator. In this order of ideas, in the Riemannian case we obtained the following general result: Theorem 1.1 [26] Let M m be a non-minimal CMC hypersurface in a Riemannian space form N m+1 (c) and assume that ‖A‖ The first goal of this paper is to establish a version of Theorem 1.1 when the ambient is a pseudo-Riemannian space form. Indeed, we shall prove: where in the last equality of (1.8) we used the fact that the shape operator A is self-adjoint.
To the purpose of a first quick comparison between Theorems 1.1 and 3.5, we observe that, if the induced metric on the hypersurface is Riemannian, then Tr A 2 = ‖A‖ 2 . The quantity Tr A 2 already appeared in [35] in the context of the study of proper biharmonic surfaces in a pseudo-Riemannian three-dimensional space. We also point out that in [11,19], with a slight abuse of notation, the right-hand side of (1.8) is denoted by ‖A‖ 2 instead of Tr A 2 .
In general, A is not diagonalizable. Indeed, since ⟨ , ⟩ is not positive definite, it is wellknown that complex eigenvalues may appear or, in some cases, A may not even be diagonalizable over ℂ (see [15]). We shall study these situations in detail in Sect. 3. In order to have a better understanding of the results presented in Sect. 3, it is important to observe that, differently from the Riemannian case, in the pseudo-Riemannian context it may happen that Tr A 2 ≤ 0 and, in particular, Tr A 2 = 0 does not imply that A = 0 , a fact which will lead us to show that the family of proper r-harmonic hypersurfaces ( r ≥ 3 ) of Theorem 3.5 with Tr A 2 ≤ 0 is not empty (see Theorem 3.18).
Next, we shall describe several new examples and some geometric applications. In order to state in detail our results, it is convenient to introduce first some basic notions concerning pseudo-Riemannian manifolds and geometry. Therefore, our paper is organized as follows.
In Sect. 2, we shall review some basic aspects of the theory of pseudo-Riemannian space forms and pseudo-Riemannian geometry. Moreover, we shall describe how to generalize the notion of r-harmonicity in this context. In Sect. 3, we shall state our main results and their geometric applications. Finally, in Sect. 4, we shall provide all the proofs. For the sake of completeness, we mention in this introduction that another possible, interesting definition of an r-order version of the energy functional, which was proposed by Eells-Sampson and Eells-Lemaire (see [12,14]), is As for a detailed discussion and comparison between definitions (1.3), (1.4) and (1.9), we refer to [6]. We believe that most of the techniques used in [6] could also be applied in the pseudo-Riemannian context, but we shall not pursue this option in this paper.

Pseudo-Riemannian geometry, pseudo-Riemannian space forms and r-harmonicity
A basic reference for pseudo-Riemannian geometry is the classical book of O'Neill (see [30]), but for the specific topics treated in this section we also refer to [1,3,11,19,35,39]. Let (M m t , g) be a pseudo-Riemannian manifold of dimension m with a non-degenerate metric of index t ( 0 ≤ t ≤ m ). In order to clarify the notion of index t, let us first recall that non-degeneracy means that the only vector X ∈ T p M satisfying g p (X, Y) = 0 for all Y ∈ T p M is X = 0 , for any p ∈ M . A local orthonormal frame field of (M m t , g) is a set of local vector fields {e i } m i=1 such that g(e i , e j ) = i ij , with 1 = … t = −1 , t+1 = … m = 1. Next, let us fix terminology and notations concerning pseudo-Riemannian space forms. The m-dimensional pseudo-Euclidean space with index t is denoted by The m-dimensional pseudo-Riemannian sphere denoted by m t (c) is defined as follows: The m-dimensional pseudo-Riemannian hyperbolic space denoted by ℍ m t (c) is defined as follows: , with the induced metric from ℝ m+1 t+1 , is a complete pseudo-Riemannian manifold with index t and constant negative sectional curvature c.
A pseudo-Riemannian space form refers to one of the three spaces ℝ m t , m t (c), ℍ m t (c) . We shall write m t and ℍ m t for m t (1) and ℍ m t (−1) , respectively. Sometimes, to provide a unified treatment, we also use N m t (c) to denote a pseudo-Riemannian space form of sectional curvature c.
The flat ( c = 0 ) pseudo-Riemannian space ℝ m t is called Minkowski space, while m t (c) and ℍ m t (c) are known as de Sitter space and anti-de Sitter space, respectively. When the index is t = 1 , these spaces are also referred to as Lorentz space forms. We also point out In particular, m m−1 (c) and ℍ m 1 (c) are not simply connected. Finally, we point out that in this paper we restrict our study to connected manifolds. Therefore, even if in some definitions formally one has two connected components (for instance, according to (2.1), (2.2), this happens to m m (c) and ℍ m 0 (c) ), we shall implicitly assume that we just work with one connected component.
Next, we recall for future use that the sectional curvature tensor field of N m t (c) is described by the following simple expression: Pseudo-Riemannian space forms have important applications in the theory of general relativity and ℝ 3 1 , 3 1 (c) and ℍ 3 1 (c) are model spaces for Minkowski, de Sitter and anti-de Sitter space-time, respectively. Now, in order to make this paper as self-contained as possible, we follow [11] and recall here how the basic operators of Riemannian geometry extend to the pseudo-Riemannian context. In a chart of (M m t , g) , for a local orthonormal frame field we have the following basic identities, which hold for any X ∈ C(TM) , f ∈ C ∞ (M) and bilinear form b ∈ C(⊙ 2 TM): In this paper, we shall focus on the study of pseudo-Riemannian (or, non-degenerate) hypersurfaces of a pseudo-Riemannian space form. We can describe such hypersurfaces by means of a smooth ∶ M m t � → N m+1 t and the hypothesis that M m t ′ is pseudo-Riemannian simply means that the metric induced by is non-degenerate. Therefore, we can assume that, locally, there always exists an orthonormal frame field {e i } m i=1 on M m t ′ . Moreover, denoting a unit normal vector field, {e i , } is a local orthonormal frame field on N m+1 t and we have two possibilities: either ⟨ , ⟩ = = 1 and t � = t , or = −1 and t � = t − 1.
In particular, when the ambient space is Lorentzian and = −1 , then the hypersurface is Riemannian and the standard terminology is to say that it is space-like. Let us denote B the second fundamental form and A = A the associated shape operator. The classical Gauss and Weingarten formulas hold as in the Riemannian case, i.e., for all tangent vector fields to M m t ′ . We observe that it is easy to deduce from (2.5) that The mean curvature vector field of M m t ′ , denoted , is defined by where, using (2.6), we deduce that the mean curvature function f is given by In this paper, we shall always say that a hypersurface of a pseudo-Riemannian space is minimal if f vanishes identically. However, we mention that, when the ambient space is Lorentzian and the hypersurface is space-like, the term maximal instead of minimal is also used in the literature, because in this case small regions are local maximizers of the volume functional.
Next, we introduce an important family of pseudo-Riemannian hypersurfaces of N = m+1 t ( m ≥ 2 , 1 ≤ t ≤ m ) and give, in Table 1, the expression for A and (see [1]). The definition of the hypersurfaces in Table 1 are given in Tables 2 and 3.
The importance of the families of hypersurfaces described in Table 1 lays in the fact that essentially any pseudo-Riemannian hypersurface of m+1 t with diagonalizable shape operator having at most two distinct constant principal curvatures is locally congruent to one of these (see Theorem 5.1 of [1], and also Theorem 3.13, for more details).

Remark 2.1
A similar family of hypersurfaces of pseudo-hyperbolic spaces is available, but we omit its description because all the properties of r-harmonic submanifolds in ℍ m+1 t * can be deduced from those of r-harmonic submanifolds in m+1 t , with t = m + 1 − t * . This is a consequence of the fact that, up to a multiplicative constant factor -1 in the metric (that determines a non-isometric transformation), irrelevant for the r-harmonic equation, we can identify m+1 t with ℍ m+1 m+1−t . For this reason, without loss of generality, in this paper all the existence and classification results for r-harmonic submanifolds that do not depend on a given signature of the metric will be stated only for the cases that the curvature of the ambient space is either c = 1 or c = 0.
Next, to prepare the ground for the study of r-harmonicity, we consider a general smooth map ∶ (M m t � , g) → (N n t , h) between two pseudo-Riemannian manifolds. Let us denote ∇  Table 2 Definition of hypersurfaces in Table 1 m Table 3 Definition of hypersurfaces in Table 1 k the induced connection on the pullback bundle −1 TN . In the pseudo-Riemannian context, the operator corresponding to the classical Riemannian rough Laplacian on sections of −1 TN , which will still be denoted Δ , becomes We are now in the right position to summarize the key points which enable us to describe the generalization of the notions of harmonicity and, more generally, r-harmonicity to the pseudo-Riemannian context.
h) between two pseudo-Riemannian manifolds is harmonic if its tension field vanishes identically, i.e., is a local orthonormal frame field on M m t ′ as above. As for papers and examples in this context, we cite [18,34].
Similarly, taking into account (2.8), we can define the r-energy for a map between two pseudo-Riemannian manifolds precisely as in (1.3), (1.4). In the Riemannian case, the explicit expression for the r-tension field associated with E r ( ) was obtained by Maeta and Wang (see [20,37]). Essentially, in the pseudo-Riemannian context, the only relevant difference appears when one has to take a trace (for instance, see the last equation in (2.4)). More specifically, we have the following expression for the r-tension field, where Δ and are given in (2.8) and (2.9), respectively: is a local orthonormal frame field tangent to M m t ′ (the sum over i is not written but understood). Similarly, From the analytic point of view, one of the major differences with respect to the Riemannian case is the fact the PDE's system r ( ) = 0 is not elliptic when 1 ≤ t � ≤ m − 1.

Statement of the results
In the Riemannian case, when the ambient space form has nonpositive sectional curvature there are several results which assert that, under suitable conditions, an r-harmonic submanifold is minimal (see [9,20,23] and [28], for instance). Things drastically change when the ambient is the Euclidean sphere m+1 . Indeed, in this case several examples of proper r-harmonic hypersurfaces have been constructed and studied (see [6,22,26,27] and references therein).
In the pseudo-Riemannian setting, there are in the literature several interesting results, but most of them are limited to the biharmonic case. One of the instances which have attracted more attention is the study of space-like biharmonic hypersurfaces in a Lorentzian space form. In this case, it seems that the sign of the curvature of the ambient produces phenomena which are, in some sense, dual with respect to the Riemannian case. In other words, positive curvature increases the rigidity of space-like biharmonic hypersurfaces. For example, the following interesting result was proved by Ouyang: In this order of ideas, we also have: . Then M must be totally geodesic, i.e., ℝ 2 or 2 .
This trend is confirmed in the r-harmonic case ( r ≥ 3 ), as we shall show in Corollary 3.8.
By contrast, when the index of the hypersurface is positive, some results in the biharmonic case are available, but their interpretation is less evident. For the purpose of comparison with the results of this paper, we report here, using our notations, the following interesting result of Liu and Du (see Theorem 1.2 of [19], where the cases c ≠ 0, 1 are also dealt with explicitly): Theorem 3.3 [19] Let M m t ′ be a pseudo-Riemannian proper biharmonic hypersurface in N m+1 t (c) . If M m t ′ has diagonalizable shape operator with at most two distinct principal curvatures, then c ≠ 0 . Furthermore, when c = 1 , then t � = t and M m t ′ is congruent to either m t (2) or

Remark 3.4
When c = 1 , our statement of Theorem 3.3 is equivalent to Theorem 1.2 of [19], but our formulation makes it easier the comparison with the classical results for biharmonic hypersurfaces in the Riemannian case.
The results of [19] in the biharmonic case were recently refined in [11]. Finally, for the sake of completeness, we also have to cite the interesting paper [35] by Sasahara, where the author classifies proper biharmonic curves and surfaces in de Sitter 3-space and anti-de Sitter 3-space. In Theorem 3.20, we shall extend this classification for surfaces in N 3 1 (c) to the case r ≥ 3.
Our investigation of this type of problems in the pseudo-Riemannian context starts with the following general result, which is a pseudo-Riemannian version of Theorem 1.1. Remark 3. 6 The special case r = 2 in Theorem 3.5 is part of the result proved in [19], where the condition of biharmonicity for a general hypersurface M m t � ↪ N m+1 t (c) was computed. In the case of surfaces, (1.5) was first obtained in [35].

Remark 3.7
It is important to point out that in Theorem 3.5 we do not require that the shape operator A be diagonalizable. By way of example, in Theorem 3.18 we shall exhibit a new family of r-harmonic surfaces in ℍ 3 1 whose shape operator is not diagonalizable.
A first, immediate consequence of Theorem 3.5 is the following In the same spirit:

Remark 3.12
The third-order polynomial P 3 (c) in (3.4) is equivalent to the one which appears in the Riemannian case. To see this, it is enough to set Then, up to a constant factor −1∕x 3 , P 3 (c) becomes the third-order polynomial obtained in Theorem 1.2 of [27]. We point out that, in the special case that k = m − k , we have For a more detailed discussion on the existence and qualitative behaviour of admissible roots of P 3 (c) in the general case, we refer to [27].
Next, we show that the examples given in Corollary 3.10 and Theorem 3.11 are the only possible ones within a certain class of hypersurfaces. More precisely, we prove the following result: Remark 3.14 Of course, it is straightforward to state explicitly the version of Corollary 3.10, Theorems 3.11 and 3.13 in the case that c = −1 . However, for the reasons explained in Remark 2.1, we omit the details.
Next, we obtain some geometric results for triharmonic surfaces which provide a version in the pseudo-Riemannian case of some facts which we proved in [26] in the Riemannian case.

Theorem 3.15 Assume r ≥ 3 . Let M 2 t ′ be a pseudo-Riemannian triharmonic CMC surface in N 3 t (c) and assume that its shape operator
Next, we focus on the case that the ambient space is 3 t . Our result is:

Theorem 3.16 Let M 2 t ′ be a CMC proper triharmonic surface in 3 t and assume that its shape operator A is diagonalizable. Then M 2 t ′ is an open part of the small pseudo-hypersphere 2 t (3).
A widely studied family of hypersurfaces in the pseudo-Riemannian setting is that of isoparametric Lorentzian hypersurfaces of a Lorentzian space form. For this specific topic, we refer to [25,36] and, for further background, to [15,29]. We recall that a Moreover, according to [30,Chapter 9] there exist bases where the shape operator A assumes one of the following Jordan canonical forms: Here b 0 is assumed to be non-zero. In cases I, II and III, the eigenvalues are real, while a 0 ± ib 0 are complex eigenvalues in case IV. A Lorentzian isoparametric hypersurface in N m+1 1 (c) is called of type I, II, III or IV according to the form of its shape operator A.
We observe that a Lorentzian isoparametric hypersurface in N m+1 1 (c) is a CMC hypersurface with Tr A 2 constant. Therefore, Theorem 3.5 applies to this type of hypersurfaces.
Now, a direct computation shows that for a Lorentzian isoparametric hypersurface of type I, II and III, with Tr A ≠ 0 , we have Tr A 2 > 0 . Differently, for hypersurfaces of type IV it is possible to have Tr A ≠ 0 and Tr A 2 ≤ 0 . If the curvature of the ambient space is c = 0 , we easily deduce from Theorem 3.5 that a non-minimal Lorentzian isoparametric hypersurface is r-harmonic ( r ≥ 2 ) if and only if Tr A 2 = 0 . These observations prove our first result in this context, that is In [25, Theorem 4.10], Magid gave a proof that there exist no isoparametric Lorentzian hypersurfaces of type IV in the flat Lorentzian space form ℝ m+1 1 . However, it was pointed out in [10] that there could be some gaps in the arguments of [25].
By contrast, in the case that the curvature of the ambient space is c = −1 , we shall construct some new examples of proper r-harmonic Lorentzian isoparametric surfaces in ℍ 3 1 with non-diagonalizable shape operator of type IV. More precisely, first, as above, we observe that Theorem 3.5 implies that a Lorentzian isoparametric hypersurface in ℍ m+1 1 can be r-harmonic only if its shape operator is non-diagonalizable and of type IV. The following result proves that complex circles (see [24]) provide a geometrically significant example of such hypersurfaces: where is a real constant and k(s) ≠ 0 (see [2,15,35] for details).

Remark 3.21
Comparison of [35,Theorem 5.4] with our Theorem 3.20 shows that the relevant differences between the biharmonic case and the case r ≥ 3 are the appearance in Theorem 3.20 of the family of solutions of type (6) for all r ≥ 3 , and of type (3) and (4) for r ≥ 5 . These facts, together with Theorem 3.18, support the following general idea: We can find geometrically interesting situations where there is no biharmonic instance, but there exist examples of r-harmonic immersions, r ≥ 3 . Moreover, comparing the result of the present paper with those of [26], we see that the notion of r-harmonicity is more flexible in the pseudo-Riemannian setup compared to the Riemannian case as it allows for a larger class of solutions. In this order of ideas, we also cite the phenomenon illustrated in Remark A.1.

Example 3.22
First, we point out that the family of r-harmonic surfaces in 3 1 obtained in Case (5) of Theorem 3.20 is very ample. This is a consequence of the fact that, given any real constant and smooth function k(s), it is always possible to determine (at least locally) a null curve such that its associated B-scroll (3.6) verifies (3.7). As proved in [17], the existence of such a null curve can be deduced by solving a suitable Cauchy problem for a first-order linear system of ordinary differential equations. In the proof of Theorem 3.20 we shall show that K = 2 + 1 and any such surface is proper r-harmonic provided that 2 = r − 1 (r ≥ 2). As a special case, following the procedure described in the Appendix of [17], here we give the explicit expression of the null curve (s) and of its associated vector field B(s) assuming that k(s) ≡ 1 . Note that when, as in our example, the function k(s) is bounded, (s) is defined on the whole ℝ . To this end, let where Then (s) is a null curve in 3 1 with associated vector field and the parametrized surface defines a B-scroll with Gauss curvature K = 1 + 2 and k(s) ≡ 1 . Choosing = ± √ r − 1, we obtain the desired r-harmonic surface in 3 1 .
For the sake of completeness, the method of computation used to determine the curve (s) and the vector field B(s) will be given at the end of this paper, in "Appendix A".

Proofs
As a preliminary step, in the following lemma we state without proof some standard facts which we shall use in this section.  Next, we perform our first computation:

c) be a pseudo-Riemannian hypersurface and denote by the unit normal vector field. Then
Proof We work with a geodesic frame field X i m i=1 around an arbitrarily fixed point p ∈ M m t � . Also, we simplify the notation writing ∇ for ∇ M . Since = f , around p we have: Then at p we have: where, for the second equality, we also used (2.6). Now we take the sum over i as in (2.8) and, using (2.4) and Lemma 4.1, we obtain (4.1). (Note that the sign convention for Δ and Δ is given in (2.4), (2.8).) ◻ Next, we assume that the mean curvature function f is constant and we obtain

c) be a pseudo-Riemannian hypersurface and assume that its mean curvature function f is equal to a constant . Then
Proof Since f is constant, according to Lemma 4.2 we have Δ = Tr A 2 and so Now, around p: At p: Next, computing as in Lemma 4.2, we find and, since f is constant, the proof ends immediately. ◻ We are now in the right position to prove our first theorem.

Proof of Theorem 3.5 As Tr
Proof of Corollary 3. 10 We use (2.7), (1.8) and Table 1 to compute Tr A 2 = m (c − 1) and = √ c − 1 . Then the thesis follows by direct substitution in (1.7), where the curvature of the ambient space is 1 and also = 1 . ◻

Proof of Theorem 3.11
The hypotheses of Theorem 3.5 are verified and so the condition for r-harmonicity is equation (1.7) where, as in the proof of Corollary 3.10, the curvature of the ambient space is 1 and = 1 . Then we compute using the explicit expression for the shape operator given in Table 1: Next, after direct substitution and simplification, we find that equation (1.7) is equivalent to where P 3 (c) is the third-order polynomial defined in (3.4). Now, since ck − m = 0 corresponds to = 0 , the conclusion of the theorem follows immediately.  Table 1. Then a case-by-case direct inspection of (1.7), using again (2.7), (1.8) and Table 1 to compute Tr A 2 and , shows that the only r-harmonic hypersurfaces in this family are those given in Corollary 3.10 and Theorem 3.11. In the case that c = 0 , Theorem 5.1 of [1] says that M m t ′ is one of the hypersurfaces listed in (R-1)-(R-6), p. 131 of [1]. Then, again, the thesis follows easily by direct inspection. ◻

Proof of Theorem 3.15
The 3-tension field is described by (2.11) with s = 1 . In the first part of the proof, for future reference, we do not make any assumption on the dimension m and the curvature c. We observe that ( ) = m and use Lemma 4.2 with f constant and (2.3).
We have: Similarly, we compute Using (4.2), (4.6) and (4.7) into (2.11) and Lemma 4.3, we obtain the explicit expression of the 3-tension field: Therefore, we conclude that M m t ′ is a triharmonic hypersurface in N m+1 t (c) if and only if either it is minimal or From now on, we assume that M m t is not minimal and we use the hypothesis c ≤ 0 . First, we analyse the case c = 0 . If Tr A 2 is constant, then it follows immediately from (4.8) that Tr A 2 2 = 0 and so, since A is diagonalizable, the hypersurface is totally geodesic, a contradiction. If Tr A 2 is not a constant, then there exists an open set U of the surface such that grad Tr A 2 ≠ 0 on U. We deduce from (4.8)(ii) that 0 is an eigenvalue of A on U. Now we use the assumption that m = 2 . Since M 2 t ′ is CMC, it is easy to conclude that necessarily Tr A 2 is again a constant on U, and this is a contradiction. Next, if c ≠ 0 , by the assumption c ≤ 0 we deduce that c and have opposite sign. Then, using this observation in (4.8) (i), the proof follows essentially the same argument as in the case c = 0 . ◻ Proof of Theorem 3. 16 First, from an argument similar to the proof of Theorem 3.15, we deduce that Tr A 2 must be a constant on M 2 t ′ and also the two principal curvatures are constant.
Then, according to Theorem 3.13, we deduce that the only possibility is that M 2 t ′ is an open part of a small pseudo-hypersphere 2 t (3) because, according to Theorem 3.11, in these dimensions there exists no generalized pseudo-Clifford torus which is proper triharmonic. Indeed, using k = 1, r = 3 and m = 2 , (3.4) becomes Now, the only real root of P 3 (c) is c = 2 and it corresponds to the minimal pseudo-Clifford torus. ◻ Proof of Theorem 3. 18 We have to find under what conditions on a, b the complex circle x(s, t) is r-harmonic. By using standard techniques of the theory of surfaces in ℍ 3 1 ⊂ ℝ 4 2 , we can compute the shape operator of the complex circle x(s, t) and obtain The condition ab ≠ 0 ensures that Tr A ≠ 0 , that is M 2 1 is not minimal. (4.8) The case r = 2 . From Theorem 3.5, we deduce that M 2 1 is biharmonic if and only if Tr A 2 + 2 = 0 . Now, a direct computation gives from which, since b 2 − 1 = a 2 , we deduce that M 2 1 cannot be proper biharmonic. The case r > 2 . In this case, using Theorem 3.5, we have two possibilities, that is: (i) either (ii) or 4a 2 b 2 − 1 ≠ 0 and, according to (1.7), Case (i), together with the hypothesis b 2 − a 2 = 1 , gives exactly point (1) of the theorem. As for Case (ii), we have a proper solution if and only if from which we deduce that the only possible value is r = 3 . Taking into account that, by assumption, b 2 − a 2 = 1 , it is easy to conclude Case (2) of the theorem. ◻ Proof of Theorem 3.20 Case (i). Here we assume that the shape operator A is diagonalizable. Since the hypersurface has dimension m = 2 and constant principal curvatures, we can apply Theorem 3.13. Moreover, taking into account Remark 2.1, the proper r-harmonic surfaces in ℍ 3 1 can be deduced from those of 3 2 . Putting these facts together and using (3.5), we obtain Cases (1) − (4) in the statement of the theorem.
Case (ii). Now, we assume that the shape operator A is non-diagonalizable. In this case, t � = 1 and so = 1 . If A has a double, real eigenvalue , then the argument of [2, Proposition 4.1] enables us to conclude that M 2 1 is a B-scroll in N 3 1 (c) over a null curve , and its shape operator A, with respect to the coordinate frame field { ∕ s, ∕ u} associated with (3.6), is given by with k(s) ≠ 0 . Therefore, Tr A 2 = 2 2 and = . Next, we assume that the B-scroll is not minimal, i.e., ≠ 0 , and we apply (1.7): If c = −1, m = 2, = 1, then there are no solutions; if c = 1, m = 2, = 1 , then it is easy to deduce that the B-scroll is proper r-harmonic if and only if Finally, observing that the Gauss curvature K of the B-scroll is given by K = det(A) + c , with c = 1 , we conclude Case (5).
It remains to analyse the case that the shape operator A has two complex eigenvalues. This will lead us to Case (6). Indeed, from [2, p. 453], we conclude that M 2 1 is a flat Lorentzian surface in ℍ 3 1 with parallel second fundamental form in ℝ 4 2 . Then, according to [24], M 2 1 is locally congruent to a complex circle, and therefore, the conclusion follows from Theorem 3.18. ◻

Appendix A
In this appendix, we describe how to construct the B-scrolls in 3 1 ⊂ ℝ 4 1 described in Example 3.22. We follow exactly the method given in Appendix of [17].
First, we introduce the following matrices: Next, we introduce a 4 × 4 matrix X(s): where the first 3 columns will be the Cartan frame  1 Finally, we point out that the method of this Appendix does not require k(s) ≠ 0 . If we assume that the smooth function k(s) vanishes at an isolated point, say k(s 0 ) = 0 , we have a very interesting phenomenon. The associated B-scroll is still a proper r-harmonic surface in 3 1 provided that 2 = r − 1 . Indeed, on the whole surface we have Tr A = 2 and Tr A 2 = 2 2 and so equation (1.7) applies and r-harmonicity follows exactly as in the proof of Theorem 3.20. But, in this case, the minimal polynomial of A is (x − ) at all points x(s 0 , u) , while it is (x − ) 2 elsewhere. Therefore, these proper r-harmonic surfaces in 3 1 are not isoparametric. By contrast, in the Riemannian case, all the known proper r-harmonic CMC hypersurfaces with Tr A 2 constant are isoparametric (see [26]).