He considered the case of general locally symmetric base manifold and have proved that the geodesic curvatures of projection of any (non-vertical) geodesic line on the tangent sphere bundle are all constant. Nevertheless it was still interesting to find a clearer description of projections of geodesics for the case of classical rank one symmetric spaces. The second author made a first step in this defection and proved that the projection to the base space of any non-vertical geodesic line on the tangent or tangent sphere bundle of a complex space form $C{P}^n$  is a curve of constant curvatures $k_1,...,k_5$  and zero curvatures $k_6,...,k_{n-1}$  In is this paper we make a contribution in more clear understanding of geometry of projected geodesics in the case of tangent (sphere) bundle of almost all classical locally symmetric spaces, namely spheres, complex and quaternionic projective spaces and their non-compact dual from a unified viewpoint using the recurrent properties of powers of the curvature operator of these spaces. This approach allows to give also a unified proof of the results from [3] , [4] and [5] We also use an easy to prove result [1] , stating that the geodesics of tangent or tangent sphere bundle with Sasaki metric have the same projections to the base manifold.

Remark on notations. Throughout the paper $<\cdot ,\cdot >$  and $|\cdot |$  mean the scalar product and the norm of vectors with respect to the corresponding metrics.

1 Summary of main results.

Let $\left(M^n\right(c),g)$  be a Riemannian manifold of constant curvature $c$  , $\left(M^{2n}\right(c);J;g)$  a Riemannian manifold with complex structure $J$  of constant holomorphic curvature $c$  and $\left(M^{4n}\right(c);J_1,J_2,J_3;g)$  a Riemannian manifold with quaternionic structure $(J_1,J_2,J_3)$  of constant quaternionic curvature $c$  . For the sake of brevity, denote by $\mathcal{ℳ}\left(c\right)$  one of these space forms with corresponding standard metrics and will refer to $\mathcal{ℳ}\left(c\right)$  just to a space form of constant curvature $c$  . The main result is the following statement.

Theorem 1.1 Let $\mathcal{ℳ}\left(c\right)$  be a space form of constant curvature $c\ne 0$  .

Let $\Gamma$  be non-vertical geodesic line on the tangent or tangent sphere bundle over $\mathcal{ℳ}\left(c\right)$  . Let $\gamma =\pi \circ \Gamma$  be the projection of $\Gamma$  to $\mathcal{ℳ}\left(c\right)$  . Then the geodesic curvatures $k_1,k_2,...$  of $\gamma$  are all constant and

• (a) $k_3=...=k_{n-1}=0$  for the real space form;
• (b) $k_6=...=k_{2n-1}=0$  for the complex space form;
• (c) $k_{10}=...=k_{4n-1}=0$  for the quaternionic space form.

As the referee remarked, the result of the Theorem  1.1 can be expressed in more clear geometrical terms, namely the projected curve $\gamma =\pi \circ \Gamma$  lies in a totally geodesic $S^3$  or $H^3$  , in a totally geodesic $C{P}^3$  or $C{H}^3$  and in a totally geodesic $Q{P}^3$  or $Q{H}^3$  for the real, complex and quaternionic space form respectively. These assertions can be derived from ( 6 ), ( 10 ) and ( 14 ).

Proof of the Theorem  1.1 is based on the recurrent property of powers of curvature operator of spaces under consideration. Let $R_{XY}$  be the curvature operator of $\mathcal{ℳ}\left(c\right)$  . Define a power of curvature operator $R_{XY}^p$  recurrently in the following way:

$$R_{XY}^{p}Z=R_{XY}^{p-1}\left(R_{XY}Z\right)p>1.$$  The basic tool for our considerations are a chain of lemmas.

Lemma 1.1 Let $R_{XY}$  be the curvature operator of the real space form $\left(M^n\right(c),g)$  . Then for any $X$  and $Y$  $$R_{XY}^p=\{\begin{array}{c}(-b^2{c}^2{)}^{s-1}{R}_{XY}\text{for p=2s-1}\\ (-b^2{c}^2{)}^{s-1}{R}_{XY}^2\text{for p=2s},\end{array}s\ge 1$$  where $b=|X\wedge Y|$  is a norm of bivector $X\wedge Y$  .

Lemma 1.2 Let $R_{XY}$  be the curvature operator of the non-flat complex space form $\left(M^n\right(c);J;g)$  . Denote by $b=|X\wedge Y|$  the norm of a bivector $X\wedge Y$  and $m=<X,JY>$  . Then for any $X$  and $Y$  $$R_{XY}^p=\{\begin{array}{c}Lin(J{R}_{XY}^2,R_{XY},J)\text{for p=2s-1}\\ Lin(R_{XY}^2,J{R}_{XY},E)\text{for p=2s},\end{array}s\ge 2$$  where $E$  is the identity operator and $Lin$  means a linear combination of corresponding operators with coefficients being polynomials in $\frac{1}{c},b,m$  .

Lemma 1.3 Let $R_{XY}$  be the curvature operator of the non-flat quaternionic space form $\left(M^n\right(c);J_1,J_2,J_3;g)$  . Denote by $b=|X\wedge Y|$  the norm of a bivector $X\wedge Y$  . Set $m_1=<X,J_{1}Y>$  , $m_2=<X,J_{2}Y>$  , $m_3=<X,J_{3}Y>$  , $m^2=m_1^2+m_2^2+m_3^2$  , $\mathcal{J}=m_1{J}_1+m_2{J}_2+m_3{J}_3$  . Then for any $X$  and $Y$  $$R_{XY}^p=\{\begin{array}{c}Lin(\mathcal{J}{R}_{XY}^4,R_{XY}^3,\mathcal{J}{R}_{XY}^2,R_{XY},\mathcal{J})\text{for p=2s-1}\\ Lin(R_{XY}^4,\mathcal{J}{R}_{XY}^3,R_{XY}^2,\mathcal{J}{R}_{XY},E)\text{for p=2s},\end{array}s\ge 3$$  where $E$  is the identity operator and $Lin$  means a linear combination of corresponding operators with coefficients being polynomials in $\frac{1}{c},b,m$  .

2 Necessary facts and proof of the main result.

Let $(M^n,g)$  be a Riemannian manifold and $T{M}^n$  be its tangent bundle.

Denote by $(u^1,...,u^n)$  a local coordinate system on $M^n$  . Then in each tangent space of $M^n$  the natural coordinate frame $\left\{\partial /\partial u^1,...,\partial /\partial u^n\right\}$  form a local basis. Let $\xi$  be any tangent vector over the given local chart.

Then $\xi$  can be decomposed as $$\xi ={\xi }^1\partial /\partial u^1+...+{\xi }^n\partial /\partial u^n.$$  The parameters $(u^1,...,u^n;{\xi }^1,...,{\xi }^n)$  form the so called natural induced coordinate system in $T{M}^n$  . The Sasaki metric line element $d{\sigma }^2$  with respect to this coordinate system is

$$\begin{array}{c}d{\sigma }^2=d{s}^2+|D\xi {|}^2,\end{array}$$

(1)

where $d{s}^2$  is a line element of $M^n$  , $D\xi$  is the covariant differential of $\xi$  with respect to Levi-Civita connection on $M^n$  and $|\cdot |$  means the norm with respect to Riemannian metric on $M^n$  .

The tangent sphere bundle $T_1{M}^n$  can be considered as a hypersurface in the tangent bundle defined by the condition $\left|\xi \right|=1$  . We will consider $T_1{M}^n$  as a submanifold in $T{M}^n$  with the induced metric.

With respect to the natural coordinate system, each curve $\Gamma$  on $T{M}^n$  can be represented as $\Gamma \left(\sigma \right)=\left\{u^1\left(\sigma \right)...,u^n\left(\sigma \right);{\xi }^1\left(\sigma \right),...,{\xi }^n\left(\sigma \right)\right\}$  with respect to the arc-length parameter $\sigma$  and can be interpreted as the vector field $\xi \left(\sigma \right)={\xi }^1\left(\sigma \right)\partial /\partial u^1+...+{\xi }^n\left(\sigma \right)\partial /\partial u^n$  along the projected curve $\gamma =\pi \circ \Gamma =\left(u^1\right(\sigma ),...,u^n(\sigma \left)\right)$  . If $\xi$  is a unit vector field then $\Gamma$  lies in $T_1{M}^n$  and represents an arbitrary curve in $T_1{M}^n$  .

Denote by ${(}^{\prime })$  the covariant derivative along $\gamma$  with respect to parameter $\sigma$  . Then $\Gamma$  is a geodesic line on $T{M}^n$  or $T_1{M}^n$  if $\gamma$  and $\xi$  satisfy respectively the system of equations

$$\begin{array}{c}\begin{array}{c}T{M}^n:\{\begin{array}{c}{\gamma }^{\prime \prime }=R_{{\xi }^{\prime }\xi }{\gamma }^{\prime },\\ {\xi }^{\prime \prime }=0;\end{array}{T}_1{M}^n:\{\begin{array}{c}{\gamma }^{\prime \prime }=R_{{\xi }^{\prime }\xi }{\gamma }^{\prime },\\ {\xi }^{\prime \prime }=-{\rho }^2\xi ,\end{array}\end{array}\end{array}$$

(2)

where ${\rho }^2=|{\xi }^{\prime }{|}^2$  and $R_{{\xi }^{\prime }\xi }$  is the curvature operator of $M^n$  based on bivector ${\xi }^{\prime }\wedge \xi$  .

From ( 2 ) it follows that $\rho =const$  in both cases. Denote by $s$  the arc-length parameter on $\gamma$  . Then from ( 1 ) it follows that

$$\begin{array}{c}\frac{ds}{d\sigma }=\sqrt{1-{\rho }^2},\end{array}$$

(3)

so that $0\le \rho \le 1$  . According to the latter inequality, the set of geodesics of $T{M}^n$  and $T_1{M}^n$  can be splitted naturally into 3 classes, namely,

• ∙ horizontal geodesics ( $\rho =0$  ) generated by parallel (unit) vector fields along the geodesics on the base manifold;
• ∙ vertical geodesics ( $\rho =1$  ) represented by geodesics on a fixed fiber;
• ∙ umbilical geodesics corresponding to $0<\rho <1$  .

In what follows we will consider the properties of projections of umbilical geodesics.

Lemma 2.1 (cf. [2] )Let $(M^n,g)$  be a locally symmetric Riemannian manifold and $R_{XY}$  its curvature operator. Let $\gamma =\pi \circ \Gamma$  be a projection of geodesic line on $T{M}^n$  or $T_1{M}^n$  to the base space. Then for the derivatives of $\gamma$  of order $p$  we have $${\gamma }^{\left(p\right)}=R_{{\xi }^{\prime }\xi }^{p-1}{\gamma }^{\prime }=R_{{\xi }^{\prime }\xi }{\gamma }^{(p-1)}$$  and as a consequence all the geodesic curvatures of $\gamma$  are constant.

Proof. The equalities follow from parallelism of curvature tensor of $M^n$  and the equations ( 2 ). Moreover, from the evident identity $$<{\gamma }^{\left(p\right)},{\gamma }^{(p-1)}>\equiv 0$$  for all $p>1$  , we conclude that $\left|{\gamma }^{\left(p\right)}\right|=const$  for all $p>1$  and therefore, by induction, all the geodesic curvatures of $\gamma$  are constant.

Proof of Theorem  1.1 . Case (a). Denote by $e_1,...,e_{n-1}$  the Frenet frame of $\gamma$  . Using the Frenet formulas for the curve with constant geodesic curvatures and keeping in mind ( 3 ), it is easy to see that

$$\begin{array}{c}\begin{array}{c}{\gamma }^{(2s-1)}=(1-{\rho }^2{)}^{s-1/2}{k}_1{k}_2...k_{2s-2}{e}_{2s-1}+Lin\left\{e_1,e_3,...,e_{2s-3}\right\},\\ {\gamma }^{\left(2s\right)}=(1-{\rho }^2{)}^s{k}_1{k}_2...k_{2s-1}{e}_{2s}+Lin\left\{e_2,e_4,...,e_{2s-2}\right\}\end{array}\end{array}$$

(4)

for all $s\ge 1$  ( with formal setting $k_0\equiv 1$  ). Setting $s=1,2$  in even derivatives, we see that

$$\begin{array}{c}\begin{array}{c}{\gamma }^{\left(2\right)}=(1-{\rho }^2)k_1{e}_2\\ {\gamma }^{\left(4\right)}=(1-{\rho }^2{)}^2{k}_1{k}_2{k}_3{e}_4+Lin(e_2).\end{array}\end{array}$$

(5)

On the other hand, applying Lemma  2.1 , Lemma  1.1 and Lemma  2.1 again, we get

$$\begin{array}{c}{\gamma }^{\left(4\right)}=R_{{\xi }^{\prime }\xi }^3{\gamma }^{\prime }=-b^2{c}^2{R}_{{\xi }^{\prime }\xi }{\gamma }^{\prime }=-b^2{c}^2{\gamma }^{\left(2\right)}.\end{array}$$

(6)

Substitution from ( 5 ) gives $$(1-{\rho }^2{)}^2{k}_1{k}_2{k}_3{e}_4+Lin(e_2)=0$$  and therefore $k_3=0$  , which completes the proof.

Remark, that $b^2$  is constant along $\gamma$  since $${\left(b^2\right)}^{\prime }=\left(|{\xi }^{\prime }\wedge \xi {|}^2\right)={\left({\rho }^2|\xi {|}^2-<{\xi }^{\prime },\xi {>}^2\right)}^{\prime }=2{\rho }^2<{\xi }^{\prime },\xi >-2<{\xi }^{\prime },\xi >{\rho }^2\equiv 0.$$  Case (b). Denote by $e_1,...,e_{2n-1}$  the Frenet frame of $\gamma$  . Similar to the case (a) considerations, the Frenet formulas give

$$\begin{array}{c}\begin{array}{c}{\gamma }^{(2s-1)}=(1-{\rho }^2{)}^{s-1/2}{k}_1{k}_2...k_{2s-2}{e}_{2s-1}+Lin\left\{e_1,e_3,...,e_{2s-3}\right\},\\ {\gamma }^{\left(2s\right)}=(1-{\rho }^2{)}^s{k}_1{k}_2...k_{2s-1}{e}_{2s}+Lin\left\{e_2,e_4,...,e_{2s-2}\right\}\end{array}\end{array}$$

(7)

for all $s\ge 1$  . Setting $s=1,2,3,4$  in odd derivatives, we get

$$\begin{array}{c}\begin{array}{c}{\gamma }^{\prime }=(1-{\rho }^2{)}^{1/2}{e}_1,\\ {\gamma }^{\left(3\right)}=(1-{\rho }^2{)}^{3/2}{k}_1{k}_2{e}_3+Lin(e_1),\\ {\gamma }^{\left(5\right)}=(1-{\rho }^2{)}^{5/2}{k}_1...k_4{e}_5+Lin(e_1,e_3),\\ {\gamma }^{\left(7\right)}=(1-{\rho }^2{)}^{7/2}{k}_1...k_6{e}_7+Lin(e_1,e_3,e_5).\end{array}\end{array}$$

(8)

On the other hand, applying Lemma  2.1 , Lemma  1.2 and Lemma  2.1 again, we get

$$\begin{array}{c}\{\begin{array}{c}{\gamma }^{\left(5\right)}=R_{{\xi }^{\prime }\xi }^4{\gamma }^{\prime }=Lin(R_{{\xi }^{\prime }\xi }^2,J{R}_{{\xi }^{\prime }\xi },E){\gamma }^{\prime }=Lin({\gamma }^{\left(3\right)},J{\gamma }^{\left(2\right)},{\gamma }^{\prime }),\\ {\gamma }^{\left(7\right)}=R_{{\xi }^{\prime }\xi }^6{\gamma }^{\prime }=Lin(R_{{\xi }^{\prime }\xi }^2,J{R}_{{\xi }^{\prime }\xi },E){\gamma }^{\prime }=Lin({\gamma }^{\left(3\right)},J{\gamma }^{\left(2\right)},{\gamma }^{\prime }),\\ \end{array}\end{array}$$

(9)

Excluding $J{\gamma }^{\left(2\right)}$  from ( 9 ), we come to the equation

$$\begin{array}{c}{\gamma }^{\left(7\right)}=Lin({\gamma }^{\left(5\right)},{\gamma }^{\left(3\right)},{\gamma }^{\prime }).\end{array}$$

(10)

Substitution from ( 8 ) imply $$(1-{\rho }^2{)}^{7/2}{k}_1...k_6{e}_7+Lin(e_1,e_3,e_5)=0$$  and we conclude that $k_6=0$  which completes the proof.

Remark, that the coefficients of all linear combinations are constants.

Indeed, by Lemma  1.2 the coefficients are polynomials in $1/c,b=|{\xi }^{\prime }\wedge \xi |$  and $m=<{\xi }^{\prime },J\xi >$  . The value $b$  is constant along $\gamma$  by the same reasons as in case (a). The value $m$  is constant along $\gamma$  since $$m^{\prime }=<{\xi }^{\prime },J\xi {>}^{\prime }=<{\xi }^{\prime \prime },J\xi >+<{\xi }^{\prime },J{\xi }^{\prime }>\equiv 0.$$  Case (c). Denote by $e_1,...,e_{4n-1}$  the Frenet frame of $\gamma$  . As above, the Frenet formulas give

$$\begin{array}{c}\begin{array}{c}{\gamma }^{(2s-1)}=(1-{\rho }^2{)}^{s-1/2}{k}_1{k}_2...k_{2s-2}{e}_{2s-1}+Lin\left\{e_1,e_3,...,e_{2s-3}\right\},\\ {\gamma }^{\left(2s\right)}=(1-{\rho }^2{)}^s{k}_1{k}_2...k_{2s-1}{e}_{2s}+Lin\left\{e_2,e_4,...,e_{2s-2}\right\}\end{array}\end{array}$$

(11)

for all $s\ge 1$  . Setting $s=1,2,3,4,5,6$  in odd derivatives, we get

$$\begin{array}{c}\begin{array}{c}{\gamma }^{\prime }=(1-{\rho }^2{)}^{1/2}{e}_1,\\ {\gamma }^{\left(3\right)}=(1-{\rho }^2{)}^{3/2}{k}_1{k}_2{e}_3+Lin(e_1),\\ {\gamma }^{\left(5\right)}=(1-{\rho }^2{)}^{5/2}{k}_1...k_4{e}_5+Lin(e_1,e_3),\\ {\gamma }^{\left(7\right)}=(1-{\rho }^2{)}^{7/2}{k}_1...k_6{e}_7+Lin(e_1,e_3,e_5),\\ {\gamma }^{\left(9\right)}=(1-{\rho }^2{)}^{9/2}{k}_1...k_8{e}_9+Lin(e_1,e_3,e_5,e_7),\\ {\gamma }^{\left(11\right)}=(1-{\rho }^2{)}^{11/2}{k}_1...k_{10}{e}_{11}+Lin(e_1,e_3,e_5,e_7,e_9).\end{array}\end{array}$$

(12)

Applying again Lemma  2.1 , Lemma  1.3 and then Lemma  2.1 , we get

$$\begin{array}{c}\{\begin{array}{cc}{\gamma }^{\left(7\right)}=R_{{\xi }^{\prime }\xi }^6{\gamma }^{\prime }=& Lin(R_{{\xi }^{\prime }\xi }^4,\mathcal{J}{R}_{{\xi }^{\prime }\xi }^3,R_{{\xi }^{\prime }\xi }^2,\mathcal{J}{R}_{{\xi }^{\prime }\xi },E){\gamma }^{\prime }=\\ & Lin({\gamma }^{\left(5\right)},\mathcal{J}{\gamma }^{\left(4\right)},{\gamma }^{\left(3\right)},\mathcal{J}{\gamma }^{\left(2\right)},{\gamma }^{\prime }),\\ {\gamma }^{\left(9\right)}=R_{{\xi }^{\prime }\xi }^8{\gamma }^{\prime }=& Lin(R_{{\xi }^{\prime }\xi }^4,\mathcal{J}{R}_{{\xi }^{\prime }\xi }^3,R_{{\xi }^{\prime }\xi }^2,\mathcal{J}{R}_{{\xi }^{\prime }\xi },E){\gamma }^{\prime }=\\ & Lin({\gamma }^{\left(5\right)},\mathcal{J}{\gamma }^{\left(4\right)},{\gamma }^{\left(3\right)},\mathcal{J}{\gamma }^{\left(2\right)},{\gamma }^{\prime }),\\ {\gamma }^{\left(11\right)}=R_{{\xi }^{\prime }\xi }^{10}{\gamma }^{\prime }=& Lin(R_{{\xi }^{\prime }\xi }^4,\mathcal{J}{R}_{{\xi }^{\prime }\xi }^3,R_{{\xi }^{\prime }\xi }^2,\mathcal{J}{R}_{{\xi }^{\prime }\xi },E){\gamma }^{\prime }=\\ & Lin({\gamma }^{\left(5\right)},\mathcal{J}{\gamma }^{\left(4\right)},{\gamma }^{\left(3\right)},\mathcal{J}{\gamma }^{\left(2\right)},{\gamma }^{\prime }).\end{array}\end{array}$$

(13)

Excluding $\mathcal{J}{\gamma }^{\left(2\right)}$  and $\mathcal{J}{\gamma }^{\left(4\right)}$  from ( 13 ), we come to the equation

$$\begin{array}{c}{\gamma }^{\left(11\right)}=Lin({\gamma }^{\left(9\right)},{\gamma }^{\left(7\right)},{\gamma }^{\left(5\right)},{\gamma }^{\left(3\right)},{\gamma }^{\prime }).\end{array}$$

(14)

Substitution from ( 12 ) imply $$(1-{\rho }^2{)}^{11/2}{k}_1...k_{10}{e}_{11}+Lin(e_1,e_3,e_5,e_7,e_9)=0$$  and we conclude that $k_{10}=0$  which completes the proof.

Remark, that the coefficients of all linear combinations are constants.

Indeed, by Lemma  1.3 the coefficients are polynomials in $1/c,b=|{\xi }^{\prime }\wedge \xi |$  and $m=\sqrt{m_1^2+m_2^2+m_3^2}$  . The value $b$  is constant along $\gamma$  by the same reasons as in case (a). The values $m_1,m_2,m_3$  are all constant along $\gamma$  since $$m_i^{\prime }=<{\xi }^{\prime },J_i\xi {>}^{\prime }=<{\xi }^{\prime \prime },J_i\xi >+<{\xi }^{\prime },J_i{\xi }^{\prime }>\equiv 0$$  for $i=1,2,3$  .

3 Proofs of basic Lemmas

Proof of Lemma  1.1 .The curvature operator $R_{XY}$  of the real space form $\left(M^n\right(c),g)$  has the following expression $$R_{XY}Z=c\left[<Y,Z>X-<X,Z>Y\right].$$  Then $R_{XY}^{2}Z=c\left[<Y,R_{XY}Z>X-<X,R_{XY}Z>Y\right]=c^2\left[<Y,<Y,Z>X-<X,Z>Y>X-<X,<Y,Z>X-<X,Z>Y>Y\right]=c^2\left[<Y,Z><X,Y>X-<X,Z>|Y{|}^{2}X-<Y,Z>|X{|}^{2}Y+<X,Z><X,Y>Y\right]=c^2\left[<Y,Z>\left(<X,Y>X-|X{|}^{2}Y\right)+<X,Z>\left(<X,Y>Y-|Y{|}^{2}X\right)\right]=c\left[<Y,Z>R_{XY}X+<X,Z>R_{YX}Y\right].$  Therefore, $R_{XY}^{3}Z=c\left[<Y,R_{XY}Z>R_{XY}X+<X,R_{XY}Z>R_{YX}Y\right]=c^3\left[\left(<Y,Z><X,Y>-<X,Z>|Y{|}^2\right)\left(<X,Y>X-|X{|}^{2}Y\right)+\left(<Y,Z>|X{|}^2-<X,Z><X,Y>\right)\left(<X,Y>Y-|Y{|}^{2}X\right)\right]=c^3\left[-<Y,Z>X\left(\left|X{|}^2\right|Y{|}^2-<X,Y{>}^2\right)+<X,Z>Y\left(\left|X{|}^2\right|Y{|}^2-<X,Y{>}^2\right)\right]=-c^2{b}^2{R}_{XY}Z,$  where, evidently, $b^2=\left|X{|}^2\right|Y{|}^2-<X,Y{>}^2$  is the square norm of $X\wedge Y$  .

Now we can find the other powers for $R_{XY}$  inductively.

Proof of Lemma  1.2  The curvature operator $R_{XY}$  of the complex space form $\left(M^{2n}\right(c);J;g)$  has the following expression $$R_{XY}Z=\frac{c}{4}\left[<Y,Z>X-<X,Z>Y+<JY,Z>JX-<JX,Z>JY+2<X,JY>JZ\right].$$  Introduce the unit sphere type operator $S$  acting as $$S\left(Z\right)\stackrel{def}{=}{S}_{XY}Z=<Y,Z>X-<X,Z>Y,$$  and the operator $\hat{S}\left(Z\right)$  acting as $$\hat{S}\left(Z\right)\stackrel{def}{=}{S}_{JXJY}Z=<JY,Z>JX-<JX,Z>JY.$$  Finally, if we denote $m=<X,JY>$  , then the curvature operator under consideration takes the form

$$\begin{array}{c}{R}_{XY}Z=\frac{c}{4}\left[S+\hat{S}+2mJ\right]Z.\end{array}$$

(15)

Since $|X\wedge Y|=\left|\right(JX)\wedge (JY\left)\right|$  , the operators $S$  and $\hat{S}$  satisfy $$S^3=-b^{2}S,{\hat{S}}^3=-b^2\hat{S},$$  where $b^2=|X\wedge Y{|}^2$  .

In what follows, we need a ”table of products” for the operators $S$  and $\hat{S}$  . Namely,

$$\begin{array}{c}\begin{array}{cccc}& S& \hat{S}& J\\ S& S^2& mJ\hat{S}& J\hat{S}\\ \hat{S}& mJS& {\hat{S}}^2& JS\\ J& JS& J\hat{S}& -E\\ \end{array}\end{array}$$

(16)

Indeed,

$$\begin{array}{c}\begin{array}{c}\left(S\hat{S}\right)\left(Z\right)=S_{XY}{\hat{S}}_{JXJY}Z=S_{XY}\left[<JY,Z>JX-<JX,Z>JY\right]=\\ <Y,<JY,Z>JX-<JX,Z>JY>X-<X,<JY,Z>JX-<JX,Z>JY>Y=\\ <Y,JX><JY,Z>X+<JX,Z><X,JY>Y=\\ mJ\left[<JY,Z>JX-<JX,Z>JY\right]=\left(mJ\hat{S}\right)\left(Z\right),\end{array}\end{array}$$

$$\begin{array}{c}\begin{array}{c}\left(\hat{S}S\right)\left(Z\right)=S_{JXJY}{S}_{XY}Z=S_{JXJY}\left[<Y,Z>X-<X,Z>Y\right]=\\ <JY,<Y,Z>X-<X,Z>Y>JX-<JX,<Y,Z>X-<X,Z>Y>JY=\\ <JY,X><Y,Z>JX+<X,Z><JX,Y>JY=\\ mJ\left[<Y,Z>X-<X,Z>Y\right]=\left(mJS\right)\left(Z\right),\end{array}\end{array}$$

$$\begin{array}{c}\begin{array}{c}\left(SJ\right)\left(Z\right)=S_{XY}JZ=<Y,JZ>X-<X,JZ>Y=\\ J\left[<JY,Z>JX-<JX,Z>JY\right]=\left(J\hat{S}\right)\left(Z\right),\end{array}\end{array}$$

$$\begin{array}{c}\begin{array}{c}\left(\hat{S}J\right)\left(Z\right)=S_{JXJY}JZ=<JY,JZ>JX-<JX,JZ>JY=\\ J\left[<Y,Z>X-<X,Z>Y\right]=\left(JS\right)\left(Z\right),\end{array}\end{array}$$

and the other entries of the table can be found in a similar way.

From ( 16 ) we see that $J(S+\hat{S})=(S+\hat{S})J$  and we have

$$\begin{array}{c}(S+\hat{S}{)}^2=S^2+{\hat{S}}^2+S\hat{S}+\hat{S}S=S^2+{\hat{S}}^2+mJ(S+\hat{S})\end{array}$$

$$\begin{array}{c}\begin{array}{cc}(S+\hat{S}{)}^3=& (S+\hat{S})[S^2+{\hat{S}}^2+mJ(S+\hat{S}\left)\right]=\\ & S^3+{\hat{S}}^3+\hat{S}{S}^2+S{\hat{S}}^2+mJ(S+\hat{S}{)}^2=\\ & -b^{2}S-b^2\hat{S}+\left(\hat{S}S\right)S+\left(S\hat{S}\right)\hat{S}+mJ(S+\hat{S}{)}^2=\\ & -b^2(S+\hat{S})+mJ(S^2+{\hat{S}}^2)+mJ(S+\hat{S}{)}^2=\\ & -b^2(S+\hat{S})+mJ\left[\right(S+\hat{S}{)}^2-mJ(S+\hat{S})]+\\ & mJ(S+\hat{S}{)}^2=(m^2-b^2\left)\right(S+\hat{S})+2mJ(S+\hat{S}{)}^2.\end{array}\end{array}$$

Thus,

$$\begin{array}{c}(S+\hat{S}{)}^3=Lin(S+\hat{S},J(S+\hat{S}{)}^2)\end{array}$$

(17)

On the other hand, setting for brevity $R_{XY}=R$  , from ( 15 ) we derive

$$\begin{array}{c}\begin{array}{c}S+\hat{S}=\frac{4}{c}R-2mJ=Lin(R,J),\\ (S+\hat{S}{)}^2=Lin(R^2,JR,E).\end{array}\end{array}$$

(18)

Comparing ( 17 ) and ( 18 ) we conclude $$(S+\hat{S}{)}^3=Lin\left[Lin(R,J),JLin(R^2,JR,E)\right]=Lin(J{R}^2,R,J).$$  On the other hand, from  18  $({)}_1$  $$(S+\hat{S}{)}^3={\left(\frac{4}{c}\right)}^3{R}^3+Lin(J{R}^2,R,J).$$  So, finally $$R^3=Lin(J{R}^2,R,J).$$  It is easy to trace that the coefficients of all linear combinations are polynomials in $\frac{1}{c},b,m$  . To complete the proof we should remark that

$$\begin{array}{cc}{R}^4=R^{3}R=& Lin(J{R}^2,R,J)R=Lin(J{R}^3,R^2,JR)=\end{array}$$

$$\begin{array}{cc}& Lin\left[JLin(J{R}^2,R,J),R^2,JR\right]=Lin(R^2,JR,E)\end{array}$$

which allows to find all powers of $R$  inductively.

Proof of Lemma  1.3  The curvature operator $R_{XY}$  of the quaternionic space form $\left(M^{4n}\right(c);J_1,J_2,J_3;g)$  has the following expression

$$\begin{array}{cc}{R}_{XY}Z=& \frac{c}{4}[<Y,Z>X-<X,Z>Y+<J_{1}Y,Z>J_{1}X-<J_{1}X,Z>J_{1}Y+\end{array}$$

$$\begin{array}{cc}& <J_{2}Y,Z>J_{2}X-<J_{2}X,Z>J_{2}Y+<J_{3}Y,Z>J_{3}X-<J_{3}X,Z>J_{3}Y+\end{array}$$

$$\begin{array}{cc}& 2<X,J_{1}Y>J_{1}Z+2<X,J_{2}Y>J_{2}Z+2<X,J_{3}Y>J_{3}Z].\end{array}$$

where $J_1,J_2,J_3$  are operators of quaternionic structure $$J_1{J}_2=J_3,J_2{J}_3=J_1,J_3{J}_1=J_2,J_i^2=-E,<X,J_{i}Y>=-<J_{i}X,Y>,i=\overline{1,3}$$  Introduce the unit sphere type operator $S$  acting as $$S\left(Z\right)\stackrel{def}{=}{S}_{XY}Z=<Y,Z>X-<X,Z>Y,$$  the operators $S_i\left(Z\right)$  acting as $$S_i\left(Z\right)\stackrel{def}{=}{S}_{J_{i}X{J}_{i}Y}Z=<J_{i}Y,Z>J_{i}X-<J_{i}X,Z>J_{i}Y,i=\overline{1,3},$$  and the operator $\hat{S}\left(Z\right)$  acting as $$\hat{S}\left(Z\right)\equiv S_1\left(Z\right)+S_2\left(Z\right)+S_3\left(Z\right).$$  Finally, denote $m_i=<X,J_{i}Y>,i=\overline{1,3}$  , $m^2=m_1^2+m_2^2+m_3^2$  , $\mathcal{J}=m_1{J}_1+m_2{J}_2+m_3{J}_3$  . Then the curvature operator under consideration takes the form

$$\begin{array}{c}{R}_{XY}Z=\frac{c}{4}\left[S+\hat{S}+2\mathcal{J}\right]Z.\end{array}$$

(19)

Since $|X\wedge Y|=\left|\right(J_{i}X)\wedge (J_{i}Y\left)\right|(i=\overline{1,3})$  , the operators $S$  and $S_i$  satisfy $$S^3=-b^{2}S,S_i^3=-b^2{S}_i(i=\overline{1,3}),$$  where $b^2=|X\wedge Y{|}^2$  .

The table of products for the operators $S$  and $\hat{S}$  is the following one.

$$\begin{array}{c}\begin{array}{cccccccc}& S& S_1& S_2& S_3& J_1& J_2& J_3\\ S& S^2& m_1{J}_1{S}_1& m_2{J}_2{S}_2& m_3{J}_3{S}_3& J_1{S}_1& J_2{S}_2& J_3{S}_3\\ S_1& m_1{J}_{1}S& S_1^2& -m_3{J}_3{S}_2& -m_2{J}_2{S}_3& J_{1}S& J_2{S}_3& J_3{S}_2\\ S_2& m_2{J}_{2}S& -m_3{J}_3{S}_1& S_2^2& -m_1{J}_1{S}_3& J_1{S}_3& J_{2}S& J_3{S}_1\\ S_3& m_3{J}_{3}S& -m_2{J}_2{S}_1& -m_1{J}_1{S}_2& S_3^2& J_1{S}_2& J_2{S}_1& J_{3}S\\ J_1& S_1{J}_1& S{J}_1& S_3{J}_1& S_2{J}_1& -E& J_3& -J_2\\ J_2& S_2{J}_2& S_3{J}_2& S{J}_2& S_1{J}_2& -J_3& -E& J_1\\ J_3& S_3{J}_3& S_2{J}_3& S_1{J}_3& S{J}_3& J_2& -J_1& -E\\ \end{array}\end{array}$$

(20)

The expressions for products $S{S}_i,S_{i}S,S{J}_i$  one can find similar to the table ( 16 ) making formal replacements $\hat{S}\to S_i$  and $J\to J_i$  . As concerns the other entries, we have

$$\begin{array}{c}\left(S_1{S}_2\right)\left(Z\right)=\end{array}$$

$$\begin{array}{c}{S}_{J_{1}X{J}_{1}Y}{S}_{J_{2}X{J}_{2}Y}Z=S_{J_{1}X{J}_{1}Y}\left[<J_{2}Y,Z>J_{2}X-<J_{2}X,Z>J_{2}Y\right]=\end{array}$$

$$\begin{array}{c}<J_{1}Y,<J_{2}Y,Z>J_{2}X-<J_{2}X,Z>J_{2}Y>J_{1}X-\end{array}$$

$$\begin{array}{c}<J_{1}X,<J_{2}Y,Z>J_{2}X-<J_{2}X,Z>J_{2}Y>J_{1}Y=\end{array}$$

$$\begin{array}{c}<J_{1}Y,J_{2}X><J_{2}Y,Z>J_{1}X+<J_{1}X,J_{2}Y><J_{2}X,Z>J_{1}Y=\end{array}$$

$$\begin{array}{c}{J}_1\left[<X,J_{3}Y><J_{2}Y,Z>X-<X,J_{3}Y><J_{2}X,Z>Y\right]=\end{array}$$

$$\begin{array}{c}-J_1{J}_2\left[m_3<J_{2}Y,Z>J_{2}X-m_3<J_{2}X,Z>J_{2}Y\right]=(-m_3{J}_3{S}_2)\left(Z\right),\end{array}$$

$$\begin{array}{c}\left(S_1{J}_1\right)\left(Z\right)=S_{J_{1}X{J}_{1}Y}{J}_{1}Z=<J_{1}Y,J_{1}Z>J_{1}X-<J_{1}X,J_{1}Z>J_{1}Y=\end{array}$$

$$\begin{array}{c}{J}_1\left[<Y,Z>X-<X,Z>Y\right]=\left(J_{1}S\right)\left(Z\right),\end{array}$$

$$\begin{array}{c}\left(S_1{J}_2\right)\left(Z\right)=S_{J_{1}X{J}_{1}Y}{J}_{2}Z=<J_{1}Y,J_{2}Z>J_{1}X-<J_{1}X,J_{2}Z>J_{1}Y=\end{array}$$

$$\begin{array}{c}{J}_1\left[<J_{3}Y,Z>X-<J_{3}X,Z>Y\right]=\end{array}$$

$$\begin{array}{c}-J_1{J}_3\left[<J_{3}Y,Z>J_{3}X-<J_{3}X,Z>J_{3}Y\right]=\left(J_2{S}_3\right)\left(Z\right)\end{array}$$

and so on.

From ( 20 ) we see that

$$\begin{array}{c}(S+\hat{S})\mathcal{J}=(S+S_1+S_2+S_3)(m_1{J}_1+m_2{J}_2+m_3{J}_3)=\end{array}$$

$$\begin{array}{c}{m}_1{J}_1{S}_1+m_2{J}_2{S}_2+m_3{J}_3{S}_3+m_1{J}_{1}S+m_2{J}_2{S}_3+m_3{J}_3{S}_2+m_1{J}_1{S}_3+\end{array}$$

$$\begin{array}{c}{m}_2{J}_{2}S+m_3{J}_3{S}_1+m_1{J}_1{S}_2+m_2{J}_2{S}_1+m_3{J}_{3}S=\end{array}$$

$$\begin{array}{cc}(m_1{J}_1+m_2{J}_2+m_3{J}_3)(S+S_1+S_2+S_3)=\mathcal{J}(S+\hat{S}).& \end{array}$$

Therefore, the operators $(S+\hat{S})$  $\mathcal{J}$  commute and hence for the operator $R=\frac{c}{4}\left\{\right(S+\hat{S})+2\mathcal{J}\}$  the usual formula for powers can be applied:

$$R^n={\left(\frac{c}{4}\right)}^n{\sum }_{l=0}^n\left(\genfrac{}{}{0ex}{}{n}{l}\right){\left(S+\hat{S}\right)}^{n-l}{2}^l{\left(\mathcal{J}\right)}^l.$$  The powers for $\mathcal{J}$  can be found trivially, since

$$\begin{array}{cc}{\mathcal{J}}^2=& m_1^2{J}_1^2+m_1{m}_2(J_1{J}_2+J_2{J}_1)+m_1{m}_3(J_1{J}_3+J_3{J}_1)+m_2^2{J}_2^2+\end{array}$$

$$\begin{array}{cc}& m_2{m}_3(J_2{J}_3+J_3{J}_2)+m_3^2{J}_3^2=-m_1^{2}E-m_2^{2}E-m_3^{2}E=-m^{2}E,\end{array}$$

where $m_1^2+m_2^2+m_3^2=m^2$  .

As concerns the powers of $S+\hat{S}$  , the following proposition gives the answer.

Proposition 3.1 The operator $S+\hat{S}$  possesses the recurrent property $$(S+\hat{S}{)}^5=-2(b^2+m^2\left)\right(S+\hat{S}{)}^3-(b^2-m^2)(S+\hat{S}),$$  where $b^2=|X\wedge Y{|}^2$  and $m^2=m_1^2+m_2^2+m_3^2=<X,J_{1}Y{>}^2+<X,J_{2}Y{>}^2+<X,J_{3}Y{>}^2$  .

Proof. The proof is technical and in what follows we will use some auxiliary operator products. Namely,

$$\begin{array}{c}\begin{array}{c}\begin{array}{cc}S\hat{S}=S\mathcal{J},& \hat{S}S=\mathcal{J}S,\\ S\mathcal{J}S=-m^{2}S,& S\hat{S}\mathcal{J}=-m^{2}S,\\ S(S_1^2+S_2^2+S_3^2)=S^2\mathcal{J},& \hat{S}{S}^2=\mathcal{J}{S}^2,\\ \hat{S}S\mathcal{J}=\mathcal{J}S\mathcal{J},& \hat{S}\mathcal{J}S=\mathcal{J}{S}^2,\\ \end{array}\\ \begin{array}{c}{\hat{S}}^2=S_1^2+S_2^2+S_3^2-\hat{S}\mathcal{J}+\mathcal{J}S,\\ \hat{S}(S_1^2+S_2^2+S_3^2)=-b^2\hat{S}-(S_1^2+S_2^2+S_3^2)\mathcal{J}+\mathcal{J}{S}^2.\end{array}\end{array}\end{array}$$

(21)

The proof is straightforward. Applying ( 20 ), we get $S\hat{S}=S(S_1+S_2+S_3)=m_1{J}_1{S}_1+m_2{J}_2{S}_2+m_3{J}_3{S}_3=m_{1}S{J}_1+m_{2}S{J}_2+m_{3}S{J}_3=S\mathcal{J}.$  In a similar way we find $\hat{S}S=(S_1+S_2+S_3)S=m_1{J}_{1}S+m_2{J}_{2}S+m_3{J}_{3}S=\mathcal{J}S,$  ${\hat{S}}^2=(S_1+S_2+S_3)(S_1+S_2+S_3)=S_1^2+S_2^2+S_3^2+S_1{S}_2+S_1{S}_3+S_2{S}_1+S_2{S}_3+S_3{S}_1+S_3{S}_2=S_1^2+S_2^2+S_3^2-m_3{S}_1{J}_3-m_2{S}_1{J}_2-m_3{S}_2{J}_3-m_1{S}_2{J}_2-m_2{S}_3{J}_2-m_1{S}_3{J}_1=S_1^1+S_2^2+S_3^2-\hat{S}\mathcal{J}+m_1{S}_1{J}_1+m_2{S}_2{J}_2+m_3{S}_3{J}_3=S_1^2+S_2^2+S_3^2-\hat{S}\mathcal{J}+\mathcal{J}S,$  $S\mathcal{J}S=S(m_1{J}_1+m_2{J}_2+m_3{J}_3)S=(m_1{J}_1{S}_1+m_2{J}_2{S}_2+m_3{J}_3{S}_3)S=-m_1^{2}S-m_2^{2}S-m_3^{2}S=-m^{2}S,$  $S\hat{S}\mathcal{J}=S\mathcal{J}\mathcal{J}=-m^{2}S,$  $S(S_1^2+S_2^2+S_3^2)=m_{1}S{J}_1{S}_1+m_{2}S{J}_2{S}_2+m_{3}S{J}_3{S}_3=m_1{S}^2{J}_1+m_2{S}^2{J}_2+m_3{S}^2{J}_3=S^2\mathcal{J},$  $\hat{S}\mathcal{J}S=(S_1+S_2+S_3)(m_1{J}_1+m_2{J}_2+m_3{J}_3)S=(m_1{J}_{1}S+m_2{J}_2{S}_3+m_3{J}_3{S}_2+m_1{J}_1{S}_3+m_2{J}_{2}S+m_3{J}_3{S}_1+m_1{J}_1{S}_2+m_2{J}_2{S}_1+m_3{J}_{3}S)S=\mathcal{J}{S}^2+m_2{J}_2{m}_3{J}_{3}S+m_3{J}_3{m}_2{J}_{2}S+m_1{J}_1{m}_3{J}_{3}S+m_3{J}_3{m}_1{J}_{1}S+m_1{J}_1{m}_2{J}_{2}S+m_2{J}_2{m}_1{J}_{1}S=\mathcal{J}{S}^2,$  $\hat{S}(S_1^2+S_2^2+S_3^2)=(S_1+S_2+S_3)(S_1^2+S_2^2+S_3^2)=S_1^3+S_2^3+S_3^3+S_1{S}_2^2+S_1{S}_3^2+S_2{S}_1^2+S_2{S}_3^2+S_3{S}_1^2+S_3{S}_2^2=-b^2\hat{S}-m_3{J}_3{S}_2^2-m_2{J}_2{S}_3^2-m_3{J}_3{S}_1^2-m_1{J}_1{S}_3^2-m_2{J}_2{S}_1^2-m_1{J}_1{S}_2^2=-b^2\hat{S}-m_3{S}_1^2{J}_3-m_2{S}_1^2{J}_2-m_3{S}_2^2{J}_3-m_1{S}_2^2{J}_1-m_2{S}_3^2{J}_2-m_1{S}_3^2{J}_1=-b^2\hat{S}-(S_1^2+S_2^2+S_3^2)\mathcal{J}+m_1{S}_1^2{J}_1+m_2{S}_2^2{J}_2+m_3{S}_3^2{J}_3=-b^2\hat{S}+(S_1^2+S_2^2+S_3^2)\mathcal{J}+m_1{J}_1{S}^2+m_2{J}_2{S}^2+m_3{J}_3{S}^2=-b^2\hat{S}+(S_1^2+S_2^2+S_3^2)\mathcal{J}+\mathcal{J}{S}^2.$  Now we are ready to find the powers of $(S+\hat{S})$  . Using ( 21 ), we get $(S+\hat{S}{)}^2=S^2+S\hat{S}+\hat{S}S+{\hat{S}}^2=S^2+S\mathcal{J}+\mathcal{J}S+S_1^2+S_2^2+S_3^2-\hat{S}\mathcal{J}+\mathcal{J}S=S^2+S\mathcal{J}+2\mathcal{J}S-\hat{S}\mathcal{J}+S_1^2+S_2^2+S_3^2.$  Multiplying the result by $S+\hat{S}$  and applying again ( 21 ), we find $(S+\hat{S}{)}^3=(S+\hat{S}\left)\right[S^2+S\mathcal{J}+2\mathcal{J}S-\hat{S}\mathcal{J}+S_1^2+S_2^2+S_3^2]=S^3+S^2\mathcal{J}+2S\mathcal{J}S-S\hat{S}\mathcal{J}+S(S_1^2+S_2^2+S_3^2)+\hat{S}{S}^2+\hat{S}S\mathcal{J}+2\hat{S}\mathcal{J}S-{\hat{S}}^2\mathcal{J}+\hat{S}(S_1^2+S_2^2+S_3^2)=-b^{2}S+S^2\mathcal{J}-2m^{2}S+m^{2}S+S^2\mathcal{J}+\mathcal{J}{S}^2+\mathcal{J}S\mathcal{J}+2\hat{S}\mathcal{J}S-[(S_1^2+S_2^2+S_3^2)-\hat{S}\mathcal{J}+\mathcal{J}S]\mathcal{J}+[-b^2\hat{S}-(S_1^2+S_2^2+S_3^2)\mathcal{J}+\mathcal{J}{S}^2]=-(b^2+m^2)S+2S^2\mathcal{J}+\mathcal{J}{S}^2+\mathcal{J}S\mathcal{J}+2\mathcal{J}{S}^2-(S_1^2+S_2^2+S_3^2)\mathcal{J}+\hat{S}{\mathcal{J}}^2-\mathcal{J}S\mathcal{J}-b^2\hat{S}-(S_1^2+S_2^2+S_3^2)\mathcal{J}+\mathcal{J}{S}^2=-(b^2+m^2)S+2S^2\mathcal{J}+4\mathcal{J}{S}^2-2(S_1^2+S_2^2+S_3^2)\mathcal{J}-m^2\hat{S}-b^2\hat{S}=-(b^2+m^2\left)\right(S+\hat{S})+2S^2\mathcal{J}+4\mathcal{J}{S}^2-2(S_1^2+S_2^2+S_3^2)\mathcal{J}.$  Continue the process.

$(S+\hat{S}{)}^4=(S+\hat{S}\left)\right[-(b^2+m^2)(S+\hat{S})+2S^2\mathcal{J}+4\mathcal{J}{S}^2-2(S_1^2+S_2^2+S_3^2)\mathcal{J}]=-(b^2+m^2\left)\right(S+\hat{S}{)}^2+2S^3\mathcal{J}+4S\mathcal{J}{S}^2-2S(S_1^2+S_2^2+S_3^2)\mathcal{J}+2\hat{S}{S}^2\mathcal{J}+4\hat{S}\mathcal{J}{S}^2-2\hat{S}(S_1^2+S_2^2+S_3^2)\mathcal{J}=-(b^2+m^2)(S+\hat{S}{)}^2-2b^{2}S\mathcal{J}-4m^2{S}^2-2S^2\mathcal{J}\mathcal{J}+2\mathcal{J}{S}^2\mathcal{J}+4\mathcal{J}{S}^3-2[-b^2\hat{S}-(S_1^2+S_2^2+S_3^2)\mathcal{J}+\mathcal{J}{S}^2]\mathcal{J}=-(b^2+m^2\left)\right(S+\hat{S}{)}^2-2b^{2}S\mathcal{J}-4m^2{S}^2+2m^2{S}^2+2\mathcal{J}{S}^2\mathcal{J}-4b^2\mathcal{J}S+2b^2\hat{S}\mathcal{J}+2(S_1^2+S_2^2+S_3^2){\mathcal{J}}^2-2\mathcal{J}{S}^2\mathcal{J}=-(b^2+m^2)(S+\hat{S}{)}^2-2b^{2}S\mathcal{J}-2m^2{S}^2-4b^2\mathcal{J}S+2b^2\hat{S}\mathcal{J}-2m^2(S_1^2+S_2^2+S_3^2)=-(b^2+m^2\left)\right(S+\hat{S}{)}^2-2m^2[S^2+S\mathcal{J}+2\mathcal{J}S-\hat{S}\mathcal{J}+(S_1^2+S_2^2+S_3^2\left)\right]+(2m^2-2b^2)(S\mathcal{J}+2\mathcal{J}S-\hat{S}\mathcal{J})=-(b^2+3m^2)(S+\hat{S}{)}^2+(2m^2-2b^2\left)\right(S\mathcal{J}+2\mathcal{J}S-\hat{S}\mathcal{J}).$  Finally, $(S+\hat{S}{)}^5=(S+\hat{S}\left)\right[-(b^2+3m^2)(S+\hat{S}{)}^2+(2m^2-2b^2\left)\right(S\mathcal{J}+2\mathcal{J}S-\hat{S}\mathcal{J}\left)\right]=-(b^2+3m^2)(S+\hat{S}{)}^3+(2m^2-2b^2\left)\right[S^2\mathcal{J}+2S\mathcal{J}S-S\hat{S}\mathcal{J}+\hat{S}S\mathcal{J}+2\hat{S}\mathcal{J}S-{\hat{S}}^2\mathcal{J}]=-(b^2+3m^2\left)\right(S+\hat{S}{)}^3+(2m^2-2b^2)[S^2\mathcal{J}-2m^{2}S-S{\mathcal{J}}^2+\mathcal{J}S\mathcal{J}+2\mathcal{J}{S}^2-(S_1^2+S_2^2+S_3^2-\hat{S}\mathcal{J}+\mathcal{J}S)\mathcal{J}]=-(b^2+3m^2)(S+\hat{S}{)}^3+(2m^2-2b^2\left)\right[S^2\mathcal{J}-m^{2}S+\mathcal{J}S\mathcal{J}+2\mathcal{J}{S}^2-(S_1^2+S_2^2+S_3^2)\mathcal{J}+\hat{S}{\mathcal{J}}^2-\mathcal{J}S\mathcal{J}]=-(b^2+3m^2\left)\right(S+\hat{S}{)}^3+(m^2-b^2)[2S^2\mathcal{J}+4\mathcal{J}{S}^2-2(S_1^2+S_2^2+S_3^2)\mathcal{J}-2m^{2}S-2m^2\hat{S}]=-(b^2+3m^2)(S+\hat{S}{)}^3+(m^2-b^2\left)\right[2S^2\mathcal{J}+4\mathcal{J}{S}^2-2(S_1^2+S_2^2+S_3^2)\mathcal{J}-(m^2+b^2)(S+\hat{S})+(b^2-m^2)(S+\hat{S})]=-(b^2+3m^2\left)\right(S+\hat{S}{)}^3+(m^2-b^2)\left[\right(S+\hat{S}{)}^3+(b^2-m^2)(S+\hat{S})]=-2(b^2+m^2\left)\right(S+\hat{S}{)}^3-(b^2-m^2{)}^2(S+\hat{S})$  which completes the proof.

Thus,

$$\begin{array}{c}(S+\hat{S}{)}^5=Lin((S+\hat{S}{)}^3,S+\hat{S})\end{array}$$

(22)

On the other hand, setting for brevity $R_{XY}=R$  , from ( 19 ) we derive

$$\begin{array}{c}S+\hat{S}=\frac{4}{c}R-2\mathcal{J}=Lin(R,\mathcal{J})\end{array}$$

(23)

Since $(S+\hat{S})$  and $\mathcal{J}$  commute, ( 19 ) implies the commutation of $R$  and $\mathcal{J}$  .

Keeping this and ${\mathcal{J}}^2=-m^{2}E$  , from ( 23 ) we derive

$$\begin{array}{c}(S+\hat{S}{)}^3={\left(\frac{4}{c}\right)}^3{R}^3+Lin(\mathcal{J}{R}^2,R,\mathcal{J})\end{array}$$

(24)

$$\begin{array}{c}(S+\hat{S}{)}^5={\left(\frac{4}{c}\right)}^5{R}^5+Lin(\mathcal{J}{R}^4,R^3,\mathcal{J}{R}^2,R,\mathcal{J})\end{array}$$

(25)

From ( 22 ), ( 23 ) and ( 24 ) $$(S+\hat{S}{)}^5=Lin\left[Lin(R^3,\mathcal{J}{R}^2,R,\mathcal{J}),Lin(R,\mathcal{J})\right]=Lin(R^3,\mathcal{J}{R}^2,R,\mathcal{J}).$$  So, finally from ( 25 ) $$R^5=Lin(\mathcal{J}{R}^4,R^3,\mathcal{J}{R}^2,R,\mathcal{J}).$$  It is easy to trace that the coefficients of all linear combinations are polynomials in $\frac{1}{c},b,m$  . To complete the proof we should remark that $R^6=R^{5}R=Lin(\mathcal{J}{R}^4,R^3,\mathcal{J}{R}^2,R,\mathcal{J})R=Lin(\mathcal{J}{R}^5,R^4,\mathcal{J}{R}^3,R^2,\mathcal{J}R)=Lin\left[\mathcal{J}Lin(\mathcal{J}{R}^4,R^3,\mathcal{J}{R}^2,R,\mathcal{J}),R^4,\mathcal{J}{R}^3,R^2,\mathcal{J}R\right]=Lin(R^4,\mathcal{J}{R}^3,R^2,\mathcal{J}R,E)$  which allows to find all powers of $R$  inductively.

References

1. Azo K. A note on the projection curves of geodesics of the tangent and tangent sphere bundles, Math. Repts. Toyama Univ., 1988
2. Nagy P. Geodesics on the tangent sphere bundle of a Riemannian manifold, Geometria Didicata 7 (1978), 2, 233-244.
3. Sasaki S. Geodesics on the tangent sphere bundles over space forms, Journ. Reine Angew. Math. 288 (1976), 106-120.
4. Sato K. Geodesics on the tangent bundles over space forms, Tensor 32 (1978), 5-10.
5. Yampolsky A. Characterization of projections of geodesics of Sasaki metric of $TC{P}^n$  and $T_{1}C{P}^n$  , Ukr. Geom. Sbornik 34(1991), 121-126.