Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
You have accessResearch article

Geometric description of a discrete power function associated with the sixth Painlevé equation

Published:https://doi.org/10.1098/rspa.2017.0312

    Abstract

    In this paper, we consider the discrete power function associated with the sixth Painlevé equation. This function is a special solution of the so-called cross-ratio equation with a similarity constraint. We show in this paper that this system is embedded in a cubic lattice with W~(3A1(1)) symmetry. By constructing the action of W~(3A1(1)) as a subgroup of W~(D4(1)), i.e. the symmetry group of PVI, we show how to relate W~(D4(1)) to the symmetry group of the lattice. Moreover, by using translations in W~(3A1(1)), we explain the odd–even structure appearing in previously known explicit formulae in terms of the τ function.

    1. Introduction

    The cross-ratio equation is a key discrete integrable equation, which connects discrete integrable systems and discrete differential geometry. One of its powerful roles is to provide a discrete analogue of the Cauchy–Riemann relation, thereby leading to a discrete theory of complex analytic functions. In this paper, we focus on the following system of partial difference equations:

    (zn,mzn+1,m)(zn+1,m+1zn,m+1)(zn+1,mzn+1,m+1)(zn,m+1zn,m)=1x1.1
    and
    (a0+a2+a4)zn,m=(na2)(zn+1,mzn,m)(zn,mzn1,m)zn+1,mzn1,m+(ma1a2a4)(zn,m+1zn,m)(zn,mzn,m1)zn,m+1zn,m1,1.2
    where (n,m)Z2 are independent variables, z is a dependent variable and x, ai, i=0,1,2,4, are parameters. Equation (1.1) is a special case of the cross-ratio equation and equation (1.2) is its similarity constraint [1]. The purpose of this paper is to provide a full geometric description of the system of equations (1.1) and (1.2) and its expression in terms of τ-variables.

    Nijhoff et al. [1] showed that equations (1.1) and (1.2) can be regarded as a part of the Bäcklund transformations of the Painlevé VI equation (PVI). In this context, x is identified as the independent variable of PVI, and ai, i=0,…,4 (a3 does not appear in the equations) are parameters corresponding to the simple roots of the affine root system of type D4(1) that relates to the symmetry group of PVI [25]. However, the complete characterization of equations (1.1) and (1.2) in relation to the Bäcklund transformations of PVI remains unclear.

    The solution of equations (1.1) and (1.2) for (n,m)Z+2 with the initial condition

    z0,0=0,z1,0=1andz0,1=erπi,1.3
    where i=1 and r=a0+a2+a4, gives rise to a discrete power function [68] when the parameters take special values:
    x=1,a2=0,a1+a2+a4=0.1.4
    (There is more than one definition of what constitutes a discrete power function. See Boole [9] for classical definitions.) Note that each of the ratios on the right-hand side of (1.2) is proportional to a harmonic mean of the forward and backward differences of zn,m, so that the continuum limit of equation (1.2) is
    ZdFdZ=XFX+YFY=2rF,1.5
    which is an equation satisfied by the power function F(Z)=Z2r, where Z=X+iYC. Equations (1.1) and (1.2) also appear as conditions for consistency of quadrilaterals on surfaces and enable a definition of discrete conformality [8,10]. It was shown in [11] that the special cases zn,0, zn,1 are given, respectively, by the Gamma function and Gauss hypergeometric functions. By applying the Bäcklund transformations, Ando et al. [12] and Hay et al. [13] found explicit formulae for the discrete power function in terms of the hypergeometric τ functions of PVI as stated in theorem 1.1.

    Theorem 1.1 [12,13].

    Define the functionτν(a,b,c;t)(cZ,νZ+)by

    τν(a,b,c;x)={det(φ(a+i1,b+j1,c;x))1i,jν(ν>0),1(ν=0),1.6
    where
    φ(a,b,c;x)=C0Γ(a)Γ(b)Γ(c)F(a,b,c;x)+C1Γ(ac+1)Γ(bc+1)Γ(2c)t1cF(ac+1,bc+1,2c;x).1.7
    Here, F(a,b,c;x) is the Gauss hypergeometric function, Γ(x) is the Gamma function, and C0and C1are arbitrary constants. Then, for(n,m)Z+2, zn,msatisfying equations (1.1) and (1.2) with
    a2=0,a1+a2+a4=01.8
    and the initial condition
    z0,0=0,z1,0=C0andz0,1=C1xr,1.9
    where r=a0+a2+a4is given as follows.
    • (1) Case where n≤m.

      zn,m={C1xrnN(r+1)N1(r+1)Nτn(N,rN+1,r;x)τn(N+1,rN+2,r+2;x),n+m is even,C1xrn(r+1)M1(r+1)M1τn(M+1,rM+1,r;x)τn(M+2,rM+2,r+2;x),n+m is odd.1.10

    • (2) Case where n≥m.

      zn,m={C0N(r+1)N1(r+1)Nτm(N+2,rN+1,r+2;x)τm(N+1,rN+2,r+2;x),n+m is even,C0(r+1)M1(r+1)M1τm(M+2,rM+1,r+1;x)τm(M+1,rM+2,r+1;x),n+m is odd.1.11

    Here, N=(n+m)/2, M=(n+m+1)/2 and (u)j=u(u+1)⋯(u+j−1) is the Pochhammer symbol.

    We remark that, in the case of generic values of parameters ai, the explicit expressions of zn,m in terms of τ-variables are given in [13]. This result was obtained by a heuristic approach and verified by the bilinear formulation of τ-variables. However, the appearance of the odd–even structure in the explicit formulae has not yet been fully understood. In §§4 and 5, we will see that this odd–even structure can be explained by the notion of projective reduction and the Weyl group symmetry on τ-variables.

    In this paper, we characterize equations (1.1) and (1.2) in relation to the Bäcklund transformations of PVI completely, and clarify the odd–even structure appearing in equations (1.10) and (1.11). The starting point is the appearance of the cross-ratio equation in a cubic lattice, as described in the ABS theory [1418]. We show how the cubic lattice is embedded in the lattice on which τ-variables of PVI are assigned. Four copies of the weight lattice of D4 will be needed to describe the lattice of τ-variables (see §3).

    In this paper, we refer to several lattices associated with affine reflection groups. A detailed description of (and notations for) these lattices is provided in the following section.

    (a) Background and notation

    For completeness, we explain the notation used to describe reflection groups, Coxeter groups and associated lattices in this paper. The notation is illustrated here for reflection groups arising from symmetries of cubes, but it applies to reflection groups of any dimension and type.

    A cube is left unchanged by certain reflections across hyperplanes. All such reflections can be expressed in terms of basic operations on vectors, which form a group, denoted by B3 [19]. In higher dimensions, we would have the reflection group Bn, n≥1, where the subscript refers to an integer number of dimensions. The group is generated by simple reflections, which we denote sj, j=1,…,n. The simple roots associated with these reflections are denoted αj, j=1,…,n, while corresponding co-roots are denoted αj.

    Repeated translation of a fundamental cube leads to a space-filling cubic lattice, on which equations (1.1) and (1.2) are iterated. Reflections across hyperplanes in the cubic lattice form an affine reflection group, which is denoted by Bn(1). We will refer also to the root lattice denoted by Q(Bn(1))=j=0nZαj. When interpreted as a Coxeter group, the reflection group is denoted by W(Bn(1)).

    In a Dynkin diagram, the simple reflections sj are represented by nodes. Obviously, such a diagram remains unchanged when the nodes are permuted. Such an operation is called a diagram automorphism. Let Ω be the group of such automorphisms. Incorporating these automorphisms into the space of allowable operations, we get a lift of the Coxeter group W(Bn(1)) to an extended Coxeter group denoted W~(Bn(1))=W(Bn(1))Ω [20].

    The classification of finite reflection groups leads to several fundamental types. We will encounter An, Bn and Dn in this paper, along with their corresponding lattices, affine Coxeter groups and extended Coxeter groups. Where direct sums of the same type of groups occur, we use an integer coefficient as an abbreviation, e.g. A1A1A1A1=:4A1. For conciseness, we replace ⊕ by +, e.g. D4A1A1=D4+2A1.

    (b) Main result and plan of the paper

    The main contribution of this paper, given in theorem 4.1, is to clarify the relation of the discrete power function to Bäcklund transformations of the PVI by including the system of defining difference equations into a lattice with W~(3A1(1)) symmetry. In general, a cubic lattice has W(B3(1)) as a symmetry group, but we chose to consider the symmetry group as W~(3A1(1)) by taking different selections of reflection hyperplanes. This idea allows us to realize W~(3A1(1)) as a subgroup of W~(D4(1)), which is the symmetry group of PVI. To achieve this realization, we use the combination of the cross-ratio equation together with its similarity relation, which is associated with the symmetry of the cubic lattice.

    The plan of this paper is as follows. In §2, we introduce a multi-dimensionally consistent cubic lattice with the symmetry W~(3A1(1)) on which the cross-ratio equation is naturally defined. In §3, we give actions of W~(D4(1)) on the parameters and the τ-variables. Based on this material, we construct a realization of W~(3A1(1)) in W~(D4(1)) and derive equations (1.1) and (1.2) in §4. We then explain in §5 the odd–even structures appearing in the explicit formula in [12,13]. Concluding remarks will be given in §6.

    2. Cubic lattice

    In this section, we describe a multi-dimensionally consistent cubic lattice in terms of a symmetry structure, which will turn out to be related to the symmetry group of PVI.

    While each cube inside the lattice is composed of six component faces, we will regard these as composed of two triple faces, i.e. a pair of three faces around a common vertex. We will see in §4 how this perspective is embedded in the symmetry group lattice of PVI. The symmetry group W~(3A1(1)) is suggested by the similarity equation

    ζ0zl1,l2=(l1β1)(zl1+1,l2zl1,l2)(zl1,l2zl11,l2)zl1+1,l2zl11,l2+(l2γ1)(zl1,l2+1zl1,l2)(zl1,l2zl1,l21)zl1,l2+1zl1,l21.2.1
    Here,
    β1=a2,γ1=a1+a2+a4andζ0=a0+a2+a4.2.2
    Indeed, the similarity equation has three parameters, β1, γ1 and ζ0, and two directions, ρ1 and ρ2, defined by zl1,l2=ρ1l1ρ2l2(z0,0), which shift the parameters as follows:
    ρ1:(β1,γ1,ζ0)(β11,γ1,ζ0)andρ2:(β1,γ1,ζ0)(β1,γ11,ζ0).2.3
    It is, therefore, natural to expect that there exists a Bäcklund transformation of the similarity equation, which shifts the parameters as follows:
    ρ0:(β1,γ1,ζ0)(β1,γ1,ζ0+1).2.4

    This implies that the three directions ρ0,ρ1,ρ2 are translations in W~(3A1(1)) and {β1,γ1,ζ0} are the parameters associated with the simple roots (or co-roots) of type 3A1(1). Consequently, we take the symmetry group of the cubic lattice to be W~(3A1(1)).

    Consider the cubic lattice constructed by the directions ρi, i=0,1,2. We place the following equations on each respective face of a triple of faces associated with each cube (referred to as the face equations):

    (ul1,l2,l0+ul1+1,l2,l0)(ul1+1,l2+1,l0+ul1,l2+1,l0)(ul1+1,l2,l0+ul1+1,l2+1,l0)(ul1,l2+1,l0+ul1,l2,l0)=κl1(1)κl2(2),2.5a
    (1ul1+1,l2,l0+1ul1,l2,l0)(ul1+1,l2,l0+1+ul1,l2,l0+1)=κl1(1)κl0(3)2.5b
    and(1ul1,l2+1,l0+1ul1,l2,l0)(ul1,l2+1,l0+1+ul1,l2,l0+1)=κl2(2)κl0(3),2.5c
    where {,κ1(i),κ0(i),κ1(i),}i=1,2,3 are parameters and
    ul1,l2,l0=ρ1l1ρ2l2ρ0l0(u0,0,0).2.6
    It is well known that this system of equations is multi-dimensionally consistent in the cubic lattice. Indeed, equation (2.5a) and equations (2.5b) and (2.5c) are called Q1 and H1, respectively, in the Adler–Bobenko–Suris (ABS) list (figure 1) [1418].
    Figure 1.

    Figure 1. The multi-dimensionally consistent cubic lattice. The face equations of bottom and top are given by Q1 and those of sides are given by H1.

    Using the geometric theory of PVI, we can realize the action of W~(3A1(1)) associated with these equations (see §4).

    3. Actions of affine Weyl group W~(D4(1))

    In this section, we provide the basic ingredients needed to describe actions of the symmetry group of PVI, i.e. W~(D4(1)). The explicit formulation of the τ function of PVI in terms of the weight lattice is provided here; to our knowledge, such an explicit formulation has not appeared in the literature. More precisely, four copies of the weight lattice of D4 will be needed to define the lattice of τ-variables. We will consider a root system of type D4+2A1 instead of D4 to give a full description.

    (a) Linear action

    In this subsection, we define the transformation group W~(D4(1)) and describe its linear actions on the weight lattice.

    We consider the following Z-modules:

    Q=Zα0+Zα1+Zα2+Zα3+Zα4+Zα5+Zα63.1
    and
    P=Zh0+Zh1+Zh2+Zh3+Zh4+Zh5+Zh6,3.2
    with the bilinear form ,:Q×PZ defined by
    αi,hj=δij,0i,j6.3.3
    The weight lattice of type D4(1) is spanned by hi, i=1,3,4 and h2h0. The generators {α0,,α4} and {α5,α6} are identified with the co-root of type D4(1) and type 2A1, respectively. Let us define the roots of type D4(1)+2A1: {α0,…,α6}, which satisfy
    αi,αj={Aij,0i,j42,(i,j)=(5,5),(6,6)0,otherwise,3.4
    by the following:
    (α0α1α2α3α4)=(Aij)i,j=04(h0h1h2h3h4),α5=2h5,α6=2h6.3.5
    Here, (Aij)i,j=04 is the generalized Cartan matrix of type D4(1),
    (Aij)i,j=04=(2010002100112110012000102).3.6
    We note that
    α0+α1+2α2+α3+α4=0.3.7

    We define the transformations si, i=0,…,4, which are reflections for the roots {α0,…,α4}, by

    si(λ)=λαi,λαi,i=0,,4,λP,3.8
    which give
    si:hih2hi,i=0,1,3,4,s2:h2h0+h1h2+h3+h4.3.9

    Note that the hi’s which are not explicitly shown in equation (3.9) remain unchanged. Also, define the transformations σi, i=1,2,3, which are the automorphisms of the Dynkin diagram of type D4(1) (figure 2) and the reflections for the simple roots α5 and α6, by

    σ1:(h0,h1,h2,h3,h4h5,h6)(h1+h5,h0+h5,h2+2h5,h4+h5,h3+h5h5,h6),3.10a
    σ2:(h0,h1,h2,h3,h4h5,h6)(h3+h6,h4+h6,h2+2h6,h0+h6,h1+h6h5,h6)3.10b
    andσ3=σ1σ2.3.10c
    From definitions (3.3), (3.5), (3.8) and (3.10), we can compute actions on the simple roots αi, i=0,…,6,
    si:(αi,α2)(αi,α2+αi),i=0,1,3,4,3.11a
    s2:(α0,α1,α2,α3,α4)(α0+α2,α1+α2,α2,α3+α2,α4+α2),3.11b
    σ1:(α0,α1,α3,α4,α5)(α1,α0,α4,α3,α5),3.11c
    σ2:(α0,α1,α3,α4,α6)(α3,α4,α0,α1,α6)3.11d
    andσ3:(α0,α1,α3,α4,α5,α6)(α4,α3,α1,α0,α5,α6).3.11e
    Under the linear actions on the weight lattice (3.8) and (3.10), W~(D4(1))=s0,,s4,σ1,σ2 forms an extended affine Weyl group of type D4(1). Indeed, the following fundamental relations hold:
    (sisj)mij=1,0i,j4,σi2=1,i=1,23.12a
    and
    σ1s{0,1,2,3,4}=s{1,0,2,4,3}σ1,σ2s{0,1,2,3,4}=s{3,4,2,0,1}σ2,σ1σ2=σ2σ1,3.12b
    where
    mij={1,i=j3,i=2,j2ori2,j=22,otherwise.3.13
    Figure 2.

    Figure 2. Dynkin diagram of type D4(1).

    Remark 3.1

    We can also define the action of W~(D4(1)) on the co-root lattice Q by replacing αi with αi in the actions (3.11). Then, the transformations in W~(D4(1)) preserve the form 〈⋅,⋅〉, that is, the following holds:

    γ,λ=w(γ),w(λ),3.14
    for arbitrary wW~(D4(1)), γQ and λ∈P.

    (b) Action on τ-variables

    It is well known that we can extend the linear action of W~(D4(1)) to the action of the τ function of the Painlevé VI equation [2,3]. In this section, we define the τ function on the weight lattice and give the action of W~(D4(1)) on it.

    Let Mk, k=0,…,4, be the orbits of hi, i=0,…,4, defined by

    Mk={w(hk)|wW~(D4(1))},k=0,,43.15
    and M be their disjoint union defined by
    M=k=04Mk.3.16
    Each of Mk is decomposed to the orbit of W(D4(1)) by
    Mk=Mk(00)Mk(10)Mk(01)Mk(11),3.17
    where
    Mk(ij)={λMk|α5,λ=i,α6,λ=j},k=0,1,3,43.18
    and
    M2(ij)={λM2|α5,λ=2i,α6,λ=2j}.3.19
    For fixed k∈{0,1,2,3,4}, Mk(ij) are transformed to each other by the action of Dynkin automorphisms σi, i=1,2,3. For example, the weight h0 lies in M0(00) and the orbit of h0 under σ1 and σ2 can be calculated by equations (3.3) and (3.10), which give
    σ1(h0)=h1+h5M0(10)andσ2(h0)=h3+h6M0(01).3.20

    Remark 3.2

    • (i) Given any wW~(D4(1)), the element w(hk), k=0,1,3,4, can be expressed as i=04nihi+lh5+mh6, where niZ, l,m∈{0,1}. By computing bilinear forms with α5 and α6, we see that w(hk) lies in Mk(lm).

    • (ii) From the relations (3.12b), we can express any element of W~(D4(1)) as one of

      w1,σ1w2,σ2w3,σ1σ2w4,
      where wiW(D4(1)), i=1,…,4. Note that the action of wi, i=1,…,4, on i=04nihi+lh5+mh6 does not change l and m but σi, i=1,2, does change them.

    • (iii) Based on these observations in (i) and (ii), it follows that if w(hk)Mk(lm), where k=0,1,3,4 and wW~(D4(1)), the superscript ‘(lm)’ of Mk(lm) is determined by the number of iterations of σ1 and σ2 occurring in w as

      l={0(if number of σ1 is even),1(if number of σ1 is odd)andm={0(if number of σ2 is even),1(if number of σ2 is odd).3.21
      This is easily verified from the actions of σi, i=1,2, on hk given in (3.10).

    • (iv) It follows that for any element w(hk)=i=04nihi+lh5+mh6, where k=0,1,3,4 and wW~(D4(1)), l and m are uniquely determined from ni, i=0,…,4. For example, suppose that wW~(D4(1)) is decomposed as σ1w1, where w1W(D4(1)). Then we have from (3.10) that

      w1(h0)=n1h0+n0h1+n2h2+n4h3+n3h4+(n0+n1+2n2+n3+n4l)h5+mh6.3.22
      Since w1(h0) should not include h5 and h6 as seen in equation (3.9), we have l=n0+n1+2n2+n3+n4 and m=0.

    Lemma 3.3.

    The following properties hold.

    • (i) Lattices Mkare mutually disjoint.

    • (ii) Each of Mk, k=0,1,3,4, is isomorphic to the weight lattice of typeD4(1).

    Proof.

    First we prove the property (i). Define

    δ=α0+α1+2α2+α3+α4,3.23
    which is invariant under the action of W~(D4(1)). From remark 3.1, it is obvious that M2 is disjoint from the others, since 〈δ,λ〉=2 for any λ∈M2 and 〈δ,λ〉=1 for any λ∈Mk, k=0,1,3,4. Let us consider the 16 lattices Mk(ij), k=0,1,3,4. It is easy to see that these are mutually disjoint if the superscripts (ij) are different from each other. Then it is sufficient to show that Mk(00), k=0,1,3,4, are mutually disjoint. We prove it by contradiction. Suppose that there exists an element of M0(00)M1(00). Since we are now restricting our attention to Mk(00), we use only the actions of W(D4(1)). Then, there exists an element wW(D4(1)) such that h1=w(h0), which is the reflection with respect to the root vector h1h0. Therefore, it contradicts the assumption.

    Next, we show the property (ii). We prove it for M0, by constructing a map defined by

    p:M0k=04Zhk,3.24
    where every element of M0 is given by i=04nihi+lh5+mh6, with the proviso that l and m are determined by {ni}, as mentioned in remark 3.2. We define the image of such an element under p as i=04nihi. It follows from this construction that p is injective, because two elements in the domain of p with the same image must therefore be identical. Since p is injective and 〈δ,p(λ)〉=1 for λ∈M0, lattice M0 is isomorphic to the weight lattice of type D4. Note that we have
    p(M0)=p(M0(00))p(M0(10))p(M0(01))p(M0(11))M0(00)M1(00)M3(00)M4(00).3.25
    In a similar manner, we can prove the property (ii) for Mk, k=1,3,4. Therefore, we have completed the proof. ▪

    We now consider the τ-variables assigned on the lattice M, and the action of W~(D4(1)) on them. As we will show in §4, the dependent variable of the system of partial difference equations (1.1) and (1.2) is defined by a ratio of τ-variables, and the system is derived from the action of the Weyl group on the τ-variables.

    First, let {τi}i=04, τ2(σ1), τ2(σ2), τ2(σ3), τ2(s2) be the variables and ai, i=0,…,4, x be the complex parameters satisfying the following conditions:

    xτ2+x1/2(x1)1/2τ2(σ1)=ix1/2τ2(σ2)=τ2i(x1)1/2τ2(σ3)3.26
    and
    a0+a1+2a2+a3+a4=1.3.27
    Therefore, the number of essential variables and parameters are 7 and 5, respectively. We now proceed to give the actions of W~(D4(1)) on these τ-variables, which were deduced in [2,3],
    s0:τ0u2v2τ2(σ1)τ0,τ2(s2)u2v2τ2(σ1)τ2(s2)a0u2τ0134τ02,s1:τ1τ2τ1,τ2(s2)τ2(s2)τ2τ12,s2:{τ2τ2(s2),τ2(s2)τ2,τ2(σ1)τ2(σ1)τ2(s2)+a2u4v2τ0134τ2,τ2(σ2)τ2(σ2)τ2(s2)+ia2u4τ0134τ2,τ2(σ3)τ2(σ3)τ2(s2)+ia2u2v2τ0134τ2,s3:τ3iv2τ2(σ3)τ3,τ2(s2)iv2τ2(σ3)τ2(s2)a3u2τ0134τ32,s4:τ4iu2τ2(σ2)τ4,τ2(s2)iu2τ2(σ2)τ2(s2)a4u2τ0134τ42,σ1:{τ0uvτ1,τ1u1v1τ0,τ3u1vτ4,τ4uv1τ3,τ2τ2(σ1),τ2(σ1)τ2,τ2(σ2)τ2(σ3),τ2(σ3)τ2(σ2),τ2(s2)s2(τ2(σ1)),σ2:{τ0eπi/4uτ3,τ1eπi/4u1τ4,τ3eπi/4u1τ0,τ4eπi/4uτ1,τ2τ2(σ2),τ2(σ1)τ2(σ3),τ2(σ2)τ2,τ2(σ3)τ2(σ1),τ2(s2)s2(τ2(σ2)),3.28
    where
    u=x1/4,v=(x1)1/4,τ0134=τ0τ1τ3τ4.3.29
    Here, we define the action on the parameters ai, i=0,…,4, by
    si(aj)=ajAijai,0i,j43.30a
    and
    σ1a{0,1,2,3,4}=a{1,0,2,4,3}andσ2a{0,1,2,3,4}=a{3,4,2,0,1},3.30b
    where Aij is given in (3.6). Note that each parameter ai is associated with the root αi, i=0,…,4, and the condition (3.27) corresponds to (3.7).

    The fundamental relations (3.12) also hold under the actions si, i=0,…,4, on the τ-variables and the parameters. Note that the extended parts σi, i=1,2, are modified to become

    σ22.(τ0,τ1,τ2,τ3,τ4τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2))=(iτ0,iτ1,τ2,iτ3,iτ4τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2)),3.31a
    σ1s2.(τ2,τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2))=s2σ1.(τ2,τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2))3.31b
    andσ1σ2.(τ0,τ1,τ2,τ3,τ4τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2))=σ2σ1.(iτ0,iτ1,τ2,iτ3,iτ4τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2)).3.31c

    Remark 3.4

    The correspondence between the notations in this paper and those in [2,3] is given by

    (a0,a1,a2,a3,a4,xs0,s1,s2,s3,s4,σ1,σ2,σ3τ0,τ1,τ2,τ3,τ4τ2(σ1),τ2(σ2),τ2(σ3),τ2(s2))(α0,α1,α2,α3,α4,ts0,s1,s2,s3,s4,s5,s6,s7τ0,τ1,τ2,τ3,τ4s5(τ2),s6(τ2),s7(τ2),s2(τ2)).3.32

    We then define a mapping τ on M by

    τ(hi)=τi,τ(σj.h2)=τ2(σj),τ(s2.h2)=τ2(s2),τ(w.hi)=w(τi),3.33
    where i=0,…,4, j=1,2,3 and wW~(D4(1)).

    4. Geometric description of discrete power function

    In this section, we provide our main result, which is based on the action of W~(D4(1)) constructed in the previous sections.

    Theorem 4.1.

    • I. There exist elementsρiW~(D4(1)), i=0,1,2, which satisfy the following properties.

      • (i) (Projective reduction) For each i=0,1,2, the element ρi2is a translation ofW~(D4(1)), while ρiitself is not a translation.

      • (ii) (Commutativity) ρiρjjρi, for i,j=0,1,2.

      • (iii) (Bilinear forms) There exists a subspaceZβ1Zγ1Zζ1, of the parameter spaceZa1Za2Za3Za4, on which the actions become translations. That is, the projections of ρi, i=0,1,2, on the subspace satisfy the following properties:

        (β1|γ1)=(γ1|ζ1)=(ζ1|β1)=0and(β1|β1)=(γ1|γ1)=(ζ1|ζ1)=2,4.1
        where the bilinear form ( | ) on the parameter space is defined by
        (ai|aj)=αi,αj=Aij.4.2

      • (iv) (Shift actions) Each ρi, i=0,1,2, shifts one of {β111} but not the others, i.e.

        ρ1(β1)β1Z0,ρ1:(γ1,ζ1)(γ1,ζ1),4.3a
        ρ2(γ1)γ1Z0,ρ2:(β1,ζ1)(β1,ζ1)4.3b
        andρ0(ζ1)ζ1Z0,ρ0:(β1,γ1)(β1,γ1).4.3c

    • II. Define zl1,l2by

      zl1,l2=(1)l1+l2ρ1l1ρ2l2(z0),4.4
      where
      z0=ρ02(τ0)τ0.4.5
      Then zl1,l2satisfies the cross-ratio equation
      (zl1,l2zl1+1,l2)(zl1+1,l2+1zl1,l2+1)(zl1+1,l2zl1+1,l2+1)(zl1,l2+1zl1,l2)=1x4.6
      and similarity equation (2.1).

    Remark 4.2

    • (a) We note that the Z-module spanned by the parameters ai, i=0,…,4, gives a co-root- or root-lattice of type D4(1) denoted by Q(D4(1)), with the pairing ( | ).

    • (b) The properties (iii) and (iv) are natural requirements since (iii) is the property that β1, γ1, ζ1 form simple roots of type 3A1(1) and (iv) comes from the idea that ρi are the fundamental translations in W~(3A1(1)).

    • (c) The choice of z0 was motivated by the observation that the solution of the cross-ratio equation is given by a ratio of two τ functions, where one of the τ functions lies in a direction orthogonal to the plane in which the equation is iterated [13]. This observation also shows that the variable must be defined in terms of a square of ρ0. The remaining directions ρ1 and ρ2 act as the two shift directions.

    • (d) By choosing τ-variables appropriately as in [12], we can reproduce the initial values (1.3) of the discrete power function solution of (1.1) and (1.2). In particular, note that the choice τ0=1 (up to gauge) and ρ02(τ0)=0 leads to the value z0=0.

    • (e) The property (i) is called a projective reduction, in the sense that we took ρi, which were not translations in W~(D4(1)), and considered them in a subgroup in which they became translations. Many discrete integrable systems can be derived by using such a procedure [2124].

    Proof.

    We begin with the construction of the subspace of parameters needed in item (iii). Choose

    β1=a2.4.7
    The parameters γ1 and ζ1 can be chosen from
    a1+a2+a4,a1+a2+a3anda2+a3+a44.8
    to satisfy the conditions in item (iii). For example, take γ1=a1+a2+a4. Then (β1|γ1)=0 and the other required conditions follow from equations (3.6) and (4.2). There is one remaining choice, which we label as μ1. To be specific, from now on, we define
    γ1=a1+a2+a4,ζ1=a1+a2+a3andμ1=a2+a3+a44.9
    and also define an additional set of four parameters
    β0=1β1,γ0=1γ1,ζ0=1ζ1andμ0=1μ1,4.10
    which are needed as simple roots of 4A1(1).

    The Z-module spanned by the parameters βi, γi, ζi, μi, i=0,1, denoted by Q(4A1(1)), gives a co-root- or root-lattice of type 4A1(1), with the pairing ( | ). The transformation group for Q(4A1(1)) forms an extended affine Weyl group of type 4A1(1). (See appendix A for more detail.) Notice that by relinquishing one pair of parameters, say μi, i=0,1, we obtain Q(3A1(1)) from Q(4A1(1)). We here focus on the sublattice spanned by the parameters βi, γi, ζi, denoted by Q(3A1(1)). Nevertheless, we continue to use the notations including μ such that

    πβμ:(β0,β1,μ0,μ1)(β1,β0,μ1,μ0),4.11a
    πγμ:(γ0,γ1,μ0,μ1)(γ1,γ0,μ1,μ0)4.11b
    andπζμ:(ζ0,ζ1,μ0,μ1)(ζ1,ζ0,μ1,μ0),4.11c
    for consistency with appendix A.

    We next construct ρi in W~(3A1(1)) by using the elements of W~(D4(1)). Recall that W~(D4(1))=s0,,s4,σ1,σ2,σ3. The extended affine Weyl group W~(3A1(1)) corresponding to the parameters chosen above is given by

    W~(3A1(1))=sβ0,sβ1,sγ0,sγ1,sζ0,sζ1,πβμ,πγμ,πζμ,4.12
    where
    sβ0=s4s3s1s0s2s4s3s1s0,sβ1=s2,sγ0=s0s3s2s0s3,sγ1=s1s4s2s1s4,4.13a
    sζ0=σ2s1s3s2s1s3σ2,sζ1=s1s3s2s1s3,πβμ=σ1s4s3s1s0,πγμ=σ34.13b
    andπζμ=σ2.4.13c
    The action of W~(3A1(1)) on the parameters βi, γi, ζi, μi, i=0,1, is given by
    sβ0:(β0,β1)(β0,β1+2β0),sβ1:(β0,β1)(β0+2β1,β1),4.14a
    sγ0:(γ0,γ1)(γ0,γ1+2γ0),sγ1:(γ0,γ1)(γ0+2γ1,γ1),4.14b
    sζ0:(ζ0,ζ1)(ζ0,ζ1+2ζ0),sζ1:(ζ0,ζ1)(ζ0+2ζ1,ζ1)4.14c
    and (4.11). We are now in a position to define ρi by using the translations in W~(3A1(1)),
    ρ1=sβ0πβμ,ρ2=sγ0πγμandρ0=sζ0πζμ,4.15
    whose actions on the parameters βi, γi, ζi, μi, i=0,1, are given by
    ρ1:(β0,β1,μ0,μ1)(β0+1,β11,μ1,μ0),4.16a
    ρ2:(γ0,γ1,μ0,μ1)(γ0+1,γ11,μ1,μ0)4.16b
    andρ0:(ζ0,ζ1,μ0,μ1)(ζ0+1,ζ11,μ1,μ0).4.16c
    We can easily verify that the transformations ρi, i=0,1,2, satisfy the properties (i), (ii) and (iv).

    For the second part of the theorem, note that, with z0 as defined in (4.5), the remaining directions ρ1 and ρ2 act as the two shift directions. Now we introduce

    z1=ρ1(z0),z2=ρ2(z0)andz12=ρ1ρ2(z0).4.17
    We can verify that these variables z0, z1, z2 and z12 are related by
    (z0+z1)(z12+z2)(z1+z12)(z2+z0)=1x,4.18
    by using the framework constructed, namely (3.28), (4.13), (4.15), (4.5) and (4.17). For example, to express z1 in terms of τ-variables, we need the composition of ρ1 with ρ0. The starting point is (4.15), which in turn are given by (4.13), and the actions on τ functions (3.28). In this way, we deduce
    z0+z1=ζ0qx1/2τ1τ3τ0τ4,z12+z2=ζ0xqx1/2(1q)τ1τ3τ0τ44.19a
    and
    z1+z12=ζ0(xq)qx1/2(1q)τ1τ3τ0τ4,z2+z0=ζ0x1/2τ1τ3τ0τ4,4.19b
    where q=ix1/2τ2(σ2)/τ2. Then we immediately obtain equation (4.18) from equations (4.19).

    Furthermore, from the actions of the inverses of ρ1 and ρ2, we find

    ζ0z0=β1(z1+z0)(z0+ρ11(z0))z1+ρ11(z0)γ1(z2+z0)(z0+ρ21(z0))z2+ρ21(z0).4.20
    Using equations (4.18) and (4.20), respectively, we then obtain the cross-ratio equation (4.6) and similarity equation (2.1) by defining the dependent variable by (4.4). ▪

    Remark 4.3

    In the proof of part II of the theorem, we used the action of W~(D4(1)) on τ functions. However, an alternative approach is available by starting with the definitions of the variables z0, z1, z2, z12, i.e. (4.5) and (4.17). For this purpose, action (3.28) is interpreted in terms of the variables z0, z1, z2, z12. The resulting expressions are very large and to simplify these we can instead introduce z12 by using (4.18), namely

    z12=xz2(z0+z1)z1(z2+z0)x(z0+z1)(z2+z0).4.21
    Complete details are provided in appendix B.

    Remark 4.4

    Let

    ul1,l2,l0=ρ1l1ρ2l2ρ0l0(z0),l1,l2,l0Z.4.22
    From equation (4.18) and the following actions of ρi, 1,2,3:
    (1z1+1z0)(ρ0(z1)+ρ0(z0))=ζ0(ζ0+1)x4.23a
    and
    (1z2+1z0)(ρ0(z2)+ρ0(z0))=ζ0(ζ0+1),4.23b
    we obtain
    (ul1,l2,l0+ul1+1,l2,l0)(ul1+1,l2+1,l0+ul1,l2+1,l0)(ul1+1,l2,l0+ul1+1,l2+1,l0)(ul1,l2+1,l0+ul1,l2,l0)=1x,4.24a
    (1ul1+1,l2,l0+1ul1,l2,l0)(ul1+1,l2,l0+1+ul1,l2,l0+1)=(ζ0+l0)(ζ0+l0+1)x4.24b
    and(1ul1,l2+1,l0+1ul1,l2,l0)(ul1,l2+1,l0+1+ul1,l2,l0+1)=(ζ0+l0)(ζ0+l0+1).4.24c
    These equations are equivalent to the lattice equations of ABS type discussed in §2. Indeed, equations (4.24) can be obtained from equations (2.5) by the following specialization of parameters:
    κl(1)=1x,κl(2)=1andκl(3)=(ζ0+l)(ζ0+l+1),4.25
    where lZ.

    5. Description by weight lattice

    In this section, we relate our results to those obtained from a different perspective, namely the determinantal structure of hypergeometric solutions. We show that our z-variable (4.4) is identical to that found in [12,13].

    To make this correspondence, we consider the actions of ρi, i=0,1,2, on the sublattice of the weight lattice. We here refer to a vector (wid).x for wW~(D4(1)), xM0, as a displacement vector corresponding to w. Since the dependent variable of the discrete power function in terms of the τ function is given by

    zl1,l2=(1)l1+l2ρ02ρ1l1ρ2l2(τ0)ρ1l1ρ2l2(τ0)5.1
    and ρi, i=0,1,2, satisfy the properties (i) and (ii) in §4, we here consider the displacement vectors of ρ1 and ρ2 on the following sublattice of M0 (3.15):
    L=L(0)L(1)L(2)L(12),5.2
    where
    L(0)={ρ12l1ρ22l2(h0)|l1,l2Z}=h0+Zv1+Zv2,5.3a
    L(1)={ρ12l1+1ρ22l2(h0)|l1,l2Z}=h1+h2+h5+Zv1+Zv2,5.3b
    L(2)={ρ12l1ρ22l2+1(h0)|l1,l2Z}=h4+h5+h6+Zv1+Zv25.3c
    andL(12)={ρ12l1+1ρ22l2+1(h0)|l1,l2Z}=h2h3+h6+Zv1+Zv2.5.3d
    Here, vi, i=0,1,2, are the displacement vectors on M0 given by
    v0=h0+h1+h3h4,v1=h0h1+2h2h3h4andv2=h0+h1h3+h4,5.4
    which correspond to ρi2, i=0,1,2, respectively. As illustrated in figure 3, the transformations ρ1 and ρ2 correspond to the vectors on L as the following:
    ρ1{h0h1+h2+h5,on L(0)L(12),h2h3h4h5,on L(1)L(2)5.5
    and
    ρ2{h0+h4+h5+h6,on L(0),h1h3h5+h6,on L(1),h1h3h5h6,on L(2),h0+h4+h5h6,on L(12).5.6
    Figure 3.

    Figure 3. The actions of ρ1 and ρ2 on the lattice L. The sublattices are distinguished by producing them in different colours, that is, L(0), L(1), L(2) and L(12) are coloured in black, red, blue and green, respectively. (Online version in colour.)

    Let us define the following translations in W~(D4(1)) by

    T^13=s1s2s0s4s2s1σ3,T^40=s4s2s1s3s2s4σ3,T^34=s3s2s0s1s2s3σ15.7a
    and
    T14=s1s4s2s0s3s2σ2.5.7b
    Using these translations, the actions of ρ02, ρ1 and ρ2 on L are expressed by
    ρ02=T^131T^401T^342,5.8a
    ρ1={T^401T14,on L(0)L(12),T^131T14,on L(1)L(2)5.8b
    andρ2={T^401,on L(0)L(12),T^131,on L(1)L(2).5.8c
    Therefore, the dependent variable of the discrete power function can be expressed by
    zl1,l2={(1)l11τ(l1+l2)/21,(l1+l2)/21,2,l1τ(l1+l2)/2,(l1+l2)/2,0,l1,if l1+l2 is even,(1)l11τ(l1+l2+1)/2,(l1+l2+1)/21,2,l1τ(l1+l21)/2,(l1+l2+1)/2,0,l1,if l1+l2 is odd,5.9
    where
    τk,l,m,n=T^13kT^40lT^34mT14n.5.10
    This expression is equivalent to that in [12,13]. Note that the difference of the coefficients between (5.1) and (5.9) arises from equations (3.31).

    6. Concluding remarks

    In this paper, we have investigated the geometric structure of the system of partial difference equations (1.1) and (1.2). The cross-ratio equation (1.1) is naturally defined on a multi-dimensionally consistent cubic lattice as a face equation from the ABS theory. On the other hand, we have shown in this paper that its similarity constraint (1.2) is derived from the symmetries of the lattice.

    However, it was known that the system (1.1) and (1.2) could be derived from the Bäcklund transformations of PVI, which form W~(D4(1)). The question of how to relate the symmetry group W~(D4(1)) to the symmetry group of the lattice remained open in the literature. We have answered this question in this paper by constructing the action of W~(3A1(1)) as a subgroup of W~(D4(1)). Moreover, considering the realization on the level of the τ function simultaneously, we gave the geometric characterizations of the discrete power function. In particular, we explained the odd–even structure appearing in the explicit formula in [12,13] in terms of the projective reduction.

    It is interesting to relate the geometric structures reported in this paper to other properties of (1.1) and (1.2) reported in the literature. For example, the discrete part of the Lax pair given in [25] can be obtained from translations on the W~(3A1(1)). However, the remaining part of the Lax pair involves a monodromy variable, which is not visible from the point of view of the Painlevé VI equation itself. We believe that this is related to the choice of gauge we took in defining the τ functions. This is an open subject for further investigation.

    Not many examples of discrete complex analytic functions are known. In [26], one was constructed on the hexagonal lattice. Further interesting directions concern the remaining Painlevé and discrete Painlevé equations. We expect our analysis to be useful for such equations, not only of second order, to give geometric characterization of known examples or to identify new examples. Our results in these directions will be reported in forthcoming publications.

    Data accessibility

    This paper does not contain any additional data.

    Authors' contributions

    All authors contributed equally to writing the paper.

    Competing interests

    We have no competing interests.

    Funding

    This research was supported by an Australian Laureate Fellowship no. FL120100094 and grant no. DP160101728 from the Australian Research Council and JSPS KAKENHI grant nos. JP16H03941, JP16K13763 and JP17J00092.

    Appendix A. Extended affine Weyl group W~(4A1(1))

    In this section, we consider an affine Weyl group for Q(4A1(1)) spanned by the parameters βi, γi, ζi, μi, i=0,1.

    Let

    sβ0=s4s3s1s0s2s4s3s1s0,sβ1=s2,sγ0=s0s3s2s0s3,sγ1=s1s4s2s1s4,A 1a
    sζ0=σ2s1s3s2s1s3σ2,sζ1=s1s3s2s1s3,sμ0=s0s1s2s0s1,A 1b
    sμ1=s3s4s2s3s4,πβγ=σ2s4s3s1s0,πβζ=σ3s4s3s1s0,πβμ=σ1s4s3s1s0A 1c
    andπγζ=σ1,πγμ=σ3,πζμ=σ2,r1=s1s4,r2=s1s3,r3=s3s4,A 1d
    whose actions on the parameters βi, γi, ζi, μi, i=0,1, are given by
    sβ0:(β0,β1)(β0,β1+2β0),sβ1:(β0,β1)(β0+2β1,β1),A 2a
    sγ0:(γ0,γ1)(γ0,γ1+2γ0),sγ1:(γ0,γ1)(γ0+2γ1,γ1),A 2b
    sζ0:(ζ0,ζ1)(ζ0,ζ1+2ζ0),sζ1:(ζ0,ζ1)(ζ0+2ζ1,ζ1),A 2c
    sμ0:(μ0,μ1)(μ0,μ1+2μ0),sμ1:(μ0,μ1)(μ0+2μ1,μ1),A 2d
    πβγ:(β0,β1,γ0,γ1)(β1,β0,γ1,γ0),πβζ:(β0,β1,ζ0,ζ1)(β1,β0,ζ1,ζ0),A 2e
    πβμ:(β0,β1,μ0,μ1)(β1,β0,μ1,μ0),πγζ:(γ0,γ1,ζ0,ζ1)(γ1,γ0,ζ1,ζ0),A 2f
    πγμ:(γ0,γ1,μ0,μ1)(γ1,γ0,μ1,μ0),πζμ:(ζ0,ζ1,μ0,μ1)(ζ1,ζ0,μ1,μ0),A 2g
    r1:(βi,γi,ζi,μi)(γi,βi,μi,ζi),r2:(βi,γi,ζi,μi)(ζi,μi,βi,γi)A 2h
    andr3:(βi,γi,ζi,μi)(μi,ζi,γi,βi),A 2i
    where i=0,1. They collectively form the extended affine Weyl group of type 4A1(1):
    W~(4A1(1))=sβ0,sβ1,sγ0,sγ1,sζ0,sζ1,sμ0,sμ1,πβγ,πβζ,πβμ,πγζ,πγμ,πζμ,r1,r2,r3.A 3
    Indeed, W~(4A1(1)) satisfies the following fundamental relations:
    sβ02=sβ12=sγ02=sγ12=sζ02=sζ12=sμ02=sμ12=1,A 4a
    (sβ0sβ1)=(sγ0sγ1)=(sζ0sζ1)=(sμ0sμ1)=1,A 4b
    sβisγj=sγisβj,sβisζj=sζisβj,sβisμj=sμisβj,sγisζj=sζisγj,A 4c
    sγisμj=sμisγj,sζisμj=sμisζj,A 4d
    πβγ2=πβζ2=πβμ2=πγζ2=πγμ2=πζμ2=1,A 4e
    πγζ=πβγπβζ,πγμ=πβγπβμ,πζμ=πβζπβμ,A 4f
    πβγπβζ=πβζπβγ,πβγπβμ=πβμπβγ,πβζπβμ=πβμπβζ,A 4g
    r12=r22=r32=1,r3=r1r2=r2r1,A 4h
    πβγs{β0,β1,γ0,γ1,ζi,μi}=s{β1,β0,γ1,γ0,ζi,μi}πβγ,A 4i
    πβζs{β0,β1,γi,ζ0,ζ1,μi}=s{β1,β0,γi,ζ1,ζ0,μi}πβζ,A 4j
    πβμs{β0,β1,γi,ζi,μ0,μ1}=s{β1,β0,γi,ζi,μ1,μ0}πβμ,A 4k
    r1s{βi,γi,ζi,μi}=s{γi,βi,μi,ζi}r1,r2s{βi,γi,ζi,μi}=s{ζi,μi,βi,γi}r2A 4l
    r1π{βγ,βζ,βμ,γζ,γμ,ζμ}=π{βγ,γμ,γζ,βμ,βζ,ζμ}r1andr2π{βγ,βζ,βμ,γζ,γμ,ζμ}=π{ζμ,βζ,γζ,βμ,γμ,βγ}r2.A 4m
    Here, the relation (ww)=1 for transformations w and w′ means that there is no positive integer N such that (ww′)N=1. Note that, from the definitions (4.13) and (A.1), W~(3A1(1))W~(4A1(1)) holds.

    Appendix B. Birational actions of W~(3A1(1)) on the z-variables

    The action of W~(3A1(1)) on the variables z0, z1, z2, z12, where z12 is related to other variables as (4.21), is given by the following:

    sβ0(z0)=2β0z1z12(β0+γ0+ζ0μ0)z0z1+(β0γ0+ζ0+μ0)z0z122β0z0+(β0γ0ζ0+μ0)z1(β0+γ0ζ0μ0)z12,B 1a
    sβ0(z2)=2β0z1z12+(β0+γ0+ζ0μ0)z1z2(β0γ0+ζ0+μ0)z2z122β0z2(β0γ0ζ0+μ0)z1+(β0+γ0ζ0μ0)z12,B 1b
    sγ0(z0)=2γ0z2z12(β0γ0ζ0μ0)z0z12(β0+γ0+ζ0μ0)z0z22γ0z0(β0+γ0ζ0μ0)z12(β0γ0+ζ0μ0)z2,B 1c
    sγ0(z1)=2γ0z2z12+(β0γ0ζ0μ0)z1z12+(β0+γ0+ζ0μ0)z1z22γ0z1+(β0+γ0ζ0μ0)z12+(β0γ0+ζ0μ0)z2,B 1d
    sζ0:(z0,z1,z2,z12)(z01,z11,z21,z121),B 1e
    sβ1=πβμsβ0πβμ,sγ1=πγμsγ0πγμ,sζ1=πζμsζ0πζμ,B 1f
    πβμ:(z0,z1,z2,z12)(z1,z0,z12,z2),πγμ:(z0,z1,z2,z12)(z2,z12,z0,z1),B 1g
    πζμ(z0)=ζ0(β0γ1+ζ0μ1)z0+(β0+γ1ζ0μ1)z12(z0+z2)(z1+z12)ζ0(β1γ1ζ0+μ0)z2+(β1+γ1ζ0+μ0)z122(z0+z2)(z1+z12),B 1h
    πζμ(z1)=ζ0(β1+γ1ζ0μ0)z0+(β1γ1+ζ0μ0)z12(z0+z2)(z1+z12)ζ0(β0+γ1ζ0+μ1)z2+(β0γ1ζ0+μ1)z122(z0+z2)(z1+z12),B 1i
    πζμ(z2)=ζ0(β1γ0ζ0+μ1)z0+(β1+γ0ζ0+μ1)z12(z0+z2)(z1+z12)ζ0(β0γ0+ζ0μ0)z2+(β0+γ0ζ0μ0)z122(z0+z2)(z1+z12)B 1j
    andπζμ(z12)=ζ0(β0+γ0ζ0+μ0)z0+(β0γ0ζ0+μ0)z12(z0+z2)(z1+z12)ζ0(β1+γ0ζ0μ1)z2+(β1γ0+ζ0μ1)z122(z0+z2)(z1+z12).B 1k
    Under the actions on the parameters (4.11) and (4.13) and z-variables (B.5), the following fundamental relations hold:
    sβi2=sγi2=sζi2=1,i=0,1,(sβ0sβ1)=(sγ0sγ1)=(sζ0sζ1)=1,B 2a
    (sβisγj)2=(sβisζj)2=(sγisζj)2=1,i,j=0,1,πβμ2=πγμ2=πζμ2=1,B 2b
    πβμs{βk,γk,ζk}=s{βk+1γk,ζk}πβμ,πγμs{βk,γk,ζk}=s{βk,γk+1,ζk}πγμB 2c
    andπζμs{βk,γk,ζk}=s{βk,γk,ζk+1}πζμ,kZ/(2Z).B 2d

    Footnotes

    Published by the Royal Society. All rights reserved.

    References