Arxiv:2103.01725V3 [Math.AG] 16 Jul 2021

arXiv:2103.01725v4 [math.AG] 28 Aug 2021 ⋆ ◦ E-mail E-mail ⋆ 00MteaisSbetClassiﬁcation Subject Mathematics 2010 roots unit transformations, Frobenius e words Key [email protected] upre npr yNFgatDMS-1954266 grant NSF by part in supported , [email protected] : [email protected] : aut fMteaisadMcais oooo ocwSta Moscow Lomonosov Mechanics, and Mathematics of Faculty ⋆ eateto ahmtc,Uiest fNrhCrln at Carolina North of University Mathematics, of Department iha pedxb tvnSperber Steven by appendix an with Zeutos euto modulo reduction equations, KZ : ⋆ OE NSLTOSO ZEQUATIONS KZ OF SOLUTIONS ON NOTES ocwCne fFnaetladApidMathematics Applied and Fundamental of Center Moscow eisieGr ,199 ocwGP1 Russia GSP-1, Moscow 119991 1, Gory Leninskiye eisieGr ,199 ocwGP1 Russia GSP-1, Moscow 119991 1, Gory Leninskiye MODULO 2 icn al inaoi,M 55,USA 55455, MN Minneapolis, Hall, Vincent 127 ◦ colo ahmtc,Uiest fMinnesota of University Mathematics, of School hplHl,N 79-20 USA 27599-3250, NC Hill, Chapel LXNE VARCHENKO ALEXANDER p s AND 33 1G5 41,3C0 32G20) 33C60, 14G10, (11G25, 13A35 : p AI LIMIT -ADIC p s , p s ◦ hpremti solutions, -hypergeometric n lxne Varchenko Alexander and ⋆ s ∞ → hplHill Chapel eUniversity te p ai limit, -adic 2 ALEXANDER VARCHENKO

Abstract. We consider the differential KZ equations over C in the case, when the hypergeometric solutions are one-dimensional hyperelliptic integrals of genus g. In this case the space of solutions of the differential KZ equations is a 2g-dimensional complex vector space. We also consider the same differential equations modulo ps, where p is an odd prime number and s is a positive integer, and over the field Qp of p-adic numbers. We describe a construction of polynomial solutions of the differential KZ equations modulo ps. These polynomial solutions have integer coefficients and are ps-analogs of the hyperelliptic integrals. We call them the ps-hypergeometric solutions. We consider the space s Mps of all p -hypergeometric solutions, which is a module over the ring of polynomial quasi- s constants modulo p . We study basic properties of Mps , in particular its natural filtration, and the dependence of Mps on s. We show that the p-adic limit of Mps as s → ∞ gives us a g-dimensional vector space of solutions of the differential KZ equations over the field Qp. The solutions over Qp are power series at a certain asymptotic zone of the KZ equations. In the appendix written jointly with Steven Sperber we consider all asymptotic zones of the KZ equations in the special case g = 1 of elliptic integrals. It turns out that in this case the p-adic limit of Mps as s → ∞ gives us a one-dimensional space of solutions over Qp at every asymptotic zone. We apply Dwork’s theory of the classical hypergeometric function over Qp and show that our germs of solutions over Qp defined at different asymptotic zones analytically continue into a single global invariant line subbundle of the associated KZ connection. Notice that the corresponding KZ connection over C does not have proper nontrivial invariant subbundles, and therefore our invariant line subbundle is a new feature of the KZ equations over Qp. Also in the appendix we follow Dwork and describe the Frobenius transformations of solutions of the KZ equations for g = 1. Using these Frobenius transformations we recover the unit roots of the zeta functions of the elliptic curves defined by the affine equations 2 × y = β x(x − 1)(x − α) over the finite field Fp. Here α, β ∈ Fp , α 6= 1. Notice that the same elliptic curves considered over C are used to construct the complex holomorphic solutions of the KZ equations for g = 1.

Contents 1. Introduction 3 2. KZ equations 6 3. Complex solutions 7 4. Solutions modulo ps 8 4.1. Leading terms 8 4.2. Quasi-constants 8 4.3. Solutions of system (2.1) modulo ps 9 4.4. ps-Hypergeometric solutions 9 4.5. Modules 10 5. Independence of modules from the choice of M 10 5.1. More general construction of solutions 10 5.2. More modules 11 6. Filtrations and homomorphisms 13 6.1. Reduction from modulo ps to modulo ps−m 13 6.2. Multiplication by pm 13 6.3. The composition of homomorphisms 13 6.4. Graded modules and homomorphisms 14 NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 3

7. Coefficients of solutions 14 7.1. Homogeneous polynomials 14 7.2. Formula for coefficients 14 8. Multiplication by p and Cartier-Manin matrix 15 8.1. Linear independence 15 8.2. Cartier-Manin matrices 15 8.3. Matrix of grcs,1 16 9. Change of variables 17 9.1. Change of the variable x 17 9.2. Change of variables z 18 9.3. Change of variables in the KZ equations 18 ˆ[lps−1] 9.4. Taylor expansion of Ips (u) 18 9.5. Taylor expansion of holomorphic solutions 20 9.6. Formal solutions over Qp and truncation 22 9.7. Example n = 3 22 10. p-Adic convergence 23 10.1. Teichmuller representatives 23 10.2. Approximation of binomial coefficients 24 10.3. Example n = 3, continuation 24 10.4. p-Adic convergence for arbitrary n 25 Appendix A. The case n = 3 and Dwork’s theory 27 A.1. Dwork on Legendre family 27 A.2. KZ equations 31 A.3. Solutions over C 31 A.4. Solutions as vectors of first derivatives 32 A.5. Six coordinate systems 33 A.6. Subbundle 35 A.7. Boundedness 35 A.8. More domains 36 A.9. Frobenius map on solutions of KZ equations 36 A.10. Eigenvalue ω(β(p−1)/2)(−1)(p−1)/2f(ω(α)) 37 ijk A.11. Approximation of analytic element (I )(i,j,k)∈S3 by rational functions 38 A.12. Further directions 40 References 40

1. Introduction 1.1. The KZ equations were introduced in [KZ] as the differential equations satisfied by conformal blocks on sphere in the Wess-Zumino-Witten model of conformal field theory. The solutions of the KZ equations in the form of multidimensional hypergeometric integrals were constructed more than 30 years ago, see [SV1]. The KZ equations and the hypergeometric solutions are related to many subjects in algebra, representation theory, theory of integrable systems, enumerative geometry.

The polynomial solutions of the KZ equations over the ﬁnite ﬁeld Fp of a prime number p of elements were constructed relatively recently in [SV2], see also [V4]-[V8], [RV1, RV2]. 4 ALEXANDER VARCHENKO

These solutions were called the Fp-hypergeometric solutions. The general problem is to un- derstand relations between the hypergeometric solutions of the KZ equations over C and the Fp-hypergeometric solutions and observe how the remarkable properties of hypergeometric solutions are reflected in the properties of the Fp-hypergeometric solutions. For example, the Fp-hypergeometric solutions inherit some determinant properties of the hypergeometric solutions and some Selberg integral properties, see [V8, RV1, RV2]. This program is in the first stages, where we consider essential examples and study the corresponding Fp-hypergeometric solutions by direct methods. In this paper we consider the differential KZ equations over C in the case, when the hypergeometric solutions are one-dimensional hyperelliptic integrals of genus g. In this case the space of solutions of the differential KZ equations is a 2g-dimensional complex vector space. We also consider the same differential equations modulo ps, where p is an odd prime number and s is a positive integer, and over the field Qp of p-adic numbers. We give a construction of polynomial solutions of the differential KZ equations modulo ps for positive integers s. We call such solutions the ps-hypergeometric solutions. This construction is a straightforward modification of the construction in [SV2] of polynomial solutions modulo p. s In this paper we consider the space Mps of all p -hypergeometric solutions, which is a module over the ring of polynomial quasi-constants modulo ps. We study basic properties of Mps , in particular its natural filtration, and dependence of Mps on s.

We show that the p-adic limit of Mps as s →∞ gives us a g-dimensional vector space of solutions of the differential KZ equations over the field Qp. The solutions over Qp are power series at a certain asymptotic zone of the KZ equations. This is the main result of the paper, see Lemma 9.5 and Theorem 10.7. 1.2. In the appendix written jointly with Steven Sperber we consider all six asymptotic zones of the KZ equations in the special case g = 1 of elliptic integrals. It turns out that in this case the p-adic limit of Mps as s →∞ gives us a one-dimensional space of solutions over Qp at every asymptotic zone. We apply Dwork’s theory of the classical hypergeometric function over Qp and show that our germs of solutions over Qp defined at different asymptotic zones analytically continue into a single global invariant line subbundle of the associated KZ connection. Notice that the corresponding KZ connection over C does not have proper nontrivial invariant subbundles, and therefore our invariant line subbundle is a new feature of the KZ equations over Qp. Following Dwork we show that our line subbundle is spanned at any point of the base by the germs of all solutions of the KZ equations bounded in their discs of convergence. This statement gives a definition of the line subbundle independent of asymptotic zones and analytic continuation. Also in the appendix we follow Dwork and describe the Frobenius transformations of solutions of the KZ equations for g = 1. Using these Frobenius transformations we recover the unit roots of the zeta functions of the elliptic curves defined by the affine equations 2 × y = β x(x − 1)(x − α) over the finite field Fp. Here α, β ∈ Fp ,α =6 1. Notice that the same elliptic curves considered over C are used to construct the complex holomorphic solutions of the KZ equations for g = 1. NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 5

In the end of Section A.10 we argue that the KZ equations for g = 1 contain more arithmetic information than the associated hypergeometric diﬀerential equation (1.2) for the hypergeometric function I(z) in (1.1), studied in [Dw].

1.3. Our p-adic limit of Mps as s →∞ is similar to the p-adic limit in the following classical example, see [Ig, Ma, Cl, BV1]. Consider the elliptic integral 1 ∞ dx ∞ −1/2 2 (1.1) I(z)= = zk . π 1 x(x − 1)(x − z) k Z Xk=0 It satisfies the hypergeometric differentialp equation (1.2) z(1 − z)I′′ + (1 − 2z)I′ − (1/4)I =0. The coefficients of the power series I(z) are p-adic integers and the power series I(z) converges p-adically for |z|p < 1, where |z|p is the p-adic norm of z ∈ Qp. One may show that for any positive integer s the polynomial

s (p −1)/2 s 2 (p − 1)/2 k (1.3) I s (z)= z (p −1)/2 k Xk=0 is a solution of the differential equation (1.2) modulo ps. Thus we get a sequence ∞ (I(ps−1)/2(z))s=1 of polynomials with integer coefficients, each of which is a solution of the differential equation (1.2) modulo ps, and the p-adic limit of the sequence, as s tends to ∞, is the p-adic power series solution I(z) of the differential equation (1.2). The ps-hypergeometric solutions of our differential KZ equations are analogs of the polynomials I(ps−1)/2(z) with an analogous p-adic limit. The difference is that the construction of the ps-hypergeometric solutions does not indicate the analogous p-adic limiting solutions I(z), and the analogous limiting p-adic power series solutions I(z) can be discovered only after rewriting the ps-hypergeometric solutions in a suitable asymptotic zone of the differential KZ equations. In the simplest example of our differential KZ equations, the p-adic solution is the 3-vector ∞ −1/2 −3/2 k +1 −1/2 (1.4) I(u ,u ) = u−3/2 , 1, uk , 1 2 1 k +1 k −1/2 − k −1/2 − k 2 Xk=0 ∞ s s while the sequence (I(ps−3)/2(u1,u2))s=1 of the p -hypergeometric solutions modulo p of the same equations is given by the formula

s (p −3)/2 s s (ps−3)/2 (p − 1)/2 (p − 3)/2 (1.5) I s (u ,u ) = u (p −3)/2 1 2 1 k +1 k Xk=0 k +1 (ps − 1)/2 × , 1, uk , (ps − 1)/2 − k (ps − 1)/2 − k 2 see Section 9.7. (ps−3)/2 ∞ The sum k=0 in (1.5) is the truncation of the sum k=0 in (1.4), similar to what happens in (1.1) and (1.3). A new feature appears when we compare the prefactor u−3/2 and P P 1 6 ALEXANDER VARCHENKO

(ps−3)/2 ∞ (ps−3)/2 ∞ the sequence of prefactors u1 s=1. As s →∞ the sequence of prefactors u1 s=1 tends p-adically to the prefactor u−3/2 multiplied by a Teichmuller constant on a suitable domain in Zp, where Zp is the ring of p-adic integers, see Section 10.3 and Theorem 10.5. 1.4. The paper is organized as follows. In Section 2 we define our system of KZ equations. In Section 3 we describe its complex solutions as hyperelliptic integrals. In Section 4 we describe the ps-hypergeometric solutions of our KZ equations modulo ps and define the filtered module s Mps of all p -hypergeometric solutions. In Section 5 we prove the independence of the module Mps from some arithmetic data involved in its definition. In Section 6 we discuss the properties of the operator Mps →Mps of multiplication by p. In Section 7 we calculate the coefficients of the Taylor expansion of the ps-hypergeometric solutions. In Section 8 we relate the operator Mps → Mps of multiplication by p and the Cartier-Manin matrix associated 2 with the hyperelliptic curve defined by the affine equation y =(x−z1) ··· (x−zn). In Section 9 we consider one of the asymptotic zones of our KZ equations. Using the coordinates in that asymptotic zone we describe the p-adic limit of the ps-hypergeometric solutions in Section 10. In Appendix A we apply Dwork’s theory in [Dw] to the case g = 1. In Section A.12 we discuss open problems related to the case of an arbitrary g. The author thanks Masha Vlasenko for numerous clarifying remarks on basics of the p-adic theory of geometric differential equations and comments on drafts of this paper. The author thanks Steven Sperber for collaboration on this project. The author thanks Pavel Etingof and Vadim Schechtman for useful discussions.

2. KZ equations Let g be a simple Lie algebra with an invariant scalar product. The Casimir element is

Ω= hi ⊗ hi ∈ g ⊗ g, i X n where (hi) ⊂ g is an orthonormal basis. Let V = ⊗i=1Vi be a tensor product of g-modules, κ ∈ C× a nonzero number. The differential KZ equations is the system of differential equations on a V -valued function I(z1,...,zn), ∂I 1 Ω = i,j I, i =1, . . . , n, ∂zi κ j6=i zi − zj X where Ωi,j : V → V is the Casimir operator acting in the ith and jth tensor factors, see [KZ, EFK]. This system is a system of Fuchsian first order linear differential equations. The equations are defined on the complement in Cn to the union of all diagonal hyperplanes. The object of our discussion is the following particular case. Let p be an odd prime number, n =2g +1 an odd positive integer, p > n > 2. We study the system of equations for a column vector I(z)=(I1(z), . . . , In(z)) :

∂I 1 Ωij (2.1) = I, i =1,...,n, I1(z)+ ··· + In(z)=0, ∂zi 2 zi − zj Xj6=i NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 7

where z =(z1,...,zn), the n × n-matrices Ωij have the form: .i .j . . i ··· −1 ··· 1 ··· (2.2) Ω = . . , ij  . .    j ··· 1 ··· −1 ···  . .   . .    and all other entries are zero. This joint system of differential and algebraic equations will be called the system of KZ equations in this paper. System (2.1) is the system of the differential KZ equations with parameter κ = 2 associated with the Lie algebra sl2 and the subspace of singular vectors of weight 2g − 1 of the tensor 2 ⊗(2g+1) power (C ) of two-dimensional irreducible sl2-modules, up to a gauge transformation, see this example in [V3, Section 1.1]. We consider system (2.1) over the field C. We also consider the same system of equations s modulo p and over the field Qp of p-adic numbers. 3. Complex solutions Consider the master function n −1/2 (3.1) Φ(x, z)= (x − za) a=1 Y and the column n-vector of hyperelliptic integrals Φ(x, z) (3.2) I(γ)(z)=(I (z),...,I (z)), I = dx . 1 n j x − z Z j The integrals Ij, are over an element γ of the first homology group of the algebraic curve with affine equation 2 y =(x − z1) ... (x − zn) . (γ) Starting from such γ, chosen for given values {z1,...,zn}, the vector I (z) can be analytically continued as a multivalued holomorphic function of z to the complement in Cn of the union of the diagonal hyperplanes zi = zj, i =6 j. Theorem 3.1. The vector I(γ)(z) is a solution of system (2.1). Theorem 3.1 is a classical statement. Much more general algebraic and differential equations satisfied by analogous multidimensional hypergeometric integrals were considered in [SV1]. Theorem 3.1 is discussed as an example in [V3, Section 1.1]. Proof. The theorem follows from Stokes’ theorem and the two identities: 1 Φ(x, z) Φ(x, z) ∂Φ (3.3) − + ··· + = (x, z) , 2 x − z1 x − zn ∂x ∂ 1 Ω Φ(x, z) Φ(x, z) ∂Ψi (3.4) − i,j ,..., = (x, z), ∂zi 2 zi − zj x − z1 x − zn ∂x Xj6=i 8 ALEXANDER VARCHENKO

where Ψi(x, z) is the column n-vector (0,..., 0, − Φ(x,z) , 0,..., 0) with the nonzero element x−zi at the i-th place. Theorem 3.2 ([V1, Formula (1.3)]). All solutions of system (2.1) have this form. Namely, the complex vector space of solutions of the form (3.2) is n − 1-dimensional. This theorem follows from the determinant formula for multidimensional hypergeometric integrals in [V1], in particular, from [V1, Formula (1.3)].

4. Solutions modulo ps

4.1. Leading terms. For a ring R denote R[z] = R[z1,...,zn]. For a positive integer t pt pt pt denote R[z ]= R[z1 ,...,zn ]. d1 dn Consider the lexicographical ordering of monomials z1 ...zn , so we have z1 > ··· > zn d1 dn and so on. For a nonzero polynomial f(z) = d1,...,dn ad1,...,dn z1 ...zn let fl(z) be the nonzero summand a zd1 ...zdn with the largest monomial zd1 ...zdn . We call f (z) the d1,...,dn 1 n P 1 n l d1 dn leading term of f(z), the coeﬃcient ad1,...,dn – the leading coeﬃcient, the monomial z1 ...zn – the leading monomial. s Let s be a positive integer. An element a ∈ Z/p Z has a unique presentation a = a0 + s−1 a1p + ··· + as−1p , where ai ∈ {0,...,p − 1}. An element a is invertible if and only if a0 =6 0.

Denote Fp = Z/pZ. s s n s n Let πs denote the homomorphisms Z → Z/p Z, Z[z] → (Z/p Z)[z], Z[z] → (Z/p Z)[z] s t s t and for t

4.2. Quasi-constants. We say that a polynomial f(z) ∈ Z[z] is a quasi-constant modulo ps if ∂f ∈ psZ[z] for i = 1,...,n. The quasi-constants modulo ps form a subring of Z[z] ∂zi ps denoted by Z[z]ps . For example (z1 + z2) ∈ Z[z]ps .

~~s−t d1 dn Lemma 4.1. As a Z-module the ring Z[z]ps is spanned by the monomials p z1 ...zn , t where t is the maximal integer such that t 6 s and p divides every d1,...,dn.~~

~~2 lps s−1 p p For example, z1 and p z1 z2 are such monomials.~~

d1 dn d1 dn Proof. Let f(z) = d cdz1 ...zn ∈ Z[z]ps . We show that each summand cdz1 ...zn is a 0 0 d1 dn multiple of a monomial of Lemma 4.1. Indeed, let cd0 z1 ...zn be the leading term of f(z). P s s−t Then all ﬁrst partial derivatives of it must lie in p Z[z] Hence cd0 ∈ p Z, where t is the t 0 0 maximal integer such that t 6 s and p divides every d1,...,dn. Subtracting the leading term from f(z) and repeating the reasoning we prove the lemma.

~~s t Lemma 4.2. Let f(z) be a quasi-constant modulo p and t ∈ Z>0. Then p f(z) is a quasi- constant modulo pr for any 1 6 r 6 s + t.~~

~~The rings of quasi-constants form a decreasing ﬁltration, Z[z]p ⊃ Z[z]p2 ⊃ .... NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 9~~

4.3. Solutions of system (2.1) modulo ps. We say that a column n-vector I(z) ∈ Z[z]n s of polynomials with integer coeﬃcients is a solution of system (2.1) modulo p , if πsI(z) ∈ (Z/psZ)[z]n satisﬁes system (2.1). Lemma 4.3. Let I(z) be a solution of system (2.1) modulo ps. t r (i) Let t ∈ Z>0. Then p I(z) is a solution of system (2.1) modulo p for any 1 6 r 6 s+t. (ii) Let f(z) be a quasi-constant modulo ps. Then f(z)I(z) is a solution of system (2.1) modulo ps. (iii) Let 1 6 t

4.4. ps-Hypergeometric solutions. Let M be the least positive integers such that 1 (4.1) M ≡ − (mod ps) . 2 We have ps − 1 p − 1 M = = 1+ p + ··· + ps−1 . 2 2 Introduce the master polynomial n M (4.2) Φps (x, z) = (x − zi) ∈ Z[x, z]. i=1 Y Let

Φps (x, z) Φps (x, z) i i (4.3) P s (x, z) = ,..., = P s (z) x , p x − z x − z p 1 n i X i where Pps (x, z) is a column n-vector of polynomials in x, z1,...,zn and Pps (z) are n-vectors of polynomials in z1,...,zn with coeﬃcients in Z. For a positive integer l, denote [lps−1] lps−1 Ips (z) = Pps (z) .

s [lp −1] n Theorem 4.4. For any positive integer l, the vector of polynomials Ips (z) ∈ Z[z] is a solution of system (2.1) modulo ps. Proof. We have the following modiﬁcations of identities (3.3), (3.4):

Φ s (x, z) Φ s (x, z) ∂Φ s (4.4) M p + ··· + p = p (x, z) , x − z1 x − zn ∂x i ∂ Ω Φ s (x, z) Φ s (x, z) ∂Ψ s (4.5) + M i,j p ,..., p = p (x, z), ∂zi zi − zj x − z1 x − zn ∂x Xj6=i i Φps (x,z) where Ψ s (x, z) is the column n-vector (0,..., 0, − , 0,..., 0) with the nonzero element p x−zi at the i-th place. The theorem follows from identities (4.4), (4.5). 10 ALEXANDER VARCHENKO

Remark. In [SV2] it was explained on how to construct polynomial solutions modulo p of an arbitrary system differential KZ equations, associated with any Kac-Moody algebra and any tensor product of highest weight representations. The same construction gives polynomial solutions modulo ps. The details will be provided elsewhere. s 6 ps−1 The range for the index l is defined by the inequalities 0 < lp − 1 n 2 . Hence [lps−1] l =1,...,g. The solutions Ips (z), l =1,...,g, given by this construction, will be called the ps-hypergeometric solutions in Z[z]n. For t = 1,...,s − 1 and l = 1,...,g, the vector t s−t [lp −1] s p Ipt (z) is a solution of system (2.1) modulo p , see Lemma 4.3. Such solutions also will be called ps-hypergeometric solutions in Z[z]n. 4.5. Modules. Consider the increasing filtration 0 1 s−1 s (4.6) 0= Mps ⊂Mps ⊂···⊂Mps ⊂Mps = Mps , where t g t s−r [lpr−1] (4.7) Mps = πs cr,l(z) p Ipr (z) | cr,l(z) ∈ Z[z]pr , r=1 n X Xl=1 o s n t = 1,...,s. We have Mps ⊂ (Z/p Z)[z] . Every element of Mps is a polynomial solution s of system (2.1) with coefficients in Z/p Z, see Lemma 4.3. The set Mps is a module over s the ring Z[z]ps of quasi-constants modulo p , where f(z) ∈ Z[z]ps acts by multiplication by t πsf(z). Each Mps is an Z[z]ps -submodule of Mps . t t Each Mps is also a module over the larger ring Z[z]pt of quasi-constants modulo p , where f(z) ∈ Z[z]pt acts by multiplication by πsf(z). s s n The elements of Mps will be called the p -hypergeometric solutions in (Z/p Z)[z] .

~~5. Independence of modules from the choice of M~~

~~5.1. More general construction of solutions. For i = 1,...,n, let Mi be a positive integer such that 1 (5.1) M ≡ − (mod ps) . i 2~~

Denote M~ =(M1,...,Mn). Consider the master polynomial n Mi (5.2) Φ(x, z, M~ )= (x − zi) ∈ Z[x, z], i=1 Y and the Taylor expansion Φ(x, z, M~ ) Φ(x, z, M~ ) P (x, z, M~ )= ,..., = P i(z, M~ ) xi, x − z x − z 1 n i X i where P (z, M~ ) are n-vectors of polynomials in z1,...,zn with coeﬃcients in Z. For a positive integer l, denote

~~s s I[lp −1](z, M~ ) = P lp −1(z, M~ ) . NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 11~~

Theorem 5.1. For any positive integers l, t, t 6 s, the vector of polynomials I[lpt−1](z, M~ ) ∈ Z[z]n is a solution of system (2.1) modulo pt. Proof. The theorem follows from straightforward modifications of identities (4.4), (4.5). 5.2. More modules. Consider the increasing filtration 0 ~ 1 ~ s−1 ~ s ~ ~ (5.3) 0= Mps (M) ⊂Mps (M) ⊂···⊂Mps (M) ⊂Mps (M)= Mps (M) , where t t ~ s−r [lpr−1] ~ (5.4) Mps (M) = πs cr,l(z) p I (z, M) | cr,l(z) ∈ Z[z]pr , r=1 > n X Xl 1 o s n t = 1,...,s. We have Mps (M~ ) ⊂ (Z/p Z)[z] . Every element of Mps (M~ ) is a polynomial s solution of system (2.1) with coefficients in Z/p Z, see Lemma 4.3. The set Mps (M~ ) is s a module over the ring Z[z]ps of quasi-constants modulo p , where f(z) ∈ Z[z]ps acts by t ~ ~ multiplication by πsf(z). Each Mps (M) is an Z[z]ps -submodule of Mps (M). t ~ t Each Mps (M) is also a module over the larger ring Z[z]pt of quasi-constants modulo p , where f(z) ∈ Z[z]pt acts by multiplication by πsf(z).

Theorem 5.2. Filtration (5.3) does not depend on the choice of M~ =(M1,...,Mn), satisfying congruences (5.1). Moreover, ﬁltration (5.3) coincides with ﬁltration (4.6). For s = 1 the statement is [SliV, Theorem 3.1].

Proof. First we show that Mps (M~ ) and filtration (5.3) do not depend on the choice of M~ . ~ ~ ′ ′ ′ Let M =(M1,...,Mn), M =(M1,...,Mn) be two vectors satisfying congruences (5.1). We ~ ′ > ~ ′ > ps−1 ps−1 say that M M if Mi Mi for all i. The vector ( 2 ,..., 2 ) is the minimal vector with respect to this partial order. To show that Mps (M~ ) and filtration (5.3) do not depend on the choice of M~ it is enough to show that the filtrations are the same for a vector M~ and ~ ′ ~ ps−1 for a vector M = M +(0,..., 0, 2 , 0,..., 0), where the nonzero element stays at the j-th position for some j. Then s p s s s p s P (x, z, M~ ′) = P (x, z, M~ ) · (x − z )p = P (x, z, M~ ) (−1)p −a xazp −a. j a j a=0 X ~ i ~ i ~ ′ i ~ ′ i Recall that P (x, z, M) = i P (z, M) x , P (x, z, M ) = i P (z, M ) x , For any r 6 s and l we have P s P p s r s p s r (5.5) P lp −1(z, M~ ′)= (−1)p −a zp −aP lp −a−1(z, M~ ). a j a=0 X We are interested in this formula, since I[lpr−1](z, M~ ′)= P lpr−1(z, M~ ′) is a solution of system (2.1) modulo pr. c s s−c Lemma 5.3. Let b, c ∈ Z>0 be such that bp 6 p , p 6| b. Then p is the maximal power of ps p dividing bpc . 6 s ps s ps−1 ps−1 Proof. For a p we have a a = p a−1 , and p 6| a−1 by Lucas’ theorem, [Lu]. 12 ALEXANDER VARCHENKO

Lemma 5.4. Let a ∈{0,...,ps}, a = bpc, p 6| b. Consider the vector s s p s r V =(−1)p −a zp −aP lp −a−1(z, M~ ) a j appearing in (5.5). If lpr 6 a, then V = 0. For lpr > a, we write lpr − a = vpu, where u = min(r, c). Then

~~u (5.6) V = d(z) pr−uI[vp −1](z, M~ ),~~

ps−a ps−a ps r−u u where d(z)= zj (−1) a /p is a quasi-constant modulo p . Proof. The lemma follows from Lemma 5.3. r ~ ′ r ~ Corollary 5.5. For any r =1,...,s, we have Mps (M ) ⊂Mps (M). r ~ ′ r ~ Lemma 5.6. For any r =1,...,s, we have Mps (M ) ⊃Mps (M). r ~ ′ s−r Proof. Let w be the greatest integer such that wp 6 degt P (t, z, M). Then w = w + p ′ r ~ ′ is the greatest integer such that w p 6 degt P (t, z, M ). Comparing the coeﬃcients in (5.5) and using Lemma 5.5, we observe that for any l =1,...,w we have

s−r r r (5.7) I[(l+p )p −1](z, M~ ′)= I[lp −1](z, M~ ) l−1 r−1 [mpr−1] r−k [mpk−1] + cr,m(z)I (z, M~ )+ ck,m(z) p I (z, M~ ), m=1 m>1 X Xk=1 X [lpr−1] ~ where ci,j(z) ∈ Z[z]pi . This triangular system of equations with respect to I (z, M), l =1,...,w, can be written as r−1 [lpr−1] ~ ′ [mpr−1] ~ ′ ′ r−k [mpk−1] ~ (5.8) I (z, M)= cr,m(z)I (z, M )+ ck,m(z) p I (z, M), m>1 m>1 X Xk=1 X ′ l = 1,...,w, for suitable ci,j(z) ∈ Z[z]pi . Applying the previous reasoning to the sum r−1 r−k [mpk−1] ~ k=1 m>1 ck,m(z) p I (z, M), we obtain P P r−1 [lpr−1] ~ ′ [mpr−1] ~ ′ ′′ r−k [mpk−1] ~ ′ (5.9) I (z, M)= cr,m(z)I (z, M )+ ck,m(z) p I (z, M ), m>1 m>1 X Xk=1 X ′′ l =1,...,w, for suitable ci,j(z) ∈ Z[z]pi . This proves the lemma.

~~Corollary 5.7. The module Mps (M~ ) and ﬁltration (5.3) do not depend on the choice of~~

~~M~ =(M1,...,Mn), satisfying congruences (5.1).~~

min ps−1 ps−1 min ~ s ~ s Lemma 5.8. Let M =( 2 ,..., 2 ). Then Mp (M )= Mp . Proof. The proof of the lemma is a straightforward modiﬁcation of the proof of Corollary 5.5 and Lemma 5.6. Theorem 5.2 is proved. NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 13

6. Filtrations and homomorphisms 6.1. Reduction from modulo ps to modulo ps−m. If I(z) is a polynomial solution of system (2.1) modulo ps, then I(z) is also a polynomial solution of system (2.1) modulo ps−m for any 1 6 m

~ − ~ (6.1) rs,s−m : Mps (M) →Mps m (M), s s s−r [lpr−1] s−r [lpr−1] πs cr,l(z) p I (z, M~ ) 7→ πs−m cr,l(z) p I (z, M~ ) , r=1 > r=m+1 > X Xl 1 X Xl 1 where M~ is a vector with coordinates satisfying congruences (5.1). See these sums in (5.4). − In the last sum we have ps−rI[lpr−1](z, M~ ) = ps−m−(r−m)I[(lpm)pr m−1](z, M~ ) and a solution − I[lpr−1](z, M~ ) modulo pr also can be considered as a solution I[(lpm)pr m−1](z, M~ ) modulo pr−m. r ~ ~ For any r = 1,...,s, the submodule Mps (M) ⊂ Mps (M) is mapped by rs,s−m to the r−m ~ submodule Mps−m (M). The induced map r ~ r−m ~ (6.2) rs,s−m : Mps (M) →Mpt (M) is a homomorphism of Z[z]pr -modules. Thus the map (6.1) is a homomorphism of filtered modules decreasing the index of filtration by m. By Theorem 5.2 we have Mps (M~ ) = Mps . Hence homomorphism (6.1) also can be considered as a homomorphism of filtered modules,

~~(6.3) rs,s−m : Mps →Mps−m , decreasing the index of ﬁltration by m. It is rather nontrivial to write a formula for this map in terms of the generators I[lpr−1](z) of these modules.~~

6.2. Multiplication by pm. If I(z) is a polynomial solution of system (2.1) modulo ps, then for any positive integer m the polynomial pmI(z) is a polynomial solution of system (2.1) modulo ps+m. In particular, multiplication by pm deﬁnes a map t t (6.4) ps,s+m : Mps →Mps+m , t g t g s−r [lpr−1] s+m−r [lpr−1] πs cr,l(z) p Ipr (z) 7→ πs+m cr,l(z) p Ipr (z) r=1 r=1 X Xl=1 X Xl=1 for any t =1,...,s. See these sums in (4.7). Clearly this map is an isomorphism of ﬁltered Z[z]pt -modules.

~~6.3. The composition of homomorphisms. For m ~~s. 14 ALEXANDER VARCHENKO~~~~

r s−r [lp −1] As we know, the module Mps is generated by the elements πs(p Ipr (z)), r =1,...,s, l =1,... , g. For l =1,...,g, we have s−1 g [lps−1] [lps−1] l,s s−r [kpr−1] (6.6) cs,1 : πs(I (z)) 7→ πs(pI (z)) = πs cr,k(z) p Ipr (z) r=1 X Xk=1 l,s for suitable coefficients cr,k(z) ∈ Z[z]pr . l,s The set of the coefficients (cr,k(z))l,s,r,k determines the homomorphisms cs,m for all s, m. l,s In what follows we shall describe the coefficients cs−1,k(z) for all l,s,k, see Theorem 8.4. 6.4. Graded modules and homomorphisms. Denote

~~s t t t t−1 (6.7) grMps = ⊕t=1 grMps , grMps = Mps Mps .~~

. t t pt Lemma 6.1. For any t = 1,...,s, the action of Z[z]pt on Mps makes grMps an Fp[z ]- t module. Multiplication by p on Mps induces a homomorphism t t−1 (6.8) grcs,1 : grMps → grMps pt of Fp[z ]-modules.

7. Coefficients of solutions s [lps−1] [lp −1] 7.1. Homogeneous polynomials. For l =1,...,g, the solution I (z)=(I1 ,..., [lps−1] In ) is a homogeneous polynomial in z of degree ps − 1 ps − 1 (7.1) δ = (2g + 1) − lps =(g − l)ps + − g . l 2 2 p−1 Notice that (−1)δl =(−1)s 2 + l .

ps−1 s 7.2. Formula for coeﬃcients. Recall that M = 2 . Projection of this integer to Z/p Z is invertible. Let s s s [lp −1] [lp −1] d1 dn [lp −1] Zn I (z)= Id1,...,dn z1 ...zn , Id1,...,dn ∈ . d1X,...,dn Lemma 7.1 ([V8, Lemma 3.1]). We have n [lps−1] M d1 dn (7.2) I = (−1)δl 1 − ,..., 1 − . d1,...,dn d M M i=1 i Y The sum of coordinates of this vector is divisible by ps. Lemma 7.2 (cf. [V8, Theorem 6.1]). For l =1,...,g, the leading term of the ps-hypergeometric solution I[lps−1](z) is

[lps−1] M l (7.3) I (z) = (−1)δl 0,..., 0, , 1,..., 1 zM ...zM zM−l , l l M 1 2g−2l 2g−2l+1 where 0 is repeated 2g − 2l times and 1 is repeated 2l times. Proof. The lemma follows from Lemma 7.1. NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 15

Z sZ M M l Lemma 7.3. The projections to /p of the integers l , l M are invertible. M Proof. The invertibility of l follows from Lucas’ theorem, [Lu]. 8. Multiplication by p and Cartier-Manin matrix 8.1. Linear independence. Lemma 8.1. The projections of the ps-hypergeometric solutions I[lps−1](z) ∈ Z[z]n, l = n 1,...,g, to Fp[z] are linearly independent over Fp[z], that is, if g [lps−1] n (8.1) cl(z)I (z) ∈ pZ[z] Xl=1 for some cl(z) ∈ Z[z], then all cl(z) ∈ pZ[z].

n n Proof. Recall that the projection Z[z] → Fp[z] is denoted by π1. By Lemma 7.3, the leading [lps−1] coefficient of π1(cl(z)I (z)) equals the product of the leading coefficient of π1(cl(z)) and [lps−1] the leading coefficient of π1(I (z)), if π1(cl(z)) is nonzero. In that case the leading [lps−1] l coefficient of π1(cl(z)I (z)) is a nonzero multiple of the nonzero vector π1((0,..., 0, M , 1, . . . , 1)). If relation (8.1) holds and some of the coefficients cl(z) have nonzero projections π1(cl(z)), then for several values of such indices l the sum of the corresponding leading coefficients has to l be equal to zero, which is impossible due to the fact that the vectors π1((0,..., 0, M , 1,..., 1)) are linear independent over Fp. Corollary 8.2. The projections of the ps-hypergeometric solutions I[lps−1](z) ∈ Z[z]n, l = n ps 1,...,g, to Fp[z] are linearly independent over Fp[z ]. Denote by t t (8.2) grt : Mps → grMps

~~s−t [lpt−1] the natural projection. Then the elements grt(πs(p I (z))), l =1,...,g, generate the pt t Fp[z ]-module grMps .~~

~~pt t Corollary 8.3. For t = 1,...,s, the Fp[z ]-module grMps is a free module of rank g with s−t [lpt−1] a basis grt(πs(p I (z))), l =1,...,g.~~

~~pt t Denote the basis vectors of the Fp[z ]-module grMps by~~

~~l s−t [lpt−1] (8.3) vs,t := grt(πs(p I (z))), l =1, . . . , g.~~

~~8.2. Cartier-Manin matrices. Let~~

~~f(x, z)=(x − z1) ... (x − zn), n =2g +1. Consider the hyperelliptic curve X deﬁned by the aﬃne equation 2 y =(x − z1) ... (x − zn). 16 ALEXANDER VARCHENKO~~

− 1 xi 1dx Consider the space Ω (X) of regular 1-forms on X with basis y , i =1,...,g. Define the Cartier map C : Ω1(X) → Ω1(X) as follows. We have xi−1dx xi−1yp−1dx xi−1f(x)(p−1)/2dx = = . y yp−1y yp i−1 (p−1)/2 j j Let x f(x, z) = j ci (z)x . Define g P xi−1dx xj−1dx C : 7→ cjp−1(z) , y i y j=1 X j g see [AH]. The map C is identified with the g × g-matrix (Ci (z))i,j=1 , j jp−1 Ci (z) = ci (z).

8.3. Matrix of grcs,1. Recall that multiplication of solutions by p defines a map t t−1 (8.4) grcs,1 : grMps → grMps , t pt l g t−1 pt−1 where grMps is a free Fp[z ]-module with a basis (vs,t)l=1 and grMps is a free Fp[z ]- l g pt module with a basis (vs,t−1)l=1, see (8.3). The map (8.4) is a homomorphism of Fp[z ]- modules. pt−1 Theorem 8.4. The matrix of grcs,1 is the Cartier-Manin matrix C(z ), g l m l pt−1 (8.5) grcs,1 : vs,t 7→ vs,t−1Cm(z ), l =1, . . . , g. m=1 X t s−t+2 n s−t [lp −1] Proof. The problem is to express modulo p Z[z] the element p·p Ipt (z) in terms of t−1 t s−t+1 [mp −1] [lp −1] the elements p Ipt−1 (z), m =1,...,g. In other words, we need to express Ipt (z) t−1 [mp −1] n in terms of Ipt−1 (z), m =1, . . . , g, modulo pZ[z] . t [lp −1] lpt−1 By definition, the vector Ipt (z) is the coefficient of x in the Taylor expansion of the polynomial n pt−1(p−1)/2 Ppt (x, z)= Ppt−1 (x, z) (x − zi) , i=1 t−1 Y [mp −1] mpt−1−1 while the vector Ipt−1 (z) is the coefficient of x in the Taylor expansion of the polynomial Ppt−1 (x, z), see notations in Section 4.4. We have n n − pt−1(p−1)/2 pt−1 pt 1 (p−1)/2 (x − zi) ≡ (x − zi ) mod p . i=1 i=1 t Y t−1 Y [lp −1] g [mp −1] l pt−1 Hence Ipt (z) ≡ m=1 Ipt−1 (z) Cm(z ) mod p. Theorem 8.4 is proved. t t−m Corollary 8.5. TheP matrix of grcs,m : grMps → grMps is the product of Cartier-Manin −1 −2 − matrices C(zpt )C(zpt ) ···C(zpt m ). Moreover, this statement, applied to the map s 1 1 grcs,s−1 : grMps → grMps = Mps , NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 17

lps−1 s can be reformulated as follows. For any l = 1,...,g, the solution Ips (z) modulo p of n n system (2.1), projected to Fp[z] , equals the projection to Fp[z] of the solution g s−2 s−1 [m1p−1] m2 p ms−1 p l p Ip (z)Cm1 (z ) ···Cms−2 (z )Cms−1 (z ) m1,...,m −1=1 Xs modulo p of system (2.1). See these sums in [V5, Section 8].

~~9. Change of variables~~

~~9.1. Change of the variable x. Change the variable x and set x = v + zn. Then n−1 (ps−1)/2 (ps−1)/2 (9.1) Φ˜ ps (v, z) := Φps (v + zn, z)= (v − (zi − zn)) v . i=1 Y Let~~

(9.2) P˜ps (v, z): = Pps (v + zn, z) ˜ ˜ ˜ Φps (v, z) Φps (v, z) Φps (v, z) i i = ,..., , = P˜ s (z) v , v − (z − z ) v − (z − z ) v p 1 2 n−1 n i X ˜i where Pps (z) are n-vectors of polynomials in z with integer coefficients. For a positive integer l, denote ˜[lps−1] ˜lps−1 (9.3) Ips (z) := Pps (z) . ˜[lps−1] The polynomial Ips (z) is nonzero if l = 1,...,g. Notice that every polynomial ˜[lps−1] Ips (z) is a function of differences zi − zn, i =1,...,n − 1. Consider the increasing filtration ˜ 0 ˜ 1 ˜ s−1 ˜ s ˜ (9.4) 0= Mps ⊂ Mps ⊂···⊂ Mps ⊂ Mps = Mps , where s g ˜ s−r ˜[lpr−1] (9.5) Mps = πs cr,l(z) p Ipr (z) | cr,l(z) ∈ Z[z]pr , n r=1 l=1 o Xt Xg ˜ t s−r ˜[lpr−1] (9.6) Mps = πs cr,l(z) p Ipr (z) | cr,l(z) ∈ Z[z]pr , r=1 n X Xl=1 o t =1,...,s.

s ˜[lp −1] n Theorem 9.1. For any l, the vector of polynomials Ips (z) ∈ Z[z] is a solution of system s ˜ t t (2.1) modulo p . For any t =1,...,s we have Mps = Mps . Proof. The proof is the same as the proof of Theorem 5.2 and the proof of [V5, Lemma 5.2]. In the proof of Theorem 9.1 the following Lemma 9.2 is used instead of Lemma 5.3.

~~mpr+lps−1 s−r Lemma 9.2. Let r =0,...,s − 1 and m6| p, then lps−1 is divisible by p . 18 ALEXANDER VARCHENKO~~

~~mpr+lps−1 mpr+lps−1 lps mpr+lps−1 Proof. We have lps−1 = mpr = mpr mpr−1 .~~

~~9.2. Change of variables z. We introduce the new variables u1,...,un by the formulas :~~

~~z2 − zn zn−1 − zn (9.7) u1 = z1 − zn, u2 = , ... un−1 = , un = z1 + ··· + zn, z1 − zn zn−2 − zn or~~

~~zi − zn = u1 ··· ui, i =1,...,n − 1, z1 + ··· + zn = un .~~

~~For any l,s we denote u =(u1,...,un−1), ˆ[lps−1] ˜[lps−1] (9.8) Ips (u) := Ips (z(u)) .~~

s ˆ[lp −1] n Each Ips (u) is an n-vector of polynomials in u with coeﬃcients in Z . s ˆ[lp −1] s Each Ips (u) is a solution of system (2.1) modulo p , in which the change of variables z = z(u) is performed.

9.3. Change of variables in the KZ equations. It is known that system (2.1) of the differential KZ equations has suitable asymptotic zones with appropriate local coordinates, in which the differential KZ equations have singularities only at the coordinate hyperplanes. See a definition of the asymptotic zones, for example, in [V2]. The coordinates u defined in (9.7) are local coordinates in one of the asymptotic zones. In these coordinates, system (2.1) takes the form,

~~∂I ∂I 1 Ωi (9.9) =0, = + Regi(u) I, i =1,...,n − 1, ∂un ∂ui 2 ui I1 + ··· + In =0,~~

~~where Ωi = i6k~~

s n−1 i (p −1)/2 (ps−1)/2 (9.10) Φˆ ps (v,u) := Φ˜ ps (v, z(u)) = v − uj v , i=1 j=1 ! Y Y ˆ ˆ ˆ Φps (v,u) Φps (v,u) Φps (v,u) i i (9.11) Pˆ s (v,u) := ,..., , = Pˆ s (u) v , p v − u v − u ··· u v p 1 1 n−1 i X NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 19

î where Pps (u) are n-vectors of polynomials in u with coefficients in Z. For a positive integer l, we have ˆlps−1 ˆ[lps−1] (9.12) Pps (u) = Ips (u) , ˆ[lps−1] where Ips (u) is defined in (9.8). For l =1,...,g, denote n−2l ps−1 l,s δl −l (9.13) u = (−1) (u1 ··· un−2l) (u1 ··· ui) 2 , i=1 Y or ps−1 ps−1 ps−1 ps−1 g,s δg 2 −g g−1,s δg−1 3 2 −g+1 2 2 −g+1 2 −g+1 u =(−1) u1 , u =(−1) u1 u2 u3 ,... ps−1 ps−1 ps−1 1,s δ1 (n−2) 2 −1 (n−3) 2 −1 2 −1 u =(−1) u1 u2 ··· un−2 , see (−1)δl in (7.1). Denote s p −1 2l (9.14) Cl,s = 2 0,..., 0, , 1,..., 1 , l ps − 1 where 0 is repeated 2g − 2l times and 1 is repeated 2l times, cf. formula (7.3). ˆ[lps−1] Theorem 9.3. For l =1,...,g, the polynomial Ips (u) has the following form, ˆ[lps−1] l,s l,s (9.15) Ips (u) = u T (u) , l,s l,s l,s T (u) = (T1 (u),...,Tn (u)), l,s with coordinates Tj defined as follows. If j =1,...,n − 1, then s n−1 s l,s,j p −3 p −1 l,s 2 2 (9.16) Tj = uj+1 ··· un−2l aj ai X i=1Y,i6=j n−2l−1 2l−1 ai an−2l+i × (ui+1 ··· un−2l) (un−2l+1 ··· un−2l+i) , i=1 i=1 Y Y l,s,j 6 6 ps−1 where the summation is over all a1,...,an−1 ∈ Z, 0 ai 2 , such that a1 + ··· + an−2l = an−2l+1 + ··· + an−1 + l − 1, if j 6 n − 2l; and such that a1 + ··· + an−2l = P an−2l+1 + ··· + an−1 + l, if n − 2l < j 6 n − 1; (9.17) n−1 s n−2l−1 2l−1 l,s,n p −1 T l,s = 2 (u ··· u )ai (u ··· u )an−2l+i , n a i+1 n−2l n−2l+1 n−1 i=1 i i=1 i=1 X Y Y Y l,s,n 6 6 ps−1 where the summation is over all a1,...,an−1 ∈ Z, 0 ai 2 , such that a1 + ··· + an−2l = an−2l+1 + ··· + an−1 + l. The constant term ofPT l,s(u) equals Cl,s.

~~Notice that the factor uj+1 ··· un−2l in (9.16) equals 1 if j > n − 2l. 20 ALEXANDER VARCHENKO~~

Proof. Make the change of variables v = wu1 ··· un−2l in (9.11), ◦ ˆ î i ◦,i i Pps (w,u) := Pps (wu1 ··· un−2l,u)= Pps (u)(wu1 ··· un−2l) =: Pps (u) w . i i X X Hence ˆ[lps−1] −(lps−1) ◦,lps−1 Ips (u) = u1 ··· un−2l Pps (u) . ◦ We transform the factors in the polynomial Pps(w,u) as follows. For any positive integers k and i 6 n − 2l we write k k k (9.18) (wu1 ··· un−2l − u1 ··· ui) = (u1 ··· ui) (wui+1 ··· un−2l − 1) k k = (u ··· u )k (−1)k−a (wu ··· u )a, 1 i a i+1 n−2l a=0 X and if i > n − 2l, we write k k k (9.19)(wu1 ··· un−2l − u1 ··· ui) = (wu1 ··· un−2) (1 − un−2l+1 ··· ui/w) k k = (wu ··· u )k (−1)a (u ··· u /w)a. 1 n−2 a n−2l+1 i a=0 X ps−1 ps−1 Notice that for factors in (9.11), we have k = 2 or k = 2 − 1. This explains the binomial coefficients in (9.16) and (9.17). We prove formula (9.16) for j = 1, the proof for other values of j is similar. −(lps−1) ◦,lps−1 Our goal is to calculate the first coordinate of the vector u1 ··· un−2 Pps (u). That is we need to calculate the coefficient of wlps−1 in s−1 −(lps−1) p −1 ps−1 ps−1 ps−1 2 2 −1 2 2 u1 ··· un−2 u1 (wu2 ··· un−2l − 1) (u1u2) (wu3 ··· un−2l − 1) ··· ps−1 ps−1 ps−1 ps−1 ··· (u1··· un−2l) 2 (w − 1) 2 (wu1 ··· un−2l) 2 (1 − un−2l+1/w) 2 ··· ps−1 ps−1 ps−1 ··· (wu1 ··· un−2l) 2 (1 − un−2l+1 ··· un−1/w) 2 (wu1 ··· un−2l) 2 , which is the same as the coefficient of wl−1 in s− s− l,s p 1 −1 p 1 (9.20) u2 ··· un−2l u (wu2 ··· un−2l − 1) 2 (wu3 ··· un−2l − 1) 2 ··· ps−1 ps−1 ps−1 ··· (w − 1) 2 (1 − un−2l+1/w) 2 ··· (1 − un−2l+1 ··· un−1/w) 2 . Expanding the binomials we obtain formula (9.16) for j = 1. The constant term of T l,s(u) is given by the summands in (9.16) and (9.17), corresponding to a1 = ··· = an−2l−1 = an−2l+1 = ··· = an−1 = 0 and j = n − 2l,...,n. Theorem 9.3 is proved.

9.5. Taylor expansion of holomorphic solutions. Recall the multivalued holomorphic solutions of system (2.1) described in Section 3, Φ(x, z) Φ(x, z) I(γ)(z)= ,..., dx . γ x − z1 x − zn Z NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 21

We make the same changes of variables in the integrals I(γ)(z) as we did in the previous sections. Namely, ﬁrst we change the integration variable x and set x = v + zn, then we make the change of variables z and set z = z(u). The resulting integral is Φ(ˆ v,u) Φ(ˆ v,u) Iˆ(γ)(u) = ,..., dv , γ v − u1 v Z −1/2 ˆ n−1 i −1/2 where Φ(v,u)= i=1 v − j=1 uj v . For l =1,...,g,Q we change Q the integration variable v and set v = wu1 ...un−2l. Then n−2l (l−n/2)πi −l −1/2 Φ(ˆ wu1 ...un−2l,u)= e (u1 ··· un−2l) (u1 ··· ui) i=1 Y × (1 − wu2 ··· un−2l)(1 − wu3 ··· un−2l) ··· (1 − w)

−1/2 −l × (1 − un−2l+1/w) ··· (1 − un−2l+1 ··· un−1/w) w . Choose the integration cycle γ = γl to be the circle |w| =1/2 oriented counter-clockwise. We assume that all the variables u2,...,un−1 lie inside the circle. We ﬁx the branch of the function

(9.21) (1 − wu2 ··· un−2l)(1 − wu3 ··· un−2l) ··· (1 − w) −1/2 × (1 − un−2l+1/w) ··· (1 − un−2l+1 ··· un−1/w) over the circle by choosing the argument of the function in (9.21) at w = 1/2, u2 = ··· = enπi/2 un−1 = 0 to be 0. We multiply the circle with the chosen branch of the integrand by 2πi . This ﬁnishes the description of γ1. The resulting integral is enπi/2 Φ(ˆ wu ...u ,u) Φ(ˆ wu ...u ,u) (γl) 1 n−2l 1 n−2l Iˆ (u) = ,..., u1 ...un−2l dw . 2πi |w|=1/2 wu1 ...un−2l − u1 wu1 ...un−2l Z Denote n−2l l −l −1/2 (9.22) u := (u1 ··· un−2l) (u1 ··· ui) . i=1 Y Theorem 9.4. For l =1,...,g, the function Iˆ(γl)(u) has the following form,

ˆ(γl) l l l l l (9.23) I (u) = u T (u) , T (u)=(T1(u),...,Tn(u)), l with coordinates Tj deﬁned as follows. If j =1,...,n − 1, then n−1 l,j −3 −1 l 2 2 (9.24) Tj = uj+1 ··· un−2l aj ai X i=1Y,i6=j n−2l−1 2l−1 ai an−2l+i × (ui+1 ··· un−2l) (un−2l+1 ··· un−2l+i) , i=1 i=1 Y Y 22 ALEXANDER VARCHENKO

l,j where the summation is over all a1,...,an−1 ∈ Z>0 such that a1 +···+an−2l = an−2l+1 + ··· + an−1 + l − 1, if j 6 n − 2l; and such that a1 + ··· + an−2l = an−2l+1 + ··· + an−1 + l, if n − 2l < j 6 n − 1; P n−1 n−2l−1 2l−1 l,n −1 (9.25) T l = 2 (u ··· u )ai (u ··· u )an−2l+i , n a i+1 n−2l n−2l+1 n−1 i=1 i i=1 i=1 X Y Y Y l,n where the summation is over all a1,...,an−1 ∈ Z>0 such that a1 +···+an−2l = an−2l+1 + ··· + an−1 + l. Pl n−1 The power series T (u) converges in the polydisc {(u2,...,un−1) ∈ C | |ui| < 1, i = 2,...,n − 1}. Proof. We prove formula (9.16) for j = 1, the proof for other values of j is similar. l The function T1(u) equals nπi/2 − e 3 −3 −1 −1 2 2 2 2 u1 (wu2 ··· un−2l − 1) (u1u2) (wu3 ··· un−2l − 1) ··· 2πi |w|=1/2 Z −1 −1 −1 −1 ··· (u1 ··· un−2l) 2 (w − 1) 2 (wu1 ··· un−2l) 2 (1 − un−2l+1/w) 2 ··· −1 −1 −1 ··· (wu1 ··· un−2l) 2 (1 − un−2l+1 ··· un−1/w) 2 (wu1 ··· un−2l) 2 u1 ··· un−2l dw

l−1 n−2l (−1) −1 = u−1(u ··· u )−l+1 (u ··· u ) 2 2πi 1 1 n−2l 1 i i=1 Y −3 × (1 − wu2 ··· un−2l) 2 (1 − wu3 ··· un−2l) ··· (1 − w) |w|=1/2 Z −1 2 dw × (1 − u /w) ··· (1 − u ··· u /w) . n−2l+1 n−2l+1 n−1 wl Expanding the binomials we obtain formula (9.24) for j = 1. The convergence property is clear.

(γl) 9.6. Formal solutions over Qp and truncation. For l =1,...,g, the formal series Iˆ (u) is a formal solution of system (2.1), in which the change of variables z = z(u) is performed and which is considered over the ﬁeld Qp of p-adic numbers.

Lemma 9.5. The formal series Iˆ(γl)(u), l = 1,...,g, are linear independent over the ﬁeld Qp. Proof. The proof follows from the fact that the monomials ul, l = 1,...,n, are linear independent over Qp. 9.7. Example n =3. In this case we have g = 1, u = z − z , u = z2−z3 , 1 1 3 2 z1−z3 ∞ − 3 − 1 − 1 − 3 − 1 − 1 (9.26) Iˆ(γ1)(u ,u ) = u−3/2 2 2 , 2 2 , 2 2 ua 1 2 1 a a a +1 a a +1 a 2 a=0 X ∞ − 1 − 3 a +1 −1/2 = u−3/2 2 2 , 1, ua 1 a +1 a −1/2 − a −1/2 − a 2 a=0 X NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 23

[ps−1] (9.27) Iˆ (u1,u2)= ps−1 s s s s s s s− 2 p −3 p −1 p −1 p −3 p −1 p −1 ps−3 p 3 = (−1) 2 u 2 2 2 , 2 2 , 2 2 ua 1 a a a +1 a a +1 a 2 a=0 X ps−1 s s s− 2 p −1 p −3 s ps−3 p 3 a +1 (p − 1)/2 =(−1) 2 u 2 2 2 , 1, ua 1 a +1 a (ps − 1)/2 − a (ps − 1)/2 − a 2 a=0 X 10. p-Adic convergence

Consider the ﬁeld Qp with the standard p-adic norm |t|p, t ∈ Qp. In this section we consider δl l,s ∞ ˆ(γl) the polynomial solutions ((−1) T (u))s=1 and the formal solutions I (u) as functions n−1 on Qp . Recall that Zp ⊂ Qp denotes the ring of p-adic integers.

10.1. Teichmuller representatives. For t ∈ Zp there exists the unique solution ω(t) ∈ Zp of the equation ω(t)p = ω(t) that is congruent to t modulo p. The element ω(t) is called the ps Teichmuller representative. It also can be deﬁned by ω(t) = lims→∞ t . The Teichmuller character is the homomorphism × × × Fp → Zp , α 7→ ω(α), α ∈ Fp .

~~For α ∈ Fp, r> 0, deﬁne the disc~~

~~Dα,r = {t ∈ Zp | |t − ω(α)|p~~

The space Zp is the disjoint union of the discs Dα,1, α ∈ Fp. The function ω : Zp → Zp, t 7→ ω(t), is a locally constant function equal to ω(α) on the disc Dα,1. For a subset S ⊂ Zp and a function f : S → Zp deﬁne the norm

k f k= sup |f(t)|p . t∈S ps ∞ Lemma 10.1. For any α ∈ Fp the sequence of polynomial functions (x )s=1 uniformly converges on Dα,1 to the constant function ω(α). ps Proof. We use the “fundamental inequality” from [Ro, II.4.3]: if |t|p 6 1, then |(1 + t) − s p−1 1|p 6 |t|p · max(|t|p, 1/p) . Now let t ∈ Dα,1 . Then t = 1+ t1, |t1|p 6 1/p. We have ps+1 ps ps (p−1)ps ps s+1 |t − t |p = |t |p|t − 1|p 6 |(1 + t1) − 1|p 6 1/p . ps1+s2 ps1 ps1+s2 ps1+s2−1 For positive integers s1,s2 and t ∈ Dα,1 we have |t − t |p = |t − t + s +s −1 s +1 s s p 1 2 p 1 p 1 s1+1 p ∞ t + ··· + t − t |p 6 1/p . Hence the sequence (x )s=1 is a Cauchy sequence. ps s For t ∈ D0,1, we have |t |p 6 1/p . For t1, t2 ∈ Dα,1, α =6 0, we have t1/t2 =1+ t with ps ps ps s+1 |t|p 6 1/p and |t1 − t2 |p = |(1 + t) − 1|p 6 1/p . The lemma is proved. (ps−1)/2 ∞ For α ∈ Fp consider the sequence of polynomial functions (x )s=1 on Dα,1. This sequence uniformly converges to 0 on the disc D0,1. (p−1)/2 2 2 Let α = 1. Let β ∈ Fp be such that β = α. The function Dβ,1 → Dα,1, t 7→ t , 1/2 is an analytic diﬀeomorphism. The inverse function Dα,1 → Dβ,1 will be denoted by x . 1/2 1/2 There are two square roots ±x . The root x corresponds to the chosen β ∈ Fp and the 1/2 root −x corresponds to −β ∈ Fp. 24 ALEXANDER VARCHENKO

2 (ps−1)/2 ∞ We change the variable x, set x = y , and lift the sequence (x )s=1 to the sequence ps−1 ∞ (y )s=1 of polynomial functions on Dβ,1. ps−1 ∞ Lemma 10.2. The sequence of polynomial functions (y )s=1 uniformly converges on Dβ,1 (ps−1)/2 ∞ to the function ω(β)/y. In other words, the sequence of polynomial functions (x )s=1 −1/2 uniformly converges on Dα,1 to the function ω(β)x . Proof. The lemma follows from Lemma 10.1.

−1 Let α(p−1)/2 = −1. Then ω(α)(ps−1)/2 = (−1)1+p+···+ps = (−1)s and the sequence (ps−1)/2 ∞ (x )s=1 has no limit on Dα,1. 1 3 − 2 − 2 Z 10.2. Approximation of binomial coeﬃcients. It is known that a , a ∈ p for a ∈ Z> . 0 Lemma 10.3. Let l1 > 0, l2 > 0 be integers. Then there exists an integer s0 > 0, such that > ps−1 > > for any integer s s0 and any integer a with 2 − l1 l2 + a 0 we have s − 1 − l p −1 − l (10.1) 2 1 − 2 1 6 1/ps−d−a , where d = l + l − 1/2. l + a l + a 1 2 2 2 p

~~Proof. We have~~

− 1 − l 1 l2+a 2 1 = (2(l + k) − 1), l2+a 1 l2 + a (−2) (l2 + a)! k=1 s Y p −1 − l 1 l2+a 2 1 = (2(l + k) − 1 − ps) . l2+a 1 l2 + a (−2) (l2 + a)! kY=1 1 1 ps−1 The p-adic norm of − 2 −l1 is 6 1. The diﬀerence − 2 −l1 − 2 −l1 is the sum of products l2+a l2+a l2+a 1 l2+a l +a (2(l1 + k) − 1) in each of which at least one of the factors 2(l1 + k) − 1 is (−2) 2 (l2+a)! k=1 replaced with −ps. We prove that even one such replacement implies that this summand of Q the diﬀerence has p-adic norm 6 1/ps+1/2−l1−l2−a. c 6 6 ps−1 Indeed, let 2(l1+k)−1= bp , p 6|b. We have 2(l1+k)−1 2(l1+l2+a)−1 2(l1+ 2 )−1. c s Hence bp 6 2l1 − 2+ p . Hence for any s, large enough, we have c 6 s, and the replacement s s−c of 2(l1 + k) − 1 with −p makes the norm of that summand 6 1/p . c 6 6 c−1 pc 6 bpc 6 1 We also have bp 2(l1 + l2 + a) − 1. Hence c p < 2 2 a + l1 + l2 − 2 . This shows that each summand has the p-adic norm 1/ 6 ps−c 6 1/ps+1/2−l1−l2−a. The lemma is proved. 10.3. Example n =3, continuation. Consider the formal power series ∞ − 3 − 1 − 1 − 3 − 1 − 1 (10.2) T 1(x) = 2 2 , 2 2 , 2 2 xa a a a +1 a a +1 a a=0 X in (9.26) and the sequence of polynomials

~~ps−1 2 ps−3 ps−1 ps−1 ps−3 ps−1 ps−1 (10.3) T 1,s(x)= 2 2 , 2 2 , 2 2 xa a a a +1 a a +1 a a=0 X NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 25~~

~~in (9.27) as functions on Zp.~~

~~1 Proposition 10.4. The power series T (x) uniformly convergence on D0,1. The sequence 1,s ∞ 1 of polynomial functions (T (x))s=1 uniformly converges on D0,1 to the function T (x).~~

1 3 − 2 − 2 Proof. The fact that the binomials a , a are p-adic integers implies the uniform con- 1 vergence of the power series T (x) on D0,1. ps−1 1 ∞ 1 a 1,s 2 1,s a 1 1,s 3 Let us write T (x) = a=0 Ta x and T (x) = a=0 Ta x , where Ta , Ta ∈ Zp. Then ps−1 1 1,s ∞ 1 a 2 1 1,s a T (x) − T (x)= psP+1 Ta x + a=0 (Ta − Ta P) x . Clearly the p-adic norm of the ﬁrst a= 2 ps−1 sum is 6 1/p 2 .P By Lemma 10.3 Pif s is big enough, then for each summand of the second 1 1,s a s−d sum and t ∈ D0,1 we have |(Ta − Ta ) t |p 6 1/p for some d independent of s and of the summand. This proves the proposition.

~~Consider the formal series~~

~~∞ − 3 − 1 − 1 − 3 − 1 − 1 (10.4) Iˆ(γ1)(u ,u ) = u−3/2 2 2 , 2 2 , 2 2 ua 1 2 1 a a a +1 a a +1 a 2 a=0 X and the sequence of polynomials~~

~~s− p 3 [ps−1] (10.5) (−1) 2 Iˆ (u1,u2)= ps−1 2 s s s s s s ps−3 p −3 p −1 p −1 p −3 p −1 p −1 = u 2 2 2 , 2 2 , 2 2 ua 1 a a a +1 a a +1 a 2 a=0 X ~~

2 1/2 as functions on Dα,1 × D0,1, where α = β for some β ∈ Fp. Then the function u1 : Dα,1 → (γ1) Dβ,1 is well-deﬁned and the series Iˆ (u1,u2) is a well-deﬁned function on Dα,1 × D0,1.

s p −3 s ∞ The sequence of polynomial functions 2 ˆ[p −1] uniformly Theorem 10.5. (−1) I (u1,u2) s=1 converges on to the function ˆ(γ1) . Dα,1 × D0,1 ω(β)I (u1,u 2) Proof. The theorem follows from Lemma 10.2 and Proposition 10.4 .

10.4. p-Adic convergence for arbitrary n. Given l, 1 6 l 6 g, consider the sequence of [lps−1] δl ˆ ˆ(γl) polynomials (−1) Ips (u) and the series I (u). Here u =(u1,...,un−1). l l We multiply the polynomials and the series by the same factor (u1 ··· un−2l) =(zn−2l −zn) [lps−1] l,s l δl ˆ and study the convergence of the sequence of polynomials J := (u1 ··· un−2l) (−1) Ips (u) 26 ALEXANDER VARCHENKO

~~l l (γl) to the series J := (u1 ··· un−2l) Iˆ (u). Introduce new variables :~~

n−2l i n−2l z − z x = u = (z − z ) , x = u ··· u = n−2l n , 1 j i n 2 2 n−2l z − z i=1 j=1 i=1 1 n Y Y Y zn−2l − zn zn−2l − zn x3 = u3 ··· un−2l = , ... xn−2l = un−2l = , z2 − zn zn−2l−1 − zn zn−2l+1 − zn zn−2l+2 − zn xn−2l+1 = un−2l+1 = , xn−2l+2 = un−2l+1un−2l+2 = , zn−2l − zn zn−2l − zn zn−1 − zn ... xn−1 = un−2l+1 ··· un−1 = . zn−2l − zn

~~Let x =(x1,...,xn−1). Then~~

1 1 l − 2 l − 2 l l J (x)= x1 Q (x2,...,xn−1)= x1 (Q1(x2,...,xn−1),...,Qn(x2,...,xn−1)), l with coordinates Qj deﬁned as follows. If j =1,...,n − 1, then (10.6) n−1 n−2l n−1 l,j −3 −1 Ql = x 2 2 xai−1 xai , j =1,...,n − 2l − 1, j j+1 a a i i j i i=2 X i=1Y,i6=j Y i=nY−2l+1 n−1 n−2l n−1 l,j −3 −1 Ql = 2 2 xai−1 xai , j = n − 2l,...,n − 1, j a a i i j i i=2 X i=1Y,i6=j Y i=nY−2l+1 l,j where the summation is over all a1,...,an−1 ∈ Z>0 such that a1 + ··· + an−2l = an−2l+1 +···+an−1 +l−1, if j 6 n−2l; and such that a1 +···+an−2l = an−2l+1 +···+an−1 +l, if n − 2l < j 6 n − 1; P

n−1 n−2l n−1 l,n −1 (10.7) Ql = 2 xai−1 xai , n a i i i=1 i i=2 X Y Y i=nY−2l+1 l,n where the summation is over all a1,...,an−1 ∈ Z>0 such that a1 + ··· + an−2l = an−2l+1 + ··· + an−1 + l. We also have P ps−1 ps−1 l,s 2 l,s 2 l,s l,s J (x)= x1 Q (x2,...,xn−2)= x1 (Q1 (x2,...,xn−2),...,Qn (x2,...,xn−2)), l,s with coordinates Qj (x) deﬁned as follows. If j =1,...,n − 1, then (10.8) s n−1 s n−2l n−1 l,j p −3 p −1 Ql,s = x 2 2 xai−1 xai , j =1,...,n − 2l − 1, j j+1 a a i i j i i=2 X i=1Y,i6=j Y i=nY−2l+1 s n−1 s n−2l n−1 l,j p −3 p −1 Ql,s = 2 2 xai−1 xai , j = n − 2l,...,n − 1, j a a i i j i i=2 X i=1Y,i6=j Y i=nY−2l+1 NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 27

l,s,j Z 6 6 ps−1 where the summation is over all a1,...,an−1 ∈ , 0 ai 2 , such that a1 + ··· + an−2l = an−2l+1 + ··· + an−1 + l − 1, if j 6 n − 2l; and such that a1 + ··· + an−2l = P an−2l+1 + ··· + an−1 + l, if n − 2l < j 6 n − 1; n−1 s n−2l n−1 l,n p −1 (10.9) Ql,s = 2 xai−1 xai , n a i i i=1 i i=2 X Y Y i=nY−2l+1 l,s,n 6 6 ps−1 where the summation is over all a1,...,an−1 ∈ Z, 0 ai 2 , such that a1 + ··· + an−2l = an−2l+1 + ··· + an−1 + l. P l n−2 Proposition 10.6. The power series Q (x2,...,xn−1) uniformly convergens on D0,1 . The l,s ∞ n−2 sequence of polynomial functions (Q (x2,...,xn−1))s=1 uniformly converges on D0,1 to the l function Q (x2,...,xn−1).

1 3 − 2 − 2 Proof. The fact that the binomials a , a are p-adic integers implies the uniform conver- l n−2 gence of the power series Q (x2,...,xn−1) on D0,1 . The proof of the uniform convergence l,s ∞ l of (Q (x2,...,xn−1))s=1 to Q (x2,...,xn−1) follows from Lemma 10.3 in the same way as the uniform convergence in the proof of Proposition 10.4. 1 l − 2 l Consider the formal series J (x) = x1 Q (x2,...,xn−1) and the sequence of polynomials ps−1 l,s 2 l,s n−2 2 J (x)= x1 Q (x2,...,xn−2) as functions on Dα,1 ×D0,1 , where α = β for some β ∈ Fp. 1/2 l Then the function x1 : Dα,1 → Dβ,1 is well-deﬁned and the series J (x) is a well-deﬁned n−2 function on Dα,1 × D0,1 . The sequence of polynomial functions l,s ∞ uniformly converges on Theorem 10.7. J (x) s=1 D × Dn−2 to the function ω(β)J l(x) . α,1 0,1 Proof. The theorem follows from Lemma 10.2 and Proposition 10.6.

~~Appendix A. The case n =3 and Dwork’s theory by Steven Sperber and Alexander Varchenko~~

In this appendix we consider only the special case n = 3 of previous considerations and show how this special case is related to Dwork’s theory in the classical paper [Dw]. A.1. Dwork on Legendre family. A.1.1. Consider the Legendre family of elliptic curves E(λ) defined by the affine equation y2 = x(x − 1)(x − λ). Let γ = γ(λ) be a family of 1-cycles on the curves E(λ) flat under the Gauss-Manin connection. Then the function dx (A.1) h(γ)(λ)= y Zγ satisfies the hypergeometric differential equation (A.2) λ(1 − λ)h′′ + (1 − 2λ)h′ − (1/4)h =0. 28 ALEXANDER VARCHENKO

~~All solutions of this equation are obtained in this way. One of the solutions h(δ1)(λ), for a suitable δ1, equals~~

~~1 1 ∞ −1/2 2 (A.3) F (λ)= F , ; 1; λ = λk , 2 1 2 2 k Xk=0 where 2F1 is the classical hypergeometric function.~~

~~A.1.2. If h(λ) is the elliptic integral in (A.1), then its derivative is 1 dx (A.4) h′ = . 2 (x − λ)y Zγ~~

~~A.1.3. Equation (A.2) can be written as a system of ﬁrst order linear diﬀerential equations for column 2-vectors I =(h, h′),~~

~~dI 0 1 (A.5) = B(λ)I, B(λ)= . dλ 1 2λ−1 4λ(1−λ) λ(1−λ) ~~

~~A.1.4. In [Dw] Dwork considers equation (A.2) over the ﬁeld Qp of p-adic numbers and studies analytic properties of the solution F (λ). We remind these properties below.~~

A.1.5. Let D ⊂ Qp be an open subset. A function f : D → Qp is called analytic at α ∈ D ∞ k if f can be presented by a power series k=0 ck(λ − α) with nonzero radius of convergence. A function f : D → Q is called analytic on D if f is the uniform limit on D of a sequence p P of rational function regular on D. In that case f is analytic at every point of D.

For i = 1,...,N let fi : Di → Qp be an analytic function on some domain Di. Assume that Di ∩ Di+1 =6 ∅ for i = 1,...,N − 1, and fi = fi+1 on Di ∩ Di+1, then we say that the N collection of functions fi deﬁnes an analytic element on ∪i=1Di. Cf. [Dw, Section 0].

~~A.1.6. Igusa noted in [Ig] that modulo p the polynomial~~

~~(p−1)/2 −1/2 2 (A.6) g(λ)= λj j j=0 X is the unique polynomial solution of equation (A.2) of degree less than p up to multiplication by a constant. Deﬁne~~

~~−1 D1 = {λ ∈ Zp | |g(λ)|p =1}, D2 = {λ | λ ∈ D1}, D = D1 ∪ D2.~~

~~Notice that D1, D2 are open and D1 ∩ D2 =6 ∅. More precisely,~~

~~D1 ∩ D2 = {λ ∈ Zp | |g(λ)|p =1, |λ|p =1}. NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 29~~

~~A.1.7. Dwork considers the functions F (λ) F ′(λ) (A.7) f(λ)= , η(λ)= , F (λp) F (λ)~~

~~deﬁned in a neighborhood of 0 ∈ D1 as ratios of the corresponding convergent power series expansions.~~

Dwork proves that f(λ) can be analytically continued to the domain D1. For that, he indicates a sequence of regular rational functions on D1, that sequence uniformly converges p on D1, and its limit equals F (λ)/F (λ ) in a neighborhood of 0, see [Dw, Lemma 3.4].

~~From that Dwork deduces that η(λ) has analytic continuation to the domain D1 in the same sense, see [Dw, Lemma 3.1].~~

~~Since η(λ) is analytic on D1, the function η(1/λ) is analytic on D2.~~

Using the properties of equation (A.2) Dwork shows that η(1 − λ) = −η(λ) on D1 and 2 shows that η(λ) = −η(1/λ)/λ − 1/(2λ) on D1 ∩ D2. Hence the function η(λ) on D1 and 2 the function −η(1/λ)/λ − 1/(2λ) on D2 deﬁne an analytic element on D. We will use the formulas

~~(A.8) η(1 − λ)= −η(λ), η(1/λ)= −λ2η(λ) − λ/2 ,~~

~~in Section A.5.~~

~~A.1.8. For α ∈ Qp let Vα be the space of germs at α of holomorphic solutions of equation (A.2). For α =06 , 1 we have dim Vα = 2 and for α =0, 1 we have dim Vα = 1.~~

~~For α ∈ D1 let Uα be the space of germs at α of analytic functions deﬁned by the equation du (A.9) = η(λ)u . dλ~~

By [Dw, Lemma 3.2], Uα is a subspace of Vα. We have Uα = Vα for α =0, 1. ′ For α ∈ D2 let Uα be the space of germs at α of analytic functions deﬁned by the equation du (A.10) =(−η(1/λ)/λ2 − 1/(2λ))u . dλ ′ ′ By [Dw, Lemma 3.2], Uα is a subspace of Vα. For α ∈ D1 ∩ D2 we have Uα = Uα.

~~A.1.9. By [Dw, Lemma 4.2], for α ∈ D the subspace Uα ⊂ Vα also can be characterized as the subspace of germs at α of holomorphic functions bounded in their disc of convergence.~~

More precisely, let Cp be the metric completion of the algebraic closure Q¯ p of the ﬁeld Qp. ∞ k Let α ∈ D. Let u(λ) = k=0 ck(λ − α) be an element of Vα. Consider u(λ) as a germ at α ∈ Qp ⊂ Cp of an analytic function on Cp. The germ u(λ) is called bounded on its disc of P convergence if u(λ) is bounded on its disc of convergence in Cp. ∞ k k Let r be the radius of convergence of u(λ)= k=0 ck(λ−α) . Deﬁne |u(λ)|0 = supk |ck|r . Then u(λ) is bounded in its disc of convergence if and only if |u(λ)|0 < ∞. P 30 ALEXANDER VARCHENKO

A.1.10. The function F (λ) is a holomorphic solution of equation (A.2) on the disc D0,1. A second solution of (A.5) is of the form G(λ)= F (λ) log λ + H(λ), where the function H(λ) is holomorphic on D0,1. Dwork speciﬁes H(λ) by [Dw, Equation (4.19)]. Then F G (A.11) F(λ)= F ′ G′ is a fundamental matrix of solutions of equation (A.5). In [Dw] Dwork introduces a 2 × 2-matrix function A(λ), and then proves the formula (−1)(p−1)/2 b (A.12) A(λ)= F(λ) M F(λp)−1, M = , 0 (−1)(p−1)/2p where b is a suitable number, see [Dw, Lemma 6.2] and formulas on page 72 in [Dw]. Dwork shows that A(λ) extends to an analytic function on the domain D3 ∪ D4, where

~~(A.13) D3 = {λ ∈ Zp | |λ|p =1, |λ − 1|p =1}, D4 = {λ ∈ Zp | ε< |λ|p < 1},~~

D3 ∪ D4 = {λ ∈ Zp | ε< |λ|p 6 1, |λ − 1|p =1}. Here ε is some explicit number, 0 <ε< 1. See the bottom of page 62 in [Dw] and the ﬁrst sentence of the proof of Theorem 6 in [Dw]. Formula (A.12) immediately implies that (A.14) A(λ)F(λp)= F(λ)M

on D4. The matrix A(λ) is called the matrix of the Frobenius transformation of solutions of equation (A.5) relative to the fundamental matrix F(λ). It follows from formula (A.14) that F (λp) F (λ) (A.15) A(λ) =(−1)(p−1)/2 F ′(λp) F ′(λ) on D4. This can be reformulated as the relation 1 1 (A.16) A(λ) =(−1)(p−1)/2f(λ) η(λp) η(λ) on D4. By the already formulated analytic properties of η(λ) and A(λ), relation (A.16) can be analytically continued to the domain D1 ∩ D3. × Equation (A.16) implies that for any α ∈ Fp −{1} such that ω(α) ∈ D1, the vector (1, η(ω(α))) is an eigenvector of the Frobenius matrix A(ω(α)) with eigenvalue (−1)(p−1)/2f(ω(α)), 1 1 (A.17) A(λ) =(−1)(p−1)/2f(ω(α)) . η(ω(α)) η(ω(α)) It is known that the zeta function of the elliptic curve deﬁned over Fp by the equation y2 = x(x − 1)(x − α) has two zeros, which are 1/((−1)(p−1)/2f(ω(α))), (−1)(p−1)/2f(ω(α))/p. (p−1)/2 It is also known that |f(ω(α))|p = 1. The number (−1) f(ω(α)) is called the unit root. See [Dw] and also (A.53). NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 31

A.2. KZ equations. The KZ equations (2.1) for n = 3 is the following system of diﬀerential and algebraic equations for a column 3-vector I = (I1,I2,I3) depending on variables z = (z1, z2, z3): ∂I 1 Ω Ω ∂I 1 Ω Ω (A.18) = 12 + 13 I, = 21 + 23 I, ∂z1 2 z1 − z2 z1 − z3 ∂z2 2 z2 − z1 z2 − z3 ∂I 1 Ω31 Ω32 = + I, 0= I1 + I2 + I3, ∂z3 2 z3 − z1 z3 − z2 where Ωij = Ωji and −1 1 0 −10 1 0 0 0 Ω12 = 1 −1 0 , Ω13 = 0 0 0 , Ω23 = 0 −1 1 .  0 00  1 0 −1 0 1 −1 We introduce new variables     z2 − z3 (A.19) u1 = z1 − z3, u2 = , u3 = z1 + z2 + z3, z1 − z3 see (9.7). Then system (A.18) takes the form ∂I 1 Ω + Ω + Ω ∂I 1 Ω Ω (A.20) = 12 13 23 I, = 12 + 23 I, ∂u1 2 u1 ∂u2 2 u2 − 1 u2 ∂I = 0, 0= I1 + I2 + I3. ∂u3 The variables in system (A.20) are separated, cf. (9.9). Denote W˜ = {(I1,I2,I3) | I1 + I2 + I3 =0}. Then

~~(A.21) (Ω12 + Ω13 + Ω23)|W˜ = −3 Id . Hence all solutions of system (A.20) have the form −3/2 (A.22) I = u1 (J1(u2),J2(u2),J3(u2)), J1 + J2 + J3 =0,~~

where the column vector J(u2) is a solution of the diﬀerential equation ∂J 1 Ω Ω (A.23) = 12 + 23 J. ∂u2 2 u2 − 1 u2 A.3. Solutions over C. Any solution of system (A.18) has the form 1 1 1 dx (A.24) I(γ)(z)= , , . γ x − z1 x − z2 x − z3 (x − z1)(x − z2)(x − z3) Z where γ is a ﬂat family of 1-cycles on the ellipticp curves of our family of curves. We change x and z in this integral by setting x = (z1 − z3)w + z3 and z = z(u) as in (A.19). Then integral (A.24) takes the form

(γ) −3/2 1 1 1 dw (A.25) I (u1,u2)= u1 , , . γ w − 1 w − u2 w (w − 1)(w − u2) w Z We take γ = γ1 to be the circle |w| = 1/2 orientedp counter-clockwise. We assume that u2 lies in this circle. We ﬁx the branch of (w − 1)(w − u2) w over the circle by choosing the argument of (w − 1)(w − u ) w at w =1/2, u =0 to be π/2. We multiply the circle 2 p 2 p 32 ALEXANDER VARCHENKO

1 with the chosen branch of the integrand by − 2π . This ﬁnishes the description of γ1. See the deﬁnition of cycles γl in Section 9.5. (γ1) We expand the integral I (u1,u2) as a power series in u2 and obtain ∞ − 3 − 1 − 1 − 3 − 1 − 1 (A.26) I(γ1)(u ,u ) = u−3/2 2 2 , 2 2 , 2 2 ua 1 2 1 a a a +1 a a +1 a 2 a=0 X ∞ − 1 − 3 a +1 −1/2 = u−3/2 2 2 , 1, ua , 1 a +1 a −1/2 − a −1/2 − a 2 a=0 X see Theorem (9.4). Denote a ∞ − 1 3 −1/2 (A.27) I := I(γ1) =(z − z )−3/2 2 − 2 a+1 , 1, z2−z3 1 3 a=0 a+1 a −1/2−a −1/2−a z1−z3

−3/2 ∞ 1 3 −1/2 P − 2 − 2 a+1 a (A.28) = u1 a=0 a+1 a −1/2−a , 1, −1/2−a u2 . This series is a solution of systemP (A.18 ). Remark. Formulas (A.25) and (A.26) imply that

−3/2 1 3 1 3 3 1 1 3 (A.29) I(u1,u2)= u1 2F1 2 , 2 ; 1; u2 , − 2 2F1 2 , 2 ; 2; u2 , − 2 2F1 2 , 2 ; 2; u2 , where 2F1(a, b; c; λ) is the classical hypergeometric function. A.4. Solutions as vectors of ﬁrst derivatives. Introduce the function dx (A.30) ℓ(γ)(z)= . γ (x − z )(x − z )(x − z ) Z 1 2 3 Then p ∂ℓ(γ) ∂ℓ(γ) ∂ℓ(γ) (A.31) I(γ)(z)= 2 , , . ∂z1 ∂z2 ∂z3 Changing the variable x = w(z1 − z3)+ z3 we write

~~(γ) −1/2 dw (A.32) ℓ (z) = (z1 − z3) , γ (w − 1)(w − z2−z3 )w Z z1−z3~~

−1/2 (γ)qz2 − z3 = (z1 − z3) h z1 − z3 −1/2 (γ) = u1 h (u2), where h(γ)(λ) is the elliptic integral in (A.1) and h(γ)(λ) is a solution of equation (A.2). Denote h(λ) := h(γ)(λ). Then (A.33) I(γ) =(z − z )−3/2 −h z2−z3 −2h′ z2−z3 z2−z3 , 2h′ z2−z3 , h z2−z3 +2h′ z2−z3 z2−z1 1 3 z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 −3/2 ′ ′ ′ (A.34) = u1 (−h(u2) − 2h (u2)u2, 2h (u2), h(u2)+2 h(u2)(u2 − 1)) . NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 33

~~Formula (A.34) relates solutions of system (A.5) and solutions of system (A.20). If (h, h′) is a solution of system (A.5), then equation~~

I1(u1,u2) −1 −2u2 −3/2 h(u2) (A.35) I2(u1,u2) = u1 0 2 ′     h (u2) I3(u1,u2) 1 2u2 − 2     gives a solution (I1,I2,I3) of system (A.20). Conversely, if (I1,I2,I3) is a solution of system (A.20) then formula

~~I1(u1,u2) h(u2) 3/2 −1 −u2 0 (A.36) ′ = u1 I2(u1,u2) h (u2) 0 1/2 0   I3(u1,u2)   gives a solution (h, h′) of system (A.5). Using the cycle γ1 we can evaluate~~

(γ1) −1/2 z2 − z3 (A.37) ℓ = (z1 − z3) F z1 − z3 1 1 = u−1/2F (u ), where F (λ)= F , ; 1; λ . 1 2 2 1 2 2 Denote I 2 ∂ℓ ∂ℓ ∂ℓ (A.38) ℓ := ℓ(γ1), I := = , , , ℓ ℓ ∂z1 ∂z2 ∂z3 where I is deﬁned in (A.27). Formulas (A.33) and (A.34) imply (A.39) I = 1 − 1 − 2η z2−z3 z2−z3 , 2η z2−z3 , 1+2η z2−z3 z2−z1 , z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 (A.40) = 1 − 1 − 2η(u )u , 2η(u ), 1+2η(u )(u − 1) , u1 2 2 2 2 2 where the function η(λ) is deﬁned in (A.7).

A.5. Six coordinate systems. System (A.18) of KZ equations has 6 distinguished coordinate systems (asymptotic zones). They are labeled by permutations σ =(i, j, k) ∈ S3. The σ σ σ σ coordinate system u =(u1 ,u2 ,u3 ) is deﬁned by the formulas

σ σ zj − zk σ (A.41) u1 = zi − zk, u2 = , u3 = z1 + z2 + z3. zi − zk For the identity element id = (1, 2, 3) the corresponding coordinate system is defined in (A.19). Having one of these coordinate systems we repeat the constructions of Sections A.3-A.4 and construct a scalar function ℓσ(z) and vector-valued functions Iσ(z), Iσ(z), such that Iσ(z) = Iσ(z)/ℓσ(z). For the identity element id = (1, 2, 3) these functions are ℓ(z), I(z), I(z) in (A.37), (A.27), (A.39). Notice that the functions ℓ(z) and I(z) are defined as integrals over γ1, and that γ1 is defined with the help of coordinates u1,u2,u3. For any σ the function Iσ(z) is a power series solution of system (A.18) in the chart with coordinates uσ, see (A.26) and (A.27). 34 ALEXANDER VARCHENKO

Below we list the functions Iσ : (A.42) I123 = 1 − 1 − 2η z2−z3 z2−z3 , 2η z2−z3 , 1+ η z2−z3 z2−z1 , z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 z1−z3 I321 = 1 1+2η z2−z1 z2−z3 , 2η z2−z1 , −1 − 2η z2−z1 z2−z1 , z3−z1 z3−z1 z3−z1 z3−z1 z3−z1 z3−z1 I213 = 1 2η z1−z3 , −1 − 2η z1−z3 z1−z3 , 1+2η z1−z3 z1−z2 , z2−z3 z2−z3 z2−z3 z2−z3 z2−z3 z2−z3 I132 = 1 − 1 − η z3−z2 z3−z2 , 1+2η z3−z2 z3−z1 , 2η z3−z2 , z1−z2 z1−z2 z1−z2 z1−z2 z1−z2 z1−z2 I231 = 1 1+2η z3−z1 z3−z2 , −1 − 2η z3−z1 z3−z1 , 2η z3−z1 , z2−z1 z2−z1 z2−z1 z2−z1 z2−z1 z2−z1 312 1 z1−z2 z1−z2 z1−z3 z1−z2 z1−z2 I = 2η , 1+2η , −1 − 2η 2 . z3−z2 z3−z2 z3−z2 z3−z2 z3−z2 z3−z2) ijk kji jik kji Theorem A.1. For (i, j, k) ∈ S3 consider the three functions I , I , I . Then I is transformed to Iijk by application of formula η(1 − λ) = −η(λ) and Ijik is transformed to Iijk by application of formula η(1/λ)= −λ2η(λ) − λ/2. Proof. The proof is straightforward. For example we check the statement for (i, j, k) = (1, 2, 3). In this case the functions I123, I321, I132 are 1 − 1 − 2η(u )u , 2η(u ), 1+2η(u )(u − 1) , u1 2 2 2 2 2 1 − 1+2η(1 − u )u , −2η(1 − u ), 1+2η(1 − u )(1 − u ) , u1 2 2 2 2 2 1 2η 1 , −1 − 2η 1 1 , 1+2η 1 1−u2 , u1u2 u2 u2 u2 u2 u2 where u1,u2 are defined in (A.19 ). Then formula η(1 − u2)= −η(u2) transforms the second 2 function to the first and the formula η(1/u2)= −u2η(u2)−u2/2 transforms the third function to the first.

~~Deﬁne ˜ 3 D0 = {(z1, z2, z3) ∈ Qp | zi =6 zj ∀i =6 j}.~~

~~For any σ =(i, j, k) ∈ S3 deﬁne~~

˜ σ ˜ zj − zk zj − zk D1 = (z1, z2, z3) ∈ D0 ∈ Zp, g =1 , zi − zk zi − zk p n o ˜ σ ˜ zi − zk ˜ σ ˜ σ ˜ σ ˜ σ ˜ ˜ σ D2 = (z1, z2, z3) ∈ D0 ∈ D1 , D = D1 ∪ D 2 , D = D , zj − zk n o σ∈S3 X where the function g is deﬁned in (A.6). ijk kji jik σ For any any (i, j, k) ∈ S3 the functions I , I , I deﬁne an analytic element on D˜ , see Section A.1.7 and [Dw]. Theorem A.1 implies the following corollary.

~~ijk ˜ Corollary A.2. The functions (I )(i,j,k)∈S3 deﬁne an analytic element on D.~~

Remark. Dwork’s formulas (A.8) present the S3-symmetries of the analytic element (η(λ), 2 −η(1 − λ), −η(1/λ)/λ − 1/(2λ)). Dwork’s S3-symmetries reformulated as S3-symmetries of ijk the analytic element (I )(i,j,k)∈S3 look even more well-rounded. NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 35

˜ 3 A.6. Subbundle. Denote W = {(I1,I2,I3) ∈ Qp | I1 + I2 + I3 =0}. System (A.18) of KZ equations defines a flat connection on the trivial bundle W˜ × D˜ 0 → D˜ 0. The flat sections of that bundle are solutions of system (A.18) of KZ equations. For any α ∈ D˜ such that α ∈ D˜ σ the vector Iσ(α) spans a one-dimensional subspace σ U˜α ⊂ W˜ . That subspace does not depend on σ such that α ∈ D˜ . The union of these subspaces defines a one-dimensional subbundle U˜ → D˜ of the trivial bundle W˜ × D˜ → D˜ . Theorem A.3. The subbundle U˜ → D˜ is invariant with respect to the KZ connection on W˜ × D˜ → D˜ .

σ Proof. For any σ ∈ S3 the subbundle U˜ → D˜ is generated by the ﬂat section I near the σ ˜ ˜ points where u2 = 0. Hence the subbundle U → D is generated by a ﬂat section near any point of D˜ , see Section A.1.8 and [Dw, Lemma 3.1].

σ Remark. For any σ ∈ S3 the flat section I generates the subbundle U˜ → D˜ near the σ σ points where u2 = 0. The power series I considered over C is the expansion of an integral σ over a cycle vanishing at the points where u2 = 0. The analytic continuation over C of that integral over that vanishing cycle could not generate a one-dimensional subbundle of the trivial bundle W˜ × D˜ → D˜ since the monodromy representation of the complex KZ equations in this case is irreducible. In contrast with this fact over C, the p-adic power σ series solutions I , σ ∈ S3, defined at different points glue together into a single line bundle U˜ → D˜ . This line bundle is what Dwork calls a p-adic cycle. This p-adic phenomenon was stressed by Dwork in [Dw] who titled his paper P -adic Cycles. Remark. The invariant subbundles of the KZ connection over C usually are related to some additional conformal block constructions, see [FSV1, FSV2, SV2, V7]. Apparently the subbundle U˜ → D˜ is of a different p-adic nature, cf. [V7].

σ A.7. Boundedness. Let σ ∈ S3 and α ∈ D˜ . For w ∈ W let I(z; w) be the germ at α of the solution of the KZ equations with initial condition I(α,w)= w. By formula (A.34), the coordinates of I(z; w) have the form

σ −3/2 σ ′ σ σ σ −3/2 ′ σ σ −3/2 σ ′ σ σ (u1 ) (−h(u2 ) − 2h (u2 )u2 ), (u1 ) 2h (u2 ), (u1 ) (h(u2 )+2h (u2 )(u2 − 1)), σ σ where h is the germ at the point u2 = u2 (α) of a solution of equation (A.2). We say that σ ′ σ σ ′ σ σ the germ I(z; w) is bounded if each of the germs −h(u2 ) − 2h (u2 )u2 , 2h (u2 ), h(u2 ) + ′ σ σ 2h (u2 )(u2 − 1) is bounded in its disc of convergence.

~~Theorem A.4. The germ I(z; w) is bounded if and only if w ∈ U˜α.~~

˜ σ Proof. Let w ∈ Uα. Then the germ h belongs to the corresponding subspace Uu2 (α) deﬁned in Section A.1.8. By [Dw, Lemma 4.2] the germ h is bounded in its disc of convergence, σ ′ σ σ ′ σ σ see Section A.1.9. Hence each of the three germs −h(u2 ) − 2h (u2 )u2 , 2h (u2 ), h(u2 )+ ′ σ σ 2h (u2 )(u2 − 1) is bounded in its disc of convergence. ˜ σ If w 6∈ Uα , then h 6∈ Uu2 (α). By [Dw, Lemma 4.2] the germ h is unbounded in its σ ′ σ σ ′ σ disc of convergence. Then at least one of the three germs −h(u2 ) − 2h (u2 )u2 , 2h (u2 ), σ ′ σ σ h(u2 )+2h (u2 )(u2 − 1) is unbounded in its disc of convergence. 36 ALEXANDER VARCHENKO

~~A.8. More domains. Denote (A.43) D˜ = (z , z , z ) ∈ D˜ z2−z3 ∈ Z , z2−z3 =1, z2−z3 − 1 =1 , 3 1 2 3 0 z1−z3 p z1−z3 p z1−z3 p~~

~~z2−z3 z2−z3 D˜ 4 = (z1, z2, z3) ∈ D˜ 0 ∈ Zp, ε< < 1 , z1−z3 z1 −z3 p where ε is the same number as in (A.13).~~

A.9. Frobenius map on solutions of KZ equations. Formula (A.14) describes the Frobenius map on solutions of equation (A.5). Solutions of equation (A.5) are identified with solutions of the KZ system (A.20) by formulas (A.35) and (A.36). That allows us to define the Frobenius map on solutions of the KZ system (A.20). Denote −1 −2u −1 −u 0 2 (A.44) B(u ,u )= u3/2 2 , C(u ,u )= u−3/2 0 2 . 1 2 1 0 1/2 0 1 2 1   1 2u2 − 2  1 0 0  1 0 we have B(u ,u )C(u ,u ) = , C(u ,u )B(u ,u ) = 0 1 0 . The second 1 2 1 2 0 1 1 2 1 2   −1 −1 0 ˜ 3 matrix defines the identity operator on the space W = {(I1,I2,I3) ∈ Qp | I1+ I2 + I3 =0}.

~~Recall the matrix F(u2) deﬁned in (A.11). By formula (A.35) the matrix~~

F˜(u1,u2)= C(u1,u2)F(u2) is a fundamental matrix of solutions of system (A.20). Recall the matrices A(λ), M in (A.12). Denote p p (A.45) A˜(u1,u2)= C(u1,u2) A(u2) B((u1) , (u2) ). ˜ 3 This is a 3 × 3 matrix valued function, whose values preserve the subspace W ⊂ Qp. Theorem A.5. We have p p (A.46) A˜(u1,u2) F˜((u1) , (u2) ) = F˜(u1,u2) M.

~~The matrix A˜(u1,u2) extends to an analytic function on the domain D˜ 3 ∪ D˜ 4. Proof. The theorem is a corollary of formula (A.14) and Dwork’s statements listed in Section A.1.10.~~

We call A˜(u1,u2) the matrix of the Frobenius transformation of solutions of system(A.20) relative to the fundamental matrix F˜(u1,u2) on the domain D˜ 3 ∪ D˜ 4. Recall the distinguished solution −3/2 ′ ′ ′ (A.47) I(u1,u2)= u1 (−F (u2) − 2F (u2)u2, 2F (u2), F (u2)+2F (u2)(u2 − 1))

of system (A.20) deﬁned near the points where u2 =0, see (A.34). By (A.15) we have p p (p−1)/2 (A.48) A˜(u1,u2) I((u1) , (u2) )=(−1) I(u1,u2) ˜ −1/2 on D4. Recall ℓ(u1,u2) = u1 F (u2) in (A.37). Dividing both sides in (A.47) by p p ℓ((u1) , (u2) ) we can reformulate (A.48) as NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 37

~~˜ p p (p−1)/2 (p−1)/2 (A.49) A(u1,u2) I((u1) , (u2) )=(−1) u1 f(u2) I(u1,u2)~~

on D˜ 4, see I(u1,u2) in (A.40) and f(u2) in (A.7). As in Section A.1.10 we conclude with ˜ (1,2,3) ˜ Dwork that relation (A.49) can be analytically continued to the domain D1 ∩ D3. × × Equation (A.49) implies that for any α ∈ Fp −{1}, β ∈ Fp such that ω(α) ∈ D1, the vector I(ω(β),ω(α)) is an eigenvector of the Frobenius matrix A˜(ω(β),ω(α)) with eigenvalue ω(β(p−1)/2)(−1)(p−1)/2f(ω(α)), (A.50) A˜(ω(β),ω(α)) I(ω(β),ω(α)) = ω(β(p−1)/2)(−1)(p−1)/2f(ω(α)) I(ω(β),ω(α)) .

In this Section A.9 we described the matrix A˜(u1,u2) of the Frobenius transformation of solutions of system (A.18) written in coordinates u1,u2,u3 corresponding to the chart labeled by the identify permutation (1, 2, 3) ∈ S3, In the same way we may start with the chart corresponding to any permutation σ ∈ S3 and describe the matrix of the Frobenius σ σ σ transformation of solutions of system (A.18) written in coordinates u1 ,u2 ,u3 . A.10. Eigenvalue ω(β(p−1)/2)(−1)(p−1)/2f(ω(α)). Theorem A.6. The number ω(β(p−1)/2)(−1)(p−1)/2f(ω(α)) is the unit root of the elliptic curve E(α, β) deﬁned over Fp by the aﬃne equation (A.51) w2 = β v(v − 1)(v − α).

× 2 Proof. Assume that β ∈ Fp is a square, β = γ for some γ ∈ Fp. Then on the one hand the change of the variablew ˜ = w/γ makes E(α, β) isomorphic to E(α, 1). On the other hand β(p−1)/2 = 1 and ω(β(p−1)/2)(−1)(p−1)/2f(ω(α))=(−1)(p−1)/2f(ω(α)), where the last number is the unit root of the elliptic curve E(α, 1) by [Dw]. × Assume that β ∈ Fp is not a square. Denote by N1,β the number of points on E(α, β). Then

(A.52) N1,1 + N1,β =4+4+2(p − 3)=2p +2. Indeed the number 4 + 4 corresponds to the points (0, 0), (0, 1), (0,α), ∞ on E(α, β) and on E(α, 1). The number 2(p − 3) corresponds to p − 3 elements of Fp −{0, 1,α}. Namely if v0 ∈ Fp −{0, 1,α}, then exactly one of the two elements βv0(v0 −1)(v0 −α), v0(v0 −1)(v0 −α) is a square in Fp and exactly one of the two elliptic curves has two points over v = v0, while the other curve does not have points over v0. It is known that the zeta function of the curve E(α, β) has the form ∞ N (1 − R T )(1 − (p/R )T ) (A.53) exp s,β T s = β β . s (1 − T )(1 − pT ) s=1 X Here Ns,β is the number points on E(α, β) considered over the ﬁeld Fps , while the number Rβ has |Rβ|p = 1 and is called the unit root, for example see [Mo]. Equation (A.53) implies s s s that for any s we have Ns,β =1+ p − Rβ − (p/Rβ) . In particular for s = 1 we have

~~(A.54) N1,β =1+ p − Rβ − p/Rβ. 38 ALEXANDER VARCHENKO~~

~~From (A.52) and (A.54) we obtain~~

~~0= Rβ + p/Rβ + R1 + p/R1 =(Rβ + R1)(1 + p/RβR1).~~

Since the second factor is nonzero we conclude that Rβ = −R1. By [Dw] we have R1 = (p−1)/2 (p−1)/2 (p−1)/2 (p−1)/2 (−1) f(ω(α)). Hence Rβ = − (−1) f(ω(α)) = ω(β )(−1) f(ω(α)). The theorem is proved. The relation between the eigenvector I(ω(β),ω(α)) and the elliptic curve E(α, β), indi- cated in Theorem A.6, can be explained as follows. Over C the vector I is given by integrals 2 over cycles on elliptic curves with equation y =(x − z1)(x − z2)(x − z3). After the change of variables z2 − z3 x =(z1 − z3)w + z3, u1 = z1 − z3, u2 = , y =(z1 − z3)v, z1 − z3 the equation takes the form 2 (A.55) v = u1(w − 1)(w − u2)w.

The eigenvector I(ω(β),ω(α)) corresponds to the curve in (A.55) with (u1,u2)=(ω(β),ω(α)), and (ω(β),ω(α)) ≡ (β,α) mod p. It is more surprising that system (A.18) of KZ equations gives a bit more arithmetic information than the hypergeometric equation (A.5), despite the fact system (A.18) and equation (A.5) are equivalent by (A.35) and (A.36). Indeed Dwork’s eigenvectors in (A.16) give unit roots of elliptic curves E(α, 1) while the eigenvectors in (A.50) coming from the KZ equations give unit roots of more general elliptic curves E(α, β). ijk A.11. Approximation of analytic element (I )(i,j,k)∈S3 by rational functions. Let s be a positive integer.

~~s 3 M A.11.1. Let M =(p − 1)/2, Φps (x, z)= i=1(x − zi) ,~~

Φps (x, z) Φps (x, z) Φps (x, z) i i P s (x, z)= , Q , = P s (z) x . p x − z x − z x − z p 1 2 3 i X [ps−1] ps−1 Denote Ips (z) = Pps (z). The functions −1 2 [ps−1] [ps −1] s−2 [p −1] s−1 [p−1] Ips (z), pIps−1 (z), ..., p Ip2 (z), p Ip (z), are solutions of system (A.18) modulo ps by Theorem 4.4. i i Let Φps (x, z)= i Φps (z)x . Denote ps−1 (A.56) P ℓps (z) = Φps (z). A.11.2. For k =1, 2, 3 let i i Pk,ps(v, z) := Pps (v + zk, z)= Pk,ps (z) v . i X [ps−1] ps−1 Denote Ik,ps (z) = Pk,ps (z). The functions −1 2 [ps−1] [ps −1] s−2 [p −1] s−1 [p−1] Ik,ps (z), pIk,ps−1 (z), ..., p Ik,p2 (z), p Ik,p (z), are solutions of system (A.18) modulo ps by Theorem 9.1. NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 39

i i Let Φk,ps (v, z) := Φps (v + zk, z)= i Φk,ps (z)v . Denote [ps−1] ps−1 (A.57) ℓk,pPs (z) = Φk,ps (z). A.11.3. Recall the homomorphisms Z → Z/psZ, Z[z] → (Z/psZ)[z], Z[z]3 → (Z/psZ)[z]3 r denoted by πs. Recall the subring Z[z]pr ⊂ Z[z] of quasi-constants modulo p . A.11.4. Define the filtration 0 1 s−1 s (A.58) 0= Mps ⊂Mps ⊂···⊂Mps ⊂Mps = Mps , r t t s−r [p −1] where Mps = πs r=1 cr(z) p Ipr (z) | cr(z) ∈ Z[z]pr , t = 1,...,s. Every element s of Mps is a polynomial solution of system (A.18) with coefficients in Z/p Z. P Define the filtration 0 1 s−1 s (A.59) 0= Mk,ps ⊂Mk,ps ⊂···⊂Mk,ps ⊂Mk,ps = Mk,ps , r t t s−r [p −1] where Mk,ps = πs r=1 cr(z) p Ik,pr (z) | cr(z) ∈ Z[z]pr , t =1,...,s. Every element s of M s is a polynomial solution of system (A.18) with coefficients in Z/p Z by Theorem 9.1. p P t t By Theorem 9.1 filtrations (A.58) and (A.59) coincide, Mps = Mk,ps for any k, t. A.11.5. Define the filtration 0 1 s−1 s (A.60) 0= Lps ⊂Lps ⊂···⊂Lps ⊂Lps = Lps , r t t s−r [p −1] where Lps = πs r=1 cr(z) p ℓpr (z) | cr(z) ∈ Z[z]pr , t =1,...,s. Define the filtration P 0 1 s−1 s (A.61) 0= Lk,ps ⊂Lk,ps ⊂···⊂Lk,ps ⊂Lk,ps = Lk,ps , r t t s−r p −1] where Lk,ps = πs r=1 cr(z) p ℓk,pr (z) | cr(z) ∈ Z[z]pr , t =1,...,s. t t It is easy to see that P filtrations (A.60) and (A.61) coincide, Lps = Lk,ps for any k, t.

A.11.6. Let u =(u1,u2,u3) be the coordinates in the chart corresponding to (1, 2, 3) ∈ S3, [ps−1] [ps−1] see (A.19). Consider the functions I3,ps (z) ∈Mps , ℓ3,ps (z) ∈Lps . Denote ˆ[ps−1] [ps−1] ˆ[ps−1] [ps−1] (A.62) I (u) := I3,ps (z(u)), ℓ (u) := ℓ3,ps (z(u)). We have s ps−3 ps−1 s− s− s− s− s− s− p −3 s p 3 p 1 p 1 p 3 p 1 p 1 2 ˆ[p −1] 2 2 2 2 2 2 2 2 a (A.63) (−1) I = u1 a=0 a a , a+1 a , a+1 a u2,

s ps−1 ps−1 s− 2 p −1 s p 1 2 Pˆ[p −1] 2 2 2 a (A.64) (−1) ℓ = u1 a=0 a u2 , see formula (A.63) and (10.5). P s− p 1 ps−1 2 2 2 a s Notice that a=0 a u2 is a solution of eqution (A.2) modulo p . s p −3 s ∞ 2 ˆ[p −1] As s →∞ theP sequence (−1) I s=1 of vector-valued polynomials uniformly converges to the series I in (A.28) near the points where u2 = 0, see Theorem 10.5. Similarly as s p −3 s ∞ 2 ˆ[p −1] s → ∞ the sequence (−1) ℓ s=1 of scalar polynomials uniformly converges to the series ℓ in (A.37) near the points where u = 0. 2 40 ALEXANDER VARCHENKO

As the sequence ˆ[ps−1] ˆ[ps−1] ∞ of vector-valued rational Corollary A.7. s → ∞ − I ℓ s=1 functions uniformly converges to the function I in (A.40) near the points where u =0. 2 σ A.11.7. Let σ =(i, j, k) ∈ S3. Similarly to Section A.11.6 consider the coordinates u and [ps−1] [ps−1] the functions Ik,ps (z) ∈Mps , ℓk,ps (z) ∈Lps . Denote ˆ[σ,ps−1] σ [ps−1] σ ˆ[σ,ps−1] σ [ps−1] σ (A.65) I (u ) := Ik,ps (z(u )), ℓ (u ) := Ik,ps (z(u )). Similarly to Section A.11.6 we obtain the following corollary. As the sequence ˆ[σ,ps−1] ˆ[σ,ps−1] ∞ of vector-valued rational Corollary A.8. s → ∞ − I ℓ s=1 functions uniformly converges to the function Iσ in (A.40) near the points where uσ =0. 2 A.11.8. This Appendix A is devoted to the relation between the analytic element ijk (I )(i,j,k)∈S3 and Dwork’s theory in [Dw].

As additional information, for any σ =(i, j, k) ∈ S3, Corollary A.8 indicates the sequence [ps−1] [ps−1] of polynomials Ik,ps (z) ∈Mps , ℓk,ps (z) ∈Lps whose ratio p-adically tends to the function ijk σ ijk I near the points u2 = 0 where the function I is initially defined. A.12. Further directions. In Sections 2 - 10 we considered system (2.1) of KZ equations with parameter n =2g + 1 and constructed polynomial solutions of system (2.1) modulo ps. We defined the module Mps of the constructed solutions and studied the limit of Mps as s →∞. Namely we considered a special coordinate system u = u(z) in (9.7) associated with one of the asymptotic zones of the KZ equations and showed that in this coordinate system the limit of Mps as s →∞ produces a g-dimensional space of solutions of system (2.1) over p-adic numbers Qp in the neighborhood of the point u = 0. Constructions in this appendix for g = 1 and Dwork’s theory in [Dw] suggest the following project. Consider all asymptotic zones of system (2.1), see their definition for example in [V2]. The asymptotic zones are labeled by suitable trees T . These trees are analogs of the elements σ ∈ S3 in the appendix. Each asymptotic zone has a distinguished system of T coordinates u . Probably, for every asymptotic zone the limit of Mps as s → ∞ produces a g-dimensional space VT of solutions of system (2.1) considered over Qp in a neighborhood T of the point u = 0. Probably the spaces VT of p-adic solutions, defined at different places uT = 0, analytically continue into a single global invariant g-dimensional vector subbundle of the associated KZ connection on the trivial vector bundle of rank 2g. Following Dwork and Theorem A.4 we may expect that this subbundle is spanned at any point of the base by the germs of all solutions of the KZ equations bounded in their polydiscs of convergence. This subbundle would give a generalization of the line subbundle generated by the analytic ijk element (I )(i,j,k)∈S3 constructed in this appendix. Probably, that g-dimensional subbundle 2 n will determine the set of unit roots of the curves with equation y = i=1(x − zi) over the field Fp similarly to how it is done in Sections A.9 and A.10 for the elliptic curves. Q References [AH] J. Achter, E. Howe, Hasse-Witt and Cartier-Manin matrices: A warning and a re- quest, arXiv:1710.10726, 1–14 [BV1] F. Beukers, M. Vlasenko, Dwork Crystals I, arXiv:1903.11155, 1–27 NOTES ON SOLUTIONS OF KZ EQUATIONS MODULO ps 41

[BV2] F. Beukers, M. Vlasenko, Dwork Crystals II, arXiv:1907.10390, 1–12 [Cl] H.C. Clemens, A scrapbook of complex curve theory, Second edition, Graduate Stud- ies in Mathematics, 55, AMS, Providence, RI, 2003. xii+188 pp [Dw] B. Dwork, P -adic Cycles, Publications Mathematiques de lHES, 37 (1969), 27–115 [EFK] P. Etingof, I. Frenkel, A.Kirillov, Lectures on representation theory and Knizhnik- Zamolodchikov equations, Mathematical Surveys and Monographs, 58, AMS, Prov- idence, RI, 1998.l xiv+198 pp. ISBN: 0-8218-0496-0 [FSV1] B. Feigin, V. Schechtman, and A. Varchenko, On algebraic equations satisfied by hypergeometric correlators in WZW models, I, Comm. Math. Phys. 163 (1994), 173– 184 [FSV2] B. Feigin, V. Schechtman, and A. Varchenko, On algebraic equations satisfied by hypergeometric correlators in WZW models, II, Comm. In Math. Phys. 70 (1995), 219–247 [Ig] J. Igusa, Class number of a definite quaternion with prim discriminant, Proc. Natl. Acad. Sci. U S A, 44(4) (1958), 312–314 [KZ] V. Knizhnik and A. Zamolodchikov, Current algebra and the Wess-Zumino model in two dimensions, Nucl. Phys. B247 (1984), 83–103 [Lu] E. Lucas, Theorie des Fonctions Numeriques Simplement Periodiques, American Journal of Mathematics. 1 (2) (1878) 184–196, doi:10.2307/2369308, JSTOR 2369308, MR 1505161 [Ma] Y.I. Manin, The Hasse-Witt matrix of an algebraic curve, (Russian) Izv. Akad. Nauk SSSR Ser. Mat. 25 (1961), 153–172 [Mo] P. Monsky, p-adic analysis and zeta functions, Lectures in Mathematics, Kyoto Uni- versity, Kinokuniya Book-Store Co., 1970, 123 pp [RV1] R. Rimányi, A. Varchenko, The Fp-Selberg integral, arXiv:2011.14248, 1–19 [RV2] R. Rimányi, A. Varchenko, The Fp-Selberg integral of type An, arXiv:2012.01391, 1–21 [Ro] Alain M. Robert, A Course in p-adic Analysis, (Graduate Texts in Mathematics), Springer, 2000 [SV1] V. Schechtman, A. Varchenko, Arrangements of Hyperplanes and Lie Algebra Ho- mology, Invent. Math. 106 (1991), 139–194 [SV2] V. Schechtman, A. Varchenko, Solutions of KZ differential equations modulo p, The Ramanujan Journal, 48 (3), 2019, 655–683, https://doi.org/10.1007/s11139-018-0068-x, arXiv:1707.02615 [SliV] A. Slinkin, A. Varchenko, Hypergeometric Integrals Modulo p and Hasse–Witt matrices, arXiv:2001.06869, 1–36 [V1] A. Varchenko, Beta-Function of Euler, Vandermonde Determinant, Legendre Equa- tion and Critical Values of Linear Functions of Configuration of Hyperplanes, I. Izv. Akademii Nauk USSR, Seriya Mat., 53:6 (1989), 1206–1235; II, Izv. Akademii Nauk USSR, Seriya Mat. 54:1 (1990), 146–158 [V2] A. Varchenko, Asymptotic solutions to the KZ equation and crystal base, Comm. in Math. Phys. 171 (1995) 99–137 [V3] A. Varchenko, Special functions, KZ type equations, and Representation theory, CBMS, Regional Conference Series in Math., n. 98, AMS (2003) 42 ALEXANDER VARCHENKO

[V4] A. Varchenko, Solutions modulo p of Gauss-Manin diﬀerential equations for multidimensional hypergeometric integrals and associated Bethe ansatz, arXiv:1709.06189, Mathematics 2017, 5(4), 52; doi:10.3390/math5040052, 1–18 [V5] A. Varchenko, Hyperelliptic integrals modulo p and Cartier-Manin matrices, arXiv:1806.03289, Pure and Applied Math Quarterly, 16 (2020), n. 3, 315–336 [V6] A. Varchenko, Remarks on the Gaudin model modulo p, arXiv:1708.06264, Journal of Singularities, 18 (2018), 486–499 [V7] A. Varchenko, An invariant subbundle of the KZ connection mod p and reducibility

~~of sl2 Verma modules mod p, arXiv:2002.05834 , 1–14 [V8] A. Varchenko, Determinant of Fp-hypergeometric solutions under ample reduction, arXiv:2010.11275c , 1–22~~