A Reduction Method for Higher Order Variational Equations of Hamiltonian Systems

Let $\mathbf{k}$ be a differential field and let $[A]\,:\,Y'=A\,Y$ be a linear differential system where $A\in\mathrm{Mat}(n\,,\,\mathbf{k})$. We say that $A$ is in a reduced form if $A\in\mathfrak{g}(\bar{\mathbf{k}})$ where $\mathfrak{g}$ is the Lie algebra of $[A]$ and $\bar{\mathbf{k}}$ denotes the algebraic closure of $\mathbf{k}$. We owe the existence of such reduced forms to a result due to Kolchin and Kovacic \cite{Ko71a}. This paper is devoted to the study of reduced forms, of (higher order) variational equations along a particular solution of a complex analytical hamiltonian system $X$. Using a previous result \cite{ApWea}, we will assume that the first order variational equation has an abelian Lie algebra so that, at first order, there are no Galoisian obstructions to Liouville integrability. We give a strategy to (partially) reduce the variational equations at order $m+1$ if the variational equations at order $m$ are already in a reduced form and their Lie algebra is abelian. Our procedure stops when we meet obstructions to the meromorphic integrability of $X$. We make strong use both of the lower block triangular structure of the variational equations and of the notion of associated Lie algebra of a linear differential system (based on the works of Wei and Norman in \cite{WeNo63a}). Obstructions to integrability appear when at some step we obtain a non-trivial commutator between a diagonal element and a nilpotent (subdiagonal) element of the associated Lie algebra. We use our method coupled with a reasoning on polylogarithms to give a new and systematic proof of the non-integrability of the H\'enon-Heiles system. We conjecture that our method is not only a partial reduction procedure but a complete reduction algorithm. In the context of complex Hamiltonian systems, this would mean that our method would be an effective version of the Morales-Ramis-Sim\'o theorem.Comment: 15 page

A Reduced Form for Linear Differential Systems and its Application to Integrability of Hamiltonian Systems

Let $[A]: Y'=AY$ with $A\in \mathrm{M}_n (k)$ be a differential linear system. We say that a matrix $R\in {\cal M}_{n}(\bar{k})$ is a {\em reduced form} of $[A]$ if $R\in \mathfrak{g}(\bar{k})$ and there exists $P\in GL_n (\bar{k})$ such that $R=P^{-1}(AP-P')\in \mathfrak{g}(\bar{k})$. Such a form is often the sparsest possible attainable through gauge transformations without introducing new transcendants. In this article, we discuss how to compute reduced forms of some symplectic differential systems, arising as variational equations of hamiltonian systems. We use this to give an effective form of the Morales-Ramis theorem on (non)-integrability of Hamiltonian systems.Comment: 28 page

A Characterization of Reduced Forms of Linear Differential Systems

A differential system $[A] : \; Y'=AY$, with $A\in \mathrm{Mat}(n, \bar{k})$ is said to be in reduced form if $A\in \mathfrak{g}(\bar{k})$ where $\mathfrak{g}$ is the Lie algebra of the differential Galois group $G$ of $[A]$. In this article, we give a constructive criterion for a system to be in reduced form. When $G$ is reductive and unimodular, the system $[A]$ is in reduced form if and only if all of its invariants (rational solutions of appropriate symmetric powers) have constant coefficients (instead of rational functions). When $G$ is non-reductive, we give a similar characterization via the semi-invariants of $G$. In the reductive case, we propose a decision procedure for putting the system into reduced form which, in turn, gives a constructive proof of the classical Kolchin-Kovacic reduction theorem.Comment: To appear in : Journal of Pure and Applied Algebr

Non-integrability of the generalised spring-pendulum problem

We investigate a generalisation of the three dimensional spring-pendulum system. The problem depends on two real parameters $(k,a)$, where $k$ is the Young modulus of the spring and $a$ describes the nonlinearity of elastic forces. We show that this system is not integrable when $k\neq -a$. We carefully investigated the case $k= -a$ when the necessary condition for integrability given by the Morales-Ramis theory is satisfied. We discuss an application of the higher order variational equations for proving the non-integrability in this case.Comment: 20 pages, 1 figur

Effective descent for differential operators

A theorem of N. Katz \cite{Ka} p.45, states that an irreducible differential operator $L$ over a suitable differential field $k$, which has an isotypical decomposition over the algebraic closure of $k$, is a tensor product $L=M\otimes_k N$ of an absolutely irreducible operator $M$ over $k$ and an irreducible operator $N$ over $k$ having a finite differential Galois group. Using the existence of the tensor decomposition $L=M\otimes N$, an algorithm is given in \cite{C-W}, which computes an absolutely irreducible factor $F$ of $L$ over a finite extension of $k$. Here, an algorithmic approach to finding $M$ and $N$ is given, based on the knowledge of $F$. This involves a subtle descent problem for differential operators which can be solved for explicit differential fields $k$ which are $C_1$-fields.Comment: 21 page

Diagonals of rational functions and selected differential Galois groups

International audienceWe recall that diagonals of rational functions naturally occur in lattice statistical mechanics and enumerative combinatorics. In all the examples emerging from physics, the minimal linear differential operators annihilating these diagonals of rational functions have been shown to actually possess orthogonal or symplectic differential Galois groups. In order to understand the emergence of such orthogonal or symplectic groups, we analyze exhaustively three sets of diagonals of rational functions, corresponding respectively to rational functions of three variables, four variables and six variables. We impose the constraints that the degree of the denominators in each variable is at most one, and the coefficients of the monomials are 0 or $\;\pm 1,$ so that the analysis can be exhaustive. We find the minimal linear differential operators annihilating the diagonals of these rational functions of three, four, five and six variables. We find that, even for these sets of examples which, at first sight, have no relation with physics, their differential Galois groups are always orthogonal or symplectic groups. We discuss the conditions on the rational functions such that the operators annihilating their diagonals do not correspond to orthogonal or symplectic differential Galois groups, but rather to generic special linear groups

Efficient Algorithms for Computing Rational First Integrals and Darboux Polynomials of Planar Polynomial Vector Fields

International audienceWe present fast algorithms for computing rational first integrals with bounded degree of a planar polynomial vector field. Our approach builds upon a method proposed by Ferragut and Giacomini, whose main ingredients are the calculation of a power series solution of a first order differential equation and the reconstruction of a bivariate polynomial annihilating this power series. We provide explicit bounds on the number of terms needed in the power series. This enables us to transform their method into a certified algorithm computing rational first integrals via systems of linear equations. We then significantly improve upon this first algorithm by building a probabilistic algorithm with arithmetic complexity $\~O(N^{2 \omega})$ and a deterministic algorithm solving the problem in at most $\~O(d^2N^{2 \omega+1})$ arithmetic operations, where~$N$ denotes the given bound for the degree of the rational first integral, and where $d \leq N$ is the degree of the vector field, and $\omega$ the exponent of linear algebra. We also provide a fast heuristic variant which computes a rational first integral, or fails, in $\~O(N^{\omega+2})$ arithmetic operations. By comparison, the best previous algorithm uses at least $d^{\omega+1}\, N^{4\omega +4}$ arithmetic operations. We then show how to apply a similar method to the computation of Darboux polynomials. The algorithms are implemented in a Maple package RationalFirstIntegrals which is available to interested readers with examples showing its efficiency

Differential algebra on lattice green and calabi-yau operators

International audienceWe revisit miscellaneous linear differential operators mostly associated with lattice Green functions in arbitrary dimensions, but also Calabi-Yau operators and order-seven operators corresponding to exceptional differential Galois groups. We show that these irreducible operators are not only globally nilpotent, but are such that they are homomorphic to their (formal) adjoints. Considering these operators, or, sometimes, equivalent operators, we show that they are also such that, either their symmetric square or their exterior square, have a rational solution. This is a general result: an irreducible linear differential operator homomorphic to its (formal) adjoint is necessarily such that either its symmetric square, or its exterior square has a rational solution, and this situation corresponds to the occurrence of a special differential Galois group. We thus define the notion of being "Special Geometry" for a linear differential operator if it is irreducible, globally nilpotent, and such that it is homomorphic to its (formal) adjoint. Since many Derived From Geometry n-fold integrals ("Periods") occurring in physics, are seen to be diagonals of rational functions, we address several examples of (minimal order) operators annihilating diagonals of rational functions, and remark that they also seem to be, systematically, associated with irreducible factors homomorphic to their adjoint