A Geometric Approach to the Sundman Transformation and Its Applications to Integrability

: A geometric approach to the integrability and reduction of dynamical systems, both when dealing with systems of differential equations and in classical physics, is developed from a modern perspective. The main ingredients of this analysis are infinitesimal symmetries and tensor fields that are invariant under the given dynamics. A particular emphasis is placed on the existence of alternative invariant volume forms and the associated Jacobi multiplier theory, and then the Hojman symmetry theory is developed as a complement to the Noether theorem and non-Noether constants of motion. We also recall the geometric approach to Sundman infinitesimal time-reparametrisation for autonomous systems of first-order differential equations and some of its applications to integrability, and an analysis of how to define Sundman transformations for autonomous systems of second-order differential equations is proposed, which shows the necessity of considering alternative tangent bundle structures. A short description of alternative tangent structures is provided, and an application to integrability, namely, the linearisability of scalar second-order differential equations under generalised Sundman transformations, is developed.


Introduction: The Meaning of Integrability by Quadratures
Obtaining the solutions of autonomous systems of differential equations appearing in many physical problems is not an easy task.In geometric terms, they are represented by vector fields, whose integral curves are the solutions of the system.
In fact, in order to study an autonomous system of first-order differential equations, use is made of symmetry and reduction techniques [1].The main parts of these techniques are based on the transformation of the given system into a related one.From a local point of view, sometimes, transformations are seen as coordinate changes (passive point of view).Or, more explicitly, a diffeomorphism F : M / / M of a differentiable manifold changes each local chart (U, ϕ) into a new one (F(U), ϕ • F −1 ).
In order to find their solutions, i.e., the flow of the vector fields, as in a generic case, there is no way of writing them in an explicit way, i.e., using fundamental functions, and we are happy if, at least, we can express the solutions in terms of quadratures.More specifically, integrability by quadratures means that one can determine the solutions (i.e., the flow of X) by means of a finite number of algebraic operations (including inversion of functions) and the quadratures of some functions [2].
The characterisation of integrable systems admitting such a type of solution has received a lot of attention, and the answer is always based on the use of Lie algebras of vector fields containing the given one [2][3][4].In order to study such systems, use is made of symmetry and reduction techniques.The knowledge of some particular solutions of them or related systems may be useful.But symmetry and constants of motion emerge from the existence of compatible geometric structures provided by special tensor fields.
Both concepts, reduction and integrability, appear very often in physics and mathematics and with very different meanings.For instance, students in physics relate the reduction process to infinitesimal symmetries or constants of motion, generally in the framework of Lagrangian or Hamiltonian mechanics, because several reductions in the context of symplectic geometry are well known.Note, however, that reduction processes can be developed in a more general framework and need not be related to symmetry properties, but these only aim to find simpler solvable systems providing at least partial information on the solution of the original problem.The existence of compatible geometric structures on M provided by special tensor fields is also very useful.
The problems in which reduction processes are relevant are the following ones: • Given a problem that is difficult to solve, we can look for related simpler problems whose solutions provide at least partial information on those of the first one.• Sometimes, we find geometric structures difficult to handle, for instance, because they have some kind of degeneracy.We then look for a related structure of the same kind that is easier to manage, for instance, a nondegenerate one.

•
When the original problem admits some geometric structures, then the reduction procedure tries to preserve such structures.
These are the types of problems we are going to deal with.It is remarked that there is no systematic way of selecting the related reduced systems, and sometimes, there are different possibilities.Of course, in order for the method to be efficient, we should have, in advance, some simple problems for which we know the solutions.
The problem to be analysed here is that of (systems of) first-order differential equations.The higher-order case can be related to a first-order one.
In principle, one tends to assume that a system with fewer unknown variables or a lower-order differential equation is simpler than one with more unknown variables or a higher-order differential equation.The reducibility to a system solvable by quadratures is the final objective.Sometimes, the problem of reduction is related to the so-called coordinate separability.
Usually, reduction processes starting from linear systems may produce nonlinear interacting ones, which, from the physical viewpoint, means that interaction terms can be provided by reduction [5].Conversely, given a nonlinear system, is it possible to obtain it as a reduction from a simpler linear one?

The Reduction of a System of Ordinary Differential Equations
If F : M / / N is a differentiable map, the vector field X ∈ X(M) is F-related to the vector field Y ∈ X(N) if, for each point p ∈ M, F * p (X p ) = Y F(p) , i.e., TF • X = Y • F. We say that X is F-projectable.The following diagram is commutative: Note that this implies that if c : I / / M is an integral curve of X, then F • c is an integral curve of Y (see, e.g., [5]).
The reduction of the dynamics given by the vector field X ∈ X(M) means finding a manifold N, a differentiable map F : M / / N and an F-related vector field Y ∈ X(N).
The integral curves of Y are F-projections of the integral curves of X.Another problem remains: the reconstruction of the integral curves of X ∈ X(M) from the integral curves of Y ∈ X(N).
In the reduction process, the cases to be considered are either when i : N / / M is an immersed submanifold of the n-dimensional manifold M or when F : M / / N is a surjective submersion.
(a) In the first case, the reduction implies that X |i(N) must be tangent to i(N), and therefore, there exists a vector field X on N such that it is i-related to X.Then, the integral curves of X provide us with the integral curves of X.We can suppose that the (n − k)submanifold N is (at least locally) defined by k functions Φ 1 , . . ., Φ k , which should be first integrals for X because of the assumed tangency condition.Note that only some integral curves of X will be obtained in this way, and we need a family of such submanifolds depending on k parameters in order to recover all integral curves of X ∈ X(M).
This procedure is very commonly used in mechanics: under very general conditions, determining k functionally independent constants of motion reduces the original problem to another problem involving only (n − k) variables but depending on k parameters.
(b) Finding a manifold N and a surjective submersion F : M / / N such that X is a projectable vector field is not an easy task.The map F defines an equivalence relation on M, and N is the space of equivalence classes.The starting point may be the equivalence relation, and then F is the projection map.Usually, one considers a transformation Lie group G of M that is a symmetry group for X and the equivalence relation associated with the action of G on M: we denote the space of orbits by N, assumed to be a differentiable manifold, with the function F being the natural projection of each element of M on its orbit.The symmetry condition implies that X is a projectable vector field.In this way, the symmetry condition of X is incorporated into the reduction process.
In other situations, the equivalence relation is defined by a foliation (involutive distribution), and the equivalence classes are the leaves.The set of leaves is assumed to have a manifold structure.If the original system has an invariance group, we can consider invariant foliations in such a way that the leaves are preserved.
The above reduction procedures can be used in an iterative way.In summary, there are two main techniques in the process of solving the system: • The determination of first integrals for X: First integrals provide us with foliations such that the vector field X is tangent to the leaves, and in this way, the problem reduces to that of a family of lower-dimensional problems, one on each leaf.

•
Search for symmetries of the vector field X: The knowledge of infinitesimal groups of symmetries of the vector field (i.e., of the system of differential equations) suggests that we use adapted local coordinates, with the system then decoupling into lowerdimensional subsystems.
More specifically, the knowledge of r functionally independent (i.e., such that dF 1 ∧ • • • ∧ dF r ̸ = 0) first integrals, F 1 , . . ., F r , allows us to reduce the problem to that of a family of vector fields X c , with c ∈ R r , defined on the n − r-dimensional submanifolds M c given by the level sets of the vector function of rank r, (F 1 , . . ., F r ) : M / / R r .
Of course, the best situation is when r = n − 1: the leaves are one-dimensional, giving us the solutions to the problem, up to a reparametrisation.
There is another way to reduce the problem.Given an infinitesimal symmetry (i.e., a vector field Y such that [Y, X] = 0) in a neighbourhood of a point where Y is different from zero, we can choose adapted coordinates, (y 1 , . . ., y n ), for which Y is written as (Straightening out Theorem, see, e.g., [6]) Then, the symmetry condition [Y, X] = 0 implies that X has the form and its integral curves are obtained by solving the system We have reduced the problem to a subsystem involving only the first n − 1 equations, and once this has been solved, the last equation is used to obtain the function y n (t) by means of one quadrature.Note that the new coordinates, y 1 , . . ., y n−1 , are first integrals for Y, and therefore, we cannot easily find such coordinates in a general case.
Moreover, the information provided by two different symmetry vector fields cannot be used simultaneously in the general case, because it is not possible to find local coordinates (y 1 , . . ., y n ) such that In terms of adapted coordinates for X, i.e., such that X = ∂/∂y n , the integration is immediate, the solution being This proves that the concept of integrability by quadratures depends on the choice of initial coordinates.In adapted coordinates, it is always integrable.However, it can be proved that when X is part of a family of vector fields satisfying appropriate conditions, it is integrable by quadratures for any choice of initial coordinates (see, e.g., [2,3]).
We also remark that the set is an R-linear space, and that the flows of such vector fields Y ∈ X X (M) preserve the set of integral curves of X.
Both constants of motion and infinitesimal symmetries can be used simultaneously if some compatibility conditions are satisfied.We can say that a system admitting r < n − 1 functionally independent first integrals for X, F 1 , . . ., F r , is integrable when we know that there are s commuting infinitesimal symmetries Y 1 , . . ., Y s , with r and The first integrals determine an s-dimensional foliation, and the latter conditions mean that the restrictions of the s vector fields Y a to the leaves are tangent to such leaves.
Sometimes, we have additional geometric structures that are compatible with the dynamics.For instance, an interesting example to be studied in Section 4 is that of a symplectic structure ω on a 2n-dimensional manifold M. Such a 2-form relates, by contraction, vector fields and 1-forms in a one-to-one manner.The vector fields X F associated with exact 1-forms dF are said to be Hamiltonian vector fields.The compatibility of ω means that the vector field X is a Hamiltonian vector field X = X H .The regular Poisson bracket defined by ω (i.e., {F 1 , F 2 } = X F 2 F 1 ) allows us to express the tangency conditions as i.e., the first integral functions F a must be in involution, and the corresponding Hamiltonian vector fields commute.This leads to the particularly interesting Arnold-Liouville definition of (Abelian) complete integrability (r = s = n).The vector fields are X a = X F a and, for instance, F 1 = H.
In summary, the problem of integrability by quadratures depends on determining, by quadratures, the necessary first integrals and infinitesimal symmetries and finding adapted coordinates, or, in other words, as shown in the next sections, finding a sufficient number of invariant tensors.
An appropriate set of infinitesimal symmetries of the vector field X provides constants of motion.More specifically, let {X 1 , . . ., X d } be a set of d vector fields that take linearly independent values at every point and that are infinitesimal symmetries of the dynamical vector field X.If they generate an involutive distribution (i.e., there exist functions , then, for each triple of numbers i, j, k, the functions f ij k are first integrals for X, i.e., X( f ij k ) = 0.In fact, the Jacobi identity for the vector fields X, X i , X j , i.e., [[X, X i ], Moreover, for any other index l, X l ( f ij k ) is also a first integral for X, because as X l is The constants of motion so obtained are not functionally independent, but at least this proves the usefulness of finding these families of vector fields when looking for first integrals.

Invariant Tensor Fields and the Integrability of a Vector Field
The main point is that the integrability of a given vector field X on a manifold M is usually studied by means of its compatible geometric structures.Geometric structures are defined by tensor fields t, and they are said to be compatible with X when they are invariant under X, i.e., L X t = 0.
As an example mentioned before, a function F ∈ C ∞ (M) invariant under X defines a constant of motion, also called the first integral for X.The corresponding exact 1-form dF defines an (n − 1)-dimensional foliation (integrable distribution), and the restriction of the vector field X to each leaf is tangent to the leaf, thus defining a family of vector fields on lower-dimensional manifolds, which can be considered reduced systems.This property, where a first integral, i.e., an invariant function, can be used to reduce the integrability of the vector field to a family of lower-dimensional problems, is a particular example of a more general case: the mentioned role of tensor fields invariant under a vector field X on an n-dimensional manifold M in its integrability (see, e.g., [4]).
Another example that we have also considered is that of invariant vector fields Y ∈ X(M), i.e., such that Y ∈ X X (M) = {Y ∈ X(M) | [Y, X] = 0}, then giving rise to infinitesimal symmetries of X, whose flows transform solutions into solutions.As the flows of such vector fields Y ∈ X X (M) preserve the set of integral curves of X, we can also consider a set of vector fields Y ∈ X(M) whose flows preserve each integral curve of X up to reparametrisations; i.e., they preserve the one-dimensional distribution generated by X, and then the vector field Y is a symmetry of the one-dimensional distribution spanned by X: there exists a function h ∈ C ∞ (M) such that [Y, X] = h X.Such a set of vector fields generating flows preserving the set of integral curves of X up to a reparametrisation is also a real linear space containing X X (M), usually denoted by One can check that X X (M) is a real Lie algebra containing X X (M) as a real Lie subalgebra, but X X (M) is not an ideal in X X (M).
Moreover, if F is a first integral for X, then when Y is a symmetry of the vector field X or, more generally, of the one-dimensional distribution spanned by X, the function Strictly positive symmetric forms g on M invariant under X ∈ X(M) imply that the given vector field is a Killing vector field whose integrability we have additional information on (see, e.g., [7]), while a nondegenerate closed 2-form ω on M is invariant under X ∈ X(M), and such a vector field is a locally Hamiltonian vector field on the symplectic manifold (M, ω).This is a particularly interesting case that will be studied in the next section.
Later on, we will see several examples of how to find new first integrals when additional tensor fields enter the game.Of course, some invariant multivector fields are also interesting, but the case of Poisson bivector fields is particularly relevant.Differential forms that are invariant under X give rise to absolutely integral invariants, a theory developed by Poincaré in [8].In particular, invariant 2-forms, for instance, the relevant case in the mechanics of symplectic forms, are very relevant in the study of the integrability of X, and Arnold-Liouville integrability is based on the existence of an appropriate invariant symplectic form.
We now remark that when we have a (1, 1)-tensor field R invariant under X, it may be used as a generator of symmetry, in the sense that if Y is a symmetry of X, then L X R = 0 implies that L X (R(Y)) = 0, and in this sense, R generates symmetry; for instance, as X itself is a symmetry, the vector fields R k (X), k ∈ N, are infinitesimal symmetries, too.
Furthermore, if we know such an invariant (1,1)-tensor field R and choose a local basis of vector fields {X 1 , . . ., X n }, then if the matrices A and B are the matrix representatives of R and L X in such a basis, respectively, then the invariance condition L X R = 0 is written as where Ȧi j denotes X(A i j ), because if R(X i ) = A i j X j and L X X i = B i j X j , as we have and, therefore, to the matrix Equation ( 2), which is called the Lax equation [9][10][11].Also, note that as the matrix The importance of these equations is that, as for any pair of matrices, the trace of the commutator is zero, and we have d dt TrA k = 0, and consequently, a (1, 1)-tensor field R invariant under X provides us with n constants of motion; in principle, not all of them are functionally independent, as indicated by the Hamilton-Cayley theorem.Another remark is that, as the coefficients of the characteristic equation det(A − λI) = 0 can be reconstructed from the traces of powers of A (the Le Verrier method of determining the characteristic equation of a matrix; see, e.g., [12]), the roots of such a characteristic equation, the eigenvalues of A, are constants of motion.
The search for these invariant (1,1)-tensor fields R is not an easy task and used to come from the existence of alternative structures [5,[13][14][15], with the associated constants of motion being called non-Noether constants of motion (see, e.g., [13,16,17]), but they may have a different origin [10,18].We will come back to this point in Section 5.
Another interesting case is when there exists a volume form invariant under X.This case and an alternative method to find the first integrals associated with symmetries of a one-dimensional distribution spanned by the vector field X will be analysed in Section 5.

Invariant Symplectic Structures
We recall that (M, ω) is a symplectic manifold if M is a finite-dimensional differentiable manifold and ω is a nondegenerate 2-form that satisfies dω = 0 (that is, ω ∈ Z 2 (M)).
The dimension of M is necessarily even, dim M = 2n.General results can be found in reference textbooks, for instance [19,20].By the nondegeneracy of ω, we mean that for every point u ∈ M, the map ω u : has a maximal rank, i.e., (ω) ∧n ̸ = 0.Such a map ω is a base-preserving fibred map, and hence, it induces a mapping between the linear spaces of sections, which, with a slight abuse of notation, we will also write as ω : X(M) / / 1 (M).
The following well-known result completely characterises these symplectic manifolds from the local point of view: Such coordinates are called Darboux coordinates.Since closed and, in particular, exact 1-forms are distinguished elements in 1 (M), the corresponding vector fields are called locally Hamiltonian vector fields and Hamiltonian vector fields, respectively.If H ∈ C ∞ (M), the Hamiltonian vector field X H associated with the Hamiltonian H is the unique vector field satisfying ω(X H ) = i(X H )ω = dH.The set of Hamiltonian vector fields will be denoted by X H (M, ω), and that of locally Hamiltonian vector fields by Observe that ω −1 is an isomorphism of real vector spaces.
The Cartan homotopy identity, i.e., L X α = i(X)dα + d(i(X)α), for any α ∈ p (M), shows that, given X ∈ X(M), L X ω = 0 if and only if i(X)ω is a closed 1-form, i.e., X ∈ X LH (M, ω).In particular, L X H ω = 0.In Darboux coordinates, the Hamiltonian vector field X H corresponding to the function H is given by and therefore, the equations determining its integral curves are similar to Hamilton equations.
Consequently, the set X LH (M, ω) is a Lie algebra, and As an important property, if (M, ω) is a symplectic manifold, we can define the Poisson bracket of two functions F, G ∈ C ∞ (M) as being the function {F, G} given by In Darboux coordinates for ω, the expression for {F, G} is the usual one: The aforementioned property (4) This and taking into account that H}, as well as the corresponding expressions for each cyclic reordering, we find that Finally, we recall that at each point of a neighbourhood of every point of M, the values of the set of Hamiltonian vector fields generate the tangent space, and therefore, dω = 0 if and only if the Jacobi identity holds.Hence, the Poisson bracket endows C ∞ (M) with a real Lie algebra structure, with the Jacobi identity being a consequence of the closedness of ω, and property (7) A Hamiltonian system is a triplet (M, ω, H), where (M, ω) is a symplectic manifold and H ∈ C ∞ (M).Infinitesimal symmetries of a Hamiltonian system (M, ω, H) are vector fields X whose flows are symplectomorphisms, i.e., L X ω = 0, such that they are also symmetries of H, i.e., XH = 0, and, therefore, they are symmetries of X H , [X, It is remarked that, in the case of a Hamiltonian dynamical system, the constants of motion F play a dual role in the reduction process, either as constants of motion (i.e., first integrals for X H ) or as generating infinitesimal symmetries X F , because if F is a constant of motion, {F, H} = 0, then the Hamiltonian vector field X F is an infinitesimal symmetry of H.
Moreover, both X H and X F are tangent to the level sets of F, because X H F = {F, H} = 0, and X F F = {F, F} = 0.The restriction of X F on each leaf of the foliation F H , whose leaves are the level sets of H, can be used to determine adapted coordinates, because [X F , X H ] = X {H,F} = 0, in such a way that the problem of determining the integral curves of X H is reduced not just in one but also in two degrees of freedom.
In order to be able to simultaneously use the information given by different constants of motion, F 1 , . . ., F r , it is sufficient that {F i , F j } = 0, ∀i, j = 1, . . ., r, because [X F i , X F j ] = X {F j ,F i } .If the condition is satisfied, then [X F i , X F j ] = 0 for any pair of indices i and j, and we can find adapted coordinates such that X F i = ∂/∂y i for i = 1, . . ., r.The condition is not necessary because if, for instance, {F i , F j } is constant for each pair of indices, then it is also true that [X F i , X F j ] = 0.
The simplest case will be that of Hamiltonian systems in a space of dimension N = 2 n for which such n functionally independent first integrals, F 1 , . . ., F n , in involution are known: they are called completely integrable systems.
We are now ready to recall the notion of Liouville-Arnold integrability: A completely integrable system that admits more than n functionally independent first integrals is said to be superintegrable, and when the number is the maximum number, 2n − 1, the system is called maximally superintegrable [21].Many such systems can be found in the physics literature (see, e.g., the recent papers [22,23] and references therein).
The two main examples of Hamiltonian dynamical systems are those of Hamiltonian systems on the cotangent bundle T * Q of a manifold Q and those defined by regular Lagrangians on the tangent bundle TQ, to be studied in this section.In fact, the cotangent bundle π : The geometric framework for the study of Lagrangian mechanics is that of tangent bundles [5,19,20,24,25].The tangent bundle τ : TQ / / Q is characterised by two geometric tensors, the vertical endomorphism S, a (1,1)-tensor field on TQ, also called a tangent structure, which satisfies Im S = ker S and an integrability condition, the vanishing of the Nijenhuis tensor, N S = 0, and the Liouville vector field ∆ generating dilations along fibres in TQ [26].If (U, φ) is a local chart on Q and pr i : R n / / R are the natural projections on the i-th factor and q i = pr i • φ, we define the coordinate system (U, q 1 , . . ., q n ) on Q, and the corresponding chart in U = τ −1 (U), given by (U , φ, φ * ), defines a coordinate system (q 1 , . . ., q n , v 1 , . . ., v n ) on the open set U = τ −1 (U) of TQ.Correspondingly, we consider the coordinate basis of X(U), usually denoted by {∂/∂q j | j = 1, . . ., n}, and its dual basis for 1 (U), ).With this notation, the coordinate expressions of the vertical endomorphism S and the Liouville vector field ∆ are [24,25] Similarly, we can introduce a coordinate system (q 1 , . . ., q n , p 1 , . . ., p n ) in the open set Ū = π −1 (U) of the cotangent bundle π : T * Q / / Q, and the coordinate expression of the which shows that ω = −dθ is a canonical symplectic structure on T * Q.
The vector field ∆ and the tensor field S can be used to select a special kind of vector field whose integral curves are given by the lifting solutions of second-order differential equations.These vector fields, called second-order differential equation vector fields D ∈ X(TQ) (hereafter shortened to SODE vector fields), are characterised by S(D) = ∆.
Given a Lagrangian L ∈ C ∞ (TQ), we can define a 1-form θ L = dL • S and the exact 2form ω L = −dθ L .When ω L is regular, the Lagrangian L is said to be regular: i.e., (TQ, ω L ) is a symplectic manifold.Then, the dynamics are given by the uniquely defined SODE vector field Γ L such that where the energy function This shows that (TQ, ω L , E L ) is a Hamiltonian dynamical system.The contraction of both sides of the first equation in (10) shows that the energy function is a first integral.Moreover, this implies that In the usual local tangent bundle coordinates, the expressions of θ L and E L are while that of ω L is It may be of interest to know under what conditions two regular Lagrange functions L, L ′ ∈ C ∞ (TQ), in TQ produce the same symplectic structure and the same energy function; i.e., L 0 = L − L ′ is such that ω L 0 ≡ 0 and Then, a function L 0 ∈ C ∞ (TQ) is such that ω L 0 ≡ 0 if and only if there exist a closed 1-form α ∈ Z 1 (Q) and a function h ∈ C ∞ (Q) such that L 0 = α + h • τ (see, e.g., [27]).But, as for such a function L 0 , we have E L 0 = −h, we see that L, L ′ ∈ C ∞ (TQ) in TQ produces the same Hamiltonian dynamical system, i.e., ω L = ω L ′ and E L = E L ′ , if and only if there exists a closed 1-form α ∈ Z 1 (Q) such that L ′ = L + α, and both Lagrangians are then said to be gauge-equivalent.
It is also possible to identify a complete lift vector field X c ∈ X(TQ) that is a symmetry of the Hamiltonian dynamical system (TQ, ω L , E L ) in terms of symmetries of L itself.Recall that if X ∈ X(Q), its complete lift, denoted by X c , is the vector field on TQ, whose flow is given by ϕ t * , where ϕ t is the local flow of the vector field X.If the local coordinate expression of the vector field X in a chart of Q is X(q) = X i (q)∂/∂q i , then the corresponding expression for X c in the associated chart of TQ is while, as ∆ is Φ-related to itself, (∆L) At the infinitesimal level, ( 13) and ( 14) mean that if All of these properties can be used to derive Noether's theorem in the Lagrangian formalism: Proof.First, the vector field X c is Hamiltonian because and then and then if Γ L is the dynamical vector field determined by L, i.e., i(Γ L )ω L = dE L , we have Particularly interesting cases are the geodesic Lagrangians defined through a Riemann structure g by L = 1 2 g(v, v) (see, e.g., [7,28,29]) and the so-called natural Lagrangians defined by means of a Riemann structure g and a function V in Q as follows: L = 1 2 g(v, v) − V(q).A remarkable fact is the possibility of the existence of alternative compatible structures of the same type.For instance, a Hamiltonian vector field X H on the symplectic manifold (M, ω) can also be a Hamiltonian vector field with respect to an alternative symplectic structure ω ′ ; i.e., there exists H ′ ∈ C ∞ (M) such that i(X H )ω ′ = dH ′ .This also implies that L X H ω ′ = 0, and we see in this way that ω ′ is an alternative X H -invariant 2-form for X H .But if ω ′ is a 2-form such that L X H ω ′ = 0, then L X H (ω ′ − λ ω) = 0, λ ∈ R, also implies that L X H (ω ′ − λ ω) ∧n = 0, and as (ω ′ − λ ω) ∧n = f λ (ω) ∧n , we see that L X H f λ = 0, and therefore, the function coefficients of powers of λ are first integrals.They are usually called non-Noether first integrals [13,16].Such constants are the same ones as those associated with the invariant (1,1)-tensor field R = ω −1 • ω ′ in Section 3.
We now mention two particular methods for finding such an X H -invariant 2-form ω ′ .First, if a vector field Y is such that [Y, X H ] = 0 but is not a locally Hamiltonian vector field, then is a noncanonical symmetry of the Hamiltonian vector field defined by the Hamiltonian dynamical system (M, ω, H), from the relation i ) * H), and therefore, L X H ((ϕ −1 ) * (ω)) = 0, and equivalently, L X H (ϕ * (ω)) = 0.As it was assumed that ω ′ = ϕ * (ω) ̸ = ω, we can choose such an X H -invariant closed 2-form ω ′ to define, together with ω, the pencil.We also recall that ϕ ∈ Diff(M) is said to be a canonoid transformation of the Hamiltonian dynamical system (M, ω, H) when ϕ * (X H ) is also a Hamiltonian vector field or, equivalently, when X H is Hamiltonian with respect to the transformed 2-form ϕ * (ω), and in this case, L X H (ϕ * (ω)) = 0 [14,15]; i.e., we can also choose ω ′ = ϕ * (ω) as the invariant X H -form.

Invariant Volume Forms and Jacobi Multipliers
Another interesting case with many applications, both in the theory of differential equations and in classical mechanics, is that of the existence of a volume form invariant under a given vector field X on an oriented manifold (M, Ω).Each volume form is of the form R Ω, with R ∈ C ∞ (M), and the invariance condition is L X (R Ω) = 0, and if we take into account that L Also, recall that for each vector field X ∈ X(M), the volume form L X Ω is proportional to Ω, and the proportionality function is called the divergence of X, i.e., and if local coordinates are chosen such that Vector fields X such that div(X) = 0 are called divergence-free vector fields and possess interesting properties (see, e.g., [4]).
Recalling the properties of Lie derivatives, we find that, as for each Nonvanishing functions R such that L X (R Ω) = L R X (Ω) = 0, i.e., div(R X) = 0, are called Jacobi multipliers, and (17) shows that the nonvanishing function R is a Jacobi multiplier if and only if (see [30,31] and references therein) Locally defined Jacobi multipliers, obtained as particular solutions of ( 18), always exist, but in some particular cases, global solutions, giving rise to invariant volume forms, do not exist [32,33].
Note that the relation L )Ω) means that, for each positive function f , the function R is a Jacobi multiplier for X on the oriented manifold (M, Ω) if and only if f R is a Jacobi multiplier for X in the oriented manifold (M, f −1 Ω).Moreover, if two Jacobi multipliers are known, there exist two invariant volume forms, Ω 1 = R 1 Ω and Ω 2 = R 2 Ω, such that Ω 2 = (R 2 /R 1 )Ω 1 , and hence, as An interesting case is that of a SODE vector field Γ ∈ X(TQ) admitting a Lagrangian formulation L; i.e., there exists a regular Lagrange function L ∈ C ∞ (TQ) such that i(Γ)ω L = dE L or, equivalently, L Γ θ L = dL.The tangent bundle TQ is orientable, and a local chart of TQ induced by one on Q induces a (local) orientation.As L Γ ω L = 0, the volume form (ω L ) ∧n is Γ-invariant, and therefore, if a volume form Ω was previously fixed, there will be a Jacobi multiplier R such that (ω L ) ∧n = R Ω. If, for instance, Ω is the volume form determined by a tangent bundle local chart, we find, according to (12), that the determinant of the Hessian matrix W with elements [30] and references therein).Alternative Lagrangian formulations provide different Jacobi multipliers giving rise to constants of motion, which allows the system to be reduced.
The main result is very general and a simple consequence of a geometric relation that, when particularised to specific cases, gives rise to many results scattered throughout the physics literature.So, different relations among vector fields and their divergences can be used to establish first integrals and integral invariants for vector fields on a manifold M.
The following geometric relation plays a fundamental role: if X, Y is an arbitrary pair of vector fields on an oriented manifold (M, Ω), then because as the Lie derivatives of a p-form and from the difference between and the corresponding relation we obtain (19).Consequently, from relation (19), we see that if the vector field X ∈ X(M) on an oriented manifold (M, Ω) is divergence-free, then if the vector field Y is an infinitesimal symmetry of X, i.e., [X, Y] = 0, we find that div (Y) is a constant of motion for X.Something similar happens when the vector field Y is an infinitesimal symmetry of the one-dimensional distribution generated by X: i.e., there exists a function h such that [Y, X] = h X, because then, in this case, the function div (Y) + h is a first integral for X.Actually, if div (X) = 0, then L X Ω = 0, and hence, as we obtain L X (div (Y) + h) = 0, and therefore, the following function I is a constant of motion for X: When the vector field X ∈ X(M) is not divergence-free, because, for any nonvanishing function R, the constants of motion of X coincide with those of R X, if we choose R as a Jacobi multiplier for X, then the vector field X = R X is such that when [Y, X] = h X, we have [Y, X] = h X, with h = (Y(R)/R) + h, and hence, the preceding result shows that the function I = div (Y) + h = div (Y) + Y(log R) + h is a constant of motion for X and, therefore, for X [43].

A Geometric Approach to a Generalised Sundman Transformation
An infinitesimal time-reparametrisation, usually called a Sundman transformation [44], introduced when looking for an analytic solution to the three-body problem [45][46][47], is very intriguing, at least from a geometric perspective, because its geometric meaning is not clear.So, in the nonautonomous approach, t is not a coordinate but rather a parameter of the integral curve.It was introduced for a second-order differential equation, but it admits a clear geometric generalisation to systems of first-order differential equations.
The geometric approach was developed in [28,48].This infinitesimal timereparametrisation is written as It has allowed researchers to find solutions to many different problems in the theory of differential equations and related physical problems, for instance, the linearisation of differential equations, or to regularise some equations of motion and to avoid collision singularities.More generally, it is used to transform a given equation into another one of a more appropriate form.Moreover, as a SODE vector field on Q can be seen as a particular case of a vector field on TQ, it is enough to consider, for the time being, Sundman transformations for systems of first-order differential equations.
The meaning of infinitesimal time-reparametrisation ( 22) is based on the following property: if γ(t) is a curve, we carry out the mentioned Sundman transformation by introducing a new parameter τ into such a curve such that and we obtain the reparametrised curve γ(τ), given by γ(τ(t)) = γ(t), and then Consequently, if γ(t) is an integral curve of X, i.e., dγ/dt = X γ(t) , then i.e., the curve γ is an integral curve of the vector field f X.This implies that we can consider the given Sundman transformation to be a change from the vector field X to the conformally related one f X [28,48], but note that the integral curves of the vector fields X ∈ X(M) and f X ∈ X(M) coincide up to their respective reparametrisations.Of course, such infinitesimal time-reparametrisation implies a new velocity v with respect to the new time, and then it is related to the old one by vi = f v i .
It was proven in [49] that for a given incomplete vector field X, it is possible to find a function f ∈ C ∞ (M) such that the vector field f X is a complete vector field, and this property shows the usefulness of Sundman transformations in the regularisation process.The relevance of the property of replacing a vector field with a conformally related one in the Kustaanheimo-Stiefel regularisation of the Kepler problem was also demonstrated in [49].Conformally related vector fields have also been considered in the study of superintegrable systems (see, e.g., [22,23]).Moreover, two conformally related vector fields X and f X on a manifold M have the same local first integrals, but, as indicated in [49], this is not true for general tensor fields T, because L f X T ̸ = f L X T. This is an important fact because it provides us with a way to obtain, from a vector field X that does not preserve T, vector fields f X under which a given tensor field T is invariant.
As an example, we can focus on the linearisation of a system of first-order differential equations.

• Linearisation
Recall that a manifold M can be endowed with a linear structure if there exists a complete vector field ∆, playing the role of a Liouville vector field, with only one nondegenerate critical point, such that F (0) and F (1) ∆ separates derivations [5].Here, F (k) ∆ denotes the set of functions on the manifold M defined by Linear functions (with respect to ∆) are those defined by ∆ , and linear vector fields X ∈ X(M) are those preserving the R-linear subspace F (1) This happens when there exists h such that [X, ∆] = h X and the function f is such that h = ∆(log f ).
A chart (U, φ) of M identifies U with an open set of R n , and then U is endowed with the complete vector field ∆ = x i ∂/∂x i , and if X = X i (x) ∂/∂x i , we say that the vector field X is linear when [∆, X] = 0, and as [∆, X] = 0 is equivalent to X i (x) = A i j x j , with integral curves for X solutions of ẋi = A i j x j .Similarly, vector fields X whose integral curves are solutions of an inhomogeneous linear system, ẋi = A i j x j + B i , are characterised by The preceding results can be summarised as follows [48]: Theorem 3. Let X be a vector field on a differentiable manifold M endowed with a linear structure ∆.

1.
L ∆ X = h X, for some function h, if and only if there exists a nonvanishing positive function f such that the vector field f X satisfies L ∆ ( f X) = 0.

2.
L ∆ X = 0 if and only if, for any coordinate system (U, φ) such that ∆| U = x i ∂ ∂x i , the local expression of X is linear.

3.
L 2 ∆ X = −L ∆ X if and only if for any coordinate system such that ∆| U = x i ∂ ∂x i , the local expression of X is affine in x i .

Divergence-free case
If a vector field X ∈ X(M) on an oriented manifold (M, Ω) is such that L X Ω ̸ = 0, there may be a positive function f ∈ C ∞ (M) such that L f X Ω = 0.These functions are called Jacobi multipliers and play a relevant role in its integrability by quadratures.
Recall that as L X Ω is also a volume form, we can define the divergence of X by L X Ω = div X Ω, and then X is divergence-free if and only if L X Ω = 0. Note that as L f X Ω = L X ( f Ω), the search for such a Jacobi multiplier is equivalent to the determination of an X-invariant volume form f Ω (more details can be found in [30,31]).

Geodesic case
Another possibility would be to consider a Riemannian structure on the manifold M given by a nondegenerate symmetric two-times covariant tensor field g.As L f X g ̸ = f L X g, some non-Killing vector fields X, i.e., such that L X g ̸ = 0, can give rise to Killing ones [7] just by multiplication by a convenient function f .

Hamiltonian case
Given a closed 2-form ω that is non-invariant under the vector field X on a manifold, i.e., L X ω ̸ = 0, if there is a positive function f such that L f X ω = 0, i.e., if f i(X)ω is a closed form, then when the nondegeneracy condition d f ∧ ω = 0 is satisfied, the 2-form f ω is symplectic, and the vector field X is locally Hamiltonian with respect to this symplectic form f ω.Recall that, in the case where X is also locally Hamiltonian with respect to the original symplectic form ω, such a function f must be a constant of motion for X, X f = 0.Moreover, we can make use of the alternative symplectic structure to find non-Noether constants of motion, as explained in Section 4.
The same can be said about bivector fields Λ, particularly those of Poisson structures, i.e., such that [Λ, Λ] S = 0, where [•, •] S denotes the Schouten bracket.In this case, defines a Poisson bracket.Hamiltonian vector fields X f on M are those of the form X f g = {g, f }.Given a vector field X and a Poisson structure Λ, sometimes, there exists a function f such that f X is Hamiltonian with respect to the Poisson structure Λ: a Hamiltonisation process.

Sundman Transformations for Systems of Second-Order Differential Equations
As indicated above, the Sundman transformation was first introduced to deal with systems of second-order differential equations, such as Newton's equations of motion.It is well known that an autonomous system is associated with a system of 2n first-order differential equations whose solutions define the integral curves of the SODE vector field Γ on TQ i.e., a system of second-order differential equations can be dealt with using a SODE vector field on the tangent bundle, Γ ∈ X(TQ), such that S(Γ) = ∆.
The theory of Sundman transformations developed for systems of first-order differential equations suggests that we proceed by similarity and define the transformed vector field as obtained by multiplication by the function defining the Sundman transformation.But the appearance of second-order derivatives suggests that this approach may need to be modified.On the other hand, when multiplying by the function f , the SODE characteristic property is lost.This leads us to examine this definition more carefully.It will be shown that we can overcome these two problems in the particularly important case of the function f defining the Sundman transformation, being a constant-sign basic function.
As remarked in [48], when applying a Sundman transformation (22), the new velocities with respect to the new time, v, are related to the previous ones by vi = f v i , because and therefore, This suggests the use of non-natural coordinates on TQ, the so-called quasi-velocities [50].With the change of coordinates to a new coordinate system on the manifold TQ (we keep the notation xi for clarity), we see, by direct computation, that Consequently, the new coordinate expression of the preceding vector field Γ is Under a generalised Sundman transformation of the given system, we can use, as suggested in [48], the operational relations and hence, if x(t) is a solution of the system, and the new parameter τ is given by then x(τ) such that x(τ(t)) = x(t) satisfies the SODE system together with vi = f v i , i.e., d xi /dτ = f dx i /dt.The system can be rewritten as This shows that the images under the generalised Sundman transformation of the integral curves of the considered vector field Γ are the integral curves of the vector field Γ( x, v) and a simple comparison shows that

Alternative Tangent Bundle Structures
The relevance and usefulness of the existence of alternative geometric structures in the description of mechanical systems, particularly the tensorial characterisation of linear structures and vector bundle structures, are also well known [5].
We first review the work by Crampin and Thompson on almost-tangent structures on M [26].An almost-tangent structure on a manifold M is a (1,1)-tensor field S on M such that, at each point p ∈ M, the kernel of the linear map S p : T p M / / T p M coincides with its image.It follows that S 2 = 0.The almost-tangent structure S is said to be integrable if its Nijenhuis tensor N S vanishes.We recall that the Nijenhuis tensor N T of a (1,1)-tensor field T on a manifold M is a (1,2)-tensor field given by The N S = 0 condition implies that the vertical distribution defined by ker S = Im S is involutive, because if X, Y ∈ ker S = Im S, then there exist X, Y ∈ X(M) such that S( X) = X and S( Y) = Y, and then i.e., [X, Y] lies in the image of S, which coincides with ker S, and therefore, the distribution defined by ker S is involutive and hence integrable in the Darboux sense.The integral submanifolds of such a foliation are the fibres of a fibration π : M / / Q, where Q is an n-dimensional manifold.The vertical vector fields are those taking values in ker S.Then, if X ∈ X(M) projects onto a vector field X ∈ X(Q), then S(X) only depends on X.
Crampin and Thompson showed [26,51] that there are local coordinates (x i , u i ) on M such that S = (∂/∂u i ) ⊗ dx i .Here, (x 1 , . . ., x n ) are coordinates on Q, and if X 1 , . . ., X n are vector fields on M projecting onto ∂/∂x 1 , . . ., ∂/∂x n , respectively, as [S(X i ), S(X j )] = 0, then local coordinates on the fibres can be chosen such that S(X i ) = ∂/∂u i .The coordinates u i are determined only up to an additive constant on each fibre, because they are solutions of the system S(X i )u j = δ j i , and we can change u i by ūi = u i + f i (x).The ambiguity is fixed by the choice of the zero section.
Such a choice uniquely determines, for each coordinate system (x i ) on the base manifold Q, a system of local coordinates (x i , u i ) such that Furthermore, there is a uniquely defined vertical vector field ∆ such that L ∆ S(X) = −S(X) and ∆ = 0 on the zero section.The explicit expression of ∆ is Such a vector field ∆ provides a linear space structure to every fibre of π : M / / Q.
The vector fields associated with systems of second-order differential equations with respect to the new integrable almost-tangent structure are such that S(D) = ∆, i.e., of the form ∂ ∂u i .It was shown in [48] that, starting with the standard tangent bundle structure given by the vertical endomorphism S and a constant-sign basic function τ * Q h on TQ, with h ∈ C ∞ (Q), we can introduce a new integrable almost-tangent structure (τ * Q h) S on TQ, simply denoted by S = h S, because ker(h S) = Im(h S), and N h S = 0.The vertical distribution defined by ker(h S) coincides with the usual one defined by ker S, and the set of leaves can be identified with the base manifold Q.
If we choose a local set of coordinates for Q, (q 1 , . . ., q n ), and if X 1 , . . ., X n are τ Qprojectable vector fields on TQ and τ Q -related to ∂/∂q 1 , . . ., ∂/∂q n , respectively, then the vector fields h S(X 1 ), . . ., h S(X n ) are pairwise commuting, and there exist local coordinates v1 , . . ., vn such that S(h X i ) = h S(X i ) = ∂/∂ vi , and as S(X i ) = ∂/∂v i , with the same choice for the zero section, we can see that S = (τ * Q h) S has the local expression with the fibre coordinates given by h vi = v i .The vector field ∆ coincides with ∆.
Note that vi values are not coordinates adapted to the original tangent bundle structure because the 1-forms α i defining such coordinates of the vector v by vi = α i (v) are not exact, with α i = (1/h)dq i , and correspond to the so-called quasi-velocities [50].In summary, we obtain a new tangent bundle structure on TQ, with its corresponding Liouville vector field ∆ = ∆ and a new concept of a SODE vector field.When h = f −1 , we see that Γ = f Γ is a SODE vector field with respect to S = (1/ f )S.
If we successively apply two Sundman transformations to a system characterised by the functions f 1 ∈ C ∞ (Q) and f 2 ∈ C ∞ (Q), we obtain a new system by applying the Sundman transformation defined by the product f 2 f 1 ∈ C ∞ (Q), and we see that the transformed vector field is a SODE vector field with respect to the new tangent bundle structure S = f −1 2 f −1 1 S.As the set of positive real functions is an Abelian multiplicative Lie group, the Sundman transformation corresponding to f −1 is the inverse of the Sundman transformation defined by f .More details can be found in [48].

Linearisability of Scalar Second-Order Differential Equations
A particularly interesting case is that of systems of second-order differential equations that are linear in a given coordinate system: i.e., the functions X i are of the form where A i j and B i j are real numbers.In geometric terms, this can be written as [ ∆, Γ] = 0, with ∆ = ∆ Q + ∆, where Note that the vector field ∆ is not intrinsic, but it depends on the choice of the chart on U.More explicitly, the vector field ∆ is the complete lift of ∆ Q .Our aim is to study under which circumstances a given system of second-order differential equations can be transformed by an appropriate Sundman transformation into a linear one.
Consider the simplest case of a one-dimensional autonomous SODE, i.e., n = 1, and X 1 = F is a basic function.The equation of the integral curves of the vector field Γ transformed under a Sundman transformation of Γ = v ∂/∂q + F(q) ∂/∂v is and, in the particular case of a system defined by a potential function V, where F = −V ′ , the expressions of the SODE vector field Γ in coordinates (q, v) and its transformed vector field Γ = f Γ in coordinates ( q, v) are given, respectively, by But as the energy function V is a conserved quantity, if we restrict ourselves to studying the motions for a given energy E, This is an inhomogeneous linear differential equation in the variable q if and only if there exist constants A, B and C such that f 2 (E − V ) = A q 2 + B q + C, from where In the relevant case of the radial equation for a given fixed angular momentum ℓ for the Coulomb-Kepler problem, for which q is the radial variable r and V(r and from f 2 (E − V ) = Ar 2 + B + C, we obtain , and therefore, the general solution is a multiple of f (r) = r, with and then the transformed second-order differential equation is [52] The Sundman transformation is independent of ℓ, because f (r) is proportional to r.
In the opposite direction, if we start from the linear SODE vector field Γ of the onedimensional harmonic oscillator, i.e., Γ = v (∂/∂x) − ω 2 x (∂/∂v), which corresponds to F(x) = −ω 2 x, i.e., V = 1 2 ω 2 x 2 , under the Sundman transformation defined by a given function f , it becomes and the projections onto the base manifold of its integral curves are solutions of i.e., it appears that an additional quadratic damping term and linearity are lost.
Conversely, an equation of this last type can be reduced to a harmonic oscillator by the Sundman transformation defined by f −1 .If, for instance, f (x) = x 2 , the former equation reduces to i.e., it is the Ermakov-Pinney equation under the action of a damping quadratic term.This is a prototype for linearisable examples.As indicated in [48], the linear structure on Q = R depends on the choice of the coordinate y, ∆ Q = y ∂/∂y, and then the linearisability character of a given SODE depends on the choice of y.Therefore, in the study of the linearisation process, we need to first find a coordinate transformation y = φ(x) from the original coordinate to the adapted one and then a Sundman transformation, dτ = h(x) dt, that transforms our original SODE into a linear one in the new coordinate y and its velocity.This will be implemented in several steps by finding several coordinate transformations and Sundman transformations that simplify the form of the SODE.Then, we describe the composition of ordinary changes of coordinates with the Sundman transformation considered so far.We will refer to a transformation of the form y = x, dτ = h(x) dt as a pure Sundman transformation and to y = φ(x), dτ = dt as a pure coordinate transformation.By composing such types of transformations, we obtain a group of generalised Sundman transformations (h, φ) defined as with the composition law This set of generalised Sundman transformations is a semidirect product group whose neutral element is (1, id), and the inverse of (h, φ) is ((1/h) • φ −1 , φ −1 ).Pure Sundman transformations are those of the form (h, id) and constitute an Abelian invariant subgroup, while the usual coordinate transformations correspond to those of the form (1, φ) and also make up a subgroup.Finally, each transformation can be factorised as Definition 2. A SODE is fibre-linearisable (or linearisable in velocities) under a generalised Sundman transformation if it can be transformed to a SODE of the form ẍ + A(x) ẋ + b(x) = 0, where A and b are real functions, and a SODE is linearisable under generalised a Sundman transformation if it can be transformed into a SODE of the form ẍ + α ẋ + Bx + C = 0 for some real numbers α, B, C ∈ R. Note that when C ̸ = 0, the transformed equation is an inhomogeneous linear equation.
The group properties of the set of generalised Sundman transformations can be used to study the linearisability of a scalar SODE: The possible linearising transformations will be factorisable as a composition of a pure coordinate transformation and a pure Sundman transformation.Inverting the process, we can first see the form of the image under a pure Sundman transformation of the prototype linear equation.We can see that its image under the general transformation rule for the one-dimensional case is an equation of the form i.e., the function X(x, v) is a polynomial of degree two at most in the variable v = dx/dt, ∂ 3 X/∂v 3 = 0.Such a class is invariant under changes of coordinates because if x = φ(x), and therefore, x satisfies an equation of the same type.Therefore, we obtain the following: Theorem 4. A necessary condition for fibre-linearisability under a generalised Sundman transformation is that the function X(x, v) be a polynomial of degree 2 at most in the variable v = dx/dt, or in other words, ∂ 3 X/∂v 3 = 0, and consequently, second-order differential equations linearisable under a generalised Sundman transformation must be of the previous form.
Such a class of equations is invariant under a generalised Sundman transformation and contains a subset of linear equations that correspond to the null function γ 0 = 0. We must determine which of these equations are linearisable.
Note that the quadratic in the velocity term can always be eliminated by an appropriate pure Sundman transformation, x 1 = x, dt 1 = h(x) dt.Under such a transformation, the SODE becomes We can choose the function h as a nontrivial solution h 0 (x) of the linear ODE dh dx + γ 0 (x)h = 0, with the general solution h(x) = K exp − γ 0 (x) dx , with K being a constant, which is irrelevant for our purposes.With this choice, under the pure Sundman transformation defined by h, the second-order differential equation becomes the fibre-linear equation with A 1 = A 0 /h 0 and b 1 = b 0 /h 2 0 .Therefore, any SODE of the given class is fibrelinearisable by a pure Sundman transformation.If both A 1 and b 1 vanish identically (i.e., A 0 and b 0 vanish identically), then we have already obtained a special linear equation, more specifically, d 2 x 1 /dt 2 1 = 0. Note that, in the case where at least one of such functions is not the zero function, once we have obtained an inhomogeneous fibre-linear SODE, we can use only the Sundman transformation preserving the set of inhomogeneous fibre-linear SODEs.These are more general than affine transformations of coordinates, as stated in the following lemma.
The value of the constant c does not have any influence on the linearisability, and without losing generality, we can take c = 1 (or c = −1 when appropriate).
This shows that we can only transform the fibre-linear SODE by means of a Sundman transformation x 2 = φ(x 1 ), and dt 2 = h(x 1 ) dt 1 , with φ(x) = x h(ζ) dζ.As for such a transformation, dx 2 /dt 2 = dx 1 /dt 1 , the expression of the transformed SODE is where the substitution of x 1 by φ −1 (x 2 ) in the last two terms is understood.This SODE is linear if and only if the function A 1 /h is constant and the function b 1 /h is affine in the variable x 2 .There are two different situations: (a) If A 1 vanishes identically (also recall that A 0 vanishes identically), we can choose the function h as h = |b 1 | in the Sundman transformation, and the transformed SODE is The Sundman transformation is x 2 = where we assume that x 1 is replaced by its corresponding value in terms of x 2 .
A generalised Sundman transformation preserving the form of the above SODE, i.e., inhomogeneous fibre-linear with a constant coefficient α, must be of the form x = m x + n, d t = m dt, with m and n as constants, and therefore, it is either linear or, otherwise, not linearisable.Obviously, the new SODE is linear if and only if db 2 /dx 2 is constant, db 2 /dx 2 = B ∈ R.
One can find the conditions in the original data γ 0 , A 0 and b 0 for b 2 (x 2 ) to be an affine function.The following linearisability condition can be proven [48]: the SODE is linearisable if and only if the functions A 0 , b 0 and γ 0 satisfy Notice that this condition is also satisfied in the first two cases (A 0 = b 0 = 0 and A 0 = / 0, b 0 = 0), and moreover, we can prove that such a condition is invariant under a generalised Sundman transformation.
All of these results can be summarised by the following statement [48]: (ii) A second-order differential equation is linearisable if and only if it is of the preceding form and the coefficients γ(x), A(x) and b(x) satisfy More specifically, we find the following: To illustrate the theory and for the sake of completeness, we mention several examples, which have been more developed in [48].
2. Geodesics on the sphere: The differential equation is ẍ = 2 ẋ2 cot x + sin x cos x, so that γ(x) = −2 cot x, A(x) = 0 and b(x) = − sin x cos x.It is Sundman-linearisable into its canonical form by means of the transformation given by y = Notice that this equation is not linearisable by point transformations, as it does not satisfy Lie criteria, but it is Sundman-linearisable.
4. As a generalisation of the preceding example, we consider the differential equation which, in our particular case, turns out to be and therefore, Then, given the function f , the function g is any function of the linear space of solutions of the linear second-order differential equation in the variable g.As g = f is a solution of such an equation, we introduce the change of the variable g = f ζ, and the new equation is f ζ ′′ − f ′ ζ ′ = 0, which shows that the general solution is and therefore, the linearisability condition implies that the function g is Finally, the Liénard-type equation containing a dissipative term studied in [57] ẍ + f (x) ẋ2 + g(x) ẋ + h(x) = 0 (53) can be reduced by a pure Sundman transformation dτ = F(x) dt to a Liénard Equation ( 51), as more explicitly shown in [48].

Lemma 1 .
If a generic Sundman transformation x = φ(x), d t = h(x) dt preserves the set of inhomogeneous fibre-linear SODEs, then there exists a real constant c = / 0 such that dφ/dx = c h.Moreover, a Sundman transformation x = φ(x), d t = h(x) dt, such that the condition holds, transforms every inhomogeneous fibre-linear SODE d 2 x dt 2 + A(x) dx dt + b(x) = 0 (41) into an analogous one with the coefficients Ā( x) and b( x), where