Search

This paper is a plea for diagonals and telescopers of rational, or algebraic, functions using creative telescoping, in a computer algebra experimental mathematics learn-by-examples approach. We show that diagonals of rational functions (and this is also the case with diagonals of algebraic functions) are left invariant when one performs an infinite set of birational transformations on the rational functions. These invariance results generalize to telescopers. We cast light on the almost systematic property of homomorphism to their adjoint of the telescopers of rational, or algebraic, functions. We shed some light on the reason why the telescopers, annihilating the diagonals of rational functions of the form P/Q^k and 1/Q, are homomorphic. For telescopers with solutions (periods) corresponding to integration over non-vanishing cycles, we have a slight generalization of this result. We introduce some challenging examples of generalization of diagonals of rational functions, like diagonals of transcendental functions, yielding simple $_2F_1$ hypergeometric functions associated with elliptic curves, or (differentially algebraic) lambda-extension of correlation of the Ising model., Comment: 48 pages

This paper provides several illustrations of the numerous remarkable properties of the lambda-extensions of the two-point correlation functions of the Ising model, sheding some light on the non-linear ODEs of the Painlev\'e type. We first show that this concept also exists for the factors of the two-point correlation functions focusing, for pedagogical reasons, on two examples namely C(0,5) and C(2,5) at $\nu = -k$. We then display, in a learn-by-example approach, some of the puzzling properties and structures of these lambda-extensions: for an infinite set of (algebraic) values of $ \lambda$ these power series become algebraic functions, and for a finite set of (rational) values of lambda they become D-finite functions, more precisely polynomials (of different degrees) in the complete elliptic integrals of the first and second kind K and E. For generic values of $ \lambda$ these power series are not D-finite, they are differentially algebraic. For an infinite number of other (rational) values of $ \lambda$ these power series are globally bounded series, thus providing an example of an infinite number of globally bounded differentially algebraic series. Finally, taking the example of a product of two diagonal two-point correlation functions, we suggest that many more families of non-linear ODEs of the Painlev\'e type remain to be discovered on the two-dimensional Ising model, as well as their structures, and in particular their associated lambda extensions. The question of their possible reduction, after complicated transformations, to Okamoto sigma forms of Painlev\'e VI remains an extremely difficult challenge., Comment: 23 pages

We revisit, with a pedagogical heuristic motivation, the lambda extension of the low-temperature row correlation functions C(M,N) of the two-dimensional Ising model. In particular, using these one-parameter series to understand the deformation theory around selected values of $\lambda$, namely $\lambda = \cos(\pi \, m/n)$ with m and n integers, we show that these series yield perturbation coefficients, generalizing form factors, that are D-finite functions. As a by-product these exact results provide an infinite number of highly non-trivial identities on the complete elliptic integrals of the first and second kind. These results underline the fundamental role of Jacobi theta functions and Jacobi forms, the previous D-finite functions being (relatively simple) rational functions of Jacobi theta functions, when rewritten in terms of the nome of elliptic functions., Comment: 35 pages

We study the factorizations of Ising low-temperature correlations C(M,N) for $\nu=-k$ and M+N odd, $M \le N$, for both the cases $M\neq 0$ where there are two factors, and $M=0$ where there are four factors. We find that the two factors for $ M \neq 0$ satisfy the same non-linear differential equation and, similarly, for M=0 the four factors each satisfy Okamoto sigma-form of Painlev\'e VI equations with the same Okamoto parameters. Using a Landen transformation we show, for $M\neq 0$, that the previous non-linear differential equation can actually be reduced to an Okamoto sigma-form of Painlev\'e VI equation. For both the two and four factor case, we find that there is a one parameter family of boundary conditions on the Okamoto sigma-form of Painlev\'e VI equations which generalizes the factorization of the correlations C(M,N) to an additive decomposition of the corresponding sigma's solutions of the Okamoto sigma-form of Painlev\'e VI equation which we call lambda extensions. At a special value of the parameter, the lambda-extensions of the factors of C(M,N) reduce to homogeneous polynomials in the complete elliptic functions of the first and second kind. We also generalize some Tracy-Widom (Painlev\'e V) relations between the sum and difference of sigma's to this Painlev\'e VI framework., Comment: 46 pages

We present Painlev{\'e} VI sigma form equations for the general Ising low and high temperature two-point correlation functions $ C(M,N)$ with $M \leq N $ in the special case $\nu = -k$ where $\nu = \, \sinh 2E_h/k_BT/\sinh 2E_v/k_BT$. More specifically four different non-linear ODEs depending explicitly on the two integers $M $ and $N$ emerge: these four non-linear ODEs correspond to distinguish respectively low and high temperature, together with $ M+N$ even or odd. These four different non-linear ODEs are also valid for $M \ge N$ when $ \nu = -1/k$. For the low-temperature row correlation functions $ C(0,N)$ with $ N$ odd, we exhibit again for this selected $\nu = \, -k$ condition, a remarkable phenomenon of a Painlev\'e VI sigma function being the sum of four Painlev\'e VI sigma functions having the same Okamoto parameters. We show in this $\nu = \, -k$ case for $ T < T_c $ and also $ T > T_c$, that $ C(M,N)$ with $ M \leq N $ is given as an $ N \times N$ Toeplitz determinant.

We show that the unresolved examples of Christol's conjecture $ \, _3F_{2}\left([2/9,5/9,8/9],[2/3,1],x\right)$ and $_3F_{2}\left([1/9,4/9,7/9],[1/3,1],x\right)$, are indeed diagonals of rational functions. We also show that other $\, _3F_2$ and $\, _4F_3$ unresolved examples of Christol's conjecture are diagonals of rational functions. Finally we give two arguments that show that it is likely that the $\, _3F_2([1/9, 4/9, 5/9], \, [1/3,1], \, 27 \cdot x)$ function is a diagonal of a rational function., Comment: 13 pages

We provide a set of diagonals of simple rational functions of three and four variables that are squares of Heun functions. These Heun functions obtained through creative telescoping, turn out to be either pullbacked $_2F_1$ hypergeometric functions and in fact classical modular forms. We also obtain Heun functions that are Shimura curves as solutions of telescopers of rational functions., Comment: 64 pages

We recall that diagonals of rational functions naturally occur in lattice statistical mechanics and enumerative combinatorics. We find that a seven-parameter rational function of three variables with a numerator equal to one (reciprocal of a polynomial of degree two at most) can be expressed as a pullbacked 2F1 hypergeometric function. This result can be seen as the simplest non-trivial family of diagonals of rational functions. We focus on some subcases such that the diagonals of the corresponding rational functions can be written as a pullbacked 2F1 hypergeometric function with two possible rational functions pullbacks algebraically related by modular equations, thus showing explicitely that the diagonal is a modular form. We then generalise this result to eight, nine and ten parameters families adding some selected cubic terms at the denominator of the rational function defining the diagonal. We finally show that each of these previous rational functions yields an infinite number of rational functions whose diagonals are also pullbacked 2F1 hypergeometric functions and modular forms., Comment: 39 pages

We show that non-linear Schwarzian differential equations emerging from covariance symmetry conditions imposed on linear differential operators with hypergeometric function solutions, can be generalized to arbitrary order linear differential operators with polynomial coefficients having selected differential Galois groups. For order three and order four linear differential operators we show that this pullback invariance up to conjugation eventually reduces to symmetric powers of an underlying order-two operator. We give, precisely, the conditions to have modular correspondences solutions for such Schwarzian differential equations, which was an open question in a previous paper. We analyze in detail a pullbacked hypergeometric example generalizing modular forms, that ushers a pullback invariance up to operator homomorphisms. We expect this new concept to be well-suited in physics and enumerative combinatorics. We finally consider the more general problem of the equivalence of two different order-four linear differential Calabi-Yau operators up to pullbacks and conjugation, and clarify the cases where they have the same Yukawa couplings., Comment: 47 pages

We discuss the implications of studies of partition function zeros and equimodular curves for the analytic properties of the Ising model on a square lattice in a magnetic field. In particular we consider the dense set of singularities in the susceptibility of the Ising model at $H=0$ found by Nickel and its relation to the analyticity of the field theory computations of Fonseca and Zamolodchikov., Comment: 21 pages, 13 figures

We give examples of infinite order rational transformations that leave linear differential equations covariant. These examples are non-trivial yet simple enough illustrations of exact representations of the renormalization group. We first illustrate covariance properties on order-two linear differential operators associated with identities relating the same $_2F_1$ hypergeometric function with different rational pullbacks. We provide two new and more general results of the previous covariance by rational functions: a new Heun function example and a higher genus $_2F_1$ hypergeometric function example. We then focus on identities relating the same hypergeometric function with two different algebraic pullback transformations: such remarkable identities correspond to modular forms, the algebraic transformations being solution of another differentially algebraic Schwarzian equation that emerged in a paper by Casale. Further, we show that the first differentially algebraic equation can be seen as a subcase of the last Schwarzian differential condition, the restriction corresponding to a factorization condition of some associated order-two linear differential operator. Finally, we also explore generalizations of these results, for instance, to $_3F_2$, hypergeometric functions, and show that one just reduces to the previous $_2F_1$ cases through a Clausen identity. In a $_2F_1$ hypergeometric framework the Schwarzian condition encapsulates all the modular forms and modular equations of the theory of elliptic curves, but these two conditions are actually richer than elliptic curves or $_2F_1$ hypergeometric functions, as can be seen on the Heun and higher genus example. This work is a strong incentive to develop more differentially algebraic symmetry analysis in physics., Comment: 43 pages

We study the class of non-holonomic power series with integer coefficients that reduce, modulo primes, or powers of primes, to algebraic functions. In particular we try to determine whether the susceptibility of the square-lattice Ising model belongs to this class, and more broadly whether the susceptibility is a solution of a differentially algebraic equation. Initial results on Tutte's non-linear ordinary differential equation (ODE) and other simple quadratic non-linear ODEs suggest that a large set of differentially algebraic power series solutions with integer coefficients might reduce to algebraic functions modulo primes, or powers of primes. Here we give several examples of series with integer coefficients and non-zero radius of convergence that reduce to algebraic functions modulo (almost) every prime (or power of a prime). These examples satisfy differentially algebraic equations with the encoding polynomial occasionally possessing quite high degree (and thus difficult to identify even with long series). Additionally, we have extended both the high- and low-temperature Ising square-lattice susceptibility series to 5043 coefficients. We find that even this long series is insufficient to determine whether it reduces to algebraic functions modulo $3$, $5$, etc. This negative result is in contrast to the comparatively easy confirmation that the corresponding series reduce to algebraic functions modulo powers of $2$., Comment: 36 pages

We present the reduction of the correlation functions of the Ising model on the anisotropic square lattice to complete elliptic integrals of the first, second and third kind, the extension of Kramers-Wannier duality to anisotropic correlation functions, and the linear differential equations for these anisotropic correlations. More precisely, we show that the anisotropic correlation functions are homogeneous polynomials of the complete elliptic integrals of the first, second and third kind. We give the exact dual transformation matching the correlation functions and the dual correlation functions. We show that the linear differential operators annihilating the general two-point correlation functions are factorised in a very simple way, in operators of decreasing orders., Comment: 22 pages

We consider the isotropic perimeter generating functions of three-choice, imperfect, and 1-punctured staircase polygons, whose 8th order linear Fuchsian ODEs are previously known. We derive simple relationships between the three generating functions, and show that all three generating functions are joint solutions of a common 12th order Fuchsian linear ODE. We find that the 8th order differential operators can each be rewritten as a direct sum of a direct product, with operators no larger than 3rd order. We give closed-form expressions for all the solutions of these operators in terms of $_2F_1$ hypergeometric functions with rational and algebraic arguments. The solutions of these linear differential operators can in fact be expressed in terms of two modular forms, since these $_2F_1$ hypergeometric functions can be expressed with two, rational or algebraic, pullbacks., Comment: 28 pages

We previously reported on a recursive method to generate the expansion of the lattice Green function of the $d$-dimensional face-centred 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 2-dimensional array $ T_d(n,j)$ where only the coefficients $ t_d(n,0)$ correspond to the coefficients of the lattice Green function of the d-dimensional fcc. The coefficients $ t_d(n,j)$ are associated to D-finite bivariate series annihilated by linear partial differential equations that we analyze., Comment: 28 pages

We recall that the full susceptibility series of the Ising model, modulo powers of the prime 2, reduce to algebraic functions. We also recall the non-linear polynomial differential equation obtained by Tutte for the generating function of the q-coloured rooted triangulations by vertices, which is known to have algebraic solutions for all the numbers of the form $2 +2 \cos(j\pi/n)$, the holonomic status of the q= 4 being unclear. We focus on the analysis of the q= 4 case, showing that the corresponding series is quite certainly non-holonomic. Along the line of a previous work on the susceptibility of the Ising model, we consider this q=4 series modulo the first eight primes 2, 3, ... 19, and show that this (probably non-holonomic) function reduces, modulo these primes, to algebraic functions. We conjecture that this probably non-holonomic function reduces to algebraic functions modulo (almost) every prime, or power of prime numbers. This raises the question to see whether such remarkable non-holonomic functions can be seen as ratio of diagonals of rational functions, or algebraic, functions of diagonals of rational functions., Comment: 27 pages

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

We study the full susceptibility of the Ising model modulo powers of primes. We find exact functional equations for the full susceptibility modulo these primes. Revisiting some lesser-known results on discrete finite automata, we show that these results can be seen as a consequence of the fact that, modulo 2^r, one cannot distinguish the full susceptibility from some simple diagonals of rational functions which reduce to algebraic functions modulo 2^r, and, consequently, satisfy exact functional equations modulo 2^r. We sketch a possible physical interpretation of these functional equations modulo 2^r as reductions of a master functional equation corresponding to infinite order symmetries such as the isogenies of elliptic curves. One relevant example is the Landen transformation which can be seen as an exact generator of the Ising model renormalization group. We underline the importance of studying a new class of functions corresponding to ratios of diagonals of rational functions: they reduce to algebraic functions modulo powers of primes and they may have solutions with natural boundaries., Comment: 21 pages

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 two-point 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., Comment: 26 pages

We revisit an order-six 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 order-two self-adjoint operators. From these results, we deduce a new form of decomposition of operators for this selected order-six linear differential operator in terms of three order-two self-adjoint operators. We then generalize the previous decomposition to decompositions in terms of an arbitrary number of self-adjoint 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 self-adjoint 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 4-polytopes defining Calabi-Yau 3-folds, obtained recently by P. Lairez, correspond to a particular form of the decomposition detailed in this paper., Comment: 40 pages

In this paper we compare the integrable hard hexagon model with the non-integrable hard squares model by means of partition function roots and transfer matrix eigenvalues. We consider partition functions for toroidal, cylindrical, and free-free boundary conditions up to sizes $40\times40$ and transfer matrices up to 30 sites. For all boundary conditions the hard squares roots are seen to lie in a bounded area of the complex fugacity plane along with the universal hard core line segment on the negative real fugacity axis. The density of roots on this line segment matches the derivative of the phase difference between the eigenvalues of largest (and equal) moduli and exhibits much greater structure than the corresponding density of hard hexagons. We also study the special point $z=-1$ of hard squares where all eigenvalues have unit modulus, and we give several conjectures for the value at $z=-1$ of the partition functions., Comment: 46 pages

We show that the globally nilpotent G-operators 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 self-adjoint 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 self-adjoint (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 order-twelve and order-23 operators. We show that the order-twelve 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 order-23 operator is shown to factorize in an order-two operator and an order-21 operator. The symmetric square of this order-21 operator has a rational solution. Its differential Galois group is, thus, in the orthogonal group $ SO(21, \mathbb{C})$., Comment: 33 pages

We study the analyticity of the partition function of the hard hexagon model in the complex fugacity plane by computing zeros and transfer matrix eigenvalues for large finite size systems. We find that the partition function per site computed by Baxter in the thermodynamic limit for positive real values of the fugacity is not sufficient to describe the analyticity in the full complex fugacity plane. We also obtain a new algebraic equation for the low density partition function per site., Comment: 49 pages, IoP styles files, lots of figures (png mostly) so using PDFLaTeX. Some minor changes added to version 2 in response to referee reports

We show that the n-fold integrals $\chi^{(n)}$ of the magnetic susceptibility of the Ising model, as well as various other n-fold integrals of the "Ising class", or n-fold 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 Calabi-Yau ODEs, and, more generally, Picard-Fuchs 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., Comment: This paper is the short version of the larger (100 pages) version, available as arXiv:1211.6031 , where all the detailed proofs are given and where a much larger set of examples is displayed

We show that the n-fold integrals $\chi^{(n)}$ of the magnetic susceptibility of the Ising model, as well as various other n-fold integrals of the "Ising class", or n-fold 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 Calabi-Yau ODEs, and more generally of MUM ODEs, is, almost always, diagonal of rational functions. We revisit Christol's conjecture that globally bounded series of G-operators 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 Calabi-Yau 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., Comment: 100 pages

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 Picard-Fuchs systems of two-variables "above" Calabi-Yau 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 n-fold 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 two-parameter family of elliptic curves. We then find a first set of non-Nickelian 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 Calabi-Yau manifolds on such problems. We discuss the accumulation of these singular curves for the non-holonomic anisotropic full susceptibility., Comment: 36 pages

Understanding the relationship which integrable (solvable) models, all of which possess very special symmetry properties, have with the generic non-integrable models that are used to describe real experiments, which do not have the symmetry properties, is one of the most fundamental open questions in both statistical mechanics and quantum field theory. The importance of the two-dimensional Ising model in a magnetic field is that it is the simplest system where this relationship may be concretely studied. We here review the advances made in this study, and concentrate on the magnetic susceptibility which has revealed an unexpected natural boundary phenomenon. When this is combined with the Fermionic representations of conformal characters, it is suggested that the scaling theory, which smoothly connects the lattice with the correlation length scale, may be incomplete for $H \neq 0$., Comment: 33 pages

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 Calabi-Yau 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 order-six 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" Ising-operators, are "special" ones, reducing, in fact systematically, to (selected, k-balanced, ...) $_{q+1}F_q$ hypergeometric functions, or correspond to the more general solutions of Calabi-Yau equations., Comment: 35 pages

We present a general method for analytically factorizing the n-fold form factor integrals $f^{(n)}_{N,N}(t)$ for the correlation functions of the Ising model on the diagonal in terms of the hypergeometric functions $_2F_1([1/2,N+1/2];[N+1];t)$ which appear in the form factor $f^{(1)}_{N,N}(t)$. New quadratic recursion and quartic identities are obtained for the form factors for n=2,3. For n= 2,3,4 explicit results are given for the form factors. These factorizations are proved for all N for n= 2,3. These results yield the emergence of palindromic polynomials canonically associated with elliptic curves. As a consequence, understanding the form factors amounts to describing and understanding an infinite set of palindromic polynomials, canonically associated with elliptic curves. From an analytical viewpoint the relation of these palindromic polynomials with hypergeometric functions associated with elliptic curves is made very explicitly, and from a differential algebra viewpoint this corresponds to the emergence of direct sums of differential operators homomorphic to symmetric powers of a second order operator associated with elliptic curve., Comment: 37 pages

We show that almost all the linear differential operators factors obtained in the analysis of the n-particle 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 order-four 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 Calabi-Yau 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 pull-back, the corresponding Calabi-Yau fourth-order 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., Comment: 46 pages

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$., Comment: 19 pages, 1 figure

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., Comment: 25 pages, 2 figures, IoP style files

This paper deals with $\tilde{\chi}^{(6)}$, the six-particle 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 left-most factor of order six being equivalent to the symmetric fifth power of the linear differential operator corresponding to the elliptic integral $E$. The right-most 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., Comment: 23 pages

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

We consider some two-dimensional 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 two-form. 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 two-form. 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., Comment: 26 pages, 11 figures

We consider the Fuchsian linear differential equation obtained (modulo a prime) for $\tilde{\chi}^{(5)}$, the five-particle 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 left-most 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 $(n-1)$-th power of $L_E$ occurs as a left-most factor in the minimal order linear differential operators for all $\tilde{\chi}^{(n)}$'s., Comment: 33 pages

We recall various multiple integrals related to the isotropic square Ising model, and corresponding, respectively, to the n-particle contributions of the magnetic susceptibility, to the (lattice) form factors, to the two-point correlation functions and to their lambda-extensions. These integrals are holonomic and even G-functions: 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 p-curvature. 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 weight-1 modular form emerging in the three-particle contribution $ \chi^{(3)}$ of the magnetic susceptibility of the square Ising model. In the case where we do not have G-functions, but Hamburger functions (one irregular singularity at 0 or $ \infty$) that correspond to the confluence of singularities in the scaling limit, the p-curvature is also found to verify new structures associated with simple deformations of the nilpotent property., Comment: 55 pages

We calculate very long low- and high-temperature series for the susceptibility $\chi$ of the square lattice Ising model as well as very long series for the five-particle contribution $\chi^{(5)}$ and six-particle 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 diff-Pad\'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 series-solutions 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 power-law 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$., Comment: 54 pages, 2 figures

We study birational transformations of the projective space originating from lattice statistical mechanics, specifically from various chiral Potts models. Associating these models to \emph{stable patterns} and \emph{signed-patterns}, we give general results which allow us to find \emph{all} chiral $q$-state spin-edge Potts models when the number of states $q$ is a prime or the square of a prime, as well as several $q$-dependent family of models. We also prove the absence of monocolor stable signed-pattern with more than four states. This demonstrates a conjecture about cyclic Hadamard matrices in a particular case. The birational transformations associated to these lattice spin-edge models show complexity reduction. In particular we recover a one-parameter family of integrable transformations, for which we give a matrix representation, Comment: 22 pages 0 figure The paper has been reorganized, splitting the results into two sections : results pertaining to Physics and results pertaining to Mathematics

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., Comment: 39 pages, 7 figures

We use the recently derived form factor expansions of the diagonal two-point 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 russian-doll 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$., Comment: 18 pages

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 self-dual 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 Kramers-Wannier duality. We revisit the singular behavior (J. Phys. A {\bf 38} (2005) 9439-9474) 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$)., Comment: 34 pages, 1 figure

We briefly recall the Fuchs-Painlev\'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 Toeplitz-like determinants, are equivalent to the second order operator $ L_E$ which has $ E$ as solution (or, for off-diagonal 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 singled-out second order linear differential operator $ L_E$ being replaced by an isomonodromic system of two third-order 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 non-linear 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., Comment: Dedicated to the : Special issue on Symmetries and Integrability of Difference Equations, SIDE VII meeting held in Melbourne during July 2006

We study the Ising model two-point 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 two-point 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 singled-out 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., Comment: 39 pages

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}., Comment: 18 pages, Dedicated to the centenary of the publication of the Painleve VI equation in the Comptes Rendus de l'Academie des Sciences de Paris by Richard Fuchs in 1905

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., Comment: 10 pages 2005 Nankai conference on differential geometry

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 deduce all the critical behaviors of the solutions $\chi^{(3)}$ and $\chi^{(4)}$, as well as the asymptotic behavior of the coefficients in the corresponding series expansions. We confirm that the newly found quadratic number singularities of the Fuchsian ODE associated to $\chi^{(3)}$ are not singularities of the particular solution $\chi^{(3)}$ itself. We use the previous connection matrices to get the exact 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. As far as the physics implications of the solutions are concerned, we find challenging qualitative differences when comparing the corrections to scaling for the full susceptibillity $\chi$ at high temperature (resp. low temperature) and the first two terms $\chi^{(1)}$ and $\chi^{(3)}$ (resp. $\chi^{(2)}$ and $\chi^{(4)}$) ., Comment: 38 pages, no figure Note added in the proofs

Using an expansion method in the variables $ x_i$ that appear in the $(n-1)$-dimensional integrals representing the $n$-particle contribution to the Ising square lattice model susceptibility $\chi$, we generate a long series of coefficients for the 3-particle contribution $\chi^{(3)}$, using a $ N^4$ polynomial time algorithm. We give the Fuchsian differential equation of order seven for $\chi^{(3)}$ that reproduces all the terms of our long series. An analysis of the properties of this Fuchsian differential equation is performed., Comment: 15 pages, no figures, submitted to J. Phys. A

We perform a Random Matrix Theory (RMT) analysis of the quantum four-state chiral Potts chain for different sizes of the chain up to size L=8. Our analysis gives clear evidence of a Gaussian Orthogonal Ensemble statistics, suggesting the existence of a generalized time-reversal invariance. Furthermore a change from the (generic) GOE distribution to a Poisson distribution occurs when the integrability conditions are met. The chiral Potts model is known to correspond to a (star-triangle) integrability associated with curves of genus higher than zero or one. Therefore, the RMT analysis can also be seen as a detector of ``higher genus integrability''., Comment: 23 pages and 10 figures

We will try to sketch Professor F. Y. Wu's contributions in lattice statistical mechanics, solid state physics, graph theory, enumerative combinatorics and so many other domains of physics and mathematics. We will recall F. Y. Wu's most important and well-known classic results and we will also sketch his most recent researches dedicated to the connections of lattice statistical mechanical models with deep problems in pure mathematics. Since it is hard to provide an exhaustive list of all his contributions, to give some representation of F. Y. Wu's "mental connectivity" we will concentrate on the interrelations between the various results he has obtained in so many different domains of physics and mathematics. Along the way we will also try to understand Wu's motivations and his favorite concepts, tools and ideas., Comment: One LaTex file, 63 pages In honor of F.Y. Wu on the occasion of his 70th birthday. Statphys-Taiwan 2002. The second APCTP and sixth Taiwan International Symposium on Statistical Physics. Table of contents + minor correction