The aim of this paper is to prove that PHH harmonic submersions are actually (weakly) stable, yet another property which relates maps in this class to holomorphic maps (compare to [3] ).

We start by recalling some basic notions and facts from [1] , and [5] , see also [3] . We start with a map $\phi :(M^m,g_M)\to (N^{2n},J_N,g_N)$  , defined on a Riemann manifold with value in a Kähler manifold. For any point $x\in M$  , we denote $d{\phi }_x^{*}:T_{\phi \left(x\right)}N\to T_{x}M$  the adjoint of the tanget map $d{\phi }_x:T_{x}M\to T_{\phi \left(x\right)}N$  .

The map $\phi$  is called PHWC at $x$  (pseudo-horizontally weakly conformal) if and only if $d{\phi }_x\circ d{\phi }_x^{*}$  commutes with $J_{N,\phi \left(x\right)}$  . Naturally, $\phi$  is called PHWC if it is PHWC at any point $x\in M$  . If $\phi$  is PHWC at $x$  , we say that $\phi$  is PHH at $x$  (pseudo-horizontally homothetic) if and only if $$d{\phi }_x\left(\left({\nabla }_v^{M}d{\phi }_x^{*}\right(J_{N}Y){)}_x\right)=J_{N,\phi \left(x\right)}d{\phi }_x\left(\left({\nabla }_v^{M}d{\phi }_x^{*}\right(Y){)}_x\right)$$  for any horizontal tangent vector $v\in T_{x}M$  , and any vector field $Y$  , locally defined in a neighbourhood of $\phi \left(x\right)$  , where ${\nabla }^M$  is the Levi-Civita connection on $M$  (the Levi-Civita connections on $N$  is denoted by ${\nabla }^N$  , and the induced connection in the bundle ${\phi }^{-1}TN$  is denoted by $\tilde{\nabla }$  ). The map $\phi$  is called PHH if it is PHH at any point $x$  of $M$  . It is easy to see that a PHWC map is PHH if and only if $$d\phi \left({\nabla }_X^{M}d{\phi }^{*}\right(J_{N}Y\left)\right)=J_{N}d\phi \left({\nabla }_X^{M}d{\phi }^{*}\right(Y\left)\right),$$  for any horizontal vector field $X$  on $M$  and any vector field $Y$  on $N$  .

The two conditions PHWC and PHH can be seen in terms of the almost complex structure on the horizontal bundle, defined by $J_H=d{\phi }^{-1}\circ J_N\circ d\phi$  . $J_H$  is the restriction of an $f$  -structure on $M$  , also denoted by $J_H$  , which vanishes on the vertical distribution. The PHWC condition is the compatibility condition of $J_H$  with the metric on $M$  , while pseudo-horizontally homothetic condition is equivalent to $J_H$  being parallel in horizontal directions.

Suppose next that the map $\phi$  is harmonic and submersive, and $M$  is compact. In this case, we know from Theorem 2.1 (a), Proposition 3.1 and Proposition 3.3 in [1] that the fibres of $\phi$  are minimal submanifolds.

Recall next that the stability of harmonic maps is controlled by a condition on the Hessian of the energy-functional:

$$H\left(E{)}_{\phi }\right(V,V)\ge 0,$$  for all $V$  a section of the bundle ${\phi }^{-1}TN$  . The Hessian is computed by (see [7] , pp. 155) $$H\left(E{)}_{\phi }\right(V,W)={\int }_M{g}_N(J_{\phi }V,W)v_M,$$  for all $V,W$  sections in ${\phi }^{-1}TN$  , where $v_M$  is the volume form on $M$  , and $J_{\phi }$  is a second order selfadjoint elliptic differential operator acting on sections of ${\phi }^{-1}TN$  in the following way. Denote by ${}^{N}R$  the curvature tensor field on $N$  , and let $\{{\varepsilon }_1,...,{\varepsilon }_m\}$  be an orthogonal vector frame on $M$  , and $V$  be a section in ${\phi }^{-1}TN$  . Then $$J_{\phi }V:=-{\sum }_{i=1}^m\frac{1}{\left|\right|{\varepsilon }_i|{|}^2}\left({\tilde{\nabla }}_{{\varepsilon }_i}{\tilde{\nabla }}_{{\varepsilon }_i}-{\tilde{\nabla }}_{{\nabla }_{{\varepsilon }_i}{\varepsilon }_i}\right)V-{\sum }_{i=1}^m{\frac{1}{\left|\right|{\varepsilon }_i|{|}^2}}^{N}R(V,d\phi ({\varepsilon }_i\left)\right)d\phi \left({\varepsilon }_i\right).$$  The second-order elliptic differential operator $${\overline{\Delta }}_{\phi }:=-{\sum }_{i=1}^m\frac{1}{\left|\right|{\varepsilon }_i|{|}^2}\left({\tilde{\nabla }}_{{\varepsilon }_i}{\tilde{\nabla }}_{{\varepsilon }_i}-{\tilde{\nabla }}_{{\nabla }_{{\varepsilon }_i}{\varepsilon }_i}\right)$$  is called the rough Laplacian of $\phi$  , cf. [7] , pp. 155. The second sum which appears in the formula defining the Jacobi operator $J_{\phi }$  is denoted by ${\mathcal{ℛ}}_{\phi }$  , so $J_{\phi }={\overline{\Delta }}_{\phi }-{\mathcal{ℛ}}_{\phi }$  . One of the useful properties of the rough Laplacian, which will be constantly used in the sequel is the following, cf. [7] , pp. 156.

Proposition 1 The rough Laplacian ${\overline{\Delta }}_{\phi }$  satisfies $${\int }_M{g}_N({\overline{\Delta }}_{\phi }V,W)v_M={\int }_M{g}_N(\tilde{\nabla }V,\tilde{\nabla }W)v_M={\int }_M{g}_N(V,{\overline{\Delta }}_{\phi }W)v_M,$$  where $V$  and $W$  are sections on ${\phi }^{-1}TN$  , and $$g_N(\tilde{\nabla }V,\tilde{\nabla }W)={\sum }_{i=1}^m\frac{1}{\left|\right|{\varepsilon }_i|{|}^2}{g}_N({\tilde{\nabla }}_{{\varepsilon }_i}V,{\tilde{\nabla }}_{{\varepsilon }_i}W).$$ 

After these preparations, we arrive at the statement of the main result of this paper.

Theorem 2 Let $(M^m,g_M)$  be a compact Riemann manifold, $(N^{2n},J_N,g_N)$  be a Kähler manifold, and $\phi :M\to N$  be a harmonic PHH submersion. Then $\phi$  is (weakly) stable.

Proof. As in the proof of Theorem 4.1 of [1] , we choose a (local) frame $\{e_1,...,e_n,$  $J_N{e}_1,...,J_N{e}_n\}$  in ${\phi }^{-1}TN$  such that the system $\left\{d{\phi }^{*}\right(e_1),...,d{\phi }^{*}(e_n),$  $d{\phi }^{*}\left(J_N{e}_1\right),...,d{\phi }^{*}\left(J_N{e}_n\right)\}$  is an orthogonal frame in the horizontal distribution. We also choose $\{u_1,...,u_s\}$  an orthonormal basis for the vertical distribution. We denote $E_i=d{\phi }^{*}\left(e_i\right)$  , and $E_i^{\prime }=d{\phi }^{*}\left(J_N{e}_i\right)$  , for all $i=1,...,n$  .

With this notation, we apply the same strategy of proof as in [7] , pp. 172, Theorem 3.2.

For $V$  a section in ${\phi }^{-1}TN$  , we apply Proposition 1, and compute:

$$\begin{array}{ccc}H\left(E{)}_{\phi }\right(V,V)& =& {\int }_M{g}_N(\tilde{\nabla }V,\tilde{\nabla }V)v_M-{\int }_M{g}_N({\mathcal{ℛ}}_{\phi }V,V)v_M.\end{array}$$

By definition

$$\begin{array}{ccc}{g}_N(\tilde{\nabla }V,\tilde{\nabla }V)& =& {\sum }_{i=1}^n\left(\frac{1}{\left|\right|E_i|{|}^2}{g}_N({\tilde{\nabla }}_{E_i}V,{\tilde{\nabla }}_{E_i}V)+\frac{1}{\left|\right|E_i^{\prime }|{|}^2}{g}_N({\tilde{\nabla }}_{E_i^{\prime }}V,{\tilde{\nabla }}_{E_i^{\prime }}V)\right)\end{array}$$

$$\begin{array}{ccc}& & +{\sum }_{j=1}^s{g}_N({\tilde{\nabla }}_{u_j}V,{\tilde{\nabla }}_{u_j}V).\end{array}$$

Analogous to the operator used in the proof of Theorem 3.2, Chapter 5, [7] , we define, for any $V\in \Gamma \left({\phi }^{-1}TN\right)$  , the operator $DV\in \Gamma \left({\phi }^{-1}TN\otimes {\mathcal{\mathscr{H}}}^{*}\right)$  , where $\mathcal{\mathscr{H}}$  is the horizontal distribution on $M$  , by $$DV\left(X\right):={\tilde{\nabla }}_{J_{H}X}V-J_N{\tilde{\nabla }}_{X}V,$$  for any $X$  a horizontal vector field on $M$  .

Next, we compute

$$\begin{array}{ccc}{g}_N(DV,DV)& =& {\sum }_{i=1}^n\{\frac{1}{\left|\right|E_i|{|}^2}{g}_N\left(DV\right(E_i),DV(E_i\left)\right)\end{array}$$

$$\begin{array}{ccc}& & +\frac{1}{\left|\right|E_i^{\prime }|{|}^2}{g}_N\left(DV\right(E_i^{\prime }),DV(E_i^{\prime }\left)\right)\}.\end{array}$$

Since $J_H{E}_i=E_i^{\prime }$  , $J_H{E}_i^{\prime }=-E_i$  , and $\left|\right|E_i\left|\right|=\left|\right|E_i^{\prime }\left|\right|$  , we obtain

$$\begin{array}{ccc}{g}_N(DV,DV)& =& 2{\sum }_{i=1}^n\frac{1}{\left|\right|E_i|{|}^2}(g_N({\tilde{\nabla }}_{E_i}V,{\tilde{\nabla }}_{E_i}V)+g_N({\tilde{\nabla }}_{E_i^{\prime }}V,{\tilde{\nabla }}_{E_i^{\prime }}V)\end{array}$$

$$\begin{array}{ccc}& & -2g_N({\tilde{\nabla }}_{E_i^{\prime }}V,J_N{\tilde{\nabla }}_{E_i}V))\end{array}$$

Therefore

$$\begin{array}{ccc}{\int }_M\left(g_N(J_{\phi }V,V)-\frac{1}{2}{g}_N(DV,DV)\right)v_M& =& {\sum }_{i=1}^n{\int }_M\frac{1}{\left|\right|E_i|{|}^2}(2g_N({\tilde{\nabla }}_{E_i^{\prime }}V,J_N{\tilde{\nabla }}_{E_i}V)\end{array}$$

$$\begin{array}{ccc}& & -g_N{(}^{N}R(V,d\phi (E_i\left)\right)d\phi \left(E_i\right),V)\end{array}$$

$$\begin{array}{ccc}& & -g_N{(}^{N}R(V,d\phi (E_i^{\prime }\left)\right)d\phi \left(E_i^{\prime }\right),V))v_M.\end{array}$$

Next, taking into account the identities $d\phi \left(E_i^{\prime }\right)=J_{N}d\phi \left(E_i\right)$  , $d\phi \left(E_i\right)=-J_{N}d\phi \left(E_i^{\prime }\right)$  , and the basic properties of the curvature tensor field ${}^{N}R$  , we obtain $${}^{N}R(V,d\phi (E_i\left)\right)d\phi \left(E_i\right){+}^{N}R(V,d\phi (E_i^{\prime }\left)\right)d\phi \left(E_i^{\prime }\right)={J_N}^{N}R\left(d\phi \right(E_i),d\phi (E_i^{\prime }\left)\right)V,$$  and thus $${\int }_M\left(g_N(J_{\phi }V,V)-\frac{1}{2}{g}_N(DV,DV)\right)v_M$$  $$={\sum }_{i=1}^n{\int }_M\frac{1}{\left|\right|E_i|{|}^2}\left(2g_N({\tilde{\nabla }}_{E_i^{\prime }}V,J_N{\tilde{\nabla }}_{E_i}V)-g_N\left({J_N}^{N}R\right(d\phi \left(E_i\right),d\phi \left(E_i^{\prime }\right))V,V)\right)v_M.$$  We compute

$$\begin{array}{ccc}-g_N\left({J_N}^{N}R\right(d\phi \left(E_i\right),d\phi \left(E_i^{\prime }\right))V,V)=g_N{(}^{N}R\left(d\phi \right(E_i),d\phi (E_i^{\prime }\left)\right)V,J_{N}V)& & \end{array}$$

$$\begin{array}{ccc}& =g_N({\tilde{\nabla }}_{E_i}{\tilde{\nabla }}_{E_i^{\prime }}V-{\tilde{\nabla }}_{E_i^{\prime }}{\tilde{\nabla }}_{E_i}V-{\tilde{\nabla }}_{[E_i,E_i^{\prime }]}V,J_{N}V)& \end{array}$$

$$\begin{array}{ccc}& =E_i{g}_N({\tilde{\nabla }}_{E_i^{\prime }}V,J_{N}V)-E_i^{\prime }{g}_N({\tilde{\nabla }}_{E_i}V,J_{N}V)& \end{array}$$

$$\begin{array}{ccc}& -g_N({\tilde{\nabla }}_{{\nabla }_{E_i}{E}_i^{\prime }}V,J_{N}V)+g_N({\tilde{\nabla }}_{{\nabla }_{E_i^{\prime }}{E}_i}V,J_{N}V)& \end{array}$$

$$\begin{array}{ccc}& -g_N({\tilde{\nabla }}_{E_i^{\prime }}V,{\tilde{\nabla }}_{E_i}{J}_{N}V)+g_N({\tilde{\nabla }}_{E_i}V,{\tilde{\nabla }}_{E_i^{\prime }}{J}_{N}V).& \end{array}$$

Similarly to [7] , pp. 180, we define a ${\mathcal{C}}^{\infty }$  function $\phi$  on $M$  by the formula:

$$\phi :={\sum }_{i=1}^n\frac{1}{\left|\right|E_i|{|}^2}(E_i{g}_N({\tilde{\nabla }}_{E_i^{\prime }}V,J_{N}V)-E_i^{\prime }{g}_N({\tilde{\nabla }}_{E_i}V,J_{N}V)$$  $$-g_N({\tilde{\nabla }}_{{\nabla }_{E_i}{E}_i^{\prime }}V,J_{N}V)+g_N({\tilde{\nabla }}_{{\nabla }_{E_i^{\prime }}{E}_i}V,J_{N}V)).$$  Since $$g_N({\tilde{\nabla }}_{E_i}V,{\tilde{\nabla }}_{E_i^{\prime }}{J}_{N}V)=-g_N({\tilde{\nabla }}_{E_i}{J}_{N}V,{\tilde{\nabla }}_{E_i^{\prime }}V),$$  we have $${\int }_M\left(g_N(J_{\phi }V,V)-\frac{1}{2}{g}_N(DV,DV)\right)v_M={\int }_M\phi v_M.$$  The proof of the Theorem will be concluded if we prove $${\int }_M\phi v_M=0.$$  For this, we use Green's formula. We choose $X$  a horizontal vector field on $M$  defined by the property:

$$g_M(X,Y)=g_N({\tilde{\nabla }}_{J_{H}Y}V,J_{N}V),$$  for any vector field $Y$  on $M$  , and we prove $div\left(X\right)=\phi$  . Indeed, since the fibres of $\phi$  are minimal, and $X$  is horizontal, it follows:

$$div\left(X\right)={\sum }_{i=1}^n\frac{1}{\left|\right|E_i|{|}^2}\left(g_M(E_i,{\nabla }_{E_i}X)+g_M(E_i^{\prime },{\nabla }_{E_i^{\prime }}X)\right).$$  Next, $$div\left(X\right)={\sum }_{i=1}^n\frac{1}{\left|\right|E_i|{|}^2}\left(E_i{g}_M(E_i,X)-g_M({\nabla }_{E_i}{E}_i,X)+E_i^{\prime }{g}_M(E_i^{\prime },X)-g_M({\nabla }_{E_i^{\prime }}{E}_i^{\prime },X)\right)$$  $$={\sum }_{i=1}^n\frac{1}{\left|\right|E_i|{|}^2}\left(E_i{g}_N({\tilde{\nabla }}_{J_H{E}_i}V,J_{N}V)-g_N({\tilde{\nabla }}_{J_H{\nabla }_{E_i}{E}_i}V,J_{N}V)\right)$$  $$+{\sum }_{i=1}^n\frac{1}{\left|\right|E_i|{|}^2}\left(E_i^{\prime }{g}_N({\tilde{\nabla }}_{J_H{E}_i^{\prime }}V,J_{N}V)-g_N({\tilde{\nabla }}_{J_H{\nabla }_{E_i^{\prime }}{E}_i^{\prime }}V,J_{N}V)\right).$$  By the PHH condition, we have $J_H{\nabla }_{E_i^{\prime }}{E}_i^{\prime }=-({\nabla }_{E_i^{\prime }}{E}_i{)}^h$  , and $J_H{\nabla }_{E_i}{E}_i=({\nabla }_{E_i}{E}_i^{\prime }{)}^h$  , where by $(.{)}^h$  we denoted the horizontal component of $(.)$  , so, $$div\left(X\right)=\phi .$$  We proved $${\int }_M{g}_N(J_{\phi }V,V)v_M=\frac{1}{2}{g}_N(DV,DV)v_M\ge 0.$$ 

Remark 3 Our result improves the main result of [6] , provided that the source manifold is compact (condition which is not needed in [6] ).

Acknowledgements. This work was financed by a NATO fellowship. The author is grateful to the Fourier Institute in Grenoble for hospitality, and to J. C. Wood for some useful remarks on an early version of the manuscript.

References

1. M. A. Aprodu, M. Aprodu, V. Brînzănescu, A class of harmonic maps and minimal submanifolds. Int. J. Math. 11 (2000), 1177-1191.
2. M. A. Aprodu, M. Aprodu, Implicitly defined harmonic PHH submersions. Manuscripta Math. 100 (1999), 103-121.
3. D. Burns, F. Burstall, P. De Bartolomeis, J. Rawnsley, Stability of harmonic maps of Kähler manifolds. J. Differential Geom. 30 (1989), 579-594.
4. A. Lichnerowicz, Applications harmoniques et veriétés kählériennes. Symp. Math. III (Bologna 1970), 341-402.
5. E. Loubeau, Pseudo Harmonic Morphisms. Int. J. Math. 7 (1997), 943-957.
6. S. Montaldo, Stability of harmonic morphisms to a surface, Int. J. Math. 9 (1998), 865-875.
7. K. Urakawa, Calculus of variations and harmonic maps. Transl. Math. Monographs Vol. 132. AMS Providence, Rhode Island: 1993.