Applications of Nijenhuis geometry: Nondegenerate singular points of Poisson-Nijenhuis structures

We study and completely describe pairs of compatible Poisson structures near singular points of the recursion operator satisfying natural non-degeneracy condition.

The paper [2] suggests a programme and general strategy for studying such operators (and demonstrates that the programme is realistic by proving a series of non-trivial results). Important part of the strategy is the study of singular points of Nijenhuis operators.
The goal of this paper is to illustrate how methods, ideas and results of [2] work in describing singular points of more complicated geometric structures involving Nijenhuis operators as one of their ingredients (cf. [2, Sections 6.2 and 6.3]) and also to give a proof for Theorem 6.3 announced in [2].
The geometric structure we consider in this paper came from the theory of integrable systems, where Nijenhuis operators serve as recursion operators for bihamiltonian structures [7,12,14,15], see also [4,5].
Recall that two Poisson structures P and P are called compatible, if their sum P + P is also a Poisson structure. If P and P are non-degenerate and therefore come from symplectic forms ω = P −1 and ω = P −1 , then the compatibility condition is equivalent to the property that the recursion operator L given by the relation ω( · , · ) = ω( L · , ·) (1) is Nijenhuis. Since ω can be recovered from ω and L, we can reformulate the compatibility condition by saying that a symplectic structure ω and a Nijenhuis operator L are compatible if (I) the form ω( · , · ) = ω( L · , ·) is skew-symmetric, i.e., is a differential 2-form, (II) d ω = 0, i.e., this form is closed.
In the case when L is non-degenerate, a compatible pair (ω, L) defines a Poisson-Nijenhuis structure in the sense of [12,15]. However, non-degeneracy of L is not very essential as we can replace it with L + c · Id, c ∈ R and here we will think of compatible pairs (ω, L) as a natural subclass of Poisson-Nijenhuis structures.
It is natural to ask about local simultaneous canonical form for compatible ω and L. In the case when L is algebraically generic at a point p ∈ M 2n , i.e., its algebraic type (see the next section for details) remains the same in a certain neighbourhood U(p), the answer was obtained in [15], [18], [7]. We briefly recall some of these results.
Notice, first of all, that condition (I) imposes natural algebraic restrictions on the algebraic type of L: in the Jordan decomposition of L all the blocks can be partitioned into pairs of equal blocks (i.e. of the same size and with the same eigenvalue). In particular, the characteristic polynomial of L is a full square and each eigenvalue has even multiplicity.
If L has n = 1 2 dim M distinct real eigenvalues, each of multiplicity 2, and in addition their differentials are linearly independent at a point p ∈ M, then there exists a local symplectic coordinate system x 1 , . . . , x n , p 1 , . . . , p n in which L is diagonal with x i being its i-th eigenvalue [15]: ω = i dp i ∧ dx i , L = diag(x 1 , x 2 , . . . , x n , x 1 , x 2 , . . . , x n ).
Equivalently, one can say that the pair (ω, L) is the direct product of two-dimensional blocks For further use and without pretending that it is new, let us slightly generalise this result (cf. [15], [18], [7] Type 2. L has one real constant eigenvalue of multiplicity 2k; Type 3. L has one pair of non-constant complex conjugate eigenvalues of multiplicity 2: where z = x + i y, w = u + i v, J denotes the corresponding complex structure and α(z, w) + i β(z, w) is a holomorphic function in z and w.
Type 4. L has one pair of constant complex conjugate eigenvalues of multiplicity 2k: where z j = x j + i y j , w j = u j + i v j , J denotes the corresponding complex structure and α + i β ∈ C, β = 0.
Let ω and L be compatible. Suppose that L is semisimple and algebraically generic in a neighbourhood of p ∈ M. Then the pair (ω, L) locally splits into a direct product of 'elementary blocks' of 4 types described above. If dλ(p) = 0 or d(α(p) + i β(p)) = 0 for some real or complex eigenvalue, then in the corresponding 'elementary block' we may set λ = x (see Type 1) and α + i β = z (see Type 3).
Notice that to admit a compatible symplectic partner, a semisimple algebraically generic Nijenhuis operator L should satisfy the following additional condition: non-constant eigenvalues of L must be all of multiplicity two.
If L is algebraically generic but not necessarily semisimple, the description (rather nontrivial) of compatible pairs (ω, L) was obtained by Turiel [18] under some additional assumptions on the differentials of tr L k , k = 1, . . . , n. These assumptions basically mean that each eigenvalue is either constant or its differential does not vanish.
The goal of the present paper is to give a simultaneous canonical form for (ω, L) at typical singular points, where the algebraic type of L changes.
If we want to study singularities of Poisson-Nijenhuis structures (cf. Problem 5.17 in [1]), then it is natural to ask which Nijenhuis operators admit compatible symplectic structures and what is a simultaneous canonical form for ω and L near a (possibly singular) point p ∈ M?
Let us explain in more detail what kind of singularities we have in mind. For a symplectic structure ω, all points p ∈ M are obviously equivalent to each other (in other words, there are no singular points). However it is not the case for an operator L since its algebraic properties may essentially vary from point to point. More precisely, at each point p ∈ M, we can define the algebraic type (or Segre characteristic, see e.g. [6]) of the operator L(p) : T p M → T p M as the structure of its Jordan canonical form which is characterised by the sizes of Jordan blocks related to each eigenvalue λ i of L(p) (the specific values of λ i 's are not important here). Following [2, Definition 2.7], we call a point p ∈ M algebraically generic, if the algebraic type of L does not change in some neighbourhood U(p) ⊂ M. Otherwise, we say that p is singular.
Study of singularities of Nijenhuis and other related structures is an important part of the programme suggested in [2] (see also [11]). Actually, in most branches of geometry understanding of singular points of relevant geometric structures plays a central role.
Among singular points, we distinguish (see [2,Definition 2.10]) an important subclass of the so-called differentially non-degenerate singular points p ∈ M which are characterised by the property that the differentials of the functions tr L k , k = 1, . . . , dim M, are linearly independent at p. Instead of the traces of powers of L one can equivalently consider the coefficients σ 1 , . . . , σ m , m = dim M, of the characteristic polynomial of L: Relation (I), however, implies that the eigenvalues of L have always even multiplicity, which in turn implies that L cannot be differentially non-degenerate (since the characteristic polynomial is a full square and has at most 1 2 dim M functionally independent coefficients). As a natural analog of differential non-degeneracy in this setup we suggest the condition that the differentials of the functions tr L, tr L 2 ,. . . , tr L n , n = 1 2 dim M, are linearly independent at a given (singular) point p ∈ M. In terms of the characteristic polynomial of L which is a full square (of a polynomial with smooth coefficients) in our case: this condition amounts to linear independence of the differentials dh 1 (p), . . . , dh n (p). Our goal is to study local canonical forms for ω and L at such points. This problem is clearly important for finite dimensional integrable systems and we do hope that it might also be important for understanding of bihamiltonian structures in infinite dimension appearing in the theory of integrable ODEs and PDEs. The main result of this paper is the following theorem whose weaker version is announced in [2].
Theorem 2.1. Let ω and L be compatible (i.e., define a Poisson-Nijenhuis structure on M 2n ) and real analytic. Suppose that at a point p ∈ M 2n , the differentials d tr L, . . . , d tr L n are linearly independent. Then there exists a local coordinate system x 1 , . . . , x n , p 1 , . . . , p n such that ω = i dx i ∧ dp i and L is given by the matrix where 0 n is the zero n × n-matrix, and A ⊤ denotes the transposed of A.
Remark 2.1. We would like to point out that this theorem mainly concerns singular points of Poisson-Nijenhuis structures, i.e., those at which the algebraic type of L changes (recall that the case of algebraically generic points has basically been covered by the results of [15], [18], [7]). The typical and most important example of such a situation is a "perturbation of a Jordan block" when all the eigenvalues of L simultaneously vanish at the point p ∈ M 2n . In this situation, the x-coordinates of p are (0, 0, . . . , 0). However, the formulas (3) and (4) are universal in the sense that they are applicable to both algebraically generic and singular points of any algebraic kind under the only condition that the differentials d tr L, . . . , d tr L n are linearly independent.
More precisely, if a Nijenhuis-Poisson structure (L, ω) satisfies the conditions of Theorem 2.1 and in addition, at the point p ∈ M 2n , the operator L has s distinct eigenvalues λ 1 , . . . , λ s (some of which can be complex conjugate) of multiplicities 2k 1 , . . . , 2k s , then the local canonical form for L, ω is given by (3) and (4) where the x-coordinates x 0 1 , . . . , x 0 n of the point p are the coefficients of the polynomial that is, the square root of the characteristic polynomials χ L(p) .
The rest of the paper is devoted to the proof of Theorem 2.1. First we discuss one cohomological property of gl-regular operators serving as a key step in our proof and also interesting on its own. This property uses some existence and uniqueness results from general theory of PDEs that we clarify in Appendix A. The proof of Theorem 2.1 is given in Section 4. In Appendix B we formulate an alternative version of Theorem 2.1.
All the objects we are dealing with (functions, manifolds, maps) are assumed to be real analytic.
3 One cohomology-like property of gl-regular Nijenhuis operators.
Let V be a finite-dimensional vector space. We will say that an operator A : V → V is gl-regular, if the dimension of its conjugacy class O(A) = {P AP −1 : P ∈ GL(V )} is maximal (i.e., equals n 2 − n). Equivalently, this regularity condition means that in an appropriate basis the matrix of A takes the so-called companion form as in (9). Typical examples of gl-regular operators include a semisimple operator with distinct eigenvalues and a single Jordan block.
To clarify the relation between the gl-regularity condition and our main Theorem 2.1, it is worth mentioning that the linear independence of the differentials d tr A(p), . . . , d tr A n (p) automatically implies that A(p) is gl-regular 1 . Notice that A and L from Theorem 2.1 are related by (3), hence tr L k = 2 tr A k so that the assumptions of this theorem imply gl-regularity of the block A. These two operators A and L are related via the Turiel extension construction discussed in Appendix B. However in this section, the reader should think of A just as an operator defined on an n-dimensional manifold Q and not related to L in any sense. Proposition 3.1. Let A be a (real-analytic) Nijenhuis operator which is gl-regular at a point p ∈ Q. Assume that Ω is a closed (real-analytic) 2-form such that the form is also closed. Then locally in some neighbourhood of p ∈ Q there exists a (real-analytic) function U such that Remark 3.1. Note that the gl-regularity condition is essential. It can be easily checked that the statement of Proposition 3.1 fails for the operator A = where the eigenvalues collide, i.e., Remark 3.2. Proposition 3.1 can be interpreted in terms of the Nijenhuis differential [12], [16]. Denote by Λ k (Q) the module of differential k−forms on a manifold Q, Λ 0 are functions on Q (definition works for both smooth and analytic cases) and d : Λ k → Λ k+1 is standard differential of k−forms.
Given an operator A, one can define a new differential d A : Λ k → Λ k+1 which is completely determined by the following properties (see [16,Appendix] If N A = 0, we have that d 2 A = 0 (see [12, Corollary 6.1] for E = T Q, N = A and µ being standard commutator of vector fields).
The goal of the rest of this section is to prove Proposition 3.1. Recall that we are mainly interested in those points p ∈ Q where A changes its algebraic type. The key example of such a situation is a perturbation of a Jordan block, i.e., A(p) is similar to a single Jordan block with zero eigenvalue, but generically the operator A(q), q ∈ V (p) is semisimple with distinct eigenvalues (as in (4)).
We however first consider the case when A(p) is regular and semisimple. Lemma 3.1. In the assumptions of Proposition 3.1, let A(p) be semisimple with distinct eigenvalues. Then the statement of Proposition 3.1 holds. Moreover, the general solution U of (6) depends on (is parametrised by) n arbitrary functions of one variable.
Proof. Consider a local coordinate system y 1 , . . . , y n in which A takes diagonal form Since A is Nijenhuis, such a coordinate system exists [9,17].
We need to solve (6) which takes, in coordinates y 1 , . . . , y n , the following form: or, equivalently,

This equation is consistent if and only if
which can also be rewritten as These two linear relations follow from (and in fact, are equivalent to) Hence, (8) (or equivalently (6)), admits a solution U 0 . Since (8) is linear in U, its general solution is defined modulo solutions of the corresponding homogeneous system that is, U = U 0 (y 1 , . . . , y n ) + u 1 (y 1 ) + · · · + u n (y n ).
In other words, U depends on n functions in one variables, as stated.
Remark 3.3. On can similarly check that Proposition 3.1 holds under the additional condition that p ∈ Q is algebraically generic, i.e. the algebraic type of A remains locally constant in a neighbourhood of p ∈ Q. The splitting theorem [2, Section 3.2] allows one to assume without loss of generality that A has either a single real eigenvalue or single pair of complex conjugate eigenvalues. In this case, due to gl-regularity, A will be similar to a single Jordan block (real or complex) and the general solution of (6) can be found explicitly by using the canonical forms for such operators from [2]. We omit details of this computation here, although the explicit description of functions U satisfying the condition d (A * d U) = 0 for A being a Nijenhuis Jordan block, could be interesting for several reasons. We will provide such a description in another paper.
The following statement is a simple fact from Linear Algebra. It is well known that every gl-regular operator A can be reduced (at one single point p ∈ Q) to companion form 2 , i.e., where σ i will be automatically the coefficients of the characteristic polynomial χ A (t) = det(t · Id − A).
Lemma 3.2. Let x 1 , . . . , x n be a local coordinate system in a neighbourhood of p ∈ Q such that in these coordinates A(p) = A comp . Then PDE system (6) can be resolved w.r.t. all second order derivatives of U except those of the form U x i xn (i = 1, . . . , n), i.e., can be rewritten in the form . , x n , U x 1 , . . . , U xn , U x 1 xn , . . . , U xnxn ), i, j = 1, . . . , n − 1. (10) Proof. At the point p ∈ Q itself, the proof is a simple exercise in Linear Algebra. In all neighbouring points, representation (10) for system (6) holds by continuity.
In order to prove Proposition 3.1 we will show that PDE system (10) possesses sufficiently many solutions (namely, the space of solutions is parametrised by n functions in one variable).
(ii) The following compatibility conditions identically hold in virtue of (10): Proof. The proof will be given in Appendix A.
In our case, the functions h ij are linear (but not homogeneous) in the derivatives U x i and U x i xn , i.e., For the given operator A and 2-form Ω, the compatibility conditions (11) take the form where c α are some functions in coordinates (x 1 , . . . , x n ) that can be explicitly written in terms of the components of the operator A and 2-form Ω and their partial derivatives, i.e., Recall that both A and Ω are rather special and satisfy very strong conditions. However, straightforward verification of compatibility conditions (12) (and even explicit description of them in terms of A, Ω and their derivatives) is not an easy task. Instead, we will use the fact that c α (x) is real-analytic in x and for this reason it is sufficient to verify these conditions only for some open subset V 0 ⊂ V (p). As such a subset we consider a small neighbourhood U(q) ⊂ V (p) of a point q ∈ V (p) at which A(q) is semisimple (similar verification can be done for q being algebraically generic 3 , cf. Remark 3.3). According to Lemma 3.1, the general solution U of PDE system (6) (which is equivalent to (10)) depends on n arbitrary real analytic functions of one variable. Basically, this property is equivalent to condition (i) from Lemma 3.3 and therefore the compatibility conditions (11) from (ii) hold.
However, the collections of n functions parametrising the space of solutions of (6) in Lemma 3.1 and in Lemma 3.3 are formally different. So we need to justify the fact that by choosing appropriate functions u 1 (y 1 ), . . . , u n (y n ) (see Lemma 3.1) we can fulfil arbitrary initial conditions v(x n ), v 1 (x n ), . . . , v n−1 (x n ) from Lemma 3.3.
It can be easily seen that each of these functions φ i is locally invertible, i.e., dφ i dt = 0 at t = 0. Indeed, this property means that the entries of the last column ξ of the Jacobi matrix J = ∂y ∂x (at the origin) are all different from zero. This matrix satisfies the relation A comp (x) = J −1 A diag (y)J or, equivalently, JA comp = A diag J. The latter can be interpreted as follows: the columns of J are (in the reverse order) of the form ξ, A diag ξ, A 2 diag ξ, . . . A n−1 diag ξ.
In particular, these vectors are linearly independent. For a diagonal matrix A diag this condition holds if and only if all the components of ξ ∈ R n are different from zero, as stated.
The initial conditions (i) from Lemma 3.3 can be rewritten in the form It is more convenient to rewrite the first condition in the form U xn (γ(t)) = v n (t) = v ′ (t), U(γ(0)) = c (where c ∈ R is an arbitrary constant). In other words, we replace v(t) with its derivative v ′ (t) and the value c ∈ R at t = 0.
Using the standard relations U x j = α ∂yα ∂x j U yα we can rewrite these conditions in coordinates y 1 , . . . , y n as follows: and using U(y 1 , . . . , y n ) = U 0 (y 1 , . . . , y n ) + u 1 (y 1 ) + · · · + u n (y n ): Taking into account the fact that φ α (t) is locally invertible, we see that the derivatives of the functions u 1 , . . . , u n can be uniquely reconstructed from v 1 , . . . , v n and vice versa as required.
We have shown that in a neighbourhood of any point q at which A(q) is diagonalisable over R, we can find a solution of (6) with arbitrary initial conditions (i) from Lemma 3.3. This guarantees the fulfilment of the stong compatibility conditions (11) from (ii) in this neighbourhood and therefore (due to real analyticity) in the whole neighbourhood V (p). The latter, in turn, guarantees the existence of local solutions of (6) in a neighbourhood of p ∈ Q, completing the proof of Proposition 3.1.

Proof of the main theorem
Recall that we consider ω and L to be compatible and real analytic. We also assume that at a point p ∈ M 2n all the eigenvalues of L vanish, but the differentials d tr L, . . . , d tr L n are linearly independent.
Since the bilinear form ω(·, ·) = ω(L·, ·) is skew-symmetric, the characteristic polynomial of L is a full square: i.e., the square of a certain polynomial of degree n. The condition that the differentials of the functions tr L, . . . , tr L n are linearly independent implies that the differentials of the functions h 1 , . . . , h n are linearly independent too and we can take them as the first n coordinates x 1 , . . . , x n of a local coordinate system. It is well known that these functions commute with respect to the Poisson bracket related to ω (see [15,7]) which, by Darboux theorem, implies local existence of functions p 1 , . . . , p n such that (x 1 , . . . , x n , p 1 , . . . , p n ) is a canonical coordinate system for ω.
In these coordinates, the operator L takes an "almost required" form. To see this we use the following formula proved in [2]. Let σ 1 , . . . , σ m , m = dim M, be the coefficients of the characteristic polynomial of L: considered as smooth functions on M. Then in any coordinate system y 1 , . . . , y m , the matrix L(y) satisfies the following relation and J = ∂σ i ∂y j is the Jacobi matrix of the functions σ 1 , . . . , σ m w.r.t. y 1 , . . . , y m . In our case, the characteristic polynomial χ L (t) is a full square as in (13) so that its coefficients σ i (i = 1, . . . , 2n) are some explicit polynomials in h α = x α (α = 1, . . . , n) and do not depend on the other half of coordinates p α (α = 1, . . . , n). It is now straightforward to check that in coordinates (x, p) the first n rows of the operator L are as in (3,4). Furthermore, the algebraic compatibility condition (I) implies that the lower right n × n-block of L is transposed to the upper left n × n-block, so that L takes the form with A as in (4) and S being just a skew-symmetric matrix (again due to condition (I)) whose components may a priori depend on all variables.
Next we use the differential compatibility condition (II) saying that the 2-form ω(·, ·) = ω(L·, ·) is closed. This form is given by and its differential is Substituting S = S + T in the formula above, where S is as in (4), we observe that the first sum in (16) cancels with the terms coming from S and we obtain Since ω is closed, we see that T ij does not depend on the p-variables, so that the form n i,j=1 T ij (x)dx i ∧ dx j can be viewed as a closed 2-form T on the local n-dimensional coordinate chart (x 1 , . . . , x n ).
To complete the proof, we will show that the form T can be "killed" by a suitable canonical transformation of the form (X 1 , . . . , X n , P 1 , . . . , P n ) = x 1 , . . . , x n , p 1 + ∂U ∂x 1 , . . . , p n + ∂U ∂xn , where U is a function of x 1 , . . . , x n (called sometimes generating function). This transformation preserves the (canonical) symplectic structure ω so that we only need to look after the change of L.
The Jacobi matrix of coordinate transformation (18) is J = id 0 d 2 U id , where id denotes the identity n × n matrix and d 2 U = ∂ 2 U ∂x i ∂x j is the Hessian matrix of U. Hence, after transformation (18) the matrix of L given by (14) takes the form Taking into account the transformation of the p-coordinates in S, we see that L new takes the required form (3,4) if and only if the generating function U satisfies the following equation: which can be equivalently rewritten in a more invariant form: that coincides with the PDE system (6) treated in the previous section. In order to check the existence of a generating function U solving (20), it remains to verify two conditions from Proposition 3.1 imposed on the 2-form T .
They both follow from the fact that L = A 0 S A ⊤ with S = S + T and A and S as in (4), is a Nijenhuis operator. It is easily checked by a straightforward computation that vanishing Nijenhuis torsion of L is exactly equivalent to the two geometric conditions we need, namely that the 2-forms are closed. This allows us to apply Proposition 3.1 that now guarantees the existence of U solving (20) and hence completes the proof of Theorem 2.1.
Lemma 3.3 is not a new fact and can be understood as a special case of the Cartan-Kähler theorem (see e.g. [8]). We suppose that people working in the Cartan-Kähler theory would find it trivial. However, deriving Lemma 3.3 from a very general and non-trivial Cartan-Kähler theorem seems to require more work than proving it directly. Unfortunately, we have not found this lemma in the required form in the literature, so we will prove it here in order to make the paper self-contained.
As one usually does in the theory of differential equations, we set . . , f n = U xn and rewrite the above system in the form In our case the number of unknown functions f j equals the number of variables x i , i.e., k = n but in general this assumption is irrelevant.
In this view, we will prove an analogue of Lemma 3.3 for a more general PDE system (22) (from which Lemma 3.3 can be easily derived by using the above substitution f i = U x i ). Recall that all the functions involved are assumed to be real analytic.
Lemma 5.1. The following properties of (22) are equivalent: • for any initial straight line x 1 = a 1 , . . . , x n = a n and any real analytic initial conditions on it f 1 (a 1 , . . . , a n−1 , x n ) = v 1 (x n ) f 2 (a 1 , . . . , a n−1 , x n ) = v 2 (x n ) . . . f k (a 1 , . . . , a n−1 , x n ) = v k (x n ) there locally exists a unique real analytic solution F = • The following compatibility conditions are identically fulfilled in virtue of (22): Proof. We give a proof in the simplest non-trivial case 4 , namely, for three variables x, y, z and one function f so that (22) becomes f x = h 1 (x, y, z, f, f z ) and f y = h 2 (x, y, z, f, f z ).
Assume that the condition D x h 2 = D y h 1 holds identically in virtue of the system. In more details: And similarly In these two expressions we see an additional formal variable f zz . In this simples case the corresponding coefficients in the both take the form ∂h 1 ∂fz ∂h 2 ∂fz and therefore coincide automatically. However in general it is not necessarily the case so that the equality these coefficients should be considered as an extra condition additional to The identical fulfilment of this condition means that the functions in the right hand side and the left hand side coincide as functions of 5 independent variables x, y, z, f and f z .
Without loss of generality, we will assume that the parameters a 1 and a 2 in the initial condition vanish. For given initial condition f (0, 0, z) = v(z) we construct the desired solution in two steps. We first construct a unique function of two variables u(y, z) that is a local realanalytic solution of the Cauchy problem u y = h 2 (0, y, z, u, u z ), u(0, z) = v(z).
The existence of u(y, z) follows from the classical Cauchy-Kovalevskaya theorem. If the desired solution f of (23) exists then necessarily f (x, 0, z) = u(x, z).

The next step is solving the Cauchy problem
Again, a unique local real-analytic solution f (x, y, z) exists by Cauchy-Kovalevskaya theorem.
Thus we have constructed a function that satisfies the initial condition and solves the first equation of (23), as for the the second equation, it is fulfilled only on the plane x = 0. Our goal is to show that D x h 2 = D y h 2 guarantees that the second equation holds at all the points. Then if the compatibility condition D x h 2 = D y h 1 is identically fulfilled, then f y (x, y, z) = h 2 (x, y, z, f, f z ) for all x, y, z. where By [19], this operator is Nijenhuis and compatible with the canonical symplectic structure i dp i ∧ dx i on T * Q. Example 6.1. By direct calculations we see that the "normal form" from Theorem 2.1 is in fact the Turiel extension of the Nijenhuis operator in the (first) companion form (9) with σ i = x i : indeed, S and A from (4) are related by (25).
Observe that the Turiel construction is evidently geometric and does not depend on the choice of a coordinate system x 1 , . . . , x n in which we apply it. Indeed, L is the recursion operators for two 2-forms ω and ω each of which is geometric (see [10]): the first of them is canonical ω = −dθ, θ = p α dx α , the second is defined by the following invariant formula ω = −d θ with θ = A * θ = A i j (x)p i dx j . One therefore may look for coordinate systems on Q such that in the corresponding coordinate system on T * Q the extended operator L takes a more convenient form. The interpretation of "convenient" may depend on the problem we are solving and on the situation in which the recursion Nijenhuis operator appeared. As an illustration, let us show that one can choose a coordinate system on Q in such a way that S in (24) vanishes.
We say that a Nijenhuis operator A is in the second companion form if its matrix, in local coordinates y 1 , . . . , y n , has the form (cf. (9)) Not every operator of the form (26) is Nijenhuis. The next lemma gives necessary and sufficient condition for this property to hold.
Condition (27) simply means that the differential 1-form given by the last row of A is closed, i.e., ∂σ n−i+1 ∂y j = ∂σ n−j+1 ∂y i for any i, j; it is a linear condition on σ i . The other condition (28) is nonlinear.
If A is Nijenhuis and both α and A * α are closed, this formula implies that also A * 2 α is closed. For our A given by (26), the forms dy 1 and dy 2 = A * dy 1 are closed. Thus, if A is Nijenhuis, also the form (A * ) n dy 1 = A * dy n is closed, which gives us (27). Similarly, the form (A * ) n+1 dy 1 = A * 2 dy n is closed, which gives us (28).
In order to prove Lemma 6.1 in the other direction, let us observe that since N A is a tensor, it is sufficient to check that (29) vanishes for the basis forms α = dy i , i = 1, . . . , n. For them each term in (29) is zero.
Example 6.2. The Turiel extension of the Nijenhuis operator A in the second companion form is where 0 n is the zero n × n matrix. Indeed, for fixed i the 1-form A α i dx α is just A * (i−1) dx 1 ; so the corresponding sum in (25) vanishes for i < n trivially and for i = n by (27).
We use this observation to give another canonical form for Poisson-Nijenhuis structures from Theorem 2.1. To that end we use a coordinate transformation x → y that reduces the Nijenhuis operator A from Theorem 2.1 (being in the first companion form): to the second companion form (26). Suitable new coordinates (y 1 , . . . , y n ) are easy to find by using [2, Proposition 2.1] which states that for any Nijenhuis operator A and any k ∈ N we have Let us set: y 1 = − tr A, y 2 = − 1 2 tr A 2 , . . . , y n = − 1 n tr A n .
Notice that the transition formulas x → y express the traces y 1 , . . . , y n of powers of A (with some coefficients) in terms of the coefficients of its characteristic polynomial x 1 , . . . , x n . Such a transformation is known to be invertible and polynomial in both directions (see below the formulas for the inverse transformation y → x).
From (32) we see that A * dy i = dy i+1 (for i ∈ 1, . . . , n − 1) meaning that the first n − 1 rows of A in the y-coordinate system are as we want, i.e., in these coordinates A takes the second companion form (26). The corresponding functions σ k (y) in the last row of (26) can be explicitly constructed. Indeed, σ k are the coefficients of the characteristic polynomial χ A (t) = det(t · Id −A) = t n + σ 1 t n−1 + · · · + σ n so that the functions σ k (y 1 , . . . , y n ) coincides with x k and are given by the condition σ k − tr A, − 1 2 tr A 2 , . . . , − 1 n tr A n = coefficient of t n−k in χ A (t).
This condition is well known and the functions σ i satisfying (33) can be obtained from the famous Newton-Girard polynomials (by appropriate rescaling and changing sings). They depend neither on the matrix A nor on the dimension and are given by the following universal formula, see e.g.
Theorem 6.1. Let L and ω be compatible (i.e., define a Poisson-Nijenhuis structure on M 2n ) and real analytic. Suppose that at a point p ∈ M 2n the differentials d tr L, . . . , d tr L n are linearly independent. Then there exists a local coordinate system y 1 , . . . , y n , p 1 , . . . , p n such that ω = i dy i ∧ dp i and L is given by the matrix where A is given by (26) with σ k given by (34).