Proper r-harmonic functions from Riemannian manifolds

We introduce a new method for constructing complex-valued r-harmonic functions on Riemannian manifolds. We then apply this for the important semisimple Lie groups SO(n)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SO }(n)$$\end{document}, SU(n)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SU }(n)$$\end{document}, Sp(n)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{Sp }(n)$$\end{document}, SLn(R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SL }_{n}({\mathbb {R}})$$\end{document}, Sp(n,R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{Sp }(n,{\mathbb {R}})$$\end{document}, SU(p,q)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SU }(p,q)$$\end{document}, SO(p,q)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SO }(p,q)$$\end{document}, Sp(p,q)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{Sp }(p,q)$$\end{document}, SO∗(2n)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SO }^*(2n)$$\end{document} and SU∗(2n)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{SU }^*(2n)$$\end{document}.


Introduction
Biharmonic functions are important in physics and arise in areas of continuum mechanics, including elasticity theory and the solution of Stokes flows. Here the literature is vast, but usually the domains are either surfaces or open subsets of flat Euclidean space, with only very few exceptions, see for example [1]. The development of the very last years has changed this and can be traced, e.g., in the following publications [2][3][4][5]8]. There the authors develop methods for constructing explicit r -harmonic functions on classical Lie groups and even some symmetric spaces.
In this paper, we introduce a general scheme for constructing complex-valued r -harmonic functions φ : (M, g) → C on Riemannian manifolds. We then show that this can be employed in the important cases of the classical compact simple Lie groups SO(n), SU(n), Sp(n), and furthermore in the non-compact cases of SL n (R), Sp(n, R), SU( p, q), SO( p, q), Sp( p, q), SO * (2n) and SU * (2n).

Preliminaries
Let (M, g) be a smooth manifold equipped with a Riemannian metric g. We complexify the tangent bundle T M of M to T C M and extend the metric g to a complex-bilinear form on Then the gradient ∇φ of a complex-valued function φ : (M, g) → C is a section of T C M. In this situation, the well-known linear Laplace-Beltrami operator (alt. tension field) τ on (M, g) acts locally on φ as follows: For two complex-valued functions φ, ψ : (M, g) → C, we have the following well-known relation where the conformality operator κ is given by For a positive integer r , the iterated Laplace-Beltrami operator τ r is defined by

Definition 2.1
For a positive integer r , we say that a complex-valued function φ : It should be noted that the harmonic functions are exactly the 1-harmonic and the biharmonic functions are the 2-harmonic ones. In some texts, the r -harmonic functions are also called polyharmonic of order r . We also note that, if a function is r -harmonic, then it is also p-harmonic for any p ≥ r . Hence, one is usually interested in studying functions which are proper r -harmonic.

Definition 2.3
Let (M, g) be a Riemannian manifold. Then a complex-valued function φ : M → C is said to be an eigenfunction if it is eigen both with respect to the Laplace-Beltrami operator τ and the conformality operator κ, i.e., there exist complex numbers λ, μ ∈ C such that

Remark 2.4
It is clear that any element φ ∈ E of an eigenfamily is an eigenfunction in the sense of Definition 2.3.

A new construction method
Let φ : (M, g) → C be a complex-valued function and f : U → C be a holomorphic function defined on an open subset U of C containing the image φ(M) of φ. Then it is a direct consequence of the chain rule that the composition In general, this differential equation seems difficult, if not impossible, to solve. But as we will now see, this problem can actually be solved in the case when φ : (M, g) → C is an eigenfunction in the sense of Definition 2.3. In that case, equation (3.1) takes the following form Thus, will be harmonic if and only if f satisfies the complex ordinary differential equation Having constructed a harmonic function, we can now use it to construct a biharmonic function.
Then we can make 2 proper biharmonic by requiring that it solves the Poisson equation Employing the chain rule, we see that this equation is equivalent to i.e., we want f 2 to satisfy the complex ordinary differential equation which can again be solved using elementary methods. This procedure can be continued inductively to obtain proper r -harmonic functions. Indeed, suppose that a holomorphic function f r −1 : C\(−∞, 0] → C has been constructed so that r −1 = f r −1 • φ is proper (r − 1)-harmonic. Then we can consider the function r = f r • φ : W → C, where f r : C\(−∞, 0] → C is holomorphic, and study the Poisson equation which, as before, reduces to a complex ordinary differential equation that can easily be solved. To state the solutions furnished by the method that we have just described, we will make extensive use of the holomorphic principal logarithm log : C\(−∞, 0] → C and the standard notation z α = exp(α log(z)).
As already hinted, our proper r -harmonic functions will be defined on the open subset W of M given by We wish this domain to be as large as possible, so it might be convenient to introduce branch cuts other than ( for some λ, μ ∈ C. Then for a natural number r ≥ 1 and (c 1 , c 2 ) ∈ C 2 , any non-vanishing function

Eigenfunctions on semisimple Lie groups
In the last section, we have introduced a new method for constructing r -harmonic functions on Riemannian manifolds in terms of eigenfunctions. This method is not worth much without the existence of such objects. The goal of this section is to show that these certainly do exist in many of the important cases of the classical semisimple Lie groups. As already pointed out in Remark 2.4, the elements φ ∈ E of a given eigenfamily are automatically eigenfunctions i.e. exactly the objects that we are seeking. Eigenfamilies on classical semisimple Lie groups have already been constructed in [5,6] and [7]. We use the results from these papers to present interesting examples, see Table 1. Let us discuss our notation and recall the definitions of some of the classical semisimple Lie groups. We will denote the n × n identity matrix by I n , and we will also use the standard notations The Lie groups in this section are all assumed to be equipped with the left-invariant Riemannian metric induced by the canonical inner product on the Lie algebra gl n (C). We will consider the well-known special linear group Moreover, we will consider the important orthogonal, unitary, and quaternionic unitary Lie groups as well as the special orthogonal and the special unitary groups

SU(n) = U(n) ∩ SL n (C).
We recall that their generalizations O( p, q), U( p, q), Sp( p, q) are defined as the groups of linear transformations on R p+q , C p+q , H p+q , respectively, preserving a bilinear (or Hermitian) form induced by I p,q . More precisely, and we also put For the quaternionic unitary group Sp( p, q), we use the standard complex representation of GL p+q (H) in C 2( p+q)×2( p+q) given by We will also consider the symplectic group Sp(n, R) = {q ∈ SL 2n (R) | q · J n · q t = J n }.
Note that every element q ∈ Sp(n, R) has a unique representation as for some matrices x, y, z, w ∈ R n×n . Finally, we will consider the Lie groups SO * (2n) = {q ∈ SU(n, n) | q · I n,n · J n · q t = I n,n · J n }, Note that for any element q ∈ SO * (2n), we have a unique representation where z, w are n × n matrices.
Having recalled the definitions of these classical semisimple Lie groups, we are ready to present examples of eigenfunctions; see Table 1. The reader should note that in Table 1 we interpret the vectors as row vectors, so that if a, v ∈ C n , then a t v is an n × n matrix with entries a i v j . We also denote The machinery required for calculating the tension field and the conformality operator of the functions in Table 1 can be found in [5,6] and [7], and is therefore omitted from this paper. Most of the results presented here are immediate consequences of the calculations made there, the only exception being the results concerning the Lie groups SL n (R), SU(n), SU( p, q) and SU * (2n), which require a simple modification.
Let us describe this for the case of SU(n), the other cases follow the same pattern. Following Table 1, define φ : SU(n) → C by where a, v ∈ C n . Note that this function extends naturally to U(n). If we letτ andκ denote the tension field and the conformality operator of U(n), respectively, then as a consequence of Lemma 5.1 and Theorem 5.2 in [6], we get τ (φ) = −n φ andκ(φ, φ) = −φ 2 . Now it is easy to see that the matrix is of unit length and generates the orthogonal complement of su(n) in u(n). If we denote the tension field and the conformality operator of SU(n) by τ and κ, respectively, we therefore obtain as claimed in Table 1.