# Renormalization, Isogenies, and Rational Symmetries

Authors: A. Bostan; S. Boukraa; S. Hassani; J.-M. Maillard; J.-A. Weil; N. Zenine; N. Abarenkova
Journal: Advances in Mathematical Physics
 (2010)

---

## Abstract

We give an example of infinite-order rational transformation that
leaves a linear differential equation covariant. This example can be seen as a nontrivial
but still simple illustration of an exact representation of the renormalization
group.


---

## Body

## 1. Introduction

There is no need to underline the success of the renormalization group revisited by Wilson [1, 2] which is nowadays seen as a fundamental symmetry in lattice statistical mechanics or field theory. It contributed to promote 2d conformal field theories and/or scaling limits of second-order phase transition in lattice statistical mechanics.1 If one does not take into account most of the subtleties of the renormalization group, the simplest sketch of the renormalization group corresponds to Migdal-Kadanoff decimation calculations, where the new coupling constants created at each step of the (real-space) decimation calculations are forced2   to stay in some (slightly arbitrary) finite-dimensional parameter space. This drastic projection may be justified by the hope that the basin of attraction of the fixed points of the corresponding (renormalization) transformation in the parameter space is “large enough.”One heuristic example is always given because it is one of the very few examples ofexact renormalization, the renormalization of the one-dimensional Ising model without a magnetic field. It is a straightforward undergraduate exercise to show that performing various decimations summing over every two, three, orN spins, one gets exact generators of the renormalization group reading T
N
:
t
→
t
N, where t is (with standard notations) the high temperature variable t
=
tanh
(
K
). It is easy to see that these transformations T
N, depending on the integer N, commute together. Such an exact symmetry is associated with a covariance of the partition function per site Z
(
t
)
=
C
(
t
)
·
Z
(
t
2
). In this particular case one recovers the (very simple) expression of the partition function per site, 2
cosh
(
K
), as an infinite product of the action of (for instance) T
2 on the cofactor C
(
t
). In this very simple case, this corresponds to the using of the identity (valid for |
x
|
<
1):(1.1)
∏
n
=
0
∞
(
1
+
x
2
n
)
=
1
1
-
x
.
For T
3
:
t
→
t
3 one must use the identity(1.2)
∏
n
=
0
∞
(
1
+
x
3
n
+
x
2
·
3
n
)
=
∏
n
=
0
∞
(
1
-
x
3
n
+
1
1
-
x
3
n
)
=
1
1
-
x
,
and for T
N
:
t
→
t
N a similar identity where the 3 in the exponents is changed into N.Another simple heuristic example is the one-dimensional Ising modelwith a magnetic field. Straightforward calculations enable to get an infinite number of exact generators of the corresponding renormalization group, represented as rational transformations3
(1.3)
T
N
:
(
x
,
z
)
→
T
N
(
x
,
z
)
=
(
x
N
,
z
N
)
,
where the first two transformations T
2 and T
3 read in terms of the two (low-temperature well-suited and fugacity-like) variables x
=
e
4
K and z
=
e
2
H: (1.4)
x
2
=
(
x
+
z
)
(
1
+
x
z
)
x
·
(
1
+
z
)
2
,
z
2
=
z
·
(
1
+
x
z
)
x
+
z
,
x
3
=
x
·
(
z
2
x
+
2
z
+
1
)
(
z
2
+
2
z
+
x
)
(
z
2
x
+
z
+
x
z
+
x
)
2
,
z
3
=
z
·
z
2
x
+
2
z
+
1
z
2
+
2
z
+
x
.
One simply verifies that these rational transformations of two (complex) variables commute. This can be checked by formal calculations for T
N and T
M for any N and M less than 30, and one can easily verify a fundamental property expected for renormalization group generators:(1.5)
T
N
·
T
M
=
T
M
·
T
N
=
T
N
M
,
where the “dot” denotes the composition of two transformations. The infinite number of these rational transformations of two (complex) variables (1.3) are thus a rational representation of the positive integers together with their product. Such rational transformations can be studied “per se” as discrete dynamical systems, the iteration of any of these various exact generators corresponding to an orbit of the renormalization group.Of course these two examples of exact representation of the renormalization group are extremely degenerate since they correspond to one-dimensional models.4 Migdal-Kadanoff decimation will quite systematically yieldrational
5 transformations similar to (1.3) in two, or more, variables.6 Consequently, they are never (except “academical” self-similar models) exact representations of the renormalization group. The purpose of this paper is to provide simple (but nontrivial) examples of exact renormalization transformations that are not degenerate like the previous transformations on one-dimensional models.7 In several papers [3, 4] for Yang-Baxter integrable models with a canonical genus one parametrization [5, 6] (elliptic functions of modulus k), we underlined that the exact generators of the renormalization group must necessarily identify with the various isogenies which amount to multiplying or dividing τ, the ratio of the two periods of the elliptic curves, by an integer. The simplest example is the Landen transformation [4] which corresponds to multiplying (or dividing because of the modular group symmetry τ
↔
1
/
τ), the ratio of the two periods is(1.6)
k
→
k
L
=
2
k
1
+
k
,
τ
⟷
2
τ
.
The other transformations8 correspond to τ
↔
N
·
τ, for various integers N. In the (transcendental) variable τ, it is clear that they satisfy relations like (1.5). However, in the natural variables of the model (e
K
,
tanh
(
K
)
,
k
=
sinh
(
2
K
), not transcendental variables like τ), these transformations are algebraic transformations corresponding in fact to the fundamental modular curves. For instance, (1.6) corresponds to the genus zero fundamental modular curve
(1.7)
j
2
·
j
'
2
-
(
j
+
j
'
)
·
(
j
2
+
1487
·
j
j
'
+
j
'
2
)
+
3
·
1
5
3
·
(
16
j
2
-
4027
j
j
'
+
16
j
'
2
)
-
12
·
3
0
6
·
(
j
+
j
'
)
+
8
·
3
0
9
=
0
,
or(1.8)
5
9
v
3
u
3
-
12
·
5
6
u
2
v
2
·
(
u
+
v
)
+
375
u
v
·
(
16
u
2
+
16
v
2
-
4027
v
u
)
-
64
(
v
+
u
)
·
(
v
2
+
1487
v
u
+
u
2
)
+
2
12
·
3
3
·
u
v
=
0
,
which relates the two Hauptmoduls u
=
1
2
3
/
j
(
k
), v
=
1
2
3
/
j
(
k
L
):(1.9)
j
(
k
)
=
256
·
(
1
-
k
2
+
k
4
)
3
k
4
·
(
1
-
k
2
)
2
,
j
(
k
L
)
=
16
·
(
1
+
14
k
2
+
k
4
)
3
(
1
-
k
2
)
4
·
k
2
.
One verifies easily that (1.7) is verified with j
=
j
(
k
) and j
'
=
j
(
k
L
).The selected values ofk, the modulus of elliptic functions, k
=
0,1, are actually fixed points of the Landen transformations. The Kramers-Wannier duality k
↔
1
/
k maps k
=
0 onto k
=
∞. For the Ising (resp. Baxter) model these selected values of k correspond to the three selected subcases of the model (T
=
∞, T
=
0, and the critical temperature T
=
T
c), for which the elliptic parametrization of the model degenerates into a rational parametrization [4]. We have the same property for all the other algebraic modular curves corresponding to τ
↔
N
·
τ. This is certainly the main property most physicists expect for an exact representation of a generator of the renormalization group, namely, that it maps a generic point of the parameter space onto the critical manifold (fixed points). Modular transformations are, in fact, the only transformations to be compatible with all the other symmetries of the Ising (resp. Baxter) model like, for instance, the gauge transformations, some extended s
l
(
2
)
×
s
l
(
2
)
×
s
l
(
2
)
×
s
l
(
2
) symmetry [7], and so forth. It has also been underlined in [3, 4] that seeing (1.6) as a transformation on complex variables (instead of real variables) provides two other complex fixed points which actually correspond to complex multiplication for the elliptic curve, and are, actually, fundamental new singularities9 discovered on the χ
(
3
) linear ODE [8–10]. In general, this underlines the deep relation between the renormalization group and the theory of elliptic curves in a deep sense, namely, isogenies of elliptic curves, Hauptmoduls,
10
modular curves and modular forms.Note that an algebraic transformation like (1.6) or (1.8) cannot be obtained from any local Migdal-Kadanoff transformation which naturally yields rational transformations; an exact renormalization group transformation like (1.6) can only be deduced from nonlocal decimations. The emergence of modular transformations as representations of exact generators of the renormalization group explains, in a quite subtle way, the difficult problem of how renormalization group transformations can be compatible with reversibility
11 (iteration forward and backwards). An algebraic modular transformation (1.8) corresponds to τ
→
2
τ
and
τ
→
τ
/
2
in the same time, as a consequence of the modular group symmetry
τ
↔
1
/
τ.A simple rational parametrization12 of the genus zero modular curve (1.8) reads:(1.10)
u
=
1728
z
(
z
+
16
)
3
,
v
=
1728
z
2
(
z
+
256
)
3
=
u
(
2
12
z
)
.
Note that the previously mentioned reversibility is also associated with the fact that the modular curve (1.8) is invariant by u
↔
v, and, within the previous rational parametrization (1.10), with the fact that permuting u and v corresponds13 to the Atkin-Lehner involution z
↔
2
12
/
z.For many Yang-Baxter integrable models of lattice statistical mechanics the physical quantities (partition function per site, correlation functions, etc.) are solutions of selected14 linear differential equations. For instance, the partition function per site of the square (resp. triangular, etc.) Ising model is an integral of an elliptic integral of the third kind. It would be too complicated to show the precise covariance of these physical quantities with respect to (algebraic) modular transformations like (1.8). Instead, let us give, here, an illustration of the nontrivial action of the renormalization group on some elliptic function that actually occurs in the 2D Ising model: a weight-one modular form. This modular form actually, and remarkably, emerged [11] in a second-order linear differential operator factor denoted Z
2 occurring [8] for χ
(
3
), and that the reader can think as a physical quantity solution of a particular linear ODE replacing the too complicated integral of an elliptic integral of the third kind. Let us consider the second-order linear differential operator (D
z denotes d
/
d
z):(1.11)
α
=
D
z
2
+
(
z
2
+
56
z
+
1024
)
z
·
(
z
+
16
)
(
z
+
64
)
·
D
z
-
240
z
·
(
z
+
16
)
2
(
z
+
64
)
,
which has the (modular form) solution(1.12)
F
1
2
(
[
1
12
,
5
12
]
,
[
1
]
;
1728
z
(
z
+
16
)
3
)
=
2
·
(
z
+
256
z
+
16
)
-
1
/
4
·
F
1
2
(
[
1
12
,
5
12
]
,
[
1
]
;
1728
z
2
(
z
+
256
)
3
)
.
Do note that the two pull-backs in the arguments of the same hypergeometric function are actually related by the modular curve relation (1.8) (see (1.10)). The covariance (1.12) is thus the very expression of a modular form property with respect to a modular transformation (τ
↔
2
τ) corresponding to the modular transformation (1.8).The hypergeometric function at the rhs of (1.12) is solution of the second-order linear differential operator(1.13)
β
=
D
z
2
+
z
2
+
416
z
+
16384
(
z
+
256
)
(
z
+
64
)
z
·
D
z
-
60
(
z
+
64
)
(
z
+
256
)
2
,
which is the transformed of operator α by the Atkin-Lehner duality z
↔
2
12
/
z, and, also, a conjugation of α:(1.14)
β
=
(
z
+
16
z
+
256
)
-
1
/
4
·
α
·
(
z
+
16
z
+
256
)
1
/
4
.Along this line we can also recall that the (modular form) function15
(1.15)
F
(
j
)
=
j
-
1
/
12
·
F
2
1
(
[
1
12
,
5
12
]
,
[
1
]
;
1
2
3
j
)
,
verifies:(1.16)
F
(
(
z
+
16
)
3
z
)
=
2
·
z
-
1
/
12
·
F
(
(
z
+
256
)
3
z
2
)
.A relation like (1.12) is a straight generalization of the covariance we had in the one-dimensional model Z
(
t
)
=
C
(
t
)
·
Z
(
t
2
), which basically amounts to seeing the partition function per site as some “automorphic function” with respect to the renormalization group, with the simple renormalization group transformation t
→
t
2 being replaced by the algebraic modular transformation (1.8) corresponding to τ
↔
2
τ (i.e., the Landen transformation (1.6)).We have here all the ingredients for seeing the identification of exact algebraic representations of the renormalization group with the modular curves structures we tried so many times to promote (preaching in the desert) in various papers [3, 4]. However, even if there are no difficulties, just subtleties, these Ising-Baxter examples of exact algebraic representations of the renormalization group already require some serious knowledge of the modular curves, modular forms, and Hauptmoduls in the theory of elliptic curves, mixed with the subtleties naturally associated with the various branches of such algebraic (multivalued) transformations.The purpose of this paper is to present another elliptic hypergeometric function and other much simpler (Gauss hypergeometric) second-order linear differential operators covariant by infinite-order rational transformations.The replacement ofalgebraic (modular) transformations by simple rational transformations will enable us to display a complete explicit description of an exact representation of the renormalization group that any graduate student can completely dominate.


## 2. Infinite Number of Rational Symmetries on a Gauss Hypergeometric ODE

Keeping in mind modular form expressions like (1.12), let us recall a particular Gauss hypergeometric function introduced by Vidunas in [12](2.1)
F
2
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
=
1
4
·
z
-
1
/
4
·
∫
0
z
t
-
3
/
4
(
1
-
t
)
-
1
/
2
d
t
=
(
1
-
z
)
-
1
/
2
·
F
2
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
-
4
z
(
1
-
z
)
2
)
.
This hypergeometric function corresponds to the integral of a holomorphic form on a genus-one curve P
(
y
,
t
)
=
0:(2.2)
d
t
y
,
with
:
y
4
-
t
3
·
(
1
-
t
)
2
=
0
.
Note that the function(2.3)
ℱ
(
z
)
=
z
1
/
4
·
F
2
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
,
which is exactly an integral of an algebraic function, has an extremely simple covariance property with respect to the infinite-order rational transformation z
→
-
4
z
/
(
1
-
z
)
2:(2.4)
ℱ
(
-
4
z
(
1
-
z
)
2
)
=
(
-
4
)
1
/
4
·
ℱ
(
z
)
.
The occurrence of this specific infinite-order transformation is reminiscent of Kummer's quadratic relation (2.5)
F
2
1
(
[
a
,
b
]
,
[
1
+
a
-
b
]
;
z
)
=
(
1
-
z
)
-
a
·
F
2
1
(
[
a
2
,
1
+
a
2
-
b
]
,
[
1
+
a
-
b
]
;
-
4
z
(
1
-
z
)
2
)
,
but it is crucial to note that, relation (2.4) does not relate two different functions, but is an “automorphy” relation on the same function.It is clear from the previous paragraph that we want to see such functions as “ideal” examples of physical functions covariant by an exact (here, rational) generator of the renormalization group. The function (2.3) is actually solution of the second-order linear differential operator:(2.6)
Ω
=
D
z
2
+
1
4
3
-
5
z
z
·
(
1
-
z
)
·
D
z
=
ω
1
·
D
z
,
with
ω
1
=
D
z
+
1
4
3
-
5
z
z
·
(
1
-
z
)
=
D
z
+
1
4
·
d
ln
(
z
3
(
1
-
z
)
2
)
d
z
.
From the previous expression of ω
1 involving a log derivative of a rational function it is obvious that this second-order linear differential operator has two solutions, the constant function and an integral of an algebraic function. Since these two solutions behave very simply under the infinite-order rational transformation z
→
-
4
z
/
(
1
-
z
)
2, it is totally and utterly natural to see how the linear differential operator Ω transforms under the rational change of variable z
→
R
(
z
)
=
-
4
z
/
(
1
-
z
)
2 (which amounts to seeing how the two-order-one operators ω
1 and D
z transform). It is a straightforward calculation to see that introducing the cofactor C
(
z
) which is the inverse of the derivative of R
(
z
)
(2.7)
C
(
z
)
=
-
1
4
·
(
1
-
z
)
3
1
+
z
,
1
C
(
z
)
=
d
R
(
z
)
d
z
,
D
z and ω
1, respectively, transform under the rational change of variable z
→
R
(
z
)
=
-
4
z
/
(
1
-
z
)
2 as(2.8)
D
z
→
C
(
z
)
·
D
z
,
ω
1
→
(
ω
1
)
(
R
)
=
C
(
z
)
2
·
ω
1
·
1
C
(
z
)
,
yielding
:
Ω
→
C
(
z
)
2
·
Ω
.
Since z
→
-
4
z
/
(
1
-
z
)
2 is of infinite-order, the second-order linear differential operator (2.6) has an infinite number of rational symmetries (isogenies):(2.9)
z
→
-
4
z
(
1
-
z
)
2
→
16
·
(
1
-
z
)
2
·
z
(
1
+
z
)
4
→
-
64
·
(
1
-
z
)
2
(
1
+
z
)
4
z
(
1
-
6
z
+
z
2
)
4
→
⋯
.Once we have found a second-order linear differential operator (written in a unitary or monic form)Ω, covariant by the infinite-order rational transformation z
→
-
4
z
/
(
1
-
z
)
2, it is natural to seek for higher-order linear differential operators also covariant by z
→
-
4
z
/
(
1
-
z
)
2. One easily verifies that the successive symmetric powers of Ω are (of course) also covariant. The symmetric square of Ω,(2.10)
D
z
3
+
3
4
3
-
5
z
(
1
-
z
)
z
·
D
z
2
+
3
8
1
-
5
z
(
1
-
z
)
z
2
·
D
z
,
factorizes in simple order-one operators(2.11)
(
D
z
+
2
4
3
-
5
z
(
1
-
z
)
z
)
·
(
D
z
+
1
4
3
-
5
z
(
1
-
z
)
z
)
·
D
z
,
and, more generally, the symmetric Nth power16 of Ω reads(2.12)
(
D
z
+
N
4
3
-
5
z
z
(
1
-
z
)
)
·
(
D
z
+
N
-
1
4
3
-
5
z
z
(
1
-
z
)
)
⋯
(
D
z
+
1
4
3
-
5
z
z
(
1
-
z
)
)
·
D
z
.
The covariance of such expressions is the straight consequence of the fact that the order-one factors(2.13)
ω
k
=
D
z
+
k
4
3
-
5
z
z
·
(
1
-
z
)
,
k
=
0,1
,
…
,
N
,
transform very simply under z
→
-
4
z
/
(
1
-
z
)
2:(2.14)
ω
k
→
(
ω
k
)
(
R
)
=
(
C
(
z
)
)
k
+
1
·
ω
k
·
(
C
(
z
)
)
-
k
.
More generally, let us consider a rational transformation z
→
R
(
z
), the corresponding cofactor C
(
z
)
=
1
/
R
'
(
z
), and the order-one operator ω
1
=
D
z
+
A
(
z
). We have the identity(2.15)
C
(
z
)
·
D
z
·
(
1
C
(
z
)
)
=
D
z
-
d
ln
(
C
(
z
)
)
d
z
.
The change of variable z
→
R
(
z
) on ω
1 reads(2.16)
D
z
+
A
(
z
)
→
C
(
z
)
·
D
z
+
A
(
R
(
z
)
)
=
C
(
z
)
·
(
D
z
+
B
(
z
)
)
.
We want to impose that this rhs expression can be written (see (2.8)) as(2.17)
C
(
z
)
2
·
(
D
z
+
A
(
z
)
)
·
1
C
(
z
)
,
which, because of (2.15), occurs if(2.18)
B
(
z
)
=
A
(
z
)
-
d
ln
(
C
(
z
)
)
d
z
,
yielding a “Rota-Baxter-like” [13, 14] functional equation on A
(
z
) and R
(
z
)
(2.19)
(
d
R
(
z
)
d
z
)
2
·
A
(
R
(
z
)
)
=
d
R
(
z
)
d
z
·
A
(
z
)
+
d
2
R
(
z
)
d
z
2
.Remark 2.1.
Coming back to the initial Gauss hypergeometric differential operator the covariance ofΩ becomes a conjugation. Let us start with the Gauss hypergeometric differential operator for (2.1):
(2.20)
H
0
=
8
z
·
(
1
-
z
)
·
D
z
2
+
2
·
(
5
-
7
z
)
·
D
z
-
1
.
It is transformed by z
→
R
(
z
)
=
-
4
z
/
(
1
-
z
)
2 into
(2.21)
H
1
=
8
z
·
(
1
-
z
)
·
D
z
2
-
2
(
3
z
-
5
)
·
D
z
+
4
1
-
z
=
(
1
-
z
)
1
/
2
·
H
0
·
(
1
-
z
)
-
1
/
2
,
then by z
→
R
(
R
(
z
)
)
=
R
2
(
z
)
=
16
z
(
1
-
z
)
2
/
(
1
+
z
)
4 into
(2.22)
H
2
=
8
z
·
(
1
-
z
)
·
D
z
2
-
2
(
3
z
-
1
)
(
z
+
5
)
z
+
1
·
D
z
+
16
z
-
1
(
z
+
1
)
2
=
(
z
+
1
z
-
1
)
·
H
0
·
(
z
+
1
z
-
1
)
-
1
,
and more generally for z
→
R
N
=
R
(
R
(
R
⋯
(
R
(
z
)
⋯
)
(2.23)
H
N
=
C
N
·
H
0
·
C
N
-
1
,
where
:
C
N
=
z
1
/
4
·
R
N
-
1
/
4
.


### 2.1. A Few Remarks on the “Rota-Baxter-Like” Functional Equation

The functional equation17(2.19) is the (necessary and sufficient) condition for Ω
=
(
D
z
+
A
(
z
)
)
·
D
z to be covariant by z
→
R
(
z
).Using the chain rule formula of derivatives of composed functions:(2.24)
d
R
(
R
(
z
)
)
d
z
=
d
R
(
z
)
d
z
·
[
d
R
(
z
)
d
z
(
R
(
z
)
)
]
,
d
2
R
(
R
(
z
)
)
d
z
2
=
d
2
R
(
z
)
d
z
2
·
[
d
R
(
z
)
d
z
(
R
(
z
)
)
]
+
(
d
R
(
z
)
d
z
)
2
·
[
d
2
R
(
z
)
d
z
2
(
R
(
z
)
)
]
,
one can show that, for A
(
z
) fixed, the “Rota-Baxter-like” functional equation (2.19) is invariant by the composition of R
(
z
) by itself R
(
z
)
→
R
(
R
(
z
)
), R
(
R
(
R
(
z
)
)
)
,
…
. This result can be generalized to any composition of various R
(
z
)'s satisfying (2.19). This is in agreement with the fact that (2.19) is the condition for Ω
=
(
D
z
+
A
(
z
)
)
·
D
z to be covariant by z
→
R
(
z
) it must be invariant by composition of R
(
z
)'s (for A
(
z
) fixed).Note that we have not used here the fact that for globally nilpotent [11] operators, A
(
z
) and B
(
z
) are necessarily log derivatives of Nth roots of rational functions. For R
(
z
)
=
-
4
z
/
(
1
-
z
)
2:(2.25)
A
(
z
)
=
1
4
·
d
ln
(
a
(
z
)
)
d
z
,
B
(
z
)
=
1
4
·
d
ln
(
b
(
z
)
)
d
z
,
a
(
z
)
=
(
1
-
z
)
2
·
z
3
,
b
(
z
)
=
z
3
·
(
1
+
z
)
4
(
1
-
z
)
10
.
The existence of the underlying a
(
z
) in (2.25) consequence of a global nilpotence of the order-one differential operator, can however be seen in the following remark on the zeros of the lhs and rhs terms in the functional equation (2.19). When R
(
z
) is a rational function (e.g., -
4
z
/
(
1
-
z
)
2 or any of its iterates R
(
n
)
(
z
)), the lhs and rhs of (2.19) are rational expressions. The zeros are roots of the numerators of these rational expressions. Because of (2.25) the functional equation (2.19) can be rewritten (after dividing by R
'
(
z
)) as(2.26)
(
d
R
(
z
)
d
z
)
·
A
(
R
(
z
)
)
=
A
(
z
)
+
d
d
z
(
ln
(
d
R
(
z
)
d
z
)
)
=
1
4
·
d
d
z
(
ln
(
a
(
z
)
·
(
d
R
(
z
)
d
z
)
4
)
)
.
One easily verifies, in our example, that the zeros of the rhs of (2.26) come from the zeros of A
(
R
(
z
)
) (and not from the zeros of R
'
(
z
) in the lhs of (2.26)). The zeros of the log-derivative rhs of (2.26) correspond to a
(
z
)
·
R
'
(
z
)
4
=
ρ, where ρ is a constant to be found. Let us consider for R
(
z
) the nth iterates of -
4
z
/
(
1
-
z
)
2 that we denote R
(
n
)
(
z
). A straightforward calculation shows that the zeros of A
(
R
(
n
)
(
z
)
) or a
'
(
R
(
n
)
(
z
)
) (where a
'
(
z
) denotes the derivative of a
(
z
) namely, (
z
-
1
)
(
5
z
-
3
)
·
z
2) actually correspond to the general closed formula:(2.27)
5
5
·
a
(
z
)
·
(
d
R
(
n
)
(
z
)
d
z
)
4
-
4
·
3
3
·
(
-
4
)
n
=
0
.
More precisely the zeros of 5
·
R
(
n
)
(
z
)
-
3 verify (2.27), or, in other words, the numerator of 5
R
(
n
)
(
z
)
-
3 divides the numerator of the lhs of (2.27).In another case forT
(
z
) given by (2.45), which also verifies (2.19) (see below), the relation (2.27) is replaced by(2.28)
5
5
·
a
(
z
)
·
(
d
T
(
n
)
(
z
)
d
z
)
4
-
4
·
3
3
·
(
-
7
-
24
i
)
n
=
0
.
More generally for a rational function ρ
(
x
), obtained by an arbitrary composition of -
4
z
/
(
1
-
z
)
2 and T
(
z
), we would have(2.29)
5
5
·
a
(
z
)
·
(
d
ρ
(
z
)
d
z
)
4
-
4
·
3
3
·
λ
n
=
0
.
where λ corresponds to(2.30)
ρ
(
x
)
=
λ
·
z
+
⋯
,
λ
=
[
d
ρ
(
z
)
d
z
]
z
=
0
.


### 2.2. Symmetries ofΩ, Solutions to the “Rota-Baxter-Like” Functional Equation

Let us now analyse all the symmetries of the linear differential operatorΩ
=
(
D
z
+
A
(
z
)
)
·
D
z by analyzing all the solutions of (2.19) for a given A
(
z
). For simplicity we will restrict to A
(
z
)
=
(
3
-
5
z
)
/
z
/
(
1
-
z
)
/
4 which corresponds to R
(
z
)
=
-
4
z
/
(
z
-
1
)
2 and all its iterates (2.9). Let us first seek for other (more general) solutions that are analytic at
z
=
0:(2.31)
R
(
z
)
=
a
1
·
z
+
a
2
·
z
2
+
a
3
·
z
3
+
⋯
.
It is a straightforward calculation to get, order by order from (2.19), the successive coefficients a
n in (2.31) as polynomial expressions (with rational coefficients) of the first coefficient a
1 with(2.32)
a
2
=
-
2
5
·
a
1
·
(
a
1
-
1
)
,
a
3
=
1
75
·
a
1
·
(
a
1
-
1
)
·
(
7
a
1
-
17
)
,
a
4
=
-
2
4875
·
a
1
·
(
a
1
-
1
)
·
(
41
a
1
2
-
232
a
1
+
366
)
,
…
,
a
n
=
-
n
5
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
P
n
(
-
4
)
,
where P
n
(
a
1
) is a polynomial with integer coefficients of degree n
-
2. Since we have here a series depending on one parameter a
1 we will denote it R
a
1
(
z
). This is a quite remarkable series depending on one parameter.18 One can easily verify that this series actually reduces (as it should!) to the successive iterates (2.9) of -
4
z
/
(
1
-
z
)
2 for a
1
=
(
-
4
)
n. In other words this one-parameter family of “functions” actually reduces to rational functions for an infinite number of integer values a
1
=
(
-
4
)
n.Furthermore, one can also verify a quite essential property we expect for a representation of the renormalization group, namely, that twoR
a
1
(
z
) for different values of a
1 commute, the result corresponding to the product of these two a
1:(2.33)
R
a
1
(
R
b
1
(
z
)
)
=
R
b
1
(
R
a
1
(
z
)
)
=
R
a
1
·
b
1
(
z
)
.
The neutral element must necessarily correspond to a
1
=
1 which is actually the identity transformation R
1
(
z
)
=
z. We have an “absorbing” element corresponding to a
1
=
0, namely, R
0
(
z
)
=
0. Performing the inverse of R
a
1
(
z
) (with respect to the composition of functions) amounts to changing a
1 into its inverse 1
/
a
1. Let us explore some “reversibility” property of our exact representation of a renormalization group with the inverse of the rational transformations (2.9). The inverse of R
-
4
(
z
)
=
-
4
z
/
(
1
-
z
)
2 must correspond to a
1
=
-
1
/
4:(2.34)
R
-
1
/
4
(
z
)
=
-
1
4
·
z
-
1
8
z
2
-
5
64
z
3
-
7
128
z
4
-
21
512
z
5
+
⋯
.
However, a straight calculation of the inverse of R
-
4
(
z
)
=
-
4
z
/
(
1
-
z
)
2 gives a multivalued function, or if one prefers, two functions(2.35)
S
-
1
/
4
(
1
)
(
z
)
=
z
-
2
+
2
1
-
z
z
=
-
1
4
·
z
-
1
8
z
2
+
⋯
,
S
-
1
/
4
(
2
)
(
z
)
=
z
-
2
-
2
1
-
z
z
=
-
4
z
+
2
+
1
4
z
+
1
8
z
2
+
⋯
,
which are the two roots of the simple quadratic relation (R
-
4
(
z
'
)
=
z):(2.36)
z
'
2
-
2
·
(
1
-
2
z
)
·
z
'
+
1
=
0
,
where it is clear that the product of these two functions is equal to +
1. The radius of convergence of S
-
1
/
4
(
1
)
(
z
) is 1.Because of our choice to seek for functions analytical atz
=
0 our renormalization group representation “chooses” the unique root that is analytical at z
=
0, namely, S
-
1
/
4
(
1
)
(
z
). For the next iterate of R
-
4
(
z
)
=
-
4
z
/
(
1
-
z
)
2 in (2.9) the inverse transformation corresponds to the roots of the polynomial equation of degree four (R
16
(
z
'
)
=
z):(2.37)
z
'
4
+
(
4
-
16
z
)
·
z
′
3
+
(
6
+
32
z
)
·
z
′
2
+
(
4
-
16
z
)
·
z
'
+
1
=
0
,
which yields four roots, one of which is analytical at z
=
0 and corresponds to a
1
=
1
/
(
-
4
)
2 in our one-parameter family of (renormalization) transformations:(2.38)
S
1
/
16
(
1
)
(
z
)
=
1
16
z
+
3
128
z
2
+
53
4096
z
3
+
277
32768
z
4
+
3181
524288
z
5
+
⋯
,
its (multiplicative) inverse S
1
/
16
(
2
)
(
z
)
=
1
/
S
1
/
16
(
1
)
(
z
):(2.39)
S
1
/
16
(
2
)
(
z
)
=
16
z
-
6
-
17
16
z
-
67
128
z
2
-
1333
4096
z
3
-
7445
32768
z
4
+
⋯
,
and two (formal) Puiseux series (u
=
±
z):(2.40)
S
1
/
16
(
3
)
(
z
)
=
1
+
u
+
1
2
u
2
+
3
8
u
3
+
1
4
u
4
+
27
128
u
5
+
5
32
u
6
+
⋯
.
Many of these results are better understood when one keeps in mind that there is a special transformation J
:
z
↔
1
/
z which is also a R-solution of (2.19) and verifies many compatibility relations with these transformations (I
d denotes the identity transformation R
0
(
z
)):(2.41)
R
-
4
·
J
=
R
-
4
,
S
-
1
/
4
(
2
)
·
R
-
4
=
J
,
R
-
4
·
S
-
1
/
4
(
1
)
=
S
-
1
/
4
(
1
)
·
R
-
4
=
I
d
,
S
1
/
16
(
1
)
(
z
)
=
S
-
1
/
4
(
1
)
·
S
-
1
/
4
(
1
)
,
S
1
/
16
(
2
)
(
z
)
=
S
-
1
/
4
(
1
)
·
S
-
1
/
4
(
2
)
,
J
·
S
-
1
/
4
(
1
)
=
S
-
1
/
4
(
2
)
,
J
·
S
-
1
/
4
(
2
)
=
S
-
1
/
4
(
1
)
,
…
,
where the dot corresponds, here, to the composition of functions. These symmetries of the linear differential operator Ω correspond to isogenies of the elliptic curve (2.2).It is clear that we have another one-parameter family corresponding toJ
·
R
a
1 with an expansion of the form(2.42)
J
·
R
a
1
=
b
1
z
-
2
5
·
(
b
1
-
1
)
-
1
15
·
b
1
2
-
1
b
1
·
z
-
2
975
·
(
b
1
-
1
)
(
4
b
1
+
1
)
(
4
b
1
+
3
)
b
1
2
·
z
2
-
1
248625
·
(
b
1
-
1
)
(
4
b
1
+
1
)
(
1268
b
1
2
+
951
b
1
+
91
)
b
1
3
·
z
3
-
2
2071875
·
(
b
1
-
1
)
(
4
b
1
+
1
)
(
3688
b
1
3
+
2766
b
1
2
+
404
b
1
+
17
)
b
1
4
·
z
4
+
⋯
.
For b
1
=
-
1
/
4, b
1
=
(
-
1
/
4
)
2, b
1
=
(
-
1
/
4
)
3, this family reduces to the (multiplicative) inverse of the successive rational functions displayed in (2.9)(2.43)
-
1
4
·
(
1
-
z
)
2
z
→
1
16
·
(
1
+
z
)
4
(
1
-
z
)
2
·
z
→
-
1
64
·
(
1
-
6
z
+
z
2
)
4
(
1
-
z
)
2
(
1
+
z
)
4
·
z
→
⋯
,
which can also be written as:
(2.44)
-
1
4
·
(
z
+
1
z
)
+
1
2
,
1
16
·
(
z
+
1
z
)
+
3
8
+
z
(
1
-
z
)
2
,
-
1
64
·
(
z
+
1
z
)
+
13
32
-
z
4
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
,
1
256
·
(
z
+
1
z
)
+
51
128
+
z
16
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
+
z
(
-
4
)
n
-
6
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
,
-
1
1024
·
(
z
+
1
z
)
+
205
512
-
z
164
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
-
4
z
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
-
64
z
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
(
1
+
20
z
-
26
z
2
+
20
z
3
+
z
4
)
4
,
…
,
1
(
-
4
)
n
·
(
z
+
1
z
)
+
2
5
4
n
(
4
n
-
(
-
1
)
n
)
+
z
(
-
4
)
n
-
2
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
+
z
(
-
4
)
n
-
6
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
+
z
(
-
4
)
n
-
8
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
(
1
+
20
z
-
26
z
2
+
20
z
3
+
z
4
)
4
+
⋯
,
where we discover some “additive structure” of these successive rational functions.In fact, due to the specificity of this elliptic curve (occurrence of complex multiplication), we have another remarkable rational transformation solution of (2.19), preserving covariantly Ω. Let us introduce the rational transformation (i denotes -
1):(2.45)
T
(
z
)
=
z
·
(
z
-
(
1
+
2
i
)
1
-
(
1
+
2
i
)
·
z
)
4
,
we also have the remarkable covariance [12]:(2.46)
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
=
1
-
z
/
(
1
+
2
i
)
1
-
(
1
+
2
i
)
z
·
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
T
(
z
)
)
,
which can be rewritten in a simpler way on (2.3) (see (2.4)).It is a straightforward matter to see thatT
(
z
) actually belongs to the R
a
1
(
z
) one-parameter family:(2.47)
T
(
z
)
=
R
a
1
(
z
)
=
-
(
7
+
24
i
)
·
z
+
⋯
,
a
1
=
-
25
·
ρ
,
ρ
=
(
7
+
24
i
)
25
,
|
ρ
|
=
1
.As far as the reduction of (2.32) to a rational function is concerned, it is straightforward to see that:(2.48)
(
1
-
z
)
2
·
(
1
+
z
)
4
·
R
a
1
(
z
)
=
a
1
·
z
+
⋯
-
2
175746796875
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
4
)
·
(
a
1
-
16
)
·
P
8
(
a
1
)
·
z
8
+
⋯
-
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
4
)
·
(
a
1
-
16
)
·
P
n
(
a
1
)
·
z
n
+
⋯
,
where N
(
n
) is a large integer growing with n, and P
n is a polynomial with integer coefficients of degree n
-
4, or(2.49)
(
1
-
(
1
+
2
i
)
·
z
)
4
·
R
a
1
(
z
)
=
a
1
·
z
+
⋯
−
4
1243125
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
7
+
24
i
)
·
(
P
6
(
a
1
)
+
i
Q
6
(
a
1
)
)
·
z
6
+
⋯
+
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
7
+
24
i
)
·
(
P
n
(
a
1
)
+
i
Q
n
(
a
1
)
)
·
z
n
+
⋯
,
where P
n and Q
n are two polynomials with integer coefficients of degree, respectively, n
-
3 and n
-
4.Similar calculations can be performed forT
*
(
z
) defined by(2.50)
T
*
(
z
)
=
z
·
(
z
-
(
1
-
2
i
)
(
1
-
2
i
)
z
-
1
)
4
,
for which we also have the covariance(2.51)
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
=
1
-
z
/
(
1
-
2
i
)
1
-
(
1
-
2
i
)
z
·
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
T
*
(
z
)
)
.It is a simple calculation to check that any iterate ofT
(
z
) (resp. T
*
(
z
)) is actually a solution of (2.19) and corresponds to R
a
1
(
z
) for the infinite number of values a
1
=
(
-
7
-
24
i
)
N (resp. (
-
7
+
24
i
)
N). Furthermore, one verifies, as it should (see (2.33)), that the three rational functions R
-
4
(
z
), T
(
z
), and T
*
(
z
) commute. It is also a straightforward calculation to see that the rational function built from any composition of R
-
4
(
z
), T
(
z
), and T
*
(
z
) is actually a solution of (2.19). We thus have a triple infinity of values of a
1, namely a
1
=
(
-
4
)
M
·
(
-
7
-
24
i
)
N
·
(
-
7
+
24
i
)
P
for any integer
M
,
N
and
P, for which R
a
1
(
z
) reduces to rational functions. We are in fact describing (some subset of) the isogenies of the elliptic curve (2.2), and identifying these isogenies with a discrete subset of the renormalization group. Conversely, a functional equation like (2.19) can be seen as a way to extend the n-fold composition of a rational function R
(
z
) (namely R
(
R
(
⋯
R
(
z
)
⋯
)
)) to n
any complex number.


### 2.3. Revisiting the One-Parameter Family of Solutions of the “Rota-Baxter-Like” Functional Equation

This extension can be revisited as follows. Keeping in mind the well-known example of the parametrization of the standard mapz
→
4
z
·
(
1
-
z
) with z
=
sin
2
(
θ
), yielding θ
→
2
θ, let us seek for a (transcendental) parametrization z
=
P
(
u
) such that(2.52)
R
-
4
(
P
(
u
)
)
=
P
(
-
4
u
)
or
:
R
-
4
=
P
·
H
-
4
·
P
-
1
,
where H
a
1 denotes the scaling transformation z
→
a
1
·
z (here H
-
4
:
z
→
-
4
·
z) and P
-
1 denotes the inverse transformation of P (for the composition). One can easily find such a (transcendental) parametrization order by order(2.53)
P
(
z
)
=
z
-
2
5
z
2
+
7
75
z
3
-
82
4875
z
4
+
1078
414375
z
5
−
452
1243125
z
6
+
57311
1212046875
z
7
-
1023946
175746796875
z
8
+
⋯
,
and similarly for its inverse (for the composition) transformation(2.54)
Q
(
z
)
=
P
-
1
(
z
)
=
z
+
2
5
z
2
+
17
75
z
3
+
244
1625
z
4
+
45043
414375
z
5
+
2302
27625
z
6
+
128941
1939275
z
7
+
15365176
281194875
z
8
+
⋯
.
This approach is reminiscent of the conjugation introduced in Siegel's theorem [15–17]. It is a straightforward matter to see (order by order) that one actually has(2.55)
R
a
1
(
P
(
u
)
)
=
P
(
a
1
·
u
)
or
:
R
a
1
=
P
·
H
a
1
·
P
-
1
.
The structure of the (one-parameter) renormalization group and the extension of the composition of n times a rational function R
(
z
) (namely, R
(
R
(
⋯
R
(
z
)
⋯
)
)) to n
any complex number become a straight consequence of this relation. Along this line one can define some “infinitesimal composition” (ϵ
≃
0):(2.56)
R
1
+
ϵ
(
z
)
=
P
·
H
1
+
ϵ
·
P
-
1
(
z
)
=
z
+
ϵ
·
F
(
z
)
+
⋯
,
where one can find, order by order, the “infinitesimal composition” function F
(
z
):(2.57)
F
(
z
)
=
z
-
2
5
z
2
-
2
15
z
3
-
14
195
z
4
-
154
3315
z
5
-
22
663
z
6
-
418
16575
z
7
-
9614
480675
z
8
-
2622
160225
z
9
+
⋯
.
It is straightforward to see, from (2.33), that the function F
(
z
) satisfies the following functional equations involving a rational function R
(
z
) (in the one-parameter family R
a
1
(
z
)):(2.58)
d
R
(
z
)
d
z
·
F
(
z
)
=
F
(
R
(
z
)
)
,
d
R
(
n
)
(
z
)
d
z
·
F
(
z
)
=
F
(
R
(
n
)
(
z
)
)
,
where
:
R
(
n
)
(
z
)
=
R
(
R
(
⋯
R
(
z
)
)
⋯
)
.
F
(
z
) cannot be a rational or algebraic function. Let us consider the fixed points of R
(
n
)
(
z
). Generically d
R
(
n
)
(
z
)
/
d
z is not equal to 0 or ∞ at any of these fixed points. Therefore one must have F
(
z
)
=
0 or F
(
z
)
=
∞ for the infinite set of these fixed points: F
(
z
) cannot be a rational or algebraic function, it is a transcendental function, and similarly for the parametrization function P
(
z
). In fact, let us introduce the function(2.59)
G
(
z
)
=
(
1
-
z
)
·
F
(
z
)
,
G
(
z
)
=
z
-
7
5
z
2
+
4
15
z
3
+
4
65
z
4
+
28
1105
z
5
+
44
3315
z
6
+
44
5525
z
7
+
836
160225
z
8
+
1748
480675
z
9
+
⋯
+
g
n
·
z
n
+
⋯
.
One actually finds that the successive g
n satisfies the very simple (hypergeometric function) relation:(2.60)
g
n
+
1
g
n
=
4
n
-
9
4
n
+
1
.
The function G
(
z
) is actually the hypergeometric function solution of the homogeneous operator(2.61)
D
z
2
+
1
4
13
z
-
3
z
·
(
1
-
z
)
·
D
z
+
3
4
6
z
2
-
3
z
+
1
(
1
-
z
)
2
·
z
2
,
or of the inhomogeneous ODE(2.62)
4
z
·
(
1
-
z
)
·
d
G
(
z
)
d
z
+
(
9
z
-
3
)
·
G
(
z
)
-
z
·
(
1
-
z
)
2
=
0
.
One deduces the expression of F
(
z
) as a hypergeometric function(2.63)
F
(
z
)
=
z
·
(
1
-
z
)
1
/
2
·
F
2
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
z
)
=
∂
R
a
1
∂
a
1
|
a
1
=
1
.Finally we get the linear differential operator annihilatingF
(
z
)
(2.64)
Ω
F
=
D
z
2
+
1
4
·
5
z
-
3
z
(
1
-
z
)
·
D
z
+
1
4
·
3
-
6
z
+
5
z
2
(
1
-
z
)
2
z
2
=
D
z
·
(
D
z
-
1
4
·
3
-
5
z
z
·
(
1
-
z
)
)
,
which is, in fact, nothing but Ω
* the adjoint of linear differential operator Ω (see (2.6)). One easily checks19 that the second-order differential equation Ω
F
(
y
(
z
)
)
=
0 transforms under the change of variable z
→
-
4
z
/
(
1
-
z
)
2 into the second-order differential equation Ω
F
(
R
)
(
y
(
z
)
)
=
0 with Ω
F
(
R
)
=
C
(
z
)
2
·
ω
F
(
R
) where the unitary (monic) operator ω
F
(
R
) is the conjugate of Ω
F:(2.65)
ω
F
(
R
)
=
D
z
2
-
1
4
·
11
z
2
+
30
z
+
3
z
·
(
1
-
z
)
(
1
+
z
)
·
D
z
+
1
4
·
3
+
12
z
+
50
z
2
+
12
z
3
+
3
z
4
z
2
·
(
1
-
z
)
2
(
1
+
z
)
2
=
(
1
C
(
z
)
)
·
D
z
·
(
D
z
-
1
4
·
3
-
5
z
z
·
(
1
-
z
)
)
·
C
(
z
)
=
(
1
C
(
z
)
)
·
Ω
F
·
C
(
z
)
=
(
1
C
(
z
)
)
·
Ω
*
·
C
(
z
)
with C
(
z
)
=
1
/
R
'
(
z
) and the “dot” denotes the composition of operators. Actually, the factors in the adjoint Ω
* transform under the change of variable z
→
-
4
z
/
(
1
-
z
)
2 as follows20:(2.66)
D
z
→
C
(
z
)
·
D
z
,
ω
1
*
=
→
(
ω
1
*
)
(
R
)
=
ω
1
*
·
C
(
z
)
,
Ω
*
→
Ω
F
(
R
)
=
C
(
z
)
·
Ω
*
·
C
(
z
)
which is precisely the transformation we need to match with (2.58) and see the ODE Ω
*
(
F
(
z
)
)
=
0 compatible with the change of variable z
→
-
4
z
/
(
1
-
z
)
2:(2.67)
Ω
*
(
F
(
z
)
)
=
0
→
(
C
(
z
)
·
Ω
*
·
C
(
z
)
)
(
F
(
R
(
z
)
)
)
=
(
C
(
z
)
·
Ω
*
·
C
(
z
)
)
(
R
′
(
z
)
·
F
(
z
)
)
=
C
(
z
)
·
Ω
*
(
F
(
z
)
)
=
0
.
This is, in fact, a quite general result that will be seen to be valid in a more general (higher genus) framework (see (2.148), (2.150) in what follows).Not surprisingly one can deduce from (2.33) and the previous results, in particular (2.63), the following results for R
a
1
(
z
):(2.68)
-
4
·
∂
R
a
1
∂
a
1
|
a
1
=
-
4
=
F
(
R
(
z
)
)
,
(
-
4
)
n
·
∂
R
a
1
∂
a
1
|
a
1
=
(
-
4
)
n
=
F
(
R
(
n
)
(
z
)
)
,
where R
(
z
)
=
-
4
z
/
(
1
-
z
)
2 and R
(
n
)
(
z
) denotes R
(
R
(
⋯
R
(
R
(
z
)
)
)
). Of course we have similar relation for T
(
z
), -
4 being replaced by -
7
-
24
i. Therefore the partial derivative ∂
R
a
1
/
∂
a
1 that can be expressed in terms of hypergeometric functions for for a double infinity of values of a
1, namely, a
1
=
(
-
4
)
M
×
(
-
7
-
24
i
)
N.One can, of course, check, order by order, that (2.58) is actually verified for any function in the one-parameter family R
a
1
(
z
):(2.69)
d
R
a
1
(
z
)
d
z
·
F
(
z
)
=
F
(
R
a
1
(
z
)
)
,
which corresponds to an infinitesimal version of (2.33).From (2.56) one simply deduces(2.70)
z
·
d
P
(
z
)
d
z
=
F
(
P
(
z
)
)
,
that we can check, order by order from (2.53), the series expansion of P
(
z
), and from (2.57) the series expansion of F
(
z
), but also(2.71)
d
Q
(
z
)
d
z
·
F
(
z
)
=
Q
(
z
)
that we can, check order by order, from (2.54), the series expansion of Q
(
z
)
=
P
-
1
(
z
) and from (2.57). We now deduce that the log-derivative of the “well-suited change of variable” Q
(
z
) is nothing but the (multiplicative) inverse of a hypergeometric function F
(
z
):(2.72)
d
ln
(
Q
(
z
)
)
d
z
=
1
F
(
z
)
,
Q
(
z
)
=
λ
·
exp
(
∫
z
d
z
F
(
z
)
)
.
The function Q
(
z
) is solution of the nonlinear differential equation(2.73)
-
4
z
2
·
(
1
-
z
)
2
·
(
Q
·
Q
(
1
)
·
Q
(
3
)
+
(
Q
(
1
)
)
2
·
Q
(
2
)
-
2
Q
·
(
Q
(
2
)
)
2
)
+
z
·
(
3
-
5
z
)
(
1
-
z
)
·
Q
(
1
)
·
(
Q
·
Q
(
2
)
-
(
Q
(
1
)
)
2
)
+
(
5
z
2
-
6
z
+
3
)
·
Q
·
(
Q
(
1
)
)
2
=
0
,
where the Q
(
n
)'s denote the nth derivative of Q
(
z
). At first sight Q
(
z
) would be a nonholonomic function, however, remarkably, it is a holonomic function solution of an order-five operator which factorizes as follows:(2.74)
Ω
Q
=
(
D
z
+
3
-
5
z
(
1
-
z
)
·
z
)
·
(
D
z
+
3
4
·
3
-
5
z
(
1
-
z
)
·
z
)
·
×
(
D
z
+
2
4
·
3
-
5
z
(
1
-
z
)
·
z
)
·
(
D
z
+
1
4
·
3
-
5
z
(
1
-
z
)
·
z
)
·
D
z
,
yielding the exact expression of Q
(
z
) in terms of hypergeometric functions:(2.75)
Q
(
z
)
=
z
·
(
F
2
1
(
z
1
-
z
[
1
2
,
1
4
]
,
[
5
4
]
;
z
z
1
-
z
)
)
4
=
z
1
-
z
·
(
2
F
1
(
[
1
4
,
3
4
]
,
[
5
4
]
;
-
z
1
-
z
)
)
4
,
that is, the fourth power of (2.3), with the differential operator (2.74) being the symmetric fourth power of Ω. From (2.3) we immediately get the covariance of Q
(
z
):(2.76)
Q
(
-
4
z
(
1
-
z
)
2
)
=
-
4
·
Q
(
z
)
,
and, more generally, Q
(
R
a
1
)
=
a
1
·
Q
(
z
). Since Q
(
z
) and F
(
z
) are expressed in terms of the same hypergeometric function, the relation (2.71) must be an identity on that hypergeometric function. This is actually the case. This hypergeometric function verifies the inhomogeneous equation:(2.77)
4
·
z
·
d
ℋ
(
z
)
d
z
+
ℋ
(
z
)
-
(
1
-
z
)
-
1
/
2
=
0
,
where(2.78)
ℋ
(
z
)
=
F
2
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
.RecallingQ
(
P
(
z
)
)
=
z, one has the following functional relation on P
(
z
):(2.79)
P
(
z
)
·
F
2
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
P
(
z
)
)
4
=
z
.Noting thatQ
(
z
4
)
1
/
4
=
ℱ
(
z
4
) (see (2.3)) can be expressed in term of an incomplete elliptic integral of the first kind of argument -
1
(2.80)
z
·
2
F
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
z
4
)
=
Elliptic
F
(
z
,
-
1
)
,
one can find that (2.79) rewrites on P
(
z
) as(2.81)
Elliptic
F
(
P
(
z
)
1
/
4
,
-
1
)
=
z
1
/
4
,
from which we deduce that the function P
(
z
) is nothing but aJacobi elliptic function
21
(2.82)
P
(
z
)
=
(
sn
(
z
1
/
4
,
-
1
)
)
4
.
In Appendix B we display a set of “Painlevé-like” ODEs22 verified by P
(
z
). From the simple nonlinear ODE on the Jacobi elliptic sinus, namely, S
′
′
+
2
·
S
3
=
0, and the exact expression of P
(
z
) in term of Jacobi elliptic sinus, one can deduce other nonlinear ODEs verified by the nonholonomic function P
(
z
) (P
(
1
)
=
d
P
(
z
)
/
d
z, P
(
2
)
=
d
2
P
(
z
)
/
d
z
2):(2.83)
z
3
/
2
·
(
P
(
1
)
)
2
-
(
1
-
P
)
·
P
3
/
2
=
0
,
(2.84)
P
(
2
)
-
3
4
·
(
P
(
1
)
)
2
P
+
3
4
·
P
(
1
)
z
+
1
2
·
P
3
/
2
z
3
/
2
=
0
.


### 2.4. Singularities of the Jacobi Elliptic FunctionP
(
z
)

Most of the results of this section, and to some extent, of the next one, are straight consequences of the exact closed expression ofP
(
z
) in terms of an elliptic function. Following the pedagogical approach of this paper we will rather follow a heuristic approach not taking into account the exact result (2.82), to display simple methods and ideas that can be used beyond exact results on a specific example.From a diff-Padé analysis of the series expansion ofP
(
z
), we got the sixty (closest to z
=
0) singularities. In particular we got that P
(
z
) has a radius of convergence R
≃
11.81704500807
⋯ corresponding to the following (closest to z
=
0) singularity z
=
z
s of P
(
z
):(2.85)
z
s
=
-
11.817045008077115768316337283432582087420697
⋯
=
(
-
4
)
·
2
F
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
1
)
4
=
-
1
16
·
π
6
Γ
(
3
/
4
)
8
.
This singularity corresponds to a pole of order four: P
(
z
)
≃
(
z
-
z
s
)
-
4. The function P
(
z
) has many other singularities:(2.86)
3
4
·
z
s
,
(
161
±
240
i
)
·
z
s
,
(
-
7
±
24
i
)
·
z
s
,
(
-
119
±
120
i
)
·
z
s
,
⋯
5
4
·
z
s
,
(
41
±
840
i
)
·
z
s
,
(
-
527
±
336
i
)
·
z
s
,
(
-
1519
±
720
i
)
·
z
s
,
⋯
7
4
·
z
s
,
(
1241
±
2520
i
)
·
z
s
,
(
-
567
±
1944
i
)
·
z
s
,
(
-
3479
±
1320
i
)
·
z
s
,
⋯In fact, introducingx and y the real and imaginary part of these singularities in z
s units, one finds out that they correspond to the double infinity of points(2.87)
x
=
(
m
1
2
-
2
m
1
m
2
-
m
2
2
)
·
(
m
1
2
+
2
m
1
m
2
-
m
2
2
)
,
y
=
4
·
m
1
m
2
·
(
m
2
-
m
1
)
·
(
m
2
+
m
1
)
,
where m
1 and m
2 are two integers, and they all lie on the intersection of an infinite number of genus zero curves indexed by the fourth power of an integer M
=
m
4 (m
=
m
1 or m
=
m
2):(2.88)
2
12
M
4
-
2
11
·
x
·
M
3
-
2
7
·
(
17
y
2
+
14
x
2
)
·
M
2
-
2
5
·
x
·
(
8
x
2
+
7
y
2
)
·
M
+
y
4
=
0
.
The parametrization (2.87) describes not only the poles of P
(
z
) when m
1
+
m
2 is odd, but also the zeros of P
(
z
) when m
1
+
m
2 is even. This (infinite) proliferation of singularities confirms the nonholonomic character of P
(
z
).These results are simply inherited from (2.82). The zeros and poles of the elliptic sinus s
n
(
z
,
i
) correspond to two lattice of periods. Denoting 𝒦
1 and 𝒦
2 the two periods of the elliptic curve, the location of the poles and zeros reads, respectively,(2.89)
P
n
1
,
n
2
=
2
n
1
·
𝒦
1
+
(
2
n
2
+
1
)
·
𝒦
2
,
Z
n
1
,
n
2
=
2
n
1
·
𝒦
1
+
2
n
2
·
𝒦
2
,
𝒦
1
=
π
3
/
2
2
3
/
2
·
1
Γ
(
3
/
4
)
2
,
𝒦
2
=
(
1
-
-
1
)
·
K
1
making crystal clear the fact that we have complex multiplication for this elliptic curve. The formula (2.87) just amount to saying that the poles and zeros of s
n
(
z
1
/
4
,
i
) are located at P
n
1
,
n
2
4 and Z
n
1
,
n
2
4:(2.90)
P
n
1
,
n
2
4
=
-
z
s
4
·
(
(
2
n
1
+
2
n
2
+
1
)
+
i
·
(
2
n
2
+
1
)
)
4
,
Z
n
1
,
n
2
4
=
-
z
s
4
·
(
(
2
n
1
+
2
n
2
+
1
)
+
i
·
2
n
2
)
4
.
The correspondence with (2.87) is m
1
=
n
1
+
2
n
2
+
1, m
2
=
-
n
1 for the poles and m
1
=
n
1
+
2
n
2, m
2
=
-
n
1 for the zeros.Remark 2.2.
let us consider thea
1
→
∞ limit of the one-parameter series R
a
1 (see (2.31), (2.32)) rewriting R
a
1
(
z
) as R
̃
b
1
(
u
)
(2.91)
R
̃
b
1
(
u
)
=
R
a
1
(
z
)
,
with
:
b
1
=
1
a
1
,
u
=
z
b
1
.
In the a
1
→
∞ limit, that is the b
1
→
0 limit, one easily verifies, order by order in u, that R
̃
b
1
(
u
) becomes exactly the transcendental parametrization function (2.53):
(2.92)
R
̃
b
1
(
u
)
→
P
(
u
)
when
b
1
→
0
.
For a
1
=
(
-
4
)
n (n
→
∞), one finds that the radius of convergence of the R
a
1
(
z
) series becomes in the n
→
∞ limit R
n
≃
z
s
/
4
n, in agreement23 with (2.91).


### 2.5.P
(
z
) and an Infinite Number of Rational Transformations: The Sky Is the Limit

Note that some nonlinear ODEs associated withP
(
z
) and displayed in Appendix B, namely (B.3) and (B.10), and the functional equation (2.79), are invariant by the change of variable (
P
(
z
)
,
z
)
→
(
-
4
P
(
z
)
/
(
1
-
P
(
z
)
)
2
,
-
4
z
). In fact (B.3), (2.79), and (B.10) are invariant by (
P
(
z
)
,
z
)
→
(
-
4
P
(
z
)
/
(
1
-
P
(
z
)
)
2
,
-
4
z
), but also (
P
(
z
)
,
z
)
→
(
-
(
1
-
P
(
z
)
)
2
/
4
/
P
(
z
)
,
-
z
/
4
), and also by (
P
(
z
)
,
z
)
→
(
1
/
P
(
z
)
,
z
).The functionP
(
z
) satisfies the functional equation:(2.93)
P
(
-
4
·
z
)
=
-
4
P
(
z
)
1
-
P
(
z
)
2
,
but also(2.94)
P
(
(
-
7
-
24
i
)
·
z
)
=
T
(
P
(
z
)
)
,
P
(
(
-
7
+
24
i
)
·
z
)
=
T
*
(
P
(
z
)
)
,
and, more generally, as can be checked order by order on series expansions(2.95)
P
(
a
1
·
z
)
=
R
a
1
(
P
(
z
)
)
.
For example, considering the “good” branch (2.35) for the inverse of -
4
z
/
(
1
-
z
)
2, namely S
-
1
/
4
(
1
)
(
z
), we can even check, order by order, on the series expansions of P
(
z
) and S
-
1
/
4
(
1
)
(
z
) the functional relation(2.96)
S
-
1
/
4
(
1
)
(
P
(
z
)
)
=
P
(
-
z
4
)
,
valid for
|
P
(
z
)
|
<
1
since the radius of convergence of
S
-
1
/
4
(
1
)
(
z
)
is
1.Recalling the functional equations (2.94) it is natural to say that if P
(
z
) is singular at z
=
z
s, then, for almost all the rational functions, in particular T
(
z
) (resp. T
*
(
z
)) the T
(
P
(
z
)
) is also singular z
=
z
s, and thus, from (2.94), P
(
z
) is also singular at z
=
(
-
7
±
24
i
)
·
z
s. It is thus extremely natural to see the emergence of the infinite number of singularities in (2.87) of the form z
=
(
N
1
+
i
·
N
2
)
·
z
s, as a consequence of (2.95) together with a reduction of the one-parameter series R
a
1
(
z
) to a rational function for an infinite number of selected values of a
1, namely the N
1
+
i
·
N
2 in (2.87). This is actually the case for all the values displayed in (2.87). For instance, for a
1
=
3
4
=
81 we get the following simple rational function:(2.97)
R
81
(
z
)
=
z
·
(
z
2
+
6
z
-
3
3
z
2
-
6
z
-
1
)
4
,
for which it is straightforward to verify that this rational transformation commutes with T
(
z
), T
*
(
z
), -
4
z
/
(
1
-
z
)
2, and is a solution of the Rota-Baxter-like functional equation (2.19). The case a
1
=
5
4
=
625 in (2.87), is even simpler, since it just requires to compose T
(
z
) and T
*
(
z
)
(2.98)
R
625
(
z
)
=
T
(
T
*
(
z
)
)
=
T
*
(
T
(
z
)
)
=
z
·
(
z
2
-
2
z
+
5
5
z
2
-
2
z
+
1
)
4
·
(
1
-
12
z
-
26
z
2
+
52
z
3
+
z
4
1
+
52
z
-
26
z
2
-
12
z
3
+
z
4
)
4
,
which, again verifies (2.19) and commutes with all the other rational functions, in particular (2.97). We also obtained the rational function corresponding to a
1
=
7
4
=
2401, namely:
(2.99)
R
2401
(
z
)
=
z
·
(
N
2401
(
z
)
D
2401
(
z
)
)
4
,
with:(2.100)
N
2401
(
z
)
=
z
12
·
D
2401
(
1
z
)
,
(2.101)
D
2401
(
z
)
=
1
+
196
z
-
1302
z
2
+
14756
z
3
-
15673
z
4
-
42168
z
5
+
111916
z
6
-
82264
z
7
+
35231
z
8
-
19852
z
9
+
2954
z
10
+
308
z
11
-
7
z
12
.The polynomialN
2401
(
z
) satisfies many functional equations, like, for instance (withR
-
4
(
z
)
=
-
4
z
/
(
1
-
z
2
)):(2.102)
4
12
·
D
2401
(
1
R
-
4
(
z
)
)
=
D
2401
(
z
)
·
D
2401
(
1
z
)
and also(2.103)
(
1
-
z
)
49
·
D
2401
(
R
-
4
(
z
)
)
2
=
D
2401
(
z
)
4
-
z
49
·
D
2401
(
1
z
)
4
.We also obtained the rational function corresponding toa
1
=
1
1
4
=
14641, namely,
(2.104)
R
14641
(
z
)
=
z
·
(
N
14641
(
z
)
D
14641
(
z
)
)
4
,
with:(2.105)
N
14641
(
z
)
=
z
30
·
D
14641
(
1
z
)
,
(2.106)
D
14641
(
z
)
=
1
+
1210
z
-
33033
z
2
+
2923492
z
3
+
5093605
z
4
-
385382514
z
5
+
3974726283
z
6
-
14323974808
z
7
+
57392757037
z
8
-
291359180310
z
9
+
948497199067
z
10
-
1642552094436
z
11
+
1084042069649
z
12
+
1890240552750
z
13
-
6610669151537
z
14
+
9712525647792
z
15
-
8608181312269
z
16
+
5384207244702
z
17
-
3223489742187
z
18
+
2175830922716
z
19
-
1197743580033
z
20
+
387221579866
z
21
-
50897017743
z
22
-
7864445336
z
23
+
5391243935
z
24
-
815789634
z
25
+
28366041
z
26
-
5092956
z
27
+
207691
z
28
+
2794
z
29
-
11
z
30
,
and, of course, one can verify that R
14641
(
z
) actually commutes with R
-
4, R
81, R
625, R
2401
(
z
), and is a solution of the Rota-Baxter-like functional equation (2.19). Similarly to R
2401
(
z
) (see (2.102), (2.103)), we also have the functional equations:(2.107)
4
30
·
D
14641
(
1
R
-
4
(
z
)
)
=
D
14641
(
z
)
·
D
14641
(
1
z
)
,
and also(2.108)
(
1
-
z
)
(
4
·
30
+
1
)
·
D
14641
(
R
-
4
(
z
)
)
2
=
D
14641
(
z
)
4
-
z
(
4
·
30
+
1
)
·
D
14641
(
1
z
)
4
.Next we obtained the rational function corresponding toa
1
=
1
3
4
=
28561, which verifies (2.19) namely,(2.109)
R
28561
(
z
)
=
z
·
(
N
28561
(
z
)
D
28561
(
z
)
)
4
,
with
N
28561
(
z
)
=
z
42
·
D
28561
(
1
z
)
,
N
28561
(
z
)
=
z
42
·
D
28561
(
1
z
)
,
(2.110)
D
28561
(
z
)
=
(
1
-
22
z
+
235
z
2
-
228
z
3
+
39
z
4
+
26
z
5
+
13
z
6
)
·
D
28561
(
36
)
(
z
)
D
28561
(
36
)
(
z
)
=
1
+
2388
z
-
61098
z
2
+
19225300
z
3
+
606593049
z
4
-
1543922656
z
5
+
7856476560
z
6
-
221753896032
z
7
+
1621753072244
z
8
-
4542779886736
z
9
+
2731418674664
z
10
+
36717669656304
z
11
-
200879613202428
z
12
+
547249607666784
z
13
-
934179604482832
z
14
+
1235038888776160
z
15
-
1788854212778642
z
16
+
3018407750933816
z
17
-
4349780716415868
z
18
+
4419228090228152
z
19
-
2899766501472914
z
20
+
931940880451552
z
21
+
413258559018224
z
22
-
857795672629664
z
23
+
659989056851972
z
24
-
304241349909008
z
25
+
87636987790824
z
26
-
14593362219920
z
27
+
1073204980340
z
28
+
45138167200
z
29
-
23660433008
z
30
+
2028597792
z
31
-
29540327
z
32
+
3238420
z
33
-
73386
z
34
-
492
z
35
+
z
36
.
We get similar results, mutatis mutandis, than the ones previously obtained (commutation, functional equations like (2.107), (2.108), etc.), namely,(2.111)
4
42
·
D
28561
(
1
R
-
4
(
z
)
)
=
D
28561
(
z
)
·
D
28561
(
1
z
)
,
…
.The “palindromic” nature of (2.97) (2.98), (2.99), (2.104) and (2.109) (see (2.100), (2.105)), (2.109)) corresponds to the fact that these rational transformations commute with J:(2.112)
1
R
81
(
z
)
=
R
81
(
1
z
)
,
1
R
625
(
z
)
=
R
625
(
1
z
)
,
…
.
In fact, more generally, we have R
N
4
(
1
/
z
)
=
1
/
R
N
4
(
z
) for N any odd integer (N
=
9,21
,
…) and R
N
4
(
1
/
z
)
=
R
N
4
(
z
) for any even integer  N.From (2.87) one can reasonnably conjecture that the fourth power of any integer will provide a new example of R
a
1
(
z
) being a rational function. The simple nontrivial example corresponds to the already found rational function(2.113)
R
16
(
z
)
=
16
·
z
·
(
1
-
z
)
2
(
z
+
1
)
4
.
We already have explicit rational functions for all values of a
1 of the form N
4 for N
=
2,3
,
…
,
16 and of course, we can in principle, build explicit rational functions for all the N's product of the previous integers. Along this line it is worth noticing that the coefficients of the series R
a
1
(
z
) are all integers when a
1 is the fourth power of any integer.We are thus starting to build an infinite number of (elementary)commuting rational transformations, any composition of these (infinite number of) rational transformations giving rational transformations satisfying (2.19) and preserving the linear differential operator Ω. This set of rational transformations is a pretty large set! Actually this set of rational transformations corresponds to the isogenies of the underlying elliptic function.The proliferation of the singularities ofP
(
z
) corresponds to this (pretty large) set of rational transformations. Recalling (2.96), the previous singularity argument is not valid
24 for the (well-suited) inverse transformations (S
-
1
/
4
(
1
)
(
z
), … ) of these rational transformations because (2.96) requires |
P
(
z
)
|
<
1 (corresponding to the radius of convergence of S
-
1
/
4
(
1
)
(
z
)) and the singularity z
=
z
s corresponds precisely to “hit” the value P
(
z
)
=
1.


### 2.6. Other Examples of Selected Gauss Hypergeometric ODEs

For heuristic reasons we have focused onA
(
z
)
=
(
3
-
5
z
)
/
z
/
(
1
-
z
)
/
4, but of course, one can find many other examples and try to generalize these examples.For instance, introducing(2.114)
A
(
z
)
=
1
6
·
d
ln
(
(
1
-
z
)
3
z
5
)
d
z
=
1
6
·
5
-
8
z
(
1
-
z
)
z
the rational transformation(2.115)
R
(
z
)
=
-
27
·
z
(
1
-
4
z
)
3
,
verifies the “Rota-Baxter-like” functional relation (2.19). This example corresponds to the following covariance [12] on a Gauss hypergeometric integral (of the c
=
1
+
b type, see below):(2.116)
2
F
1
(
[
1
2
,
1
6
]
,
[
7
6
]
;
z
)
=
z
-
1
/
6
6
·
∫
0
z
t
-
5
/
6
(
1
-
t
)
-
1
/
2
·
d
t
=
(
1
-
4
z
)
-
1
/
2
·
2
F
1
(
[
1
2
,
1
6
]
,
[
7
6
]
;
-
27
z
(
1
-
4
z
)
3
)
which is associated with the elliptic curve(2.117)
y
6
-
(
1
-
t
)
3
·
t
5
=
0
.
Another example (of the c
=
1
+
a type, see below) is(2.118)
A
(
z
)
=
1
3
·
d
ln
(
(
1
-
z
)
2
z
2
)
d
z
=
2
3
·
1
-
2
z
(
1
-
z
)
z
,
where the rational transformation(2.119)
R
(
z
)
=
z
·
(
z
-
2
)
3
(
1
-
2
z
)
3
=
-
8
z
-
36
z
2
-
126
z
3
-
387
z
4
+
⋯
verifies the “Rota-Baxter-like” functional relation (2.19). This example corresponds to the following covariance [12] on a Gauss hypergeometric integral:(2.120)
2
F
1
(
[
1
3
,
2
3
]
,
[
4
3
]
;
z
)
=
z
-
1
/
3
3
·
∫
0
z
t
-
2
/
3
(
1
-
t
)
-
2
/
3
·
d
t
=
1
2
·
2
-
z
1
-
2
z
·
2
F
1
(
[
1
3
,
2
3
]
,
[
4
3
]
;
z
·
(
z
-
2
)
3
(
1
-
2
z
)
3
)
,
which is associated with the elliptic curve(2.121)
y
3
-
(
1
-
t
)
2
·
t
2
=
0
.
Note that, similarly to the main example of the paper, there exist many rational transformations25 satisfying (2.19) that cannot be reduced to iterates of (2.119), for instance,(2.122)
T
(
z
)
=
-
27
·
z
·
(
1
-
z
)
(
z
2
-
z
+
1
)
3
(
z
3
+
3
z
2
-
6
z
+
1
)
3
=
-
27
z
-
378
z
2
-
3888
z
3
-
34074
z
4
-
271620
z
5
-
2032209
z
6
+
⋯
.
One verifies immediately that (2.122) actually verifies (2.19) with (2.118). Not surprisingly, the two rational transformations (2.119) and (2.122) commute.Another simple example with rational symmetries corresponds toΩ
=
(
D
z
+
A
(
z
)
)
·
D
z with(2.123)
A
(
z
)
=
-
1
2
·
3
z
-
1
z
(
1
-
z
)
=
1
2
·
d
ln
(
z
·
(
1
-
z
)
2
)
d
z
.
It has the simple (genus zero) hypergeometric solution26:(2.124)
ℱ
(
z
)
=
z
1
/
2
·
2
F
1
(
[
1
,
1
2
]
,
[
3
2
]
;
z
)
=
arctanh
(
z
1
/
2
)
.
The linear differential operator Ω is covariant under the change of variable z
→
1
/
z and z
→
R
(
z
), where27
(2.125)
R
(
z
)
=
4
z
(
1
+
z
)
2
.
One can easily check that (2.123) and (2.125) satisfy the functional equation (2.19). One also verifies that (2.123) and z
→
1
/
z or the iterates of (2.125) satisfy the functional equation (2.19). The solution of the adjoint operator are (
1
-
z
)
·
z
1
/
2 and(2.126)
F
(
z
)
=
z
·
(
1
-
z
)
·
2
F
1
(
[
1
,
1
2
]
,
[
3
2
]
;
z
)
=
z
1
/
2
·
(
1
-
z
)
·
arctanh
(
z
1
/
2
)
=
z
-
2
3
z
2
-
2
15
z
3
-
2
35
z
4
-
2
63
z
5
-
2
99
z
6
-
2
143
z
7
+
⋯
.
One verifies, again, that (2.126) and (2.125) commute, (2.126) corresponding to the “infinitesimal composition” of (2.125) (see (2.56)).A first natural generalization amounts to keeping the remarkable factorization (2.6) which will, in fact, reduce the covariance of a second-order operator to the covariance of a first-order operator.28 Such a situation occurs for Gauss hypergeometric functions2
F
1
(
[
a
,
b
]
,
[
1
+
a
]
;
z
) solution of the (
a
,
b
)-symmetric linear differential operator(2.127)
z
·
(
1
-
z
)
·
D
z
2
+
(
c
-
(
a
+
b
+
1
)
·
z
)
·
D
z
-
a
·
b
,
as soon as29
c
=
1
+
a. For instance(2.128)
ℱ
(
z
)
=
z
a
·
2
F
1
(
[
a
,
b
]
,
[
1
+
a
]
;
z
)
,
is an integral of a simple algebraic function and is solution with the constant function of the second-order operator(2.129)
Ω
=
(
D
z
+
(
a
-
b
-
1
)
z
+
1
-
a
z
·
(
1
-
z
)
)
·
D
z
=
(
D
z
+
d
ln
(
(
1
-
z
)
b
·
z
1
-
a
)
d
z
)
·
D
z
,
yielding a new A
(
z
):(2.130)
A
(
z
)
=
(
1
-
a
)
+
(
a
-
b
-
1
)
z
(
1
-
z
)
·
z
=
1
-
a
z
-
b
1
-
z
.
The adjoint of (2.129) has the simple solution z
1
-
a
·
(
1
-
z
)
b:(2.131)
F
(
z
)
=
z
·
(
1
-
z
)
b
·
2
F
1
(
[
a
,
b
]
,
[
1
+
a
]
;
z
)
.Due to the(
a
,
b
)-symmetry of (2.127) we have a similar result for c
=
1
+
b. The function ℱ
(
z
)
=
z
b
·
2
F
1
(
[
a
,
b
]
,
[
1
+
b
]
;
z
) is solution of (2.129) where a and b have been permuted:(2.132)
(
D
z
+
(
b
-
a
-
1
)
z
+
1
-
b
z
·
(
1
-
z
)
)
·
D
z
,
yielding another A
(
z
)
(2.133)
A
(
z
)
=
(
1
-
b
)
+
(
b
-
a
-
1
)
z
(
1
-
z
)
·
z
,
The adjoint of (2.132) has the solution (
1
-
z
)
a
·
z
1
-
b together with the hypergeometric function:(2.134)
F
(
z
)
=
z
·
(
1
-
z
)
a
·
2
F
1
(
[
a
,
b
]
,
[
1
+
b
]
;
z
)
,
where one recovers the previous result (2.126).We are seeking for (Gauss hypergeometric) second-order differential equations30 with an infinite number of (hopefully rational, if not algebraic) symmetries: this is another way to say that we are not looking for generic Gauss hypergeometric differential equations, but Gauss hypergeometric differential equations related to elliptic curves, and thus having an infinite set of such isogenies. We are necessarily in the framework where the two parameters a and b of the Gauss hypergeometric are rational numbers in order to have integral of algebraic functions (yielding globally nilpotent [11] second-order differential operators). Let us denote by D the common denominator of the two rational numbers a
=
N
a
/
D and b
=
N
b
/
D, function (2.128) is associated to a period of the algebraic curve(2.135)
y
D
=
(
1
-
t
)
N
b
·
t
D
-
N
a
.
We just need to restrict to triplets of integers (
N
a
,
N
b
,
D
) such that the previous curve is an elliptic curve.Let us give an example (of thec
=
1
+
b type) that does not correspond to a genus one curve, with(2.136)
2
F
1
(
[
1
3
,
1
6
]
,
[
7
6
]
;
z
)
=
1
6
·
z
-
1
/
6
·
∫
0
z
t
-
5
/
6
(
1
-
t
)
-
1
/
3
·
d
t
,
which corresponds to the genus two curve:(2.137)
y
6
-
(
1
-
t
)
2
·
t
5
=
0
.
Again one introduces A
(
z
)
(2.138)
A
(
z
)
=
1
6
·
d
ln
(
(
1
-
z
)
2
z
5
)
d
z
=
1
6
·
5
-
7
z
z
·
(
1
-
z
)
and seeks for R
(
z
) as series expansions analytical at z
=
0. One gets actually, order by order, a one-parameter family(2.139)
R
a
1
(
z
)
=
a
1
·
z
-
2
7
a
1
·
(
a
1
-
1
)
·
z
2
+
1
637
a
1
·
(
a
1
-
1
)
·
(
17
a
1
-
87
)
·
z
3
-
2
84721
a
1
·
(
a
1
-
1
)
·
(
113
a
1
2
-
856
a
1
+
3438
)
·
z
4
-
1
38548055
a
1
·
(
a
1
-
1
)
·
(
3674
a
1
3
+
121194
a
1
2
-
552261
a
1
+
2095059
)
·
z
5
+
⋯
+
1
+
ϵ
n
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
·
z
n
+
⋯
,
where ϵ
n
=
0 for n odd and ϵ
n
=
1 for n even, and N
(
n
) is a (large) integer depending on n, and P
n
(
a
1
) is a polynomial with integer coefficients of degree n
-
2. One easily verifies, order by order, that one gets a one-parameter family of transformations commuting for different values of the parameter:(2.140)
R
a
1
(
R
b
1
(
z
)
)
=
R
b
1
(
R
a
1
(
z
)
)
=
R
a
1
b
1
(
z
)
.As far as the “algorithmic complexity” of this series (2.139) is concerned it is worth noticing that the degree growth [18] of the series coefficients is actually linear and not exponential as we could expect [19] at first sight. Even if this series is transcendental, it is not a “wild” series.Seeking for selected values ofa
1 such that the previous series (2.139) reduces to a rational function one can try to reproduce the simple calculations (2.48), (2.49), but unfortunately “shooting in the dark” because we have no hint of a well-suited denominator (if any!) like the polynomials in the lhs of (2.48), (2.49).It is also worth noticing that if we slightly changeA
(
z
) into(2.141)
A
(
z
)
=
1
N
·
d
ln
(
(
1
-
z
)
2
z
5
)
d
z
=
1
N
·
5
-
7
z
z
·
(
1
-
z
)
,
the algebraic curve (2.137) becomes y
N
-
(
1
-
t
)
2
·
t
5
=
0 which has, for instance genus five for N
=
11, but genus zero for N
=
7. For any of these cases of (2.141) one can easily get, order by order, a one-parameter series R
a
1 totally similar to (2.139) with, again, polynomials P
n
(
a
1
) of degree n
-
2.The first coefficienta
2 is in general(2.142)
a
2
=
-
2
2
N
-
5
·
a
1
·
(
a
1
-
1
)
.For the genus zero case,N
=
7
(2.143)
a
2
=
-
2
9
·
a
1
·
(
a
1
-
1
)
,
a
3
=
-
1
1296
·
a
1
·
(
a
1
-
1
)
·
(
127
-
a
1
)
,
a
4
=
-
1
134136
·
a
1
·
(
a
1
-
1
)
·
(
254
a
1
2
+
185
a
1
+
7499
)
,
…
,
a
n
=
-
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
,
which corresponds to the solution(2.144)
2
7
·
∫
0
z
z
-
5
/
7
·
(
1
-
z
)
-
2
/
7
·
d
t
=
z
2
/
7
·
2
F
1
(
[
2
7
,
2
7
]
,
[
9
7
]
;
z
)
.
Using the parametrization of the genus zero curve
(2.145)
y
=
-
(
u
+
1
)
2
·
u
5
(
u
+
1
)
7
-
u
7
,
t
=
-
u
7
(
u
+
1
)
7
-
u
7
,
one can actually perform the integration (2.144) of d
t
/
y and get an alternative form of the hypergeometric function (2.144):(2.146)
∫
0
z
z
-
5
/
7
·
(
1
-
z
)
-
2
/
7
·
d
t
=
∫
0
u
ρ
(
u
)
·
d
u
=
∫
0
v
v
1
-
v
7
·
d
v
,
where
:
z
=
-
u
7
(
u
+
1
)
7
-
u
7
,
ρ
(
u
)
=
(
u
+
1
)
4
·
u
(
u
+
1
)
7
-
u
7
,
and
:
v
=
u
1
+
u
,
z
=
v
7
v
7
-
1
.
Except transformations like v
→
ω
·
v (with ω
7
=
1) which have no impact on z, it seems difficult to find rational symmetries in this genus zero case.ForN
=
11 (genus five) the first successive coefficients read:(2.147)
a
2
=
-
2
17
·
a
1
·
(
a
1
-
1
)
,
a
3
=
-
1
8092
·
a
1
·
(
a
1
-
1
)
·
(
143
a
1
+
367
)
,
a
4
=
-
1
206346
·
a
1
·
(
a
1
-
1
)
·
(
1186
a
1
2
+
2473
a
1
+
5011
)
,
…
,
a
n
=
-
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
.
The “infinitesimal composition” function F
(
z
) (see (2.56), (2.57), and (2.58)) reads,(2.148)
F
(
z
)
=
∂
R
a
1
∂
a
1
|
a
1
=
1
=
z
-
2
17
z
2
-
15
238
z
3
-
5
119
z
4
-
37
1190
z
5
-
888
36295
z
6
-
2183
108885
z
7
-
4366
258213
z
8
-
58941
4045337
z
9
-
1807524
141586795
z
10
-
46543743
4106017055
z
11
-
5305986702
521464165985
z
12
+
⋯
and again we can actually check that this is actually the series expansion of the hypergeometric function(2.149)
z
·
(
1
-
z
)
·
2
F
1
(
[
1
,
15
11
]
,
[
17
11
]
;
z
)
,
solution of Ω
* the adjoint of the Ω linear differential operator corresponding to this (genus 5) N
=
11 case:(2.150)
Ω
*
=
D
z
·
(
D
z
+
1
11
·
7
z
-
5
z
·
(
1
-
z
)
)
.We have similar results for (2.128), (2.129), (2.130). As far as these one-parameter families of transformations R
a
1 are concerned, the only difference between the generic cases corresponding to arbitrary genus and genus one cases like (2.118) is that in the generic higher genus case, only a finite number of values of the parameter a
1 can correspond to rational functions. Note that this higher genus result generalizes to the arbitrary genus Gauss hypergeometric functions (2.128) and associated operators (2.129) and function (2.130). In this general case one can also get order by order a one-parameter family of transformations R
a
1 satisfying a commutation relation (2.141).Note thatR
(
z
)
=
1
/
z is actually a solution of (2.19) for this genus-two example (2.139). Along this line of selected R
(
z
) solutions of (2.19) many interesting subcases of this general case (2.128), (2.129), (2.130) are given in Appendix C.1.In our previous genus-one examples, with this close identification between therenormalization group and the isogenies of elliptic curves, we saw that, in order to obtain linear differential operators covariant by an infinite number of transformations (rational or algebraic), we must restrict our second-order Gauss hypergeometric differential operator to Gauss hypergeometric associated to elliptic curves (see Appendices C and D). Beyond this framework we still have one-parameter families (see (2.141)) but we cannot expect an infinite number of rational (and probably algebraic) transformations to be particular cases of such families of transcendental transformations.


## 2.1. A Few Remarks on the “Rota-Baxter-Like” Functional Equation

The functional equation17(2.19) is the (necessary and sufficient) condition for Ω
=
(
D
z
+
A
(
z
)
)
·
D
z to be covariant by z
→
R
(
z
).Using the chain rule formula of derivatives of composed functions:(2.24)
d
R
(
R
(
z
)
)
d
z
=
d
R
(
z
)
d
z
·
[
d
R
(
z
)
d
z
(
R
(
z
)
)
]
,
d
2
R
(
R
(
z
)
)
d
z
2
=
d
2
R
(
z
)
d
z
2
·
[
d
R
(
z
)
d
z
(
R
(
z
)
)
]
+
(
d
R
(
z
)
d
z
)
2
·
[
d
2
R
(
z
)
d
z
2
(
R
(
z
)
)
]
,
one can show that, for A
(
z
) fixed, the “Rota-Baxter-like” functional equation (2.19) is invariant by the composition of R
(
z
) by itself R
(
z
)
→
R
(
R
(
z
)
), R
(
R
(
R
(
z
)
)
)
,
…
. This result can be generalized to any composition of various R
(
z
)'s satisfying (2.19). This is in agreement with the fact that (2.19) is the condition for Ω
=
(
D
z
+
A
(
z
)
)
·
D
z to be covariant by z
→
R
(
z
) it must be invariant by composition of R
(
z
)'s (for A
(
z
) fixed).Note that we have not used here the fact that for globally nilpotent [11] operators, A
(
z
) and B
(
z
) are necessarily log derivatives of Nth roots of rational functions. For R
(
z
)
=
-
4
z
/
(
1
-
z
)
2:(2.25)
A
(
z
)
=
1
4
·
d
ln
(
a
(
z
)
)
d
z
,
B
(
z
)
=
1
4
·
d
ln
(
b
(
z
)
)
d
z
,
a
(
z
)
=
(
1
-
z
)
2
·
z
3
,
b
(
z
)
=
z
3
·
(
1
+
z
)
4
(
1
-
z
)
10
.
The existence of the underlying a
(
z
) in (2.25) consequence of a global nilpotence of the order-one differential operator, can however be seen in the following remark on the zeros of the lhs and rhs terms in the functional equation (2.19). When R
(
z
) is a rational function (e.g., -
4
z
/
(
1
-
z
)
2 or any of its iterates R
(
n
)
(
z
)), the lhs and rhs of (2.19) are rational expressions. The zeros are roots of the numerators of these rational expressions. Because of (2.25) the functional equation (2.19) can be rewritten (after dividing by R
'
(
z
)) as(2.26)
(
d
R
(
z
)
d
z
)
·
A
(
R
(
z
)
)
=
A
(
z
)
+
d
d
z
(
ln
(
d
R
(
z
)
d
z
)
)
=
1
4
·
d
d
z
(
ln
(
a
(
z
)
·
(
d
R
(
z
)
d
z
)
4
)
)
.
One easily verifies, in our example, that the zeros of the rhs of (2.26) come from the zeros of A
(
R
(
z
)
) (and not from the zeros of R
'
(
z
) in the lhs of (2.26)). The zeros of the log-derivative rhs of (2.26) correspond to a
(
z
)
·
R
'
(
z
)
4
=
ρ, where ρ is a constant to be found. Let us consider for R
(
z
) the nth iterates of -
4
z
/
(
1
-
z
)
2 that we denote R
(
n
)
(
z
). A straightforward calculation shows that the zeros of A
(
R
(
n
)
(
z
)
) or a
'
(
R
(
n
)
(
z
)
) (where a
'
(
z
) denotes the derivative of a
(
z
) namely, (
z
-
1
)
(
5
z
-
3
)
·
z
2) actually correspond to the general closed formula:(2.27)
5
5
·
a
(
z
)
·
(
d
R
(
n
)
(
z
)
d
z
)
4
-
4
·
3
3
·
(
-
4
)
n
=
0
.
More precisely the zeros of 5
·
R
(
n
)
(
z
)
-
3 verify (2.27), or, in other words, the numerator of 5
R
(
n
)
(
z
)
-
3 divides the numerator of the lhs of (2.27).In another case forT
(
z
) given by (2.45), which also verifies (2.19) (see below), the relation (2.27) is replaced by(2.28)
5
5
·
a
(
z
)
·
(
d
T
(
n
)
(
z
)
d
z
)
4
-
4
·
3
3
·
(
-
7
-
24
i
)
n
=
0
.
More generally for a rational function ρ
(
x
), obtained by an arbitrary composition of -
4
z
/
(
1
-
z
)
2 and T
(
z
), we would have(2.29)
5
5
·
a
(
z
)
·
(
d
ρ
(
z
)
d
z
)
4
-
4
·
3
3
·
λ
n
=
0
.
where λ corresponds to(2.30)
ρ
(
x
)
=
λ
·
z
+
⋯
,
λ
=
[
d
ρ
(
z
)
d
z
]
z
=
0
.


## 2.2. Symmetries ofΩ, Solutions to the “Rota-Baxter-Like” Functional Equation

Let us now analyse all the symmetries of the linear differential operatorΩ
=
(
D
z
+
A
(
z
)
)
·
D
z by analyzing all the solutions of (2.19) for a given A
(
z
). For simplicity we will restrict to A
(
z
)
=
(
3
-
5
z
)
/
z
/
(
1
-
z
)
/
4 which corresponds to R
(
z
)
=
-
4
z
/
(
z
-
1
)
2 and all its iterates (2.9). Let us first seek for other (more general) solutions that are analytic at
z
=
0:(2.31)
R
(
z
)
=
a
1
·
z
+
a
2
·
z
2
+
a
3
·
z
3
+
⋯
.
It is a straightforward calculation to get, order by order from (2.19), the successive coefficients a
n in (2.31) as polynomial expressions (with rational coefficients) of the first coefficient a
1 with(2.32)
a
2
=
-
2
5
·
a
1
·
(
a
1
-
1
)
,
a
3
=
1
75
·
a
1
·
(
a
1
-
1
)
·
(
7
a
1
-
17
)
,
a
4
=
-
2
4875
·
a
1
·
(
a
1
-
1
)
·
(
41
a
1
2
-
232
a
1
+
366
)
,
…
,
a
n
=
-
n
5
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
P
n
(
-
4
)
,
where P
n
(
a
1
) is a polynomial with integer coefficients of degree n
-
2. Since we have here a series depending on one parameter a
1 we will denote it R
a
1
(
z
). This is a quite remarkable series depending on one parameter.18 One can easily verify that this series actually reduces (as it should!) to the successive iterates (2.9) of -
4
z
/
(
1
-
z
)
2 for a
1
=
(
-
4
)
n. In other words this one-parameter family of “functions” actually reduces to rational functions for an infinite number of integer values a
1
=
(
-
4
)
n.Furthermore, one can also verify a quite essential property we expect for a representation of the renormalization group, namely, that twoR
a
1
(
z
) for different values of a
1 commute, the result corresponding to the product of these two a
1:(2.33)
R
a
1
(
R
b
1
(
z
)
)
=
R
b
1
(
R
a
1
(
z
)
)
=
R
a
1
·
b
1
(
z
)
.
The neutral element must necessarily correspond to a
1
=
1 which is actually the identity transformation R
1
(
z
)
=
z. We have an “absorbing” element corresponding to a
1
=
0, namely, R
0
(
z
)
=
0. Performing the inverse of R
a
1
(
z
) (with respect to the composition of functions) amounts to changing a
1 into its inverse 1
/
a
1. Let us explore some “reversibility” property of our exact representation of a renormalization group with the inverse of the rational transformations (2.9). The inverse of R
-
4
(
z
)
=
-
4
z
/
(
1
-
z
)
2 must correspond to a
1
=
-
1
/
4:(2.34)
R
-
1
/
4
(
z
)
=
-
1
4
·
z
-
1
8
z
2
-
5
64
z
3
-
7
128
z
4
-
21
512
z
5
+
⋯
.
However, a straight calculation of the inverse of R
-
4
(
z
)
=
-
4
z
/
(
1
-
z
)
2 gives a multivalued function, or if one prefers, two functions(2.35)
S
-
1
/
4
(
1
)
(
z
)
=
z
-
2
+
2
1
-
z
z
=
-
1
4
·
z
-
1
8
z
2
+
⋯
,
S
-
1
/
4
(
2
)
(
z
)
=
z
-
2
-
2
1
-
z
z
=
-
4
z
+
2
+
1
4
z
+
1
8
z
2
+
⋯
,
which are the two roots of the simple quadratic relation (R
-
4
(
z
'
)
=
z):(2.36)
z
'
2
-
2
·
(
1
-
2
z
)
·
z
'
+
1
=
0
,
where it is clear that the product of these two functions is equal to +
1. The radius of convergence of S
-
1
/
4
(
1
)
(
z
) is 1.Because of our choice to seek for functions analytical atz
=
0 our renormalization group representation “chooses” the unique root that is analytical at z
=
0, namely, S
-
1
/
4
(
1
)
(
z
). For the next iterate of R
-
4
(
z
)
=
-
4
z
/
(
1
-
z
)
2 in (2.9) the inverse transformation corresponds to the roots of the polynomial equation of degree four (R
16
(
z
'
)
=
z):(2.37)
z
'
4
+
(
4
-
16
z
)
·
z
′
3
+
(
6
+
32
z
)
·
z
′
2
+
(
4
-
16
z
)
·
z
'
+
1
=
0
,
which yields four roots, one of which is analytical at z
=
0 and corresponds to a
1
=
1
/
(
-
4
)
2 in our one-parameter family of (renormalization) transformations:(2.38)
S
1
/
16
(
1
)
(
z
)
=
1
16
z
+
3
128
z
2
+
53
4096
z
3
+
277
32768
z
4
+
3181
524288
z
5
+
⋯
,
its (multiplicative) inverse S
1
/
16
(
2
)
(
z
)
=
1
/
S
1
/
16
(
1
)
(
z
):(2.39)
S
1
/
16
(
2
)
(
z
)
=
16
z
-
6
-
17
16
z
-
67
128
z
2
-
1333
4096
z
3
-
7445
32768
z
4
+
⋯
,
and two (formal) Puiseux series (u
=
±
z):(2.40)
S
1
/
16
(
3
)
(
z
)
=
1
+
u
+
1
2
u
2
+
3
8
u
3
+
1
4
u
4
+
27
128
u
5
+
5
32
u
6
+
⋯
.
Many of these results are better understood when one keeps in mind that there is a special transformation J
:
z
↔
1
/
z which is also a R-solution of (2.19) and verifies many compatibility relations with these transformations (I
d denotes the identity transformation R
0
(
z
)):(2.41)
R
-
4
·
J
=
R
-
4
,
S
-
1
/
4
(
2
)
·
R
-
4
=
J
,
R
-
4
·
S
-
1
/
4
(
1
)
=
S
-
1
/
4
(
1
)
·
R
-
4
=
I
d
,
S
1
/
16
(
1
)
(
z
)
=
S
-
1
/
4
(
1
)
·
S
-
1
/
4
(
1
)
,
S
1
/
16
(
2
)
(
z
)
=
S
-
1
/
4
(
1
)
·
S
-
1
/
4
(
2
)
,
J
·
S
-
1
/
4
(
1
)
=
S
-
1
/
4
(
2
)
,
J
·
S
-
1
/
4
(
2
)
=
S
-
1
/
4
(
1
)
,
…
,
where the dot corresponds, here, to the composition of functions. These symmetries of the linear differential operator Ω correspond to isogenies of the elliptic curve (2.2).It is clear that we have another one-parameter family corresponding toJ
·
R
a
1 with an expansion of the form(2.42)
J
·
R
a
1
=
b
1
z
-
2
5
·
(
b
1
-
1
)
-
1
15
·
b
1
2
-
1
b
1
·
z
-
2
975
·
(
b
1
-
1
)
(
4
b
1
+
1
)
(
4
b
1
+
3
)
b
1
2
·
z
2
-
1
248625
·
(
b
1
-
1
)
(
4
b
1
+
1
)
(
1268
b
1
2
+
951
b
1
+
91
)
b
1
3
·
z
3
-
2
2071875
·
(
b
1
-
1
)
(
4
b
1
+
1
)
(
3688
b
1
3
+
2766
b
1
2
+
404
b
1
+
17
)
b
1
4
·
z
4
+
⋯
.
For b
1
=
-
1
/
4, b
1
=
(
-
1
/
4
)
2, b
1
=
(
-
1
/
4
)
3, this family reduces to the (multiplicative) inverse of the successive rational functions displayed in (2.9)(2.43)
-
1
4
·
(
1
-
z
)
2
z
→
1
16
·
(
1
+
z
)
4
(
1
-
z
)
2
·
z
→
-
1
64
·
(
1
-
6
z
+
z
2
)
4
(
1
-
z
)
2
(
1
+
z
)
4
·
z
→
⋯
,
which can also be written as:
(2.44)
-
1
4
·
(
z
+
1
z
)
+
1
2
,
1
16
·
(
z
+
1
z
)
+
3
8
+
z
(
1
-
z
)
2
,
-
1
64
·
(
z
+
1
z
)
+
13
32
-
z
4
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
,
1
256
·
(
z
+
1
z
)
+
51
128
+
z
16
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
+
z
(
-
4
)
n
-
6
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
,
-
1
1024
·
(
z
+
1
z
)
+
205
512
-
z
164
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
-
4
z
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
-
64
z
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
(
1
+
20
z
-
26
z
2
+
20
z
3
+
z
4
)
4
,
…
,
1
(
-
4
)
n
·
(
z
+
1
z
)
+
2
5
4
n
(
4
n
-
(
-
1
)
n
)
+
z
(
-
4
)
n
-
2
·
17
-
60
z
+
102
z
2
-
60
z
3
+
17
z
4
(
1
-
z
)
2
(
1
+
z
)
4
+
z
(
-
4
)
n
-
6
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
+
z
(
-
4
)
n
-
8
·
(
1
-
z
)
2
(
1
+
z
)
4
(
z
2
-
6
z
+
1
)
4
(
1
+
20
z
-
26
z
2
+
20
z
3
+
z
4
)
4
+
⋯
,
where we discover some “additive structure” of these successive rational functions.In fact, due to the specificity of this elliptic curve (occurrence of complex multiplication), we have another remarkable rational transformation solution of (2.19), preserving covariantly Ω. Let us introduce the rational transformation (i denotes -
1):(2.45)
T
(
z
)
=
z
·
(
z
-
(
1
+
2
i
)
1
-
(
1
+
2
i
)
·
z
)
4
,
we also have the remarkable covariance [12]:(2.46)
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
=
1
-
z
/
(
1
+
2
i
)
1
-
(
1
+
2
i
)
z
·
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
T
(
z
)
)
,
which can be rewritten in a simpler way on (2.3) (see (2.4)).It is a straightforward matter to see thatT
(
z
) actually belongs to the R
a
1
(
z
) one-parameter family:(2.47)
T
(
z
)
=
R
a
1
(
z
)
=
-
(
7
+
24
i
)
·
z
+
⋯
,
a
1
=
-
25
·
ρ
,
ρ
=
(
7
+
24
i
)
25
,
|
ρ
|
=
1
.As far as the reduction of (2.32) to a rational function is concerned, it is straightforward to see that:(2.48)
(
1
-
z
)
2
·
(
1
+
z
)
4
·
R
a
1
(
z
)
=
a
1
·
z
+
⋯
-
2
175746796875
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
4
)
·
(
a
1
-
16
)
·
P
8
(
a
1
)
·
z
8
+
⋯
-
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
4
)
·
(
a
1
-
16
)
·
P
n
(
a
1
)
·
z
n
+
⋯
,
where N
(
n
) is a large integer growing with n, and P
n is a polynomial with integer coefficients of degree n
-
4, or(2.49)
(
1
-
(
1
+
2
i
)
·
z
)
4
·
R
a
1
(
z
)
=
a
1
·
z
+
⋯
−
4
1243125
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
7
+
24
i
)
·
(
P
6
(
a
1
)
+
i
Q
6
(
a
1
)
)
·
z
6
+
⋯
+
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
(
a
1
+
7
+
24
i
)
·
(
P
n
(
a
1
)
+
i
Q
n
(
a
1
)
)
·
z
n
+
⋯
,
where P
n and Q
n are two polynomials with integer coefficients of degree, respectively, n
-
3 and n
-
4.Similar calculations can be performed forT
*
(
z
) defined by(2.50)
T
*
(
z
)
=
z
·
(
z
-
(
1
-
2
i
)
(
1
-
2
i
)
z
-
1
)
4
,
for which we also have the covariance(2.51)
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
=
1
-
z
/
(
1
-
2
i
)
1
-
(
1
-
2
i
)
z
·
2
F
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
T
*
(
z
)
)
.It is a simple calculation to check that any iterate ofT
(
z
) (resp. T
*
(
z
)) is actually a solution of (2.19) and corresponds to R
a
1
(
z
) for the infinite number of values a
1
=
(
-
7
-
24
i
)
N (resp. (
-
7
+
24
i
)
N). Furthermore, one verifies, as it should (see (2.33)), that the three rational functions R
-
4
(
z
), T
(
z
), and T
*
(
z
) commute. It is also a straightforward calculation to see that the rational function built from any composition of R
-
4
(
z
), T
(
z
), and T
*
(
z
) is actually a solution of (2.19). We thus have a triple infinity of values of a
1, namely a
1
=
(
-
4
)
M
·
(
-
7
-
24
i
)
N
·
(
-
7
+
24
i
)
P
for any integer
M
,
N
and
P, for which R
a
1
(
z
) reduces to rational functions. We are in fact describing (some subset of) the isogenies of the elliptic curve (2.2), and identifying these isogenies with a discrete subset of the renormalization group. Conversely, a functional equation like (2.19) can be seen as a way to extend the n-fold composition of a rational function R
(
z
) (namely R
(
R
(
⋯
R
(
z
)
⋯
)
)) to n
any complex number.


## 2.3. Revisiting the One-Parameter Family of Solutions of the “Rota-Baxter-Like” Functional Equation

This extension can be revisited as follows. Keeping in mind the well-known example of the parametrization of the standard mapz
→
4
z
·
(
1
-
z
) with z
=
sin
2
(
θ
), yielding θ
→
2
θ, let us seek for a (transcendental) parametrization z
=
P
(
u
) such that(2.52)
R
-
4
(
P
(
u
)
)
=
P
(
-
4
u
)
or
:
R
-
4
=
P
·
H
-
4
·
P
-
1
,
where H
a
1 denotes the scaling transformation z
→
a
1
·
z (here H
-
4
:
z
→
-
4
·
z) and P
-
1 denotes the inverse transformation of P (for the composition). One can easily find such a (transcendental) parametrization order by order(2.53)
P
(
z
)
=
z
-
2
5
z
2
+
7
75
z
3
-
82
4875
z
4
+
1078
414375
z
5
−
452
1243125
z
6
+
57311
1212046875
z
7
-
1023946
175746796875
z
8
+
⋯
,
and similarly for its inverse (for the composition) transformation(2.54)
Q
(
z
)
=
P
-
1
(
z
)
=
z
+
2
5
z
2
+
17
75
z
3
+
244
1625
z
4
+
45043
414375
z
5
+
2302
27625
z
6
+
128941
1939275
z
7
+
15365176
281194875
z
8
+
⋯
.
This approach is reminiscent of the conjugation introduced in Siegel's theorem [15–17]. It is a straightforward matter to see (order by order) that one actually has(2.55)
R
a
1
(
P
(
u
)
)
=
P
(
a
1
·
u
)
or
:
R
a
1
=
P
·
H
a
1
·
P
-
1
.
The structure of the (one-parameter) renormalization group and the extension of the composition of n times a rational function R
(
z
) (namely, R
(
R
(
⋯
R
(
z
)
⋯
)
)) to n
any complex number become a straight consequence of this relation. Along this line one can define some “infinitesimal composition” (ϵ
≃
0):(2.56)
R
1
+
ϵ
(
z
)
=
P
·
H
1
+
ϵ
·
P
-
1
(
z
)
=
z
+
ϵ
·
F
(
z
)
+
⋯
,
where one can find, order by order, the “infinitesimal composition” function F
(
z
):(2.57)
F
(
z
)
=
z
-
2
5
z
2
-
2
15
z
3
-
14
195
z
4
-
154
3315
z
5
-
22
663
z
6
-
418
16575
z
7
-
9614
480675
z
8
-
2622
160225
z
9
+
⋯
.
It is straightforward to see, from (2.33), that the function F
(
z
) satisfies the following functional equations involving a rational function R
(
z
) (in the one-parameter family R
a
1
(
z
)):(2.58)
d
R
(
z
)
d
z
·
F
(
z
)
=
F
(
R
(
z
)
)
,
d
R
(
n
)
(
z
)
d
z
·
F
(
z
)
=
F
(
R
(
n
)
(
z
)
)
,
where
:
R
(
n
)
(
z
)
=
R
(
R
(
⋯
R
(
z
)
)
⋯
)
.
F
(
z
) cannot be a rational or algebraic function. Let us consider the fixed points of R
(
n
)
(
z
). Generically d
R
(
n
)
(
z
)
/
d
z is not equal to 0 or ∞ at any of these fixed points. Therefore one must have F
(
z
)
=
0 or F
(
z
)
=
∞ for the infinite set of these fixed points: F
(
z
) cannot be a rational or algebraic function, it is a transcendental function, and similarly for the parametrization function P
(
z
). In fact, let us introduce the function(2.59)
G
(
z
)
=
(
1
-
z
)
·
F
(
z
)
,
G
(
z
)
=
z
-
7
5
z
2
+
4
15
z
3
+
4
65
z
4
+
28
1105
z
5
+
44
3315
z
6
+
44
5525
z
7
+
836
160225
z
8
+
1748
480675
z
9
+
⋯
+
g
n
·
z
n
+
⋯
.
One actually finds that the successive g
n satisfies the very simple (hypergeometric function) relation:(2.60)
g
n
+
1
g
n
=
4
n
-
9
4
n
+
1
.
The function G
(
z
) is actually the hypergeometric function solution of the homogeneous operator(2.61)
D
z
2
+
1
4
13
z
-
3
z
·
(
1
-
z
)
·
D
z
+
3
4
6
z
2
-
3
z
+
1
(
1
-
z
)
2
·
z
2
,
or of the inhomogeneous ODE(2.62)
4
z
·
(
1
-
z
)
·
d
G
(
z
)
d
z
+
(
9
z
-
3
)
·
G
(
z
)
-
z
·
(
1
-
z
)
2
=
0
.
One deduces the expression of F
(
z
) as a hypergeometric function(2.63)
F
(
z
)
=
z
·
(
1
-
z
)
1
/
2
·
F
2
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
z
)
=
∂
R
a
1
∂
a
1
|
a
1
=
1
.Finally we get the linear differential operator annihilatingF
(
z
)
(2.64)
Ω
F
=
D
z
2
+
1
4
·
5
z
-
3
z
(
1
-
z
)
·
D
z
+
1
4
·
3
-
6
z
+
5
z
2
(
1
-
z
)
2
z
2
=
D
z
·
(
D
z
-
1
4
·
3
-
5
z
z
·
(
1
-
z
)
)
,
which is, in fact, nothing but Ω
* the adjoint of linear differential operator Ω (see (2.6)). One easily checks19 that the second-order differential equation Ω
F
(
y
(
z
)
)
=
0 transforms under the change of variable z
→
-
4
z
/
(
1
-
z
)
2 into the second-order differential equation Ω
F
(
R
)
(
y
(
z
)
)
=
0 with Ω
F
(
R
)
=
C
(
z
)
2
·
ω
F
(
R
) where the unitary (monic) operator ω
F
(
R
) is the conjugate of Ω
F:(2.65)
ω
F
(
R
)
=
D
z
2
-
1
4
·
11
z
2
+
30
z
+
3
z
·
(
1
-
z
)
(
1
+
z
)
·
D
z
+
1
4
·
3
+
12
z
+
50
z
2
+
12
z
3
+
3
z
4
z
2
·
(
1
-
z
)
2
(
1
+
z
)
2
=
(
1
C
(
z
)
)
·
D
z
·
(
D
z
-
1
4
·
3
-
5
z
z
·
(
1
-
z
)
)
·
C
(
z
)
=
(
1
C
(
z
)
)
·
Ω
F
·
C
(
z
)
=
(
1
C
(
z
)
)
·
Ω
*
·
C
(
z
)
with C
(
z
)
=
1
/
R
'
(
z
) and the “dot” denotes the composition of operators. Actually, the factors in the adjoint Ω
* transform under the change of variable z
→
-
4
z
/
(
1
-
z
)
2 as follows20:(2.66)
D
z
→
C
(
z
)
·
D
z
,
ω
1
*
=
→
(
ω
1
*
)
(
R
)
=
ω
1
*
·
C
(
z
)
,
Ω
*
→
Ω
F
(
R
)
=
C
(
z
)
·
Ω
*
·
C
(
z
)
which is precisely the transformation we need to match with (2.58) and see the ODE Ω
*
(
F
(
z
)
)
=
0 compatible with the change of variable z
→
-
4
z
/
(
1
-
z
)
2:(2.67)
Ω
*
(
F
(
z
)
)
=
0
→
(
C
(
z
)
·
Ω
*
·
C
(
z
)
)
(
F
(
R
(
z
)
)
)
=
(
C
(
z
)
·
Ω
*
·
C
(
z
)
)
(
R
′
(
z
)
·
F
(
z
)
)
=
C
(
z
)
·
Ω
*
(
F
(
z
)
)
=
0
.
This is, in fact, a quite general result that will be seen to be valid in a more general (higher genus) framework (see (2.148), (2.150) in what follows).Not surprisingly one can deduce from (2.33) and the previous results, in particular (2.63), the following results for R
a
1
(
z
):(2.68)
-
4
·
∂
R
a
1
∂
a
1
|
a
1
=
-
4
=
F
(
R
(
z
)
)
,
(
-
4
)
n
·
∂
R
a
1
∂
a
1
|
a
1
=
(
-
4
)
n
=
F
(
R
(
n
)
(
z
)
)
,
where R
(
z
)
=
-
4
z
/
(
1
-
z
)
2 and R
(
n
)
(
z
) denotes R
(
R
(
⋯
R
(
R
(
z
)
)
)
). Of course we have similar relation for T
(
z
), -
4 being replaced by -
7
-
24
i. Therefore the partial derivative ∂
R
a
1
/
∂
a
1 that can be expressed in terms of hypergeometric functions for for a double infinity of values of a
1, namely, a
1
=
(
-
4
)
M
×
(
-
7
-
24
i
)
N.One can, of course, check, order by order, that (2.58) is actually verified for any function in the one-parameter family R
a
1
(
z
):(2.69)
d
R
a
1
(
z
)
d
z
·
F
(
z
)
=
F
(
R
a
1
(
z
)
)
,
which corresponds to an infinitesimal version of (2.33).From (2.56) one simply deduces(2.70)
z
·
d
P
(
z
)
d
z
=
F
(
P
(
z
)
)
,
that we can check, order by order from (2.53), the series expansion of P
(
z
), and from (2.57) the series expansion of F
(
z
), but also(2.71)
d
Q
(
z
)
d
z
·
F
(
z
)
=
Q
(
z
)
that we can, check order by order, from (2.54), the series expansion of Q
(
z
)
=
P
-
1
(
z
) and from (2.57). We now deduce that the log-derivative of the “well-suited change of variable” Q
(
z
) is nothing but the (multiplicative) inverse of a hypergeometric function F
(
z
):(2.72)
d
ln
(
Q
(
z
)
)
d
z
=
1
F
(
z
)
,
Q
(
z
)
=
λ
·
exp
(
∫
z
d
z
F
(
z
)
)
.
The function Q
(
z
) is solution of the nonlinear differential equation(2.73)
-
4
z
2
·
(
1
-
z
)
2
·
(
Q
·
Q
(
1
)
·
Q
(
3
)
+
(
Q
(
1
)
)
2
·
Q
(
2
)
-
2
Q
·
(
Q
(
2
)
)
2
)
+
z
·
(
3
-
5
z
)
(
1
-
z
)
·
Q
(
1
)
·
(
Q
·
Q
(
2
)
-
(
Q
(
1
)
)
2
)
+
(
5
z
2
-
6
z
+
3
)
·
Q
·
(
Q
(
1
)
)
2
=
0
,
where the Q
(
n
)'s denote the nth derivative of Q
(
z
). At first sight Q
(
z
) would be a nonholonomic function, however, remarkably, it is a holonomic function solution of an order-five operator which factorizes as follows:(2.74)
Ω
Q
=
(
D
z
+
3
-
5
z
(
1
-
z
)
·
z
)
·
(
D
z
+
3
4
·
3
-
5
z
(
1
-
z
)
·
z
)
·
×
(
D
z
+
2
4
·
3
-
5
z
(
1
-
z
)
·
z
)
·
(
D
z
+
1
4
·
3
-
5
z
(
1
-
z
)
·
z
)
·
D
z
,
yielding the exact expression of Q
(
z
) in terms of hypergeometric functions:(2.75)
Q
(
z
)
=
z
·
(
F
2
1
(
z
1
-
z
[
1
2
,
1
4
]
,
[
5
4
]
;
z
z
1
-
z
)
)
4
=
z
1
-
z
·
(
2
F
1
(
[
1
4
,
3
4
]
,
[
5
4
]
;
-
z
1
-
z
)
)
4
,
that is, the fourth power of (2.3), with the differential operator (2.74) being the symmetric fourth power of Ω. From (2.3) we immediately get the covariance of Q
(
z
):(2.76)
Q
(
-
4
z
(
1
-
z
)
2
)
=
-
4
·
Q
(
z
)
,
and, more generally, Q
(
R
a
1
)
=
a
1
·
Q
(
z
). Since Q
(
z
) and F
(
z
) are expressed in terms of the same hypergeometric function, the relation (2.71) must be an identity on that hypergeometric function. This is actually the case. This hypergeometric function verifies the inhomogeneous equation:(2.77)
4
·
z
·
d
ℋ
(
z
)
d
z
+
ℋ
(
z
)
-
(
1
-
z
)
-
1
/
2
=
0
,
where(2.78)
ℋ
(
z
)
=
F
2
1
(
[
1
2
,
1
4
]
,
[
5
4
]
;
z
)
.RecallingQ
(
P
(
z
)
)
=
z, one has the following functional relation on P
(
z
):(2.79)
P
(
z
)
·
F
2
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
P
(
z
)
)
4
=
z
.Noting thatQ
(
z
4
)
1
/
4
=
ℱ
(
z
4
) (see (2.3)) can be expressed in term of an incomplete elliptic integral of the first kind of argument -
1
(2.80)
z
·
2
F
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
z
4
)
=
Elliptic
F
(
z
,
-
1
)
,
one can find that (2.79) rewrites on P
(
z
) as(2.81)
Elliptic
F
(
P
(
z
)
1
/
4
,
-
1
)
=
z
1
/
4
,
from which we deduce that the function P
(
z
) is nothing but aJacobi elliptic function
21
(2.82)
P
(
z
)
=
(
sn
(
z
1
/
4
,
-
1
)
)
4
.
In Appendix B we display a set of “Painlevé-like” ODEs22 verified by P
(
z
). From the simple nonlinear ODE on the Jacobi elliptic sinus, namely, S
′
′
+
2
·
S
3
=
0, and the exact expression of P
(
z
) in term of Jacobi elliptic sinus, one can deduce other nonlinear ODEs verified by the nonholonomic function P
(
z
) (P
(
1
)
=
d
P
(
z
)
/
d
z, P
(
2
)
=
d
2
P
(
z
)
/
d
z
2):(2.83)
z
3
/
2
·
(
P
(
1
)
)
2
-
(
1
-
P
)
·
P
3
/
2
=
0
,
(2.84)
P
(
2
)
-
3
4
·
(
P
(
1
)
)
2
P
+
3
4
·
P
(
1
)
z
+
1
2
·
P
3
/
2
z
3
/
2
=
0
.


## 2.4. Singularities of the Jacobi Elliptic FunctionP
(
z
)

Most of the results of this section, and to some extent, of the next one, are straight consequences of the exact closed expression ofP
(
z
) in terms of an elliptic function. Following the pedagogical approach of this paper we will rather follow a heuristic approach not taking into account the exact result (2.82), to display simple methods and ideas that can be used beyond exact results on a specific example.From a diff-Padé analysis of the series expansion ofP
(
z
), we got the sixty (closest to z
=
0) singularities. In particular we got that P
(
z
) has a radius of convergence R
≃
11.81704500807
⋯ corresponding to the following (closest to z
=
0) singularity z
=
z
s of P
(
z
):(2.85)
z
s
=
-
11.817045008077115768316337283432582087420697
⋯
=
(
-
4
)
·
2
F
1
(
[
1
4
,
1
2
]
,
[
5
4
]
;
1
)
4
=
-
1
16
·
π
6
Γ
(
3
/
4
)
8
.
This singularity corresponds to a pole of order four: P
(
z
)
≃
(
z
-
z
s
)
-
4. The function P
(
z
) has many other singularities:(2.86)
3
4
·
z
s
,
(
161
±
240
i
)
·
z
s
,
(
-
7
±
24
i
)
·
z
s
,
(
-
119
±
120
i
)
·
z
s
,
⋯
5
4
·
z
s
,
(
41
±
840
i
)
·
z
s
,
(
-
527
±
336
i
)
·
z
s
,
(
-
1519
±
720
i
)
·
z
s
,
⋯
7
4
·
z
s
,
(
1241
±
2520
i
)
·
z
s
,
(
-
567
±
1944
i
)
·
z
s
,
(
-
3479
±
1320
i
)
·
z
s
,
⋯In fact, introducingx and y the real and imaginary part of these singularities in z
s units, one finds out that they correspond to the double infinity of points(2.87)
x
=
(
m
1
2
-
2
m
1
m
2
-
m
2
2
)
·
(
m
1
2
+
2
m
1
m
2
-
m
2
2
)
,
y
=
4
·
m
1
m
2
·
(
m
2
-
m
1
)
·
(
m
2
+
m
1
)
,
where m
1 and m
2 are two integers, and they all lie on the intersection of an infinite number of genus zero curves indexed by the fourth power of an integer M
=
m
4 (m
=
m
1 or m
=
m
2):(2.88)
2
12
M
4
-
2
11
·
x
·
M
3
-
2
7
·
(
17
y
2
+
14
x
2
)
·
M
2
-
2
5
·
x
·
(
8
x
2
+
7
y
2
)
·
M
+
y
4
=
0
.
The parametrization (2.87) describes not only the poles of P
(
z
) when m
1
+
m
2 is odd, but also the zeros of P
(
z
) when m
1
+
m
2 is even. This (infinite) proliferation of singularities confirms the nonholonomic character of P
(
z
).These results are simply inherited from (2.82). The zeros and poles of the elliptic sinus s
n
(
z
,
i
) correspond to two lattice of periods. Denoting 𝒦
1 and 𝒦
2 the two periods of the elliptic curve, the location of the poles and zeros reads, respectively,(2.89)
P
n
1
,
n
2
=
2
n
1
·
𝒦
1
+
(
2
n
2
+
1
)
·
𝒦
2
,
Z
n
1
,
n
2
=
2
n
1
·
𝒦
1
+
2
n
2
·
𝒦
2
,
𝒦
1
=
π
3
/
2
2
3
/
2
·
1
Γ
(
3
/
4
)
2
,
𝒦
2
=
(
1
-
-
1
)
·
K
1
making crystal clear the fact that we have complex multiplication for this elliptic curve. The formula (2.87) just amount to saying that the poles and zeros of s
n
(
z
1
/
4
,
i
) are located at P
n
1
,
n
2
4 and Z
n
1
,
n
2
4:(2.90)
P
n
1
,
n
2
4
=
-
z
s
4
·
(
(
2
n
1
+
2
n
2
+
1
)
+
i
·
(
2
n
2
+
1
)
)
4
,
Z
n
1
,
n
2
4
=
-
z
s
4
·
(
(
2
n
1
+
2
n
2
+
1
)
+
i
·
2
n
2
)
4
.
The correspondence with (2.87) is m
1
=
n
1
+
2
n
2
+
1, m
2
=
-
n
1 for the poles and m
1
=
n
1
+
2
n
2, m
2
=
-
n
1 for the zeros.Remark 2.2.
let us consider thea
1
→
∞ limit of the one-parameter series R
a
1 (see (2.31), (2.32)) rewriting R
a
1
(
z
) as R
̃
b
1
(
u
)
(2.91)
R
̃
b
1
(
u
)
=
R
a
1
(
z
)
,
with
:
b
1
=
1
a
1
,
u
=
z
b
1
.
In the a
1
→
∞ limit, that is the b
1
→
0 limit, one easily verifies, order by order in u, that R
̃
b
1
(
u
) becomes exactly the transcendental parametrization function (2.53):
(2.92)
R
̃
b
1
(
u
)
→
P
(
u
)
when
b
1
→
0
.
For a
1
=
(
-
4
)
n (n
→
∞), one finds that the radius of convergence of the R
a
1
(
z
) series becomes in the n
→
∞ limit R
n
≃
z
s
/
4
n, in agreement23 with (2.91).


## 2.5.P
(
z
) and an Infinite Number of Rational Transformations: The Sky Is the Limit

Note that some nonlinear ODEs associated withP
(
z
) and displayed in Appendix B, namely (B.3) and (B.10), and the functional equation (2.79), are invariant by the change of variable (
P
(
z
)
,
z
)
→
(
-
4
P
(
z
)
/
(
1
-
P
(
z
)
)
2
,
-
4
z
). In fact (B.3), (2.79), and (B.10) are invariant by (
P
(
z
)
,
z
)
→
(
-
4
P
(
z
)
/
(
1
-
P
(
z
)
)
2
,
-
4
z
), but also (
P
(
z
)
,
z
)
→
(
-
(
1
-
P
(
z
)
)
2
/
4
/
P
(
z
)
,
-
z
/
4
), and also by (
P
(
z
)
,
z
)
→
(
1
/
P
(
z
)
,
z
).The functionP
(
z
) satisfies the functional equation:(2.93)
P
(
-
4
·
z
)
=
-
4
P
(
z
)
1
-
P
(
z
)
2
,
but also(2.94)
P
(
(
-
7
-
24
i
)
·
z
)
=
T
(
P
(
z
)
)
,
P
(
(
-
7
+
24
i
)
·
z
)
=
T
*
(
P
(
z
)
)
,
and, more generally, as can be checked order by order on series expansions(2.95)
P
(
a
1
·
z
)
=
R
a
1
(
P
(
z
)
)
.
For example, considering the “good” branch (2.35) for the inverse of -
4
z
/
(
1
-
z
)
2, namely S
-
1
/
4
(
1
)
(
z
), we can even check, order by order, on the series expansions of P
(
z
) and S
-
1
/
4
(
1
)
(
z
) the functional relation(2.96)
S
-
1
/
4
(
1
)
(
P
(
z
)
)
=
P
(
-
z
4
)
,
valid for
|
P
(
z
)
|
<
1
since the radius of convergence of
S
-
1
/
4
(
1
)
(
z
)
is
1.Recalling the functional equations (2.94) it is natural to say that if P
(
z
) is singular at z
=
z
s, then, for almost all the rational functions, in particular T
(
z
) (resp. T
*
(
z
)) the T
(
P
(
z
)
) is also singular z
=
z
s, and thus, from (2.94), P
(
z
) is also singular at z
=
(
-
7
±
24
i
)
·
z
s. It is thus extremely natural to see the emergence of the infinite number of singularities in (2.87) of the form z
=
(
N
1
+
i
·
N
2
)
·
z
s, as a consequence of (2.95) together with a reduction of the one-parameter series R
a
1
(
z
) to a rational function for an infinite number of selected values of a
1, namely the N
1
+
i
·
N
2 in (2.87). This is actually the case for all the values displayed in (2.87). For instance, for a
1
=
3
4
=
81 we get the following simple rational function:(2.97)
R
81
(
z
)
=
z
·
(
z
2
+
6
z
-
3
3
z
2
-
6
z
-
1
)
4
,
for which it is straightforward to verify that this rational transformation commutes with T
(
z
), T
*
(
z
), -
4
z
/
(
1
-
z
)
2, and is a solution of the Rota-Baxter-like functional equation (2.19). The case a
1
=
5
4
=
625 in (2.87), is even simpler, since it just requires to compose T
(
z
) and T
*
(
z
)
(2.98)
R
625
(
z
)
=
T
(
T
*
(
z
)
)
=
T
*
(
T
(
z
)
)
=
z
·
(
z
2
-
2
z
+
5
5
z
2
-
2
z
+
1
)
4
·
(
1
-
12
z
-
26
z
2
+
52
z
3
+
z
4
1
+
52
z
-
26
z
2
-
12
z
3
+
z
4
)
4
,
which, again verifies (2.19) and commutes with all the other rational functions, in particular (2.97). We also obtained the rational function corresponding to a
1
=
7
4
=
2401, namely:
(2.99)
R
2401
(
z
)
=
z
·
(
N
2401
(
z
)
D
2401
(
z
)
)
4
,
with:(2.100)
N
2401
(
z
)
=
z
12
·
D
2401
(
1
z
)
,
(2.101)
D
2401
(
z
)
=
1
+
196
z
-
1302
z
2
+
14756
z
3
-
15673
z
4
-
42168
z
5
+
111916
z
6
-
82264
z
7
+
35231
z
8
-
19852
z
9
+
2954
z
10
+
308
z
11
-
7
z
12
.The polynomialN
2401
(
z
) satisfies many functional equations, like, for instance (withR
-
4
(
z
)
=
-
4
z
/
(
1
-
z
2
)):(2.102)
4
12
·
D
2401
(
1
R
-
4
(
z
)
)
=
D
2401
(
z
)
·
D
2401
(
1
z
)
and also(2.103)
(
1
-
z
)
49
·
D
2401
(
R
-
4
(
z
)
)
2
=
D
2401
(
z
)
4
-
z
49
·
D
2401
(
1
z
)
4
.We also obtained the rational function corresponding toa
1
=
1
1
4
=
14641, namely,
(2.104)
R
14641
(
z
)
=
z
·
(
N
14641
(
z
)
D
14641
(
z
)
)
4
,
with:(2.105)
N
14641
(
z
)
=
z
30
·
D
14641
(
1
z
)
,
(2.106)
D
14641
(
z
)
=
1
+
1210
z
-
33033
z
2
+
2923492
z
3
+
5093605
z
4
-
385382514
z
5
+
3974726283
z
6
-
14323974808
z
7
+
57392757037
z
8
-
291359180310
z
9
+
948497199067
z
10
-
1642552094436
z
11
+
1084042069649
z
12
+
1890240552750
z
13
-
6610669151537
z
14
+
9712525647792
z
15
-
8608181312269
z
16
+
5384207244702
z
17
-
3223489742187
z
18
+
2175830922716
z
19
-
1197743580033
z
20
+
387221579866
z
21
-
50897017743
z
22
-
7864445336
z
23
+
5391243935
z
24
-
815789634
z
25
+
28366041
z
26
-
5092956
z
27
+
207691
z
28
+
2794
z
29
-
11
z
30
,
and, of course, one can verify that R
14641
(
z
) actually commutes with R
-
4, R
81, R
625, R
2401
(
z
), and is a solution of the Rota-Baxter-like functional equation (2.19). Similarly to R
2401
(
z
) (see (2.102), (2.103)), we also have the functional equations:(2.107)
4
30
·
D
14641
(
1
R
-
4
(
z
)
)
=
D
14641
(
z
)
·
D
14641
(
1
z
)
,
and also(2.108)
(
1
-
z
)
(
4
·
30
+
1
)
·
D
14641
(
R
-
4
(
z
)
)
2
=
D
14641
(
z
)
4
-
z
(
4
·
30
+
1
)
·
D
14641
(
1
z
)
4
.Next we obtained the rational function corresponding toa
1
=
1
3
4
=
28561, which verifies (2.19) namely,(2.109)
R
28561
(
z
)
=
z
·
(
N
28561
(
z
)
D
28561
(
z
)
)
4
,
with
N
28561
(
z
)
=
z
42
·
D
28561
(
1
z
)
,
N
28561
(
z
)
=
z
42
·
D
28561
(
1
z
)
,
(2.110)
D
28561
(
z
)
=
(
1
-
22
z
+
235
z
2
-
228
z
3
+
39
z
4
+
26
z
5
+
13
z
6
)
·
D
28561
(
36
)
(
z
)
D
28561
(
36
)
(
z
)
=
1
+
2388
z
-
61098
z
2
+
19225300
z
3
+
606593049
z
4
-
1543922656
z
5
+
7856476560
z
6
-
221753896032
z
7
+
1621753072244
z
8
-
4542779886736
z
9
+
2731418674664
z
10
+
36717669656304
z
11
-
200879613202428
z
12
+
547249607666784
z
13
-
934179604482832
z
14
+
1235038888776160
z
15
-
1788854212778642
z
16
+
3018407750933816
z
17
-
4349780716415868
z
18
+
4419228090228152
z
19
-
2899766501472914
z
20
+
931940880451552
z
21
+
413258559018224
z
22
-
857795672629664
z
23
+
659989056851972
z
24
-
304241349909008
z
25
+
87636987790824
z
26
-
14593362219920
z
27
+
1073204980340
z
28
+
45138167200
z
29
-
23660433008
z
30
+
2028597792
z
31
-
29540327
z
32
+
3238420
z
33
-
73386
z
34
-
492
z
35
+
z
36
.
We get similar results, mutatis mutandis, than the ones previously obtained (commutation, functional equations like (2.107), (2.108), etc.), namely,(2.111)
4
42
·
D
28561
(
1
R
-
4
(
z
)
)
=
D
28561
(
z
)
·
D
28561
(
1
z
)
,
…
.The “palindromic” nature of (2.97) (2.98), (2.99), (2.104) and (2.109) (see (2.100), (2.105)), (2.109)) corresponds to the fact that these rational transformations commute with J:(2.112)
1
R
81
(
z
)
=
R
81
(
1
z
)
,
1
R
625
(
z
)
=
R
625
(
1
z
)
,
…
.
In fact, more generally, we have R
N
4
(
1
/
z
)
=
1
/
R
N
4
(
z
) for N any odd integer (N
=
9,21
,
…) and R
N
4
(
1
/
z
)
=
R
N
4
(
z
) for any even integer  N.From (2.87) one can reasonnably conjecture that the fourth power of any integer will provide a new example of R
a
1
(
z
) being a rational function. The simple nontrivial example corresponds to the already found rational function(2.113)
R
16
(
z
)
=
16
·
z
·
(
1
-
z
)
2
(
z
+
1
)
4
.
We already have explicit rational functions for all values of a
1 of the form N
4 for N
=
2,3
,
…
,
16 and of course, we can in principle, build explicit rational functions for all the N's product of the previous integers. Along this line it is worth noticing that the coefficients of the series R
a
1
(
z
) are all integers when a
1 is the fourth power of any integer.We are thus starting to build an infinite number of (elementary)commuting rational transformations, any composition of these (infinite number of) rational transformations giving rational transformations satisfying (2.19) and preserving the linear differential operator Ω. This set of rational transformations is a pretty large set! Actually this set of rational transformations corresponds to the isogenies of the underlying elliptic function.The proliferation of the singularities ofP
(
z
) corresponds to this (pretty large) set of rational transformations. Recalling (2.96), the previous singularity argument is not valid
24 for the (well-suited) inverse transformations (S
-
1
/
4
(
1
)
(
z
), … ) of these rational transformations because (2.96) requires |
P
(
z
)
|
<
1 (corresponding to the radius of convergence of S
-
1
/
4
(
1
)
(
z
)) and the singularity z
=
z
s corresponds precisely to “hit” the value P
(
z
)
=
1.


## 2.6. Other Examples of Selected Gauss Hypergeometric ODEs

For heuristic reasons we have focused onA
(
z
)
=
(
3
-
5
z
)
/
z
/
(
1
-
z
)
/
4, but of course, one can find many other examples and try to generalize these examples.For instance, introducing(2.114)
A
(
z
)
=
1
6
·
d
ln
(
(
1
-
z
)
3
z
5
)
d
z
=
1
6
·
5
-
8
z
(
1
-
z
)
z
the rational transformation(2.115)
R
(
z
)
=
-
27
·
z
(
1
-
4
z
)
3
,
verifies the “Rota-Baxter-like” functional relation (2.19). This example corresponds to the following covariance [12] on a Gauss hypergeometric integral (of the c
=
1
+
b type, see below):(2.116)
2
F
1
(
[
1
2
,
1
6
]
,
[
7
6
]
;
z
)
=
z
-
1
/
6
6
·
∫
0
z
t
-
5
/
6
(
1
-
t
)
-
1
/
2
·
d
t
=
(
1
-
4
z
)
-
1
/
2
·
2
F
1
(
[
1
2
,
1
6
]
,
[
7
6
]
;
-
27
z
(
1
-
4
z
)
3
)
which is associated with the elliptic curve(2.117)
y
6
-
(
1
-
t
)
3
·
t
5
=
0
.
Another example (of the c
=
1
+
a type, see below) is(2.118)
A
(
z
)
=
1
3
·
d
ln
(
(
1
-
z
)
2
z
2
)
d
z
=
2
3
·
1
-
2
z
(
1
-
z
)
z
,
where the rational transformation(2.119)
R
(
z
)
=
z
·
(
z
-
2
)
3
(
1
-
2
z
)
3
=
-
8
z
-
36
z
2
-
126
z
3
-
387
z
4
+
⋯
verifies the “Rota-Baxter-like” functional relation (2.19). This example corresponds to the following covariance [12] on a Gauss hypergeometric integral:(2.120)
2
F
1
(
[
1
3
,
2
3
]
,
[
4
3
]
;
z
)
=
z
-
1
/
3
3
·
∫
0
z
t
-
2
/
3
(
1
-
t
)
-
2
/
3
·
d
t
=
1
2
·
2
-
z
1
-
2
z
·
2
F
1
(
[
1
3
,
2
3
]
,
[
4
3
]
;
z
·
(
z
-
2
)
3
(
1
-
2
z
)
3
)
,
which is associated with the elliptic curve(2.121)
y
3
-
(
1
-
t
)
2
·
t
2
=
0
.
Note that, similarly to the main example of the paper, there exist many rational transformations25 satisfying (2.19) that cannot be reduced to iterates of (2.119), for instance,(2.122)
T
(
z
)
=
-
27
·
z
·
(
1
-
z
)
(
z
2
-
z
+
1
)
3
(
z
3
+
3
z
2
-
6
z
+
1
)
3
=
-
27
z
-
378
z
2
-
3888
z
3
-
34074
z
4
-
271620
z
5
-
2032209
z
6
+
⋯
.
One verifies immediately that (2.122) actually verifies (2.19) with (2.118). Not surprisingly, the two rational transformations (2.119) and (2.122) commute.Another simple example with rational symmetries corresponds toΩ
=
(
D
z
+
A
(
z
)
)
·
D
z with(2.123)
A
(
z
)
=
-
1
2
·
3
z
-
1
z
(
1
-
z
)
=
1
2
·
d
ln
(
z
·
(
1
-
z
)
2
)
d
z
.
It has the simple (genus zero) hypergeometric solution26:(2.124)
ℱ
(
z
)
=
z
1
/
2
·
2
F
1
(
[
1
,
1
2
]
,
[
3
2
]
;
z
)
=
arctanh
(
z
1
/
2
)
.
The linear differential operator Ω is covariant under the change of variable z
→
1
/
z and z
→
R
(
z
), where27
(2.125)
R
(
z
)
=
4
z
(
1
+
z
)
2
.
One can easily check that (2.123) and (2.125) satisfy the functional equation (2.19). One also verifies that (2.123) and z
→
1
/
z or the iterates of (2.125) satisfy the functional equation (2.19). The solution of the adjoint operator are (
1
-
z
)
·
z
1
/
2 and(2.126)
F
(
z
)
=
z
·
(
1
-
z
)
·
2
F
1
(
[
1
,
1
2
]
,
[
3
2
]
;
z
)
=
z
1
/
2
·
(
1
-
z
)
·
arctanh
(
z
1
/
2
)
=
z
-
2
3
z
2
-
2
15
z
3
-
2
35
z
4
-
2
63
z
5
-
2
99
z
6
-
2
143
z
7
+
⋯
.
One verifies, again, that (2.126) and (2.125) commute, (2.126) corresponding to the “infinitesimal composition” of (2.125) (see (2.56)).A first natural generalization amounts to keeping the remarkable factorization (2.6) which will, in fact, reduce the covariance of a second-order operator to the covariance of a first-order operator.28 Such a situation occurs for Gauss hypergeometric functions2
F
1
(
[
a
,
b
]
,
[
1
+
a
]
;
z
) solution of the (
a
,
b
)-symmetric linear differential operator(2.127)
z
·
(
1
-
z
)
·
D
z
2
+
(
c
-
(
a
+
b
+
1
)
·
z
)
·
D
z
-
a
·
b
,
as soon as29
c
=
1
+
a. For instance(2.128)
ℱ
(
z
)
=
z
a
·
2
F
1
(
[
a
,
b
]
,
[
1
+
a
]
;
z
)
,
is an integral of a simple algebraic function and is solution with the constant function of the second-order operator(2.129)
Ω
=
(
D
z
+
(
a
-
b
-
1
)
z
+
1
-
a
z
·
(
1
-
z
)
)
·
D
z
=
(
D
z
+
d
ln
(
(
1
-
z
)
b
·
z
1
-
a
)
d
z
)
·
D
z
,
yielding a new A
(
z
):(2.130)
A
(
z
)
=
(
1
-
a
)
+
(
a
-
b
-
1
)
z
(
1
-
z
)
·
z
=
1
-
a
z
-
b
1
-
z
.
The adjoint of (2.129) has the simple solution z
1
-
a
·
(
1
-
z
)
b:(2.131)
F
(
z
)
=
z
·
(
1
-
z
)
b
·
2
F
1
(
[
a
,
b
]
,
[
1
+
a
]
;
z
)
.Due to the(
a
,
b
)-symmetry of (2.127) we have a similar result for c
=
1
+
b. The function ℱ
(
z
)
=
z
b
·
2
F
1
(
[
a
,
b
]
,
[
1
+
b
]
;
z
) is solution of (2.129) where a and b have been permuted:(2.132)
(
D
z
+
(
b
-
a
-
1
)
z
+
1
-
b
z
·
(
1
-
z
)
)
·
D
z
,
yielding another A
(
z
)
(2.133)
A
(
z
)
=
(
1
-
b
)
+
(
b
-
a
-
1
)
z
(
1
-
z
)
·
z
,
The adjoint of (2.132) has the solution (
1
-
z
)
a
·
z
1
-
b together with the hypergeometric function:(2.134)
F
(
z
)
=
z
·
(
1
-
z
)
a
·
2
F
1
(
[
a
,
b
]
,
[
1
+
b
]
;
z
)
,
where one recovers the previous result (2.126).We are seeking for (Gauss hypergeometric) second-order differential equations30 with an infinite number of (hopefully rational, if not algebraic) symmetries: this is another way to say that we are not looking for generic Gauss hypergeometric differential equations, but Gauss hypergeometric differential equations related to elliptic curves, and thus having an infinite set of such isogenies. We are necessarily in the framework where the two parameters a and b of the Gauss hypergeometric are rational numbers in order to have integral of algebraic functions (yielding globally nilpotent [11] second-order differential operators). Let us denote by D the common denominator of the two rational numbers a
=
N
a
/
D and b
=
N
b
/
D, function (2.128) is associated to a period of the algebraic curve(2.135)
y
D
=
(
1
-
t
)
N
b
·
t
D
-
N
a
.
We just need to restrict to triplets of integers (
N
a
,
N
b
,
D
) such that the previous curve is an elliptic curve.Let us give an example (of thec
=
1
+
b type) that does not correspond to a genus one curve, with(2.136)
2
F
1
(
[
1
3
,
1
6
]
,
[
7
6
]
;
z
)
=
1
6
·
z
-
1
/
6
·
∫
0
z
t
-
5
/
6
(
1
-
t
)
-
1
/
3
·
d
t
,
which corresponds to the genus two curve:(2.137)
y
6
-
(
1
-
t
)
2
·
t
5
=
0
.
Again one introduces A
(
z
)
(2.138)
A
(
z
)
=
1
6
·
d
ln
(
(
1
-
z
)
2
z
5
)
d
z
=
1
6
·
5
-
7
z
z
·
(
1
-
z
)
and seeks for R
(
z
) as series expansions analytical at z
=
0. One gets actually, order by order, a one-parameter family(2.139)
R
a
1
(
z
)
=
a
1
·
z
-
2
7
a
1
·
(
a
1
-
1
)
·
z
2
+
1
637
a
1
·
(
a
1
-
1
)
·
(
17
a
1
-
87
)
·
z
3
-
2
84721
a
1
·
(
a
1
-
1
)
·
(
113
a
1
2
-
856
a
1
+
3438
)
·
z
4
-
1
38548055
a
1
·
(
a
1
-
1
)
·
(
3674
a
1
3
+
121194
a
1
2
-
552261
a
1
+
2095059
)
·
z
5
+
⋯
+
1
+
ϵ
n
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
·
z
n
+
⋯
,
where ϵ
n
=
0 for n odd and ϵ
n
=
1 for n even, and N
(
n
) is a (large) integer depending on n, and P
n
(
a
1
) is a polynomial with integer coefficients of degree n
-
2. One easily verifies, order by order, that one gets a one-parameter family of transformations commuting for different values of the parameter:(2.140)
R
a
1
(
R
b
1
(
z
)
)
=
R
b
1
(
R
a
1
(
z
)
)
=
R
a
1
b
1
(
z
)
.As far as the “algorithmic complexity” of this series (2.139) is concerned it is worth noticing that the degree growth [18] of the series coefficients is actually linear and not exponential as we could expect [19] at first sight. Even if this series is transcendental, it is not a “wild” series.Seeking for selected values ofa
1 such that the previous series (2.139) reduces to a rational function one can try to reproduce the simple calculations (2.48), (2.49), but unfortunately “shooting in the dark” because we have no hint of a well-suited denominator (if any!) like the polynomials in the lhs of (2.48), (2.49).It is also worth noticing that if we slightly changeA
(
z
) into(2.141)
A
(
z
)
=
1
N
·
d
ln
(
(
1
-
z
)
2
z
5
)
d
z
=
1
N
·
5
-
7
z
z
·
(
1
-
z
)
,
the algebraic curve (2.137) becomes y
N
-
(
1
-
t
)
2
·
t
5
=
0 which has, for instance genus five for N
=
11, but genus zero for N
=
7. For any of these cases of (2.141) one can easily get, order by order, a one-parameter series R
a
1 totally similar to (2.139) with, again, polynomials P
n
(
a
1
) of degree n
-
2.The first coefficienta
2 is in general(2.142)
a
2
=
-
2
2
N
-
5
·
a
1
·
(
a
1
-
1
)
.For the genus zero case,N
=
7
(2.143)
a
2
=
-
2
9
·
a
1
·
(
a
1
-
1
)
,
a
3
=
-
1
1296
·
a
1
·
(
a
1
-
1
)
·
(
127
-
a
1
)
,
a
4
=
-
1
134136
·
a
1
·
(
a
1
-
1
)
·
(
254
a
1
2
+
185
a
1
+
7499
)
,
…
,
a
n
=
-
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
,
which corresponds to the solution(2.144)
2
7
·
∫
0
z
z
-
5
/
7
·
(
1
-
z
)
-
2
/
7
·
d
t
=
z
2
/
7
·
2
F
1
(
[
2
7
,
2
7
]
,
[
9
7
]
;
z
)
.
Using the parametrization of the genus zero curve
(2.145)
y
=
-
(
u
+
1
)
2
·
u
5
(
u
+
1
)
7
-
u
7
,
t
=
-
u
7
(
u
+
1
)
7
-
u
7
,
one can actually perform the integration (2.144) of d
t
/
y and get an alternative form of the hypergeometric function (2.144):(2.146)
∫
0
z
z
-
5
/
7
·
(
1
-
z
)
-
2
/
7
·
d
t
=
∫
0
u
ρ
(
u
)
·
d
u
=
∫
0
v
v
1
-
v
7
·
d
v
,
where
:
z
=
-
u
7
(
u
+
1
)
7
-
u
7
,
ρ
(
u
)
=
(
u
+
1
)
4
·
u
(
u
+
1
)
7
-
u
7
,
and
:
v
=
u
1
+
u
,
z
=
v
7
v
7
-
1
.
Except transformations like v
→
ω
·
v (with ω
7
=
1) which have no impact on z, it seems difficult to find rational symmetries in this genus zero case.ForN
=
11 (genus five) the first successive coefficients read:(2.147)
a
2
=
-
2
17
·
a
1
·
(
a
1
-
1
)
,
a
3
=
-
1
8092
·
a
1
·
(
a
1
-
1
)
·
(
143
a
1
+
367
)
,
a
4
=
-
1
206346
·
a
1
·
(
a
1
-
1
)
·
(
1186
a
1
2
+
2473
a
1
+
5011
)
,
…
,
a
n
=
-
1
N
(
n
)
·
a
1
·
(
a
1
-
1
)
·
P
n
(
a
1
)
.
The “infinitesimal composition” function F
(
z
) (see (2.56), (2.57), and (2.58)) reads,(2.148)
F
(
z
)
=
∂
R
a
1
∂
a
1
|
a
1
=
1
=
z
-
2
17
z
2
-
15
238
z
3
-
5
119
z
4
-
37
1190
z
5
-
888
36295
z
6
-
2183
108885
z
7
-
4366
258213
z
8
-
58941
4045337
z
9
-
1807524
141586795
z
10
-
46543743
4106017055
z
11
-
5305986702
521464165985
z
12
+
⋯
and again we can actually check that this is actually the series expansion of the hypergeometric function(2.149)
z
·
(
1
-
z
)
·
2
F
1
(
[
1
,
15
11
]
,
[
17
11
]
;
z
)
,
solution of Ω
* the adjoint of the Ω linear differential operator corresponding to this (genus 5) N
=
11 case:(2.150)
Ω
*
=
D
z
·
(
D
z
+
1
11
·
7
z
-
5
z
·
(
1
-
z
)
)
.We have similar results for (2.128), (2.129), (2.130). As far as these one-parameter families of transformations R
a
1 are concerned, the only difference between the generic cases corresponding to arbitrary genus and genus one cases like (2.118) is that in the generic higher genus case, only a finite number of values of the parameter a
1 can correspond to rational functions. Note that this higher genus result generalizes to the arbitrary genus Gauss hypergeometric functions (2.128) and associated operators (2.129) and function (2.130). In this general case one can also get order by order a one-parameter family of transformations R
a
1 satisfying a commutation relation (2.141).Note thatR
(
z
)
=
1
/
z is actually a solution of (2.19) for this genus-two example (2.139). Along this line of selected R
(
z
) solutions of (2.19) many interesting subcases of this general case (2.128), (2.129), (2.130) are given in Appendix C.1.In our previous genus-one examples, with this close identification between therenormalization group and the isogenies of elliptic curves, we saw that, in order to obtain linear differential operators covariant by an infinite number of transformations (rational or algebraic), we must restrict our second-order Gauss hypergeometric differential operator to Gauss hypergeometric associated to elliptic curves (see Appendices C and D). Beyond this framework we still have one-parameter families (see (2.141)) but we cannot expect an infinite number of rational (and probably algebraic) transformations to be particular cases of such families of transcendental transformations.


## 3. Conclusion

We have shown that several selected Gauss hypergeometric linear differential operators associated to elliptic curves, and  factorised into order-one linear differential operators, actually present aninfinite number of rational symmetries that actually identify with the isogenies of the associated elliptic curves that are perfect illustrations of exact representations of the renormalization group. We actually displayed all these calculations, results, and structures because they are perfect examples of exact renormalization transformations. For more realistic models (corresponding to Yang-Baxter models with elliptic parametrizations), the previous calculations and structures become more involved and subtle, the previous rational transformations being replaced by algebraic transformations corresponding to modular curves. For instance, in our models of lattice statistical mechanics (or enumerative combinatorics, etc.), we are often getting globally nilpotent linear differential operators [11] of quite high orders [20–24] that, in fact, factor into globally nilpotent operators of smaller orders31 which, for Yang-Baxter integrable models with a canonical elliptic parametrization, must necessarily “be associated with elliptic curves.” Appendix D provides some calculations showing that the integral for χ
(
2
), the two-particle contribution of the susceptibility of the Ising model [25–27], is clearly and straightforwardly associated with an elliptic curve.We wanted to highlight the importance ofexplicit constructions in answering difficult or subtle questions.All the calculations displayed in this paper are elementary calculations given explicitly for heuristic reasons. The simple calculations (in particular with the introduction of a simple Rota-Baxter like functional equation) should be seen as some undergraduate training to more realistic renormalization calculations that will require a serious knowledge of fundamental modular curves,modular forms, Hauptmoduls, Gauss-Manin or Picard-Fuchs structures [28, 29] and, beyond, some knowledge of mirror symmetries [30–34] of Calabi-Yau manifolds, these mirror symmetries generalizing32  the Hauptmodul structure for elliptic curves.



---

*Source: 941560-2010-04-29.xml*