the integrability problem for systems of first order odes is indeed a relevant issue in the theory of dynamical systems ( see , for instance ,  @xcite and references therein ) , and the same question for coupled systems of odes could be also meaningful from the viewpoint of applications ( for instance , regarding synchronization problems  @xcite ) . in this context , the aim of this paper is two fold : on one hand , to present a novel approach to the explicit construction of liouville  integrable coupled systems of odes endowed with lie  poisson symmetries and , on the other hand , to show how integrable deformations of the previous systems can be systematically obtained .
the construction here presented establishes a ( to the best of our knowledge ) novel connection between poisson
lie groups and integrable deformations of dynamical systems .
such a link is based on the well - known result by drinfeld that establishes the one - to - one correspondence between poisson  lie groups and lie bialgebras  @xcite , as well as on the general construction of integrable hamiltonian systems from poisson coalgebras that was introduced in  @xcite , thus giving rise to the so - called  coalgebra symmetry approach " to finite dimensional integrable systems ( see  @xcite for applications of this method to hnon  heiles and lotka  volterra systems , respectively , and  @xcite for the coalgebra  based construction of nonlinear superposition rules arising in nonautonomous lie
hamilton systems  @xcite ) .
in particular , we will apply the formalism here introduced in order to obtain new integrable deformations of certain lorenz  @xcite and rssler  @xcite systems , as well as of some coupled versions of them .
nevertheless , we stress that the method presented in the paper is completely general and could be applied to any other interesting hamiltonian systems of odes endowed with a lie  poisson symmetry .
the basics of the poisson coalgebra approach to integrability will be summarized in the next section . in short
, we will consider a finite  dimensional integrable hamiltonian system of odes that is defined through a lie  poisson algebra @xmath3 , where @xmath4 is the algebra of smooth functions on @xmath5 and the lie
poisson bracket is @xmath6 .
a set of functions in involution will be given by a single hamiltonian function @xmath7 and a number of casimir functions @xmath8
, @xmath9 . then , since lie
poisson algebras are endowed with a ( primitive or non  deformed ) coalgebra structure , the coalgebra symmetry approach will straightforwardly provide , by following  @xcite , a coupled integrable generalization of the system .
this approach will be illustrated in section 3 through the explicit construction of an integrable coupled rssler system .    a detailed technical presentation of the construction of poisson
lie groups as deformations of lie bialgebras will be given in section 4 .
if we denote by @xmath10 the abstract lie algebra associated with the lie
poisson coalgebra @xmath3 , we will assume that we are able to find a cocommutator map @xmath11 compatible with @xmath10 , so that @xmath12 defines a lie bialgebra structure  @xcite .
then , it is well  known that the dual of the cocommutator @xmath11 will define a second lie algebra structure ( that we will denote by @xmath13 ) , and the dual of the structure tensor for @xmath14 will provide a compatible cocommutator map @xmath15 for @xmath13 , in such a way that @xmath16 is just the dual lie bialgebra to @xmath12 .
once a given @xmath17 is identified , the ( connected component containing the identity ) of the lie group @xmath18 with lie algebra @xmath13 can be constructed by exponentiation , and the unique poisson  lie structure having @xmath19 as its linearization at the identity can be computed , thus providing the deformed poisson structure @xmath20 we were looking for .
simultaneously , the group law for @xmath18 will give us the explicit deformed coproduct map @xmath21 on the algebra of smooth functions @xmath22 .
moreover , the deformation nature of this construction arises in a more transparent way if we consider the one - parametric cocommutator @xmath23 , since the parameter @xmath24 will give us information concerning the orders in the deformation process that arises from the ` exponentiation '  of the tangent lie bialgebra structure .
it is worth to be stressed that the procedure that we have just sketched is algorithmic and can be made fully explicit by making use of symbolic computation packages .
therefore , as a poisson coalgebra , @xmath25 can be interpreted as a smooth @xmath24deformation of @xmath26 on which the initial dynamical system with hamiltonian function @xmath27 is defined . as a consequence ,
if we consider the same hamiltonian function @xmath27 on @xmath25 , then the associated dynamical system is just an integrable @xmath24deformation of the initial system of odes .
moreover , the coalgebra symmetry approach  @xcite guarantees that the set of involutive functions for the deformed system will be given by the deformed coproducts of the hamiltonian @xmath27 and of the deformed casimir functions @xmath28 for the poisson structure @xmath20 , with @xmath29 .
sections 5 and 6 are devoted to use this poisson - deformation approach in order to obtain integrable deformations of two coupled rssler and lorenz systems , respectively .
it is worth stressing that , as a result of the exponentiation process that underlies the full construction , such integrable deformations of dynamical systems are provided by non
polynomial vector fields , which is a non
standard finding from the viewpoint of integrability . also , the underlying poisson  lie groups are responsible of the fact that the coupling between the systems is quite involved , since such coupling is generated by their non  abelian group law .
a final section including some remarks and open problems closes the paper .
let @xmath14 be a finite dimensional lie algebra . as it is well known , the associated lie
poisson algebra @xmath4 is the algebra of @xmath30 functions on @xmath5 , equipped with the poisson bracket @xmath31,x \rangle , \qquad x \in g^ * , \quad f , g \in { \cal f}(g^ * ) , \label{lpintrinsic}\ ] ] where @xmath32 is the non - degenerate pairing between @xmath14 and @xmath5 .
let us assume that @xmath14 is generated by the set @xmath33 , with commutation relations @xmath34=c_{ij}^k\ , x_k , \label{structurecoeffs}\ ] ] and let us denote with @xmath35 the corresponding dual basis in @xmath19 @xmath36 then the lie
poisson bracket ( [ lpintrinsic ] ) for the coordinate functions on @xmath5 is given by the fundamental poisson brackets @xmath37    now , by defining the following coproduct map @xmath38 for the coordinate functions @xmath39 we can say that @xmath26 is endowed with a poisson coalgebra structure  @xcite , since the coassociativity condition @xmath40 is fulfilled , and the coproduct map @xmath41 can be straightforwardly shown to define a poisson map between @xmath42 and @xmath43 , namely @xmath44 with respect to the poisson structure on @xmath45 induced from  .
the coalgebra approach to integrability is based on the fact that if we start from a liouville integrable system defined on @xmath42 , we can use the coproduct map ( [ primitive ] ) in order to extend it to an integrable system defined on @xmath46^n:= { \cal{f}}(g^ * ) \otimes   \dots^{n}\otimes { \cal{f}}(g^*)$ ] , with @xmath47 being an arbitrary number of copies of @xmath4 .
the keystone for this result is given by the @xmath48-th coproduct maps , which are defined recursively as @xmath49 and , by construction , satisfy the poisson homorphism property ( [ homorphism ] ) @xmath50    the simplest possible case in this framework occurs when the integrable system on @xmath42 is defined by a single hamiltonian function @xmath27 and an arbitrary number of casimir functions @xmath51 , @xmath9 .
in this case all the involutive functions defining an integrable system on @xmath46^n$ ] are obtained as images of the hamiltonian and casimir functions under the @xmath52 maps @xmath53 . in this paper
we will consider only this case , but we recall that this approach can be generalized for poisson algebras in which integrability involves either more than a single hamiltonian or even non - coassociative poisson maps ( see  @xcite ) .
in order to illustrate this integrability approach , let us consider the specific lie  poisson algebra that underlies the integrability of a dynamical system of rssler type .
we recall that the so  called rssler system @xmath54 was introduced in @xcite as a prototype of 3d system of odes exhibiting the phenomenon of continuous chaos .
when @xmath55 , the system ( [ s1 ] ) is known to be integrable in the liouville sense @xcite , and gives rise to @xmath56 in particular , this system is known to be hamiltonian  @xcite with respect to the poisson bracket @xmath57 provided that we consider the hamiltonian function @xmath58 the bracket   can be thought of as a particular case of the lie
poisson bracket @xmath59 where we have introduced a fourth  central " poisson algebra generator @xmath60 .
it is immediate to show that this is just the lie
poisson algebra corresponding to the 4d real solvable lie algebra @xmath61 ( here we follow the classification of four  dimensional real lie algebras given in  @xcite ) .
the casimir functions for ( [ la ] ) are given by @xmath62 indeed , if the hamiltonian of the system is again given by ( [ rh ] ) , the rssler system ( [ rs ] ) is recovered when @xmath63 .
poisson coalgebra @xmath64 is defined through the primitive coproduct map @xmath65 where in the second equality we have identified the coordinate functions on each copy of @xmath42 through a different subscript .
now , by using the standard poisson coalgebra structure   , we are able to define integrable coupled generalizations of the rssler system ( examples of nonintegrable couplings can be found in  @xcite ) .
explicitly , if we consider the coproduct of the hamiltonian  , we get @xmath66^{2}+[\delta(y)]^{2}\right)+\delta(z)= \dfrac{1}{2}(x_1^{2}+y_1^{2})+z_1+\dfrac{1}{2}(x_2^{2}+y_2^{2})+z_2+x_1x_2+y_1y_2,\ ] ] and the corresponding 6d equations of motion lead to a coupled rssler system @xmath67 \dot{x}_{2}=-w_{2}(y_{1}+y_{2})-z_{2},\\ \dot{y}_{2}=w_{2}(x_{1}+x_{2}),\\ \dot{z}_{2}=z_{2}(x_{1}+x_{2 } ) , \end{array } \label{coro}\ ] ] together with @xmath68 .
notice that the coupling between both systems comes from the non  linear terms of the hamiltonian  , and the restriction of the dynamics to the submanifold @xmath69 ( or , alternatively , to @xmath70 ) leads us to the uncoupled rssler system ( [ rs ] ) .
the 6d hamiltonian system   is integrable since , due to the underlying coalgebra symmetry , we have three additional ( besides the hamiltonian ) independent integrals of the motion that are in involution . two of them are given by the ( nontrivial ) casimir functions on each of the copies of the lie
poisson algebra , namely _
1=z_1 e^- , _ 2=z_2 e^- , and the the third one is just the coproduct of the casimir @xmath71 , which reads @xmath72    the @xmath47dimensional generalization of this result is straightforward by considering the @xmath47-th coproduct map @xmath73 of the hamiltonian and of the casimir functions , that is obtained as the pull - back through the @xmath47-th coproduct of the coordinate functions  @xcite & & ^(n)(x)=x11  ^n-1)1 & & + 1x1
^n-2 ) 1 +  & & + 11  ^n-1)1x , [ fo ] ( and the same for @xmath74 and @xmath75 ) .
explicitly , @xmath76^{2}+[\delta^{(n)}(y)]^{2}\right)+\delta^{(n)}(z)=
\dfrac{1}{2}\left(\sum_{i=1}^n x_i\right)^{2 } + \dfrac{1}{2}\left(\sum_{i=1}^n y_i \right)^{2}+\sum_{i=1}^n z_i , \label{nr}\ ] ] and the @xmath47-coupled rssler system reads @xmath77 this system of @xmath78 coupled odes has the following ( @xmath79 ) integrals of the motion : _
i = z_i e^- , i=1,
 ,n @xmath80 by construction , all these functions are in involution and poisson - commute with the hamiltonian @xmath81 .
the aim of this paper is to construct integrable deformations of coalgebra symmetric systems of the type   for any number @xmath82 of coupled constituents .
as we will show in the sequel , this can be achieved by considering the lie
poisson algebra @xmath19 as the linearization of a poisson  lie group structure on a certain _ different _ lie group @xmath18 . in his way , the construction of the corresponding full poisson  lie group structure @xmath25 will provide a poisson coalgebra @xmath24deformation of @xmath64 , and the application of the coalgebra approach to a hamiltonian defined on the poisson
lie group @xmath18 will automatically provide an integrable @xmath24deformation of the initial system defined on @xmath42 .      given a lie algebra @xmath14 with basis @xmath33 and lie bracket @xmath34=c_{ij}^k\ , x_k , \label{structurecoeffs}\ ] ] a lie bialgebra structure @xmath83
is given by a skew - symmetric cocommutator map @xmath84 @xmath85 such that    * \i ) @xmath11 is a 1-cocycle , _
i.e. _ , ( [ x , y])=[(x ) , y1 + 1y ] + [ x1 + 1x , ( y ) ] , x , y . [ 1cocycle ] * \ii ) the dual map @xmath86 is a lie bracket .
therefore , the dual @xmath87 of the cocommutator map defines a second ( in general , different ) lie algebra structure , that we will call @xmath13 , and whose lie brackets are @xmath88= f^{ij}_k \ , \hat x^k .
\label{geta}\ ] ] here the duality relation is established through the canonical pairing @xmath89 .
note that , given the initial lie algebra @xmath14 , the dual lie algebra @xmath13 has to fulfill the 1-cocycle condition , which implies the following compatibility equations between the two lie algebras @xmath14 and @xmath13 : f^ab_k c^k_ij = f^ak_i c^b_kj+f^kb_i c^a_kj + f^ak_j c^b_ik + f^kb_j c^a_ik .
[ compatfc ]    by taking the components of the structure tensor @xmath90 as initial data , the solutions to the ( now linear ) equations   together with the ( quadratic ) equations coming from the jacobi identity for @xmath91 provide all possible lie bialgebra structures for @xmath14 .
finally , if the solutions for the tensor @xmath92 are classified into equivalence classes under the action of all possible automorphisms for @xmath14 , the classification of all possible lie bialgebra structures of @xmath14 is obtained .
two remarks are in order :    * it is immediate to prove that if @xmath83 is a lie bialgebra with structure tensors @xmath93 , then @xmath94 with @xmath95 being the cocommutator map @xmath96 defines a lie bialgebra structure on @xmath13 with structure tensors @xmath97 , that is called the dual lie bialgebra to @xmath83 . in other words , @xmath13 and @xmath11
are defined through the same structure tensor @xmath91 , and we can say that @xmath98 . in the same manner , @xmath14 and @xmath15
are defined by the same structure tensor @xmath90 and we can say that @xmath99 . * for our purposes it will be convenient to multiply all the components @xmath92 of the cocommutator map by a real parameter @xmath24 , that will play the role of a deformation parameter . the resulting cocommutator will be denoted by @xmath100 , and the dual lie algebra structure @xmath101 will be @xmath88=\eta f^{ij}_k \,\hat x^k . \label{geta}\ ] ] obviously , @xmath13 is just @xmath102 and these two lie algebras are isomorphic under the change of basis @xmath103 .
nevertheless , as we will see in the sequel , the explicit appearance of @xmath24 will allow us to control the full deformation / exponentiation procedure .
a poisson  lie ( pl ) group is a lie group @xmath18 endowed with a poisson structure on the algebra @xmath104 of functions on the group manifold , such that the lie group multiplication is a poisson map . therefore , by defining the coproduct @xmath105 as the pull back of a function @xmath106 by the group multiplication , we obtain a natural structure of poisson coalgebra on the algebra of smooth functions on the poisson
lie group @xmath22 . in algebraic terms , if @xmath18 is a lie group , the group multiplication on @xmath18 induces a coproduct map : c^(d)c^(d)c^(d ) , which is coassociative ( since the group law is associative ) .
therefore , pl groups are instances of poisson coalgebras , since the coproduct map is -by definition- a poisson algebra homomorphism : =( ) a , bc^(d ) , where \ { a b , c d}=\{a , c } bd + ac \
{ b , d } , [ pois2 ] is the natural poisson structure on @xmath107 .
now , the relevance of lie bialgebras arises due to their one - to - one correspondence with poisson
lie groups , a result that can be explicitly stated as follows .    * theorem  @xcite*.    _ let @xmath18 be a lie group with lie algebra @xmath13 : _
\a ) if @xmath108 is a poisson  lie group , then @xmath13 has a natural lie bialgebra structure @xmath94 , called the tangent lie bialgebra of @xmath108 .
\b ) conversely , if @xmath18 is connected and simply connected , every lie bialgebra structure @xmath94 is the tangent lie bialgebra of a unique poisson structure on @xmath18 which makes @xmath18 into a poisson  lie group .
explicitly , if @xmath108 is a pl group , the canonical lie algebra structure on @xmath109 is given by _
d^:=(d\{f_1,f_2})_e=^ ( ^1 ^2 ) , ^1,^2 d^where @xmath110 are such that @xmath111 . in other words ,
the linearization of the pl bracket in terms of the local canonical coordinates @xmath112 on the @xmath18 group manifold is just @xmath113 , the dual of the cocommutator map that defines the associated tangent lie bialgebra structure .
a novel construction of coalgebra deformations of lie  poisson algebras can be envisaged if drinfeld theorem is revisited by taking the tangent lie bialgebra @xmath94 as the initial data for the construction of the associated pl group .
let us assume that we want to construct a poisson coalgebra deformation of a given lie
poisson algebra @xmath19 .
then , we can proceed as follows :    1 .
take the lie algebra @xmath14 , find a non  trivial lie bialgebra structure @xmath83 and construct its dual lie bialgebra @xmath94 .
now , recall that drinfeld theorem ensures that @xmath94 will be the tangent lie bialgebra structure of a certain poisson  lie group @xmath108 such that @xmath114 .
2 .   introduce the isomorphic lie algebra @xmath101   in order to have an explicit deformation parameter in the tangent lie bialgebra .
3 .   construct the pl group @xmath108 , that will provide the deformation of the poisson coalgebra @xmath64 . in order to obtain it
, we have firstly to find a faithful representation @xmath115 of @xmath101 , and use it to construct the matrix element of the ( connected component ) of the lie group @xmath18 through the usual exponentiation : @xmath116 4 .   with an abuse of notation
, we will make use the same symbol @xmath117 for the coordinate functions on @xmath5 and for the local coordinates of the @xmath18 group in ( [ matrixgroup ] ) . in this way
, the pullback of the group multiplication will provide a coproduct map _ : ( d ) ( d ) ( d ) , where , due to the coordinate identification that we have previously made , the coproduct @xmath21 is a @xmath24 deformation of the primitive one for @xmath64 , namely @xmath118 note that in the ` non - deformed case '   @xmath119 we have from   that the dual lie algebra @xmath13 is @xmath120=0 $ ] , and then the group @xmath18 is abelian .
therefore , its group law is just the addition of group coordinates , which is algebraically expressed through the primitive coproduct  .
5 .   therefore , by drinfeld theorem
, there exists a unique poisson bracket @xmath20 on @xmath22 and such that @xmath121 is a poisson  lie group @xmath122 it turns out that the poisson bracket @xmath20 is also an @xmath24 deformation of the lie
poisson bracket @xmath6 on @xmath19 , in the sense that @xmath123 indeed , the poisson
lie bracket @xmath20 has to be explicitly found by solving simultaneously the equations   and  , but this problem turns out to be computationally feasible .    due to the properties ( [ deltadeformation ] ) and ( [ bracketdeformation ] )
we can conclude that , as a poisson coalgebra , the poisson
lie group @xmath121 is an @xmath24deformation of the lie
poisson algebra @xmath26 .
indeed , there will be as many poisson - coalgebra deformations of @xmath19 as different pl groups @xmath18 such that @xmath19 induces an admissible tangent lie bialgebra cocommutator @xmath15 for @xmath114 .
in the sequel we will make use of the fact that the integrable rssler system is an outstanding example of lie  poisson coalgebra - symmetric system .
this will allow the construction of an @xmath24deformed coproduct map that will give rise to integrable coupled systems of odes of rssler type .
let us start by considering the real lie algebra a@xmath124 in the basis given in  @xcite ( note that this is exactly the four dimensional algebra  ) @xmath125=-w , \quad & [ x , z]=-z , \quad & [ y , z]=0 , \quad&[w,\cdot ] = 0 .
\label{laconm } \end{array}\ ] ] it is straightforward to check that the following cocommutator map _ ( x)= x z , _ ( y)= z w , _ ( z)=_(w)= 0 , is compatible with the lie  poisson algebra ( [ la ] ) in the sense of  , and thus defines a lie bialgebra structure @xmath126 . the dual lie algebra @xmath101 obtained from the dual cocommutator map is @xmath127=0 , \qquad [ \hat x,\hat z]= \eta\ , \hat x , \qquad [ \hat x,\hat w]=0 ,   \nonumber \\ & & [ \hat y,\hat z]= 0 , \qquad [ \hat y,\hat w]= 0 , \qquad [ \hat z,\hat w]= \eta\ , \hat y. \label{pda}\end{aligned}\ ] ] it can be checked that the dual lie algebra   is also isomorphic to the @xmath128 algebra . therefore , we have a self  dual lie bialgebra structure .
now we have to construct the lie group @xmath18 whose lie algebra is @xmath101 . if we define the @xmath129 matrix @xmath130 as the one with the only nonvanishing entry @xmath131 in the @xmath132th row and @xmath133th column , a faithful representation @xmath115 of the algebra   is given by @xmath134 by using @xmath115 , we construct the matrix group element @xmath135 where @xmath136 are the local coordinate functions on the group @xmath18 . at this point
it is worth stressing that the parameter @xmath24 appears as as multiplicative factor within all the nonvanishing matrix entries of the fundamental representation of the lie algebra @xmath13  .
therefore , when we exponentiate in order to obtain the group element @xmath18 , the powers of @xmath24 will give us the corresponding contributions coming from different powers of the lie algebra generators .    by solving the functional equations @xmath137 through the algorithm presented in  @xcite , we get the group law , _
i.e. _ , the coproduct maps for the coordinate functions of @xmath18 , which read @xmath138 that are indeed @xmath24-deformations of the primitive coproduct ( [ coproducts ] ) such that @xmath139 . again , by following the approach described in  @xcite , the unique poisson
lie bracket @xmath140 for which the coproduct @xmath21 is a poisson map is found to be @xmath141 as expected , @xmath142 and we have thus obtained a @xmath24deformation of the lie
 poisson algebra ( [ la ] ) . finally , the casimir functions for the algebra   are found to be @xmath143 and , as expected , they are again @xmath24deformations of the casimir functions  .    therefore , the poisson
lie group @xmath121 has been fully constructed .
now , by taking the same hamiltonian   of the rssler system , the pl bracket   gives rise to the system of odes @xmath144 , \cr \displaystyle \dot{y}=x w e^{-\eta z}=w x-(w x z ) \eta + \left(\frac{1}{2 } w x z^2\right ) \eta ^2+o[\eta^3 ] , \\
\displaystyle \dot{z}=x\frac{\left ( 1-e^{-\eta z }   \right ) } { \eta } = x z- \left(\frac{1}{2 } x z^2\right ) \eta + \left(\frac{1}{6 } x z^3\right ) \eta ^2 + o[\eta^3 ] , \cr \displaystyle \dot{w}=0 , \end{array }   \label{req}\ ] ] that , for @xmath145 , provides an integrable deformation of the rssler system  , with deformed integrals of the motion given by   and .
the preservation of the closed nature of the trajectories under deformation is clearly appreciated in figure 1 .
( 140,66)(0,0 ) [ figure2 ]      once the poisson
hopf algebra @xmath121 has been fully determined , the construction of an integrable deformation for the @xmath47coupled system   is straightforward .
namely , in the @xmath146 case the deformed coproduct @xmath21 on the hamiltonian of the rssler system ( [ rh ] ) will be _
( ) = ( [ _ ( x)]^2+[_(y)]^2)+_(z ) , and from   we get @xmath147 this hamiltonian gives rise to the following integrable deformation of the @xmath146 coupled rssler system @xmath148 the three integrals of the motion are given by the coproduct of the casimir @xmath149 _
( c^(2)_)= ( - ) , together with the ` one  site '  deformed casimirs ^(1)_= ( - ) , ^(2)_= ( - ) .
the generalization of this construction to an arbitrary number @xmath47 of integrable coupled rssler systems is straightformard by making use of the @xmath47-th order coproduct map , that can be explicitly obtained through the recurrence  .
the approach described above can be applied to any other hamiltonian system of odes endowed with lie  poisson symmetry . in particular , the celebrated lorenz equations @xcite are given by the system @xmath150 here we will consider a conservative limit of ( [ lorenz ] ) ( see @xcite ) , which is obtained through the following rescaling @xmath151 and by taking the limit @xmath152 . as a consequence of this transformation ,
the system ( [ lorenz ] ) is mapped into @xmath153 which is a hamiltonian system with respect to the lie
poisson algebra @xcite @xmath154 and provided that we have considered the hamiltonian function @xmath155    in order to proceed in a similar way as for the rssler system , let us consider the 4d lie algebra @xmath156=-\dfrac{1}{2 } ( z+w ) , \qquad [ x , z]=\dfrac{1}{2 } y , \qquad [ y , z]=x , \qquad [ w,\cdot]=0 .
\label{liel}\ ] ] this lie algebra is isomorphic to the @xmath157 lie algebra , where @xmath158 is a central generator . obviously , the lie  poisson algebra @xmath19 arising from   is a generalization of the poisson algebra   through a new central coordinate function @xmath60 @xmath159 now , if we consider the hamiltonian @xmath160 the equations of motion obtained through   are @xmath161 and we recover the conservative lorenz equations when @xmath145 . finally , the conserved quantities for the system are given by the casimir functions for ( [ clt ] ) , namely @xmath162    with this prerrequisites we can proceed with the formalism presented in section 4 . a cocommutator map compatible with the lie algebra ( [ liel ] )
is found to be @xmath163 the dual lie algebra @xmath101 obtained from the dual cocommutator map @xmath164 is given by @xmath127=0 , \qquad [ \hat x,\hat z]= \eta\ , \hat x , \qquad [ \hat x,\hat w]= \eta\ , \hat x + \eta\ , \hat y ,   \nonumber \\ & & [ \hat y,\hat z]= \eta\ , \hat y , \qquad [ \hat y,\hat w]= \frac{\eta}{2}\ , \hat x + \eta\ , \hat y , \qquad [ \hat z,\hat w]= 0 .
\label{pd}\end{aligned}\ ] ] this lie algebra @xmath101 can be shown to be isomorphic to the direct sum of two copies of the two  dimensional real lie algebra @xmath165=e_1 $ ] .
a four dimensional faithful representation @xmath115 of @xmath101 is given by the matrices @xmath166 by using @xmath115 , we construct the matrix group element @xmath167 as a consequence , the coproduct map arising from the group multiplication of two @xmath18 matrices is found to be @xmath168,\\ \delta_\eta(y)=y_{1}+e^{-\eta\ , ( z_{1}+w_{1})}\left[y_{2}\cosh \left(\tfrac{\eta } { \sqrt{2}}\ , w_{1}\right)-\sqrt{2}x_{2}\sinh \left(\tfrac{\eta}{\sqrt{2}}\ , w_{1}\right)\right ] , \\
\delta_\eta(z)=z_1+z_2,\\ \delta_\eta(w)=w_1+w_2 . \end{array } \label{colo}\ ] ]    now , by applying again the algorithmic procedure introduced in  @xcite , we get the following unique solution for the poisson  lie structure on the group @xmath18 that has @xmath19 as its linearization and is compatible with  : @xmath169 this is indeed an @xmath24deformation of the @xmath170 lie  poisson algebra  , which is recovered in the @xmath171 limit , and the ( deformed ) casimir functions for the poisson algebra ( [ dt ] ) are found to be @xmath172 - 2 e^{z \eta } [ 1+w \eta   ( 1+w \eta ) ] \right)}{\eta ^2}. \label{casl}\ ] ] note that , as expected , @xmath173    therefore , an integrable deformation of the integrable lorenz system   is obtained by considering the same hamiltonian   and the poisson  lie brackets  , and reads @xmath174 , \\[2pt ] \dot{z}= & x \ , y , \end{array } \label{cl2def}\ ] ] where the only deformed equation occurs for the @xmath74 variable , and the conserved quantities are given by the hamiltonian and the casimir functions  .
several periodic trajectories for this system are plotted in figure 2 , and are compared with the underformed lorenz ones .
( 140,66)(0,0 ) [ figure2 ]      the undeformed lie
poisson coalgebra @xmath64 is again defined through the primitive coproduct map  .
therefore , the @xmath146 coupled lorenz hamiltonian will be given by @xmath175^{2},\ ] ] where the coproduct @xmath41 is the undeformed one  .
the corresponding equations of motion read @xmath176 \dot{x}_{2}= y_2,\\ \dot{y}_{2}= 2x_2-(x_1+x_2)(w_2+z_2),\\ \dot{z}_{2}= y_2(x_1+x_2 ) , \end{array } \label{couplo}\ ] ] together with @xmath68 ( see  @xcite for some instances of nonintegrable coupled lorenz oscillators ) .
again , three additional independent integrals of the motion in involution for this system are provided by the coalgebra symmetry : namely , the casimir functions on each of the two copies of the lie  poisson algebra _
1= -2x_1 ^ 2+y_1 ^ 2+z_1(z_1 + 2 w_1 ) , _
2= -2x_2 ^ 2+y_2 ^ 2+z_2(z_2 + 2 w_2 ) , together with the coproduct of the casimir @xmath71 , which reads @xmath177    the @xmath47dimensional integrable generalization of this system arises from the @xmath47-th coproduct map @xmath73 of the hamiltonian , namely @xmath178 in this way the following @xmath47-coupled lorenz system is obtained as @xmath179 and the ( @xmath79 ) integrals of the motion ( in involution between themselves and with respect to the hamiltonian @xmath81 ) are given by  @xcite the ` one  site ' integrals _
i=-2x_i^2+y_i^2+z_i(z_i+2 w_i ) , i=1,
,n together with the @xmath48th order coproducts of the casimir functions @xmath180 , \quad k=2,\dots , n.\ ] ] note that the coupling between the 3d dynamical systems is induced by the non  linear @xmath181 term in the hamiltonian .      by following the very same approach as for the rssler system
, the deformed coproduct map @xmath21   acting on the hamiltonian function ( [ h1 ] ) gives rise to _ ( ) = 2(_(z)-_(w))- [ _ ( x)]^2 , and by taking into account the explicit expressions for @xmath21 we get ( ) = 2(z_1+z_2-w_1-w_2)-^2 .
this hamiltonian gives rise to the following ( quite involved ) deformation of the @xmath146 coupled lorenz system @xmath182\notag\\ % % % % % % % % % % % % % % % % % % % % & & + y_{1}\dfrac{\eta}{4}e^{-2\eta ( w_{1}+z_{1})}\left [ ( 2x_{2}^{2}+y_{2}^{2})\cosh(\sqrt{2}\,\eta\ , w_{1})-2\sqrt{2}\,x_{2}\,y_{2}\sinh(\sqrt{2}\,\eta\ , w_{1 } ) \right],\notag\\ % % % % % % % % % % % % % % % % % % % % % % % % \notag\\ % % % % % % % % \dot{y}_{1}&= & \dfrac{x_{1}}{2\eta}\left ( -1+\eta ( 4+\eta(2x_{1}^{2}-y_{1}^{2 } ) ) \right)\notag\\ % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % & & + \dfrac{x_{1}}{2\eta}e^{-2\eta ( w_{1}+z_{1 } ) } \left ( 1+\eta^{2}(2x_{2}^{2}+y_{2}^{2})\cosh(\sqrt{2}w_{1}\eta ) + \eta^{2}(2x_{2}^{2}-y_{2}^{2}-2\sqrt{2}x_{2}y_{2}\sinh(\sqrt{2}w_{1}\eta ) ) \right)\notag\\ % % % % % % % & & + \dfrac{1}{4\eta}e^{-3\eta ( w_{1}+z_{1})}\left [ \left ( 2x_{2}\cosh\left(\tfrac{\eta}{\sqrt{2}}w_1\right)-\sqrt{2}y_{2}\sinh\left(\tfrac{\eta}{\sqrt{2}}w_1\right ) \right)\left(1+e^{2\eta ( w_{1}+z_{1})}(-1+\eta^{2}(6x_{1}^{2}-y_{1}^{2 } ) ) \right ) \right],\notag\\ % % % % % % % % % % % % % % % % % % % % % \notag\\ % % % % % % % % \dot{z}_{1}&= &   x_{1}y_{1}+\dfrac{1}{2}y_{1}e^{-\eta(w_{1}+z_{1})}\left ( 2x_{2}\cosh \left(\tfrac{\eta}{\sqrt{2}}w_1\right)-\sqrt{2}y_{2}\sinh\left(\tfrac{\eta}{\sqrt{2}}w_1\right ) \right ) , \notag\\ % % % % % % % \notag\\ % % % % % % % % \dot{x}_{2}&= & y_{2}+\dfrac{1}{8\eta}e^{-2\eta(w_{1}+w_{2}+z_{1}+z_{2})}\left ( 1+e^{2\eta ( w_{2}+z_{2})}(-1+\eta^{2}(2x_{2}^{2}-y_{2}^{2 } ) ) \right)\times \label{deflor2}\\ % % % % % % % & & \times \left ( y_{2}-y_{2}\cosh(\sqrt{2}w_{1}\eta)+2\sqrt{2}x_{1}e^{\eta(w_{1}+z_{1})}\sinh\left(\tfrac{\eta}{\sqrt{2}}w_1\right)+\sqrt{2}x_{2}\sinh(\sqrt{2}w_{1}\eta ) \right),\notag\\ % % % % % % % \notag\\ % % % % % % % % % \dot{y}_{2}&= & 2x_{2}+\dfrac{1}{4\eta}e^{-2\eta ( w_{1}+w_{2}+z_{1}+z_{2})}\cosh\left(\tfrac{\eta}{\sqrt{2}}w_1\right )     \left(1+e^{2\eta(w_{2}+z_{2})}(-1+\eta^{2}(2x_{2}^{2}-y_{2}^{2}))\right)\times\notag\\ % % % % % % % % % % & & \times \left ( 2x_{1}e^{\eta ( w_{1}+z_{1})}+2x_{2}\cosh\left(\tfrac{\eta}{\sqrt{2}}w_1\right)-\sqrt{2}y_{2}\sinh\left(\tfrac{\eta}{\sqrt{2}}w_1\right ) \right),\notag\\ % % % % % % \notag\\ % % % % % % % \dot{z}_{2}&= & e^{-2\eta(w_{1}+z_{1})}\left ( x_{1}y_{2}e^{\eta(w_{1}+z_{1})}\cosh\left(\tfrac{\eta}{\sqrt{2}}w_1\right)+x_{2}y_{2}\cosh(\sqrt{2}\eta w_{1 } ) \right)\notag\\ % % % % % % % % % % % % % % & & -\dfrac{\sqrt{2}}{4}e^{-2\eta(w_{1}+z_{1})}\left ( 4x_{1}x_{2}e^{\eta(w_{1}+z_{1})}\sinh\left(\tfrac{\eta}{\sqrt{2}}w_1\right)+(2x_{2}^{2}+y_{2}^{2})\sinh\left(\sqrt{2}w_{1}\eta\right ) \right ) .
\notag\end{aligned}\ ] ] by taking a power series expansion we get a first order deformation in @xmath24 of the coupled lorenz system   given by @xmath183,\\[2pt ] \dot{y}_{1}=2x_{1}-(x_{1}+x_{1})(w_{1}+z_{1 } ) \\[2pt ] % % % % % % % % % % \quad\quad+ \left ( ( w_{1}^{2}+z_{1}^{2}+2 w_{1}z_{1})(x_{1}+2x_{2})-\dfrac{y_{1}^{2}}{2}(x_{1}+x_{2})+\dfrac{w_{1}y_{2}}{2}(w_{1}+z_{1})+x_{1}(x_{1}^{2}+3 x_{1}x_{2}+2x_{2}^{2 } ) \right)\eta + o[\eta^2],\\ \dot{z}_{1}=y_1(x_1+x_2)-\dfrac{y_{1}}{2}(w_{1}(2x_{2}+y_{2})+2z_{1}x_{2})\eta + o[\eta^2 ] \end{array}\ ] ] @xmath184,\\[2pt ] % % % % % % % % % % % % % % \dot{y}_{2}=(2 { x_2}-({x_1}+{x_2 } ) ( { w_2}+{z_2 } ) ) \\ % % % % % % % % % \quad\quad+\dfrac{1}{2 } \left ( ( x_{1}+x_{2})(2x_2 ^ 2-y_2 ^ 2 + 2(w_2+z_2)^2)+(w_2+z_2)(2 ( x_1 + 2x_2)(w_1+z_1)+w_1 y_2)\right)\eta + o[\eta^2],\\[2pt ]   % % % % % % % % % % % % % % % % \dot{z}_{2}=y_2({x_1}+{x_2})-\dfrac{1}{2 } \left ( 2 y_2 ( w_1+z_1)(x_1 + 2x_2)+w_1 y_2 ^ 2 + 2 w_1 x_2 ( x_1+x_2 )
\right)\eta   + o[\eta^2 ] .
\end{array}\ ] ]    finally , the three integrals of the motion for   are given by the coproduct @xmath21   of the casimir @xmath149 given in  , namely    l _
( c^(2)_)= e^-_(z ) ( 1+e^2 ( _ ( w)+_(z ) ) - + - 2 e^ _ ( z ) [ 1+_(w ) ( 1+_(w ) ) ] ) ,    together with the ` one  site '  deformed casimirs    l ^(1)_= , + ^(2)_= .
of course , since _
0\{_i,_()}_=\{_i , ( ) } , _
i=\{x_i , y_i , z_i } , all these expressions define an integrable @xmath24-deformation of the @xmath146 coupled lorenz system , and the generalization for arbitrary @xmath47 can be straightforwardly obtained through the @xmath47-th coproduct map that arises from   and from the recurrence  .
in this paper we have established a novel ( and somehow unexpected ) link between the theory of pl groups and the integrability problems for ( coupled ) odes .
in fact , since poisson  lie groups can be thought of as deformations of lie  poisson coalgebras , this property can be used in order to construct integrable deformations of ( coupled ) systems of odes endowed with a lie  poisson hamiltonian structure .
the applicability of this statement is based on the fact that the explicit construction of the pl groups that deform a given lie
 poisson algebra can be systematically obtained through the exponentiation of the associated lie bialgebra structures ( therefore , by making use of the poisson version of the well  known ` quantum duality principle '  for quantum groups  @xcite ) .
one of the main advantages of the approach here presented relies on its algorithmic nature , that makes it amenable to be computerized in a nearly complete way ( see  @xcite ) . in order to illustrate the method , in this paper
we have worked out two specific integrable deformations of ( coupled and uncoupled ) rssler and lorenz systems .
it is important to emphasize that other integrable deformations for these systems do exist .
in fact , there will be as many of them as the number of possible lie bialgebra structures associated to the initial lie
poisson algebra that provides the hamiltonian structure for the undeformed dynamical system ( namely , the a@xmath124 algebra   for the rsler system and the @xmath185 algebra   for the lorenz one ) .    indeed , the very same approach can be applied to many other hamiltonian systems of odes .
moreover , there also exists the possibility of generalizing this method in order to consider integrable systems with more than one hamiltonian , either by deforming bihamiltonian poisson pencils or by using the loop coproduct approach introduced in  @xcite .
finally it is worth mentioning that the analytic solution of coupled integrable systems with coalgebra symmetry is facilitated by the use of the so called ` cluster '   variables ( see  @xcite ) , and this fact could shed some light on the dynamical significance of the rssler and lorenz couplings here obtained .
work on all these lines is in progress and will be presented elsewhere .
this work was partially supported by the spanish ministerio de economa y competitividad ( mineco ) under grant mtm2013 - 43820-p , by junta de castilla y len under grant bu278u14 and by infn
sezione roma tre .