
We previously reported on a recursive method to generate the expansion of the
lattice Green function of the $d$dimensional facecentred cubic lattice (fcc).
The method was used to generate many coefficients for d=7 and the corresponding
linear differential equation has been obtained. In this paper, we show the
strength and the limit of the method by producing the series and the
corresponding linear differential equations for d=8, 9, 10, 11, 12. The
differential Galois groups of these linear differential equations are shown to
be symplectic for d=8, 10, 12 and orthogonal for d= 9, 11. The recursion
relation naturally provides a 2dimensional array $ T_d(n,j)$ where only the
coefficients $ t_d(n,0)$ correspond to the coefficients of the lattice Green
function of the ddimensional fcc. The coefficients $ t_d(n,j)$ are associated
to Dfinite bivariate series annihilated by linear partial differential
equations that we analyze.

We show and give the linear differential operators ${\cal L}^{scal}_q$ of
order q= n^2/4+n+7/8+(1)^n/8, for the integrals $I_n(r)$ which appear in the
twopoint correlation scaling function of Ising model $ F_{\pm}(r)=
\lim_{scaling} {\cal M}_{\pm}^{2}
< \sigma_{0,0} \, \sigma_{M,N}> = \sum_{n} I_{n}(r)$. The integrals $
I_{n}(r)$ are given in expansion around r= 0 in the basis of the formal
solutions of $\, {\cal L}^{scal}_q$ with transcendental combination
coefficients. We find that the expression $ r^{1/4}\,\exp(r^2/8)$ is a solution
of the Painlev\'e VI equation in the scaling limit. Combinations of the
(analytic at $ r= 0$) solutions of $ {\cal L}^{scal}_q$ sum to $ \exp(r^2/8)$.
We show that the expression $ r^{1/4} \exp(r^2/8)$ is the scaling limit of the
correlation function $ C(N, N)$ and $ C(N, N+1)$. The differential Galois
groups of the factors occurring in the operators $ {\cal L}^{scal}_q$ are
given.

We revisit an ordersix linear differential operator having a solution which
is a diagonal of a rational function of three variables. Its exterior square
has a rational solution, indicating that it has a selected differential Galois
group, and is actually homomorphic to its adjoint. We obtain the two
corresponding intertwiners giving this homomorphism to the adjoint. We show
that these intertwiners are also homomorphic to their adjoint and have a simple
decomposition, already underlined in a previous paper, in terms of ordertwo
selfadjoint operators. From these results, we deduce a new form of
decomposition of operators for this selected ordersix linear differential
operator in terms of three ordertwo selfadjoint operators. We then generalize
the previous decomposition to decompositions in terms of an arbitrary number of
selfadjoint operators of the same parity order. This yields an infinite family
of linear differential operators homomorphic to their adjoint, and, thus, with
a selected differential Galois group. We show that the equivalence of such
operators is compatible with these canonical decompositions. The rational
solutions of the symmetric, or exterior, squares of these selected operators
are, noticeably, seen to depend only on the rightmost selfadjoint operator in
the decomposition. These results, and tools, are applied on operators of large
orders. For instance, it is seen that a large set of (quite massive) operators,
associated with reflexive 4polytopes defining CalabiYau 3folds, obtained
recently by P. Lairez, correspond to a particular form of the decomposition
detailed in this paper.

We show that the globally nilpotent Goperators corresponding to the factors
of the linear differential operators annihilating the multifold integrals
$\chi^{(n)}$ of the magnetic susceptibility of the Ising model ($n \le 6$) are
homomorphic to their adjoint. This property of being selfadjoint up to
operator homomorphisms, is equivalent to the fact that their symmetric square,
or their exterior square, have rational solutions. The differential Galois
groups are in the special orthogonal, or symplectic, groups. This selfadjoint
(up to operator equivalence) property means that the factor operators we
already know to be Derived from Geometry, are special globally nilpotent
operators: they correspond to "Special Geometries".
Beyond the small order factor operators (occurring in the linear differential
operators associated with $ \chi^{(5)}$ and $ \chi^{(6)}$), and, in particular,
those associated with modular forms, we focus on the quite large ordertwelve
and order23 operators. We show that the ordertwelve operator has an exterior
square which annihilates a rational solution. Then, its differential Galois
group is in the symplectic group $ Sp(12, \mathbb{C})$. The order23 operator
is shown to factorize in an ordertwo operator and an order21 operator. The
symmetric square of this order21 operator has a rational solution. Its
differential Galois group is, thus, in the orthogonal group
$ SO(21, \mathbb{C})$.

We show that the nfold integrals $\chi^{(n)}$ of the magnetic susceptibility
of the Ising model, as well as various other nfold integrals of the "Ising
class", or nfold integrals from enumerative combinatorics, like lattice Green
functions, correspond to a distinguished class of function generalising
algebraic functions: they are actually diagonals of rational functions. As a
consequence, the power series expansions of the, analytic at x=0, solutions of
these linear differential equations "Derived From Geometry" are globally
bounded, which means that, after just one rescaling of the expansion variable,
they can be cast into series expansions with integer coefficients. We also give
several results showing that the unique analytical solution of CalabiYau ODEs,
and, more generally, PicardFuchs linear ODEs, with solutions of maximal
weights, are always diagonal of rational functions. Besides, in a more
enumerative combinatorics context, generating functions whose coefficients are
expressed in terms of nested sums of products of binomial terms can also be
shown to be diagonals of rational functions. We finally address the question of
the relations between the notion of integrality (series with integer
coefficients, or, more generally, globally bounded series) and the modularity
of ODEs.

We show that the nfold integrals $\chi^{(n)}$ of the magnetic susceptibility
of the Ising model, as well as various other nfold integrals of the "Ising
class", or nfold integrals from enumerative combinatorics, like lattice Green
functions, are actually diagonals of rational functions. As a consequence, the
power series expansions of these solutions of linear differential equations
"Derived From Geometry" are globally bounded, which means that, after just one
rescaling of the expansion variable, they can be cast into series expansions
with integer coefficients. Besides, in a more enumerative combinatorics
context, we show that generating functions whose coefficients are expressed in
terms of nested sums of products of binomial terms can also be shown to be
diagonals of rational functions. We give a large set of results illustrating
the fact that the unique analytical solution of CalabiYau ODEs, and more
generally of MUM ODEs, is, almost always, diagonal of rational functions. We
revisit Christol's conjecture that globally bounded series of Goperators are
necessarily diagonals of rational functions. We provide a large set of examples
of globally bounded series, or series with integer coefficients, associated
with modular forms, or Hadamard product of modular forms, or associated with
CalabiYau ODEs, underlying the concept of modularity. We finally address the
question of the relations between the notion of integrality (series with
integer coefficients, or, more generally, globally bounded series) and the
modularity (in particular integrality of the Taylor coefficients of mirror
map), introducing new representations of Yukawa couplings.

Lattice statistical mechanics, often provides a natural (holonomic) framework
to perform singularity analysis with several complex variables that would, in a
general mathematical framework, be too complex, or could not be defined.
Considering several PicardFuchs systems of twovariables "above" CalabiYau
ODEs, associated with double hypergeometric series, we show that holonomic
functions are actually a good framework for actually finding the singular
manifolds. We, then, analyse the singular algebraic varieties of the nfold
integrals $ \chi^{(n)}$, corresponding to the decomposition of the magnetic
susceptibility of the anisotropic square Ising model. We revisit a set of
Nickelian singularities that turns out to be a twoparameter family of elliptic
curves. We then find a first set of nonNickelian singularities for $
\chi^{(3)}$ and $ \chi^{(4)}$, that also turns out to be rational or ellipic
curves. We underline the fact that these singular curves depend on the
anisotropy of the Ising model. We address, from a birational viewpoint, the
emergence of families of elliptic curves, and of CalabiYau manifolds on such
problems. We discuss the accumulation of these singular curves for the
nonholonomic anisotropic full susceptibility.

We give the exact expressions of the partial susceptibilities $\chi^{(3)}_d$
and $\chi^{(4)}_d$ for the diagonal susceptibility of the Ising model in terms
of modular forms and CalabiYau ODEs, and more specifically,
$_3F_2([1/3,2/3,3/2],\, [1,1];\, z)$ and $_4F_3([1/2,1/2,1/2,1/2],\, [1,1,1];
\, z)$ hypergeometric functions. By solving the connection problems we
analytically compute the behavior at all finite singular points for
$\chi^{(3)}_d$ and $\chi^{(4)}_d$. We also give new results for $\chi^{(5)}_d$.
We see in particular, the emergence of a remarkable ordersix operator, which
is such that its symmetric square has a rational solution. These new exact
results indicate that the linear differential operators occurring in the
$n$fold integrals of the Ising model are not only "Derived from Geometry"
(globally nilpotent), but actually correspond to "Special Geometry"
(homomorphic to their formal adjoint). This raises the question of seeing if
these "special geometry" Isingoperators, are "special" ones, reducing, in fact
systematically, to (selected, kbalanced, ...) $_{q+1}F_q$ hypergeometric
functions, or correspond to the more general solutions of CalabiYau equations.

We show that almost all the linear differential operators factors obtained in
the analysis of the nparticle contribution of the susceptibility of the Ising
model for $\, n \le 6$, are operators "associated with elliptic curves". Beyond
the simplest factors which are homomorphic to symmetric powers of the second
order operator associated with the complete elliptic integral E, the second and
third order differential operators can actually be interpreted as modular forms
of the elliptic curve of the Ising model. A last orderfour globally nilpotent
operator is not reducible to this elliptic curve, modular forms scheme. It is
shown to actually correspond to a natural generalization of this elliptic
curve, modular forms scheme, with the emergence of a CalabiYau equation,
corresponding to a selected $_4F_3$ hypergeometric function which can also be
seen as a Hadamard product of the complete elliptic integral K, with a
remarkably simple algebraic pullback, the corresponding CalabiYau
fourthorder differential operator having a symplectic differential Galois
group SP(4,C). The associated mirror maps and higher order Schwarzian ODEs has
an exact (isogenies) representation of the generators of the renormalization
group, extending the modular group SL(2,Z) to a GL(2, Z) symmetry group.

We obtain in exact arithmetic the order 24 linear differential operator
$L_{24}$ and right hand side $E^{(5)}$ of the inhomogeneous
equation$L_{24}(\Phi^{(5)}) = E^{(5)}$, where $\Phi^{(5)}
=\tilde{\chi}^{(5)}\tilde{\chi}^{(3)}/2+\tilde{\chi}^{(1)}/120$ is a linear
combination of $n$particle contributions to the susceptibility of the square
lattice Ising model. In Bostan, et al. (J. Phys. A: Math. Theor. {\bf 42},
275209 (2009)) the operator $L_{24}$ (modulo a prime) was shown to factorize
into $L_{12}^{(\rm left)} \cdot L_{12}^{(\rm right)}$; here we prove that no
further factorization of the order 12 operator $L_{12}^{(\rm left)}$ is
possible. We use the exact ODE to obtain the behaviour of $\tilde{\chi}^{(5)}$
at the ferromagnetic critical point and to obtain a limited number of analytic
continuations of $\tilde{\chi}^{(5)}$ beyond the principal disk defined by its
high temperature series. Contrary to a speculation in Boukraa, et al (J. Phys.
A: Math. Theor. {\bf 41} 455202 (2008)), we find that $\tilde{\chi}^{(5)}$ is
singular at $w=1/2$ on an infinite number of branches.

We review developments made since 1959 in the search for a closed form for
the susceptibility of the Ising model. The expressions for the form factors in
terms of the nome $q$ and the modulus $k$ are compared and contrasted. The
$\lambda$ generalized correlations $C(M,N;\lambda)$ are defined and explicitly
computed in terms of theta functions for $M=N=0,1$.

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.

This paper deals with $\tilde{\chi}^{(6)}$, the sixparticle contribution to
the magnetic susceptibility of the square lattice Ising model. We have
generated, modulo a prime, series coefficients for $\tilde{\chi}^{(6)}$. The
length of the series is sufficient to produce the corresponding Fuchsian linear
differential equation (modulo a prime). We obtain the Fuchsian linear
differential equation that annihilates the "depleted" series
$\Phi^{(6)}=\tilde{\chi}^{(6)}  {2 \over 3} \tilde{\chi}^{(4)} + {2 \over 45}
\tilde{\chi}^{(2)}$. The factorization of the corresponding differential
operator is performed using a method of factorization modulo a prime introduced
in a previous paper. The "depleted" differential operator is shown to have a
structure similar to the corresponding operator for $\tilde{\chi}^{(5)}$. It
splits into factors of smaller orders, with the leftmost factor of order six
being equivalent to the symmetric fifth power of the linear differential
operator corresponding to the elliptic integral $E$. The rightmost factor has
a direct sum structure, and using series calculated modulo several primes, all
the factors in the direct sum have been reconstructed in exact arithmetics.

A detailed study is presented of the expected performance of the ATLAS
detector. The reconstruction of tracks, leptons, photons, missing energy and
jets is investigated, together with the performance of btagging and the
trigger. The physics potential for a variety of interesting physics processes,
within the Standard Model and beyond, is examined. The study comprises a series
of notes based on simulations of the detector and physics processes, with
particular emphasis given to the data expected from the first years of
operation of the LHC at CERN.

We consider some twodimensional birational transformations. One of them is a
birational deformation of the H\'enon map. For some of these birational
mappings, the post critical set (i.e. the iterates of the critical set) is
infinite and we show that this gives straightforwardly the algebraic covariant
curves of the transformation when they exist. These covariant curves are used
to build the preserved meromorphic twoform. One may have also an infinite post
critical set yielding a covariant curve which is not algebraic (transcendent).
For two of the birational mappings considered, the post critical set is not
infinite and we claim that there is no algebraic covariant curve and no
preserved meromorphic twoform. For these two mappings with non infinite post
critical sets, attracting sets occur and we show that they pass the usual tests
(Lyapunov exponents and the fractal dimension) for being strange attractors.
The strange attractor of one of these two mappings is unbounded.

We consider the Fuchsian linear differential equation obtained (modulo a
prime) for $\tilde{\chi}^{(5)}$, the fiveparticle contribution to the
susceptibility of the square lattice Ising model. We show that one can
understand the factorization of the corresponding linear differential operator
from calculations using just a single prime. A particular linear combination of
$\tilde{\chi}^{(1)}$ and $\tilde{\chi}^{(3)}$ can be removed from
$\tilde{\chi}^{(5)}$ and the resulting series is annihilated by a high order
globally nilpotent linear ODE. The corresponding (minimal order) linear
differential operator, of order 29, splits into factors of small orders. A
fifth order linear differential operator occurs as the leftmost factor of the
"depleted" differential operator and it is shown to be equivalent to the
symmetric fourth power of $L_E$, the linear differential operator corresponding
to the elliptic integral $E$. This result generalizes what we have found for
the lower order terms $\tilde{\chi}^{(3)}$ and $\tilde{\chi}^{(4)}$. We
conjecture that a linear differential operator equivalent to a symmetric
$(n1)$th power of $L_E$ occurs as a leftmost factor in the minimal order
linear differential operators for all $\tilde{\chi}^{(n)}$'s.

We recall various multiple integrals related to the isotropic square Ising
model, and corresponding, respectively, to the nparticle contributions of the
magnetic susceptibility, to the (lattice) form factors, to the twopoint
correlation functions and to their lambdaextensions. These integrals are
holonomic and even Gfunctions: they satisfy Fuchsian linear differential
equations with polynomial coefficients and have some arithmetic properties. We
recall the explicit forms, found in previous work, of these Fuchsian equations.
These differential operators are very selected Fuchsian linear differential
operators, and their remarkable properties have a deep geometrical origin: they
are all globally nilpotent, or, sometimes, even have zero pcurvature. Focusing
on the factorised parts of all these operators, we find out that the global
nilpotence of the factors corresponds to a set of selected structures of
algebraic geometry: elliptic curves, modular curves, and even a remarkable
weight1 modular form emerging in the threeparticle contribution $ \chi^{(3)}$
of the magnetic susceptibility of the square Ising model. In the case where we
do not have Gfunctions, but Hamburger functions (one irregular singularity at
0 or $ \infty$) that correspond to the confluence of singularities in the
scaling limit, the pcurvature is also found to verify new structures
associated with simple deformations of the nilpotent property.

We calculate very long low and hightemperature series for the
susceptibility $\chi$ of the square lattice Ising model as well as very long
series for the fiveparticle contribution $\chi^{(5)}$ and sixparticle
contribution $\chi^{(6)}$. These calculations have been made possible by the
use of highly optimized polynomial time modular algorithms and a total of more
than 150000 CPU hours on computer clusters. For $\chi^{(5)}$ 10000 terms of the
series are calculated {\it modulo} a single prime, and have been used to find
the linear ODE satisfied by $\chi^{(5)}$ {\it modulo} a prime.
A diffPad\'e analysis of 2000 terms series for $\chi^{(5)}$ and $\chi^{(6)}$
confirms to a very high degree of confidence previous conjectures about the
location and strength of the singularities of the $n$particle components of
the susceptibility, up to a small set of ``additional'' singularities. We find
the presence of singularities at $w=1/2$ for the linear ODE of $\chi^{(5)}$,
and $w^2= 1/8$ for the ODE of $\chi^{(6)}$, which are {\it not} singularities
of the ``physical'' $\chi^{(5)}$ and $\chi^{(6)},$ that is to say the
seriessolutions of the ODE's which are analytic at $w =0$.
Furthermore, analysis of the long series for $\chi^{(5)}$ (and $\chi^{(6)}$)
combined with the corresponding long series for the full susceptibility $\chi$
yields previously conjectured singularities in some $\chi^{(n)}$, $n \ge 7$.
We also present a mechanism of resummation of the logarithmic singularities
of the $\chi^{(n)}$ leading to the known powerlaw critical behaviour occurring
in the full $\chi$, and perform a power spectrum analysis giving strong
arguments in favor of the existence of a natural boundary for the full
susceptibility $\chi$.

We introduce some multiple integrals that are expected to have the same
singularities as the singularities of the $ n$particle contributions
$\chi^{(n)}$ to the susceptibility of the square lattice Ising model. We find
the Fuchsian linear differential equation satisfied by these multiple integrals
for $ n=1, 2, 3, 4$ and only modulo some primes for $ n=5$ and $ 6$, thus
providing a large set of (possible) new singularities of the $\chi^{(n)}$. We
discuss the singularity structure for these multiple integrals by solving the
Landau conditions. We find that the singularities of the associated ODEs
identify (up to $n= 6$) with the leading pinch Landau singularities. The second
remarkable obtained feature is that the singularities of the ODEs associated
with the multiple integrals reduce to the singularities of the ODEs associated
with a {\em finite number of one dimensional integrals}. Among the
singularities found, we underline the fact that the quadratic polynomial
condition $ 1+3 w +4 w^2 = 0$, that occurs in the linear differential equation
of $ \chi^{(3)}$, actually corresponds to a remarkable property of selected
elliptic curves, namely the occurrence of complex multiplication. The
interpretation of complex multiplication for elliptic curves as complex fixed
points of the selected generators of the renormalization group, namely
isogenies of elliptic curves, is sketched. Most of the other singularities
occurring in our multiple integrals are not related to complex multiplication
situations, suggesting an interpretation in terms of (motivic) mathematical
structures beyond the theory of elliptic curves.

We use the recently derived form factor expansions of the diagonal twopoint
correlation function of the square Ising model to study the susceptibility for
a magnetic field applied only to one diagonal of the lattice, for the isotropic
Ising model.
We exactly evaluate the one and two particle contributions $\chi_{d}^{(1)}$
and $\chi_{d}^{(2)}$ of the corresponding susceptibility, and obtain linear
differential equations for the three and four particle contributions, as well
as the five particle contribution ${\chi}^{(5)}_d(t)$, but only modulo a given
prime. We use these exact linear differential equations to show that, not only
the russiandoll structure, but also the direct sum structure on the linear
differential operators for the $ n$particle contributions $\chi_{d}^{(n)}$ are
quite directly inherited from the direct sum structure on the form factors $
f^{(n)}$.
We show that the $ n^{th}$ particle contributions $\chi_{d}^{(n)}$ have their
singularities at roots of unity. These singularities become dense on the unit
circle $\sinh2E_v/kT \sinh 2E_h/kT=1$ as $ n\to \infty$.

We briefly recall the FuchsPainlev\'e elliptic representation of Painlev\'e
VI. We then show that the polynomiality of the expressions of the correlation
functions (and form factors) in terms of the complete elliptic integral of the
first and second kind,
$ K$ and $ E$, is a straight consequence of the fact that the differential
operators corresponding to the entries of Toeplitzlike determinants, are
equivalent to the second order operator $ L_E$ which has $ E$ as solution (or,
for offdiagonal correlations to the direct sum of $ L_E$ and $ d/dt$). We show
that this can be generalized, mutatis mutandis, to the anisotropic Ising model.
The singledout second order linear differential operator $ L_E$ being replaced
by an isomonodromic system of two thirdorder linear partial differential
operators associated with $ \Pi_1$, the Jacobi's form of the complete elliptic
integral of the third kind (or equivalently two second order linear partial
differential operators associated with Appell functions, where one of these
operators can be seen as a deformation of $ L_E$). We finally explore the
generalizations, to the anisotropic Ising models, of the links we made, in two
previous papers, between Painlev\'e nonlinear ODE's, Fuchsian linear ODE's and
elliptic curves. In particular the elliptic representation of Painlev\'e VI has
to be generalized to an ``Appellian'' representation of Garnier systems.

We consider families of multiple and simple integrals of the ``Ising class''
and the linear ordinary differential equations with polynomial coefficients
they are solutions of. We compare the full set of singularities given by the
roots of the head polynomial of these linear ODE's and the subset of
singularities occurring in the integrals, with the singularities obtained from
the Landau conditions. For these Ising class integrals, we show that the Landau
conditions can be worked out, either to give the singularities of the
corresponding linear differential equation or the singularities occurring in
the integral. The singular behavior of these integrals is obtained in the
selfdual variable $w= s/2/(1+s^2)$, with $s= \sinh(2K)$, where $K=J/kT$ is the
usual Ising model coupling constant. Switching to the variable $s$, we show
that the singularities of the analytic continuation of series expansions of
these integrals actually break the KramersWannier duality. We revisit the
singular behavior (J. Phys. A {\bf 38} (2005) 94399474) of the third
contribution to the magnetic susceptibility of Ising model $\chi^{(3)}$ at the
points $1+3w+4w^2= 0$ and show that $\chi^{(3)}(s)$ is not singular at the
corresponding points inside the unit circle $ s =1$, while its analytical
continuation in the variable $s$ is actually singular at the corresponding
points $ 2+s+s^2=0$ oustside the unit circle ($ s  > 1$).

We study the Ising model twopoint diagonal correlation function $ C(N,N)$ by
presenting an exponential and form factor expansion in an integral
representation which differs from the known expansion of Wu, McCoy, Tracy and
Barouch. We extend this expansion, weighting, by powers of a variable
$\lambda$, the $j$particle contributions, $ f^{(j)}_{N,N}$. The corresponding
$ \lambda$ extension of the twopoint diagonal correlation function, $ C(N,N;
\lambda)$, is shown, for arbitrary $\lambda$, to be a solution of the sigma
form of the Painlev{\'e} VI equation introduced by Jimbo and Miwa. Linear
differential equations for the form factors $ f^{(j)}_{N,N}$ are obtained and
shown to have both a ``Russian doll'' nesting, and a decomposition of the
differential operators as a direct sum of operators equivalent to symmetric
powers of the differential operator of the elliptic integral $ E$. Each $
f^{(j)}_{N,N}$ is expressed polynomially in terms of the elliptic integrals $
E$ and $ K$. The scaling limit of these differential operators breaks the
direct sum structure but not the ``Russian doll'' structure. The previous $
\lambda$extensions, $ C(N,N; \lambda)$ are, for singledout values $ \lambda=
\cos(\pi m/n)$ ($m, n$ integers), also solutions of linear differential
equations. These solutions of Painlev\'e VI are actually algebraic functions,
being associated with modular curves.

The sigma form of the Painlev{\'e} VI equation contains four arbitrary
parameters and generically the solutions can be said to be genuinely
``nonlinear'' because they do not satisfy linear differential equations of
finite order. However, when there are certain restrictions on the four
parameters there exist one parameter families of solutions which do satisfy
(Fuchsian) differential equations of finite order. We here study this phenomena
of Fuchsian solutions to the Painlev{\'e} equation with a focus on the
particular PVI equation which is satisfied by the diagonal correlation function
C(N,N) of the Ising model. We obtain Fuchsian equations of order $N+1$ for
C(N,N) and show that the equation for C(N,N) is equivalent to the $N^{th}$
symmetric power of the equation for the elliptic integral $E$.
We show that these Fuchsian equations correspond to rational algebraic curves
with an additional Riccati structure and we show that the Malmquist Hamiltonian
$p,q$ variables are rational functions in complete elliptic integrals. Fuchsian
equations for off diagonal correlations $C(N,M)$ are given which extend our
considerations to discrete generalizations of Painlev{\'e}.

We present a simple, but efficient, way to calculate connection matrices
between sets of independent local solutions, defined at two neighboring
singular points, of Fuchsian differential equations of quite large orders, such
as those found for the third and fourth contribution ($\chi^{(3)}$ and
$\chi^{(4)}$) to the magnetic susceptibility of square lattice Ising model.
We use the previous connection matrices to get the exact explicit expressions
of all the monodromy matrices of the Fuchsian differential equation for
$\chi^{(3)}$ (and $\chi^{(4)}$) expressed in the same basis of solutions. These
monodromy matrices are the generators of the differential Galois group of the
Fuchsian differential equations for $\chi^{(3)}$ (and $\chi^{(4)}$), whose
analysis is just sketched here.