Elliptic Painlev\'e equations from next-nearest-neighbor translations on the $E_8^{(1)}$ lattice
Nalini Joshi, Nobutaka Nakazono

TL;DR
This paper introduces a new elliptic discrete Painlevé equation derived from next-nearest-neighbor translations on the $E_8^{(1)}$ lattice, expanding the class of known elliptic difference equations.
Contribution
It presents a novel elliptic Painlevé equation from next-nearest-neighbor lattice translations, generalizing previous models and connecting to Adler's Q4 equation.
Findings
New elliptic Painlevé equation with 8 parameters
Reduction of Sakai's equation to a projective form
Connection to Adler's Q4 equation
Abstract
The well known elliptic discrete Painlev\'e equation of Sakai is constructed by a standard translation on the lattice, given by nearest neighbor vectors. In this paper, we give a new elliptic discrete Painlev\'e equation obtained by translations along next-nearest-neighbor vectors. This equation is a generic (8-parameter) version of a 2-parameter elliptic difference equation found by reduction from Adler's partial difference equation, the so-called Q4 equation. We also provide a projective reduction of the well known equation of Sakai.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Elliptic Painlevé equations from next-nearest-neighbor translations on the lattice
Nalini Joshi
and
Nobutaka Nakazono
School of Mathematics and Statistics F07, The University of Sydney, New South Wales 2006, Australia.
Abstract.
The well known elliptic discrete Painlevé equation of Sakai is constructed by a standard translation on the lattice, given by nearest neighbor vectors. In this paper, we give a new elliptic discrete Painlevé equation obtained by translations along next-nearest-neighbor vectors. This equation is a generic (8-parameter) version of a 2-parameter elliptic difference equation found by reduction from Adler’s partial difference equation, the so-called Q4 equation. We also provide a projective reduction of the well known equation of Sakai.
Key words and phrases:
discrete Painlevé equation; elliptic Painlevé equation; affine Weyl group; elliptic function; projective reduction
2010 Mathematics Subject Classification:
14H70, 33E05, 33E17,
1. Introduction
At the head of the list of discrete Painlevé equations described by Sakai sits an elliptic difference equation, which has attracted a great deal of attention in recent times. This equation [22, 17] has the affine Weyl group symmetry of type and a perennial question is whether it is unique as the only elliptic-difference-type discrete Painlevé equation. While two other candidate equations are now known [18, 20], these are related to Sakai’s elliptic difference equation. In this paper, we deduce elliptic-difference equations for the first time, which are not related to Sakai’s equation by either a Miura transformation or by projective reduction. We do this by considering different (non-conjugate) translations on the lattice and describing these translations under the action of Jacobian elliptic functions.
A discrete Painlevé equation is an ordinary difference equation, which is iterated by translation in an affine Weyl group, where translation corresponds to vectors in the root- or (equivalently) weight-lattice [4, 11]. Sakai’s elliptic difference equation is realized by translation on the lattice.
To describe the construction, we first explain how such translations are characterized. Fix a point in the lattice. Then there are 240 nearest neighbors of this point in the lattice, lying at a distance whose squared length is equal to 2. We refer to the 120 vectors between this fixed point and its possible nearest neighbors as nearest-neighbor-connecting vectors (NVs). Similarly, there are 2160 next-nearest neighbors, lying at a distance whose squared length is 4. The 1080 vectors between the fixed point and such next-nearest neighbors will be referred to as next-nearest-neighbor-connecting vectors (NNVs). In [22] and [18], elliptic difference equations were constructed as translations expressed in terms of NVs.
In a previous study [2], we considered another elliptic difference equation (Equation (1.1)), which was obtained by reduction from a partial difference equation [20]. Curiously, its symmetry group turns out to be , a sub-group of . Moreover, its time iteration turns out not to be given by translation on the lattice. However, its square (i.e., composition with itself) is a translation. We have discovered that this translation is expressed in terms of NNVs, which makes it very different to the elliptic difference equations in earlier papers [22, 18]. Note that independent study by Carstea et al. [3] has also led to equations with translations given by NNVs.
Generalizing this surprising insight led us to the discovery of a new elliptic Painlevé equation (3.36), which can be regarded as the generic version of Equation (1.1). The term generic here refers to the fact that it contains the largest number of parameters possible as an equation with symmetry group . It has 8 parameters in addition to the independent variable. Moreover, we also obtain a projective reduction of Sakai’s elliptic difference equation (3), which contains 4 parameters in addition to the independent variable.
All of these elliptic difference equations have the same space of initial values, regularized by blowing up 8 points in arbitrary position in . The space of initial values is labelled as . We emphasize here that this characterization of the space of initial values is not enough to distinguish the four different elliptic difference equations described above. We provide evidence that different (non-conjugate) translations in \widetilde{W}\big{(}E_{8}^{(1)}\big{)} lead to distinct elliptic difference equations. All four elliptic difference equations are described in this paper in terms of Jacobian elliptic functions for the sake of uniformity, because this is the form originally given for Equation (1.1).
1.1. Background
Discrete Painlevé equations are nonlinear ordinary difference equations of second order, which include discrete analogues of the six Painlevé equations: PI, , PVI. Together with the Painlevé equations, the discrete Painlevé equations are now regarded as one of the most important classes of equations in the theory of integrable systems (see, e.g., [8]).
Sakai’s geometric description of discrete Painlevé equations, based on types of space of initial values, is well known[22]. This picture relies on compactifying and regularizing space of initial values. The spaces of initial values are constructed by the blow up of at eight base points (i.e. points where the system is ill defined because it approaches ) and are classified into 22 types according to the configuration of the base points as follows:
[TABLE]
In each case, the root system characterizing the surface forms a subgroup of the 10-dimensional Picard group. The symmetry group of each equation, formed by Cremona isometries, arises from the orthogonal complement of this root system inside the Picard group. Its birational actions give the discrete Painlevé equation of interest in each case.
Recently, the following elliptic Painlevé equation was obtained from the reduction of Adler’s partial difference equation (or, Q4 equation)[20]:
[TABLE]
where is the modulus of the elliptic sine,
[TABLE]
and
[TABLE]
See Appendix A for standard results about Jacobian elliptic functions.
In [2], the geometry of Equation (1.1), i.e., its space of initial values and corresponding Cremona isometries, was investigated. The space of initial values of Equation (1.1) was identified with the elliptic -surface and its Cremona isometries collectively form an affine Weyl transformation group of type , denoted by W\big{(}E_{8}^{(1)}\big{)}. Moreover, it was shown that Equation (1.1) cannot be derived from a translation of W\big{(}E_{8}^{(1)}\big{)} but can be derived by using a projective reduction. The process of deriving discrete dynamical systems of Painlevé type from elements of affine Weyl groups that are of infinite order (but that are not necessarily translations) by taking a projection on an appropriate subspace of parameters is called a projective reduction [13, 12]. Note that this process is motivated by taking symmetric versions of systems of discrete Painlevé equations, which has been known since [21].
Although the geometry of Equation (1.1) has been clarified, the realization of Equation (1.1) from the action of Cremona isometries was missing since its base points are parametrized by the Jacobian elliptic function (Jacobi’s setting) and birational actions of Cremona isometries on such setting were not explicitly known.
The present study fills this gap (see Theorem 3.1), that is, we provide the realization of the Equation (1.1) as a half-translation of the extended affine Weyl group \widetilde{W}\big{(}E_{8}^{(1)}\big{)}.
In the remainder of the paper, we refer to Sakai’s elliptic Painlevé equation [22] as the MSY elliptic Painlevé equation [17] because the former was obtained in while the latter was provided by resolution in . We work throughout the paper in .
1.2. Plan of the paper
This paper is organized as follows. In §2, we construct Cremona isometries for base points (2.4) that generalize those of Equation (1.1) and show that these form an affine Weyl group of type , denoted by W\big{(}E_{8}^{(1)}\big{)}, under the linear actions on the Picard group. In §3, using intersection theory, we obtain the birational action of W\big{(}E_{8}^{(1)}\big{)} on the coordinates and parameters of the base points, and we prove that these birational actions also satisfy the fundamental relations of an affine Weyl group of type . Adding the identity mappings on the Picard group, we obtain the extended affine Weyl group \widetilde{W}\big{(}E_{8}^{(1)}\big{)}. Equation (1.1) and three other elliptic Painlevé equations are then deduced from the action of this resulting group \widetilde{W}\big{(}E_{8}^{(1)}\big{)}. Finally, we give some concluding remarks in §4.
2. Cremona isometries
In this section, we construct Cremona isometries for a generalization of the base points of Equation (1.1) and show that these isometries collectively form an affine Weyl group of type .
Recall that Equation (1.1) has the following eight base points (see [2]):
[TABLE]
where and are complete elliptic integrals and
[TABLE]
which lie on the elliptic curve
[TABLE]
Note that these base points contain three arbitrary parameters , , .
We now generalize the number of arbitrary parameters. Take a new set of base points given by
[TABLE]
where , , and are non-zero complex parameters. These generalized base points lie on the elliptic curve
[TABLE]
These base points can be reduced to the ones given in Equation (2.1) under a specialization of the parameters. Indeed, the points (2.4) and elliptic curve (2.5) can be respectively reduced to the points (2.1) and the curve (2.3) by assuming
[TABLE]
Remark 2.1**.**
To distinguish between Weierstrass’s and Jacobi’s setting, we here denote the surface characterized by the elliptic curve (2.5) as -surface. The relation between - and -surface is given in Appendix B.
Let denote the blow up of at the points (2.4). Moreover, let the linear equivalence classes of the total transform of vertical and horizontal lines in be denoted respectively by and .
The Picard group of , denoted by Pic, is given by
[TABLE]
where , , are exceptional divisors. The intersection form is given by a symmetric bilinear form with
[TABLE]
where is the Kronecker delta. The anti-canonical divisor of is given by
[TABLE]
For later convenience, let
[TABLE]
We define the root lattice by the elements of Pic that are orthogonal to the anti-canonical divisor . The simple roots , , are given by
[TABLE]
where
[TABLE]
We can easily verify that
[TABLE]
Representing intersecting and by a line between nodes and , we obtain the Dynkin diagram of shown in Figure 1.
Definition 2.2**.**
[7, 6, 15]* An automorphism of Pic is called a Cremona isometry if it preserves*
**(i): **
the intersection form on Pic;
**(ii): **
the canonical divisor ;
**(iii): **
effectiveness of each effective divisor of Pic.
It is well-known that the reflections are Cremona isometries. In this case we define the reflections , , by the following linear actions:
[TABLE]
for all .
They collectively form an affine Weyl group of type , denoted by W\big{(}E_{8}^{(1)}\big{)}. Namely, we can easily verify that under the action on the the following fundamental relations hold:
[TABLE]
where
[TABLE]
3. Birational actions of the Cremona isometries for the Jacobi’s setting
In this section, we give the birational actions of the Cremona isometries on the coordinates and parameters of the base points (2.4). By using these birational actions, we derive various elliptic Painlevé equations.
We focus on a particular example first to explain how to deduce such birational actions. Recall and are given by the linear equivalence classes of vertical lines and horizontal lines , respectively. Applying the reflection operator given by (2.14) to , we find
[TABLE]
which means that can be described by the curve of bi-degree passing through base points and with multiplicity . (See [14] for for more detail.) This result leads us to the birational action given below in Equation (3.2b).
Similarly, from the linear actions of , , we obtain their birational actions on the coordinates and parameters of the base points (2.4) as follows. The actions of the generators of W\big{(}E_{8}^{(1)}\big{)} on the coordinates are given by
[TABLE]
Note that
[TABLE]
is invariant under the action of W\big{(}E_{8}^{(1)}\big{)}.
For Jacobi’s elliptic function it is well known that shifts by half periods give the following relations:
[TABLE]
These identities motivate our search for the transformations that are identity mappings on the . Indeed, we define such transformations , , by the following actions:
[TABLE]
Adding the transformations , we extend W\big{(}E_{8}^{(1)}\big{)} to
[TABLE]
In general, for a function , we let an element w\in\widetilde{W}\big{(}E_{8}^{(1)}\big{)} act as , that is, acts on the arguments from the left.
Theorem 3.1**.**
Under the birational actions (3.2) and (3.5), the generators of \widetilde{W}\big{(}E_{8}^{(1)}\big{)}=\langle\iota_{1},\iota_{2},\iota_{3},\iota_{4}\rangle\rtimes\langle s_{0},\dots,s_{8}\rangle satisfy the following fundamental relations:
[TABLE]
where
[TABLE]
Proof.
Let us define the function by
[TABLE]
The action (3.2b) can be expressed as
[TABLE]
By replacing the Jacobian elliptic function to the Jacobi theta-functions using (A.4) and then applying the addition formula (A.9), the function can be rewritten as
[TABLE]
which gives
[TABLE]
Therefore, the action (3.2b) can be also expressed as
[TABLE]
which leads the proof of the relation . The other relations can be directly verified by using the actions (3.2) and (3.5). Therefore, we have completed the proof. ∎
Now we are in a position to derive Equation (1.1) from the Cremona transformations associated with -surface. Note that for convenience we use the following notations for the composition of the reflections and for the summation of the parameters :
[TABLE]
respectively. Let
[TABLE]
The actions of on the root lattice and parameter space are not translational as the following shows:
[TABLE]
where is defined by (2.2). However, when the parameters take special values (2.6), the action of becomes the translational motion in the parameter subspace:
[TABLE]
and then the action on the coordinates:
[TABLE]
gives Equation (1.1). Here, we also consider the translation
[TABLE]
whose actions on the root lattice and parameter space are given by
[TABLE]
The action on the coordinates:
[TABLE]
gives the following elliptic Painlevé equation:
[TABLE]
Here, for conciseness, we provide the terms , and in Appendix A.
Remark 3.2**.**
The periodicity of Jacobian elliptic functions in Equations (1.1) and (3.36) allows us to take the actions (3.30) and (3.34) to be
[TABLE]
respectively, without loss of generality.
Remark 3.3**.**
The relation (3.32) is the key to the reduction of Equation (3.36) to Equation (1.1). Under the conditions (2.6), it gives the projective reduction of Equation (3.36) to Equation (1.1). (See [13, 12] for more detail.) We therefore refer to Equation (1.1) as the projectively-reduced equation of Equation (3.36). Conversely, we refer to Equation (3.36) as the generic version of Equation (1.1). Note that the word “generic” here has a specific meaning. An equation is called generic only when it has the same number of parameters as the corresponding surface-type discrete Painlevé equation in Sakai’s list [22].
Remark 3.4**.**
Often a discrete Painlevé equation can be extended up to a generic version by singularity confinement [9, 21]. However, Equation (1.1) has a reduced number of parameters after applying the singularity confinement criterion (see [20]). The reason is some of the coefficients in Equation (3.36) become zero under the condition (2.6). In such a case, singularity confinement cannot be used to extend the equation up to the generic version. We note that a similar observation has also been reported for another example in [10].
Next, we provide a new expression for the MSY elliptic Painlevé equation investigated in [22, 16, 17, 19, 18]. Moreover, we also construct its projectively-reduced equation.
Define a translation of \widetilde{W}\big{(}E_{8}^{(1)}\big{)} by
[TABLE]
where
[TABLE]
whose action on the root lattice is given by
[TABLE]
The translation actually corresponds to the translation in [17]. However, it was expressed in terms of Weierstrass elliptic functions. We now convert this translation to Jacobian elliptic functions by using the birational actions (3.2). Consider the following action of :
[TABLE]
Now consider the special choices of parameters:
[TABLE]
Then, the transformation acts on this parameter subspace and variables , as
[TABLE]
This leads to a scalar second-order mapping for alone, given by
[TABLE]
which is the projectively-reduced equation of the MSY elliptic Painlevé equation (3.42).
Remark 3.5**.**
Recall that a Kac translation is defined by
[TABLE]
and, therefore, its squared length is given by .
In [18], Ohta et al. discovered that NVs are equivalent to each other under the operations of the Weyl group W\big{(}E_{8}^{(1)}\big{)}=\langle s_{0},\dots,s_{8}\rangle and the same holds for NNVs. and are expressed by Kac transformations as the following:
[TABLE]
Since
[TABLE]
* and are a NNV and a NV, respectively. In this sense, the translation differs from .*
Remark 3.5 shows that Equation (3.36) forms a new system of elliptic difference equations, which differs from the MSY elliptic difference equation given originally by Sakai. Since the generic versions of equations are different, Equation (3) differs from Equation (1.1). In this sense, Equation (3) is also a new elliptic difference equation.
4. Concluding remarks
In this paper, we have constructed the birational actions of the Cremona isometries \widetilde{W}\big{(}E_{8}^{(1)}\big{)} for the -surface. Equation (1.1) is realized in terms of birational actions of Cremona isometries for the first time. Using birational actions, we further derived the elliptic Painlevé equations (3.36), (3.42) and (3).
Equation (3.42) is a new expression of the important elliptic Painlevé equation studied in [22, 16, 17, 19, 18]. For a long time, this equation remained the unique elliptic Painlevé equation known with 8 parameters. In this paper, we provided a new elliptic Painlevé equation which also has 8 parameters. We showed that it must be different because it is realized in terms of translations in \widetilde{W}\big{(}E_{8}^{(1)}\big{)}, which are not conjugate to those that give rise to the previously known elliptic difference equation. Furthermore, we provided a projectively-reduced version of the MSY elliptic Painlevé equation, namely Equation (3). We also showed that Equation (3.36) is a generic version of Equation (1.1).
These results lead to more questions. The search for Lax pairs of these new equations remains open. Moreover, special function solutions of the projectively-reduced equations (3) and (1.1) remain to be found.
Acknowledgment
The authors would like to express their sincere thanks to Profs M. Noumi, Y. Ohta, T. Takenawa and Y. Yamada for inspiring and fruitful discussions. This research was supported by an Australian Laureate Fellowship # FL120100094 and grant # DP160101728 from the Australian Research Council and JSPS KAKENHI Grant Number JP17J00092. Nakazono gratefully acknowledges support from Erwin Schrödinger International Institute for Mathematics and Physics (ESI) that enabled his attendance at the conference on Elliptic Hypergeometric Functions in Combinatorics, Integrable Systems and Physics where this work was presented in March 20-24, 2017.
Appendix A Jacobian elliptic functions
To make this paper self-contained, we recall needed information about Jacobian elliptic functions and define notations (see [5, Chapter 22] and [23]).
The Jacobian elliptic functions are given by
[TABLE]
where the functions and are even, while the functions and are odd. The Jacobi theta-functions are given by
[TABLE]
where the functions , and are even, while the function is odd (see [23, §21.62] and [1, §16.31] for definitions of these standard functions). Here, and are the modulus and the complementary modulus of the elliptic sine, respectively, which satisfy
[TABLE]
and are complete elliptic integrals and also periodical parameters of Jacobian elliptic functions, e.g. the periods of the function are and . We note that the following formulae hold:
[TABLE]
Throughout this paper, we also use the functions , and defined by
[TABLE]
where
[TABLE]
Appendix B Relation between the - and -surface
In this section, we give the one-to-one correspondence between the - and -surface, which are characterized by the base points (B.7) and (2.4), respectively.
Using Landen transformations:
[TABLE]
and defining Weierstrass function in terms of the Jacobian elliptic function by
[TABLE]
where, , , satisfy
[TABLE]
we obtain the following one-to-one correspondence between Jacobian elliptic function and Weierstrass function:
[TABLE]
where
[TABLE]
By defining new coordinates and parameters , , and as the following:
[TABLE]
the eight base points (2.4) can be rewritten as
[TABLE]
which is the Weierstrass function’s setting (-surface) investigated in [17]. Therefore, the relations (B.6) give the one-to-one correspondence between the - and -surface.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] M. Abramowitz and I. Stegun. Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables . Dover Books on Mathematics. Dover Publications, 2012.
- 2[2] J. Atkinson, P. Howes, N. Joshi, and N. Nakazono. Geometry of an elliptic difference equation related to Q 4. J. Lond. Math. Soc. (2) , 93(3):763–784, 2016.
- 3[3] A. S. Carstea, A. Dzhamay, and T. Takenawa. Fiber-dependent deautonomization of integrable 2D mappings and discrete Painlevé equations. ar Xiv:1702.04907 .
- 4[4] J. Conway, N. Sloane, E. Bannai, R. Borcherds, J. Leech, S. Norton, A. Odlyzko, R. Parker, L. Queen, and B. Venkov. Sphere Packings, Lattices and Groups . Grundlehren der mathematischen Wissenschaften. Springer New York, 2013.
- 5[5] NIST Digital Library of Mathematical Functions . http://dlmf.nist.gov/, Release 1.0.14 of 2016-12-21. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller and B. V. Saunders, eds.
- 6[6] I. Dolgachev and D. Ortland. Point sets in projective spaces and theta functions. Astérisque , (165):210 pp. (1989), 1988.
- 7[7] I. V. Dolgachev. Classical algebraic geometry . Cambridge University Press, Cambridge, 2012. A modern view.
- 8[8] B. Grammaticos and A. Ramani. Discrete Painlevé equations: a review. In Discrete integrable systems , volume 644 of Lecture Notes in Phys. , pages 245–321. Springer, Berlin, 2004.
