arxiv: v2 [math.nt] 10 Sep 2018

Transcription

1 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS arxiv: v [math.nt] 10 Sep 018 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB Abstract. It is well-known that the Apéry sequences which arise in the irrationality proofs for ζ( and ζ(3 satisfy many intriguing arithmetic properties and are related to the pth Fourier coefficients of modular forms. In this paper, we prove that the connection to modular forms persists for sequences associated to Brown s cellular integrals and state a general conjecture concerning supercongruences. 1. Introduction Period integrals on the moduli space M 0,N of curves of genus 0 with N marked points have featured prominently in a variety of mathematical and physical contexts. These period integrals are of particular importance as, for example, it is now known that they are Q-linear combinations of multiple zeta values [1] thereby proving a conjecture of Goncharov and Manin [], provide an effective computation of a large class of Feynman integrals [10], occur in superstring theory [11], [41] and have potential connections to higher Frobenius limits in relation to the Gamma conjecture [1]. Recently, Brown [13] introduced a program where such period integrals play a central role in understanding irrationality proofs of values of the Riemann zeta function. Before discussing this program, we first briefly review the classical situation and some subsequent developments. Inspired by Apéry, Beukers [7] gave another proof of the irrationality of ζ( and ζ(3 by considering the integrals and He showed that J n := 1 I n := ( 1 n x n (1 x n y n (1 y n (1 xy n+1 dxdy (1.1 x n (1 x n y n (1 y n w n (1 w n (1 (1 xyw n+1 dxdydw. (1. I n = a(nζ(+ã(n and J n = b(nζ(3+ b(n where ã(n and b(n are explicit rational numbers and n ( n ( n+k n ( n ( n+k a(n =, b(n = (1.3 k k k k k=0 k=0 Date: September 11, Mathematics Subject Classification. Primary: 11F03; Secondary: 11M3, 14H10. Key words and phrases. Modular forms, period integrals, supercongruences. 1

2 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB are the Apéry numbers. Since their introduction, there has been substantial interest in both the intrinsic arithmetic properties of the Apéry numbers and their relationship to modular forms. For example, consider η 6 (4z =: α(nq n, n=1 the unique newform in S 3 (Γ 0 (16,( 4, where η(z = q 1/4 n 1 (1 qn is the Dedekind etafunction, q = e πiz and z H. Ahlgren [1] showed that, for all primes p 5, ( a α(p (mod p, (1.4 thus confirming a conjecture in [4]. In [3], Ahlgren and Ono proved a conjecture of Beukers [9] which stated that, if p is an odd prime, then ( b β(p (mod p, (1.5 where η 4 (zη 4 (4z =: β(nq n is the unique newform in S 4 (Γ 0 (8. Both of these influential results required the existence of an underlying modular Calabi-Yau variety, specifically a K3-surface for (1.4 and a Calabi-Yau 3-fold for (1.5 in order to relate the pth Fourier coefficients to finite field hypergeometric series. Concerning arithmetic properties, Coster [17] proved the supercongruences and n=1 a(mp r a(mp r 1 (mod p 3r (1.6 b(mp r b(mp r 1 (mod p 3r (1.7 for primes p 5 and integers m, r 1. Other supercongruences of this type have been studied in numerous works (for example, see [5], [8], [15], [16], [0], [5], [31] [34], [43]. In order to view (1.4 (1.7 from a general perspective, we discuss the setup from [13]. Recall that M 0,N, N 4, is the moduli space of genus zero curves (Riemann spheres with N ordered marked points (z 1,...,z N. It is the set of N-tuples of distinct points (z 1,...,z N modulo the equivalence relation given by the action of PSL. This action is triply-transitive and so there is a unique representative of each equivalence class such that z 1 = 0, z N 1 = 1 and z N =. We introduce simplicial coordinates on M 0,N by setting t 1 = z, t = z 3,..., t N 3 = z N. This yields the identification M 0,N = {(t1,...,t N 3 (P 1 \{0,1, } N 3 t i t j for all i j}. A typical period integral on M 0,N can be given as S N t a i i (1 t j b j (t i t j c ij dt 1...dt N 3, (1.8 where a i, b j and c ij Z are such that (1.8 converges and the simplex S N = {(t 1,...,t N 3 R N 3 : 0 < t 1 <... < t N 3 < 1}

3 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 3 is a connected component of M 0,N (R. It was proved by Brown [1] that the integrals (1.8 are Q-linear combinations of multiple zeta values of weight up to and including N 3. That these typically involve multiple zeta values of all weights is an obstruction to irrationality proofs. For example, generic period integrals on M 0,6 yield linear forms in 1, ζ( and ζ(3. To ensure the vanishing of coefficients, we consider a variant of the classical dinner table problem [6], [38]. Suppose we have N 5 people sitting at a round table. We permute them such that each person has two new neighbors. We represent the new seating arrangement (non-uniquely by a permutation σ = σ N on {1,,...,N} and write σ = (σ(1,...,σ(n. A permutation σ is called convergent if no set of k elements in {1,...,N} are simultaneously consecutive for {1,...,N} and σ for all k N. If N 7, then this is equivalent to the dinner table problem. For N 8, this is a new condition. The main idea of [13] is to associate a rational function f σ and a differential (N 3-form ω σ to a given σ as follows. Formally, we define f σ = i z i z i+1 z σ(i z σ(i+1 and ω σ = dt 1...dt N 3 i z σ(i z σ(i+1, (1.9 where the product is over all indices i modulo N, z 1 = 0 and z N 1 = 1. We remove factors with z N and then let z i+1 t i for i = 1,,...,N 3. If we consider the cellular integral I σ (n = I N,σ (n := fσ n ω σ, (1.10 S N then I σ (n converges if and only if σ is convergent. For n = 0, we then obtain the cell-zeta values ζ σ (N 3 = I σ (0 studied in [14], which are multiple zeta values of weight N 3. More generally, by [13, Corollary 8.], I σ (n is a Q-linear combination of multiple zeta values of weight less than or equal to N 3. Suppose that this linear combination is of the form A σn (nζ σ (N 3, with A σn (n Q, plus a combination of multiple zeta values of weight less than N 3 (note that A σn (n is necessarily unique if multiple zeta values of different weight are linearly independent over Q; alternatively, since this independence remains open, a unique value can be selected, unconditionally, by working motivically [13]. We then say that A σ (n = A σn (n is the leading coefficient of the cellular integral I σ (n. By construction, A σ (0 = 1. Conjecturally, we can define leading coefficients in this fashion for all convergent σ. Since we will be concerned with specific permutations σ, in which case the linear combinations of multiple zeta values can be made explicit, this issue will not disturb our discussion. For example, if N = 5, then σ 5 = (1,3,5,,4 is the unique (up to dihedral symmetry, see Section convergent permutation, I σ5 (n recovers Beukers integral (1.1 after a change of variables, and the leading coefficients A σ5 (n are the Apéry numbers a(n. Similarly, for N = 6, one obtains (1. upon considering σ 6 = (1,5,3,6,,4 and the leading coefficients A σ6 (n are the Apéry numbers b(n. This general framework raises some natural questions. Is there an analogue of (1.4 and (1.5 for N 7? Do the leading coefficients A σ (n satisfy supercongruences akin to (1.6? The first main result generalizes (1.4 to all odd weights greater than or equal to 3. Theorem 1.1. For each odd positive integer N 5, there exists a convergent σ N and a modular form f k (z =: n=1 γ k(nq n of weight k = N such that, for all primes p 5, ( A σn γ k (p (mod p. (1.11

4 4 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB As the product of cellular integrals is not necessarily a cellular integral (see Example.4, powers of the Apéry numbers a(n are not automatically leading coefficients. We prove Theorem 1.1 by carefully constructingan explicit family of convergent σ N suchthat A σn (n = a(n (N 3/ and the required modular forms f k (z. For example, if N = 5, then A σ5 (n = a(n and f 3 (z = η 6 (4z. Thus, we recover (1.4. Now, consider the convergent permutation σ 8 = (8,3,6,1,4,7,,5 and η 1 (z =: γ(nq n, the unique newform in S 6 (Γ 0 (4. Our second main result is a higher weight version of (1.5. Theorem 1.. Let p be an odd prime. Then ( A σ8 γ(p (mod p. (1.1 In Section 10.1 of [13], Brown provided the following table for the number C N of convergent σ N (up to symmetries: n=1 N C N Table 1. C N for 5 N 11 Using the techniques from Section 3, we have obtained explicit binomial sum expressions for all 898 of the leading coefficients A σn (n for 7 N 10. For example, if N = 7, then the convergent permutations and associated leading coefficients are given in Table. The ordinary generating functions of these five sequences satisfy fourth order differential equations of Calabi Yau type, and thus appear as sequences in the table [4] (where they are numbered as 193, 43, 101, 198, and 7. The remaining list of convergent permutations and binomial sum expressions for the leading coefficients are available upon request. Based on numerical evidence, we make the following conjecture which suggests that (1.6 is a generic property of leading coefficients. Conjecture 1.3. For each N 5 and convergent σ N, the leading coefficients A σn (n satisfy for all primes p 5 and integers m, r 1. A σn (mp r A σn (mp r 1 (mod p 3r (1.13 Given Theorems 1.1 and 1., it is natural to wonder if such supercongruences hold for all leading coefficients A σn (n. Do they arise from L-series attached to Galois representations? It appears that a result similar to (1.5 and (1.1 also holds for a convergent σ 10 and modular form of weight 8. This and other observations will be the subject of forthcoming work. Finally, it would be of interest to examine Theorems 1.1 and 1. and Conjecture 1.3 from the recent motivic perspective of [40]. The paper is organized as follows. In Section, we first discuss the necessary background on dihedral symmetries, the multiplicative structure of cellular integrals and modular forms with complex multiplication, then prove Theorem 1.1. In Section 3, we explain how to derive multiple binomial sum representations for the leading coefficients. In Section 4, we prove some preliminary results on combinatorial congruences, recall finite field hypergeometric series and then prove Theorem 1..

5 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 5 (7,,4,1,6,3,5 (7,,4,6,1,3,5 (7,,5,1,3,6,4 (7,,5,1,4,6,3 (7,3,6,,5,1,4 j,k=0 n ( n ( n ( ( n+k +j k +j k j n k j,k=0 n ( ( ( ( ( ( n n n+k n+j n+k+j n k j n n n k j j,k=0 j,k=0 [ n ( n ( ] n+k k n k=0 n ( n ( n ( ( n+k k +j k j n n n ( n ( n ( ( n+k n k+j k j n n Table. Convergent permutations σ 7 and leading coefficients A σ7 (n. Proof of Theorem Dihedral symmetries and multiplicative structures. Let Σ N bethesymmetric group on {1,,...,N}. A dihedral structure on {1,,...,N} is an equivalence class of permutations σ Σ N, where the equivalences are generated by (σ 1,σ,...,σ N (σ,σ 3,...,σ N,σ 1 (σ 1,σ,...,σ N (σ N,σ N 1,...,σ 1. A configuration on {1,,...,N} is an equivalence class [δ,δ ] of pairs (δ,δ of dihedral structures modulo the equivalence relations (δ,δ (σδ,σδ for σ Σ N. Associated to every permutation σ Σ N is the configuration [σ] := [id,σ]. Clearly, every configuration can be thus represented. Indeed, the configurations on {1,,...,N} can be identified with the double cosets D N \Σ N /D N, where D N are the dihedral permutations of Σ N. Henceforth, we will not make a distinction between permutations and configurations. The dual of a configuration [δ,δ ] is the configuration [δ,δ]. Equivalently, the dual of [σ] is [σ] := [σ 1 ]. For further details, we refer to [13, Section 3.1]. The notion of a configuration is important for our considerations because, up to a possible factor of ( 1 n, the cellular integral I σ (n only depends on the configuration [σ]. In [13, Section 6], Brown describes the following partial multiplication on pairs of dihedral structures. A pair of dihedral structures (δ,δ on {1,,...,N} is multipliable along the triple t = (t 1,t,t 3, where t 1,t,t 3 {1,,...,N} are distinct, if the elements t 1,t,t 3 are consecutive in δ and the elements t 1,t 3 are consecutive in δ. Let (α,α and (β,β be dihedral structures on X = {1,,...,N} and Y = {1,,...,M}, respectively. If (α,α is multipliable along s, and (β,β (the dual of (β,β is multipliable

6 6 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB along t, then the product (γ,γ := (α,α s,t (β,β is a pair of dihedral structures on the disjoint union Z of X and Y, where each of the three elements s 1,s,s 3 from X gets identified with the corresponding element t 1,t,t 3 from Y. The dihedral structure γ is the unique structure such that its restriction to X (respectively, Y coincides with α (respectively, β. We say that γ was obtained by shuffling α and β (more general such shuffles are described in [14, Section.3.1]. γ is likewise obtained by shuffling α and β. In order to work explicitly, we identify this disjoint union of X and Y with the set Z = {1,,...,M +N 3} in such a way that Y is the natural subset of Z, while the N 3 elements 1,,...,N of X, with s 1,s,s 3 removed, are identified, in that order, with the elements M + 1,M +,...,M +N 3 of Z. We note that different choices for this identification lead to the same configuration [γ,γ ]. Example.1. Forillustration, letusmultiply(α,α = ((1,,3,4,5,(1,3,5,,4 with(β,β = (α,α along s = (1,,3, t = (4,,5. Here, X = {1,,3,4,5}, Y = {1,,3,4,5} and Z = {1,,...,7}. The elements 1,,3,4,5 of Y are identified with the elements 1,,3,4,5 of Z, and the elements 1,,3,4,5 of X are identified with the elements 4,,5,6,7 of Z. In other words, we replace (α,α with (α,α = ((4,,5,6,7,(4,5,7,,6. In order to obtain γ, we shuffle α = (4,,5,6,7 and β = (1,,3,4,5 such that both dihedral structures are preserved (because of the conditions of multipliability there is a unique such shuffle. The result is γ = (1,,3,4,7,6,5. Likewise, shuffling α = (4,5,7,,6 and β = (1,3,5,,4 results in γ = (1,3,5,7,,6,4. We refer to [13, Example 6.3] for another example accompanied by a helpful illustration. The main interest in this partial multiplication stems from the following result. For its statement, recall that every configuration [γ,γ ] can be represented as [σ] = [id,σ]. Since, up to a sign, the cellular integral I σ only depends on the configuration [σ], we may write I [γ,γ ] = I σ. Likewise, we write A [γ,γ ] = A σ for the corresponding leading coefficients. Proposition. ([13], Proposition 6.5. Suppose that (γ,γ = (α,α s,t (β,β for some choice of s,t. Then, for all n 0, possibly up to a sign, In particular, A [γ,γ ](n = A [α,α ](na [β,β ](n. I [γ,γ ](n = I [α,α ](ni [β,β ](n. Example.3. Recall that, for σ 5 = (1,3,5,,4, the leading terms A σ5 (n of the corresponding cellular integrals are the Apéry numbers a(n. The configuration [σ 5 ] is represented by the pairs of dihedral structures (α,α = (β,β from Example.1. Let (γ,γ be the product of these two pairs along s,t as chosen in Example.1. Observe that, as a configuration, [γ,γ ] = [σ 7 ] with σ 7 = (1,3,7,5,,6,4. By Proposition., A σ7 (n = a(n. Similarly, many leading coefficients A σ (n are products of leading coefficients of lower order. However, it is not the case that products of leading coefficients are always leading coefficients.

7 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 7 Example.4. Consider the configuration σ 8 = (8,3,6,1,4,7,,5 featured in Theorem 1. and Section 3. This configuration is self-dual and not multipliable along any choice of triple t. Indeed, we confirm numerically that the product a(na σ8 (n is not a leading coefficient A σ10 (n for any configuration σ 10. Nevertheless, the next result guarantees that all positive integer powers of the Apéry numbers a(n are leading coefficients. Proposition.5. If ρ = (N 1,ρ 0,N,ρ 1,ρ,...,ρ N 3 Σ N, then A τ (n = a(na ρ (n, where τ = (N +1,ρ 0 +1,N +,N,ρ N 3 +1,ρ N 4 +1,...,ρ 1 +1,1. Note that τ is of the same shape of ρ, so that the result can be iterated to conclude that a(n λ A ρ (n is a leading coefficient for any integer λ 0. Proof. Let (α,α = ((1,,3,4,5,(1,3,5,,4 and (β,β = ((1,,...,N,ρ. Note that (α,α is multipliable along s = (1,,3, and that the dual of (β,β is multipliable along t = (N 1,ρ 0,N. Let (γ,γ = (α,α s,t (β,β. To compute this product, we proceed as in Example.1 and replace (α,α with (α,α = ((N 1,ρ 0,N,N+1,N+,(N 1,N,N+,ρ 0,N+1. Then, shuffling α = (N 1,ρ 0,N,N+ 1,N + and β = (1,,...,N such that both dihedral structures are preserved, we obtain γ = (1,,...,N 1,N +,N +1,N. Likewise, shuffling α = (N 1,N,N +,ρ 0,N + 1 and β = (N 1,ρ 0,N,ρ 1,ρ,...,ρ N 3 results in γ = (N 1,N +1,ρ 0,N +,N,ρ 1,ρ,...,ρ N 3. Since [α,α ] = [σ 5 ] and [β,β ] = [ρ], Proposition. implies that A [γ,γ ](n = a(na ρ (n, so that it only remains to observe that [γ,γ ] = [τ]. First, swapping N and N +, we find that Finally, as configurations, we have which is [τ], as claimed. [γ,γ ] = [(N 1,N +1,ρ 0,N,N +,ρ 1,ρ,...,ρ N 3 ]. [(N 1,N +1,ρ 0,N,N +,ρ 1,ρ,...,ρ N 3 ] = [(N,N +,ρ 0 +1,N +1,1,ρ 1 +1,ρ +1,...,ρ N 3 +1] = [(N +1,ρ 0 +1,N +,N,ρ N 3 +1,ρ N 4 +1,...,ρ 1 +1,1], Example.6. Let us illustrate with ρ = (4,,5,3,1. Since A ρ (n = a(n, Proposition.5 implies that A σ7 (n = a(n, σ 7 = (6,3,7,5,,4,1. Note that this configuration agrees with the one obtained in Example.3. Iterating Proposition.5, we find A σ9 (n = a(n 3, σ 9 = (8,4,9,7,,5,3,6,1, A σ11 (n = a(n 4, σ 11 = (10,5,11,9,,7,4,6,3,8,1.

8 8 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB In fact, the family of configurations in Example.6 can be made explicit as follows. For a positive integer M and configuration (a 1,a,...,a N, we write M +(a 1,a,...,a N for (M + a 1,M +a,...,m +a N. Corollary.7. Let M be an integer. Then where the configuration σ M+1 is A σm+1 (n = a(n M 1, σ M+1 = M +(M,0,M +1,M 1, (M,M 3,...,±1, 1,±,..., (M 1. Proof. The statement is clearly true for M =, in which case σ 5 = (4,,5,3,1. The claim then follows from Proposition.5 by induction... Modular forms with complex multiplication and Hecke characters. In this section, we recall some properties of modular forms with complex multiplication and Hecke characters. For more details, see [39]. Suppose ψ is a nontrivial real Dirichlet character with corresponding quadratic field K. A newform f(z = n=1 γ(nqn, where q := e πiz, has complex multiplication (CM by ψ, or by K, if γ(p = ψ(pγ(p for all primes p in a set of density one. By the work of Hecke and Shimura we can construct CM newforms using Hecke characters. Let K = Q( d be an imaginary quadratic field with discriminant D, and let O K be its ring of integers. For an ideal f O K, let I(f denote the group of fractional ideals prime to f. A Hecke character of weight k and modulo f is a homomorphism Φ : I(f C, satisfying Φ(αO K = α k 1 when α 1 (mod f. Let N(a denote the norm of the ideal a. Then, f(z := Φ(aq N(a = γ(nq n, a where the sum is over all ideals a in O K prime to f, is a Hecke eigenform of weight k on Γ 0 ( D N(f with Nebentypus χ(n = ( D Φ(nOK n. Here, ( a n n is the Kronecker symbol. Furthermore, f has CM by K. We call f the conductor of Φ if f is minimal, i.e., if Φ is defined k 1 modulo f then f f. If f is the conductor of Φ then f(z is a newform. From [39], we also know that every CM newform comes from a Hecke character in this way. We will see that Theorem 1.1 is a consequence of Corollary.7 and Corollary.11. Theorem.8. Let k be a positive integer. Then there exists a weight k CM newform S k (Γ 0 (3, if s even, f k (z =: γ k (nq n S k (Γ 0 (4,( 4, if k 1 (mod 4, n=1 S k (Γ 0 (16,( 4, if k 3 (mod 4, such that, for any odd prime p, {( 1 (x+y 1(k 1 [ (x+iy γ k (p = k 1 +(x iy k 1], if p 1 (mod 4, p = x +y, x odd, 0, if p 3 (mod 4. Proof. For n each equivalence class modulo 4, we will define a Hecke character Ψ k and construct the required CM newform f k, using the methodology outlined above. n=1

9 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 9 For an ideal f O K, let I(f denote the group of fractional ideals prime to f, and let J(f be the subset of principal fractional ideals whose generator is multiplicatively congruent to 1 modulo f, i.e., J(f = {(α I(f α 1 (mod f}. LetK = Q( 1 whichhasdiscriminant D = 4andwhoseringofintegers iso K = Z[ 1], which is a principal ideal domain. Therefore, all fractional ideals of K are also principal, and are of the form 1 m (α where m Z\{0} and α O K. Case 1: k 1 (mod 4. Let f = (1. Then I(f = J(f is the set of all fractional ideals. We define the Hecke character Φ k : I((1 C of weight k and conductor (1 by ( Φ 1 k m (α = ( α k 1. m Therefore, f k (z := Φ k (aq N(a = a is a CM newform and, for p an odd prime, γ k (p = n=1 γ k (nq n S k (Γ 0 (4,( 4 { (x+iy k 1 +(x iy k 1, if p 1 (mod 4, p = x +y, x odd, 0, if p 3 (mod 4. Case : k 3 (mod 4. Let f = (. Then and I(( = { 1 m (α m Z odd,α = x+iy Z[ 1],x and y different parity} J(( = { 1 m (α m Z odd,α = x+iy Z[ 1],x odd, y even}. We define the Hecke character Φ k : I(( C of weight k and conductor ( by ( k 1 Φ 1 k m (x+iy = ( 1 y. Therefore, f k (z := Φ k (aq N(a = a is a CM newform, and for p an odd prime, γ k (p = ( x+iy m n=1 γ k (nq n S k (Γ 0 (16,( 4 { (x+iy k 1 +(x iy k 1, if p 1 (mod 4, p = x +y, x odd, 0, if p 3 (mod 4. Case 3: k even. Let f = (+i. Then and I((+i = { 1 m (α m Z odd, α = x+iy Z[ 1], x and y different parity} J((+i = { 1 m (α m Z odd, α = x+iy Z[ 1], x odd, y even, such that if m x (mod 4 then y 0 (mod 4, otherwise y (mod 4}. We define the Hecke character Φ k : I((+i C of weight k and conductor (+i by ( Φ 1 k 1χ(x+iy k m (x+iy k 1 = ( 1 y ( 1 m 1, ( x+iy m

10 10 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB where χ(x+iy = ( 1 x+y 1 { 1, if x odd, y even, i, if x even, y odd. Therefore, f k (z := Φ k (aq N(a = γ k (nq n S k (Γ 0 (3 a n=1 is a CM newform and, for p an odd prime, {( 1 x+y 1 [ (x+iy k 1 +(x iy k 1], if p 1 (mod 4, p = x +y, x odd, γ k (p = 0, if p 3 (mod 4. Remark.9. When s odd, an alternative expression for the CM newform f k (z constructed in Theorem.8 is given by the binary theta series f k (z = 1 ( 1 m(k 1/ (n im k 1 q n+m. (.1 4 (n,m Z The formulation (.1 was recently used to relate the L-series of f k at k 1 to an interpolated version of the leading coefficients in Corollary.7. For further details, see [36]. Corollary.10. Let k be a positive integer and let f k (z =: n=1 γ k(nq n be the weight k CM newform described in Theorem.8. Then, for any odd prime p and integer m 1, γ k (p m = m 1 t=0 ( m t p m (k 1, if p 1 (mod 4 and m even, ( m p t(k 1 m γ (m t(k 1+1 (p+ 0, otherwise. Proof. This follows from a simple application of the binomial theorem on the expression for γ k (p given in Theorem.8, while noting that ( 1 x+y 1 (x+iy ( 1 x+y 1 (x iy = p. Corollary.11. Let l 1be apositive integer and define k := l+1. Let f k (z =: n=1 γ k(nq n be the weight k CM newform described in Theorem.8. Let a(n = n n ( n+j j=0( j j be the Apéry numbers introduced in (1.3. Then, for primes p 5, a( l γ k (p (mod p. Proof. By (1.4 and Corollary.10, We can now prove Theorem 1.1. a( l γ 3 (p l γ k (p (mod p. Proof of Theorem 1.1. Let l 1 be a positive integer and define k = l + 1 and N = l + 3. Consider f k (z =: n=1 γ k(nq n as in Theorem.8. By Corollary.7, there exists a convergent σ N whose leading coefficient A σn (n satisfies Thus, the result then follows from (. and Corollary.11. A σn (n = a(n l. (.

11 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS Multiple binomial sums 3.1. Proving multiple binomial sum representations for A σ (n. In this section, we outline how to obtain multiple binomial sum representations for the leading coefficients A σ (n. Our discussion applies to any configuration, but we proceed with the configuration σ 8 = (8,3,6,1,4,7,,5 for N = 8 which is relevant to Theorem 1.. This is the self-dual configuration 8 π 6 in Brown s notation [13, Section ]. The leading coefficients A σ8 (n have initial terms 1, 33, 899, , , ,... and the following particularly symmetric binomial sum representation. Proposition 3.1. For the configuration σ 8 = (8,3,6,1,4,7,,5, the leading coefficients are n 4 ( ( n n+ki A σ8 (n =. (3.1 k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 Proof. For σ = σ 8, consider the cellular integral I σ (n = S 8 f n σ ω σ where f σ = (z 1 z (z z 3 (z 3 z 4 (z 4 z 5 (z 5 z 6 (z 6 z 7 (z 7 z 8 (z 8 z 1 (z 8 z 3 (z 3 z 6 (z 6 z 1 (z 1 z 4 (z 4 z 7 (z 7 z (z z 5 (z 5 z 8. (3. As in Section 1, we let z 1 = 0, z 7 = 1, z 8 =. Then, in the coordinates t 1 = z,t = z 3,t 3 = z 4,t 4 = z 5,t 5 = z 6, we have f σ = ( t 1(t 1 t (t t 3 (t 3 t 4 (t 4 t 5 (t 5 1, (t t 5 (t 5 ( t 3 (t 3 1(1 t 1 (t 1 t 4 and dt 1 dt dt 3 dt 4 dt 5 ω σ = (t t 5 (t 5 ( t 3 (t 3 1(1 t 1 (t 1 t 4. The domain S 8 of integration then consists of all (t 1,t,...,t 5 R 5 such that 0 < t 1 < t <... < t 5 < 1. Algorithmic approaches to computing explicit period integrals, such as I σ (n for specific values of n, are described in [10] or [37]. In particular, Panzer implemented his symbolic integration approach [37] using hyperlogarithms in a Maple package called HyperInt. Using this package, we explicitly evaluate I σ (n in terms of multiple zeta values for several small values of n, and obtain: I σ (0 = 16ζ(5 8ζ(3ζ( I σ (1 = 33I σ (0 43ζ(3+316ζ( 6 I σ ( = 899I σ ( ζ( ζ( I σ (3 = I σ ( ζ( ζ( I σ (4 = I σ ( ζ( ζ( As proven by Brown in [13, Section 4], the integrals I σ (n satisfy a linear recurrence with polynomial coefficients. Slightly more specifically, the ordinary generating function F σ (x := I σ (nx n n=0

12 1 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB satisfies a Picard Fuchs differential equation. Again, in each specific instance, this differential equation can be obtained algorithmically. A particularly efficient such approach is an extension of the Griffiths Dwork reduction method due to Lairez [8] for computing periods of rational functions. To apply this method, we observe that where F σ (x = S 8 ω σ 1 f σ x = dt 1 dt dt 3 dt 4 dt 5 S 8 A Bx A = (t t 5 (t 5 ( t 3 (t 3 1(1 t 1 (t 1 t 4, B = ( t 1 (t 1 t (t t 3 (t 3 t 4 (t 4 t 5 (t 5 1. Lairez s method (implemented in Magma then successfully determines a Fuchsian differential equation of order 7 satisfied by F σ (x. This differential equation has a two-dimensional space of analytic solutions around x = 0. As a consequence, this differential equation together with the two values I σ (0, I σ (1, explicitly obtained above, determines the values of the cellular integrals I σ (n for all n. Alternatively, the differential equation for the generating function F σ (x translates directly into a recurrence of order 1 for the coefficients I σ (n and, hence, for the leading coefficients A σ (n. That is, we find that 151n(n+1 5 (n+1a σ (n+1 = 11 j=0 q j (na σ (n j for certain polynomials q j (x Z[x] of degree 7. To complete the proof, it therefore only remains to verify that the numbers n 4 ( ( n n+ki B σ (n := n k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 satisfy the same recurrence with matching initial values. This can be done algorithmically using, for instance, creative telescoping. In practice, the fact that B σ (n is a triple sum makes the computation of the recurrence rather challenging. Yet, Koutschan s Mathematica package HolonomicFunctions [7] is able to determine a fourth order linear recurrence for B σ (n (see below for more information on this recurrence. We then verify that this recurrence is a right factor of the earlier recurrence of order 1 for A σ (n. Since the first 1 initial values match (in fact, additional reflection shows that two matching initial values suffice, we conclude that A σ (n = B σ (n. The proof of Proposition 3.1 demonstrates that any individual evaluation of leading coefficients in terms of binomial sums can, in principle, be algorithmically proven due to recent advances in symbolic computation. It is curious to note that we had to use Maple, Magma and Mathematica in that computation. Remark 3.. We find that A σ8 (n, with σ 8 as in Proposition 3.1, is the unique solution of a fourth order recurrence p 4 (na n+4 +p 3 (na n+3 +p (na n+ +p 1 (na n+1 +p 0 (na n = 0,,

13 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 13 with initial conditions a 0 = 1,a 1 = 33 and a j = 0 for j < 0. Here, the coefficients p j (n are polynomials of degree 15, satisfying p j (n = p 4 j ( 5 n, for j {1,,3,4}. The latter relation is a consequence of the fact that the configuration σ 8 is self-dual (see [13, Section 4]. 3.. Finding multiple binomial sum representations for A σ (n. In this section, we illustrate how binomial sum representations can be found for any convergent configuration σ. As in the previous section, we proceed with the configuration σ = σ 8 = (8,3,6,1,4,7,,5. Consider the cellular integral I σ (n = fσ n ω σ, S 8 where f σ is as in (3.. Because the action of PGL is triply transitive, we may also make the convenient choice z σ(6 = z 7 = 1, z σ(7 = z = 0, z σ(8 = z 5 =. Then, in the coordinates t 1 = z 1,t = z 3,t 3 = z 4,t 4 = z 6,t 5 = z 8, we have and We then substitute f σ = (t 1( t (t t 3 (t 4 1(1 t 5 (t 5 t 1 (t 5 t (t t 4 (t 4 t 1 (t 1 t 3 (t 3 1 ω σ = dt 1 dt dt 3 dt 4 dt 5 (t 5 t (t t 4 (t 4 t 1 (t 1 t 3 (t 3 1. x 1 = t 5 t, x = t t 4, x 3 = t 4 t 1, x 4 = t 1 t 3, x 5 = t 3 1. Observe that t 3 = x 5 +1, t 1 = x 4 +x 5 +1 and, likewise, t 4 = x 3 +x 4 +x 5 +1, t = x +x 3 +x 4 +x 5 +1, t 5 = x 1 +x +x 3 +x 4 +x 5 +1, so that, up to a sign in ω σ, f σ = x 4,6x,6 x,4 x 3,5 x 1,5 x 1,3 x 1 x x 3 x 4 x 5, ω σ = dx 1dx dx 3 dx 4 dx 5, x 1 x x 3 x 4 x 5 where x i,j = x i +x i x j with x 6 = 1. At this point, it is natural to also consider the integral J σ (n = 1 (πi 5 x i =ε f n σω σ = 1 (πi 5 x i =ε ( x4,6 x,6 x,4 x 3,5 x 1,5 x 1,3 x 1 x x 3 x 4 x 5 n dx 1 dx dx 3 dx 4 dx 5 x 1 x x 3 x 4 x 5, where ε is chosen sufficiently small so that the integrals converge. By the residue theorem, this integral evaluates to J σ (n = [(x 1 x x 3 x 4 x 5 n ](x 4,6 x,6 x,4 x 3,5 x 1,5 x 1,3 n. In particular, the values J σ (n are nonnegative integers. By the principle of creative telescoping, we can derive a linear recurrence which is satisfied by both I σ (n and J σ (n. In fact, we are going to show that the leading coefficients A σ (n of I σ (n are equal to J σ (n. Likely, one can prove this equality in a general uniform fashion (for instance, following the approach of [18], as suggested by Dupont. For our purposes, it suffices to observe that both sequences satisfy a

14 14 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB common recurrence and agree to sufficiently many terms. It remains to express J σ (n in terms of a multiple binomial sum. Here, it will be convenient to not specialize x 6. In order to find an explicit formula for J σ (n, we expand Λ = (x 1,5 x,6 x 1,3 x,4 x 3,5 x 4,6 n using the binomial theorem and then extract the coefficient of (x 1 x x 3 x 4 x 5 x 6 n. Some care needs to be applied at this stage, because the order in which terms are expanded can have a considerable influence on the final binomial sum. For instance, the number of summations can vary substantially. In the present case, it is natural to first expand x n 1,5 = n ( n k 1 =0 k 1 x k 1 1,3 xn k 1 4,5, x,6 = n ( n k =0 and, in a second step, x n+k 1 1,3 as well as x n+k 4 4,6, so that Λ equals ( ( n n x n+k 1 1,3 x n,4 xn 3,5 xn+k 4 4,6 x n k 4,3 x n k 1 4,5 = k 1,k 4 k 1 k 1,k 4,l 1,l 4 k 1 k 4 k 4 l 1 l 4 k x k 4,6 xn k,3, ( ( ( ( n n n+k1 n+k4 x n,4x n 3,5x n k 4+l 1,3 x n k 1+l 4 4,5 x n+k 1 l 1 1 x n+k 4 l 4 6. Since J σ (n is the coefficient of (x 1 x x 3 x 4 x 5 x 6 n in Λ, we conclude that J σ (n is the coefficient of (x x 3 x 4 x 5 n in ( ( ( ( n n n+k1 n+k4 x n,4 xn 3,5 xn k 4+k 1,3 x n k 1+k 4 4,5. (3.3 k 1,k 4 k 1 k 4 We next expand x n,4 and xn 3,5 to obtain k 1 x n,4 xn 3,5 xn k 4+k 1,3 x n k 1+k 4 4,5 = ( ( n n x n k 4+k 1 +k,3 x n k 1+k 4 +k 3 4,5 x n k 3 3 x n k 4 k k 3 k,k 3 = k,k 3,l,l 3 k 4 ( ( ( ( n n n k4 +k 1 +k n k1 +k 4 +k 3 k k 3 l x l x n k 3+n k 4 +k 1 +k l 3 x n k +n k 1 +k 4 +k 3 l 3 4 x l 3 5. The coefficient of (x x 3 x 4 x 5 n in that sum is ( ( ( ( n n n+k n+k3, n n k,k 3 k k 3 subject to the constraint k 1 +k = k 3 +k 4. Finally, combined with (3.3, we conclude that n 4 ( ( n n+ki J σ (n =, n k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 l 3

15 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 15 which is the binomial sum of Proposition 3.1. A similar procedure has been carried out to find explicit binomial expressions for the remaining 897 leading coefficients A σn (n for 7 N 10. Again, these expressions are available upon request. 4. Proof of Theorem Congruences for binomial coefficients and harmonic sums. For a nonnegative integer n, we define the harmonic sum H n by n 1 H n := j and H 0 := 0. We first recall some elementary congruences (see Section 7.7, Theorems 133, 13 and 116 in [4]: For p an odd prime, we have ( ( j=1 4! 1 (mod p, (4.1 ( 1 p (mod p, (4. ( 1 p (mod p (4.3 p and, for primes p > 3, H 0 (mod p. (4.4 We will also need the following result, which follows easily from (4.3 and (4.4. Lemma 4.1. For a prime p > 3, p 1 ph (mod p. For nonnegative integers n, let (a n := a(a+1(a+ (a+n 1 denote the rising factorial, with (a 0 := 1. Let p be an odd prime. We note that, for 0 k (/, ( +k (1+k(+k ( +k (k +1 (mod p, k!! and, similarly, ( ( ( +k ( 1 k (k+1 k k! ( as well as ( ( k ( 1 +k k Lemma 4.. For p an odd prime, ( (k +1 1 (mod p. k=0 (mod p, (4.5 (mod p. (4.6

16 16 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB Proof. Let m := and write (k +1 m = b 0 + m t=1 b t k t for appropriate integers b t. Recall that, for positive integers s, { j s 1 (mod p, if ( s, 0 (mod p, otherwise. Therefore, j=1 (k +1 m = pb 0 + k=0 Since (k +1 m is monic, b = 1. Lemma 4.3. For p an odd prime, [ 4 k 1,k,k 3,k 4 =0 ki = m k=0 t=1 ( +1 b t k t = pb 0 + m b t t=1 k=1 k t b (mod p. (4.7 ] ( H +k 4 H k 4 0 (mod p. (4.8 Proof. For brevity, we again write m =. Let 0 k and consider Similarly, Let (m k +1 m = ( k +1( k + (p k( k H m k H m k = ( k( +1 k ( k( 1 k ( 1 m (k+1 m (mod p. 1 k + 1 p k k + 1 k (H m+k H k (mod p. S = m k 1,k,k 3,k 4 =0 ki = [ 4 ] ( +1 m (H m+k4 H k4 k+ + 1 k+1 k k be the left-hand side of (4.8. By replacing with m, for i {1,,3,4}, we see that [ m 4 ] S = (m +1 m (H m k4 H m k4 S (mod p. k 1,k,k 3,k 4 =0 ki = Hence, the sum S must vanish modulo p.

17 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 17 Lemma 4.4. For p an odd prime, k 1,k,k 3,k 4 =0 ki = 4 ( ( + 4 k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 ( ( + (mod p. Proof. Once more, m :=. Replacing k 1 and k with m k 1 and m k, respectively, yields m 4 ( ( m m+ki m ( ( 4 ( ( m m ki m m+ki =. m k 1,k,k 3,k 4 =0 ki = k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 Let k := (k 1,k,k 3,k 4 and define ( ( 4 ( ( m m ki m m+ki f(p,k := m and g(p,k := i=3 4 ( ( m m+ki. Then it suffices to prove m m f(p,k g(p,k (mod p. k 1,k,k 3,k 4 =0 k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 k 1 +k =k 3 +k 4 It follows from (4.6 that ( ( m k m k = m k m ( m+k ( 1 m k i=3 (mod p. Consequently, letting m k := (m k 1,m k,m k 3,m k 4, we have that f(p,k+f(p,m k g(p,k g(p,m k [ 4 ( ][ m ( m ki ( ][ m+ki 4 ( m+ki = m 0 (mod p. i=3 4 ( ] m ki i=3 m Therefore, m m [f(p,k g(p,k] = [f(p,k+f(p,m k g(p,k g(p,m k] k 1,k,k 3,k 4 =0 k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 k 1 +k =k 3 +k 4 0 (mod p.

18 18 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB 4.. Multiplicative characters and finite field hypergeometric functions. Let F p denote the group of multiplicative characters of F p. We extend the domain of χ F p to F p, by defining χ(0 := 0 (including the trivial character ε p and denote χ as the inverse of χ. When p is odd we denote the character of order of F p by φ p. We will drop the subscript p if it is clear from the context. We recall the following orthogonality relation. For χ F p, we have {, if χ = ε, χ(x = 0, if χ ε. x F p (4.9 For A, B F p, let ( A := B( 1 B p x F p A(xB(1 x. Then, for A 0,A 1,...,A n, B 1,...,B n F p and x F p, the finite field hypergeometric function of Greene [3] is defined as ( A0, A 1,..., A n ( A0 χ n n+1f n x B 1,..., B n p := p χ F p We consider the case where A i = φ p for all i and B j = ε p for all j and write ( φp, φ n+1f n (x = n+1 F p,..., φ p n ε p,..., ε p x χ p ( Ai χ χ(x. B i χ for brevity. Let Z p denote the ring of p-adic integers and Z p its group of units. We define the Teichmüller character to be the primitive character ω : F p Z p satisfying ω(x x (mod p for all x {0,1,...,}. In fact, ω(x x pn 1 (mod p n (4.10 for all n 1. Theorem 1. is now implied by Proposition 3.1 and the following result. Theorem 4.5. Let η 1 (z =: γ(nq n n=1 be the unique newform in S 6 (Γ 0 (4. Then, for p an odd prime, γ(p 4 k 1,k,k 3,k 4 =0 k 1 +k =k 3 +k 4 ( ( + (mod p. (4.11 Proof. We confirm that (4.11 holds for p = 3 and assume p 5 henceforth. For the Legendre family of elliptic curves E λ, given by E λ : y = x(x 1(x λ, λ {0,1},

19 we define SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 19 a(p,λ := p+1 #E λ (F p. Then, from [19, Proposition.1] and using the fact that dim(s 6 (Γ 0 (4 = 1, we get that for p an odd prime, From [, Lemma.1] we have ( γ(p = 3 a(p,λ 4 3pa(p,λ +p λ= = p 3 +p ( 3 a(p,λ 4 3pa(p,λ. λ= a(p,λ = p p 3, λ= so that, combining these two equations, In [6, Section 4], Koike showed that for λ 0,1, γ(p = p 3 4p 9p 3 a(p,λ 4. (4.1 λ= a(p,λ = φ( 1 p F 1 (λ. We now recall a couple of facts about finite field hypergeometric functions from [30, Proposition 3] and [3, Theorem 4.] respectively: and, for λ 0, Therefore, p F 1 (1 = φ( 1, F 1 (λ = φ(λ F 1 ( 1 λ a(p,λ 4 = [ φ( 1 p F 1 (λ] 4 λ= λ= = [ φ( 1 p F 1 (λ] 4 [ φ( 1 p F 1 (1] 4 λ=1. = [ φ( 1 p F 1 (λ] 4 1 λ=1 = λ=1 [ ( 1 ] 4 φ( λ p F 1 1. (4.13 λ

20 0 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB By Theorem.4 in [9] (see also Theorem 1.1 in [35] and (4.10, we get that, for λ 0, φ( λ p F 1 ( 1 λ = φ( 1 λ p F 1 ( 1 λ (p+1 (p+1 (/ j=0 (/ j=0 ( ( j ( ( j +j j +j j For brevity, we again let m = (/ and write f(p,k 1,k,k 3,k 4 := Combining (4.13 and (4.14 yields λ= λ=1 ( 1 j λ jp( ( 1+jp H +j H j ( ( ( 1 j 1+jp H +j H j ω(λ j (mod p. ( ( ( m m+ki ( 1 (1+p (H m+ki H ki. j=0 m ( ( m m+j 1+ a(p,λ 4 (p+1 ( 1 j (1+jp(H m+j H j ω j (λ j j (4p+1 (3p+1 m k 1,k,k 3,k 4 =0 m k 1,k,k 3,k 4 =0 ki 0 (mod f(p,k 1,k,k 3,k 4 ω k 1+k +k 3 +k 4 (λ λ=1 f(p,k 1,k,k 3,k 4 (mod p, ( where we have used (4.9 and the facts that (p+1 4 4p+1 (mod p and (4p+1( (3p+1 (mod p. Accounting for (4.15 in (4.1, we now have γ(p 9p +(3p+1 m k 1,k,k 3,k 4 =0 ki 0 (mod f(p,k 1,k,k 3,k 4 (mod p. In order that 0 (mod, the sum must be either 0, or (. In the case = 0, we necessarily have k 1 = k = k 3 = k 4 = 0, which contributes f(p,0,0,0,0 = 1. In the third case, that is = (, we necessarily have k 1 = k = k 3 = k 4 = = m.

21 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 1 Applying (4., (4.4 and Lemma 4.1, we evaluate [( ( m f (p,m,m,m,m = m m [( m ( 1 m (1+p((H H m (1 p(h H m [ ( 1 m p (1+pH m ] 4 [( 1 m (1 ph m (1+pH m ] 4 [( 1 m ] 4 1 (mod p. The only other possibility is =. Therefore, m γ(p 3p+(3p+1 f(p,k 1,k,k 3,k 4 (mod p. k 1,k,k 3,k 4 =0 ki = Noting that H m+j H j Z p for 0 j m, we have 4 ( ( m m+ki f(p,k 1,k,k 3,k 4 = ( 1 (1+p (H m+ki H ki and so γ(p 3p + p [ 4 ( ( ]( m m+ki 1+p m k 1,k,k 3,k 4 =0 ki = g(p,k 1,k,k 3,k 4 + ] 4 ] 4 4 (H m+ki H ki m k 1,k,k 3,k 4 =0 ki = where, for brevity, [ 4 ( ( ]( m m+ki g(p,k 1,k,k 3,k 4 := 3+ We will now show that m k 1,k,k 3,k 4 =0 ki = 4 ( ( m m+ki 4 (H m+ki H ki. (mod p, (mod p, (4.16 g(p,k 1,k,k 3,k 4 3 (mod p. (4.17 Using (4.5 and (4.1, we note that, when =, [ 4 ( ]( g(p,k 1,k,k 3,k 4 ( 1 k (ki +1 4 i m 3+ (H m+ki H ki (m! [ 4 ]( 4 ( +1 m 3+ (H m+ki H ki (mod p. (4.18

22 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB Let us define g 1 (p,k 1,k,k 3,k 4 [ 4 ]( 3 := ( +1 m 3+ (H m+ki H ki (k 1 +k +k 3 (H m+k4 H k4 (4.19 as well as [ 4 ] g (p,k 1,k,k 3,k 4 := ( +1 m ((H m+k4 H k4. Then, using (4.18 and Lemma 4.3, we get that m k 1,k,k 3,k 4 =0 ki = g(p,k 1,k,k 3,k 4 m k 1,k,k 3,k 4 =0 ki = m k 1,k,k 3,k 4 =0 ki = g 1 (p,k 1,k,k 3,k 4 + m k 1,k,k 3,k 4 =0 ki = g 1 (p,k 1,k,k 3,k 4 (mod p. g (p,k 1,k,k 3,k 4 If m < k < p, then H m+k H k 1 p Z p and (k+1 m pz p. Therefore, to establish (4.17 it suffices to prove k 1,k,k 3,k 4 =0 ki = g 1 (p,k 1,k,k 3,k 4 3 (mod p. (4.0 Let k 4 := k 1 k k 3 and consider, for 1 i 3, P( := d [ ki ( +1 m (k 4 +1 ] 4m m = bi,0 + b i,t ki t d t=1 = ( +1 ] m +(k4 +1 m d d [ (k4 +1 m d [ ki ( +1 ] m. d

23 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 3 Now, d [ ki ( +1 ] m = (ki +1 m dk +k d [ i (ki +1 ] m i d [ ] = ( +1 m +k d i ( +1 m ( +1 m d = ( +1 m +k m m i ( +1 m ( +r = ( +1 m [ 1+ m s=1 s=1 r=1 r s ] 1 ( +s and, similarly, So, = ( +1 m[1+ (H m+ki H ki ], d [ (k4 +1 ] d m = (k4 +1 m [(k 4 +1 m ] d d d = (k 4 +1 m [( k 1 k k 3 +1 m ] d m m = (k 4 +1 m ( k 1 k k 3 +r s=1 r=1 r s = (k 4 +1 m (H m+k 4 H k4. P( = ( +1 m (k 4 +1 m[1+ (H m+ki H ki (H m+k4 H k4 ]. (4.1 Noting (4.7, we see that k 1 =0 P( = pb i,0 + 4m k 1 =1 t=1 By definition of P(, we have that b i,t k t i 4m t=1 b i,t k 1 =1 ( +1 m (k 4 +1 m = b i,0 + t b i, b i,( (mod p. 4m t=1 b i,t t+1 kt i, which is monic with integer coefficients. Thus, p b i, and b i,( =. Therefore, =0 P( b i, b i,( 1 (mod p. (4.

24 4 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB Considering (4.19, accounting for (4.1 and (4., and applying Lemma 4., we get that 3 3 g 1 (p,k 1,k,k 3,k 4 = (k j +1 m P( k 1,k,k 3,k 4 =0 ki = = k 1,k,k 3 =0 3 3 k j =0 j i k j =0 j i j=1 j i 3 (k j +1 m j=1 j i =0 P( 3 (k j +1 3 m 1 = 3 (mod p. This establishes (4.0, which in turn establishes (4.17. Now, accounting for (4.17 in (4.16, we see that 4 ( ( γ(p + (mod p. k 1,k,k 3,k 4 =0 ki = j=1 j i Applying Lemma 4.4 completes the proof. Acknowledgements The first and third authors are supported by a grant from the Simons Foundation (#35339, Dermot McCarthy; #514645, Armin Straub. The second author would like to thank both the Institut des Hautes Études Scientifiques and the Max-Planck-Institut für Mathematik for their support and Francis Brown for his encouragement during the initial stages of this project. He also thanks the organizers of the workshop Hypergeometric motives and Calabi-Yau differential equations, January 8 7, 017 at the MATRIX institute, The University of Melbourne at Creswick, for the opportunity to discuss these results. The authors are grateful to Clément Dupont for several insightful and helpful comments on cellular integrals. References [1] S. Ahlgren, Gaussian hypergeometric series and combinatorial congruences, Symbolic computation, number theory, special functions, physics and combinatorics (Gainesville, FL, 1999, 1 1, Dev. Math., 4, Kluwer Acad. Publ., Dordrecht, 001. [] S. Ahlgren, The points of a certain fivefold over finite fields and the twelfth power of the eta function, Finite Fields Appl. 8 (00, no. 1, [3] S. Ahlgren, K. Ono, A Gaussian hypergeometric series evaluation and Apéry number congruences, J. Reine Angew. Math. 518 (000, [4] G. Almkvist, C. van Enckevort, D. van Straten and W. Zudilin, Tables of Calabi-Yau equations, preprint available at [5] T. Amdeberhan, R. Tauraso, Supercongruences for the Almkvist-Zudilin numbers, Acta Arith. 173 (016, no. 3, [6] B. Aspvall, F. Liang, The dinner table problem, Technical Report STAN-CS-80-89, Computer Science Department, Stanford University, Stanford, California, 1980.

25 SEQUENCES, MODULAR FORMS AND CELLULAR INTEGRALS 5 [7] F. Beukers, A note on the irrationality of ζ( and ζ(3, Bull. London Math. Soc. 11 (1979, no. 3, [8] F. Beukers, Some congruences for the Apéry numbers, J. Number Theory 1 (1985, no., [9] F. Beukers, Another congruence for the Apéry numbers, J. Number Theory 5 (1987, no., [10] C. Bogner, F. Brown, Feynman integrals and iterated integrals on moduli spaces of curves of genus zero, Commun. Number Theory Phys. 9 (015, no. 1, [11] J. Broedel, O. Schlotterer and S. Stieberger, Polylogarithms, multiple zeta values and superstring amplitudes, Fortschr. Phys. 61 (013, no. 9, [1] F. Brown, Multiple zeta values and periods of moduli spaces M 0,n, Ann. Sci. Éc. Norm. Supér. (4 4 (009, no. 3, [13] F. Brown, Irrationality proofs for zeta values, moduli spaces and dinner parties, Mosc. J. Comb. Number Theory 6 (016, no. 3, [14] F. Brown, S. Carr and L. Schneps, The algebra of cell-zeta values, Compos. Math. 146 (010, no. 3, [15] H.H. Chan, S. Cooper and F. Sica, Congruences satisfied by Apéry-like numbers, Int. J. Number Theory 6 (010, no. 1, [16] S. Cooper, Sporadic sequences, modular forms and new series for 1/π, Ramanujan J. 9 (01, no. 1-3, [17] M. Coster, Supercongruences, Ph.D. thesis, Universiteit Leiden, [18] C. Dupont, Odd zeta motive and linear forms in odd zeta values, with a joint appendix with Don Zagier. Compos. Math. 154 (018, no., [19] S. Frechette, K. Ono and M. Papanikolas, Gaussian hypergeometric functions and traces of Hecke operators, Int. Math. Res. Not. 004, no. 60, [0] I. Gessel, Some congruences for Apéry numbers, J. Number Theory 14 (198, no. 3, [1] V. V. Golyshev, D. Zagier, Proof of the gamma conjecture for Fano 3-folds of Picard rank 1, Izv. Math. 80 (016, no. 1, [] A.B. Goncharov, Y.I. Manin, Multiple ζ-motives and moduli spaces M 0,n, Compos. Math. 140 (004, no. 1, [3] J. Greene, Hypergeometric functions over finite fields, Trans. Amer. Math. Soc. 301 (1987, no. 1, [4] G.H. Hardy, E.M. Wright, An introduction to the theory of numbers, fifth edition, The Clarendon Press, Oxford University Press, New York, [5] G. Kalita, A. Chetry, Congruences for generalized Apéry numbers and Gaussian hypergeometric series, Res. Number Theory 3 (017, Art. 5, 15 pp. [6] M. Koike, Orthogonal matrices obtained from hypergeometric series over finite fields and elliptic curves over finite fields, Hiroshima Math. J. 5 (1995, [7] C. Koutschan, Advanced Applications of the Holonomic Systems Approach, PhD thesis, RISC, Johannes Kepler University, Linz, Austria, 009. [8] P. Lairez, Computing periods of rational integrals, Math. Comp. 85 (016, no. 300, [9] P. Loh, R. Rhodes, p-adic and combinatorial properties of modular form coefficients, Int. J. Number Theory (006, no., [30] K. Ono, Values of Gaussian hypergeometric series, Trans. Amer. Math. Soc. 350 (1998, no. 3, [31] R. Osburn, B. Sahu, Congruences via modular forms, Proc. Amer. Math. Soc. 139 (011, no. 7, [3] R. Osburn, B. Sahu, Supercongruences for Apéry-like numbers, Adv. in Appl. Math. 47 (011, no. 3, [33] R. Osburn, B. Sahu, A supercongruence for generalized Domb numbers, Funct. Approx. Comment. Math. 48 (013, part 1, [34] R. Osburn, B. Sahu, A. Straub, Supercongruences for sporadic sequences, Proc. Edinb. Math. Soc. ( 59 (016, no., [35] R. Osburn, C. Schneider, Gaussian hypergeometric series and supercongruences, Math. Comp. 78 (009, no. 65, [36] R. Osburn, A. Straub, Interpolated sequences and critical L-values of modular forms, preprint available at [37] E. Panzer, Algorithms for the symbolic integration of hyperlogarithms with applications to Feynman integrals, Comput. Phys. Commun. 188 (015,

26 6 DERMOT MCCARTHY, ROBERT OSBURN, ARMIN STRAUB [38] P. Poulet, Permutations, L Intermédiaire des Mathématiciens, 6 (1919, [39] K. Ribet, Galois representations attached to eigenforms with Nebentypus, Modular functions of one variable, V (Proc. Second Internat. Conf., Univ. Bonn, Bonn, 1976, pp Lecture Notes in Math., Vol. 601, Springer, Berlin, [40] D. P. Roberts, F. Rodriquez-Villegas and M. Watkins, Hypergeometric motives, preprint (017. [41] O. Schlotterer, S. Stieberger, Motivic multiple zeta values and superstring amplitudes, J. Phys. A 46 (013, no. 47, , 37 pp. [4] J. Stienstra, F. Beukers, On the Picard-Fuchs equation and the formal Brauer group of certain elliptic K3-surfaces, Math Ann. 71 (1985, no., [43] A. Straub, Multivariate Apéry numbers and supercongruences of rational functions, Algebra Number Theory 8 (014, no. 8, Department of Mathematics & Statistics, Texas Tech University, Lubbock, TX , USA School of Mathematics and Statistics, University College Dublin, Belfield, Dublin 4, Ireland Department of Mathematics and Statistics, University of South Alabama, 411 University Blvd N, Mobile, AL address: dermot.mccarthy@ttu.edu address: robert.osburn@ucd.ie address: straub@southalabama.edu