Arxiv:Math-Ph/0305039V2 22 Oct 2003

arXiv:math-ph/0305039v2 22 Oct 2003 n hysoe htteivratcnb eadda h Eich the as regarded be can invariant weight the that showed they and inva Witten–Reshetikhin–Turaev the of Ref expansion In ymptotic weight. integral half with especially forms” modular h Dedekind the h iwon fteqatmivrat Zagier’s invariant, quantum the of viewpoint the o h ou link torus the for n twsson[]ta ahe’ nain o h ou kn torus the for invariant Kashaev’s that [3] shown was it w and correspondence This knot. trefoil coincide the to for function, [14] Jones shown was which 6], [5, invariant Kashaev’s fteAndrews–Gordon the of oua rpry tcnb eadda limit a as regarded be can it property; modular fve ftenme theory, a number they the but of complement, view knot of the of volume hyperbolic the tor with asymptotic the an of that studies states hyperbolic, which not [5,14] Conjecture” is “Volume knot torus the As space. n ae twsfudta Zagier’s that ind found linearly was of number it the later of and bound upper an of function ating 2000 Date eetsuisrva nitmt oncinbtentequ the between connection intimate an reveal studies Recent oiae yorpeiu euto h ou knot torus the on result previous our by Motivated a 5 03 eie nOtbr1,20.T peri Commu in appear To 2003. 17, October on Revised 2003. 15, May : ahmtc ujc Classification. Subject Mathematics 3 ucinfrtetrslink torus the for function hwta h nain srgre salmtof limit a as regarded is invariant the that show oua omo weight of form modular a A UNU NAIN O OU IKADMDLRFORMS MODULAR AND LINK TORUS FOR INVARIANT QUANTUM / BSTRACT 2 nRf 9 airfrhrsuida“tag dniy rel identity” “strange a studied further Zagier 19, Ref. In . η fnto,adcaie oeo h ihe nerlwt h with integral Eichler the of role a clarified and -function, ecnie naypoi xaso fKsavsinvarian Kashaev’s of expansion asymptotic an consider We . T (2 , 2 m q sre,wihi ht eiswt weight with series theta is which -series, ) seFg )i hsatce esalso htteinvariant the that show shall We article. this in 1) Fig. (see 3 T q / (2 sre n oua forms. modular and -series 2 , hc srltdt the to related is which 2 m q ) esalgive shall We . sre with -series .I 1. 16,5M7 53,1F3. 11F23 05A30, 57M27, 11B65, AUIOHIKAMI KAZUHIRO NTRODUCTION q 1 being q q q q sre a rgnlycnetdwt gener- a with connected originally was -series en h oto nt fteEclrintegral Eichler the of unity of root the being su sre dniyrltdt hs nains and invariants, these to related identity -series c en the being T N (2) (2 t oto nt fteEclritga of integral Eichler the of unity of root -th m , sfrhrivsiae o h ou knot, torus the for investigated further as − 2 in ftePicrehmlg sphere, Poincaré homology the of riant arneadZge tde nas- an studied Zagier and Lawrence 9 . 2 m ii fKsavsivratcoincides invariant Kashaev’s of limit sko a o eatatv o the for attractive be not may knot us character. ihaseicvleo h colored the of value specific a with 1) + ot pnetVsiivivrat [17], invariants Vassiliev ependent N .Mt.Phys. Math. n. erte bobn rmtepoint the from absorbing rather re e nerlo oua omof form modular a of integral ler nu ntivratad“nearly and invariant knot antum t oto nt onie with coincides unity of root -th T esuyKsavsinvariant Kashaev’s study we , (2 1 tdt h afdrvtvsof half-derivatives the to ated ro h ooe Jones colored the of or t / , 2 2 m l-nerlwih.From weight. alf-integral spanning 1) + lohsanearly a has also m -dimensional snwre- now is FIGURE 1. Hopf link T (2, 2) and torus link T (2, 4)

garded as the half-integration or the Eichler integral of a modular form of weight 3/2. Remarkable is that this modular form is related to the su(2)m−2 character. It is noted that recent studies [20,21] reveal a relation with Ramanujan’s mock theta functions. We also propose a q-series identity, which is new as far as we know, and studyc an asymptotic expansion thereof.

This paper is organized as follows. In section 2 we construct the colored Jones polynomial for the torus link T (2, 2 m). Using the Jones–Wenzl idempotent, we give an explicit formula of the invariant. It is known [14] that Kashaev’s invariant coincides with a speciﬁc value of the colored Jones polynomial. This correspondence enables us to give an integral form of Kashaev’s invariant for the torus link T (2, 2 m) in section 3. We further give an asymptotic expansion of the invariant, and see that the invariant for T (2, 2 m) also has a nearly modular property. We give an explicit form of Kashaev’s invariant for this torus link using the enhanced Yang–Baxter operator. Combining these results we obtain an asymptotic expansion of a certain ω-series. In section 4 we introduce the q-series related to Kashaev’s invariant for torus link, and prove a new q-series identity. We study the modular property of these q-series, and discuss how Kashaev’s invariant for T (2, 2 m) may be regarded as the Eichler integral of a modular form with weight 3/2 which is the afﬁne su(2)m−2 character, in section 5. In the last section, we collect some examples.

2. COLORED JONES POLYNOMIAL FOR TORUS LINK (2, 2 m)

The N-colored Jones polynomial for torus knot T (m, p) was studied in Refs. 12,15. Following these methods, we compute the colored Jones polynomial for torus link T (2, 2 m) in this section. 2 PSfragWe use replacements the Jones–Wenzl idempotent, andPSfrag use following replacements formulae (see, e.g., Ref. 10);

a b a b

a2 b2−c2 (a+b−c)/2 a+b−c+ + =( 1) A 2 (2.1) − PSfrag replacements PSfrag replacements c c

a a a ∆ = c c (2.2) θ(a, b, c) b c: (a, b, c)Xis admissible b b where each label denotes a color, and we mean that A2(n+1) A−2(n+1) ∆ =( 1)n − , n − A2 A−2 − a b + c, ≤ (a, b, c) is admissible a + b + c is even, and b c + a, PSfrag replacements⇔  ≤ c a + b. ≤ We have a θ-net 

a θ(a, b, c)= b (2.3) c

which is given as ∆ ! ∆ ! ∆ ! ∆ ! θ(a, b, c)= x+y+z x−1 y−1 z−1 , ∆y+z−1! ∆z+x−1! ∆x+y−1! with

a = y + z, b = z + x, c = x + y.

Proposition 1. The N-colored Jones polynomial J (h; ) for the torus link = T (2, 2 m) is N K K given by N−1 N h JN (h; ) − 1 m(N 2−1)h mhj2+(m+ε) hj+ 1 hε 2 sh K = e 2 ε e 2 , (2.4) 2 J (h; ) N ε=±1 j=0 O X X where a parameter q is set to be q = A4 = eh, 3 and denotes unknot whose invariant is given by O sh(N h/2) J (h; )= . N O sh(h/2)

Proof. We ﬁrst apply eq. (2.2) in the torus link T (2, 2 m) (see Fig. 2), and untangle crossings recursively using eq. (2.1). We see that θ(a, b, c) vanishes at the end. We have

c 2m N−1− −2(N−1)+c−(N−1)2 + 1 c2 J (h; )= ∆ ( 1) 2 A 2 N K c − − − c:(N 1,N X1,c) is admissible 2 − A−2m(N −1) N 1 = A4mj(j+1) A2(2j+1) A−2(2j+1) . A2 A−2 − − j=0 X This proves eq. (2.4).

= ⇒ PSfrag replacements

n n n n

FIGURE 2. We apply eq. (2.2) to the torus link T (2, 4). We set n = N 1. −

3. KASHAEV INVARIANT AND ASYMPTOTIC EXPANSION

It is known [14] that Kashaev’s invariant is given from the colored Jones polynomial at a speciﬁc value h 2 π i/N. By use of a result of Prop. 1 we obtain an integral form of the invariant as → follows.

Proposition 2. The Kashaev invariant for torus link T (2, 2 m) is given by an integral form as

π i (m− 1 ) 3 −πi/4 T (2, 2 m) = e 2N m 8√2(m N) 2 e h iN 2 2 sh(4 mπNw) sh(4 πw) dww2 e8mπiNw +4mπNw , (3.1) × sh(4 mπw) ZC where denotes a path passing through the origin in the steepest descent direction. C 4 Proof. A proof is essentially same with those given in Refs. 2,7,8,16 for studies in asymptotic behavior of the quantum invariants.

To rewrite eq. (2.4) into an integral form, we use an integral formula

2 2 1 z ehw = dz e− h +2wz, √π h ZC where an integration path is an inﬁnite line passing through the origin in the steepest descent C direction. Then we get N−1 m2 m ε N h JN (h; ) − 1 mh(N 2−1)− +1 h mh(j+ + )2 2 sh K = e 2 4m ε e 2m 2 J (h; ) N ε=±1 j=0 O X X N−1 m2 z2 m ε 1 − 1 mh(N 2−1)− +1 h − +2(j+ + )z = e 2 4m ε dz e mh 2m √ π m h ε=±1 j=0 X X ZC m h z2 1 − 1 mhN 2+ h− − +Nz 2 sh(N z) sh(z/m) = e 2 4 4m dz e mh . √π m h sh(z) ZC Now we take a limit 2 π i h , → N to compute Kashaev’s invariant. As LHS vanishes in this limit, we take a derivative of both sides to have 3/2 1 N πi 1 πi iz2 2 sh(N z) sh(z/m) N (m− m )− 2 N( mπ +z) T (2, 2 m) N = e 2 4 dz z e 2 . h i 4 √2 π3 m sh(z) ZC Rescaling an integral variable, we obtain eq. (3.1).

Theorem 3. Kashaev’s invariant for the torus link = T (2, 2 m) has an asymptotic expan- hKiN K sion in N as →∞ m−1 m− 2 k2 3 πi− ( 1) πi 3/2 2 k k − πiN T (2, 2 m) e 4 2mN N ( 1) (k m) sin π e 2m h iN ≃ m − − m r k=1 X ∞ (m;0) k m− 2 − ( 1) πi Ek π i + e 2mN N , (3.2) k! 2 m N Xk=0 (m;0) where Ek is deﬁned from a generating function (see also eq. (3.16)) as ∞ m sh(z) E(m;0) = k z2k (3.3) sh(m z) (2 k)! Xk=0 1 m2 (3 7 m2)(1 m2) =1+ − z2 + − − z4 + . 6 360 ··· 5 Proof. When we decompose sh(4 mπNw) in the integrand (3.1) into (e4mπNw e−4mπNw)/2, the − integral reduces to I I up to constant where I and I are 1 − 2 1 2 2 sh(4 πw) I = dww2 e8mπN(iw +w) , 1 sh(4 mπw) ZC

2 sh(4 πw) I = dww2 e8mπNiw . 2 sh(4 mπw) ZC In I we deform an integration path + i . In this deformation we have contributions from 1 C → C 2 residues at w = k i for k =1, 2,..., 2 m 1. Then we get 4 m −

2 sh(4 πw) I =2 π i (Residues)+ dww2 e8mπN(iw +w) 1 · sh(4 mπw) CZ i + 2 2 m−1 2 k2 1 k 2 k − πiN i 8mπiNz2 sh(4 π z) = ( 1) k sin π e 2m + dz z + e 32 m3 − m 2 sh(4 m π z) Xk=1 ZC m−1 k2 1 k k − πiN 1 8mπiNz2 sh(4 π z) = ( 1) (k m) sin π e 2m + I dz e . 8 m2 − − m 2 − 4 sh(4 m π z) Xk=1 ZC In the last equality, we have used a symmetry z z of the integrand. Substituting a series ↔ − expansion (3.3) into above expression, we obtain an assertion of theorem.

Asymptotic expansion for the torus knot T (2, 2 m + 1) was studied in Ref. 3. In view of these results, a tail of asymptotic expansion of Kashaev’s invariant is given in an inﬁnite series of N −1 with coefﬁcients, whose generating function seems to be related to 1 s − . ∆K(s)

Here ∆K(s) is the Alexander polynomial for knot , and in a case of the torus link = T (2, 2 m) K K we have 1 s2m ∆K(s)= − . 1+ s

We have shown that the volume conjecture [5,14] is fulﬁlled for the torus link T (2, 2 m), 2 π lim log T (2, 2 m) N =0, N→∞ N h i because the torus link T (2, 2 m) is not hyperbolic. Rather we have interests in an asymptotic expansion of the q-series. To this aim, we compute the quantum invariant by another method; 6 explicit form of Kashaev’s invariants can be directly computed from a set of the enhanced Yang– Baxter operators [5] (see also Refs. 13, 14, 18);

1−(k−j+1)(ℓ−i) ij Nω i j Rkℓ = ∗ ∗ θ , (3.4a) (ω) − − (ω) (ω) − (ω) · k ℓ [ℓ k 1] [j−ℓ] [i j] [k−i] −1+(ℓ−i−1)(k−j) −1 ij Nω i j (R )kℓ = ∗ ∗ θ , (3.4b) (ω) (ω) − (ω) (ω) − · k ℓ [ℓ−k−1] [j ℓ] [i−j] [k i] µk = δ ω1/2, (3.4c) ℓ − k,ℓ+1 where we have deﬁned the N-th root of unity as 2 π i ω = exp , (3.5) N and means a complex conjugation. We have also used [x] 0, 1,...,N 1 modulo N, and ∗ ∈{ − } i k<ℓ j, ≤ ≤ i j j i k < ℓ, θ =1, if and only if  ≤ ≤ k ℓ ℓ j i k (with ℓ

n (ω) = (1 ωk), (3.6) n − Yk=1 (ω)m m , if m n 0, PSfrag replacements = (ω)n (ω)PSfragm−n replacements≥ ≥ (3.7) n  0, others. These operators are assigned to a projection of knot as follows;

i j i j

ij −1 ij Rkℓ = R kℓ = PSfrag replacementsk ℓ PSfrag replacements k ℓ

k ℓ

k −1 k µℓ = µ ℓ = k ℓ

7 Proposition 4. Kashaev’s invariant for the torus link = T (2, 2 m) is explicitly given by hKiN K T (2, 2 m) h iN m−2 1 cm− cm− (cm− +1) ci(ci+1) ci+1 = N ( 1) 1 ω 2 1 1 ω . (3.8) − ci N−1≥cm− ≥···≥c ≥c ≥0 i=1 ! X1 2 1 Y Proof. We get this result from a direct computation using R-matrix (3.4) for (1, 1)-tangle of (2, 2 m)-torus link. See Refs. 13,18.

We note that Kashaev’s invariant for the Hopf link T (2, 2) is given by

T (2, 2) = N. h iN

Combining Prop. 4 with Thm. 3, we obtain an asymptotic expansion of ω-series.

Corollary 5.

m−2 1 cm− cm− (cm− +1) ci(ci+1) ci+1 ( 1) 1 ω 2 1 1 ω − ci N−1≥cm− ≥···≥c ≥c ≥0 i=1 ! X1 2 1 Y m−1 m− 2 k2 3 πi− ( 1) πi 2 k k − πiN √N e 4 2mN ( 1) (k m) sin π e 2m ≃ m − − m r k=1 X ∞ (m;0) k m− 2 − ( 1) πi Ek π i + e 2mN . (3.9) k! 2 m N Xk=0 In the rest of this paper, we shall reveal a meaning of this asymptotic expansion from a point of view of the modular form. As a generalization of ω-series deﬁned by Kashaev’s invariant we introduce

N−1 (a) cm− 1 cm− (cm− +1) Y (ω)= ( 1) 1 ω 2 1 1 m − c ,...,cm− =0 1 X 1 m−2 2 ··· 2 a ··· m− ci+1 + δi,a ωc1 + +cm−2+c +1+ +c 2 , (3.10) × ci i=1 ! Y for m 2 and a = 0, 1,...,m 2. See that Kashaev’s invariant for T (2, 2 m) corresponds to a ≥ − case of a =0, T (2, 2 m) = N Y (0)(ω). (3.11) h iN · m (a) It is unclear whether the ω-series Ym (ω) for a = 0 represent the quantum invariant for any three 6 manifolds. 8 (a) Conjecture 1. Let Ym (ω) be deﬁned by eq. (3.10). An asymptotic expansion of this ω-series in N is given by →∞ π i (m−1−a)2 (a) e 2mN Y (ω) · m m−1 k2 3 πi 2 k k − πiN √N e 4 ( 1) (k m) sin (a + 1) π e 2m ≃ · m − − m r k=1 X ∞ E(m;a) π i k + k , (3.12) k! 2 m N Xk=0 (m;a) where a generalized Euler number Ek is given from a generating function as ∞ m sh (a + 1) z E(m;a) = k z2k (3.13) sh(m z) (2 k)! k=0 X 1 =(a +1)+ (1 + a) (1 + a)2 m2 z2 6 − 1 + (1 + a) (1 + a)2 m2 3(1+ a)2 7 m2 z4 360 − − 1 + (1 + a) (1 + a)2 m2 3(1+ a)4 18(1+a)2 m2 + 31 m4 z6 + . 15120 − − ··· The case a =0 of this conjecture is proved in Corollary 5.

Note that we have sh (a + 1)z ∞ = χ(a) (n) e−nz, (3.14) sh(m z) 2m n=0 X (a) where the odd periodic function χ2m(n) is written as n mod 2 m m 1 a m +1+ a others − − (3.15) χ(a) (n) 1 1 0 2m − Applying the Mellin transformation to eqs. (3.13) and (3.14), we have an expression of the gener- (a) alized Euler number in terms of the L-function associated to χ2m(n);

E(m;a) = m L( 2 k, χ(a) ) k · − 2m (2 m)2k m 1 a m +1+ a = m B − − B , (3.16) − 2 k +1 2k+1 2 m − 2k+1 2 m where Bn(x) is the n-th Bernoulli polynomial. It should be remarked that the colored Jones polynomial (2.4) for the torus link = T (2, 2 m) is rewritten using the periodic function as K 2mN 2 2 N h JN (h; ) − 1 2− − m +1 (0) k 2 m(N 1)h 4m h 4m h 2 sh K = e χ2m(k) e . (3.17) 2 JN (h; ) − O k=0 9 X Based on this expression, we ﬁnd that Kashaev’s invariant is given by

2mN m− 2 k2 1 − ( 1) πi 2 (0) πi T (2, 2 m) = e 2mN k χ (k) e 2mN . (3.18) h iN −4 m N 2m Xk=0 Later we shall clarify a relationship between above conjecture and the modular form.

4. q-SERIES IDENTITY

(a) In this section we study a q-series identity, which is closely related with Ym (ω) deﬁned in eq. (3.10). We use standard notation as in eqs. (3.5) – (3.7), but in this section we replace ω, the N-th primitive root of unity, with generic q.

We deﬁne the q-series

∞ (a) cm− 1 cm− (cm− +1) c +···+cm− K (x)= ( 1) 1 q 2 1 1 x 1 1 m − c ,...,cm− =0 1 X 1 m−2 c 2+···+c 2 +ca +···+cm− ci+1 + δi,a q 1 m−2 +1 2 , (4.1) × · ci i=1 ! Y for m 2 and a = 0, 1, 2,...,m 2. We simply replace N−1 in eq. (3.10) with an inﬁnite ≥ − cm−1=0 sum ∞ , though we have introduced an additional variable x. cm−1=0 P

P (a) Theorem 6. Let the q-series Km (x) be deﬁned in eq. (4.1). We have

∞ 2 2 (a) n −(m−1−a) n−(m−1−a) (a) 4m 2 Km (x)= χ2m(n) q x , (4.2) n=0 X (a) where χ2m(n) is a periodic function in eq. (3.15).

Proof. We prove this statement by showing that both sides satisfy the same q-difference equation.

(a) It is easy to see from a periodicity of the function χ2m(n) that RHS of eq. (4.2) solves a difference equation (see, e.g., Refs. 1,3,19),

∞ n2− m− −a 2 n− m− −a (a) a+1 a+1 (a) ( 1 ) ( 1 ) K (x)=1 q x + χ (n) q 4m x 2 m − 2m n=2m X =1 qa+1 xa+1 + xm q2m−1−a K(a)(q2 x). (4.3) − m 10 We shall show that LHS of eq. (4.2) also fulﬁlls a same difference equation. To this aim, we introduce

∞ 1 (a) cm−1 cm−1(cm−1+1) cm−1 K (x ,...,x − )= ( 1) q 2 x m 1 m 1 − m−1 c ,...,cm− =0 1 X 1 a−1 2 2 ci ci ci+1 ca ca ca+1 +1 q xi q xa × ci · ca i=1 ! Y m−2 2 ci +ci ci ci+1 q xi . (4.4) × ci i=a+1 ! Y See that by deﬁnition we have

(a) (a) Km (x)= Km (x,...,x). (4.5) m−1

(a) We use same symbol Km , but we believe there is no| confusion.{z }

To prove the assertion of the theorem, we use formulae for the q-binomial coefﬁcients; n +1 n n = qc + (4.6a) c c c 1 − n n = + qn+1−c . (4.6b) c c 1 − c +1 Applying eq. (4.6a) to a+1 in eq. (4.4), we get c a (a) (0) −1 −1 Km (x1,...,xm−1)= Km (q x1,...,q xa−1, xa,...,xm−1) (a−1) + q x K (x ,...,x − , q x , x ,...,x − ). (4.7) a · m 1 a 1 a a+1 m 1 On the other hand, using eq. (4.6b), we have

(a) (0) −1 −1 Km (x1,...,xm−1)= Km (q x1,...,q xa, xa+1,...,xm−1) (a−1) + q x K (x ,...,x , q x , x ,...,x − ) (4.8) a · m 1 a a+1 a+2 m 1 Another difference equation is given as follows;

(0) Km (x1,...,xm−1) ∞ ∞ m−2 1 m− m− m− cm−1 i i i c c 1 2 c 1(c 1+1) c (c +1) c i+1 =1+ ( 1) q xm−1 q xi − ci cm− =1 c ,...,cm− =0 i=1 ! X1 1 X 2 Y (m−2) =1 q x − K (q x ,...,qx − ). (4.9) − m 1 · m 1 m 1 11 We can prove eq. (4.3) by use of eqs. (4.7)—(4.9) as follows. Recursive use of eq. (4.7) gives

(a) Km (x1,...,xa, q xa+1,...,qxm−1) (0) −1 −1 = Km (q x1,...,q xa−1, xa, q xa+1,...,qxm−1) (0) −1 −1 + q x K (q x ,...,q x − , x − , q x ,...,qx − ) a · m 1 a 2 a 1 a m 1 2 (0) −1 −1 + q x − x K (q x ,...,q x − , x − , q x − ,...,qx − ) a 1 a · m 1 a 3 a 2 a 1 m 1 a (0) + + q x x K (q x ,...,qx − ). (4.10) ··· 1 ··· a · m 1 m 1 Substituting above equation for a = m 2 into eq. (4.9), we get − (0) −1 −1 Km (q x1,...,q xm−2, xm−1) (0) −1 −1 + q x − K (q x ,...,q x − , x − , q x − ) m 1 · m 1 m 3 m 2 m 1 2 (0) −1 −1 + q x − x − K (q x ,...,q x − , x − , q x − , q x − ) m 2 m 1 · m 1 m 4 m 3 m 2 m 1 m−2 (0) + + q x x − K (x , q x ,...,qx − ) ··· 2 ··· m 1 · m 1 2 m 1 m−1 (0) =1 q x x − K (q x ,...,qx − ). (4.11) − 1 ··· m 1 · m 1 m 1 In the same way, iterated use of eq. (4.8) gives

(a) Km (x1,...,xa+1, q xa+2,...,qxm−1) (0) −1 −1 = Km (q x1,...,q xa, xa+1, q xa+2,...,qxm−1) (0) −1 −1 + q x K (q x ,...,q x − , x , q x ,...,qx − ) a · m 1 a 1 a a+1 m 1 2 (0) −1 −1 + q x − x K (q x ,...,q x − , x − , q x ,...,qx − ) a 1 a · m 1 a 2 a 1 a m 1 a (0) + + q x x K (x , q x ,...,qx − ). (4.12) ··· 1 ··· a · m 1 2 m 1 Substituting above equation for a = m 2 into eq. (4.9), we ﬁnd − (0) −1 −1 Km (q x1,...,q xm−2, xm−1) (0) −1 −1 + q x − K (q x ,...,q x − , x − , q x − ) m 2 · m 1 m 3 m 2 m 1 2 (0) −1 −1 + q x − x − K (q x ,...,q x − , x − , q x − , q x − ) m 3 m 2 · m 1 m 4 m 3 m 2 m 1 m−2 (0) + + q x x − K (x , q x ,...,qx − ) ··· 1 ··· m 2 · m 1 2 m 1 1 (0) −1 −1 = 1 Km (q x1,...,q xm−1) . (4.13) xm−1 −

Combining two equations (4.11) and (4.13) with x x = = x − , we obtain ≡ 1 ··· m 1 K(0)(x)=1 q x + q2m−1 xm K(0)(q2 x). m − · m This proves eq. (4.2) for a =0. 12 Other cases can be shown by using eqs. (4.11) and (4.13) with eqs. (4.10) and (4.12).

(a) (a) Corollary 7. We deﬁne the q-series Φm (τ) by Km (x = 1) up to constant, i.e.,

∞ 2 (m−1−a) 1 (a) e cm−1 cm−1(cm−1+1) Φ (τ)= m q 4m ( 1) q 2 m − c ,...,cm− =0 1 X 1 e m−2 c 2+···+c 2 +ca +···+cm− ci+1 + δi,a q 1 m−2 +1 2 , (4.14) × · ci i=1 ! Y where m 2 and a =0, 1,...,m 2, and we mean ≥ − q = e2πiτ .

Then we have ∞ (a) 1 2 (a) 4m n Φm (τ)= m χ2m(n) q . (4.15) n=0 X e Factor m in above deﬁnition is merely for our later convention.

(a) Corollary 8. Let the q-series Φm (τ) be deﬁned by

(a) 1 2 (a) 4m n Φm (τ)= n χ2m(n) q , (4.16) ∈Z Xn for m 2 and a =0, 1,...,m 2. Then we have ≥ − (a) Φm (τ) ∞ 2 (m−1−a) 1 m 1 a cm−1 cm−1(cm−1+1) =4 q 4m c + c + + c − + − − ( 1) q 2 1 2 ··· m 1 2 − c ,...,cm− =0 1 X 1 m−2 2 ··· 2 a ··· m− ci+1 + δi,a qc1 + +cm−2+c +1+ +c 2 . (4.17) × · ci i=1 ! Y Proof. We differentiate eq. (4.2) with respect to x and substitute x 1. →

5. MODULAR PROPERTY

(a) (a) We shall reveal the modular property of the q-series Φm (τ) and Φm (τ) deﬁned by eq. (4.16) (a) and eq. (4.15) respectively. The theta series Φm (τ) have weight 3/2, and span (m 1)-dimensional e − space; it is straightforward to get

a2 (m−1−a) πi (m−1−a) Φ (τ +1)=e 2m Φ (τ), (5.1) m · m 13 and from the standard method using the Poisson summation formula we have i 3/2 1 Φ (τ)= M Φ , (5.2) m τ m · m − τ where (m−2) Φm (τ) . . Φm(τ)=   , Φ(1)(τ)  m   (0) τ   Φm ( )  and M is an (m 1) (m 1) matrix,   m − × − 2 a b M = sin π . (5.3) m 1≤a,b≤m−1 m m r Remarkable is that the theta series defined by (m−λ) m Φm+2 (τ) chλ (τ)= , (5.4) 2 η(τ) 3 where η(τ) is the Dedekind η-function (6.3), is the affinesu(2)m character (see Examples in Sec- tion 6) [4,11]. c As studied in Ref. 9, we have interests in the Eichler integral of the modular form Φm(τ). Generally when the q-series ∞ n F (τ)= an q , n=1 X is a modular form with weight k Z≥ , the Eichler integral defined as k 1 integrations of F (τ) ∈ 2 − with respect to τ, or explicitly defined by ∞ a F (τ)= n qn, nk−1 n=1 X satisfies e (c τ + d)k−2 F (γ(τ)) F (τ)= G (τ), (5.5) · − γ a b where γ = SL(2; Z), and G (τ) is the period polynomial c d ∈ γ e e k−1 i∞ (2 π i) k−2 Gγ (z)= F (τ)(z τ) dτ. (k 2)! −1 − − Zγ (i ∞) (a) In our case, the modular form Φm (τ) in eq. (4.16) has a half-integral weight, and above story does not work any more. But the Eichler integral as an infinite q-series can be defined in a naive sense, and we may find that ∞ (a) 1 2 (a) 4m n Φm (τ)= m χ2m(n) q , n=0 X14 e which is nothing but a definition (4.15). We note that a prefactor in above definition is for our (a) (a) convention. We can regard Φm (τ) as the Eichler integral of the modular form Φm (τ) with weight 3/2. e To study a nearly modular property (see eq. (5.5)) of this Eichler integral of the half-integral weight modular form, we first recall a following result.

Proposition 9. Let Cf (n) be a periodic function with mean value 0 and modulus f. Then we have an asymptotic expansion as t 0 ց ∞ ∞ k 2 ( t) C (n) e−n t L( 2 k,C ) − , f ≃ − f k! n=1 X Xk=0 where L(k,Cf ) is the L-function associated with Cf (n), and is given by f f k n L( k,C )= C (n) B . − f −k +1 f k+1 f n=1 X Proof. It is a standard result using the Mellin transformation. See, e.g., Ref. 9.

(a) Using this property, we obtain the Eichler integral Φm (τ) near at a root of unity as follows.

Proposition 10. The Eichler integral (4.15) for τ =e M Q (N > 0, and M, N are coprime N ∈ integers) reduces to mN n2 (a) M (a) n Mπi Φ = m χ (n) 1 e 2mN . (5.6) m N 2m − m N n=0 X e Proof. We have from eq. (4.15) ∞ n2 (a) M t − t Φ + i = m C (n) e 4m , m N 2 π 2mN n=0 X where e 2 (a) Mn 2mN πi C2mN (n)= χ2m(n) e . We see that C (n +2 m N) = C (n) and C (2 m N n) = C (n), and we can 2mN 2mN 2mN − − 2mN apply Prop. 9 to get an asymptotic expansion in t 0 as ∞ ց M t L( 2 k,C ) t k Φ(a) + i − 2mN . m N 2 π ≃ k! −4 m Xk=0 Then we obtain a limitinge value M Φ(a) = m L(0,C ). m N · 2mN e 15 Recalling an explicit form of the Bernoulli polynomial B (x)= x 1 and a property C (2mN 1 − 2 2mN − n)= C (n), we obtain eq. (5.6). − 2mN

(a) Conjecture 2. Let Φm (τ) be the Eichler integral deﬁned by eq. (4.15) or (4.14). When q is the N- (a) th root of unity, Φm (1/N), which was computed as eq. (5.6), coincides with an expression (3.10) e up to constant, i.e., e m− −a 2 (a) 1 ( 1 ) πi (a) Φ = e 2mN Y (ω). (5.7) m N · m e Proof for a case of a =0. As a case of a =0 is related to Kashaev’s invariant for the torus link as in eq. (3.11), this case can be directly proved by using eq. (3.18) as follows;

m− 2 m− 2 ( 1) πi (0) 1 ( 1) πi e 2mN Y (ω)= e 2mN T (2, 2 m) m N h iN 2mN k2 1 2 (0) πi = k χ (k) e 2mN −4 m N 2 2m Xk=0 2mN k2 (0) m k πi = χ (k) e 2mN 2m 2 − 2 N Xk=0 mN k2 (0) k πi (0) 1 = χ (k) m e 2mN = Φ . 2m − N m N Xk=0 e In the third equality, we have summed an expression with k 2 m N k. As a result of eq. (5.6) → − the statement of Conjecture is true for a =0.

We now discuss how Conjectures 1 follows from Conjecture 2. We ﬁrst recall from eq. (5.6) that for N Z ∈ (m−1−a)2 (a) 2m πiN Φm (N)=(1+ a) e . (5.8)

Following Ref. 9 (see also Ref. 21)e we deﬁne the period function

m ∞ Φ(a)(τ) r(a)(z; α)= m dτ, (5.9) m 8 i √τ z r Zα − where α Q. It is defined for z in the lower half plane, z H−, but it is analytically continued to ∈ ∈ − (a) R = ∂H . To see a modular property of Φm (α + i y) in y 0, we further define ց ∞ (a) (a) e m Φm (τ) Φm (z)= dτ, (5.10) 8 i z∗ √τ z r 16Z − b where z H−. We can find that this function is nearly modular of weight 1/2 by ∈ m−1 m−1 0 (m−1−b) (m−1−b) 1 m Φm ( 1/s) ds Mm Φm = Mm − a,b − z a,b 8 i ∗ √ s−1 + z−1 s2 b=1 b=1 r Zz X X − b = √ i z Φ(m−1−a)(z) r(m−1−a)(z;0) . (5.11) m − m On the other hand, we have for z = α + i y b

n2 ∞ m πiτ (a) m (a) e 2 Φm (z)= n χ2m(n) dτ 8 i ∗ √τ z r n∈Z Zz − ∞ X b n2 (a) πiz πy = m χ (n) e 2m erfc n , 4m − m n=1 r X where erfc(x) is the complementary error function ∞ 2 2 erfc(x)= e−t dt. √π Zx (a) As erfc(0) = 1 and we know eq. (5.6) from a deﬁnition of Φm (α), we get

(a) (a) Φm (α)= Φm (α),e (5.12)

− for α Q. We stress that LHS is a limitb from thee lower half plane H while RHS is analytically ∈ continued from the upper half plane H.

Taking a limit z 1/N for N Z in eq. (5.11), we obtain → ∈ − 1 m 1 1 Φ(m−1−a) = √ i N M Φ(m−1−b) ( N)+ r(m−1−a) ;0 . (5.13) m N − m a,b m − m N b=1 X b b (a) 1 When we recall eq. (5.8) and use an asymptotic expansion of rm ;0 in N , we obtain a N → ∞ following proposition.

(a) Proposition 11. Let the Eichler integral Φm (τ) be deﬁned by eq. (4.15). Asymptotic expansion in N is then given by →∞ e m−1 b2 (a) 1 2 (m 1 a) b − πiN Φ √ i N (m b) sin − − π e 2m m N ≃ − m − m b=1 r X ∞ e E(m;a) πi k + k . (5.14) k! 2 m N Xk=0 Based on this proposition, we get a conjecture (3.12) assuming a conjecture (5.7). 17 6. EXAMPLES

6.1 (2, 4)-Torus Link: m = 2

q-series: • The modular form with weight 3/2 is

(0) (0) 1 n2 n 1 (n2+n+ 1 ) Φ (τ)= n χ (n) q 8 = ( 1) (2 n + 1) q 2 4 (6.1) 2 4 − ∈Z ∈Z Xn Xn 1 3 6 10 15 =2 q 8 1 3 q +5 q 7 q +9 q 11 q + , − − − ··· (0) where χ4 (n) is the primitive character modulo 4, n mod 4 1 3 others χ(0)(n) 1 1 0 4 − It is known by Jacobi that we can write

(0) 3 Φ2 (τ)=2 η(τ) , (6.2)

where η(τ) is the Dedekind η-function,

1/24 η(τ)= q (q)∞. (6.3) · From a modular property of the η-function we see that

(0) 1 (0) 4 πi Φ2 (τ +1)=e Φ2 (τ), (6.4) τ 3/2 Φ(0)( 1/τ)= Φ(0)(τ). (6.5) 2 − i 2 The Eichler integral is given by ∞ ∞ (0) (0) 1 n2 1 k 1 k(k+1) Φ (τ)=2 χ (n) q 8 =2 q 8 ( 1) q 2 (6.6) 2 4 − n=0 X Xk=0 1 3 6 10 15 e =2 q 8 1 q + q q + q q + . − − − ··· root of unity: • As a limit of q being the N-th root of unity, the Eichler integral (6.6) coincides with Kashaev’s invariant for torus link,

N−1 π π (0) 1 i c 1 c(c+1) 1 i Φ = e 4N ( 1) ω 2 = e 4N T (2, 4) . (6.7) 2 N − N h iN c=0 X We see thate the Eichler integral (6.6) with q being the N-th root of unity takes a same form with the original Eichler integral up to constant, only an inﬁnite sum reduces to a ﬁnite 18 (0) πiN/4 sum. Using Φ2 (N) = e , we have the nearly modular property (see eq. (3.9)),

∞ (2;0) n e (0) 1 (0) En π i Φ √ i N Φ ( N)+ , (6.8) 2 N ≃ − 2 − n! 4 N n=0 X (2;0) e e where En is the Euler number deﬁned by

24n+1 1 3 E(2;0) = B B , n −2 n +1 2n+1 4 − 2n+1 4 some of which are given as

∞ 1 E(2;0) = n x2n ch(x) (2 n)! n=0 X 1 5 61 =1 x2 + x4 x6 + . − 2 24 − 720 ···

6.2 (2, 6)-Torus Link: m = 3

q-series: • A set of the theta series is given by

(0) (0) 1 2 12 n Φ3 (τ)= n χ6 (n) q (6.9a) ∈Z Xn 1 5 8 16 21 =4 q 3 1 2 q +4 q 5 q +7 q 8 q + , − − − ··· (1) (1) 1 2 12 n Φ3 (τ)= n χ6 (n) q (6.9b) n∈Z X 1 2 4 10 14 24 =2 q 12 1 5 q +7 q 11 q + 13 q 17 q + , − − − ··· where n mod 6 2 4 others n mod 6 1 5 others χ(0)(n) 1 1 0 χ(1)(n) 1 1 0 6 − 6 − We note [11] that these series can be written in terms of the Dedekind η-function as

5 2 (0) 1 πi η(2 τ) η(τ) η(4 τ) Φ (τ)=4e 12 =4 , (6.10a) 3 1 2 η(2 τ) η(τ + 2) 5 (1) η(2 τ) Φ3 (τ)=2 2 . (6.10b) η(4 τ) 19 The modular property is written as

(1) 1 πi (1) Φ3 (τ + 1) e 6 0 Φ3 (τ) = , (6.11) (0) 2 (0) 3 πi Φ3 (τ + 1)! 0 e ! Φ3 (τ)!

(1) (1) Φ3 ( 1/τ) τ 3/2 1 1 1 Φ3 (τ) − = . (6.12) Φ(0)( 1/τ)! i · √2 1 1! · Φ(0)(τ)! 3 − − 3 The Eichler integrals are then deﬁned by

∞ a (0) 1 a 1 a(a+1) b(b+1) a Φ (τ)=3 q 3 ( 1) q 2 q 3 − b a=0 b=0 ∞ X X e (0) 1 2 12 n =3 χ6 (n) q (6.13a) n=0 X 1 5 8 16 21 =3 q 3 1 q + q q + q q + , − − − ··· ∞ a+1 (1) 1 a 1 a(a+1) b2 a +1 Φ (τ)=3 q 12 ( 1) q 2 q 3 − b a=0 b=0 ∞ X X e (1) 1 2 12 n =3 χ6 (n) q (6.13b) n=0 X 1 2 4 10 14 24 =3 q 12 1 q + q q + q q + . − − − ··· We note that Zagier’s identity [19] leads us to ﬁnd

∞ (1) 1 n 2 Φ (τ)=3 q 12 ( 1) ( 1; q ) . 3 − − n+1 n=0 X root of unity: e • From Prop. 10, the Eichler integrals (6.13) reduce in a case of q being N-th root of unity to 3N n2 (a) 1 (a) n πi Φ =3 χ (n) 1 e 6N . (6.14) 3 N 6 − 3 N n=0 X This coincides withe

N−1 π π (0) 1 2 i a 1 a(a+1)+b(b+1) a 1 2 i Φ = e 3N ( 1) ω 2 = e 3N T (2, 6) , (6.15a) 3 N − b N h iN a,bX=0 e N−1 π (1) 1 i a 1 a(a+1)+b2 a +1 Φ = e 6N ( 1) ω 2 . (6.15b) 3 N − b a,bX=0 e 20 where the second identity remains to be proved (the ﬁrst identity was established from an asymptotic expansion due to a discussion in the previous section). We have checked numerically a validity of this identity. As we have

(1) 1 πiN Φ3 (N) 2 e 6 = , (0) 2 πiN Φ (N)! e 3 ! e3 e modular property and eqs. (6.15) supports that

(1) 1 (1) ∞ (3:1) n Φ3 1 1 1 Φ3 ( N) 1 En π i N √ i N − + , (6.16) (0) 1 ≃ − · √ · (0) n! (3:0) 6 N Φ ! 2 1 1! Φ ( N)! n=0 En ! e3 N − e3 − X e e where the generalized Euler number is deﬁned by

3 62n 2 a 4+ a E(3;a) = · B − B , n −2 n +1 2n+1 6 − 2n+1 6

or some of them are as follows;

∞ (3;1) 3 sh(2 x) E x2k = k sh(3 x) sh(x) ! E(3;0)! (2 k)! Xk=0 k 2 10 x2 34 x4 910 x6 = 3 + + . − 8 2 24 − 720 ··· 1! 3 ! 32! 896!

6.3 (2, 8)-Torus Link: m = 4

q-series: • 21 A set of theta series is given by

(0) (0) 1 2 16 n Φ4 (τ)= n χ8 (n) q (6.17a) ∈Z Xn 3 3 η(τ) τ 3 = τ η η η(2 τ)! − 2 2 9 7 10 22 27 =2 q 16 3 5 q + 11 q 13 q + 19 q 21 q + , − − − ··· (1) (1) 1 2 16 n Φ4 (τ)= n χ8 (n) q (6.17b) ∈Z Xn =4 η(2 τ) 3

1 2 6 12 20 30 =4 q 4 1 3 q +5 q 7 q +9 q 11 q + , − − − ··· (2) (2) 1 2 16 n Φ4 (τ)= n χ8 (n) q (6.17c) ∈Z Xn 3 3 η(τ) τ 3 = τ + η η( 2 ) η(2 τ)! 2 1 3 5 14 18 33 =2 q 16 1 7 q +9 q 15 q + 17 q 23 q + , − − − ···

which are modular coinvariant;

(2) 1 πi (2) Φ4 (τ + 1) e 8 0 0 Φ4 (τ) (1) 1 πi (1) Φ4 (τ + 1) =  0 e 2 0  Φ4 (τ) , (6.18) (0) 9 πi (0) Φ (τ + 1)  0 0 e 8  Φ (τ)  4     4        (2) (2) Φ4 ( 1/τ) 1 √2 1 Φ4 (τ) − τ 3/2 1 Φ(1)( 1/τ) = √2 0 √2 Φ(1)(τ) . (6.19) 4 − i · 2 − · 4 Φ(0)( 1/τ)  1 √2 1  Φ(0)(τ)  4 −   −   4    22    We then have the Eichler integral as

∞ a b (0) 9 a 1 a(a+1) b(b+1) a c(c+1) b Φ (τ)=4 q 16 ( 1) q 2 q q 4 − b c a=0 b=0 c=0 ∞ X X X e (0) 1 2 16 n =4 χ8 (n) q (6.20a) n=0 X 9 7 10 22 27 =4 q 16 1 q + q q + q q + , − − − ··· ∞ a b+1 (1) 1 a 1 a(a+1) b(b+1) a c2 b +1 Φ (τ)=4 q 4 ( 1) q 2 q q 4 − b c a=0 b=0 c=0 ∞ X X X e (1) 1 2 16 n =4 χ8 (n) q (6.20b) n=0 X 1 2 6 12 20 30 =4 q 4 1 q + q q + q q + , − − − ··· ∞ a+1 b (2) 1 a 1 a(a+1) b2 a +1 c2 b Φ (τ)=4 q 16 ( 1) q 2 q q 4 − b c a=0 b=0 c=0 ∞ X X X e (2) 1 2 16 n =4 χ8 (n) q (6.20c) n=0 X 1 3 5 14 18 33 =4 q 16 1 q + q q + q q + . − − − ··· One sees that (1) (0) Φ4 (τ)=2 Φ2 (2 τ). root of unity: • e e Limiting value of the Eichler integrals when q goes to the N-th primitive root of unity is given by 4N n2 (a) 1 (a) n πi Φ =4 χ (n) 1 e 8N , (6.21) 4 N 8 − 4 N n=0 X which is rewrittene as N−1 π π (0) 1 9 i a 1 a(a+1)+b(b+1)+c(c+1) a b 1 9 i Φ = e 8N ( 1) ω 2 = e 8N T (2, 8) , (6.22a) 4 N − b c N h iN a,b,cX=0 e N−1 π (1) 1 i a 1 a(a+1)+b(b+1)+c2 a b +1 Φ = e 2N ( 1) ω 2 , (6.22b) 4 N − b c a,b,cX=0 e N−1 π (2) 1 i a 1 a(a+1)+b2 +c2 a +1 b Φ = e 8N ( 1) ω 2 . (6.22c) 4 N − b c a,b,cX=0 e 23 Though we have checked these three equalities numerically, we proved only eq. (6.22a) in this article. The nearly modular properties are written as

(2) 1 (2) Φ4 N 1 √2 1 Φ4 ( N) 1 − Φ(1) 1  √ i N √2 0 √2 Φ(1)( N) e4 N ≃ − · 2 − · e4 − Φ(0) 1   1 √2 1  Φ(0)( N) e4 N   −  e4 −        (4:2) En e e ∞ n 1 (4:1) π i + En  , (6.23) n! 8 N n=0 X E(4:0)  n    where we have (2) 1 πiN Φ4 (N) 3 e 8 (1) 1 πiN Φ (N) = 2 e 2  , e4 (0) 9 πiN Φ (N)  e 8  e4        and the generalized Euler number is deﬁned by e 26n+2 3 a 5+ a E(4;a) = B − B . n −2 n +1 2n+1 8 − 2n+1 8 Some of them are explicitly given as follows;

E(4;2) sh(3 x) ∞ k 4 x2k sh(2 x) = E(4;1) sh(4 x)   k (2 k)! sh(x) Xk=0 (4;0)   Ek      3  7 119 5587 x2 x4 x6 = 2 8 + 160 7808 + .   −   2   24 −   720 ··· 1 5 109 5465                

ACKNOWLEDGMENTS

The author would like to thank Anatol N. Kirillov for suggesting to study torus links. He also thanks Hitoshi Murakami for communications in constructing quantum knot invariants. Thanks are to Don Zagier for comments on nearly modular form. He thanks George Andrews for his interests. This work is supported in part by the Sumitomo Foundation, and Grant-in-Aid for Young Scientists from the Ministry of Education, Culture, Sports, Science and Technology of Japan. 24 REFERENCES

[1] G.E.ANDREWS, The Theory of Partitions, Addison-Wesley, London, 1976. [2] K. HIKAMI, Volume conjecture and asymptotic expansion of q-series, Exp. Math. (2003). to appear. [3] , q-series and L-functions related to half-derivatives of the Andrews–Gordon identity, Ramanujan J., (2003). to appear. [4] V. G. KAC, Infinite Dimensional Lie Algebras, Cambridge Univ. Press, Cambridge, 3 ed., 1990. [5] R.M.KASHAEV, A link invariant from quantum dilogarithm, Mod. Phys. Lett. A, 10 (1995), pp. 1409–1418. [6] , The hyperbolic volume of knots from quantum dilogarithm, Lett. Math. Phys., 39 (1997), pp. 269–275. [7] R.M.KASHAEV AND O. TIRKKONEN, Proof of the volume conjecture for torus knots, Zap. Nauch. Sem. POMI, 269 (2000), pp. 262–268. [8] R. LAWRENCE AND L. ROZANSKY, Witten–Reshetikhin–Turaev invariants of Seifert manifolds, Commun. Math. Phys., 205 (1999), pp. 287–314. [9] R. LAWRENCE AND D. ZAGIER, Modular forms and quantum invariants of 3-manifolds, Asian J. Math., 3 (1999), pp. 93–107. [10] R. LICKORISH, An Introduction to Knot Theory, Springer, New York, 1997. [11] I. G. MACDONALD, Affine root systems and Dedekind’s η-function, Invent. Math., 15 (1972), pp. 91–143. [12] H.R.MORTON, The coloured Jones function and Alexander polynomial for torus knots, Proc. Cambridge Philos. Soc., 117 (1995), pp. 129–135. [13] H. MURAKAMI, Kashaev’s invariant and the volume of a hyperbolic knot after Yokota, in “Physics and Combi- natorics”, A. N. Kirillov, A. Tsuchiya, and H. Umemura, eds., World Scientific, Singapore, 2001, pp. 244–272. [14] H. MURAKAMI AND J. MURAKAMI, The colored Jones polynomials and the simplicial volume of a knot, Acta Math., 186 (2001), pp. 85–104. [15] M. ROSSO AND V. JONES, On the invariants of torus knots derived from quantum groups, J. Knot Theory Ramifications, 2 (1993), pp. 97–112. [16] L. ROZANSKY, Higher order terms in the Melvin–Morton expansion of the colored Jones polynomial, Commun. Math. Phys., 183 (1997). [17] A. STOIMENOW, Enumeration of chord diagrams and an upper bound for Vassiliev invariants, J. Knot Theory Ramifications, 7 (1998), pp. 93–114. [18] Y. YOKOTA, On the volume conjecture for hyperbolic knots, math.QA/0009165, (2000). [19] D.ZAGIER, Vassiliev invariants and a strange identity related to the Dedekind eta-function, Topology, 40 (2001), pp. 945–960. [20] S. P. ZWEGERS, Mock ϑ-functions and real analytic modular forms, in“q-Series with Applications to Combi- natorics, Number Theory, and Physics”, B. C. Berndt and K. Ono, eds., Amer. Math. Soc., Providence, 2001, pp. 269–277. [21] , Mock Theta Functions, PhD thesis, Universiteit Utrecht, 2002.

DEPARTMENT OF PHYSICS, GRADUATE SCHOOL OF SCIENCE, UNIVERSITY OF TOKYO, HONGO 7–3–1, BUNKYO, TOKYO 113–0033, JAPAN. E-mail address: [email protected]