Next Article in Journal
Burr XII Distribution for Disease Data Analysis in the Presence of a Partially Observed Failure Mode
Next Article in Special Issue
An Asymptotic Expansion for the Generalized Gamma Function
Previous Article in Journal
Stable Exponential Cosmological Type Solutions with Three Factor Spaces in EGB Model with a Λ-Term
Previous Article in Special Issue
Sequence Spaces and Spectrum of q-Difference Operator of Second Order
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diagonals of Rational Functions: From Differential Algebra to Effective Algebraic Geometry

by
Youssef Abdelaziz
1,
Salah Boukraa
2,
Christoph Koutschan
3 and
Jean-Marie Maillard
1,*
1
LPTMC, UMR 7600 CNRS, Sorbonne Université, Tour 23, 5ème étage, Case 121, 4 Place Jussieu, CEDEX 05, 75252 Paris, France
2
LSA, IAESB, Université de Blida 1, Blida 09000, Algeria
3
Johann Radon Institute for Computational and Applied Mathematics, RICAM, Altenberger Strasse 69, A-4040 Linz, Austria
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1297; https://doi.org/10.3390/sym14071297
Submission received: 2 May 2022 / Revised: 7 June 2022 / Accepted: 19 June 2022 / Published: 22 June 2022

Abstract

:
We show that the results we had previously obtained on diagonals of 9- and 10-parameter families of rational functions in three variables x , y , and z , using creative telescoping, yielding modular forms expressed as pullbacked 2 F 1 hypergeometric functions, can be obtained much more efficiently by calculating the j -invariant of an elliptic curve canonically associated with the denominator of the rational functions. These results can be drastically generalized by changing the parameters into arbitrary rational functions of the product p = x y z . In other cases where creative telescoping yields pullbacked 2 F 1 hypergeometric functions, we extend this algebraic geometry approach to other families of rational functions in three or more variables. In particular, we generalize this approach to rational functions in more than three variables when the denominator can be associated to an algebraic variety corresponding to products of elliptic curves, or foliations in elliptic curves. We also extend this approach to rational functions in three variables when the denominator is associated with a genus-two curve such that its Jacobian is a split Jacobian, corresponding to the product of two elliptic curves. We sketch the situation where the denominator of the rational function is associated with algebraic varieties that are not of the general type, having an infinite set of birational automorphisms. We finally provide some examples of rational functions in more than three variables, where the telescopers have pullbacked 2 F 1 hypergeometric solutions, because the denominator corresponds to an algebraic variety that has a selected elliptic curve.

1. Introduction

In previous papers [1,2], using creative telescoping [3], we have obtained diagonals (for the introduction of the concept of diagonals of rational functions, see [4,5,6,7,8,9,10,11]) of 9- and 10-parameter families of rational functions, given by (classical) modular forms expressed as pullbacked 2 F 1 hypergeometric functions [12]. The natural emergence of diagonals of rational functions (the lattice Green functions are the simplest examples of such diagonals of rational functions [13,14,15,16,17,18]) in lattice statistical mechanics is explained in [19,20]. This can be seen as the reason for the frequent occurrence of modular forms and Calabi–Yau operators in lattice statistical mechanics [21,22,23,24,25,26,27]. In other previous papers [17,18] dedicated to Heun functions that are either diagonals of simple rational functions, or only solutions of telescopers [28,29] of simple rational functions of three or four variables, we have obtained many order-three telescopers having squares of Heun functions as solutions that turn out to be squares of pullbacked 2 F 1 hypergeometric functions corresponding to classical modular forms and even Shimura automorphic forms [30,31], strongly reminiscent of periods of extremal rational surfaces [32,33], and other foliations of K3 surfaces in elliptic curves. In other words, one finds experimentally that the 2 F 1 hypergeometric functions emerging in the calculation of the diagonals of rational functions, or of the solutions of the telescopers of rational functions, seem to be only special 2 F 1 ( a , b ; c ; x ) hypergeometric functions with a selected set of parameters a , b , c (see the list (B.1) in Appendix B of [17], corresponding to classical modular forms) (see Felix Klein’s connection of the 2 F 1 ( 1 / 12 , 5 / 12 ; 1 ; x ) Gauss hypergeometric function with modular forms; for instance, in the very pedagogical and heuristic paper, [12]), together with a finite set of parameters, such as 7 / 24 , 11 / 24 , 5 / 4 , corresponding to Shimura automorphic forms [30,31]), pullbacked by selected pullbacks. This last paper [17] also underlined the difference between the diagonal of a rational function Diag ( R ) , and the solutions of the telescoper of the same rational function.
These results strongly suggested to find an algebraic geometry interpretation for all of these exact results, and more generally, they suggested to provide an alternative algebraic geometry approach of the results emerging from creative telescoping (the reader may refer to [34] for an extensive survey of “creative telescoping” approaches).
This is the purpose of the present paper. In particular, we are going to show that most of these pullbacked 2 F 1 hypergeometric functions can be obtained efficiently through algebraic geometry calculations, thus providing a more intrinsic algebraic geometry interpretation of the creative telescoping calculations that are typically differential algebra calculations [28,29,34,35].
Creative telescoping [28,29,34,36] is a methodology to deal with parametrized symbolic sums and integrals that yields differential/recurrence equations for such expressions. This methodology became popular in computer algebra in the past 25 years. By the “telescoper” of a rational function, say R ( x , y , z ) , we here refer to the output of the creative telescoping program [3], applied to the transformed rational function R ˜ = R ( x / y , y / z , z ) / ( y z ) . Such a telescoper is a linear differential operator T in x and x , such that T + y · U + z · V annihilates R ˜ , where the so-called “certificates” U , V are rational functions in x , y , z . In other words, the telescoper T represents a linear ODE that is satisfied by Diag ( R ) .
The paper is essentially dedicated to the solutions of telescopers of rational functions that are not necessarily diagonals of rational functions. These solutions correspond to periods [37] of algebraic varieties over some cycles that are not necessarily vanishing [38] cycles (in French: “cycles évanescents” [5,39]) as with the case of the diagonals of rational functions. The reader who is interested in the connection between the process of taking diagonals, calculating telescopers, and the notion of periods, de Rham cohomology (i.e., differential forms) and other Picard–Fuchs equations can read in detail the thesis of Pierre Lairez [35] (see also [40]). We sketch just some of these ideas in Appendix A.
The purpose of this paper is not to give an introduction to creative telescoping [28,29], but to provide many pedagogical (non-trivial) examples of telescopers using (one can obtain these telescopers using Chyzak’s algorithm [41] or Koutschan’s semi-algorithm [3,42] (the termination is not proven). For the examples displayed in this paper, Koutschan’s package [3] is more efficient) the HolonomicFunctions Mathematica package extensively [3].
The paper is organized as follows. We first recall in Section 2 the exact results of [1,2] for the 9- and 10-parameter families of rational functions using creative telescoping, yielding modular forms expressed as pullbacked 2 F 1 hypergeometric functions. We show in Section 3 that these exact results can be obtained, much more efficiently, by calculating the j -invariant of an elliptic curve that is canonically associated with the denominator of the rational function, and we underline the fact that one can drastically generalize these results, the parameters becoming quite arbitrary rational functions. Section 4 generalizes the previous calculations to denominators of the rational functions of more than three variables, corresponding to the products (or foliations) of elliptic curves. In Section 5, we show how modular forms expressed as pullbacked 2 F 1 hypergeometric functions occur for rational functions in three variables when the denominator is associated with a genus-two curve such that its Jacobian is a split Jacobian corresponding to the product of two elliptic curves. In Section 6, we sketch the situation where the denominator of the rational function is associated with algebraic varieties of low Kodaira dimension, having an infinite set of birational automorphisms. We finally provide some examples of rational functions in more than three variables, where the telescopers have pullbacked 2 F 1 hypergeometric solutions, the denominator corresponding to an algebraic variety having a selected elliptic curve in the variety explaining these pullbacked 2 F 1 solutions.

2. Classical Modular Forms and Diagonals of 9- and 10-Parameter Families of Rational Functions

In previous papers [1,2], using creative telescoping [3], we have obtained diagonals of 9- and 10-parameter families of rational functions, given by (classical) modular forms expressed as pullbacked 2 F 1 hypergeometric functions. Let us recall these results.

2.1. Nine-Parameter Rational Functions Giving Pullbacked 2 F 1 Hypergeometric Functions for Their Diagonals

Let us recall the nine-parameter rational function in three variables, x , y and z :
1 a + b 1 x + b 2 y + b 3 z + c 1 y z + c 2 x z + c 3 x y + d y 2 z + e z x 2 .
Calculating (using the HolonomicFunctions Mathematica package [3]) the telescoper of this rational function (1), one obtains an order-two linear differential operator annihilating the diagonal of the rational function (1). The diagonal of the rational function (1) can be written [1,2] as a pullbacked hypergeometric function:
1 P 4 ( x ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 1 P 6 ( x ) 2 P 4 ( x ) 3 ,
where P 4 ( x ) and P 6 ( x ) are two polynomials of degree four and six in x , respectively. The Hauptmodul pullback in (2) has the form:
H = 1728 j = 1 P 6 ( x ) 2 P 4 ( x ) 3 = 1728 · x 3 · P 8 ( x ) P 4 ( x ) 3 ,
where j is the j-invariant and P 8 ( x ) is a polynomial of degree eight in x . Such a pullbacked 2 F 1 hypergeometric function (2) corresponds to a classical modular form [1,2].

2.2. Ten-Parameter Rational Functions Giving Pullbacked 2 F 1 Hypergeometric Functions for Their Diagonals

Let us recall the 10-parameter rational function in three variables, x , y and z :
R ( x , y , z ) = 1 a + b 1 x + b 2 y + b 3 z + c 1 y z + c 2 x z + c 3 x y + d 1 x 2 y + d 2 y 2 z + d 3 z 2 x .
Calculating the telescoper of this rational function (4), one obtains an order-two linear differential operator annihilating the diagonal of the rational function (4). The diagonal of the rational function (4) can be written [1,2] as a pullbacked hypergeometric function:
1 P 3 ( x ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 1 P 6 ( x ) 2 P 3 ( x ) 3 ,
where P 3 ( x ) and P 6 ( x ) are two polynomials of degree three and six in x , respectively. Furthermore, the Hauptmodul pullback in (5) is seen to be of the form:
H = 1728 j = 1 P 6 ( x ) 2 P 3 ( x ) 3 = 1728 · x 3 · P 9 ( x ) P 3 ( x ) 3 .
where P 9 ( x ) is a polynomial of degree nine in x . Again, (5) corresponds to a classical modular form [1,2].

3. Deducing Creative Telescoping Results from Effective Algebraic Geometry

Obtaining the previous pullbacked hypergeometric results (2) and (5) required [1,2] an accumulation of creative telescoping calculations, and a lot of “guessing”, using all of the symmetries of the diagonals of these rational functions (1) and (4). We are looking for a more efficient and intrinsic way of obtaining these exact results. These two pullbacked hypergeometric results (2) and (5), are essentially “encoded” by their Hauptmodul pullbacks (3) and (6), or, equivalently, their corresponding j -invariants. The interesting question, which will be addressed in this paper, is whether it is possible to canonically associate elliptic curves whose j -invariants correspond precisely to these Hauptmoduls H = 1728 j .

3.1. Revisiting the Pullbacked Hypergeometric Results in an Algebraic Geometry Perspective

One expects such an elliptic curve to correspond to the singular part of the rational function, namely, the denominator of the rational function. Let us recall that the diagonal of a rational function in (for example) three variables is obtained through its multi-Taylor expansion [19,20]:
R ( x , y , z ) = m n l a m , n , l · x m y n z l ,
by extracting the “diagonal” terms, i.e., the powers of the product p = x y z :
D i a g R ( x , y , z ) = m a m , m , m · x m .
Consequently, it is natural to consider the algebraic curve corresponding to the intersection of the surface defined by the vanishing condition D ( x , y , z ) = 0 of the denominator D ( x , y , z ) of these rational functions (1) and (4), with the hyperbola p = x y z (where p is seen, here, as a constant). This amounts, for instance, to eliminating the variable z , substituting z = p x y in D ( x , y , z ) = 0 .

3.1.1. Nine-Parameter Case

In the case of the rational functions (1), this corresponds to the (planar) algebraic curve:
a + b 1 x + b 2 y + b 3 p x y + c 1 y p x y + c 2 x p x y + c 3 x y + d y 2 p x y + e p x y x 2 = 0 ,
which can be rewritten as a (general, nine-parameter) biquadratic:
a x y + b 1 x 2 y + b 2 x y 2 + b 3 p + c 1 p y + c 2 p x + c 3 x 2 y 2 + d p y 2 + e p x 2 = 0 .
Using formal calculations (namely, using with(algcurves) in Maple, and in particular, the command j_invariant) one can easily calculate the genus of the planar algebraic curve (10), and find that it is actually an elliptic curve (of genus one). Furthermore, one can (almost instantaneously) find the exact expression of the j -invariant of this elliptic curve as a rational function of the nine parameters a , b 1 , b 2 , , e in (1). One actually finds that this j -invariant is precisely the j , such that the Hauptmodul H = 1728 j is the exact expression (3). In other words, the classical modular form result (2) could have been obtained, almost instantaneously, by calculating the j -invariant of an elliptic curve that is canonically associated with the denominator of the rational function (1). The algebraic planar curve (10) corresponds to the most general biquadratic of two variables, which depends on nine homogeneous parameters. Such a general biquadratic is known to be an elliptic curve for generic values of the nine parameters. (So many results in integrable models correspond to this most general biquadratic: the Bethe ansatz of the Baxter model [43,44], the elliptic curve foliating the 16-vertex model [44], so many QRT birational maps [45], etc.)
Thus, the nine-parameter exact result (2) can be seen as a simple consequence of the fact that the most general nine-parameter biquadratic is an elliptic curve.

3.1.2. Ten-Parameter Case

In the case of the rational function (4), substituting z = p x y in D ( x , y , z ) = 0 , one obtains the 10-parameter bicubic:
a x y 2 + b 1 x 2 y 2 + b 2 x y 3 + b 3 p y + c 1 p y 2 + c 2 p x y + c 3 x 2 y 3 + d 1 x 3 y 3 + d 2 y 3 + d 3 p 2 = 0 .
As before, we find that this planar algebraic curve is actually an elliptic curve (Generically, the most general planar bicubic is not a genus-one algebraic curve. It is a genus-four curve.) and that the exact expression of its j -invariant is precisely the j of the Hauptmodul H = 1728 / j in (6).
Thus, this 10-parameter result (5) can again be seen as a simple consequence of the fact that there exists a family of 10-parameter bicubics (see (11)) which are elliptic curves for generic values of the 10 parameters.
These preliminary calculations are a strong incentive to attempt to replace the differential algebra calculations of creative telescoping with more intrinsic algebraic geometry calculations, or at least, to perform effective algebraic geometry calculations to provide an algebraic geometry interpretation of the exact results obtained from creative telescoping.

3.2. Finding Creative Telescoping Results from j -Invariant Calculations

One might think that these results are a consequence of the simplicity of the denominators of the rational functions (1) or (4), being associated with biquadratics or selected bicubics. In fact, these results are very general. Let us, for instance, consider a nine-parameter family of planar algebraic curves that are not biquadratics or (selected) bicubics:
a 1 x 4 + a 2 x 3 + a 3 x 2 + a 4 x + a 5 + a 6 x 2 y + a 7 y 2 + a 8 y + a 9 x y = 0 .
One can easily calculate the genus of this planar curve and see that this genus is actually one for arbitrary values of the a n ’s. Thus the planar curve (12) is an elliptic curve for generic values of the nine parameters a 1 , , a 9 . It is straightforward to see that the algebraic surface S ( x , y , z ) = 0 , corresponding to:
z · ( a 1 x 4 + a 2 x 3 + a 3 x 2 + a 4 x + a 5 + a 6 x 2 y + a 7 y 2 + a 8 y ) + a 9 p = 0 ,
will automatically be such that its intersection with the hyperbola p = x y z gives back the elliptic curve (12).
Using this kind of “reverse engineering” suggests that we should consider the rational function in three variables x , y , and z :
R ( x , y , z ) = 1 1 + z · ( a 1 x 4 + a 2 x 3 + a 3 x 2 + a 4 x + a 5 + a 6 x 2 y + a 7 y 2 + a 8 y ) ,
which will be such that its denominator is canonically associated with an elliptic curve. Again, we can immediately calculate the j -invariant of that elliptic curve. If one calculates the telescoper of this eight-parameter family of rational functions (14), one finds that this telescoper is an order-two linear differential operator with pullbacked hypergeometric solutions of the form:
A ( x ) · 2 F 1 1 12 , 5 12 ; 1 ; H ,
where A ( x ) is an algebraic function, and where again, the pullback-Hauptmodul, H = 1728 / j , precisely corresponds to the j -invariant of the elliptic curve.
More generally, in seeking for planar elliptic curves, one can, for the given values of two integers M and N , look for planar algebraic curves:
n = 0 n = N m = 0 m = M a m , n · x n y m = 0 ,
defined by the set of a m , n ’s that are equal to zero, apart from N homogeneous parameters a m , n being, as in (10), (11), or (13), independent parameters. Finding such an N -parameter family of (planar) elliptic curves automatically provides an N -parameter family of rational functions such that their telescopers have a pullbacked 2 F 1 hypergeometric solution, which we can simply deduce from the j -invariant of that elliptic curve.
Recalling the results of Section 2.2, the natural question to ask now is whether it is possible to find families of such (planar) elliptic curves that depend on more than 10 independent parameters.
Before addressing this question, let us recall the concept of birationally equivalent elliptic curves. Let us consider the monomial transformation:
( x , y ) ( x M y N , x P y Q ) ,
where M , N , P , Q are integers, such that M · Q P · N = 1 ; then, its compositional inverse is the monomial transformation:
( x , y ) ( x Q y N , y M x P ) .
This monomial transformation (17) is thus a birational (This transformation is rational and its compositional inverse is also rational (here, monomial)) transformation. A birational transformation transforms an elliptic curve, such as (12), into another elliptic curve with the same j -invariant: these two elliptic curves are said to be birationally equivalent. In the case of the birational and monomial transformation (17), the elliptic curve (12) is changed into (one can easily verify for particular values of the M , N , P , Q , and a k ’s, using with(algcurves) in Maple, that the j -invariants of (12) and (19) are actually equal.):
a 1 · x 4 M y 4 N + a 2 · x 3 M y 3 N + a 3 · x 2 M y 2 N + a 4 · x M y N + a 5 + a 6 · x 2 M + P y 2 N + Q + a 7 · x 2 P y 2 Q + a 8 · x P y Q + a 9 · x M + P y N + Q = 0 .
With this kind of birational monomial transformation (17), we see that one can obtain families of elliptic curves (19) of arbitrary large degrees in x and y . Consequently, one can find 9- or 10-parameter families of rational functions of arbitrary large degrees yielding pullbacked 2 F 1 hypergeometric functions. There is no constraint on the degree of the planar algebraic curves (19): the only relevant question concerns the maximum number of (linearly) independent parameters in families of planar elliptic curves, which is shown to be 10. The demonstration (we thank Josef Schicho for providing this demonstration) is sketched in Appendix B.

3.3. Pullbacked 2 F 1 Functions for Higher Genus Curves: Monomial Transformations

Let us recall another important point. We have already remarked in [1,2] that once we have an exact result for a diagonal of a rational function of three variables R ( x , y , z ) , we immediately obtain another exact result for the diagonal of the rational function R ( x n , y n , z n ) for any positive integer n . As a result, we obtain a new expression for the diagonal changing x into x n . In fact, this is also a result on the telescoper of the rational function R ( x , y , z ) : the telescoper of the rational function R ( x n , y n , z n ) is the x x n pullback of the telescoper of the rational function R ( x , y , z ) . Having a pullbacked 2 F 1 solution for the telescoper of the rational function R ( x , y , z ) (respectively, the diagonal of the rational function R ( x , y , z ) ), we will immediately deduce a pullbacked 2 F 1 solution for the telescoper of the rational function R ( x n , y n , z n ) (respectively, the diagonal of the rational function R ( x n , y n , z n ) ).
Along this line, let us change, in the rational function (1), ( x , y , z ) into ( x 2 , y 2 , z 2 ) :
R 2 ( x , y , z ) = 1 a + b 1 x 2 + b 2 y 2 + b 3 z 2 + c 1 y 2 z 2 + c 2 x 2 z 2 + c 3 x 2 y 2 + d y 4 z 2 + e z 2 x 4 .
The diagonal of this new rational function (20) will be the pullbacked 2 F 1 exact expression (2), where we change x x 2 . The intersection of the algebraic surface corresponding to the vanishing condition of the denominator of the new rational function (20), with the hyperbola p = x y z (i.e., z = p x y ), is nothing but Equation (10), where we have changed ( x , y ; p ) into ( x 2 , y 2 ; p 2 ) :
a x 2 y 2 + b 1 x 4 y 2 + b 2 x 2 y 4 + b 3 p 2 + c 1 p 2 y 2 + c 2 p 2 x 2 + c 3 x 4 y 4 + d p 2 y 4 + e p 2 x 4 = 0 ,
which is no longer (if we perform the same calculations with the 10-parameter rational function, (4) we get an algebraic curve of genus 10 instead of 9 ) an elliptic curve, but a curve of genus 9 .
With that example, we see that classical modular form results, or pullbacked 2 F 1 exact expressions such as (2), can actually emerge from higher genus curves like (21). As far as these diagonals, or telescopers, of rational function calculations are concerned, higher genus curves like (21) must in fact be seen as “almost” elliptic curves up to an x x n covering.
Such results for monomial transformations like ( x , y , z ) ( x n , y n , z n ) can, in fact, be generalized to more general (not birational (in contrast with transformations such as (17))) monomial transformations. This is sketched in Appendix C.

3.4. Changing the Parameters into Functions of the Product p = x y z

All these results for the parameterized families of rational functions can be drastically generalized when one remarks that allowing any of these parameters to be a rational function of the product p = x y z also yields the previous pullbacked 2 F 1 exact expression, as in (2), where the parameter is changed into that rational function of x (see [1]). Let us consider a simple (two-parameter) illustration of this general result. Let us consider a subcase of the previous 9- or 10-parameter families, introducing, for example, the two-parameter rational function:
1 1 + 2 x + b 2 · y + 5 y z + x z + c 3 · x y .
The diagonal of this rational function (22) is the pullbacked hypergeometric function:
1 P 2 ( x ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 43200 · x 4 · P 4 ( x ) P 2 ( x ) 3 ,
where:
P 2 ( x ) = 1 8 · ( b 2 + 10 ) · x + 8 · ( 2 b 2 2 20 b 2 + 15 c 3 + 200 ) · x 2 ,
and:
P 4 ( x ) = 675 c 3 4 · x 4 + 4 c 3 2 · ( b 2 + 10 ) · ( 4 b 2 2 100 b 2 + 45 c 3 + 400 ) · x 3 + ( 64 b 2 4 32 b 2 3 c 3 8 b 2 2 c 3 2 1280 b 2 3 + 1280 b 2 2 c 3 460 b 2 c 3 2 5 c 3 3 + 6400 b 2 2 3200 b 2 c 3 800 c 3 2 ) · x 2 ( b 2 + 10 ) · ( 32 b 2 2 16 b 2 c 3 c 3 2 ) · x + 2 b 2 · ( 2 b 2 c 3 ) ,
Let us now consider the previous rational function (22) where the two parameters b 2 and c 3 become some rational functions of the product, p = x y z , for instance:
b 2 ( p ) = 1 + 3 p 1 + 7 p 2 , c 3 ( p ) = 1 + p 2 1 + 2 p where : p = x y z .
The new corresponding rational function becomes more involved, but one can easily calculate the telescoper of this new rational function of three variables x , y and z , and find that it is, again, an order-two linear differential operator having the pullbacked hypergeometric solution (23) where b 2 and c 3 are now replaced by ( p is now x ) the functions:
b 2 ( x ) = 1 + 3 x 1 + 7 x 2 , c 3 ( x ) = 1 + x 2 1 + 2 x .
In that case, one obtains a diagonal that is the pullbacked hypergeometric solution:
( 1 + 2 x ) 1 / 4 · ( 1 + 7 x 2 ) 1 / 4 · q 8 1 / 4 × 2 F 1 1 12 , 5 12 ; 1 ; 43200 · x 4 · ( 1 + 7 x 2 ) 2 · q 20 ( 1 + 2 x ) · q 8 3 ,
where q 8 and q 20 are two polynomials with integer coefficients of degree 8 and 20 in x . The exact expression (28) is nothing but (23) (with (24) and (25)), where b 2 and c 3 have been replaced by the rational functions (27). Similar calculations can be performed for the more general rational functions (1) or (4), when all the (9 or 10) parameters are more involved rational functions.
When performing our creative telescoping symbolic calculations using the HolonomicFunctions package [3], such results may look quite impressive. From the algebraic geometry viewpoint, it is almost tautological (An algebraic geometer will probably see this as a trivial remark: diagonalization is an algebraic procedure and nothing really happens to the coefficients. Therefore if one replaces the coefficients by anything else, one will find those replaced coefficients in the end result.), if one takes for granted the result of our previous Section 3.1 and Section 3.2, namely, that the pullbacked hypergeometric solution of the telescoper corresponds to the Hauptmodul 1728 / j , where j is the j -invariant of the elliptic curve corresponding to the intersection of the algebraic surface corresponding to the vanishing condition of the denominator, with the hyperbola p = x y z : this calculation of the j -invariant is performed for p fixed, and arbitrary (9 or 10) parameters a , b 1 , . It is clearly possible to force the parameters to be functions (The functions should be rational functions if one wants to stick with the diagonals and telescopers of rational functions, but the result remains valid for algebraic functions, or even transcendental functions with reasonable Taylor series expansions at x = 0 : for instance, for 2 F 1 hypergeometric functions, one obtains a differentially algebraic function corresponding to the composition of 2 F 1 hypergeometric functions.) of p , the j -invariant being changed accordingly. Of course, in that case, the parameters in the rational function are the same functions, but of the product p = x y z .
One thus obtains pullbacked hypergeometric solutions (classical modular forms) for an (unreasonably) large set of rational functions in three variables, namely, the families of the rational functions (1) or (4), but now where the 9 or 10 parameters are 9 or 10 totally arbitrary rational functions (with Taylor series expansions) of the product p = x y z .
We see experimentally that changing the parameters of the rational function into functions actually works for the diagonals of rational functions, as well as for the solutions of telescopers of rational functions, depending on the parameters.

4. Creative Telescoping on Rational Functions of More than Three Variables Associated with Products or Foliations of Elliptic Curves

Let us show that such an algebraic geometry approach to creative telescoping can be generalized to the rational functions of more than three variables, when the vanishing condition of the denominator can be associated with the products of elliptic curves, or more generally, algebraic varieties with foliations in elliptic curves.
The telescoper of the rational function in the four variables x , y , z and w :
x y z ( 1 + z ) 2 x · ( 1 x ) · ( x x y z w ) · y · ( 1 y ) · ( y x y z w ) ,
gives an order-three self-adjoint linear differential operator which is, thus, the symmetric square of an order-two linear differential operator. The latter has the pullbacked hypergeometric solution:
S 1 ( x ) = ( 1 x + x 2 ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 27 4 · x 2 · ( 1 x ) 2 ( x 2 x + 1 ) 3 = 2 F 1 1 2 , 1 2 ; 1 ; x .
In [18], we underlined the difference between the diagonal of a rational function and solutions of the telescoper of the same rational function. In this case, the diagonal of the rational function (29) is zero (The reason for this is that the integration takes place over a cycle that is homologically equivalent to the trivial cycle. The cycle becomes trivial after taking the limit p 0 . Integrals over non-vanishing cycles usually give logarithms of p , as with the second solution to the hypergeometric function 2 F 1 ( 1 / 2 , 1 / 2 ; 1 ; x ) .), and is thus different from the pullbacked hypergeometric solution (30), which is a “period” [37] of the algebraic variety corresponding to the denominator over some (non-vanishing (diagonals of rational functions correspond to periods over vanishing cycles [5,38])) cycle. From now on, we will have a similar situation in most of the following examples of this paper.
This example is a simple illustration of what we expect for the products of elliptic curves, or of algebraic varieties with foliations in elliptic curves. Introducing the product p = x y z w , the vanishing condition of the denominator of the rational function (29) yields the surface S ( x , y , z ) = 0 :
( 1 + z ) 2 x · ( 1 x ) · ( x p ) · y · ( 1 y ) · ( y p ) = 0 .
For fixed p and fixed y , Equation (31) can be seen as an algebraic curve:
( 1 + z ) 2 λ · x · ( 1 x ) · ( x p ) = 0 with : λ = y · ( 1 y ) · ( y p ) .
For fixed p and fixed y , λ can be considered as a constant, the algebraic curve (32) being an elliptic curve with an obvious Weierstrass form:
Z 2 x · ( 1 x ) · ( x p ) = 0 where : Z = 1 + z λ .
The j -invariant of (32), or (a shift z z + 1 or a rescaling z 2 z 2 λ does not change the j -invariant of the Weierstrass elliptic form) (33), is well-known and yields the Hauptmodul H :
H = 1728 j = 27 4 · p 2 · ( 1 p ) 2 ( p 2 p + 1 ) 3
For fixed p and fixed x , Equation (31) can be seen as an algebraic curve:
( 1 + z ) 2 μ · y · ( 1 y ) · ( y p ) = 0 for : μ = x · ( 1 x ) · ( x p ) ,
which is also an elliptic curve with an obvious Weierstrass form and the same Hauptmodul (34). This Hauptmodul is precisely the one that occurs in the pullbacked hypergeometric solution (30).
More generally, the rational function of the four variables x , y , z and w :
x y z ( 1 + z ) 2 x · ( 1 x ) · ( x R 1 ( p ) ) · y · ( 1 y ) · ( y R 2 ( p ) ) ,
where p = x y z w , and where R 1 ( p ) and R 2 ( p ) are two arbitrary rational functions (with Taylor series expansions) of the product p = x y z w , yields a telescoper that has an order-four linear differential operator that is the symmetric product (As the present paper belongs to the literature on symbolic computation and not on pure mathematics for algebraic geometers, we use the standard Maple (DEtools) terminology of symmetric powers and symmetric products of linear differential operators [46]. Note that “symmetric product” is not a proper mathematical name for this construction on the solution space; it is a homomorphic image of the tensor product. The (Maple/DEtools) reason for choosing the name symmetric_product is the resemblance with the function symmetric_power.) of two order-two linear differential operators having, respectively, the pullbacked hypergeometric solutions (30) where x is replaced by R 1 ( x ) and R 2 ( x ) . These two hypergeometric solutions thus have the two Hauptmodul pullbacks:
H 1 = 1728 j 1 = 27 4 · R 1 ( p ) 2 · ( 1 R 1 ( p ) ) 2 ( R 1 ( p ) 2 R 1 ( p ) + 1 ) 3 ,
H 2 = 1728 j 2 = 27 4 · R 2 ( p ) 2 · ( 1 R 2 ( p ) ) 2 ( R 2 ( p ) 2 R 2 ( p ) + 1 ) 3 ,
obtained by calculations similar to the ones previously performed on (31) but, now, for the Weierstrass form corresponding to the denominator (36).
A solution of the telescoper of (36) is thus the product of these two pullbacked hypergeometric functions. Let us give a simple illustration of this general result, with the next example.
The telescoper of the rational function in the four variables x , y , z and w :
x y z ( 1 + z ) 2 x · ( 1 x ) · ( x x y z w ) · y · ( 1 y ) · ( y 3 x y z w ) ,
corresponding to (36), with R 1 ( p ) = p and R 2 ( p ) = 3 p , gives an order-four linear differential operator that is the symmetric product of two order-two operators having, respectively, the pullbacked hypergeometric solution (30) and the solution (30) where the variable x has been changed into 3 x :
S 2 ( x ) = S 1 ( 3 x ) = ( 1 3 x + 9 x 2 ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 243 4 · x 2 · ( 1 3 x ) 2 ( 1 3 x + 9 x 2 ) 3 .

Creative Telescoping on Rational Functions of Five Variables Associated with Products or Foliations of Three Elliptic Curves

Let us now introduce the rational function in five variables x , y , z , v and w :
x y z v D ( x , y , z , v , w ) ,
where the denominator D ( x , y , z , v , w ) reads:
D p = ( 1 + v ) 2 x · ( 1 x ) · ( x p ) · y · ( 1 y ) · ( y 3 p ) · z · ( 1 z ) · ( z 5 p ) , where : p = x y z v w .
The telescoper of the rational function (41) of five variables gives (Such a creative telescoping calculation requires “some” computing time to achieve the result.) an order-eight linear differential operator, which is the symmetric product of three order-two linear differential operators having, respectively, the pullbacked hypergeometric solution (30), which is the solution (30) where x has been changed into 3 x , namely (40), and the solution (30), where x has been changed into 5 x :
S 3 ( x ) = S 1 ( 5 x ) = ( 1 5 x + 25 x 2 ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 675 4 · x 2 · ( 1 5 x ) 2 ( 1 5 x + 25 x 2 ) 3 .
In other words, the order-eight telescoper of the rational function (41) has the product S = S 1 · S 2 · S 3 , of (30), (40), and (43) as a solution. From an algebraic geometry viewpoint, this is a consequence of the fact that, for fixed p , the algebraic variety D p = 0 , where D p is given by (42), can be seen, for fixed y and z , as an elliptic curve E 1 of equation D y , z , p ( v , x ) = 0 , for fixed x and z , as an elliptic curve E 2 of equation D x , z , p ( v , y ) = 0 , and also, for fixed x and y , as an elliptic curve E 3 of equation D x , y , p ( v , z ) = 0 , the j -invariants j k , k = 1 , 2 , 3 , of these three elliptic curves E k yielding (in terms of p ), precisely, the three Hauptmoduls H k = : 1728 j k
27 4 · x 2 · ( 1 x ) 2 ( x 2 x + 1 ) 3 , 243 4 · x 2 · ( 1 3 x ) 2 ( 1 3 x + 9 x 2 ) 3 , 675 4 · x 2 · ( 1 5 x ) 2 ( 1 5 x + 25 x 2 ) 3 ,
occurring as pullbacks in the three S k ’s of the solution S = S 1 · S 2 · S 3 , of the telescoper of (41).
Remark 1.
Other examples of rational functions of three, four, and five variables where the denominators also correspond to Weierstrass (respectively, Legendre) forms, are displayed in Appendix D. They provide simple illustrations of rational functions where the denominator is associated with K3 surfaces (See the emergence of the product of elliptic curves from a Shioda–Inose structure on surfaces with Picard number 19 in [47]. In [47], Ling Long considers one-parameter families of K3 surfaces with generic Picard number 19. The existence of a Shioda–Inose structure implies that there is a one-parameter family of elliptic curves.), or Calabi–Yau threefolds. In these cases, the algebraic varieties have simple foliations in terms of two or three families of elliptic curves, and the solutions of the corresponding telescopers can be selected 3 F 2 and 4 F 3 hypergeometric functions (see (A28) in Appendix D), naturally associated with K3 surfaces and Calabi–Yau operators [27].

5. Creative Telescoping of Rational Functions in Three Variables Associated with Genus-Two Curves with Split Jacobians

In papers [17,18], dedicated to Heun functions that are solutions of telescopers of simple rational functions of three and four variables, we have obtained (See Equation (83) in Section 2.2 of [18].) an order-four telescoper of a rational function of three variables, which is the direct sum of two order-two linear differential operators, each having classical modular form solutions that can be written as pullbacked 2 F 1 hypergeometric solutions. Unfortunately, the intersection of the algebraic surface corresponding to the denominator of the rational function with the p = x y z hyperbola yields a genus-two algebraic curve.
Let us try to understand, in this section, how a genus-two curve can yield two classical modular forms. Let us first recall the results in Section 2.2 of [18].

5.1. Periods of Extremal Rational Surfaces

Let us recall the rational function in just three variables [18]:
R ( x , y , z ) = 1 1 + x + y + z + x y + y z x 3 y z .
Its telescoper is actually an order-four linear differential operator L 4 , which not only factorizes into two order-two linear differential operators, but is actually the direct sum (LCLM) of two (These two order-two linear differential operators, L 2 and M 2 , are not homomorphic) order-two linear differential operators, L 4 = L 2 M 2 . These two (non-homomorphic) order-two linear differential operators contain, respectively, the two pullbacked hypergeometric solutions:
S 1 = ( 1 + 9 x ) 1 / 4 · ( 1 + 3 x ) 1 / 4 · ( 1 + 27 x 2 ) 1 / 4 × 2 F 1 1 12 , 5 12 ; 1 ; 1728 · x 3 · ( 1 + 9 x + 27 x 2 ) 3 ( 1 + 3 x ) 3 · ( 1 + 9 x ) 3 · ( 1 + 27 x 2 ) 3 ,
and:
S 2 = 1 ( 1 + 4 x 2 x 2 36 x 3 + 81 x 4 ) 1 / 4 × 2 F 1 1 12 , 5 12 ; 1 ; 1728 · x 5 · ( 1 + 9 x + 27 x 2 ) · ( 1 2 x ) 2 ( 1 + 4 x 2 x 2 36 x 3 + 81 x 4 ) 3 .
The diagonal of (45) is actually the half-sum of the two series (46) and (47):
D i a g R ( x , y , z ) = S 1 + S 2 2 .
As far as our algebraic geometry approach is concerned, the intersection of the algebraic surface corresponding to the denominator of the rational function (45) with the hyperbola p = x y z gives the planar algebraic curve (corresponding to the elimination of the z variable by the substitution z = p x y ):
1 + x + y + p x y + x y + y p x y x 3 y p x y = 0 .
One easily finds that this algebraic curve is (for p fixed) a genus-two curve, and that this higher genus situation does not correspond to the “almost elliptic curves” described in Section 3.2, namely, an elliptic curve transformed by a monomial transformation. How can a “true” genus-two curve yield two j -invariants, namely, a telescoper with two Hauptmodul pullbacked 2 F 1 solutions? We are going to see that the answer is that the Jacobian of this genus-two curve (An algebraic geometer will probably recall that it is very well-known that a genus-two curve may have Jacobian isogenous to a product of elliptic curves. This is not the case in general. The genus-two curves that have a (non-constant) map to an elliptic curve have this property. Our purpose in Section 5.3 is to perform a creative telescoping calculation in such a selected situation.) is in fact isogenous to a product E × E of two elliptic curves (split Jacobian).

5.2. Split Jacobians

Let us first recall the concept of split Jacobian with a simple example. In [48], one has a crystal-clear example of a genus-two curve C :
y 2 ( x 3 + 420 x 5600 ) · ( x 3 + 42 x 2 + 1120 ) = 0 ,
such that its Jacobian J ( C ) is isogenous to a product of elliptic curves with j -invariants j 1 = 2 7 · 7 2 = 6272 and j 2 = 2 5 · 7 · 17 3 = 1100512 , corresponding to the following two values of the Hauptmodul, H = 1728 j : H 1 = 27 / 98 and H 2 = 54 / 34391 . Let us consider the genus-one elliptic curve:
v 2 = u 3 + 4900 u 2 + 7031500 u + 2401000000 ,
of j -invariant j = j 2 = 2 5 · 7 · 17 3 . We consider the following transformation (This transformation is rational but not birational. If it were birational, then it would preserve the genus. Here, one goes from genus one to genus two):
u = 882000 · ( x 14 ) x 3 + 420 x 5600 , v = 49000 · ( x 3 21 x 2 140 ) ( x 3 + 420 x 5600 ) 2 · y .
This change of variables (52) actually transforms the elliptic curve (51) in u and v into the genus-two curve (50) in x and y . This provides a simple example of a genus-two curve with split Jacobian through K3 surfaces.
More generally, let us consider the Jacobian of a genus-two curve C . The Jacobian is simple if it does not contain a proper abelian subvariety, otherwise the Jacobian is reducible, or decomposable or “split”. For this latter case, the only possibility for a genus-two curve is that its Jacobian is isogenous to a product E × E of two elliptic curves (along these lines, see also the concepts of Igusa–Clebsch invariants and Hilbert modular surfaces [48,49,50,51]). Equivalently, there is a degree-n map C E to some elliptic curves. Classically, such pairs (One also has an anti-isometry Galois invariant E E under Weil pairing. The decomposition corresponds to real multiplication by the quadratic ring of discriminant n 2 ) C , E arose in the reduction of hyperelliptic integrals to elliptic ones [48]. The j -invariants correspond, here, to the two elliptic subfields: see [48].

5.3. Creative Telescoping on Rational Functions in Three Variables Associated with Genus-Two Curves with Split Jacobians: A Two-Parameter Example

Let us now consider the example with two parameters, a and b, given in Section 4.5 on page 12 of [48]. Let us substitute the rational parameterization (see also [52], Section 6, page 48):
u = x 2 x 3 + a x 2 + b x + 1 , v = y · ( x 3 b x 2 ) ( x 3 + a x 2 + b x + 1 ) 2 ,
in the elliptic curve:
R · v 2 = R · u 3 + 2 · ( a b 2 6 a 2 + 9 b ) · u 2 + ( 12 a b 2 ) · u 4 ,
where:
R = 4 · ( a 3 + b 3 ) a 2 b 2 18 a b + 27 .
This gives the genus-two curve C a , b ( x , y ) = 0 , with:
C a , b ( x , y ) = R · y 2 + ( 4 x 3 + b 2 x 2 + 2 b x + 1 ) · ( x 3 + a x 2 + b x + 1 ) .
The j -invariant of the elliptic curve (54) gives the following exact expression for the Hauptmodul H = 1728 j :
H = 108 · ( b 3 ) 3 · ( 4 a 3 + 4 b 3 a 2 b 2 18 a b + 27 ) 2 · ( b 2 + 3 b + 9 ) 3 ( a 2 b 4 + 12 b 5 126 a b 3 + 216 b a 2 + 405 b 2 972 a ) 3 .
Let us consider the telescoper of the rational function of three variables x y / D a ( x , y , z ) , where the denominator D a ( x , y , z ) is C a , b ( x , y ) , given by (56), but for b = 3 + x y z :
D a ( x , y , z ) = C a , 3 + x y z ( x , y ) = x 6 y 3 z 3 + x 7 y 2 z 2 + 4 x 3 y 5 z 3 + 9 x 5 y 2 z 2 + 6 x 6 y z + 3 x 4 y 2 z 2 + 36 y 4 x 2 z 2 + 6 x 5 y z + 4 x 6 + 27 x 4 y z + 9 x 5 + 18 x 3 y z + 108 x y 3 z + 18 x 4 + 3 x 2 y z + 32 x 3 + 27 x 2 + 135 y 2 + 9 x + 1 + ( x 6 y 2 z 2 + 6 x 5 y z + 2 x 4 y z + 4 x 5 18 x y 3 z + 9 x 4 + 6 x 3 + x 2 54 y 2 ) · a y 2 · ( x y z + 3 ) 2 · a 2 + 4 y 2 · a 3 .
This telescoper of the rational function:
R a ( x , y , z ) = x y D a ( x , y , z ) ,
is an order-four linear differential operator L 4 , which is actually the direct sum, L 4 = L C L M ( L 2 , M 2 ) = L 2 M 2 , of two order-two linear differential operators, having two pullbacked hypergeometric solutions. One finds out that one of the two pullbacks precisely corresponds to the Hauptmodul H given by (57) for b = 3 + x .
Let us consider the a = 3 subcase (the discriminant in b of 4 a 3 + 4 b 3 a 2 b 2 18 a b + 27 reads: ( a 3 ) 3 · ( a 2 + 3 a + 9 ) 3 ; consequently, the exact expressions are simpler at a = 3 ). For a = 3 , the Hauptmodul H = 1728 j , given by (57) becomes for b = 3 + x :
H = 4 · x · ( 27 + 4 x ) 2 · ( x 2 + 9 x + 27 ) 3 ( 9 + x ) 3 · ( 4 x 2 + 27 x + 27 ) 3 .
The telescoper of the rational function (59), with D a ( x , y , z ) given by (58) for a = 3 , is an order-four linear differential operator, which is the direct sum of two order-two linear differential operators L 4 = L C L M ( L 2 , M 2 ) = L 2 M 2 , with these two order-two linear differential operators having the pullbacked hypergeometric solutions:
( 27 + 4 x ) 1 / 2 · x 5 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 1 + 27 4 x ,
for L 2 , and:
3 + x ( 9 + x ) 1 / 4 · ( 4 x 2 + 27 x + 27 ) 1 / 4 · x 3 / 2 · ( 27 + 4 x ) 1 / 2 × 2 F 1 1 12 , 5 12 ; 1 ; 4 · x · ( 27 + 4 x ) 2 · ( x 2 + 9 x + 27 ) 3 ( 9 + x ) 3 · ( 4 x 2 + 27 x + 27 ) 3 ,
for M 2 , where we see clearly that the Hauptmodul in (62) is precisely the Hauptmodul (60). The Jacobian of the genus-two curve is a split Jacobian corresponding to the product E 1 × E 2 of two elliptic curves; the j -invariant of the second elliptic curve corresponds to the Hauptmodul H = 1728 j given by (57) when the j -invariant of the first elliptic curve reads:
j 1 = 6912 x 27 + 4 x ,
corresponding to the Hauptmodul 1728 j 1 = 1 + 27 4 x in (61). This second invariant is, as it should, exactly the j -invariant of the second elliptic curve E , given page 48 in [52]:
j ( E ) = 256 · ( 3 b a 2 ) 3 4 a 3 c a 2 b 2 18 a b c + 4 b 3 + 27 c 2 ,
for c = 1 , a = 3 and b = 3 + x .

5.4. Creative Telescoping on Rational Functions of Three Variables Associated with Genus-Two Curves with Split Jacobians: A Simple Example

Another simpler example of a genus-two curve with pullbacked 2 F 1 solution (not a product of pullbacked 2 F 1 ) of the telescoper can be given if one considers the genus-two algebraic curve C p ( x , y ) = 0 given in Lemma 7 of [53] (see also [54,55]):
C p ( x , y ) = x 5 + x 3 + p · x y 2 .
Let us introduce the rational function x y / D ( x , y , z ) , where the denominator D ( x , y , z ) is given by:
D ( x , y , z ) = C ( p = x y z ) ( x , y ) = x 5 + x 3 + x 2 y z y 2 .
The telescoper of this rational function is an order-two linear differential operator that has the two hypergeometric solutions:
x 1 / 4 · 2 F 1 1 8 , 5 8 ; 3 4 ; 4 x
which is a Puiseux series at x = 0 and:
x 1 / 4 · 2 F 1 1 8 , 5 8 ; 1 ; 1 4 x .
These two hypergeometric solutions can be rewritten as (the fact that 2 F 1 1 8 , 5 8 ; 1 ; z can be rewritten as 2 F 1 1 12 , 5 12 ; 1 ; H ( z ) where the Hauptmodul H ( z ) is solution of a quadratic equation is given in Equation (H.14) of Appendix H of [18]):
A ( x ) · 2 F 1 1 12 , 5 12 ; 1 ; 1728 J ,
where the j -invariant J , in the Hauptmodul 1728 J in (69), corresponds exactly to the degree-two elliptic subfields:
J 2 128 · ( 2000 x 2 + 1440 x + 27 ) ( 1 4 x ) 2 · J 4096 · ( 100 x 9 ) 3 ( 1 4 x ) 3 = 0 ,
given in the first equation of page 6 of [53].
Remark 2.
In contrast to the previous example of Section 5.3, where we had two j -invariants corresponding to the two order-two linear differential operators L 2 and M 2 of the direct-sum decomposition of the order-four telescoper, here, we have just one order-two telescoper, which is enough to “encapsulate” the two j -invariants (70), since they are Galois-conjugate.

6. Rational Functions with Tri-Quadratic Denominator and N -Quadratic Denominator

We attempt to find telescopers of rational functions corresponding to (factors of) linear differential operators of “small” orders; for instance, order-two linear differential operators with pullbacked 2 F 1 hypergeometric functions, classical modular forms, or their modular generalizations (order-four Calabi–Yau linear differential operators [27], etc.). As we saw in the previous sections, this corresponds to the fact that the denominator of these rational functions is associated with an elliptic curve, or the products of elliptic curves, with K3 surfaces or with threefold Calabi–Yau manifolds corresponding to algebraic varieties with foliations in elliptic curves (Even though K3 surfaces or threefold Calabi-Yau manifolds are not abelian varieties, the Weierstrass–Legendre forms introduced in Appendix D amount to saying that K3 surfaces can be “essentially viewed” (as far as creative telescoping is concerned) as foliations in two elliptic curves, and threefold Calabi–Yau manifolds as foliations in three elliptic curves), since this paper attempts to reduce the differential algebra creative telescoping calculations to effective algebraic geometry calculations. (One has birational automorphisms in projective spaces [56,57], but since this paper is dedicated to (efficient) formal calculations, we work exclusively in affine coordinates (see for instance (A41), (A42), and (A43) below.) For algebraic geometers, an elliptic curve is a smooth complete genus-one curve with a choice of a base point. Here our elliptic curves are, in fact, an affine piece of a genus-one curve with no base point, but this does not really matter, because the j -invariant, which is all we care about in these kinds of creative telescoping calculations, is determined by that much information and structures, we want to focus on rational functions with denominators that correspond to selected algebraic varieties [44,58], beyond algebraic varieties corresponding to products of elliptic curves or foliations in elliptic curves (K3 surfaces, threefold Calabi–Yau manifolds, and higher curves with split Jacobian corresponding to products of elliptic curves,…), namely, algebraic varieties with an infinite number of birational automorphisms (the best explicit illustration of this situation emerges in integrable models [44,58,59,60]). This infinite number of birational symmetries excludes algebraic varieties of the “general type” (with a finite number. (There are even precise bounds for the number of automorphisms. The upper bound is 84 ( g 1 ) for curves of genus g , and these bounds have been extensively studied in higher dimensions [61,62,63]) of birational symmetries.) For algebraic surfaces, this amounts to discarding the surfaces of “general type” that have Kodaira dimension two, focusing on Kodaira dimension one (elliptic surfaces), Kodaira dimension zero (abelian surfaces, hyperelliptic surfaces, K3 surfaces, and Enriques surfaces), or even Kodaira dimension (ruled surfaces and rational surfaces).
In contrast with algebraic curves, where one can easily and very efficiently calculate the genus of the curves to discard algebraic curves of higher genus and, in the case of genus one, obtain the j -invariant using formal calculations (use with(algcurves) in Maple, and the commands “genus” and “j_invariant”), it is, in practice, quite difficult to see for higher dimensional algebraic varieties, that the algebraic variety is not of the “general type”, because it has an infinite number of birational symmetries. For these (low Kodaira dimension) “selected cases” that we are interested in, calculating the generalization of the j -invariant (Igusa–Shioda invariants, etc.) is quite hard.
Along this line, we want to underline that there exists a remarkable set of algebraic surfaces, namely the algebraic surfaces corresponding to tri-quadratic equations:
m = 0 , 1 , 2 n = 0 , 1 , 2 l = 0 , 1 , 2 a m , n , l · x m y n z l = 0 ,
depending on 27 = 3 3 parameters a m , n , l . More generally, one can introduce algebraic varieties corresponding to N -quadratic equations:
m 1 = 0 , 1 , 2 m 2 = 0 , 1 , 2 m N = 0 , 1 , 2 a m 1 , m 2 , , m N · x 1 m 1 x 2 m 2 x N m N = 0 .
With these tri-quadratic (71), or N -quadratic (72) equations, we will see, in Appendix E.1 and Appendix E.2, that we have automatically (selected) algebraic varieties that are not of the “general type”, having an infinite number of birational symmetries, which is precisely our requirement for the denominator of rational functions with remarkable telescopers. (Telescopers with factors of “small enough” order, possibly yielding classical modular forms, Calabi–Yau operators, etc. Rational functions with denominators of the “general type” will yield telescopers of very large orders.)
Let us first, as a warm-up, consider, in the next subsection, a remarkable example of tri-quadratic (71), where the underlying foliation in elliptic curves is crystal clear.

6.1. Rational Functions with Tri-Quadratic Denominator Simply Corresponding to Elliptic Curves

Let us first recall the tri-quadratic equation in three variables x , y and z :
x 2 y 2 z 2 2 · M · x y z · ( x + y + z ) + 4 · M · ( M + 1 ) · x y z + M 2 · ( x 2 + y 2 + z 2 ) 2 M 2 · ( x y + x z + y z ) = 0 ,
already introduced in Appendix C of [64]. This algebraic surface, symmetric in x , y , and z , can be seen for z (respectively, x or y ) fixed, as an elliptic curve whose j -invariant is independent of z , yielding the corresponding Hauptmodul:
H = 1728 j = 27 · M 2 · ( M 1 ) 2 4 · ( M 2 M + 1 ) 3 .
This corresponds to the fact that this algebraic surface (73) can be seen as a product of two instances of the same elliptic curve with the Hauptmodul (74). This is a consequence of the fact that, introducing x = s n ( u ) 2 , y = s n ( v ) 2 and z = s n ( u + v ) 2 , and M = 1 / k 2 , this algebraic surface (73) corresponds to the well-known formula for the addition on an elliptic sine (see Equation (C.3) in Appendix C of [64]):
s n ( u + v ) = s n ( u ) c n ( v ) d n ( v ) + s n ( v ) c n ( u ) d n ( u ) 1 k 2 s n ( u ) 2 s n ( v ) 2 .
For M = x y z w , the LHS of the tri-quadratic Equation (73) yields a polynomial of four variables x , y , z , and w , which we denote T ( x , y , z , w ) :
T ( x , y , z , w ) = x 2 y 2 z 2 2 · x 2 y 2 z 2 w · ( x + y + z ) + 4 · ( x y z w + 1 ) · x 2 y 2 z 2 w + x 2 y 2 z 2 w 2 · ( x 2 + y 2 + z 2 ) 2 x 2 y 2 z 2 w 2 · ( x y + x z + y z ) .
The telescoper of the rational function in four variables x , y , z , and w ,
x y z T ( x , y , z , w ) ,
is an order-three (self-adjoint) linear differential operator that is the symmetric square of the order-two linear differential operator, having the following pullbacked 2 F 1 hypergeometric solution:
x 1 / 2 · ( x 2 x + 1 ) 1 / 4 × 2 F 1 1 12 , 5 12 ; 1 ; 27 · x 2 · ( x 1 ) 2 4 · ( x 2 x + 1 ) 3 .
As it should, the Hauptmodul in (78) is the same as the Hauptmodul (74). The algebraic surface (73) can be seen as the product of two instances of the same elliptic curve with the Hauptmodul (74): As expected, the solution of the order-three telescoper is the square of the pullbacked 2 F 1 hypergeometric function (78) with that Hauptmodul.
More generally, we can also consider another tri-quadratic equation of three variables x , y , and z , and two parameters M and N :
x 2 y 2 z 2 2 M · x y z · ( x + y + z ) + N · x y z + M 2 · ( x 2 + y 2 + z 2 ) 2 M 2 · ( x y + x z + y z ) = 0 .
This surface, symmetric in x , y , and z , can be seen for z (respectively, x or y ) fixed as an elliptic curve whose j -invariant is again independent of z , yielding the corresponding Hauptmodul:
H = 1728 j = 1728 · M 6 · ( 64 M 3 N 2 ) ( 48 M 3 N 2 ) 3 .
Let us consider the following change of variables M = m 2 and N = 8 · m 3 + p in (79). For p = x y z w , the LHS of the tri-quadratic Equation (79) yields a polynomial in four variables x , y , z , and w , which we denote T m ( x , y , z , w ) :
T m ( x , y , z , w ) = x 2 y 2 z 2 2 m 2 · x y z · ( x + y + z ) + ( 8 · m 3 + x y z w ) · x y z + m 4 · ( x 2 + y 2 + z 2 ) 2 m 4 · ( x y + x z + y z ) .
For z (respectively, x or y ) fixed, the corresponding Hauptmodul (80) reads:
H = 1728 · m 12 · p · ( 16 m 3 + p ) ( 16 m 6 + 16 m 3 · p + p 2 ) 3 .
The telescoper of the rational function in four variables x , y , z , and w ,
x y z T m ( x , y , z , w ) ,
is an order-three (self-adjoint) linear differential operator that is the symmetric square of an order-two linear differential operator, having the following pullbacked 2 F 1 hypergeometric solution:
( 16 m 6 + 16 m 3 · x + x 2 ) 1 / 4 · × 2 F 1 1 12 , 5 12 ; 1 ; 1728 · m 12 · x · ( 16 m 3 + x ) ( 16 m 6 + 16 m 3 · x + x 2 ) 3 .
As it should, the Hauptmodul in (84) is the same as the Hauptmodul (82). The algebraic surface (79) can be seen as the product of two instances of the same elliptic curve with the Hauptmodul (80) (or (82)). As expected, the solution of the order-three telescoper is the square of the pullbacked 2 F 1 hypergeometric function (84) with the Hauptmodul (82).
Remark 3.
Let us perform some (slight) deformations of the rational function (77), changing the first 2 coefficient in (76) into a 3 coefficient. One thus considers the polynomial T ( x , y , z , w ) :
T ( x , y , z , w ) = x 2 y 2 z 2 3 · x 2 y 2 z 2 w · ( x + y + z ) + 4 · ( x y z w + 1 ) · x 2 y 2 z 2 w + x 2 y 2 z 2 w 2 · ( x 2 + y 2 + z 2 ) 2 · x 2 y 2 z 2 w 2 · ( x y + x z + y z ) .
The telescoper of the rational function in four variables,
x y z T ( x , y , z , w ) ,
is an (irreducible) linear differential operator L 4 of (only) order four, which is non-trivially homomorphic to its adjoint. (Its exterior square has a rational solution. However this order-four linear differential operator is not MUM (maximum unipotent monodromy [27,65,66]).) A priori, we cannot exclude the fact that L 4 could be homomorphic to the symmetric cube of a second-order linear differential operator, or to a symmetric product of two second-order operators. Furthermore, it could also be, in principle, that these second-order operators admit classical modular forms as solutions (pullbacks of special 2 F 1 hypergeometric functions). However, these options can both be excluded by using some results from differential Galois theory [67], specifically from [68] (Prop. 7, p. 50) for the symmetric cube case, and from [68] (Prop. 10, p. 69) for the symmetric product case; see also [69] (§3). Indeed, if L 4 were either a symmetric cube or a symmetric product of order-two operators, then its symmetric square would contain a (direct) factor of order 3 or 1. This is ruled out by a factorization procedure which shows that the symmetric square of L 4 is (LCLM-)irreducible.
This example does not correspond to an addition formula such as (75), but the polynomial T ( x , y , z , w ) still corresponds to a tri-quadratic. Consequently, it is an algebraic variety with an infinite number of birational automorphisms, as shown in Appendix E.1.

6.2. Rational Functions with Tri-Quadratic Denominator: Fricke Cubics Examples Associated with Painlevé VI Equations

Let us consider other simple examples of tri-quadratic surfaces that occur in different domains of mathematics and physics.
Among the Fricke families of cubic surfaces, the family [70,71,72]:
x y z + x 2 + y 2 + z 2 + b 1 x + b 2 y + b 3 z + c = 0 ,
of affine cubic surfaces parameterized by the four constants ( b 1 , b 2 , b 3 , c ) is known [71] to be a deformation of a D 4 singularity that occurs at the symmetric (Manin’s) case b 1 = b 2 = b 3 = 8 , c = 28 .
Among the symmetric b 1 = b 2 = b 3 cases, some selected sets of the four constants ( b 1 , b 2 , b 3 , c ) emerge: the Markov cubic b 1 = b 2 = b 3 = c = 0 , Cayley’s nodal cubic b 1 = b 2 = b 3 = 0 , c = 4 , Clebsch’s diagonal cubic b 1 = b 2 = b 3 = 0 , c = 20 , and Klein’s cubic b 1 = b 2 = b 3 = 1 , c = 0 .
Some of these symmetric cubics can be seen as the monodromy manifold of the Painlevé VI equation (see Equation (1.7) in [73], see also Equations (1.2) and (1.4) in [72]): the Picard–Hitchin cases ( 0 , 0 , 0 , 4 ) , ( 0 , 0 , 0 , 4 ) , ( 0 , 0 , 0 , 32 ) , Kitaev’s cases ( 0 , 0 , 0 , 0 ) , ( 8 , 8 , 8 , 64 ) , and especially Manin’s case ( 8 , 8 , 8 , 28 ) .
Let us consider the (symmetric) rational function in three variables, x , y , and z [71]:
R ( x , y , z ) = 1 x 2 + y 2 + z 2 + x y z + c ,
that takes into account the other Picard–Hitchin cases (as well as the Markov cubic b 1 = b 2 = b 3 = c = 0 , Cayley’s nodal cubic b 1 = b 2 = b 3 = 0 , c = 4 , and Clebsch’s diagonal cubic b 1 = b 2 = b 3 = 0 , c = 20 cases) ( 0 , 0 , 0 , 4 ) , ( 0 , 0 , 0 , 4 ) , ( 0 , 0 , 0 , 32 ) . The rational function (88) has an order-two telescoper that has a simple pullbacked hypergeometric solution:
1 x + c · 2 F 1 1 3 , 2 3 ; 1 ; 27 x 2 ( x + c ) 3 = ( x + c ) 1 / 4 · q 3 ( x ) 1 / 4 · 2 F 1 1 12 , 5 12 , 1 ; 1728 · x 6 · p 3 ( x ) ( x + c ) 3 · q 3 ( x ) 3 ,
where (the values c = 0 and c = 4 are the only values such that the discriminant in x of p 3 ( x ) can be zero):
p 3 ( x ) = x 3 + 3 · ( c + 9 ) · x 2 + 3 · c 2 · x + c 3 , q 3 ( x ) = x 3 + 3 · ( c + 8 ) · x 2 + 3 · c 2 · x + c 3 ,
Eliminating z = p x y in the denominator of (88) gives the genus-four algebraic curve:
x 2 y 2 · ( x 2 + y 2 ) + ( p + c ) · x 2 y 2 + p 2 = 0 .
Again, the question is to see whether the Jacobian of this genus-four algebraic curve (88) could also correspond to a split Jacobian, with a j -invariant corresponding to the Hauptmodul in (89).

7. Telescopers of Rational Functions of Several Variables

Let us consider the rational function in four variables x , y , z , u :
R ( x , y , z , u ) = 1 1 + 3 y + z + 9 y z + 11 z 2 y + 3 u x .
The telescoper of this rational function of four variables is an order-two linear differential operator L 2 that has the pullbacked hypergeometric solution:
( 1 2592 x 2 ) 1 / 4 × 2 F 1 1 12 , 5 12 ; 1 ; 419904 · x 3 · ( 5 12 x 19440 x 2 + 2665872 x 3 ) ( 1 2592 x 2 ) 3 .
The diagonal of (91) is the expansion of this pullbacked hypergeometric function (92):
1 + 648 x 2 72900 x 3 + 1224720 x 4 330674400 x 5 + 23370413220 x 6 1276733858400 x 7 + 180019474034400 x 8 12013427240614800 x 9 +
If one considers the intersection of the vanishing condition of the denominator of (91) with the hyperbola p = x y z u , eliminating, for instance, u = p x y z in the vanishing condition of the denominator of (91), one gets a condition, independent of x , that corresponds to a genus-one curve:
11 y 2 z 3 + 9 y 2 z 2 + 3 y 2 z + y z 2 + y z + 3 p = 0 .
The Hauptmodul of this elliptic curve (94) reads:
H = 419904 · p 3 · ( 5 12 p 19440 p 2 + 2665872 p 3 ) ( 1 2592 p 2 ) 3 ,
which corresponds precisely to the Hauptmodul pullback in (92).
Remark 4.
The expansion (93) of (92) is not only the diagonal of the rational function R ( x , y , z , u ) in four variables (91); it is also the diagonal of the rational function of three variables R ( x , y , z , 1 ) . Actually, using Section 3, one sees easily that eliminating x = p y z in the vanishing condition of the denominator of R ( x , y , z , 1 ) gives exactly the same elliptic curve (94).
Let us now generalize the rational function (91) of four variables x , y , z , u , by introducing the rational function of N + 3 variables x , y , z , u 1 , u 2 , , u N :
R ( x , y , z , u 1 , u 2 , , u N ) = 1 1 + 3 y + z + 9 y z + 11 z 2 y + 3 x · u 1 u 2 u N .
The telescoper of this rational function of N + 3 variables is the same order-two telescoper as for (91), which has the pullbacked hypergeometric solution (92). Again, one can verify that the diagonal of (96) is the expansion (93) of the pullbacked hypergeometric function (a pure algebraic geometer will probably consider this result as trivial from the computational point of view, saying that the variety is a fiber bundle over a family of elliptic curves with constant fiber (see also below)) (92). If one considers the intersection of the vanishing condition of the denominator of (96) with the hyperbola p = x y z u 1 u 2 u N , eliminating, for instance, u N = p x y z u 1 u N 1 in the vanishing condition of the denominator of (96), one again obtains a condition, independent of x but also of u 1 , , u N , that corresponds to the genus-one curve (94):
11 y 2 z 3 + 9 y 2 z 2 + 3 y 2 z + y z 2 + y z + 3 p = 0 .
The Hauptmodul of this elliptic curve (97) is again equal to the Hauptmodul (95), which corresponds precisely to the Hauptmodul pullback in (92).
Other examples, corresponding to the simple polynomial deformations of (91), such that their diagonal is the pullbacked 2 F 1 hypergeometric function (92), are displayed in Appendix F. This (infinite) family of rational functions corresponds to a different algebraic geometry scenario: the “canonical” algebraic surface corresponding to the intersection of the vanishing condition of the denominator of the rational function with the hyperbola p = x y z , is foliated in (generically high genus) algebraic curves depending on the variable x . One sees (experimentally) that the Hauptmodul of the pullbacked 2 F 1 hypergeometric functions corresponds to the Hauptmodul of the x = 0 algebraic curve, which is an elliptic curve (the algebraic curves for other values of x are not necessarily elliptic curves, they can be algebraic curves of quite large genus). In contrast to the other examples and results of this paper, we have no algebraic geometry interpretation of this experimental result yet.

8. Conclusions

Diagonals of rational functions emerge quite naturally in lattice statistical mechanics [19,20]. This explains the frequent occurrence of modular forms, represented as pullbacked 2 F 1 hypergeometric functions [1,2] in lattice statistical mechanics [21,22,23,24,25,26,27].
We have shown that the results that we had obtained on diagonals of 9- and 10-parameter families of rational functions in three variables, using creative telescoping yielding classical modular forms expressed as pullbacked 2 F 1 hypergeometric functions [1,2], can be obtained much more efficiently by calculating the j -invariant of an elliptic curve that is canonically associated with the denominator of the rational functions. In the case where creative telescoping yields pullbacked 2 F 1 hypergeometric functions, we generalize this result to other families of rational functions of three, and even more than three, variables, when the denominator can be associated with products of elliptic curves or foliations in terms of elliptic curves, or when the denominator is associated with a genus-two curve with a split Jacobian corresponding to the products of elliptic curves.
We have seen different scenarios. In the first cases, we have considered denominators corresponding to products of elliptic curves: in these cases, the solutions of the telescoper were products of pullbacked 2 F 1 hypergeometric functions. We have also considered denominators corresponding to genus-two curves with split Jacobians isogenous to the products of two elliptic curves, and in these cases, the solutions of the telescoper were the sums of two pullbacked 2 F 1 hypergeometric functions, sometimes with one pullbacked 2 F 1 hypergeometric function being enough to describe the two Galois-conjugate j -invariants (see Section 5.4). We also considered the denominators corresponding to algebraic varieties where the Hauptmodul pullback in the pullbacked 2 F 1 hypergeometric functions emerges from a selected ( x = 0 , see Appendix F.1 and Appendix F.2) elliptic curve of the algebraic variety. We also encountered denominators corresponding to algebraic manifolds with an infinite set of birational automorphisms and elliptic curves foliation, no longer yielding classical modular forms represented as pullbacked 2 F 1 hypergeometric functions, but more general modular structures associated with selected linear differential operators such as Calabi–Yau linear differential operators [27,65] and their generalizations.
The creative telescoping method on a rational function is an efficient way to find the periods of an algebraic variety over all possible cycles (not only the vanishing cycles [5,38] corresponding to the diagonals of rational functions). The fact that the solution of the telescoper corresponds to “periods” [37] over all possible cycles is a simple consequence of the fact that creative telescoping corresponds to purely differential algebraic manipulations on the integrand independently of the cycles, thus being blind to analytical details. In this paper, we show that the final result emerging from such differential algebra procedures (which can be cumbersome when the result depends on 9 or 10 parameters) can be obtained almost instantaneously from a more fundamental intrinsic pure algebraic geometry approach, calculating, for instance, the j-invariant of some canonical elliptic curve. This corresponds to a shift Analysis Differential Algebra Algebraic Geometry. In algebraic geometry studies of more involved algebraic varieties than the product of elliptic curves, foliation in elliptic curves (Calabi–Yau manifolds, etc.) is often a tedious and/or difficult task (finding Igusa–Shioda invariants, etc.), and formal calculation tools are not always available or user-friendly. Ironically, for such involved algebraic varieties, creative telescoping may then become a simple and efficient tool to perform effective algebraic geometry studies.
Some work remains to be done. While in the present paper, we have focused on studying the relationship between the denominator of a rational function and its diagonal, it would be interesting to explore the impact that variations of the numerator have on the telescoper or on the diagonal. Then, one may expect that all these algebraic–geometric findings could be exploited in the design of novel and more efficient creative telescoping algorithms for multivariate rational functions.

Author Contributions

Formal analysis, Y.A., S.B. and C.K.; Writing—original draft, J.-M.M. All authors have read and agreed to the published version of the manuscript.

Funding

Austrian Science Fund (FWF): F5011-N15.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

J.-M.M. would like to thank G. Christol for many enlightening discussions on diagonals of rational functions. J.-M.M. would like to thank D. van Straten for several enlightening and effective algebraic geometry explanations. J.-M.M. would like to thank the School of Mathematics and Statistics of Melbourne University where part of this work has been performed. S.B. would like to thank the LPTMC and the CNRS for their kind support. We thank Josef Schicho for providing the demonstration of the results in Appendix B. We would like to thank A. Bostan for useful discussions on creative telescoping. Y.A. and C.K. were supported by the Austrian Science Fund (FWF): F5011-N15. We thank the Research Institute for Symbolic Computation for access to the RISC software packages. We thank M. Quaggetto for technical support. This work has been performed without any support of the ANR, the ERC, or the MAE; or any PES of the CNRS.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Diagonals of Rational Functions and Picard–Fuchs Equations

For simplicity, let us consider the rational functions of three variables and double integrals [74]. The diagonal of a rational function of three variables is obtained through its multi-Taylor expansion [19,20]:
R ( x , y , z ) = m n l a m , n , l · x m y n z l ,
by extracting the “diagonal” terms, i.e., the powers of the product p = x y z :
D i a g R ( x , y , z ) = m a m , m , m · p m .
Such diagonals are closely related to the integrals of rational functions. For example, D i a g R ( x , y , z ) is the constant term (in y , z ) in the infinite expansion:
R p y z , y , z = m , n , l 0 a m , n , l · p m y n m z l m ,
which can be represented by the integral [35]:
1 ( 2 π i ) 2 R p y z , y , z d y y d z z .
The diagonal (A2) is also the constant term (in y , z ) of:
R p y , y z , z = m , n , l 0 a m , n , l · p m y n m z l n ,
which is of the form:
1 ( 2 π i ) 2 N p ( y , z ) D p ( y , z ) d y y d z z ,
where the numerator N p ( y , z ) and the denominator D p ( y , z ) are polynomials. It is well-known that such integrals satisfy a linear differential equation with respect to p having rational functions in p as coefficients, called the Picard–Fuchs equation. (The order of this linear differential equation is generally equal to the rank of the algebraic de Rham cohomology of D p ( y , z ) = 0 . For curves of genus g , this rank is 2 g .) The problem of determining such linear differential equations was commenced by Griffiths [75] with the assumption that the variety D p ( y , z ) = 0 is smooth, but later techniques were developed to include examples with singular points [35,40]. The linear differential equations (Gauss–Manin systems and telescopers) occurring in integrable models [16,23,24] are of an order that is much larger than order two (Since it is well-known that the Picard–Fuchs equation corresponding to the (Weierstrass) elliptic curve corresponds to the hypergeometric function 2 F 1 ( 1 / 12 , 5 / 12 ; 1 ; 1 / J ) ) and almost never correspond to smooth varieties. Creative telescoping (for a detailed introduction to creative telescoping [36], see for instance [34]), and more specifically, the programs [3] corresponding to a fast approach to creative telescoping [42], are a powerful way to find these linear differential operators annihilating these diagonals of rational functions in the cases emerging naturally in theoretical physics, integrable models, and enumerative combinatorics, for which the order of the linear differential operators is quite large [16,23,24] and the variety D p ( y , z ) = 0 is (most of the time) not a smooth one. All the pedagogical (but non-trivial) examples of telescopers displayed in this paper can be viewed by an algebraic geometer as a presentation of examples of families of varieties and their Picard–Fuchs equations.

Appendix B. Maximum Number of Parameters for Families of Planar Elliptic Curves

We have seen in Section 3 that the previous results on the diagonals of 9- or 10-parameter families of rational functions of three variables being pullbacked 2 F 1 hypergeometric functions (and in fact, classical modular forms) can actually be seen as corresponding to the fact (well-known in integrable models and integrable mappings) that the most general biquadratic corresponding to elliptic curves is a 9-parameter family, and that the most general ternary cubic corresponding to elliptic curves is a 10-parameter family. One can, for instance, recall page 238 of [76], which amounts to a consideration of the collection of all cubic curves in C P 2 with the homogeneous equation:
a x 3 + b x 2 y + c x y 2 + d y 3 + e x 2 z + f x z 2 + g y 2 z + h y z 2 + i z 3 + j x y z = 0 ,
and the associated problems of passing through nine given points. One can also recall the ternary cubics in [77,78] and other problems of elliptic curves of high rank [79] (see the concept of Neron–Severi rank).
Since the rational functions of three variables that we consider are essentially encoded by the denominator of these rational functions, and in the cases we have considered, the emergence of pullbacked 2 F 1 hypergeometric functions (and in fact, classical modular forms) corresponds to the fact that the intersection of these denominators with the hyperbola p = x y z corresponds to elliptic curves, one sees that these rational functions are essentially classified by the possible n -parameter families P ( x , y ) = 0 of elliptic curves.
If one considers a polynomial:
P ( x , y ) = m n a m , n · x m y n ,
with generic coefficients a m , n C , then the genus of the algebraic curve defined by P is determined by the support s u p p ( P ) = { ( m , n ) N 2 : a m , n 0 } . More precisely, the genus equals the number of interior integer lattice points inside the convex hull of s u p p ( P ) [80] (see also the discussion in [81]). For example, the support of the 10-parameter family (11) consists of the following 10 points in N 2 :
( 0 , 0 ) , ( 0 , 1 ) , ( 0 , 2 ) , ( 0 , 3 ) , ( 1 , 1 ) , ( 1 , 2 ) , ( 1 , 3 ) , ( 2 , 2 ) , ( 2 , 3 ) , ( 3 , 3 )
which forms a right triangle of side length 3. Only one of these points is an interior point, namely ( 1 , 2 ) ; hence, the genus is 1.
Therefore, we may ask: which integer lattice polytopes exist that have exactly one interior point, and what is the largest such polytope? Not surprisingly, the answer is known: there are (up to transformations like translation, rotation, and shearing) exactly 16 different polytopes with a single interior point [82] (see also Figure 5, page 548 in [83]), the above-mentioned right triangle being the one with the highest total number of lattice points.
This shows that there cannot be a family of elliptic curves with more than 10 parameters.

Appendix C. Monomial Transformations Preserving Pullbacked Hypergeometric Results

More generally, recalling Section 4.2 in [2], and Section 4.2, page 17, in [1], let us consider the monomial transformation:
( x , y , z ) M ( x , y , z ) = ( x M , y M , z M ) = x A 1 · y A 2 · z A 3 , x B 1 · y B 2 · z B 3 , x C 1 · y C 2 · z C 3 ,
where the A i ’s, B i ’s, and C i ’s are positive integers, such that A 1 = A 2 = A 3 is excluded (as well as B 1 = B 2 = B 3 and C 1 = C 2 = C 3 ), and that the determinant (Note a typo in Footnote 28, page 17, of [1], as well as in the second footnote, page 18, in [2]. The sentence has been truncated. One should read: For n = 1 , the 3 × 3 matrix (A10) is stochastic and transformation (A9) is a birational transformation if the determinant of the matrix (A10) is ± 1 ) of the 3 × 3 matrix [1,2]:
A 1 B 1 C 1 A 2 B 2 C 2 A 3 B 3 C 3 ,
is not equal to zero (We want the rational function R ˜ = R ( M ( x , y , z ) ) deduced from the monomial transformation (A9) remains a rational function of three variables, and not of two, or one, variables), and that:
A 1 + B 1 + C 1 = A 2 + B 2 + C 2 = A 3 + B 3 + C 3 .
We will denote by n = A i + B i + C i the integer in these three equal sums (A11). Condition (A11) is introduced in order to impose that the product (Recall that taking the diagonal of a rational function of three variables extracts, in the multi-Taylor expansion, only the terms that are n -th powers of the product x y z ) of x M y M z M is an integer power of the product of x y z : x M y M z M = ( x y z ) n .
If we take a rational function R ( x , y , z ) in three variables and perform such a monomial transformation (A9) ( x , y , z ) M ( x , y , z ) , on this rational function R ( x , y , z ) , we obtain another rational function that we denote by R ˜ = R ( M ( x , y , z ) ) . Now, the diagonal of R ˜ is the diagonal of R ( x , y , z ) where we have changed x into x n :
Φ ( x ) = D i a g R x , y , z , D i a g R ˜ x , y , z = Φ ( x n ) .

Appendix D. Weierstrass and Legendre Forms

The telescoper of the rational function in three variables:
x y ( 1 + y ) 2 x · ( 1 x ) · ( x x y z ) ,
associated (the diagonal extracts the terms function of the product p = x y z in the multi-Taylor series) with the elliptic curve in a Weierstrass form:
( 1 + y ) 2 x · ( 1 x ) · ( x p ) = 0 ,
is the order-two linear differential operator:
L 2 = 1 + 4 · ( 1 2 x ) · D x + 4 · x · ( 1 x ) · D x 2 ,
which has the hypergeometric solution:
2 F 1 1 2 , 1 2 ; 1 ; x = ( 1 x + x 2 ) 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 27 4 · x 2 · ( 1 x ) 2 ( 1 x + x 2 ) 3 .
The elliptic curve (A14) has the Hauptmodul:
H = 27 4 · p 2 · ( 1 p ) 2 ( 1 p + p 2 ) 3 .
in agreement with the pullback in (A16).

Appendix D.1. K3 Surfaces as Products or Foliations of Two Elliptic Curves

The examples of Section 4 correspond to denominators which are algebraic varieties that can be seen as Weierstrass elliptic curves for fixed values of all the variables except two. Let us show other simple telescopers for rational functions with denominators which are algebraic varieties with some foliation in elliptic curves (such as K3 surfaces or Calabi–Yau threefolds).
The telescoper of the rational function in four variables:
x y z ( 1 + z ) 2 x · ( 1 x ) · y · ( x y ) · ( y x y z w ) ,
associated with the K 3 surface written in a Legendre form (along this line, see the first equation, page 19 of [84])
( 1 + z ) 2 x · ( 1 x ) · y · ( x y ) · ( y p ) = 0 ,
is an order-three self-adjoint (the order-three linear differential operator is thus the symmetric square of an order-two linear differential operator) linear differential operator, L 3 :
L 3 = x · ( 2 θ + 1 ) 3 8 · θ 3 ,
which has the following 3 F 2 solution (which is also, because of Clausen’s formula, the square of a 2 F 1 function):
3 F 2 1 2 , 1 2 , 1 2 ; 1 , 1 ; x = 2 F 1 1 4 , 1 4 ; 1 ; x 2 .
The K 3 surface (A19) can be seen as being associated with the product of two Weierstrass elliptic curves ( K 3 surfaces are not abelian varieties, but they are “close” to abelian varieties: from a creative telescoping viewpoint, they essentially can be seen as products of two elliptic curves) of Hauptmoduls, respectively:
H x = 27 4 · p 2 · ( 1 p ) 2 ( 1 p + p 2 ) 3 , H y = 27 4 · y 2 · ( 1 y ) 2 ( 1 y + y 2 ) 3 .
This order-three linear differential operator L 3 is the symmetric square of the order-two linear differential operator:
M 2 = 1 + 8 · ( 2 3 x ) · D x + 16 · x · ( 1 x ) · D x 2 ,
which has the hypergeometric solutions:
2 F 1 1 4 , 1 4 ; 1 ; x = 1 x 4 1 / 4 · 2 F 1 1 12 , 5 12 ; 1 ; 27 · x 2 ( x 4 ) 3 .

Appendix D.2. Calabi–Yau Threefolds as Foliation in Three Elliptic Curves

The telescoper of the rational function in five variables x , y , z , v , and w :
x y z v ( 1 + w ) 2 x · ( 1 x ) · y · ( x y ) · z · ( y z ) · ( z x y z v w ) ,
associated (the diagonal extracts the terms function of the product p = x y z v w in the multi-Taylor series) with the Calabi–Yau threefold written in a Legendre form:
( 1 + w ) 2 x · ( 1 x ) · y · ( x y ) · z · ( y z ) · ( z p ) = 0 ,
is an order-four (self-adjoint) linear differential operator L 4 :
L 4 = 16 · θ 4 x · ( 2 θ + 1 ) 4 ,
which is a Calabi–Yau operator (This linear differential operator is self-adjoint, its exterior square is of order five, and it is MUM (maximum unipotent monodromy [27,65,66])) with the 4 F 3 solution:
4 F 3 1 2 , 1 2 , 1 2 , 1 2 ; 1 , 1 , 1 ; x .
For y and z fixed, the Calabi–Yau threefold (A26) is foliated in genus-one curves:
( 1 + w ) 2 λ · x · ( 1 x ) · ( x y ) = 0 ,
where λ is the constant expression (p is fixed):
λ = y · z · ( y z ) · ( z p ) .
The Hauptmodul of these genus-one curves is independent of p and z , reading:
H y , z = 27 4 · y 2 · ( 1 y ) 2 ( 1 y + y 2 ) 3 .
Similarly for x and z fixed, the Calabi–Yau threefold (A26) is foliated in genus-one curves:
( 1 + w ) 2 μ · y · ( x y ) · ( y z ) = 0 ,
where μ is the constant expression (p is fixed):
μ = x · z · ( 1 x ) · ( z p ) .
The genus-one curves (A32) can be written in a simpler Weierstrass form:
( 1 + w ) 2 ρ · Y · 1 Y · Y z x = 0 ,
where the constant ρ reads ρ = μ · x 3 , and the variable y has been rescaled into Y = y / x . The Hauptmodul of these genus-one curves (A32) is the same as the Hauptmodul of the genus-one curves (A29), and corresponds to expression (A31) where y has been changed into z / x (see the canonical form (A34)), namely:
H x , z = 27 4 · x 2 · z 2 · ( x z ) 2 ( x 2 x z + z 2 ) 3 .
Similarly for x and y fixed, the Calabi–Yau threefold (A26) is foliated in genus-one curves,
( 1 + w ) 2 ν · z · ( y z ) · ( z p ) = 0 ,
where ν reads:
ν = x · ( 1 x ) · y · ( x y ) .
A reduction to a canonical Weierstrass form similar to (A34) gives immediately the Hauptmodul of the genus-one curve (A36), which reads:
H x , y = 27 4 · y 2 · p 2 · ( y p ) 2 ( y 2 y p + p 2 ) 3 .
The Calabi–Yau threefold (A26) thus has a foliation in a triple of elliptic curves E 1 , E 2 , and E 3 .

Appendix E. Rational Functions with Tri-Quadratic and N-Quadratic Denominators

Appendix E.1. Rational Functions with Tri-Quadratic Denominators

Let us consider the most general tri-quadratic surface:
m = 0 , 1 , 2 n = 0 , 1 , 2 l = 0 , 1 , 2 a m , n , l · x m y n z l = 0 ,
depending on 27 = 3 3 parameters a m , n , l . It can be rewritten as:
A ( x , y ) · z 2 + B ( x , y ) · z + C ( x , y ) = 0 .
It is straightforward to see that condition (A40) is preserved by the birational involution  I z :
I z : x , y , z x , y , C ( x , y ) A ( x , y ) · 1 z ,
and we have, of course, two other similar birational involutions, I x and I y , that single out x and y , respectively. The (generically) infinite-order birational transformations K x = I y · I z , K y = I z · I x , and K z = I x · I y are birational symmetries of the surface (A39) or (A40). They are related by K x · K y · K z = i d e n t i t y . Note that the birational transformation K x preserves x . The iteration of the (generically) infinite-order birational transformation K x gives elliptic curves. Since Equations (A39) or (A40) is preserved by K x , which also preserves x , the equation of the elliptic curves corresponding to the iteration (The birational transformation K x maps the elliptic curve onto itself (self-map). One can use the iteration of the birational transformation K x to actually visualize the elliptic curve [44,85]) of K x is actually (A39) for fixed values of x . Equation (A39), for fixed values of x , is a (general) biquadratic curve in y and z , and is thus an elliptic curve depending on x . Therefore, one has a canonical foliation of the algebraic surface (A39) in elliptic curves. Of course, the iteration of K y (resp. K z ) also yields elliptic curves, and similarly yields two other foliations in elliptic curves.
We have a foliation in two families of elliptic curves E and E of the surface. Consequently, this tri-quadratic surface (A39), having an infinite set of birational automorphisms and an infinite set of birational symmetries, cannot be of the “general type” (it has Kodaira dimension of less than 2 ).

Appendix E.2. Rational Functions with N-Quadratic Denominators

The calculations of Appendix E.1 can straightforwardly be generalized to N -quadratic equations, writing the N -quadratic (72) as:
A ( x 1 , x 2 , , x N 1 ) · x N 2 + B ( x 1 , x 2 , , x N 1 ) · x N + C ( x 1 , x 2 , , x N 1 ) = 0 ,
and introducing the birational involution I N :
I N : x 1 , x 2 , , x N x 1 , x 2 , , x N 1 , C ( x 1 , x 2 , , x N 1 ) A ( x 1 , x 2 , x N 1 ) · 1 x N .
Similarly to Appendix E.1, we can introduce N involutive birational transformations I m and consider the products of two such involutive birational transformations K m , n = I m · I n . These K m , n ’s are (generically) infinite-order birational transformations preserving the N 2 variables that are not x m and x n .
Using such remarkable N -variable algebraic varieties, with an infinite set of birational automorphisms, one can build rational functions of N + 1 variables, any of the parameter of the algebraic variety, becoming an arbitrary rational (or even an arbitrary algebraic) function of the product p = x 1 x 2 x N , with a Taylor series expansion at p = 0 , the diagonal of rational functions becoming diagonal of algebraic functions of the product p = x 1 x 2 x N in order to build the denominator of the rational function. The telescopers of such rational functions are seen (experimentally using creative telescoping) to be of substantially smaller order than the ones for rational functions whose denominators are, after reduction by p = x 1 x 2 x N , associated with algebraic varieties of “general type”.

Appendix F. Telescopers of Rational Functions of Several Variables: Some Examples

Let us consider here the following family of rational functions in four variables:
R ( x , y , z , u ) = 1 1 + 3 y + z + 9 y z + 11 z 2 y + 3 u x + x · P ( x , y , z ) ,
where P ( x , y , z ) is an arbitrary polynomial of the three variables x , y , and z .

Appendix F.1. Telescopers of Rational Functions of Several Variables: A Second Example with Four Variables

Let us now consider the rational function in four variables x , y , z , u :
R ( x , y , z , u ) = 1 1 + 3 y + z + 9 y z + 11 z 2 y + 3 u x + 9 x + 2 x y + 5 x z + 7 x 2 y .
which corresponds to P ( x , y , z ) = 9 + 2 y + 5 z + 7 x y . The telescoper of this rational function of four variables is the same order-two linear differential operator L 2 as for the telescoper of (91). It has the same pullbacked hypergeometric solution (92). The diagonal of the rational function (A45) is the expansion of (92), namely (93).
Performing the intersection of the codimension-one algebraic variety:
1 + 3 y + z + 9 y z + 11 z 2 y + 3 u x + 9 x + 2 x y + 5 x z + 7 x 2 y = 0 ,
corresponding to the denominator of (A45), with the hyperbola p = x y z u , amounts to eliminating, for instance, u (writing u = p x y z ). This gives P u = 0 , where P u reads:
P u = 7 x 2 y 2 z + 2 x y 2 z + 5 x y z 2 + 9 x y z + 11 y 2 z 3 + 9 y 2 z 2 + 3 y 2 z + y z 2 + y z + 3 p .
Assuming x to be constant, the previous condition P u ( y , z ) = 0 is an algebraic curve. Calculating its genus, one finds immediately that it has genus one. Calculating its j -invariant, one deduces the expression of the Hauptmodul H p , x = 1728 J as a rational expression of p and x :
H p , x = 1728 J = 46656 p 3 · ( 7 x 2 + 2 x + 3 ) 2 · N D 3 ,
where N is a polynomial expression of degree eight in w and three in p , and D is a polynomial expression of degree four in w and two in p . In the x 0 limit of the Hauptmodul H p , x = 1728 J , one finds:
H p = 419904 · p 3 · ( 5 12 p 19440 p 2 + 2665872 p 3 ) ( 1 2592 p 2 ) 3 ,
which is actually the Hauptmodul in (92). In other words, the exact expression of the diagonal of the rational function (A45), which is (92), and is essentially encapsulated in the Hauptmodul in (92), could have been obtained from the x = 0 selection of the Hauptmoduls H p , x .

Appendix F.2. Telescopers of Rational Functions of Several Variables: A Third Example with Four Variables

Let us consider the rational function in four variables x , y , z , u :
R ( x , y , z , u ) = 1 1 + 3 y + z + 9 y z + 11 z 2 y + 3 u x + x · ( y 2 z 2 + x y 3 ) ,
which corresponds to P ( x , y , z ) = y 2 z 2 + x y 3 in the family (A44). Again, the telescoper of this rational function of four variables is the same order-two linear differential operator L 2 as for the telescoper of (91). It has the same pullbacked hypergeometric solution (92). Actually, the diagonal of the rational function (91) is the expansion (93) of the pullbacked hypergeometric function (92). In this case (A49), the elimination of u = p x y z in the vanishing condition of the denominator (A49) gives the algebraic curve:
x 2 y 4 z + x y 3 z 3 + 11 y 2 z 3 + 9 y 2 z 2 + 3 y 2 z + y z 2 + y z + 3 p = 0 .
For x fixed (and of course, p fixed) this algebraic curve (A50) is a genus-five curve, but, of course, in the x = 0 case, it reduces to the same genus-one curve as for the first example (91), namely:
11 y 2 z 3 + 9 y 2 z 2 + 3 y 2 z + y z 2 + y z + 3 p = 0 .
which corresponds to the Hauptmodul (A48).
The generalization of this result is straightforward. Let us consider the rational function in four variables x , y , z , and u :
R ( x , y , z , u ) = 1 1 + 3 y + z + 9 y z + 11 z 2 y + 3 u x + x · P ( x , y , z ) ,
where P ( x , y , z ) is an arbitrary polynomial of the three variables x , y , and z . On a large set of examples, one verifies that the diagonal of (A52) is actually the expansion (93) of the pullbacked hypergeometric function (92):
1 + 648 x 2 72900 x 3 + 1224720 x 4 330674400 x 5 + 23370413220 x 6 1276733858400 x 7 + 180019474034400 x 8 12013427240614800 x 9 +
However, as far as creative telescoping calculations are concerned (using the HolonomicFunctions package [3]), the telescoper corresponding to different polynomials P ( x , y , z ) quickly becomes a quite large non-minimal linear differential operator. For instance, even for the simple polynomial P ( x , y , z ) = x + y , one obtains a quite large order-10 telescoper. Of course, since this telescoper has the pullbacked hypergeometric function (92) as a solution, it is not minimal; it is right-divisible by the order-two linear differential operator having (92) as a solution. It is straightforward to see that the previous elimination of u = p x y z in the vanishing condition of the denominator (A52) gives an algebraic curve (of arbitrary large genus for increasing degrees of the polynomial P ( x , y , z ) ):
11 y 2 z 3 + 9 y 2 z 2 + 3 y 2 z + y z 2 + y z + 3 p + y z · P ( x , y , z ) = 0 .
which reduces again, in the x = 0 case, to the same genus-one curve (A51).
With that general example (A52), we see that there is an infinite set of rational functions depending on an arbitrary polynomial P ( x , y , z ) of three variables whose diagonals are actually a pullbacked 2 F 1 hypergeometric solution, namely (92).

References

  1. Abdelaziz, Y.; Boukraa, S.; Koutschan, C.; Maillard, J.-M. Diagonals of rational functions, pullbacked 2F1 hypergeometric functions and modular forms. J. Phys. 2018, 51, 455201. [Google Scholar]
  2. Abdelaziz, Y.; Boukraa, S.; Koutschan, C.; Maillard, J.-M. Diagonals of rational functions, pullbacked 2F1 hypergeometric functions and modular forms (unabridged version). arXiv 2018, arXiv:1805.04711v1. [Google Scholar]
  3. HolonomicFunctions Package Version 1.7.1 (9-October-2013) written by Christoph Koutschan; Copyright 2007–2013; Research Institute for Symbolic Computation (RISC), Johannes Kepler University: Linz, Austria, 2018.
  4. Christol, G. Diagonales de fractions rationnelles et équations différentielles. In Study Group on Ultrametric Analysis, 10th Year: 1982/83, No. 2, Exp. No. 18; Institute Henri Poincaré: Paris, France, 1984; pp. 1–10. Available online: http://archive.numdam.org/article/GAU_1982-1983__10_2_A4_0.pdf (accessed on 10 June 2022).
  5. Christol, G. Diagonales de Fractions Rationnelles et Équations de Picard-Fuchs. In Groupe de Travail D’analyse Ultramétrique, Tome 12 (1984–1985) no. 1, Exposé no. 13, 12p. 1984. Available online: http://archive.numdam.org/article/GAU_1984-1985__12_1_A8_0.pdf (accessed on 10 June 2022).
  6. Christol, G. Diagonales de fractions rationnelles. In Séminaire de Théorie des Nombres, Paris 1986–87 (Progr. Math. vol. 75); Birkhäuser Boston: Boston, MA, USA, 1988; pp. 65–90. [Google Scholar]
  7. Christol, G. Globally Bounded Solutions of Differential Equations, Analytic Number Theory (Tokyo, 1988) (Lecture Notes in Math. 1434); Springer: Berlin, Germany, 1990; pp. 45–64. [Google Scholar] [CrossRef]
  8. Lipshitz, L.; van der Poorten, A.J. Rational functions, diagonals, automata and arithmetic. In Number Theory (Banff, AB, 1988); de Gruyter: Berlin, Germany, 1990; pp. 339–358. [Google Scholar]
  9. Lipshitz, L. The diagonal of a D-finite power series is D-finite. J. Algebra 1988, 113, 373–378. [Google Scholar] [CrossRef] [Green Version]
  10. Denef, J.; Lipshitz, L. Algebraic power series and diagonals. J. Number Theory 1987, 26, 46–67. [Google Scholar] [CrossRef] [Green Version]
  11. Bostan, A.; Boukraa, S.; Maillard, J.-M.; Weil, J.-A. Diagonal of rational functions and selected differential Galois groups. J. Phys. A Math. Theor. 2015, 48, 504001. [Google Scholar] [CrossRef] [Green Version]
  12. Maier, R.S. On rationally parametrized modular equations. J. Ramanujan Math. Soc. 2009, 24, 1–73. Available online: http://arxiv.org/abs/math/0611041 (accessed on 10 June 2022).
  13. Glasser, M.L.; Guttmann, A.J. Lattice Green function (at 0) for the 4D hypercubic lattice. J. Phys. A 1994, 27, 7011–7014. Available online: http://arxiv.org/abs/cond-mat/9408097 (accessed on 10 June 2022). [CrossRef]
  14. Guttmann, A.J. Lattice Green’s functions in all dimensions. J. Phys. A 2010, 43, 305205. [Google Scholar] [CrossRef] [Green Version]
  15. Zenine, N.; Hassani, S.; Maillard, J.-M. Lattice Green Functions: The seven-dimensional face-centred cubic lattice. J. Phys. A Math. Theor. 2015, 48, 035205. [Google Scholar] [CrossRef] [Green Version]
  16. Hassani, S.; Koutschan, C.; Maillard, J.-M.; Zenine, N. Lattice Green Functions: The d-dimensional face-centred cubic lattice, d=8,9,10,11,12. J. Phys. A Math. Theor. 2016, 49, 164003. [Google Scholar] [CrossRef] [Green Version]
  17. Abdelaziz, Y.; Boukraa, S.; Koutschan, C.; Maillard, J.-M. Heun functions and diagonals of rational functions. J. Phys. A Math. Theor. 2020, 53, 075206. [Google Scholar] [CrossRef] [Green Version]
  18. Abdelaziz, Y.; Boukraa, S.; Koutschan, C.; Maillard, J.-M. Heun functions and diagonals of rational functions (unabridged version). arXiv 2019, arXiv:1910.10761. [Google Scholar]
  19. Bostan, A.; Boukraa, S.; Christol, G.; Hassani, S.; Maillard, J.-M. Ising n-fold integrals as diagonals of rational functions and integrality of series expansions. J. Phys. A Math. Theor. 2013, 46, 185202. [Google Scholar] [CrossRef] [Green Version]
  20. Bostan, A.; Boukraa, S.; Christol, G.; Hassani, S.; Maillard, J.-M. Ising n-fold integrals as diagonals of rational functions and integrality of series expansions: Integrality versus modularity Preprint. arXiv 2012, arXiv:1211.6031. [Google Scholar]
  21. Assis, M.; Boukraa, S.; Hassani, S.; van Hoeij, M.; Maillard, J.-M.; McCoy, B.M. Diagonal Ising susceptibility: Elliptic integrals, modular forms and Calabi-Yau equations. J. Phys. A 2012, 45, 075205. [Google Scholar] [CrossRef] [Green Version]
  22. Boukraa, S.; Hassani, S.; Maillard, J.-M.; Weil, J.-A. Differential algebra on lattice Green functions and Calabi-Yau operators. J. Phys. A Math. Theor. 2014, 48, 095203. [Google Scholar] [CrossRef]
  23. Bostan, A.; Boukraa, S.; Guttmann, A.J.; Hassani, S.; Jensen, I.; Maillard, J.-M.; Zenine, N. High order Fuchsian equations for the square lattice Ising model: χ˜(5). J. Phys. A Math. Theor. 2009, 42, 275209. [Google Scholar] [CrossRef]
  24. Boukraa, S.; Hassani, S.; Jensen, I.; Maillard, J.-M.; Zenine, N. High-order Fuchsian equations for the square lattice Ising model: χ(6). J. Phys. A Math. Theor. 2010, 43, 115201. [Google Scholar] [CrossRef] [Green Version]
  25. Boukraa, S.; Hassani, S.; Maillard, J.-M.; Zenine, N. Singularities of n-fold integrals of the Ising class and the theory of elliptic curves. J. Phys. A Math. Theor. 2007, 40, 11713–11748. [Google Scholar] [CrossRef] [Green Version]
  26. Assis, M.; van Hoeij, M.; Maillard, J.-M. The perimeter generating functions of three-choice, imperfect, and one-punctured staircase polygons. J. Phys. A Math. Theor. 2016, 49, 214002. [Google Scholar] [CrossRef] [Green Version]
  27. Van Straten, D. Calabi-Yau operators. arXiv 2017, arXiv:1704.00164v1. [Google Scholar]
  28. Bostan, A.; Lairez, P.; Salvy, B. Creative telescoping for rational functions using the Griffiths-Dwork method. In Proceedings of the 38th International Symposium on Symbolic and Algebraic Computation, Boston, MA, USA, 26–29 June 2013; pp. 93–100. Available online: http://specfun.inria.fr/bostan/publications/BoLaSa13.pdf (accessed on 10 June 2022).
  29. Chen, S.; Kauers, M.; Singer, M.F. Telescopers for Rational and Algebraic Functions via Residues. In Proceedings of the 37th International Symposium on Symbolic and Algebraic Computation, Grenoble, France, 22–25 July 2012; pp. 130–137. [Google Scholar]
  30. Takeuchi, K. Commensurability classes of arithmetic triangle groups. J. Fac. Science Univ. Tokyo Sec. IA Math. 1977, 24, 201–212. [Google Scholar]
  31. Voight, J. Shimura curves of genus at most two. Math. Comp. 2009, 78, 1155–1172. [Google Scholar] [CrossRef] [Green Version]
  32. Doran, C.F.; Malmendier, A. Calabi-Yau Manifols Realizing Symplectically Rigid Monodromy Tuples. Available online: https://arxiv.org/pdf/1503.07500.pdf (accessed on 10 June 2022).
  33. Malmendier, A.; Shaska, T. (Eds.) Higher Genus Curves in Mathematical Physics and Arithmetic Geometry, Contemporary Mathematics, AMS Special Session Higher Genus Curves and Fibrations in Mathematical Physics and Arithmetic Geometry. January 2016. Available online: https://www.amazon.com/Higher-Mathematical-Physics-Arithmetic-Geometry/dp/1470428563 (accessed on 10 June 2022).
  34. Chyzak, F. The ABC of Creative Telescoping-Algorithms, Bounds, Complexity. 2014. Available online: https://tel.archives-ouvertes.fr/tel-01069831/document (accessed on 10 June 2022).
  35. Lairez, P. Périodes d’intégrales Rationnelles: Algorithmes et Applications, Thèse de Doctorat. Available online: https://pierre.lairez.fr/these.pdf (accessed on 10 June 2022).
  36. Zeilberger, D. The Method of Creative Telescoping. J. Symb. Comput. 1991, 11, 195–204. [Google Scholar] [CrossRef] [Green Version]
  37. Kontsevich, M.; Zagier, D. Periods, IHES/M/01/22. 2001. Available online: https://www.maths.ed.ac.uk/~v1ranick/papers/kontzagi.pdf (accessed on 10 June 2022).
  38. Igusa, J.I. Abstract vanishing cycle theory. Proc. Jpn. Acad. 1958, 34, 589–593. [Google Scholar] [CrossRef]
  39. Deligne, P. Intégration sur un cycle évanescent. In Inventiones Math; Springer: New York, NY, USA, 1983; Volume 76, pp. 1-29–1-43. [Google Scholar]
  40. Lairez, P. Computing periods of rational integrals. Math. Comput. 2016, 85, 1719–1752. [Google Scholar] [CrossRef] [Green Version]
  41. Chyzak, F. An extension of Zeilberger’s fast algorithm to general holonomic functions. Discret. Math. 2000, 217, 1–3. [Google Scholar] [CrossRef] [Green Version]
  42. Koutschan, C. A fast approach to creative telescoping. Math. Comput. Sci. 2010, 4, 259–266. [Google Scholar] [CrossRef] [Green Version]
  43. Boukraa, S.; Maillard, J.-M. Symmetries of lattice models in statistical mechanics and effective algebraic geometry. J. Phys. I 1993, 3, 239–258. [Google Scholar]
  44. Bellon, M.P.; Maillard, J.-M.; Viallet, C.-M. Quasi integrability of the sixteen-vertex model. Phys. Lett. 1992, 281, 315–319. [Google Scholar] [CrossRef]
  45. Boukraa, S.; Maillard, J.-M.; Rollet, G. Determinantal identities on integrable mappings. Int. J. Mod. Phys. 1994, 8, 2157–2201. [Google Scholar] [CrossRef]
  46. Bronstein, M.; Mudders, T.; Weil, J.-A. On Symmetric Powers of Differential Operators. In Proceedings of the 1997 International Symposium on Symbolic and Algebraic Computation, Maui, HI, USA, 21–23 July 1997. [Google Scholar]
  47. Long, L. On Shioda-Inose structures of one-parameter families of K3 surfaces. J. Number Theory 2004, 109, 299–318. [Google Scholar] [CrossRef] [Green Version]
  48. Kumar, A. Hilbert Modular Surfaces for square discriminants and elliptic subfields of genus 2 function fields. Res. Math. Sci. 2015, 2, 1–46. [Google Scholar] [CrossRef] [Green Version]
  49. Kumar, A. Elliptic Fibrations on a Generic Jacobian Kummer Surface. arXiv 2014, arXiv:1105.1715v3. [Google Scholar] [CrossRef] [Green Version]
  50. Kumar, A.; Mukamel, R. Algebraic Models and Arithmetic Geometry of Teichmüller Curves in Genus Two. Available online: https://arxiv.org/pdf/1406.7057.pdf (accessed on 10 June 2022).
  51. Elkies, N.; Kumar, A. K3 surfaces and equations for Hilbert modular surfaces. Algebra Number Theory 2014, 8, 2297–2411. [Google Scholar] [CrossRef] [Green Version]
  52. Kuhn, R.M. Curves of genus 2 with split Jacobian. Trans. Am. Math. Soc. 1988, 307, 41–49. [Google Scholar] [CrossRef] [Green Version]
  53. Shaska, T. Genus 2 curves with (3,3)-split Jacobian and large automorphism group. arXiv 2002, arXiv:0201008v1. [Google Scholar]
  54. Shaska, T. Genus 2 fields with degree 3 elliptic subfields. Forum Math. 2004, 16, 263–280. [Google Scholar] [CrossRef] [Green Version]
  55. Shaska, T.; Völklein, H. Elliptic subfields and automorphisms of genus 2 function fields. arXiv 2001, arXiv:0107142v1. [Google Scholar]
  56. Bedford, E.; Kim, K.; Truong, T.T.; Abarenkova, N.; Maillard, J.-M. Degree Complexity of a Family of Birational Maps. In Mathematical Physics, Analysis and Geometry; Springer: New York, NY, USA, 2007. [Google Scholar]
  57. Anglès d’Auriac, J.-C.; Maillard, J.-M.; Viallet, C.M. On the complexity of some birational transformations. J. Phys. Math. Gen. 2006, 39, 3641–3654. [Google Scholar] [CrossRef] [Green Version]
  58. Maillard, J.-M. Automorphisms of algebraic varieties and Yang-Baxter equations. J. Math. Phys. 1986, 27, 2776. [Google Scholar] [CrossRef]
  59. Bellon, M.P.; Maillard, J.-M.; Viallet, C.-M. Infinite discrete symmetry group for the Yang-Baxter equations: Spin models. Phys. Lett. 1991, 157, 343–353. [Google Scholar] [CrossRef]
  60. Bellon, M.P.; Maillard, J.-M.; Viallet, C.-M. Infinite discrete symmetry group for the Yang-Baxter equations: Vertex models. Phys. Lett. B 1991, 260, 87–100. [Google Scholar] [CrossRef]
  61. Corti, A. Polynomial bounds for the number of automorphisms of a surface of general type. Ann. Sci. Ecole Norm. Sup. 1991, 24, 113–137. [Google Scholar] [CrossRef] [Green Version]
  62. Szabó, E. Bounding automorphism groups. Math. Ann. 1996, 304, 801–811. [Google Scholar] [CrossRef]
  63. Hacon, C.D.; McKernan, J.; Xu, C. On the birational automorphisms of varieties of general type. Ann. Math. 2013, 177, 1077–1111. [Google Scholar] [CrossRef] [Green Version]
  64. Abdelaziz, Y.; Maillard, J.-M. Modular forms, Schwarzian conditions, and symmetries of differential equations in physics. J. Phys. A Math. Theor. 2017, 50, 215203. [Google Scholar] [CrossRef] [Green Version]
  65. Almkvist, G.; van Enckevort, C.; van Straten, D.; Zudilin, W. Tables of Calabi-Yau equations. arXiv 2010, arXiv:0507430v2. [Google Scholar]
  66. Bostan, A.; Boukraa, S.; Hassani, S.; van Hoeij, M.; Maillard, J.-M.; Weil, J.-A.; Zenine, N.J. The Ising model: From elliptic curves to modular forms and Calabi-Yau equations. J. Phys. A Math. Theor. 2011, 44, 045204. [Google Scholar] [CrossRef] [Green Version]
  67. Singer, M.F. Solving homogeneous linear differential equations in terms of second order linear differential equations. Am. J. Math. 1985, 107, 663–696. [Google Scholar] [CrossRef]
  68. Person, A.C. Solving Homogeneous Linear Differential Equations of Order 4 in Terms of Equations of Smaller Order. Ph.D. Thesis, Raleigh, North Carolina, 2002. Available online: http://www.lib.ncsu.edu/resolver/1840.16/3059 (accessed on 10 June 2022).
  69. Van Hoeij, M. Solving third order linear differential equations in terms of second order equations. In Proceedings of the 2007 International Symposium on Symbolic and Algebraic Computation, ACM, Waterloo, ON, Canada, 28 July–1 August 2007. [Google Scholar]
  70. Fricke, R.; Klein, F. Vorlesungen über die Theorie der automorphen Funktionen. I; Druck und Verlag von B. G. Teubner: Leipzig, Germany, 1897; p. 366. [Google Scholar]
  71. Boalch, P.; Paluba, R. Symmetric cubic surfaces and G2 character varieties. arXiv 2013, arXiv:1305.6594v2. [Google Scholar]
  72. Cantat, S.; Loray, F. Holomorphic Dynamics, Painlevé VI Equation and Character Varieties. Available online: https://hal.archives-ouvertes.fr/hal-00186558v2 (accessed on 22 May 2022).
  73. Mazzocco, M.; Vidunas, R. Cubic and Quartic Transformations of the Sixth Painlevé Equation in Terms of Riemann-Hilbert Correspondence. Stud. Appl. Math. 2013, 130, 17–48. [Google Scholar] [CrossRef] [Green Version]
  74. Picard, E. Sur les intégrales doubles de fonctions rationnelles dont les résidus sont nuls. Bulletin des Sciences Mathématiques Série 1902, 2, 26. [Google Scholar]
  75. Griffiths, P.A. On the periods of certain rational integrals I, II. Ann. Math. 1969, 90, 460–541. [Google Scholar] [CrossRef]
  76. Silverman, J.H.; Tate, J. Rational Points on elliptic Curves. In Undergraduate Texts in Mathematics; Springer: New York, NY, USA, 1992. [Google Scholar]
  77. Sadek, M. Minimal genus one curves. Funct. Approx. 2012, 46, 117–131. [Google Scholar] [CrossRef]
  78. Poonen, B. An explicit algebraic family of genus-one curves violating the Hasse principle. arXiv 2001, arXiv:9910124v1. [Google Scholar] [CrossRef]
  79. Elkies, N. Three lectures on elliptic surfaces and curves of high rank. arXiv 2007, arXiv:0709.2908v1. [Google Scholar]
  80. Khovanskii, A.G. Newton polyhedra and the genus of complete intersections. Funct. Anal. i Ego Pril. English Transl. Funct. Anal. Appl. 1978, 12, 38–46. [Google Scholar] [CrossRef]
  81. Available online: https://mathoverflow.net/questions/16615/calculating-the-genus-of-a-curve-using-the-newton-polygon (accessed on 10 June 2022).
  82. Rabinowitz, S. A census of convex lattice polygons with at most one interior lattice point. Ars Comb. 1989, 28, 83–96. [Google Scholar]
  83. Schicho, J. Simplification of surface parametrizations–a lattice polygon approach. J. Symb. Comput. 2003, 36, 535–554. [Google Scholar] [CrossRef] [Green Version]
  84. Algreen, S. The Point of a Certain Fivefold over Finite Fields and the Twelfth Power of the eta Function. Finite Fields Their Appl. 2002, 8, 18–33. [Google Scholar] [CrossRef] [Green Version]
  85. Boukraa, S.; Hassani, S.; Maillard, J.-M. Noetherian mappings. Physica 2003, 185, 3–44. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdelaziz, Y.; Boukraa, S.; Koutschan, C.; Maillard, J.-M. Diagonals of Rational Functions: From Differential Algebra to Effective Algebraic Geometry. Symmetry 2022, 14, 1297. https://doi.org/10.3390/sym14071297

AMA Style

Abdelaziz Y, Boukraa S, Koutschan C, Maillard J-M. Diagonals of Rational Functions: From Differential Algebra to Effective Algebraic Geometry. Symmetry. 2022; 14(7):1297. https://doi.org/10.3390/sym14071297

Chicago/Turabian Style

Abdelaziz, Youssef, Salah Boukraa, Christoph Koutschan, and Jean-Marie Maillard. 2022. "Diagonals of Rational Functions: From Differential Algebra to Effective Algebraic Geometry" Symmetry 14, no. 7: 1297. https://doi.org/10.3390/sym14071297

APA Style

Abdelaziz, Y., Boukraa, S., Koutschan, C., & Maillard, J. -M. (2022). Diagonals of Rational Functions: From Differential Algebra to Effective Algebraic Geometry. Symmetry, 14(7), 1297. https://doi.org/10.3390/sym14071297

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop