Arxiv:2010.16371V2 [Math.NT] 12 Jul 2021 Kronecker Limit Formulas, Which Correspond to the Positive Deﬁnite Case, in [3]

A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS

GENE S. KOPP

Abstract. We prove an analogue of Kronecker’s second limit formula for a continuous family of “indefinite zeta functions”. Indefinite zeta functions were introduced in the author’s previous paper as Mellin transforms of indefinite theta functions, as defined by Zwegers. Our formula is valid in dimension g = 2 at s = 1 or s = 0. For a choice of parameters obeying a certain symmetry, an indefinite zeta function is a differenced ray class zeta function of a real quadratic field, and its special value at s = 0 was conjectured by Stark to be a logarithm of an algebraic unit. Our formula also permits practical high-precision computation of Stark ray class invariants.

1. Introduction In a previous paper [5], we introduced indefinite zeta functions as Mellin transforms of (1) certain indefinite theta functions associated to the intermediate Siegel half-space Hg , defined below. In this paper, we obtain a formula for the values of such an indefinite zeta function at s = 1 or s = 0, in the special case of dimension g = 2. Such formulas are traditionally called Kronecker limit formulas, after Kronecker’s first and second limit formulas giving the constant term in the Laurent expansion at s = 1 of standard and twisted real analytic Eisenstein series. When our parameters are specialised appropriately, our special value is a finite linear combination of Hecke L-values at s = 1. Our formula may be used to compute values of Hecke L-functions at s = 1 (resp. s = 0) relevant to the Stark conjectures, which we discuss in Section 1.4. For imaginary quadratic fields, Stark proved his conjectures using Kronecker’s first and second limit formulas together with the theory of singular moduli [14]. The Kronecker limit formulas give the constant Laurent series coefficient at s = 1 for families of Dirichlet series continuously interpolating the ray class zeta functions ζ(s, A)—namely, standard and twisted real analytic Eisenstein series (see [13] for details). We will discuss a generalisation of these arXiv:2010.16371v2 [math.NT] 12 Jul 2021 Kronecker limit formulas, which correspond to the positive definite case, in [3]. Kronecker limit formulas applicable to real quadratic fields were developed by Hecke, Her- glotz, and Zagier (in analogy with the first Kronecker limit formula), and by Shintani (in analogy with the second Kronecker limit formula). As in the imaginary quadratic case, these formulas are obtained by continuously interpolating between ray class zeta functions using a larger family of functions. Hecke’s formula uses cycle integrals of real analytic Eisenstein series, whereas the formulas of Herglotz [2] and Zagier [21] (see also [1, 8]) and the formulas of Shintani [10, 11, 12] use partial zeta functions defined by summing over a cone. Analogues of

Date: July 13, 2021. Key words and phrases. Kronecker limit formula, real quadratic field, indefinite quadratic form, indefinite theta function, Epstein zeta function, Stark conjectures. 1 2 GENE S. KOPP the Kronecker limit formulas in other settings have been found by Liu and Masri [6], Posingies [7], and Vlasenko and Zagier [20], among others. The main theorem of this paper supplies a new real quadratic analogue of Kronecker’s second limit formula based on a new interpolation between ray class zeta functions. The interpolation is by the indefinite zeta functions introduced in [5], which are Mellin transforms of nonholomorphic indefinite theta functions. Indefinite zeta functions have a nice functional equation, but they do not have a Dirichlet series representation for general parameters. The main results on indefinite zeta functions—stated in Section 1.3—require a lot of notation, defined in Section 1.1 and Section 1.2. The proofs of the indefinite Kronecker limit formulas are provided in Section 2.

1.1. Notational conventions. We list some notational conventions used in the paper. • e(z) := exp(2πiz) is the complex exponential, and this notation is used for z ∈ C not necessarily real. •H := {τ : Im τ > 0} is the complex upper half-plane. • Non-transposed vectors v ∈ Cg are always column vectors; the transpose v> is a row vector. • If M is a g × g matrix, then M > is its transpose, and (when M is invertible) M −> is a shorthand for (M −1)>. 1 > • QM (v) denotes the quadratic form QM (v) := 2 v Mv, where M is a g × g matrix, and v is a g × 1 column vector. • f(c)|c2 := f(c ) − f(c ), where f is any function taking values in an additive group. c=c1 2 1 v1 2 2 v1 • If v = ( v2 ) ∈ C and f is a function of C , we may write f(v) as f( v2 ) rather than v1 f(( v2 )). • We often express Ω = iM + N where M,N are real g × g symmetric matrices; N and M will always have real entries even when we do not say so explicitly. We use complex logarithms throughout this paper. If f(τ) is any nonvanishing holomorphic function on the upper half plane H, there is some holomorphic function (Log f)(τ) such that exp ((Log f)(τ)) = f(τ), because H is simply connected. Specifying a single value (or the limit as τ approaches some element of R ∪ {∞}) specifies Log f uniquely. It won’t necessarily be true that (Log f)(τ) = log(f(τ)). Conventions for square roots, when not specified, follow [5]. See Section 2.3 and Section 3.2 therein for details. We recall the definition of the Siegel intermediate half-space, as defined in [5]. Definition 1.1. For 0 ≤ k ≤ g, we define the Siegel intermediate half-space of genus g and index k to be

(k) > Hg := {Ω ∈ Mg(C) : Ω = Ω and Im(Ω) has signature (g − k, k)}. (1.1)

(k) The Hg are the open orbits of the action of Sp2g(R) by fractional linear transformations (0) on the space of complex symmetric matrices. In particular, Hg is the usual Siegel upper half-space.

1.2. Indeﬁnite theta and zeta functions. We review the relevant deﬁnitions from [5]. A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 3

Definition 1.2. For any complex number α, define the function Z α E(α) := e−πu2 du, (1.2) 0 where the integral runs along any contour from 0 to α. Definition 1.3. Let Ω = iM + N be a complex symmetric matrix whose imaginary part (1) has signature (g − 1, 1); that is, Ω ∈ Hg . Define the (nonholomorphic) indefinite theta function   c2 c> Im(Ωn + z) 1 c1,c2 X > > Θ (z, Ω) := Eq  e n Ωn + n z , (1.3) g 1 > 2 n∈Z − 2 c Im(Ω)c c=c1 g g > > where z ∈ C , c1, c2 ∈ C , c1 Mc1 < 0, and c2 Mc2 < 0. Nonholomorphic indefinite theta functions were first studied by Vignéras[18, 19] and were rediscovered by Zwegers [22]. Zwegers’s theta function is defined for real cj when N is a scalar multiple of M. More precisely, if M is real symmetric matrix of signature (g − 1, 1), τ ∈ H, g c1,c2 and c1, c2 ∈ R , then Θ (Mz, τM) is equal up to an exponential factor to the function c1,c2 ϑM (z, τ) introduced by Zwegers on page 27 of [22]. Our theta functions extend Zwegers’s to the Siegel modular setting; a related generalisation has also been studied by Roehrig [9].

(1) Definition 1.4. Let Ω = iM + N ∈ Hg . Define the indefinite theta null with characteristics p, q ∈ Rg: 1 Θc1,c2 (Ω) := e q>Ωq + p>q Θc1,c2 (p + Ωq; Ω) . (1.4) p,q 2

g > > where c1, c2 ∈ C , c1 Mc1 < 0, and c2 Mc2 < 0. We define the indefinite zeta function using a Mellin transform of the indefinite theta function with characteristics.

(1) Definition 1.5. Let Ω = iM + N ∈ Hg . The completed indefinite zeta function is Z ∞ dt c1,c2 c1,c2 s ζbp,q (Ω, s) := Θp,q (tΩ)t , (1.5) 0 t g g > > where p, q ∈ R , and c1, c2 ∈ C are parameters satisfying c1 Mc1 < 0 and c2 Mc2 < 0. The completed indefinite zeta function has an analytic continuation and satisfies a functional equation, which is Theorem 1.1 of [5].

c1,c2 Theorem 1.6 (Analytic continuation and functional equation for ζbp,q (Ω, s)). The function c1,c2 ζbp,q (Ω, s) may be analytically continued to an entire function on C. It satisﬁes the functional equation g e(p>q) c1,c2 Ωc1,Ωc2 −1 ζb Ω, − s = ζb−q,p −Ω , s . (1.6) p,q 2 pdet(−iΩ) 4 GENE S. KOPP

1.3. Kronecker limit formulas for indefinite zeta functions. The Kronecker limit formula for indefinite zeta functions involves the dilogarithm function and a rapidly convergent integral of a logarithm of an infinite product. We also require the following definition of the c function κΩ(v), which is the square root of a rational function and will appear as a factor in the integrand.

(1) 2 > 2 Definition 1.7. Suppose Ω = iM + N ∈ H2 , c ∈ C satisfying c Mc < 0, v ∈ C , and s ∈ . Let Λc := Ω − i Mcc>M. Then, we define C Ω QM (c) c>Mv κc (v) := . (1.7) Ω p q 4πi −Q (c)Q (v) −2iQ c (v) M Ω ΛΩ We now state the formula. (1) Theorem 1.8 (Indefinite Kronecker limit formula at s = 1). Let Ω = iM + N ∈ H2 , p1 2 2 2 > p = ( p2 ) ∈ R \ Z , and c1, c2 ∈ C such that cj Mcj < 0. For c = c1, c2, factor the quadratic form ξ + − QΛc = α(c)(ξ − τ (c))(ξ − τ (c)), (1.8) Ω 1 where τ +(c) is in the upper half-plane and τ −(c) is in the lower half-plane. Then,

c1,c2 + − + − ζbp,0 (Ω, 1) = I (c2) − I (c2) − I (c1) + I (c1), (1.9) where 1 I±(c) := − Li (e(±p ))κc 2 1 Ω 0 Z ∞ ± ± c ± (τ (c) + it) + 2i (Log ϕp1,±p2 )(±τ (c) + it)κΩ dt. (1.10) 0 1

The function ϕp1,p2 : H → C is deﬁned by the a product expansion, ∞ Y 1 − e ((d + p1)ξ + p2) ϕp1,p2 (ξ) := (1 − e(p1ξt + p2)) , (1.11) 1 − e ((d − p1)ξ − p2) d=1 and its logarithm (Log ϕp1,p2 )(ξ) is the unique continuous branch with the property log(1 − e(p2)) if p1 = 0, lim (Log ϕp1,p2 )(ξ) = (1.12) ξ→i∞ 0 if p1 6= 0.

Here log(1 − e(p2)) is the standard principal branch. The following specialisation looks somewhat simpler and contains all of the cases of arithmetic zeta functions ZA(s) associated to real quadratic ﬁelds. Theorem 1.9 (Indeﬁnite Kronecker limit formula at s = 1, pure imaginary case). Let M be p1 2 2 2 a 2×2 real matrix of signature (1, 1), and let Ω = iM. Let p = ( p2 ) ∈ R \Z , and c1, c2 ∈ R > such that cj Mcj < 0. Then,

c1,c2 ζbp,0 (Ω, 1) = 2i Im (I(c2) − I(c1)) , (1.13) A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 5 where 1 I(c) = − Li (e(p ))κc 2 1 Ω 0 Z ∞ c τ(c) + it + 2i (Log ϕp1,p2 )(τ(c) + it)κΩ dt. (1.14) 0 1 c Here, Log ϕp1,p2 and κΩ are deﬁned as in the statement of Theorem 1.8, and ξ = τ(c) is the ξ unique root of the quadratic polynomial Q c in the upper half plane. ΛΩ 1 It is straightforward to use the functional equation for the indeﬁnite zeta function to c1,c2 rephrase Theorem 1.8 and Theorem 1.9 as formulas for ζb0,q (Ω, 0).

(1) Theorem 1.10 (Indeﬁnite Kronecker limit formula at s = 0). Let Ω = iM + N ∈ H2 , q1 2 2 2 > q = ( q2 ) ∈ R \ Z , and c1, c2 ∈ C such that cj Mcj < 0. For c = c1, c2, factor the quadratic form ξ + − QΛΩc = β(c)(ξ − ω (c))(ξ − ω (c)), (1.15) −Ω−1 1 where ω+(c) is in the upper half-plane and ω−(c) is in the lower half-plane. Then,

c1,c2 1 + − + − ζb (Ω, 1) = J (c2) − J (c2) − J (c1) + J (c1) , (1.16) 0,q pdet(−iΩ) where ± Ωc 1 J (c) := − Li (e(∓q ))κ −1 2 1 −Ω 0 Z ∞ ± ± Ωc ± (ω (c) + it) + 2i (Log ϕ−q1,∓q2 )(±ω (c) + it)κ−Ω−1 dt. (1.17) 0 1 Here, Log ϕ and κ are deﬁned as in the statement of Theorem 1.8. Theorem 1.11 (Indeﬁnite Kronecker limit formula at s = 0, pure imaginary case). Let M q1 2 2 be a 2 × 2 real matrix of signature (1, 1), and let Ω = iM. Let q = ( q2 ) ∈ R \ Z , and 2 > c1, c2 ∈ R such that cj Mcj < 0. Then,

c1,c2 2i ζb (Ω, 0) = Im (J(c2) − J(c1)) , (1.18) 0,q pdet(M) where Ωc 1 J(c) = − Li (e(−q ))κ −1 2 1 −Ω 0 Z ∞ Ωc ω(c) + it + 2i (Log ϕ−q1,−q2 )(ω(c) + it)κ−Ω−1 dt. (1.19) 0 1 Here, Log ϕ and κ are deﬁned as in the statement of Theorem 1.8, and ξ = ω(c) is the unique ξ root of the quadratic polynomial QΛΩc 1 in the upper half plane. −Ω−1 6 GENE S. KOPP

1.4. Application: indefinite zeta functions, real quadratic fields, and Stark units. The Hecke L-value LK (1, χ) contains arithmetic information that is not well-understood in general. The abelian Stark conjectures predict that this value is an algebraic number times a regulator Regχ, which is a determinant of a matrix of linear forms in logarithms of algebraic units in a particular abelian extension of the number field K [14, 15, 16, 17]. This conjecture is known when the base field K is equal to Q or an imaginary quadratic field, but not (for instance) when K is a real quadratic field. The rank 1 abelian Stark conjectures give a partial answer to Hilbert’s 12th Problem, which asked for explicit generators for the abelian extensions of a number field in terms of special values of transcendental functions. Computationally, the Stark conjectures are used to calculate class fields in the computer algebra systems Magma and PARI/GP. The Stark conjectures are most precisely formulated in the rank 1 case—that is, when

LK (s, χ) vanishes to order 1 at s = 0. The regulator Regχ in that case is a determinant of a 1 × 1 matrix. The Stark conjectures are most succinctly written as a statement about the ray class zeta function (of a ray ideal class A)

X −s ζ(s, A) := ζK (s, A) := N(A) (1.20) a∈A rather that as a statement about the Hecke L-function X LK (s, χ) = χ(A)ζ(s, A). (1.21) A Just as definite zeta functions specialise to ray class zeta functions of imaginary quadratic fields, indefinite zeta functions specialise to differenced ray class zeta functions of real quadratic fields. The full details of this specialisation are given in Section 7 of [5]. Definition 1.12 (Ray class zeta function). Let K be any number field and c an ideal of the maximal order OK . Let S be a subset of the real places of K (i.e., the embeddings K,→ R). Let A be a ray ideal class modulo cS, that is, an element of the group

{nonzero fractional ideals of OK coprime to c} ClcS(OK ) := . (1.22) {aOK : a ≡ 1 (mod c) and a is positive at each place in S} For Re(s) > 1, deﬁne the ray class zeta function of A to be X ζ(s, A) := N(a)−s. (1.23) a∈A This function has a simple pole at s = 1 with residue independent of A. The pole may be eliminated by considering the function ZA(s), deﬁned as follows.

Definition 1.13 (Differenced ray class zeta function). Let R be the element of ClcS(OK ) defined by

R := {aOK : a ≡ −1 (mod c) and a is positive at each place in S}. (1.24) For Re(s) > 1, deﬁne the diﬀerenced ray class zeta function of A to be

ZA(s) := ζ(s, A) − ζ(s, RA). (1.25)

The function ZA(s) extends to a holomorphic function on the whole complex plane. The 0 rank 1 abelian Stark conjecture says that ZA(0) is the logarithm of an algebraic unit. A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 7

Conjecture 1.14 (Stark [16]). Let K be a real quadratic ﬁeld and {ρ1, ρ2} the real embeddings 0 of K. If R is not the identity of Clc∞2 (OK ), then ZA(0) = log(ρ1(εA)) for an algebraic unit

εA generating the ray class field Lc∞2 corresponding to Clc∞2 (OK ). The units are compatible Art(A) with the Artin map: εid = εA. The specialisation of the indefinite zeta function to a differenced real quadratic zeta function is given by the following result, which is Theorem 1.3 of [5]. Theorem 1.15 (Specialisation of indefinite zeta function at s = 0). For each ray class −1 A ∈ Clc∞1∞2 (OK ) and integral ideal b ∈ A , there exists a real symmetric 2 × 2 matrix M, 2 2 vectors c1, c2 ∈ R , and q ∈ Q such that

−s c1,c2 (2πN(b)) Γ(s)ZA(s) = ζb0,q (iM, s). (1.26)

0 We may use Theorem 1.15 to compute presumptive Stark units exp(ZA(0)). Speciﬁcally, Corollary 1.16. Under the specialisation given by Theorem 1.15,

0 c1,c2 ZA(0) = ζb0,q (iM, 0). (1.27)

Proof. Take the limit of eq. (1.26) as s → 0. We give an example of such a computation in Section 3.

2. Proof of the Kronecker limit formulas The method of proof is to compute the Fourier series in ξ for an indefinite theta function ξ with respect to an action by a one-parameter unipotent subgroup {T } of SL2(R), then take a Mellin transform and specialise variables. After taking the Mellin transform, we must allow ξ to be a complex parameter and perform a fairly delicate contour integration. Unlike in the definite case, the Fourier coefficients of the indefinite theta are not elementary functions, which ultimately leads to a more complicated Kronecker limit formula. 2 > We fix the following notation for this section. Let c1, c2 ∈ C satisfying cj Mcj < 0, and c1,c2 2 consider the indefinite theta Θp,q with characteristics p, q ∈ R , as defined in Definition 1.4. (1) Let t > 0, Ω ∈ H2 , and M = Im(Ω). Write the indefinite theta of tΩ as

c1,c2 X c1,c2 > Θp,q (tΩ) = ρIm(tΩ)(n + q) e QΩ(n + q)t + p (n + q) (2.1) n∈Z2 X c1,c2 1/2 > = ρM (n + q)t e QΩ(n + q)t + p (n + q) , (2.2) n∈Z2 where     > > c1,c2 c2 Mv c1 Mv ρM (v) := E  − E  , (2.3) q 1 > q 1 > − 2 c2 Mc2 − 2 c1 Mc1 and Z z E(z) := e−πu2 du. (2.4) 0 8 GENE S. KOPP

2.1. Some lemmas about the Siegel upper half-space. The statement of our Kronecker c (0) limit formula, Theorem 1.8, involves a matrix ΛΩ in the Siegel upper half-space H2 . Its proof (0) will require a few basic lemmas about H2 .

ω11 ω12 (0) Lemma 2.1. Let Ω = ( ω12 ω22 ) ∈ H2 . Then −1 det Ω ω 2 Im Im > Im 12 . (2.5) ω11 ω11 ω11 Proof. Express Ω in terms of its real and imaginary parts, ω ω n n m m 11 12 = 11 12 + i 11 12 . (2.6) ω12 ω22 n12 n22 m12 m22 2 Note that m11 6= 0 because m11m22 − m12 = det M > 0, and thus ω11 6= 0. By an algebraic calculation, 2 2 −1 det Ω ω12 m11m22 − m12 Im Im − Im = 2 2 . (2.7) ω11 ω11 ω11 n11 + m11 2 2 2 Now, m11m22 − m12 = det M is positive, and so is n11 + m11. Thus, the inequality eq. (2.5) holds. We will also need the following inequality. √ ω11 ω12 (0) z −ω12+ det(−iΩ) Lemma 2.2. Let Ω = ( ω ω ) ∈ H . The two roots of QΩ ( 1 ) = 0 are τ1 = √ 12 22 2 ω11 −ω12− det(−iΩ) and τ2 = . Then, Im(τ1) > 0 > Im(τ2). ω11 z 2 Proof. We have QΩ ( 1 ) = ω11z + 2ω12z + ω22, and the expressions for the roots come from the quadratic formula. 2 2 2 For any complex numbers α = a1 + ia2 and β = b1 + ib2, (Im(αβ)) − Im(α ) Im(β ) = 2 2 2 2 (a1b2 − a2b1) ≥ 0. Thus, (Im(αβ)) ≥ Im(α √) Im(β ). 1 det(−iΩ) √ In particular, taking α = √ and β = √ (for any choice of −ω11), we obtain −ω11 −ω11 the inequality !!2 pdet(−iΩ) −1 det(−iΩ) Im ≥ Im Im (2.8) ω11 ω11 −ω11 −1 det(Ω) = Im Im . (2.9) ω11 ω11 Appealing to Lemma 2.1, we see by transitivity that !!2 pdet(−iΩ) ω 2 Im > Im 12 . (2.10) ω11 ω11 By subtracting the left-hand side and factoring, this inequality may be rewritten as 0 > Im(τ1) Im(τ2). So Im(τ1) and Im(τ2) are always nonzero real numbers with opposite signs. i 0 (0) In the special case Ω = ( 0 i ), τ1 = i and τ2 = −i. Since H2 is connected, we always have Im(τ1) > 0 > Im(τ2). A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 9

2.2. Some integrals involving E(u). We will now prove a few integral formulas that we will need.

Lemma 2.3. Suppose that α, β ∈ C satisfy Re (α2 − 2iβ) > 0. Then, using the standard branch of the square root function, Z ∞ −α E(αt1/2)e(βt) dt = . (2.11) p 2 0 4πiβ α − 2iβ Proof. By integration by parts, Z ∞ 1 Z ∞ d (e(βt)) E(αt1/2)e(βt) dt = E(αt1/2) dt (2.12) 0 2πiβ 0 dt Z ∞ 1 1/2 ∞ −πα2t α −1/2 = E(αt )e(βt) t=0 − e t e(βt) dt (2.13) 2πiβ 0 2 −α Z ∞ dt = exp (− (πα − 2πiβ) t) t1/2 (2.14) 4πiβ 0 t −α Z u 1/2 du = exp(−u) 2 (2.15) 4πiβ C πα − 2πiβ u −α Z du = e−uu1/2 , (2.16) 3/2 p 2 4π iβ α − 2iβ C u where the contour C is a ray from the origin through the point α2 − 2iβ. If z ∈ C with x = Re(z) > 0, s ∈ C with σ = Re(s) > 0, and [z1, z2] denotes the oriented line segment from z1 to z2, then Z du Z du Z du lim e−uus = lim e−uus + e−uus (2.17) N→∞ [0,Nz] u N→∞ [0,Nx] u [Nx,Nz] u Z du = Γ(s) + lim e−uus (2.18) N→∞ [Nx,Nz] u = Γ(s) + lim O e−NxN σ (2.19) N→∞ = Γ(s). (2.20)

R −u 1/2 du 1 1/2 Thus, in particular, C e u u = Γ 2 = π . Plugging this into eq. (2.16) gives eq. (2.11). As usual, let M = Im(Ω). Deﬁne the following auxiliary function, which will appear as a factor in the integral in the indeﬁnite Kronecker limit formula.

Deﬁnition 2.4. For v ∈ C2 and s ∈ C, set Z ∞ c c 1/2 s dt κΩ(v, s) := − ρM vt e (QΩ(v)t) t . (2.21) 0 t Also, set

c1,c2 c2 c1 κΩ (v, s) := κΩ (v, s) − κΩ (v, s) (2.22) Z ∞ c1,c2 1/2 s dt = ρM vt e (QΩ(v)t) t . (2.23) 0 t 10 GENE S. KOPP

c c c1,c2 c1,c2 In the case s = 1, we will leave out s and set κΩ(v) := κΩ(v, 1), κΩ (v) := κΩ (v, 1). c We have the following formula for κΛ(v) for a particular Λ that will come up in the calculation. Corollary 2.5. Let Λc := Ω − i Mcc>M. Note that Λc ∈ H(0) by Lemma 3.6 of [5]. Ω QM (c) Ω 2 Then, c>Mv κc (v) = . (2.24) Ω p q 4πi −Q (c)Q (v) −2iQ c (v) M Ω ΛΩ

Proof. Follows from Lemma 2.3. The following lemma will be needed to evaluate certain integrals.

Lemma 2.6. For any real number α ∈ R, Z ∞ c1,c2 1/2 2 s dt sgn(α) c1,c2 ρM vαt e QΩ(v)α t t = − 2s κΩ (v, s). (2.25) 0 t |α|

c1,c2 Proof. Follows from the deﬁnition of κΩ (v, s). 2.3. Fourier series of a unipotent transform of an indeﬁnite theta function. Con- sider the function of ξ ∈ R (although ξ will be allowed to be complex later on) and t ∈ R≥0, −ξ −ξ T c1,T c2 ξ> ξ h(ξ, t) := Θ > t T ΩT (2.26) (T ξ) p,T −ξq

X c1,c2 ξ 1/2 ξ > ξ = ρΩ T n + q t e QΩ(T n + q)t + p (T n + q) . (2.27) n∈Z2 Write this function as a Fourier series, ∞ X h(ξ, t) = bk(t)e(kξ). (2.28) k=−∞ We are ultimately interested in the Mellin transform of this function, ∞ −ξ −ξ Z T c1,T c2 ξ> ξ s dt ζ > T ΩT , s = h(ξ, t)t (2.29) b ξ −ξ (T ) p,T q 0 t ∞ X = βk(s)e(kξ), (2.30) k=−∞ where, as we will show, Z ∞ s dt βk(s) := bk(t)t . (2.31) 0 t ω11 ω12 n1 p1 q1 Express Ω = ( ω12 ω22 ), n = ( n2 ), p = ( p2 ), q = ( q2 ). Write ∞ X h(ξ, t) = hn2 (ξ, t) = h0(ξ, t) + eh(ξ, t), (2.32) n2=−∞ where hj(ξ, t) is the sum over the terms with n2 = j, and eh(ξ, t) is the sum over all the terms where n2 6= 0. Also, assume that q1 = q2 = 0. A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 11

First, calculate h0(ξ, t): ∞ X n t1/2 1 h (ξ, t) = ρc1,c2 1 e ω n2t + p n . (2.33) 0 Ω 0 2 11 1 1 1 n1=−∞

The n1 = 0 term of this sum vanishes. We write, for n2 6= 0, Z 1

hn2 (ξ, t)e(−kξ) dξ 0 Z 1 ∞ X c1,c2 n1 + n2ξ 1/2 = ρM t 0 n2 n1=−∞ n1 + n2ξ > n1 + n2ξ · e QΩ t + p e(−kξ) dξ (2.34) n2 n2 n2−1 Z ∞ X c1,c2 n1 + n2ξ 1/2 = ρM t −∞ n2 n1=0 n1 + n2ξ > n1 + n2ξ · e QΩ t + p e(−kξ) dξ (2.35) n2 n2 n2−1 Z ∞ X c1,c2 n2ξ 1/2 = ρM t −∞ n2 n1=0 n2ξ > n2ξ n1 · e QΩ t + p e −k ξ − dξ (2.36) n2 n2 n2 n2−1 ! Z ∞ X kn1 c1,c2 ξ 1/2 = e ρM n2t n2 −∞ 1 n1=0 ξ ξ · e Q n2t + p> n e(−kξ) dξ. (2.37) Ω 1 2 1 2

n2−1 X kn1 The exponential sum e evaluates to |n2| if n2|k, and to 0 otherwise. Thus, for all n2 n1=0 k ∈ Z (including k = 0), Z 1 eh(ξ, t)e(−kξ) dξ (2.38) 0 Z ∞ X c1,c2 ξ 1/2 ξ 2 > ξ = |n2| ρM n2t e QΩ n2t + p n2 e(−kξ) dξ. −∞ 1 1 1 n2|k

Our convention here is that a sum over n2|k ranges over both positive and negative n2, and over all integers when k = 0.

2.4. Shifting the contour vertically. Fix a positive real number λ to be speciﬁed later. Let C+ (C−) be the contour consisting of the horizontal line Im(z) = λ (Im(z) = −λ), + oriented towards the right half-plane. For each d1, d2 ∈ Z, d2 6= 0, let C(d1, d2) be C if 12 GENE S. KOPP

− d1d2 > 0 or d1 = 0 and d2 > 0; let C(d1, d2) be C if d1d2 < 0 or d1 = 0 and d2 < 0. The integrands in eq. (2.38) approach zero as Re(ξ) → ±∞, so we may rewrite this formula using contour integrals

Z 1 eh(ξ, t)e(−kξ) dξ (2.39) 0 Z X c ,c ξ 1 ξ ξ 1 2 2 2 > = |n2| ρM n2t e QΩ n2t + p n2 e(−kξ) dξ. k 1 1 1 C ,n2 n2|k n2

2.5. Taking Mellin transforms term-by-term. To calculate the Mellin transform of h0(ξ, t), we need to check absolute convergence to justify reversing the order of summa- tion/integration.

1 Proposition 2.7. If σ = Re(s) > 2 , then

Z ∞ ∞ 1/2 X c1,c2 n1t 1 2 σ dt ρΩ e ω11n1t + p1n1 t < ∞. (2.40) 0 0 2 t n1=−∞ Proof. We bound the integral as follows.

Z ∞ ∞ 1/2 X c1,c2 n1t 1 2 σ dt ρΩ e ω11n1t + p1n1 t (2.41) 0 0 2 t n1=−∞ Z ∞ ∞ 1/2 σ X c1,c2 t 1 t dt = ρΩ e ω11t 2 (2.42) 0 0 2 n1 t n1=−∞ ∞ ! Z ∞ 1/2 X −2σ c1,c2 t 1 σ dt = |n1| ρΩ e ω11t t (2.43) 0 0 2 t n1=−∞ < ∞. (2.44)

1 The sum converges for σ > 2 , and the integral converges for σ > 0 (because the integrand approaches a constant at t → 0 and decays exponentially as t → ∞).

Therefore, we can switch the sum and the integral. Using Lemma 2.3 and dropping the subscript on n1, Z ∞ s dt X sgn(n)e(p1n) c1,c2 1 h0(ξ, t)t = − 2s κΩ , s (2.45) 0 t |n| 0 n∈Z\{0} 1 = − (Li (e(p )) − Li (e(−p ))) κc1,c2 , s . (2.46) 2s 1 2s 1 Ω 0

Next, we’re going to calculate the Mellin transform of eh(ξ, t). We need an absolute convergence result to justify our calculation here, too. A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 13

1 Proposition 2.8. Suppose σ = Re(s) > 2 . Then,

Z ∞ Z X X c1,c2 ξ 1/2 ρM n2t k 1 0 C ,n2 k∈Z n2|k n2 n26=0 ξ 2 > ξ s dt ·e QΩ n t + p n2 e(−kξ) t dξ < ∞. (2.47) 1 2 1 t

Proof. Let

Z ∞ Z ± c1,c2 ξ 1/2 ξ σ dt K := ρM t e QΩ t t dξ (2.48) 0 C± 1 1 t < ∞. (2.49)

Set K := max{K+,K−}. We have

Z ∞ Z X X c1,c2 ξ 1/2 ρM n2t k 1 0 C ,n2 k∈Z n2|k n2 n26=0 ξ 2 > ξ s dt ·e QΩ n t + p n2 e(−kξ) t dξ 1 2 1 t Z ∞ Z X X c1,c2 ξ 1/2 ξ 2 = ρM n2t e QΩ n2t k 1 1 0 C ,n2 k∈Z n2|k n2 n26=0 dt · e−2πλktσ dξ (2.50) t Z ∞ Z X X c1,c2 ξ 1/2 ξ = ρM t e QΩ t k 1 1 0 C ,n2 k∈Z n2|k n2 n26=0 σ −2πλk t dt · e 2 dξ (2.51) n2 t X X −2πλk −2σ ≤ K e n2 (2.52) k∈Z n2|k n26=0

X X −2πλ|d1d2| −2σ = K e d2 (2.53) d1∈Z d2∈Z\{0} < ∞. (2.54)

The proposition is proved. 14 GENE S. KOPP

Now we may justify taking the Mellin transform of the Fourier series term-by-term. It follows from Proposition 2.8 that

∞ −ξ −ξ Z T c1,T c2 ξ> ξ s dt ζ > T ΩT , s = h(ξ, t)t (2.55) b ξ (T ) p,0 0 t ∞ X = βk(s)e(kξ), (2.56) k=−∞

Z ∞ Z ∞ s dt s dt where βk(s) := bk(t)t . Deﬁne βek(s) := ebk(t)t ; then, 0 t 0 t  c1,c2 1  − (Li2s(e(p1)) − Li2s(e(−p1))) κΩ , s + βe0(s) if k = 0, βk(s) = 0 (2.57)  βek(s) if k 6= 0.

Proposition 2.8 also implies that we can switch the order of integration to compute

Z ∞ Z 1 s dt βek(s) = eh(ξ, t)e(−kξ) dξ t (2.58) 0 0 t Z X > ξ −2s c1,c2 = |n2| e n2p − kξ − sgn(n2) |n2| κΩ (ξ, s) dξ (2.59) k 1 C ,n2 n2|k n2 Z X sgn(n2) c1,c2 = − 2s−1 e(n2(p1ξ + p2) − kξ) κΩ (ξ, s) dξ. (2.60) k |n2| C ,n2 n2|k n2

2.6. Series manipulations. In this subsection, we set ξ = 0 in eq. (2.56). We will manipu- late the right-hand side of this equation to prove Theorem 1.8. First of all, we have

∞ c1,c2 X ζbp,0 (Ω, s) = βk(s) (2.61) k=−∞ ∞ 1 X = − (Li (e(p )) − Li (e(−p ))) κc1,c2 , s + β (s). (2.62) 2s 1 2s 1 Ω 0 ek k=−∞

k We will rewrite the sum of the βk(s) using the substitution (d1, d2) = ( , n2). The following e n2 manipulation is legal by Proposition 2.8.

∞ Z X X X sgn(n2) c1,c2 βek(s) = − 2s−1 e(n2(p1ξ + p2) − kξ) κΩ (ξ, s) dξ (2.63) k |n2| C ,n2 k=−∞ k∈Z n2|k n2 n26=0 Z X X sgn(d2) c1,c2 = − 2s−1 e(d2(p1ξ + p2) − d1d2ξ) κΩ (ξ, s) dξ. (2.64) |d2| C(d1,d2) d1∈Z d2∈Z\{0} A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 15

Split up the series into four pieces.

∞ Z X X X e(d2p2) c1,c2 βek(s) = − 2s−1 e(−(d1 − p1)d2ξ) κΩ (ξ, s) dξ |d2| C− k=−∞ d1>0 d2>0 Z X X e(d2p2) c1,c2 + 2s−1 e(−(d1 − p1)d2ξ) κΩ (ξ, s) dξ |d2| C+ d1>0 d2<0 Z X X e(d2p2) c1,c2 − 2s−1 e(−(d1 − p1)d2ξ) κΩ (ξ, s) dξ |d2| C+ d1≤0 d2>0 Z X X e(d2p2) c1,c2 + 2s−1 e(−(d1 − p1)d2ξ) κΩ (ξ, s) dξ (2.65) |d2| C− d1≤0 d2<0 Z X X e(d2p2) c1,c2 = − 2s−1 e((d1 − p1)d2ξ) κΩ (−ξ, s) dξ |d2| C+ d1>0 d2>0 Z X X e(d2p2) c1,c2 + 2s−1 e(−(d1 − p1)d2ξ) κΩ (ξ, s) dξ |d2| C+ d1>0 d2<0 Z X X e(d2p2) c1,c2 − 2s−1 e(−(d1 − p1)d2ξ) κΩ (ξ, s) dξ |d2| C+ d1≤0 d2>0 Z X X e(d2p2) c1,c2 + 2s−1 e((d1 − p1)d2ξ) κΩ (−ξ, s) dξ (2.66) |d2| C+ d1≤0 d2<0 Z X X e(d2p2) c1,c2 = − 2s−1 e((d1 − p1)d2ξ) κΩ (−ξ, s) dξ d2 C+ d1>0 d2>0 Z X X e(−d2p2) c1,c2 + 2s−1 e((d1 − p1)d2ξ) κΩ (ξ, s) dξ d2 C+ d1>0 d2>0 Z X X e(d2p2) c1,c2 − 2s−1 e((d1 + p1)d2ξ) κΩ (ξ, s) dξ d2 C+ d1≥0 d2>0 Z X X e(−d2p2) c1,c2 + 2s−1 e((d1 + p1)d2ξ) κΩ (−ξ, s) dξ. (2.67) d2 C+ d1≥0 d2>0

Now, move the contour integral outside the sums, and rewrite the series as

∞ Z d2 d2 X X e(−p2 + p1ξ) c1,c2 X e(p2 + p1ξ) c1,c2 βek(s) = 2s−1 κΩ (−ξ, s) − 2s−1 κΩ (ξ, s) (2.68) C+ d2 d2 k=−∞ d2≥0 d2≥0

X X 1 d2 d2 c1,c2 + 2s−1 −e((d1 − p1)ξ + p2) + e((d1 + p1)ξ − p2) κΩ (−ξ, s) d2 d1>0 d2>0

d2 d2 c1,c2 + e((d1 − p1)ξ − p2) − e((d1 + p1)ξ + p2) κΩ (ξ, s) dξ. 16 GENE S. KOPP

Setting s = 1, we obtain ∞ X βek(1) (2.69) k=−∞ Z c1,c2 −ξ c1,c2 ξ = − log(1 − e(−p2 + p1ξ))κΩ + log(1 − e(p2 + p1ξ))κΩ C+ 1 1 ∞ X −ξ + (log (1 − e((d − p )ξ + p )) − log (1 − e ((d + p )ξ − p ))) κc1,c2 1 1 2 1 1 2 Ω 1 d1=1 ξ (− log (1 − e((d − p )ξ − p )) + log (1 − e((d + p )ξ + p ))) κc1,c2 dξ. 1 1 2 1 1 2 Ω 1 We want to write this sum of logarithms as a logarithm of a product, but there is the issue of the choice of branch. In order to make a clear choice, let ∞ Y 1 − e ((d + p1)ξ + p2) ϕp1,p2 (ξ) := (1 − e (p1ξ + p2)) (2.70) 1 − e ((d − p1)ξ − p2) d=1 for ξ ∈ H. This is a function on the upper half-plane which is never zero, and the upper half-plane is simply connected, so it has a choice of continuous logarithm. Let (Log ϕp1,p2 )(ξ) be the branch such that log(1 − e(p2)) if p1 = 0, lim (Log ϕp1,p2 )(ξ) = (2.71) ξ→i∞ 0 if p1 6= 0.

Here log(1 − e(p2)) is the standard principal branch. Thus, ∞ X Z −ξ β (1) = − (Log ϕ )(ξ) · κc1,c2 ek p1,−p2 Ω 1 k=−∞ C+ ξ + (Log ϕ )(ξ) · κc1,c2 dξ. (2.72) p1,p2 Ω 1

∞ c1,c2 X Adding back the other piece of β0(1) into ζbp,0 (Ω, 1) = βk(1), we obtain k=−∞ 1 ζc1,c2 (Ω, 1) = − (Li (e(p )) − Li (e(−p ))) κc1,c2 (2.73) bp,0 2 1 2 1 Ω 0 Z c1,c2 −ξ + − (Log ϕp1,−p2 )(ξ) · κΩ (2.74) C+ 1 ξ + (Log ϕ )(ξ) · κc1,c2 dξ. (2.75) p1,p2 Ω 1

2.7. Collapsing the contour onto the branch cuts. We could declare ourselves done at c1,c2 this point. Equation (2.73) is a formula for ζbp,0 (Ω, 1), as we desired, and it appears very diﬃcult to evaluate or simplify the contour integral in any way. However, eq. (2.73) is not a useful formula for computation because the integral converges slowly. The integrand decays polynomially as ξ → ±∞ along the horizontal contour C+. A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 17

We will obtain a Kronecker limit formula with rapid convergence by shifting the contour so that the integrand decays exponentially. In doing so, we will also split up the formula as a diﬀerence of a c1-piece and a c2-piece. The movement of the contour is shown in Section 2.7. c i > Let Λ := Ω − Mcc M for c = c1, c2, as we did in Corollary 2.5. Factor the quadratic Ω QM (c) ξ polynomial Q c in ξ, ΛΩ 1

ξ QΛc = α(c)(ξ − τ1(c))(ξ − τ2(c)). (2.76) Ω 1

c (0) Since ΛΩ ∈ H2 by Lemma 3.6 of [5], we know by Lemma 2.2 that we may choose τ1(c) to be in the upper half-plane and τ2(c) in the lower half-plane. c ξ The complex function ξ 7→ κΩ 1 has branch cuts along the vertical ray from τ1(c) to i∞ and the vertical ray from τ2(c) to −i∞. We check that this function is holomorphic away c ξ ξ from these branch cuts. Since κΩ 1 has simple poles at the roots ξ = r1, r2 of QΩ 1 = 0, c1,c2 ξ we must check that the residues at the poles cancel when taking the diﬀerence κΩ 1 = c2 ξ c1 ξ κΩ 1 − κΩ 1 . We have

c ξ res κΩ ξ→r1 1 ξ c>M 1 = lim (ξ − r1) s (2.77) ξ→r1 ξ ξ 2 ξ 2πiQ c>M − 2iQ (c)Q Ω 1 1 M Ω 1 ξ c>M 1 = lim s (2.78) ξ→r1 ξ 2 ξ πiω (ξ − r ) c>M − 2iQ (c)Q 11 2 1 M Ω 1 1 = , (2.79) πiω11(r1 − r2) and similarly, res κc ξ = 1 . These residues do not depend on c, so they cancel, Ω 1 πiω11(r2−r1) ξ→r2 c1,c2 ξ and κΩ 1 is holomorphic at r1 and r2. ±ξ Move the countours of integration above the zeros of QΩ 1 . Now we may safely split up the integral into a term for c1 and a term for c2. Now we retract the integral onto the branch cut. As ξ = ±τ ± + ε and ε → 0, the denominator of the integrand blows up like ε1/2, so the integral converges. The integrand changes sign when we cross the branch cut. Thus, eq. (2.73) becomes

c1,c2 + − + − ζbp,0 (Ω, 1) = I (c2) − I (c2) − I (c1) + I (c1), (2.80) 18 GENE S. KOPP

τ1(c2) τ1(c2) r1 r1

τ1(c1) τ1(c1)

C+

r2 r2

τ1(c2) r1 r1 τ1(c2)

τ1(c1) τ1(c1)

r2 r2

+ c ξ Figure 1. The contour C is moved above the poles of κΩ 1 , then collapsed onto branch cuts. where

1 I±(c) := − Li (e(±p ))κc 2 1 Ω 0 Z ∞ ± ± c ± (τ (c) + it) + 2i (Log ϕp1,±p2 )(±τ (c) + it)κΩ dt. (2.81) 0 1

We have now proven Theorem 1.8. Theorem 1.9 follows by specializing the variables, setting g Ω = iM and restricting to c1, c2 ∈ R . Theorem 1.10 and Theorem 1.11 both follow by application of the functional equation (Theorem 1.6). A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 19

3. Example We conclude with an example to show how to use the Kronecker limit formula for indeﬁnite zeta functions√ to compute Stark√ units. This example was introduced in Section 7.1 of [5]. Let K = ( 3), so O = [ 3], and let c = 5O . The ray class group Cl (O ) ∼ /8 . Q K Z √ K c∞2 K = Z Z 0 The fundamental√ unit ε = 2 + 3 is totally positive: εε = 1. It has order 3 modulo 5: ε3 = 26 + 15 3 ≡ 1 (mod 5). In this section, we use the Kronecker limit formula for 0 indeﬁnite zeta functions to compute ZI (0), where I is the principal ray class of Clc∞2 (OK ). 2 0 1/5 2 Let M = ( 0 −6 ), q = 0 , and c1 ∈ R any column vector with the property that > 0 c1 Mc1 < 0, such as c1 = ( 1 ). By Corollary 1.16 and the discussion in Section 7.1 of [5], we have

3 0 c1,P c1 ZI (0) = ζb0,q (Ω, 0), (3.1)

2 3 where Ω = iM and P = ( 1 2 ). Now we want to use Theorem 1.11 to compute the right-hand side of eq. (3.1). If we try 3 3 26 45 3 45 3 −15 ΩP c1 to do so directly, we obtain P = ( 15 26 ), P c1 = ( 26 ), ΩP c1 = 6i ( 26 ), and Λ−Ω−1 = 675 390 ξ −i . The root of Q 3 in the upper half-plane (equivalently, the branch point 390 676/3 ΛΩP c1 1 −Ω−1 3 √ ΩP c1 ξ −2340+i 3 of κ−Ω−1 1 in the upper half-plane) is ξ = 4053 , which is very close to the real axis. log(10)N That means we’d need to use about √ ≈ 857.5N terms in the product expansion of 2π 3/4053 0 ϕ−q1,q2 (ξ) to compute ZI (0) to N decimal places of accuracy. We technically have exponential decay, but it’s not very useful. It is much more practical to break up the zeta function into pieces. We can also improve −1 the rate of convergence by choosing c1 optimally; here, we will use c = ( 1 ) in place of c1. We have

0 c,P 3c ZI (0) = ζb0,q (Ω, 0) (3.2) c,P c P c,P 2c P 2c,P 3c = ζb0,q (Ω, 0) + ζb0,q (Ω, 0) + ζb0,q (Ω, 0) (3.3) c,P c c,P c c,P c = ζb0,q (Ω, 0) + ζb0,q0 (Ω, 0) + ζb0,q00 (Ω, 0), (3.4)

1 1 0 1 2 00 1 2 where q = 5 ( 0 ), q = 5 ( 1 ), and q = 5 ( 4 ) are obtained from the residues of the global units 0 1 2 ε , ε , ε modulo 5. √ √ Ωc ξ −3 6(ξ−1) ΩP c ξ 3 6(ξ+1) Now, we have κ −1 = √ and κ −1 = √ , with −Ω 1 π(3ξ2−1) 3ξ2−3ξ+1 −Ω 1 π(3ξ2−1) 3ξ2+3ξ+1 √ √ 3+i 3 −3+i 3 branch points in the upper half-plane at ξ = 6 and ξ = 6 , respectively. We log(10)N thus need to use about √ ≈ 1.269N terms in the product expansion of each of the 2π 3/6 0 functions ϕ (ξ), ϕ 0 0 (ξ), and ϕ 00 00 (ξ) to compute Z (0) to N decimal places of −q1,−q2 −q1,−q2 −q1 ,−q2 I accuracy by this method. For q, q0, q00, we computed the corresponding values of the integrals J(c),J 0(c),J 00(c) and J(P c),J 0(P c),J 00(P c) given by eq. (1.19). The computation was per- formed in Mathematica using numerical integral of the first 40 terms of the product expansion 20 GENE S. KOPP of each ϕ. For the differences of the two integrals, we obtain J(P c) − J(c) ≈ − 0.05923843917544488329354507987 + 3.65687839020311786132893850239i, (3.5) J 0(P c) − J 0(c) ≈ − 1.33733021085943469210685014899 + 0.52477812529424663387556899167i, and (3.6) J 00(P c) − J 00(c) ≈ 2.64057587271922212456484190607 + 0.52477812529424663387556899167i. (3.7) For the ray class zeta value, we thus calculate using Theorem 1.11 that 0 c,P c c,P c c,P c ZI (0) = ζb0,q (Ω, 0) + ζb0,q0 (Ω, 0) + ζb0,q00 (Ω, 0) (3.8) 2i 2i 2i = √ Im(J(P c) − J(c)) + √ (J 0(P c) − J 0(c)) + √ (J 00(P c) − J 00(c)) det M det M det M (3.9) 1 = √ Im (J(P c) − J(c) + J 0(P c) − J 0(c) + J 00(P c) − J 00(c)) (3.10) 2 3 ≈ 1.35863065339220816259511308230. (3.11) This agrees (to 30 decimal digits) with the computations described in Section 7.1 of [5]. 0 The conjectural Stark unit is exp(ZI (0)) ≈ 3.89086171394307925533764395962. This number is appears to be the root of the polynomial √ √ √ √ x8 − (8 + 5 3)x7 + (53 + 30 3)x6 − (156 + 90 3)x5 + (225 + 130 3)x4 √ √ √ − (156 + 90 3)x3 + (53 + 30 3)x2 − (8 + 5 3)x + 1, (3.12) which we have verified lies in the appropriate class field.

4. Acknowledgements This research was partially supported by National Science Foundation (USA) grants DMS- 1401224, DMS-1701576, and DMS-1045119, and by the Heilbronn Institute for Mathematical Research (UK). This paper incorporates material from the author’s PhD thesis [4]. Thank you to Jeﬀrey C. Lagarias for advising my PhD and for many helpful conversations about the content of this paper. Thank you to Marcus Appleby, Jeﬀrey C. Lagarias, and Kartik Prasanna for helpful comments and corrections.

5. Conflict of interest statement On behalf of all authors, the corresponding author states that there is no conﬂict of interest.

References [1] William Duke, Ozlem¨ Imamo˘glu, and Arpad´ Tóth.Kronecker’s first limit formula, revisited. Research in the Mathematical Sciences, 5(2):1–21, 2018. [2] Gustav Herglotz. Uber¨ die Kroneckersche grenzformel fürreelle, quadratische Körper. I, II. Berichte über die Verhandl. d. SächsischenAkad. der Wiss. zu, 75:3–14,31–37, 1923. [3] Gene S. Kopp. A generalised Kronecker second limit formula. In preparation. A KRONECKER LIMIT FORMULA FOR INDEFINITE ZETA FUNCTIONS 21

[4] Gene S. Kopp. Indefinite Theta Functions and Zeta Functions. PhD thesis, University of Michigan, Ann Arbor, Michigan, USA, August 2017. [5] Gene S. Kopp. Indefinite zeta functions. Research in the Mathematical Sciences, 8(17):1–34, 2021. [6] Sheng-Chi Liu and Riad Masri. A Kronecker limit formula for totally real fields and arithmetic applica- tions. Research in Number Theory, 1(8):1–20, 2015. [7] Anna Posingies. Kronecker limit formulas and scattering constants for Fermat curves. arXiv preprint arXiv:1111.1132, 2011. [8] Danylo Radchenko and Don Zagier. Arithmetic properties of the Herglotz function. arXiv preprint arXiv:2012.15805, 2020. [9] Christina Roehrig. Siegel theta series for indefinite quadratic forms. Research in Number Theory, 7(3):1– 28, 2021. [10] Takuro Shintani. On Kronecker limit formula for real quadratic fields. Proceedings of the Japan Academy, 52(7):355–358, 1976. [11] Takuro Shintani. On a Kronecker limit formula for real quadratic fields. J. Fac. Sci. Univ. Tokyo, 24(1):167–199, 1977. [12] Takuro Shintani. On certain ray class invariants of real quadratic fields. Journal of the Mathematical Society of Japan, 30(1):139–167, 1978. [13] Carl L. Siegel. Advanced Analytic Number Theory. Tata Intstitute of Fundamental Research, 1980. [14] Harold M. Stark. Values of L-functions at s = 1. I. L-functions for quadratic forms. Advances in Mathe- matics, 7(3):301–343, 1971. [15] Harold M. Stark. L-functions at s = 1. II. Artin L-functions with rational characters. Advances in Mathematics, 17(1):60–92, 1975. [16] Harold M. Stark. L-functions at s = 1. III. Totally real fields and Hilbert’s twelfth problem. Advances in Mathematics, 22(1):64–84, 1976. [17] Harold M. Stark. L-functions at s = 1. IV. First derivatives at s = 0. Advances in Mathematics, 35(3):197– 235, 1980. [18] Marie-France Vignéras.Séries thêtades formes quadratiques indéfinies. Séminaire Delange-Pisot-Poitou. Théoriedes nombres, 17(1):1–3, 1976. [19] Marie-France Vignéras.Sériesthêtades formes quadratiques indéfinies.In Modular Functions of One Variable VI, pages 227–239. Springer, 1977. [20] Maria Vlasenko and Don Zagier. Higher Kronecker “limit” formulas for real quadratic fields. J. Reine Angew. Math., 679:23–64, 2013. [21] Don Zagier. A Kronecker limit formula for real quadratic fields. Math. Ann., 213:153–184, 1975. [22] Sander Zwegers. Mock Theta Functions. PhD thesis, Universiteit Utrecht, Utrecht, The Netherlands, October 2002.

School of Mathematics, University of Bristol, Bristol, UK and Heilbronn Institute for Mathematical Research, Bristol, UK Email address: [email protected]