Abstract
We study the integrability of a model of elastic satellite whose centre of mass moves in a circular Keplerian orbit around a gravity centre. The satellite is modelled by two point masses connected by an extensible massless spring that obeys Hooke’s law. It is assumed that the distance between point masses is much smaller than the radius of the orbit, so the orbital motion of the satellite is not perturbed by its rotational motion. The gravity potential of the satellite is expanded into a series with respect to its size up to quadratic terms which describe the gravity gradient torque acting on the satellite. Two cases are considered with Hooke’s centre localised in the centre of mass of the dumbbell and at an arbitrary point along a line connecting both masses. It is shown that the first case appears to be integrable and super-integrable for selected values of the parameter of the system. In the second case, model depends effectively only on one parameter and is non-integrable. In the proof, differential Galois integrability obstructions are used. For the considered sysem, these obstructions are deduced thanks to the recently developed symplectic Kovacic’s algorithm in dimension 4. According to our knowledge, this is the first application of this tool to a physical model.
Similar content being viewed by others
1 Introduction
We consider a model of an elastic dumbbell satellite following Sidorenko and Celletti [10]. It consists of two point masses connected by a massless spring and moves in a central gravity field. We assume that the dumbbell is short; thus, its orbital motion decouples from its rotational motion. The dumbbell satellite models have attracted the attention of scientists because they are suitable for an investigation of the general properties of a rigid or elastic deformable body motion in a gravity field.
The aim of this article is the integrability analysis of this model. Two sets of coordinates: Cartesian and spherical are used to show different aspects of the dynamics. The complete integrability result for physically important ranges of parameters variability is obtained using the differential Galois theory and the direct method of searching for first integrals.
The plan of the paper is as follows. In Sect. 2, we derive equations of motion for the considered model. We use two sets of coordinates which parametrise the configuration space of the problem. Using the first one, we can easily notice that the first model with Hooke’s centre located at the mass centre of the dumbbell is integrable. In Sect. 3, we show that for particular values of parameters this model is super-integrable. The second set of coordinates is used in the proof of our main theorem which states that except explicitly distinguished integrable cases the system is not integrable. So, in this paper we give necessary and sufficient conditions for integrability of the considered system. In Sect. 4, we distinguish equilibria and certain invariant sets of the system. We demonstrate also its chaotic behaviour using numerical methods. Section 5 contains a proof of our main theorem. It is based on the symplectic Kovacic algorithm recently formulated in Combot and Sanabria [1]. This algorithm can be used to investigate systems of variational equations in dimension 4. As this tool is relatively new, we collect basic theoretical facts in “Appendix”.
2 Equations of motion
Following Sidorenko and Celletti [10], we assume that two point masses \(m_1\) and \(m_2\) are connected by an extensible massless spring that obeys Hooke’s law. The centre of mass of this elastic dumbbell moves in an orbit around a gravity centre located at the origin of the inertial frame. We assume that the dimensions of the dumbbell are much smaller than the dimension of the orbit.
The radius vectors of masses \(m_1\) and \(m_2\) are
where \(\varvec{r}\) is the radius vector of the centre of mass of the dumbbell, and
The kinetic energy of the dumbbell is
where
The gravitational potential energy of point masses can be expanded into the Taylor series
where \(\kappa =GM\) is the gravitational parameter of the centre, see, e.g. Sidorenko and Celletti [10].
Thus, the gravitational potential of the dumbbell reads
where \(\varvec{e}_r=\varvec{r}/|\varvec{r} |\) is the unit vector in the direction to the centre of mass.
We assume that the elastic dumbbell is permanently straight, and the potential energy of its elastic deformations is given by
where \(d_0\) is the length of the undeformed spring. If \(d_0=0\), then Hooke’s centre is localised at the centre of mass of the dumbbell. We will discuss this case later.
In a typical study of the attitude motion of a satellite, it is usually assumed that its orbital motion is not perturbed by its rotation. Up to rotation of the Cartesian coordinates, we can assume that the mass centre of the dumbbell moves in the (x, y)-plane of the inertial frame along a circular Keplerian orbit with radius \(a:=|\varvec{r} |\). Thus, the Lagrange function of the system is
It depends explicitly on time because components of \(\varvec{e}_r\) are time dependent. To remove this dependence, we pass to the orbital frame with the origin at the centre of mass of the dumbbell and with axes \(\{\varvec{s},\varvec{t},\varvec{n}\}\), where \(\varvec{n}\) and \(\varvec{t}\) are unit vectors normal and tangential to the orbit, and \(\varvec{s}=\varvec{t}\times \varvec{n}=\varvec{r}/r=\varvec{e}_r\), see Fig. 1. We set
where \(q_i\) are component of vector \(\varvec{d}\) with respect to the orbital frame,
and \(\omega ^2=\kappa /a^3 \) is the orbital angular velocity. Then,
According to our assumptions, \(\varvec{s}=\varvec{e}_r\) and \(\varvec{\omega }= \omega [0,0,1]^T\). In terms of \((\varvec{q},\varvec{\dot{q}})\) variables, the Lagrangian reads
Assuming that \(d_0\ne 0\), we can rescale variables \(q_i\rightarrow d_0 q_i\), and time \(t\rightarrow \omega t\). Then, the above Lagrangian transforms into \(L=\mu \omega ^2d_0^2{\widetilde{L}}\), where
and \(\gamma :=c/(\mu \omega ^2)\) is a new parameter. Let us note that \(\gamma \ge 0\) that will be used in further considerations. The standard Legendre transformation of the above Lagrangian gives the following Hamiltonian function
The corresponding Hamilton equations read
Now, we introduce new variables \((\xi , \varphi ,\theta )\) useful in further analysis
These variables are well defined for \(\theta \ne 0,\pi .\) The angles \(\varphi \) and \(\theta \) are defined as it is shown in Fig. 2. In these variables, Lagrangian \({\widetilde{L}} \) given in (2.13) takes the form
and making its Legendre transformation, we obtain the following Hamilton function
The corresponding Hamilton equations take the form
At this point, we explain why we derived equations of motion in two sets of coordinates in the configuration space of the system. In the Cartesian coordinates, the state of the system is described by \((\varvec{q},\varvec{\dot{q}})\) or \((\varvec{q},\varvec{p})\) which are ‘almost’ global, that is \((\varvec{q},\varvec{p})\in (\mathbb {R}^3\setminus \{\varvec{0}\})\times \mathbb {R}^3\). We have to remove \(\varvec{q}=\varvec{0}\) from the configuration space because if \(c d_0\ne 0\), then the elastic potential
is not differentiable at \(\varvec{q}=\varvec{0}\). Moreover, if \(\gamma =0\), then the Hamiltonian (2.14) is a homogeneous quadratic polynomial function of \((\varvec{q},\varvec{p})\), so equations of motion (2.15) are linear and thus, they are integrable. However, this property is not obvious in other coordinates, see Hamiltonian (2.18) and Eq. (2.19).
On the other hand, we will use coordinates \((\xi ,\varphi ,\theta )\) in our proof of the main theorem of this paper, namely
Theorem 2.1
If \(\gamma > 0\), then the system given by Hamiltonian (2.18) is not integrable in the Liouville sense with first integrals which are meromorphic in \(p_{\varphi },p_{\theta },p_{\xi },\xi ,\cos \theta ,\sin \theta ,\cos \varphi ,\sin \varphi \).
Remark that the variable change between Cartesian coordinates and coordinates (2.16) transforms meromorphic functions in \({\varvec{p}, \varvec{q},}\sqrt{q_1^2+q_2^2+q_3^2}\) to meromorphic functions in \(p_{\varphi },p_{\theta },\) \(p_{\xi },\xi ,\cos \theta ,\sin \theta ,\cos \varphi ,\sin \varphi \). Thus, Theorem 2.1 forbids the existence of additional first integrals of Hamiltonian (2.14) meromorphic in \({\varvec{p}, \varvec{q},}\sqrt{q_1^2+q_2^2+q_3^2}\).
The proof of this theorem is quite long and the proper choice of coordinates allows to avoid several analytical difficulties and complications.
In the end, let us underline that we derive the Hamiltonians mentioned in the above theorem under assumption that \(d_0\ne 0\). The case \(d_0=0\) is itself interesting, and it is considered in the next section.
3 Integrable and super-integrable cases
In order to have possibility to investigate cases with \(d_0=0\), we rescale only time variable \(t\rightarrow \omega t\) in the Lagrange function (2.12). Then, \(L=\mu \omega ^2 {\widetilde{L}}\), where now
The corresponding Hamilton function is
Let us notice that for \(d_0=0\) or \(\gamma =0\) this Hamiltonian is a homogeneous polynomial of degree two with respect to phase variables \((\varvec{q},\varvec{p})\) and thus, its equations of motion are linear in variables.
In the remaining part of this section, we will consider case \(d_0=0\) with simplified Hamiltonian
This system is integrable with commuting first integrals
The examples of time evolution of our system for two different initial conditions are presented in Figs. 3 and 4. Figures 3a and 4a show spatial motion of the dumbbell vector \(\varvec{q}\) in the orbital system, and Figs. 3b and 4b the corresponding projections of these trajectories on plane \((q_1,q_2)\). The inclination \(\theta \) oscillates between \(-\pi \) and 0, see Figs. 3c and 4c. The time changes of the azimuth angle \(\varphi \) are more complicated as the vector \(\varvec{q}\) rotate and the direction of these rotations is changing more or less periodically, see Figs. 3d and 4d. The deformation parameter of the dumbbell \(\xi \) oscillates with high frequency which illustrate Figs. 3e and 4e .
If we additionally assume that \(\gamma =0\), then one can find one more first integral. Thus, in this case the system has three additional first integrals
They are functionally independent together with H and satisfy the following commuting relations
Thus, if \(\gamma =0\), then the system is super-integrable, but \(H,I_1,I_2,I_3,I_4\) are algebraically dependent as the following relation holds true
It is natural to ask if for other values of \(\gamma \) the system is super-integrable. To answer this question, we observe that for a generic value of \(\gamma \) the system is integrable so its invariant tori are three-dimensional manifolds in the phase space. With each such a torus, we have related three periods which are independent over \(\mathbb {Z}\). If one more additional first integral appears, then its common level with the invariant torus is typically a torus of dimension two. So, if the system is super-integrable, then the three-dimensional torus is foliated by two-dimensional tori. This is why the three periods cannot be \(\mathbb {Z}\) independent. Thus, if the system is super-integrable, then a resonance between basic frequencies appears.
For the considered system, this approach is particularly simple because Hamilton equations are linear
with \(\varvec{x}=[\varvec{q},\varvec{p}]^T\) and with matrix
To identify characteristic frequencies of this system, let us check the characteristic polynomial of this matrix
After substitution \(\lambda =\mathrm {i}\omega \), characteristic equation takes the form
Its solutions are of the form \(\pm \omega _k\), \(k=1,2,3\), where
Thus, if the matrix is diagonalisable, using a canonical change of variables the Hamiltonian (3.3) can be transformed to the form
where now \((I_i,\varphi _i)\) are canonical action-angle variables and \(\sigma _{i}=\varepsilon _i\omega _i\), with \(\varepsilon _i\in \{-1,+1\}\), for \(i=1,2,3\). Clearly, \(I_1\), \(I_2\) and \(I_3\) are first integrals of the system. If there is a resonance between frequencies of the form
then
is a first integral of the system. The number \(|n |:=|n_1|+|n_2|+|n_3|\) is called the order of the resonance. Instead of \(\psi \) it is more useful to consider the first integral \(I_4=\sin \psi \) (or \(I_4=\cos \psi \)). Since \(\sin \psi \) and \(\cos \psi \) are polynomials of degree \(|m |\) in canonical variables, the additional first integral for cases corresponding to resonance conditions (3.8) can be chosen a homogeneous and polynomial function in variables \((\varvec{q},\varvec{p})\) of degree equal to the order of resonance |m|.
Below, we list triples of integers \((m_1,m_2,m_3)\in \mathbb {Z}^3\) with \(|m|\le 4\) that corresponds to \(\gamma \ge 0.\)
Let us notice that for \(\gamma =0\) we have \(\omega _1-\omega _3=0\) and \(\omega _2=0\), thus \(|m|=2\), and this is the only resonance of order 2 with \(\gamma \ge 0\). In this case, matrix \({\mathbb {A}}\) is not diagonalisable.
There are three resonances of degree 3. For the resonance, \(\omega _1+\omega _2-\omega _3=0\) corresponding to \(\gamma =4\) one can find the following first integral
Resonances \(\omega _1-2 \omega _2=0\) and \(\omega _3-2\omega _2=0\) occur for \(\gamma =\tfrac{29}{9}+\tfrac{2 \sqrt{217}}{9}\) and \(\gamma =\frac{91}{18}+\frac{5 \sqrt{337}}{18}\), respectively. We do not list the explicit forms of the corresponding first integrals because they are too long.
There are also three resonances of order four. If \(\gamma =\frac{51}{32}+\frac{7 \sqrt{57}}{32}\), then resonance \(\omega _1+2\omega _2-\omega _3=0\) occurs and the additional first integral is rather complicated
Resonances \(\omega _1-3\omega _2=0\) and \(3\omega _2-\omega _3=0\) occur for \(\gamma =\tfrac{67}{32}+\tfrac{9 \sqrt{57}}{32}\) and \(\gamma =\tfrac{21}{8}+\tfrac{15}{4 \sqrt{2}}\), respectively. Here, we do not list the form of additional first integrals for these cases.
In a generic Liouville integrable case, connected compact common levels of three independent first integrals are three dimension tori. In a super-integrable case the dimension of invariant tori is smaller. Let us consider a common level of four independent first integrals
where h, and \(\alpha _i\) are real constants. It is a two-dimensional surface in \(\mathbb {R}^6\). We can visualise it in the configuration space \(\mathbb {R}^3\) and eliminate momenta from polynomial equations defining level \({{\mathcal {M}}}\). As result, we obtain a polynomial \(P(\varvec{q},h,\alpha _1,\alpha _2,\alpha _3)\). Its zero level defines an algebraic surface which is the image of the invariant torus in the configuration space. Polynomial P is a product of irreducible polynomials \(P=P_1\cdots P_k\), \(k>1\) and \(P_i=0\) defines a connected component of the surface. In Fig. 5, we present components of these algebraic surfaces selected by the choice of the initial condition \(q_1(0)=\tfrac{1}{2}\), \(q_2(0)=\tfrac{1}{10}\), \(q_3(0)=-\frac{1}{5}\), \(p_1(0)=p_2(0)=p_3(0)=0\). The trajectory for this initial condition is shown as black line on these surfaces. Surfaces shown in Fig. 5c, f are cylinders as the polynomials defining them do not depend of \(q_3\). Thus, for these cases invariant manifolds \({{\mathcal {M}}}\) are not compact.
Time evolution of inclination \(\theta \), azimuth angle \(\varphi \) and the deformation parameter of the dumbbell \(\xi \) corresponding to these super-integrable cases is shown in Fig. 6.
If relation (3.8) holds for two non-colinear \(m_1,m_2,\) \( m_3\), the system is doubly resonant, and then it would admit two additional first integrals. As \(\omega _1=\sqrt{\gamma +1}>0\), we can divide by it relation (3.8), and noting
we obtain that both \(\rho _1,\rho _2\) should be rational.
Proposition 3.1
The system (3.5) has a double resonance only for \(\gamma =0\).
Let us recall that for \(\gamma =0\) matrix \({\mathbb {A}}\) is not diagonalisable and this is why the system is super-integrable in this case with just one additional first integral.
Proof
Taking both Eqs. (3.10) and eliminating \(\gamma \) by a resultant formula, we obtain the relation
We are looking for rational solutions of this curve. Noting \(s=\rho _1^2\), if \(\rho _1,\rho _2\) is a rational solution of (3.11), then \(s,\rho _2\) is a rational solution of
This curve is an elliptic curve, and its Weierstrass form is \(v^2=u^3+191u+99198\). Magma computational algebra system manages to compute its Mordell–Weil group, and it admits only 6 rational points, leading to the following rational solutions \((s,\rho _2)\)
For the last one, s is not a rational squared, thus does not lead to a rational solution of (3.11). Solution (0, 1) gives \(\gamma =0\), and the others do not lead to a nonnegative \(\gamma \). \(\square \)
Removing the radicals in (3.8), we find that \(\gamma \) should satisfy the following quartic polynomial
Then, substituting \(m_1,m_2,m_3\) by integers gives resonant \(\gamma \)’s as solutions of this equation. Relation (3.8) defines a straight line with rational slope in \(\rho _1,\rho _2\) coordinates for each triplet of \((m_1,m_2,m_3)\), and the ones leading to nonnegative \(\gamma \)’s are represented in Fig. 7 for \(|m|\le 4\).
Compared with a generic Hamiltonian system with three degrees of freedom not all resonances for our system are possible because the quartic equation could have no nonnegative roots (compare with Karabanov and Morozov [4]).
4 Basic properties of the system
We start this section with numerical examples showing behaviour in time of the dumbbell with \(d_0\ne 0\) obtained by integration of equations of motion in angular variables (2.19). Presented simulations were made for \(\gamma =500\). All simulations were performed using software Mathematica with working precision at least 13 so that precision of 13 digits has been maintained during internal computations. We chose four initial conditions given below for which the dynamics looks different
-
IC 1:
\(\varphi (0)=0\), \(\theta (0)=\tfrac{\pi }{2}-\tfrac{1}{20}\), \(\xi (0)=-\tfrac{1}{\gamma +1}=-\frac{1}{501}\), \(p_{\varphi }(0)=p_{\xi }(0)=p_{\theta }(0)=0\) that belong to the energy level \(E=0.49527\),
-
IC 2:
\(\varphi (0)=1\), \(\theta (0)=\tfrac{1}{10}\), \(\xi (0)=-\tfrac{1}{10}\), \(p_{\varphi }(0)=-\tfrac{1}{2}\), \(p_{\theta }(0)=p_{\xi }(0)=0\) that belong to the energy level \(E=19.8783\),
-
IC 3:
\(\varphi (0)=0\), \(\theta (0)=-\tfrac{1}{2}\), \(\xi (0)=\tfrac{1}{10}\), \(p_{\varphi }(0)=\tfrac{1}{2}\), \(p_{\theta }(0)=-\tfrac{2}{3}\), \(p_{\xi }(0)=0\) that belong to the energy level \(E=3.25553\),
-
IC 4:
\(\varphi (0)=\tfrac{\pi }{4}\), \(\theta (0)=-\tfrac{1}{2}\), \(\xi (0)=\tfrac{1}{10}\), \(p_{\varphi }(0)=\tfrac{1}{2}\), \(p_{\theta }(0)=-\tfrac{2}{3}\), \(p_{\xi }(0)=0\) that belong to the energy level \(E=3.45886\).
Spatial motion of the dumbbell vector \(\varvec{q}(t)\) in the orbital frame expressed in angular variables in (2.16) for the above initial conditions is presented in Fig. 8. The corresponding time evolutions of the inclination angle \(\theta \), the azimuth angle \(\varphi \) and the deformation parameter \(\xi \) are shown in Fig. 9. One can notice that for initial conditions IC 1 oscillations of \(\xi \) have small amplitude, so the end of vector \(\varvec{q}\) moves approximately on a sphere; its inclination oscillates and rotates around the \(q_3\)-axis. For remaining considered initial conditions, the length of the dumbbell \(\xi \) changes considerably and the motion of the dumbbell vector \(\varvec{q}\) is a superposition of oscillations of \(\theta \) and \(\xi \) and rotations of \(\varphi \).
These numerical experiments just illustrate the complexity and the variety of dynamics of the system. However, for further analytical and numerical investigations we need to identify its simplest invariant sets.
Hamilton’s equations (2.15) have the following equilibria: saddle-centre-centres \(L_{1,2}\), saddle-saddle-centres \(S_{1,2}\) and for elastic dumbbell satellite rigid enough with \(\gamma >3\) centre-centre-centres \(O_{1,2}\) with the \((\varvec{q},\varvec{p})\) coordinates
Eigenvalues of linearization of vector field (2.15) at the respective equilibria are following
-
for \(L_{1,2}:\)
$$\begin{aligned} \lambda _{1,2}=\pm \frac{\mathrm {i}}{\sqrt{2}}\sqrt{\Delta },\quad \lambda _{3,4}=\pm \sqrt{6}\sqrt{\frac{\gamma }{\Delta }},\quad \lambda _{5,6}=\pm \mathrm {i}, \end{aligned}$$where \(\Delta =\gamma +1 +\sqrt{\gamma (\gamma +14)+1}\);
-
for \(S_{1,2}\)
$$\begin{aligned} \begin{aligned}&\lambda _{1,2}{=}\pm \sqrt{\frac{1}{2} \left( 1{+}\mathrm {i}\sqrt{15}\right) } \!\approx \pm 1.11803 \pm 0.866025 \mathrm {i},\\&\lambda _{3{,}4}{=}\pm \sqrt{\frac{1}{2} \left( 1{-}\mathrm {i}\sqrt{15}\right) } \!\approx \pm 1.11803 \mp 0.866025 \mathrm {i},\\&\lambda _{5,6}=\pm \mathrm {i}\sqrt{\gamma +1}, \end{aligned} \end{aligned}$$ -
and for \(O_{1,2}\)
$$\begin{aligned} \lambda _{1,2}&=\pm \frac{\mathrm {i}}{\sqrt{2}} \sqrt{\Delta },\quad \lambda _{3,4}=\pm \mathrm {i}\sqrt{\frac{6(\gamma -3)}{\Delta }},\nonumber \\&\!\!\!\!\!\!\!\!\!\!\lambda _{5,6}=\pm 2\mathrm {i}, \end{aligned}$$where \(\Delta =\gamma +4 +\sqrt{ \gamma (\gamma -4) +52}\).
In variables \((\varphi ,\theta ,\xi ,p_{\varphi },p_{\theta },p_{\xi })\), equilibria are given by
Equilibria \(O_{1,2}\) lie on axis \(\varvec{s}=\varvec{0}\), where angular coordinates (2.16) are not defined.
The methods used to prove that the system is not integrable require that we know a non-equilibrium solution of the system. There is no general method how to find such a solution. However, the considered system admits an invariant two-dimensional manifold
on which it reduces to a one degree of freedom Hamiltonian system describing oscillations of vertically oriented satellite.
In variables \((\varphi ,\theta ,\xi ,p_{\varphi },p_{\theta },p_{\xi })\), this manifold is given by
The Hamiltonian system (2.15) admits also the invariant manifold
on which it is a Hamiltonian system with two degrees of freedom with Hamiltonian
where \(\varvec{q}=(q_1,q_2)\).
In angular coordinates (2.16), this manifold is defined by
On this manifold, variable \(\theta \) becomes the polar angle in the plane, Hamiltonian of the restricted system has the form
and the corresponding Hamilton equations read
This two degrees of freedom system on \({{\mathscr {M}}}_{4}\) does not seem to be integrable. This is suggested in the sequence of the Poincaré cross-sections shown in Fig. 10. Following Sidorenko and Celletti [10], we take \(\gamma =500\) that corresponds to considered in this paper value of parameter \(\beta =\tfrac{1}{\gamma }=0.002\). Cross-sectional plane was chosen \(\xi =0\), and points are generated when \(p_{\xi }>0\). Angle \(\theta \) is taken modulo \(\pi \). These cross-sections show that for low energies the dumbbell oscillates periodically with small amplitude around the direction to the gravitational centre, and this periodic motion is stable, see Fig. 10a. Increasing energy we can achieve that besides quasi-periodic oscillations rotations are possible. Moreover, visible chaos appears, but still the periodic solution with small amplitude is stable. However, when the energy is big enough, then it vanishes and a new hyperbolic periodic solution with small amplitude appears.
These numerical examples suggest moreover that the original system with three degrees of freedom is not integrable. In fact, if the system restricted to \({{\mathscr {M}}}_4\) is not integrable, then the original system is not integrable. We tried to prove directly that the restricted system is not integrable using differential Galois tools, however, in vain because we were not able to find a non-equilibrium particular solution.
5 Proof of Theorem 2.1
We consider the complexification of the considered system in order to apply differential Galois methods to the integrability analysis. Applied notions and results of this theory are shortly described in “Appendix”.
Hamilton equations (2.19) have complex invariant manifold given by
On manifold \({{\mathscr {M}}}\), there is a family of particular solutions defined as solutions of the following linear inhomogeneous system
It has first integral
and its solutions are
Let \([\Phi ,\Theta ,\Xi ,P_{\varphi },P_{\theta },P_{\xi }]^T\) denote variations of variables \([\varphi ,\theta ,\xi ,p_{\varphi },p_{\theta },p_{\xi }]^T\). Then, the variational equations along this particular solution will take the form
where \(\xi =\xi (t)\) satisfies (5.1).
We notice that variational equations separate into two blocks: normal variational equations for variables \([\Phi ,\Theta ,P_{\varphi },P_{\theta }]^T\) and tangential equations for variables \([\Xi ,P_{\xi }]^T\).
We will now analyse normal variational equations. This system of four equations can be rewritten as the fourth-order equations, e.g. for variable \(\Phi \). We choose a particular solution corresponding to \(h=0\) just for simplification of expressions. After the change of independent variable \(t\rightarrow z=\xi (t)\), we obtain the following equation
where \('=\frac{\mathrm {d}}{\mathrm {d}z},\) with coefficients
The companion matrix of this equation is
see Eq. (A.5) for the definition of the companion matrix. We check whether fourth-order differential operator defined by Eq. (5.2) is symplectic. To this aim, we look for a skew-symmetric matrix
which satisfies the matrix equation
see Combot and Sanabria [1]. Common factor \(w_0\) of all entries of this matrix was distinguished just for simplifying the result form. Equation (5.4) has a solution with the following entries
Because of the presence of the square root \(\sqrt{z(\gamma z+z+2)}\) matrix W defines a projective symplectic structure. This structure is not-degenerated because
This implies that differential Galois group \({{\mathscr {G}}}\) of differential Eq. (5.2) is a subgroup of projectively symplectic matrices \(\mathrm {PSP}(4,\mathbb {C})\). In our further considerations, we use positivity of parameter \(\gamma \), \(\gamma >0\).
Our proof is based on the necessary conditions formulated in Lemma A.5 contained in “Appendix”. Thus, at first, we show that Eq. (5.2) does not have any hyperexponential solution. Function f(z) is called hyperexponential if its logarithmic derivative \(f'(z)/f(z)\) is a rational function. If \(\gamma \in \mathbb {R}^{+}\setminus \{0,\tfrac{3}{5},1\}\), Eq. (5.2) has five regular singularities: \(z_1=-1\), \(z_2=\tfrac{3}{\gamma -3}\), \(z_3=0\), \(z_4=-\tfrac{2}{\gamma +1}\) and \(z_5=\infty \). The sets of exponents at the respective singularities are following
where \(c_i\) are roots of the following polynomial of degree 4
They are given by
Equation (5.2) is Fuchsian; thus, it admits a hyperexponential solution if and only if it has the form \(P(z)\prod _i(z-z_i)^{e_i}\), where \(P(z)\in \mathbb {C}[z]\), \(z_i\in \mathbb {C}\) are singularities, \(e_i\) are exponents at \(z_i\), and there exists an exponent at infinity \(e_{\infty }\) such that the sum \(\sum _ie_i+e_{\infty }\) is a non-positive integer, see, e.g. Singer and Ulmer [12]. Taking into account exponents at finite singularities one of roots of (5.5) must be a half-integer. When we substitute \(Z=\tfrac{n}{2}\) to Eq. (5.5), then we obtain quadratic equation for \(\gamma \)
In order to have at least one real root, the discriminant of this equation should be non-negative that gives condition
As \(\lim _{n\rightarrow \pm \infty } Q(z)=-\infty \), function Q(n) takes non-negative values only for a few integer n contained in the interval \(\left[ \tfrac{2}{5}(5-3\sqrt{5}),\tfrac{2}{5}(5+3\sqrt{5})\right] \approx [-0.683,4.683]\), that is for \(n\in \{0,1,2,3,4\}\). For these values, we calculate corresponding values of \(\gamma \). Taking into account that \(\gamma >0\), we obtain that
For \(\gamma =\tfrac{1}{9}(29+2\sqrt{217})\), set of exponents at infinity is
and for \(\gamma =3\)
If \(\gamma \) is different from the above-mentioned values, then the differential Eq. (5.2) does not admit a hyperexponential solution. On the other hand, for these specific values of \(\gamma \) the equation does not depend on parameters and using a computer algebra system, e.g. Maple, we can check that it does not have any hyperexponential solution.
To check the second assumption of Lemma A.5, we can transform the system (5.4) to an equation of sixth order and check how many hyperexponential solutions it has. For computation simplifications, it is better to calculate the second exterior power of Eq. (5.2) and check how many hyperexponential solutions it admits. The second exterior power of Eq. (5.2) is a sixth-order differential equation with eight regular singularities: \(z_1=-1\), \(z_2=0\), \(z_3=-\tfrac{2}{\gamma +1}\), and \(z_4,z_5,z_6,z_7\) which are roots of the fourth-order equation
and \(z_8=\infty \) provided
where \(\gamma _i\) for \(i=1,\ldots ,11\) are roots of the following equation
The respective sets of exponents at singularities are the following:
where \(d_i\) are roots of equation
For the second exterior power to have exponential solutions, the sum \(\sum _{i=1}^8e_i\), \(e_i\in E_i\) should be a non-positive integer. We have one choice of exponents
\((e_1,e_2,e_3,e_4,e_5,e_6,e_7,e_8)=(-2,-1/2,-1/2,0,0,0,0,3)\) that sum to 0. Really, the second exterior power has a solution for every \(\gamma \) of the form
In order to find other solutions, one can built them only using a root of Eq. (5.9). Looking on exponents at finite points, we deduce that such a root must be a half-integer. Substitution \(Z=\tfrac{n}{2}\) to (5.9) gives the following second-order equation for \(\gamma \)
The discriminant of this polynomial is \(1024 ((n-12) n ((n-12) n+61)+964)\ge 0\); thus, it has real roots \(\gamma _1\) and \(\gamma _2\). Both of them are negative if simultaneously two inequalities due to Vieta’s formulas hold
This happens for integer \(n\in \mathbb {R}\setminus [2,10]\). Integers \(n\in [2,10]\) give the following positive values of \(\gamma \)
It means that for
where \(\gamma _1,\ldots ,\gamma _{11}\) satisfy Eq. (5.8), the differential Galois group contains the whole \(\mathrm {SP}(4,\mathbb {C})\) that is not solvable and thus in particular is not Abelian. These selected values of \(\gamma \) correspond to confluence of singularities of normal variational Eq. (5.2) or of its second exterior power, or possibility of existence of an exponential solution of normal variational equation or of its the second exterior power. For these values \(\gamma \), the system no longer depends on parameters and can then be treated with symbolic software as, e.g. Maple. We obtain that for all these selected cases normal variational equation does not factorise thus its differential Galois group is not reducible. The second exterior power also does not factorise and has just one hyperexponential solution even for these special values of \(\gamma \). Thus, differential Galois group of (5.2) is not irreducible. Therefore, the system is not integrable for all values \(\gamma >0\).
Data availability
Data sharing is not applicable to this article as no new data were created or analysed in this study.
References
Combot, T., Sanabria, C.: A symplectic Kovacic’s algorithm in dimension 4. In: ISSAC’18—Proceedings of the 2018 ACM International Symposium on Symbolic and Algebraic Computation, pp. 143–150. ACM, New York (2018)
Compoint, É., Weil, J.A.: Absolute reducibility of differential operators and Galois groups. J. Algebra 275(1), 77–105 (2004)
Kaplansky, I.: An Introduction to Differential Algebra, 2nd edn. Hermann, Paris (1976)
Karabanov, A.A., Morozov, A.D.: On resonances in Hamiltonian systems with three degrees of freedom. Regul. Chaotic Dyn. 24(6), 628–648 (2019)
Kovacic, J.J.: An algorithm for solving second order linear homogeneous differential equations. J. Symbolic Comput. 2(1), 3–43 (1986)
Maciejewski, A.J., Przybylska, M.: Non-integrability of the problem of a rigid satellite in gravitational and magnetic fields. Celestial Mech. Dynam. Astronom. 87(4), 317–351 (2003)
Magid, A.R.: Lectures on Differential Galois Theory, University Lecture Series, vol. 7. American Mathematical Society, Providence, RI (1994)
Morales, R., Juan, J.: Differential Galois Theory and Non-integrability of Hamiltonian Systems. Progress in Mathematics, vol. 179. Birkh\(\ddot{a}\)user Verlag, Basel (1999)
Schwarz, F.: A factorization algorithm for linear ordinary differential equations. In: Proceedings of the ACM-SIGSAM 1989 International Symposium on Symbolic and Algebraic Computation, ISSAC’89, pp. 17–25. Association for Computing Machinery, New York (1989)
Sidorenko, V.V., Celletti, A.: A “spring-mass” model of tethered satellite systems: properties of planar periodic motions. Celestial Mech. Dynam. Astronom. 107(1–2), 209–231 (2010)
Singer, M.F., Ulmer, F.: Galois groups of second and third order linear differential equations. J. Symbolic Comput. 16(1), 9–36 (1993)
Singer, M.F., Ulmer, F.: Necessary conditions for Liouvillian solutions of (third order) linear differential equations. Appl. Algebra Engrg. Comm. Comput. 6(1), 1–22 (1995)
Ulmer, F.: Liouvillian solutions of third order differential equations. J. Symb. Comput. 36(6), 855–889 (2003)
van Hoeij, M., Ragot, J.-F., Ulmer, F., Weil, J.-A.: Liouvillian solutions of linear differential equations of order three and higher. J. Symb. Comput. 28, 589–609 (1998)
Acknowledgements
The research has been partially supported by the National Science Centre of Poland under Grant No. 2020/39/D/ST1/01632.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix: Theory
Appendix: Theory
Important question concerning equations (2.19) is their integrability. Because they are Hamiltonian equations, thus the most natural is the notion of the integrability in the Liouville sense.
In the nineties of the previous century, a very effective theory has appeared that formulates obstructions to the integrability in the Liouville sense in the language of the differential Galois theory of variational equations obtained as a result of linearization of Hamiltonian equations in a neighbourhood of a certain non-equilibrium particular solution.
Roughly speaking, differential Galois group of variational equations is a linear algebraic group which preserves polynomial relations between solutions of these linear equations. As a linear algebraic group in general, it consists of several connected components in the Zariski topology, and among them the one containing identity is called the identity component. For a precise definition of the differential Galois group and differential Galois theory, see, e.g. Kaplansky [3] or Magid [7].
The presence of n commuting first integrals of a Hamiltonian system with n-degrees of freedom imposes a very special property of the differential group of variational equation that gives the necessary integrability condition formulated in Morales Ruiz [8].
Theorem A.1
Assume that a Hamiltonian system is meromorphically integrable in the Liouville sense in a neighbourhood of an analytic phase curve \(\Gamma \). Then, the identity component of the differential Galois group of the variational equations along \(\Gamma \) is Abelian.
Applications of this theorem are hard because of difficulties related to
-
finding a non-equilibrium particular solution,
-
determination of differential Galois group of variational equations. Difficulty grows very quickly with dimension of the variational equations.
In applications, it happens that variational equations split into tangential and normal subsystems and one can restrict to analysis of differential Galois group of normal variational equations that have smaller dimension than full variational equations, for detailed explanation, see, e.g. in Morales Ruiz [8] or shortly in Maciejewski and Przybylska [6]. In the case, when the dimension of normal variational equations is two and coefficients of these equations are rational functions of an independent variable (sometimes after appropriate algebraic change of independent variable) there is a very powerful tool for determination of the differential Galois group. This tool called the Kovacic algorithm formulated in Kovacic [5] is a decision procedure in a finite number of steps if one can write a general solution of differential equation
in a closed form, more precisely in a field of Liouvillian functions. This algorithm is constructive, it means that it enables to determine explicitly the form of solutions. As the byproduct, it determines the differential Galois group of Eq. (A.1). In fact, it is built on the base of complete classification of algebraic subgroups of the group \( \mathrm {SL}(2,\mathbb {C})\).
Unfortunately, there is no an equivalent of the Kovacic algorithm for linear differential equations with rational coefficients of higher orders although many partial results especially for equations of the third and fourth orders are known, see, e.g. Singer and Ulmer [11, 12], van Hoeij et al. [14], Ulmer [13] and references therein.
But in the case of Hamilton equations that can be written as
where \(\mathbb {I}_n\) is n-dimensional identity matrix, their variational equations in a neighbourhood of a particular solution \(\varphi (t)\) are also Hamiltonian
Here, variables \(\varvec{Y}\) are variations of \(\varvec{z}\). Solutions of these variational equations have very special property. Namely, if \(\varvec{y}_1\) and \(\varvec{y}_2\) are arbitrary two solutions, then quantity \((\varvec{y}_1,{\mathbb {J}}\varvec{y}_2)\) is constant. But this implies that for arbitrary element G of differential Galois group holds
and this implies \({\mathbb {J}}=G^T{\mathbb {J}}G\). But this equality means that the differential Galois group of VEs is a subgroup of \( \mathrm {Sp}(2n,\mathbb {C})\). Group \(\mathrm {Sp}(2,\mathbb {C})\) is isomorphic to \(\mathrm {SL}(2,\mathbb {C})\) but in higher dimensions \(\mathrm {Sp}(2m,\mathbb {C})\subset \mathrm {SL}(2m,\mathbb {C})\), for \(m>1\) and \(\mathrm {Sp}(2m,\mathbb {C})\) is much smaller than \(\mathrm {SL}(2m,\mathbb {C})\), thus classification of its subgroups should be much simpler. Using this classification, one can try to construct the equivalent of the Kovacic algorithm for symplectic differential operators and for dimension four this was formulated by Combot and Sanabria [1].
To describe it shortly, we introduce appropriate terminology. Let L be differential operator with coefficients in \(\mathbb {C}(z)\)
and A is its corresponding companion matrix
The 2n-th-order operator is
-
symplectic if there exists an invertible skew-symmetric matrix W with elements in \(\mathbb {C}(z)\) which is a solution of following equation
$$\begin{aligned} A^TW+WA+W'=0, \end{aligned}$$(A.6) -
projectively symplectic if there exists an invertible skew-symmetric matrix W with elements in \(\mathbb {C}(z)\) which is a solution of equation
$$\begin{aligned} A^TW+WA+W'+\lambda W=0, \end{aligned}$$(A.7)for a certain \(\lambda \in \mathbb {C}(z)\).
An operator L of order 2n is symplectic (respectively, projectively symplectic) when its Galois group is isomorphic to a subgroup of symplectic matrices \(\mathrm {Sp}(2m,\mathbb {C})\) (respectively, of projectively symplectic matrices \(\mathrm {PSp}(2m,\mathbb {C})\)) defined as
If L is projectively symplectic, then up to a multiplication of a hyperexponential function, one can ensure that the operator is symplectic. Function f(z) is called hyperexponential if its logarithmic derivative \(f'(z)/f(z)\) is a rational function.
Lemma A.1
Assume that the system \(\dot{x}=Ax \) is symplectic, i.e. there exists an invertible skew-symmetric matrix W with coefficients in \(\mathbb {C}(z)\) which is a solution of Eq. (A.6). Then, x(t) is a solution of equation \(\dot{x}(t)=Ax(t)\) if and only if \( x_{*}(t)=Wx(t)\) is a solution of the adjoint equation \({\dot{x}}_{*}(t)=-A^Tx_{*}(t)\).
This lemma can be proved by a direct check. It follows that if symplectic operator L has right factor \(L_{1}\), then its adjoint \(L^{*}\) has a right factor \({\widetilde{L}}_{1}\) of the same order. Thus, L has right factor \(L_{1}\), then it has also a left factor of the same degree.
With a system \(\dot{x}=Ax, \) we can associate its external second power. It is a system of the form
for an antisymmetric matrix W, see “Appendix” in [2]. Thus, Eq. (A.6) is an equation for the external square of a dual to the system \(x'=Ax\).
The classification theorem formulated in Combot and Sanabria [1] is the following.
Lemma A.2
A Lie subgroup of \(\mathrm {Sp}(4,\mathbb {C})\) is up to conjugacy generated by elements of the form:
-
1.
upper block triangular matrices with diagonal blocks of size at most \(2\times 2\),
-
2.
\(2\times 2\) diagonal matrices and anti-diagonal matrices
$$\begin{aligned} \begin{bmatrix} *&{}\quad *&{}\quad 0&{}\quad 0\\ *&{}\quad *&{}\quad 0&{}\quad 0\\ 0&{}\quad 0&{}\quad *&{}\quad *\\ 0&{}\quad 0&{}\quad *&{}\quad *\end{bmatrix},\quad \begin{bmatrix} 0&{}\quad 0&{}\quad *&{}\quad *\\ 0&{}\quad 0&{}\quad *&{}\quad *\\ *&{}\quad *&{}\quad 0&{}\quad 0\\ *&{}\quad *&{}\quad 0&{}\quad 0 \end{bmatrix}, \end{aligned}$$ -
3.
full group \(\mathrm {Sp}_{4}(\mathbb {C})\).
This classification is constructed on the base of known classification of Lie subgroups of wider unimodular group \(\mathrm {SL}(4,\mathbb {C})\supset \mathrm {Sp}(4,\mathbb {C})\). Subgroups of the projective symplectic group are central extensions of these and so contain multiples of the identity matrix with non-unit determinant. But as these commute with all matrices, the possible structures of subgroups in items 1, 2 are unchanged.
The next two lemmas characterise reducible case.
Lemma A.3
Let us consider the following block diagonal system
where B, C and D are \(2\times 2\) matrices with rational coefficients. Then, Eq. (A.6) has the following particular solution
where r is a rational function.
Lemma A.4
Let us consider differential operator
where \(a,b,c,d\in \mathbb {C}(z)\) and let A be its companion matrix. Then, Eq. (A.6) has the particular solution
and \(\det W=0\).
Both lemmas can be proved by a direct check.
In our proof of the non-integrability, we will use the following criterion.
Lemma A.5
Assume that equation
is projectively symplectic and A is its corresponding companion matrix. If Eq. (A.13) does not admit a hyperexponential solution and Eq. (A.6) has exactly one hyperexponential solution, then the differential Galois group of (A.13) contains \(\mathrm {Sp}(4,\mathbb {C})\).
Proof
If (A.13) does not admit a hyperexponential solution, then it does not admit a right factor of order one. Thus, if it is reducible, it admits a right factor of order two. Hence, using Lemma A.4, Eq. (A.6) has a hyperexponential solution such that the corresponding antisymmetric matrix \(W_1\) is singular. Equation (A.13) is projectively simplectic, and thus, Eq. (A.6) admits a hyperexponential solution such that the corresponding antisymmetric matrix \(W_2\) is not singular. Thus, \(W_1,W_2\) cannot be the same, but this is impossible as by assumption (A.6) has exactly one hyperexponential solution. It implies that Eq. (A.13) is irreducible. If it admits Liouvillian solutions, then it admits at least two projective Poisson structures. But this is impossible as by assumption there exists exactly one hyperexponential solution of (A.6). \(\square \)
Item 1 in Lemma A.2 corresponds to reducible cases. Differential Galois group of equation
is a reducible subgroup of \(\mathrm {Sp}(4,\mathbb {C})\) if and only if the equation itself can be factorised. If a symplectic operator has a left factor, then it has a right factor of the same order. Thus for an order four symplectic operator only the following factorizations into factors with appropriate orders are possible: (i) 1,1,1,1 or (ii) 1,2,1 or (iii) 2,2.
If this equation has a factor of order one, then it is either a right factor of order one or a left factor of order one. Since a left factor of order one leads to a right factor of order one of the adjoint differential operator \(L^{*}(y)\), then one can restrict to testing of a right factor of order one of L(y) or of \(L^{*}(y)\). If an equation has a right factor of order one, then it has an hyperexponential solution. If \(L(y) = 0\) is of Fuchsian type, then any hyperexponential solution must be of the form \(P(z)\prod _i(z-z_i)^{e_i}\), where \(P(z)\in \mathbb {C}[z]\), \(z_i\in \mathbb {C}\) are singularities except the infinity, \(e_i\) are exponents at \(z_i\). Then, using lemma 3.1 in Singer and Ulmer [12] about the necessary condition for such a solution one can formulate the following proposition
Proposition A.1
If a Fuchsian symplectic equation of order four has a factor of order one, then \(L(y)=0\) has a solution of the form \(P(z)\prod _i(z-z_i)^{e_i}\), where \(P(z)\in \mathbb {C}[z]\), \(z_i\in \mathbb {C}\) are singularities, \(e_i\) are exponents at \(z_i\), and there exists an exponent at infinity \(e_{\infty }\) such that the sum \(\sum _ie_i+e_{\infty }\) is a non-positive integer.
In the case when L(y) factorises into two factors \(L=L_1L_2\) of order two to each of them one can use the Kovacic algorithm. If one of them is not solvable, then the differential Galois group is not solvable.
Item 2 in Lemma A.2 corresponds to irreducible solvable cases. Then, operator L admits a factorisation in two operators of order two with coefficients in a quadratic extension of \(\mathbb {C}(z)\). This case can be identified by the presence of two linearly independent rank two Poisson structures (requirement as for symplectic structure but without invertibility requirement) with coefficients in a quadratic extension of \(\mathbb {C}(z)\), see Combot and Sanabria [1]. But the presence of two structures is related to the presence of two hyperexponential solutions of the exterior power of L. Thus, to identify the irreducible solvable cases we use the following proposition formulated in Combot and Sanabria [1]
Proposition A.2
The necessary condition for irreducible solvable cases is that the second exterior power of differential equation has at least two hyperexponential solutions.
The second exterior power of a differential operator L with solutions \(\{y_1,\ldots ,y_n\}\) is the equation having for its space of solutions the vector space generated by the \(2\times 2\) Wronskians of solutions \(\begin{vmatrix}y_i&y_j\\ y_i'&y_j'\end{vmatrix}\), \(i<j\) of operator L, see, e.g. Schwarz [9]. To identify these solutions for Fuchsian operators, L one can use Proposition A.1.
Item 3 in Lemma A.2 corresponds to non-solvable cases.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Combot, T., Maciejewski, A.J. & Przybylska, M. Non-integrability of a model of elastic dumbbell satellite. Nonlinear Dyn 106, 125–146 (2021). https://doi.org/10.1007/s11071-021-06771-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11071-021-06771-y