L-functions for Meromorphic Modular Forms and Sum Rules in Conformal Field Theory

We define L-functions for meromorphic modular forms that are regular at cusps, and use them to: (i) find new relationships between Hurwitz class numbers and traces of singular moduli, (ii) establish predictions from the physics of T-reflection, and (iii) express central charges in two-dimensional conformal field theories (2d CFT) as a literal sum over the states in the CFTs spectrum. When a modular form has an order-$p$ pole away from cusps, its $q$-series coefficients grow as $n^{p-1} e^{2 \pi n t}$ for $t \geq \sqrt{3}/2$. Its L-function must be regularized. We define such L-functions by a deformed Mellin transform. We study the L-functions of logarithmic derivatives of modular forms.L-functions of logarithmic derivatives of Borcherds products reveal a new relationship between Hurwitz class numbers and traces of singular moduli. If we can write 2d CFT path integrals as infinite products, our L-functions confirm T-reflection predictions and relate central charges to regularized sums over the states in a CFTs spectrum. Equating central charges, which are a proxy for the number of degrees of freedom in a theory, directly to a sum over states in these CFTs is new and relies on our regularization of such sums that generally exhibit exponential (Hagedorn) divergences.

In this paper, we define L-functions attached to meromorphic modular forms that are regular at cusps. To begin, consider an integer-weight modular form f (τ ) defined on the upper half-plane H, and bounded away from cusps. Because f (τ ) = f (τ + 1), it has a Fourier decomposition, f (τ ) = n c f (n)e 2πinτ . Taking the Mellin transform of f yields the associated L-function [1], |c f (n)|. In section III, we evaluate L reg f (s) at s = 0 and s = 1, when f is a generating function for traces of singular moduli, ∂ τ log H d (j(τ )) [2], and preform a consistency check on L reg f (s). The second part of the paper applies these mathematical results to quantum field theory (QFT). In section IV, we use the special values of the L-functions found in sections II and III to show that central charges of many two-dimensional conformal field theories (2d CFTs) be written directly as a tally of states in the theory. The terms involved in these sums increase exponentially with n. Hence, this statement crucially hinges on our L-function regularization. Further, we use these L-functions to explicitly verify two sum-rules motivated by T-reflection [6,7]. The stress-tensor dictates both the central-charge sum-rule and the T-reflection sum-rule. This is striking as the T-reflection phase, and its associated sum-rule, seems to be a new gravitational anomaly [6,7]. Connecting it to the stress-tensor, which also tracks the conformal anomaly [8], seems natural. We conclude in section V.
A. Motivations from field theory to study meromorphic modular forms Three physics considerations are the driving motivation for this work. First, there is a persistent connection between the behavior of quantum chromodynamics (QCD) at low energies and the behavior of strings. Roughly, the confining potential for QCD can be modeled by a string of finite tension. As this is a model, it does not capture all of the features of QCD, particularly at high energies. Yet, there are some commonalities.
Heuristically, vacuum loops of closed QCD-strings have the topology of a two-torus. Thus, we may expect the one-loop path integral for low-energy, low-temperature QCD to be modular invariant. (See Refs. [9][10][11] for an explicit example.) However, low-energy QCD is known to have an exponential Hagedorn growth in the spectrum: d(E) ∼ E p−1 e +β H E [12,13]. This leads to poles in the path integral Z(β) when the inverse temperature, β nears β H : Juxtaposing the expected modularity of the low-energy QCD path integral with the presence of Hagedorn poles suggests that the theory of meromorphic modular forms may give an interesting set of tools to study low-energy QCD. If d(E) exhibits Hagedorn growth, then current "zeta-function" technology cannot regularize the Casimir energy sum, 1 2 E d(E) E. An L-function for a meromorphic modular form, however, might.
Second, on a related note, L-functions for meromorphic modular forms may allow us to rewrite central charges for unitary CFTs in terms of a direct sum over states in the spectrum, even if the spectrum exhibits Hagedorn growth. While central charges are a measure of the degrees of freedom in the theory, which decrease monotonically along RG flows that progressively integrate-out modes on short distance scales, it is pleasing to see a central charge written directly as a tally of states. Our L-functions allow us to do this.
Third, we would like to test a prediction of the conjectured invariance of QFT path integrals under reflecting their temperatures to negative values (T-reflection). Based on physical arguments, Refs. [6,14], conjectured that if a path integral Z(τ ) is a weight-k modular form that can be written as an infinite product of the form, where q := e 2πiτ and the d(n) are integers for all n, then we should have However, save when Z(τ ) is a quotient of Jacobi theta functions or Dedekind eta functions, these d(n)-exponents exhibit Hagedorn growth, and the regularization procedure "reg" is not known. Our L-functions for meromorphic modular forms gives us a precise context in which we may regulate these sums. Using them, we confirm these predictions from [6,14].

B. Summary of main mathematical and physical results
We now summarize our main results. First, in section II, we prove the connection between poles and Hagedorn growth of q-series coefficients, and then provide a well-defined regularized L-function that can be attached to every meromorphic modular form. The following fact establishes Hagedorn growth for meromorphic modular forms: Proposition I. 1. Let f be a meromorphic modular form with a pole of order p at some finite location (away from cusps). There exist real numbers t ≥ √ 3/2 and C > 0 such that f (τ ) = n c f (n)q n has Cn p−1 e +2πnt ≤ |c f (n)| ≤ e 2πn(t+ǫ) for infinitely many positive integers n, for any ǫ > 0.
This qualitative feature distinguishes meromorphic modular forms from their weakly holomorphic cousins and has played prominently in many papers, such as the last joint paper between Hardy and Ramanujan [15], and in the recent work of Bialek, Berndt and Ye [16][17][18] and Bringmann, Kane, Lobrich, Ono, and Rolen [19][20][21]. However, this result may have escaped broader notice. It is a very important connection between meromorphic modular forms and field theory limits of string path integrals, which often feature Hagedorn growth and Hagedorn poles. We have highlighted it for this reason.
The following statement gives the (polar part) of a well-defined regularization of the integral L * f (s) reg when f is a meromorphic modular form: Theorem I.2. Let f be a weight k meromorphic modular form that is regular at cusps. Then it's L-integral, L * f (s) reg , evaluates to where 0 < B < ∞ and "regular(s)" denotes terms that are regular for finite s ∈ C.
Theorem I.2 is a fusion of Theorem II.14 and Corrolary II. 15, and is one of the main result of this paper. In sections III and IV, we use L * f (s) reg to establish the following two results: Theorem I. 3. Suppose that F (τ ) is a meromorphic modular form of weight k. The Lfunction for T (τ ) := q∂ q log F (τ ) is,

7)
where 0 < B < ∞ and again "regular(s)" denotes terms that are regular for finite s ∈ C.
Corollary I. 4. Suppose that F (τ ) is a meromorphic modular form of weight k, and can be written as q −∆ n (1 − q n ) −d(n) . Define f (τ ) := n d(n)q n . Then, (1. 8) In section IV, we argue that if a 2d CFT path integral equals F (τ ), then it can be interpreted as the partition function for the CFT in the grand canonical ensemble, where particle number is not fixed. Further, we argue that when the CFT is free, then f (τ ) can be interpreted as the partition function for a single-particle excitation in the CFT, i.e. the partition function for the CFT in the canonical ensemble. When a 2d CFT has a path integral with modular weight k and an infinite product expansion, then T-reflection suggests the two sum-rules in Eq. (1.8) [6,7]. We use the L-functions in Eq. (1.6) to verify this "prediction" from T-reflection. This is one of the main physics-results in this paper.
The other main physics-results that come from our L-functions derive from the interpretation of the quantities ∆ and k in Eqs. (1.7) and (1.8). Namely, in the context of a unitary 2d CFT, ∆ is the Casimir energy and is directly proportional to the central charge of the theory, c: ∆ = c/24. In writing ∆ as the special value of the L-function of Z CAN (τ ), we directly equate the Casimir energy and thus the central charge to a direct tally of the number of states in the CFT. Because the number of states in the canonical ensemble has an exponential Hagedorn growth, the only way that we make this statement is by using the L-functions for meromorphic modular forms. In this precise sense, this is a new statement.
The main mathematical application of our L-functions in this paper is in the context of traces of singular moduli. All of the notation in this paper on traces of singular moduli follows Zagier's paper of the same name [2]. To state the application, let H(d) denote the Hurwitz class number of a quadratic of negative discriminant d, t n (d) be the trace of the unique modular function J m (τ ) = q −m + O(q) ∈ M ! 0 on the CM points discriminant-d called α Q , which are defined to be the solutions to discriminant −d quadratics that lie within F , w Q be the order of the stabilizer of α Q , and let f d be the unique element in the Kohnen plus- In Theorems 3 and 5 of [2], Zagier showed that the d th Hilbert polynomial H d (j(τ )) (which is a weakly holomorphic modular function that) satisfies where A(n 2 , d) is the coefficient of q n 2 in the Fourier decomposition of, This relationship between Hurwitz class numbers and traces of singular moduli, while formal, is new and crucially hinges on the definition of L-functions for meromorphic modular forms. Finally, we point-out that the exponential growth in Proposition I.1 gives an amusing way to effectively determine H(d) without direct reference to quadratics. Explicitly: Theorem I.5. Consider the two expressions for Λ d (τ ) := q∂ q H d (j(τ )) ∈ F 2 from the two different representations of H d (j(τ )) in Eq. (1.9), and let Li 0 (x) := n x n . The q-series coefficients are t n (d). The q-series coefficientst n (d) of the pole-subtracted function, 13) have exponential growth that is bounded by |t n (d)| < e πn √ 2 .
We can find H(d) by determining the number of terms terms that need to be subtractedout from Λ d (τ ) before its q-series coefficients are bounded by |t n (d)| < e πn √ 2 .

II. L-FUNCTIONS FOR MEROMORPHIC MODULAR FORMS
Write Re(z) and Im(z) for the real and imaginary parts of a complex number z ∈ C. Let H := {τ ∈ C | Im(τ ) > 0} denote the upper half-plane. A holomorphic modular form of weight k (for the full modular group, SL 2 (Z)) is a holomorphic function f : H → C such that f (τ + 1) = f (− 1 τ )τ −k = f (τ ) for all τ ∈ H, and f (τ ) = O(1) as Im(τ ) → ∞. Write M k for the vector space of modular forms of weight k. A product of modular forms of weights k 1 and k 2 is a modular form of weight k 1 + k 2 , so M * := k M k is naturally a (graded) ring.
In this work a meromorphic modular form of weight k is a quotient f = f 1 f 2 where f i ∈ M k i , the denominator f 2 is not identically zero, and k = k 1 − k 2 . We write F k for the vector space of meromorphic modular forms of weight k, and set F * := k F k . Then F k vanishes unless k is an even integer, and F * may be regarded as a "graded ring of fractions" of M * .
Define F (∞) to be the vector space of meromorphic functions on H that satisfy f (τ ) = f (τ + 1), have only finitely many poles in any compact subset of the vertical strip and satisfy the growth condition f (τ ) = O(e CIm(τ ) ) as Im(τ ) → ∞ for some C > 0. Then F * is a subspace of F (∞). Set e(x) := e 2πix and q := e(τ ). For any f ∈ F (∞) we write c f (n) for the coefficient of q n in the Fourier expansion if f is a holomorphic cusp form) then c f (n) = O(n k ) as n → ∞, and the Dirichlet series converges absolutely for Re(s) > k + 1 [1]. But as we show in Lemma II.6, general f ∈ F k have c f (n) that grow exponentially and the right-hand side of (2.3) is nowhere convergent. In this section we introduce a regularization of (2.3), L reg f (s), that is well-defined and analytic in s for an arbitrary meromorphic modular form f ∈ F k . To introduce it, recall that when f ∈ M k , the L-integral is given by: Here, when f ∈ M k , the integral and the q-series commute and yield the sum representation of Γ(s)L f (s)/(2π) s from Eq. (2.3). The analogous statement for f ∈ F k requires regularization. Generic f ∈ F k have poles within the fundamental domain F := {τ ∈ C | 0 < Im(τ ), |τ | ≥ 1, |Re(τ )| ≤ 1/2} and poles at i∞. For this case, we define the L-function, L f (s), by a deformation of the associated L-integral. Let L * f (s) reg be the following limit of the regularized contour integral: where γ(Λ, ±ǫ) begins at τ = i/Λ, moves upwards just to the left (right) of the imaginary-τ axis, and ends at τ = iΛ. Deformations away from the Im(τ )-axis allow us to deal with poles along the imaginary-τ axis. Crucially, if f ∈ M k then Eqs. (2.3), (2.4) and (2.5) match. The remainder of this section is as follows. In section II A we describe this contour and make contact with a previous definition of L reg f (s) for the special case f ∈ M ! k . In section II B we prove that if f ∈ F k has a single order-p pole at τ = s+it ∈ F , then |c f (n)| ∼ Cn p−1 e 2πnt for some C. In section II C we rewrite the integral transform in Eq. (2.5) explicitly in terms of the q-series coefficients of f ∈ F k . In section II D, we comment on the numerical convergence of L reg f (s) and on path dependence.

A. Defining the contour and defining the regularization
We define the regularized L-integral, L * f (s) reg , by the contour γ(Λ, ǫ), which begins at i/Λ and goes to iΛ, while keeping just ǫ to the right of poles on the imaginary-τ axis: Lemma II.1. Let f ∈ M k and B be a positive real number. Then the L-integral, and as Γ(s, x) ∼ x s−1 e −x for large-x, the sum in Eq. (2.7) converges absolutely for any finite s ∈ C for every f ∈ M k .
Proof. As f ∈ M k does not have any poles, then any integration contour from i/Λ to iΛ yields the same result. So the exact shape of our particular contour, γ(Λ, 0), from i/Λ to iΛ does not have any effect on the final result. As f ∈ M k is bounded along this finite contour, the contour integral converges and we may write As Λ → ∞, the integrals in the sum on the RHS evaluate to Γ(s)/(2πn) s . Thus the sum over c f (n)× integrals yields Eq. (2.6), which matches the definition of L f (s) in Eq. (2.3). Now, we can split-up the integral from t = 1/Λ to t = Λ as the sum of an integral from t = 1/Λ to t = B and an integral from t = B to t = Λ. Noting f (−1/τ ) = τ k f (τ ), we can rewrite the integral from t = 1/Λ to t = B in the following way: Note that this process is independent of the value of B. Freedom to choose 1/Λ < B < Λ plays an important role in section II C, where we derive explicit forms for L * f (s) reg for f ∈ F k . Now, for any fixed Λ ≫ 1 we may exchange the integration and the q-series summation, in Eq. (2.9). Recalling the gamma-function Γ(s, When Re(s) > 0 and Re(s − k) > 0, we find that the Λ → ∞ limit of Eq. (2.9) equals We finish the proof by analytically continuing 1/s and 1/(k − s) to the entire s-plane, and then noting that the scaling of Γ(s, x) ∼ x s−1 e −x implies the sum converges absolutely for any finite value of s ∈ C when f ∈ M k and c f (n) = O(n 2k−1 ).
The sum n c f (n)Γ(s, 2πn)/(2πn) s converges absolutely for every finite s ∈ C for any sequence of c f (n) whose exponential growth is |c f (n)| max{e 2πn/B , e 2πnB }. Now, we note: k , and Λ ≫ 1, and consider the regularized contour integral ). For any fixed Λ, it evaluates to: As there exists C > 0 for any f ∈ M ! k where |c f (n)| ≪ e C √ n for n ≫ 1, this sum in Eq. (2.12) converges absolutely for finite s ∈ C. The explicit result in Eq. (2.12) matches L * f (s) reg in Refs. [3][4][5]. Yet, the number of n < 0-terms in Eq. (2.13) diverge as Λ → ∞. Proof. In this proof, we reproduce Eq. (2.12) and make contact with the previously defined L-functions for weakly holomorphic modular forms in [3][4][5].
First, we show that our proof of Eq. (2.7) in Lemma II.1 applies without essential modification for any finite Λ ≫ 1: As f ∈ M ! k only has poles at cusps, all integration contours from τ = i/Λ to τ = iΛ yield the same result. So again we can path-deform either of the two finite-Λ integrals that appear in Eq. (2.12) to the much simpler integral iΛ i/Λ dt t s−1 (f (it) − c f (0)). We now exploit the modularity to again rewrite the nontrivial Λ-regulated integral as in Eq. (2.10). When Re(s) > 0 and Re(s − k) > 0, we again find that the integral is: Now, note that the Λ-independent terms in Eq. (2.14) exactly match those in Eq. (2.6), and exactly match those in [4,5]. However, if f has a q-series that begins with q −N for some positive integer N, then there exists a finite number of non-zero coefficients c f (n) for n < 0.
In the sum over n = 0, these terms multiply incomplete gamma functions Γ(σ, 2πnΛ), which diverge exponentially as Λ → ∞. These terms ruin the convergence of L * f (s|Λ) reg as Λ → ∞. If not for these terms, then the terms in Eq. (2.15) would dominate L * f (s|Λ) reg in the Λ → ∞ limit, and we would recover the results of [4,5]. This finishes the proof.
Crucially, this does not reproduce the L-functions for f ∈ M ! k in [3][4][5]. If we could equate the regularized L-function with the constant term of L * f (s|Λ) reg in a large-Λ expansion, then we would reproduce their results. However, we are unaware of any reason to justify such an equality. Hence, we restrict our attention to f ∈ F k that are regular at cusps. Within this space, we prove the following Lemma: Lemma II.3. Let f ∈ F k be regular at cusps. Then the Λ → ∞ and ǫ → 0 limits of exist and yield a finite result. Hence L reg f (s) as defined in Eq. (2.5) exists for generic f ∈ F k . Proof. We require the outermost contours γ(Λ, ±ǫ) to pass to left and to the right of each on-axis pole exactly once, while enclosing no poles off the axis. Now for every f ∈ F k there is a finite integer N such that and exponentially decays as τ → i∞ (and as τ → i0). In the limit ǫ → 0, each contour encloses poles on the axis but remains finite as Λ → ∞. Because the two contour integrals differ by a finite sum of residues, when f is regular at i∞, then lim ǫ→0 ( ±ǫ 1 2 lim Λ→∞ L * f (s|Λ, ǫ)) exists and yields a finite result. Thus L * f (s) reg exists and is finite when f ∈ F k is regular at cusps.
We now comment that if the regularization prescription in Lemma II.2 applied to modular forms with poles at cusps, f ∈ M ! k , then proof of Lemma II.3 would generalize to all f ∈ F k . The reason is as follows.
Let f ∈ F k have a q-series that begins with q N (1 + O(q)) with N < 0. Now, consider any g ∈ M ! k whose polar terms match those of f ∈ F k . More explicitly, f (τ ) + g(τ ) = O(q 0 ). It is relatively simple to construct an explicit function whose polar terms c g (n)q −n match those of any f ∈ F k . Let h ∈ M k be holomorphic, and let J M ∈ M ! 0 be the unique modular function that begins with q −M + O(q). Then one can choose |N| coefficients C(k) for k ∈ {1, . . . , |N|} such that g(τ ) := h(τ ) Now, f + g ∈ F k has poles only within the interior of F , while g ∈ M ! k has poles only at the cusp. Exploiting the linearity of the integral functional, we have L * f (s) reg = L * f +g (s|Λ) reg − L * g (s|Λ) reg . If the regularization prescription in Lemma II.2 worked to regularize the Lintegral for f ∈ M ! k , then we would have L * f (s|Λ, ǫ) reg equal to a sum of integrals that are regular and convergent as Λ → ∞ and ǫ → 0. This would extend Lemma II.3 to all f ∈ F k .
Yet, when f ∈ F k has poles at cusps, there are a finite number of terms in L * f (s|Λ, ǫ) reg from the leading polar powers of q in the q-series that diverge at the upper end of integration as Γ(s, −2πnΛ) as Λ → ∞. Thus, as it stands, the regularization only applies to f ∈ F k that are regular at cusps. We briefly comment on several aspects of our regularization procedure: 1. Another contour γ between i0 + and i∞ would define another L * f (s) reg γ . If f has a finite number of poles "between" γ and the Im(τ )-axis, then L * f (s) reg γ − L * f (s) reg differ by a finite sum of residues. Residues of poles to left and right of the Im(τ )-axis have opposite signs. See sections II D and V for more discussion of path dependence.
2. The contour-regularization of the pole at the cusp in our proof of Lemma II.2 is distinguished from the regularization in [3][4][5]. Here, they deformed (f (τ ) − c f (0)) → (f (τ )−c f (0))e 2πυIm(τ ) . When f ∈ M ! k , the υ-regularized L-functions can be analytically continued to υ = 0. We regularize our L-functions by contour-deformation so that we could more easily discuss the case where f ∈ F k may have poles within F .
3. Unfortunately, the regularization prescriptions here and in [4] are in tension with each other: once the integrand is non-holomorphic then the Cauchy residue theorem no longer applies. Non-holomorphic deformations are central to [4], where they define L * f (s) reg for f ∈ M ! k with poles at cusps. Residue theorems are central for our approach, where we define L * f (s) reg for f ∈ F k with poles away from cusps.
can be extended to nontrivial integrands in many ways. For instance, we could define L * f (s) reg γ by the nonholomorphic integral γ dt t s−1 (f (τ ) − c f (0)), where t = Im(τ ), or we could regularize the poles by a non-holomorphic deformation of the integrand.

5.
Defining L * f (s) reg as a meromorphic integrand integrated along contours that pass to both sides of the the Im(τ )-axis gives an unambiguous definition of the regulated integrals that occur when poles lie along the Im(τ )-axis. By using contour-deformations and residue theorems, we obtain uniquely defined, concrete, and well-behaved expressions for L * f (s) reg = reg ∞ 0 dt t s−1 (f (it) − c f (0)) for f ∈ F k that are regular at cusps. In section II C we give explicit formulae for L * f (s) reg analogous to the sum over incomplete gamma functions and q-series coefficients in Lemmas II.1 and II.2, when f ∈ F k is regular at cusps. Before this, in section II B we discuss the lockstep relationship between exponential growth of the q-series coefficients of f ∈ F k and the locations of its poles within F .

B. Poles and exponential growth
We now show the relationship between the location of the poles of f ∈ F k and the exponential growth of its q-series coefficients. To start, let R τp : F (∞) → F (∞) denote a linear operator that projects onto the subspace of functions that are regular at τ p . To explicitly define it, write f (τ ) = n∈Z r f,τp (n)(τ − τ p ) n for the Laurent series expansion of f in a (sufficiently small) deleted neighborhood of τ p , and set for m > 0. We define R τ P , and the related P τp : F (∞) → F (∞), as where Li s (z) := n>0 z n n −s is the polylogarithm function with order s. Note that the summation in (2.19) is finite since r f,τp (n) = 0 for n ≪ 0. If f is regular at τ p , so that r f,τp (n) = 0 for n < 0, then R τp f = f . When m is a positive integer, as z → 0 we have Lemma II.4. For τ p ∈ H and f ∈ F (∞) the function R τp f is regular at τ p .
Proof. If f is regular at τ p then R τp f = f and the claim is true. Otherwise, combining (2.18), (2.19) and (2.20) we see that Observe that if X is any compact subset of H then, since any f ∈ F (∞) has only finitely many poles τ p ∈ X, we obtain a well-defined projection operator R X : 22) and the functions in the image of R X are regular in X. In what follows we consider the set V Y defined by the real number 0 < Y < ∞: When Y is a finite and positive real number, X = V Y is a rectangular subset of V , and we define R Y f and P Y f to respectively be R V Y f and P V Y f . Now consider f ∈ F (∞), and note that the imaginary parts of the poles of f are bounded. The next lemma explains how the maximum of these imaginary parts bounds the growth of the Fourier coefficients of f : Proof. By Parseval's identity we have that converges for any t > t 0 . The claim follows.
Next we verify a converse to Lemma II.5, showing that the coefficients of f ∈ F (∞) grow exponentially when f has poles in H.
Lemma II.6. Let f ∈ F (∞), let t 0 be the maximal imaginary part of a pole of f , and let m 0 be the maximal order of a pole of f with maximal imaginary part. Then there exists a constant C > 0 such that |c f (n)| > Cn m 0 −1 e 2πnt 0 for infinitely many positive integers n.
Proof. Let t 0 and m 0 be as in the statement of the lemma, and let τ 0 be a pole of f with Im(τ 0 ) = t 0 and r f,τ 0 (−m 0 ) = 0 and r f,τ 0 (n) = 0 for n < −m 0 . Set g(τ ) : is regular at τ = τ 0 , so the Fourier series g(τ ) = c g (n)e 2πinτ converges absolutely at τ = τ 0 , so |c g (n)|e −2πnt 0 → 0 as n → ∞. In particular, for any ǫ > 0 we have ǫ > |c g (n)|e −2πnt 0 for n sufficiently large. Take C = 1 2 C m 0 and ǫ = 1 2 C. If |c f (n)| ≤ Cn m 0 −1 e 2πnt 0 for sufficiently large n then we have for n sufficiently large. This is a contradiction, so we must have |c f (n)| > Cn m 0 −1 e 2πnt 0 for infinitely many positive integers n.
Lemmas II.4, II.5, and II.6 establish a direct correspondence between f ∈ F (∞) having poles of order-m at τ p and terms in its q-series that exponentially as n m−1 (e −2πiτp ) n . We will find it useful to associate every function f ∈ F (∞) with two positive real numbers, which we denote U and D and define as follows: The bound of U, D ≥ √ 3/2 ultimately comes from the fact that if τ ∈ F , then Im(τ ) ≥ √ 3 2 . This bound on D, together with Lemmas II.4, II.5, and II.6, implies the following result: Proposition II.7. Let f ∈ F k have a pole of order p within H away from the cusp. Then there exist real numbers t ≥ √ 3 2 and C, p > 0 such that the c f (n) in f (τ ) = n c f (n)q n obey Cn p−1 e +2πnt ≤ |c f (n)| ≤ e 2πn(t+ǫ) for infinitely many positive integers n, for every ǫ > 0.
Proof. If f ∈ F k has an order-p pole at τ = s + it = i∞ for t > 0 and s, t ∈ R, then Lemmas II.4, II.5, and II.6 imply that a real number C > 0 exists, such that the q-series coefficients c f (n) are bounded by Cn p−1 e +2πnt ≤ |c f (n)| ≤ e 2πn(t+ǫ) for every ǫ > 0. As f is modular, it has at least one pole within F . Hence there is an order-p pole with t ≥ where ǫ is any positive non-zero number, m f is maximum pole-order of the pole below yet nearest to the line Im(τ f ) < B, and β(B, f ) is the imaginary part of the pole(s) in f with maximal imaginary part that lie below the strip V B . It is defined by (2.29) Note that 0 < β < B. Thus, there exist ǫ > 0 such that β < ǫ + β < B.

C. The regularized L-function
We now write explicit expressions for L * f (s) reg when f ∈ F k . Recall D in Eq. (2.26) and let B ≤ D. Define γ ± B (Λ, ǫ) to be the part of γ(Λ, ǫ) above or below the line Re(τ ) = B, and define Sγ + Proof. We prove the Lemma by identical steps to those in the proof of Lemma II.2.
Lemma II.9. If Re(s) > 0 and Re(s − k) > 0, then independent of ǫ it follows that Proof. This identity follows exactly as in the proofs of Lemma II.1 and Lemma II.2.
Let U, D and β be as defined in Eqs. (2.26), (2.27), and (2.29) and choose B < D such that it is not equal to the imaginary part of any pole of f . This is important in Lemma II.13, where we explicitly integrate τ s−1 Li −N (e(τ − τ p )) along the contours γ + B (Λ, ǫ) and Sγ − B (Λ, ǫ). After evaluating the constant term in Lemma II.9, the remaining part of L * f (s|Λ, ǫ) reg is a nontrivial integral over contours contained within V B . We partition this nontrivial integral into a part that is regular within the strip V B (but contains poles at ∞) and a part which contains all of the poles of f within V B . Our partitioning involves the following integrals: . Thus, the integrands and the contours of these convergent integrals match.
Lemma II.11. Let f ∈ F k be regular at cusps, and let Λ > U. Then, where c Rf (n|B) are the q-series coefficients of R B f (τ ), whose growth is bounded by |c Rf (n|B)|e 2π(ǫ−B)n → 0 for some ǫ > 0. This sum converges absolutely.
Proof. To show Eq. (2.34), note that R B f (τ ) explicitly lacks poles in V B , that L * Rf (s|Λ, ǫ|B) depends only on the end-points of the contour, and that all of the relevant contours are contained within V B . So we may replace the contours with contours along the vertical axis.
To show Eq. (2.35), we proceed as in the proof of Lemma II.2 and commute the q-series sum with the integral at finite Λ. This is possible because the integrals converge at finite Λ. Recall that by Lemmas II.4, II.5 and II.6 we know that there exists an ǫ > 0 such that Thus, the two sums exponentially converge for any s ∈ C. Finally, noting Lemma II.12. Let f ∈ F k , and let Λ > U. Then Proof. The proof follows directly from the definition of the function P B f (τ ) ∈ F (∞).
Explicit results for L * f (s) reg come from evaluating the integral of (τ /i) σ−1 Li −N (e(τ − τ p )) along the contour γ + (S)B (Λ, ǫ) as Λ → ∞ and then as ǫ → 0, for σ = s and σ = k − s. This is unambiguous when τ p is off the imaginary-τ axis. However, when f ∈ F k has poles at imaginary values, τ p = iy, we must be more careful. Here we begin by defining This is unambiguous and finite for any value of Λ ≫ y and for any value of ǫ = 0. By holomorphy, I N (y, s|Λ, ǫ|B) depends only on the sign of ±ǫ. When t ≫ y, |Li −N (e 2π(y−t+iǫ )| = O(e −2πt ), and we may safely take Λ → ∞. The resulting finite integral, I N (s, y|ǫ|B) := lim Λ→∞ I N (s, y|Λ, ǫ|B), depends only on the sign of ǫ. More precisely, where Res(N, s, iy) is the residue of −i(τ /i) s−1 Li −N (e(iy − τ )). Explicitly, it is given by: Note that each integral ±ǫ I N (s, y| ± ǫ|B) is separately finite. We now define, We may have well-defined expressions for the integrals in Eq. (2.37) for any τ P ∈ V B .
Lemma II.13. Let Im(τ p ) = y = B. Then I N (s, −iτ p |B) as defined above equals (2.42) Proof. The essence of the proof comes from evaluating I N (s, −iτ p |D) when the pole τ p = iy is on the Im(τ )-axis. There are only two cases to consider: when y < B and when y > B.
When y < B, there is no pole along the Im(τ )-axis and we may write, which converges because e 2πny Γ(s, 2πnB) ∼ e 2πn(y−B) /(2πn) s−1 exponentially decays. When y > B, then we may break-up the contour into three pieces γ + , in the limit where Λ → ∞. As emphasized above, we define I N (s, y|B) as the average of the contour integrals along γ + B (Λ, ±ǫ) is finite. Now, each contour circles the pole at τ = iy by an angle of ±π. Averaging over contours γ + B (Λ, ±ǫ) thus counts the residue at τ = iy exactly + 1 2 − 1 2 = 0 times. Because of this, we can safely evaluate I N (s, y|B) by considering the integral evaluated at the end-points t = B and at t = ∞. Further, because Li −N (e −2πt ) = O(e −2πt ) as t → ∞, we know that the integral evaluated at the upper bound of the integration contour vanishes.
So we are left with the problem of evaluating the integral at the lower end-point t = B. To compute this, we note that when N is a positive integer, we have This allows us to rewrite the integral evaluated at the lower end-point t = B as, Thus, because y > B, we have e −2πny |Γ(s, −2πnB)| ∼ e −2πn(y−B) /(2πn) s−1 for n ≫ 1. Thus this sum converges exponentially quickly for generic finite s ∈ C. We finish the proof by noting that the above manipulations apply equally well to τ p = iy + x when x = 0.
Theorem II.14. Let f ∈ F k be regular at cusps, and let B, c Rf (n|B), Γ(s, x), τ p , V B , r * f,τ P (m), and I N (s, z|A) as defined above. Then the L-integral that yields the L-function is Proof. In Lemma II.3, we defined L * f (s) reg as the Λ → ∞ and ǫ → 0 limit of the average of the integrals over contours that begin at τ = i/Λ and end at τ = iΛ and pass just ǫ the left, and just ǫ the right, of poles on the imaginary-τ axis. Each integral converges, and yields a finite result at fixed Λ > U and ǫ = 0, and each has a finite limit as Λ → ∞ and ǫ → 0.
Then in Lemma II.8, we split-up each of individual contour into a sum of three terms: two contour integrals with integrands (τ /i) σ−1 (f (τ ) − c f (0)), and two contour integrals with . Following this, in Lemma II.10, we split-up the integral over the nontrivial contours, which are entirely contained in the strip V B ⊂ V into two pieces. One piece is regular within V B but contains poles at cusps whose two integrands are proportional to R B f (τ ) − c f (0). The other piece is regular at the cusp but has poles within V B whose two integrands are proportional to P B f (τ ). We call these terms L * Rf (s|Λ, ǫ|B) reg and L * P f (s|Λ, ǫ|B) reg , respectively. By appealing to the q-series growth bounds in Lemmas II.4, II.5, and II.6 and Proposition II.7 and to Lemma II.2, we showed that the which converge exponentially quickly for finite s ∈ C. In Lemmas II.11, II.12, and II.13 we found explicit expressions for I N (−iτ p , s|B), the integrals within P B f (τ ). Combining these results yields Eq. (2.50) and completes the proof of the Theorem.
This has an important Corollary: . Thus, when evaluating L * f (s|Λ, ǫ) reg we inflect around B, and when evaluating L * f (k − s|Λ, ǫ) reg we inflect around 1/B. The results for L * f (s) reg and L * f (k − s) reg are each given wholly by terms in the first line of Eq. (2.50), and are identical up to the overall factor i k . (See similar analysis in [4] for L * f (s) reg when f ∈ M ! k .) This procedure holds for f ∈ F k : we define L * f (s|Λ, ǫ) reg by inflecting around Im(τ ) = B, and writing (the nontrivial part of) L * f (s|Λ, ǫ) reg in terms of contour integrals along . We then split-up the integral into the sum of two integrals, whose integrands are (τ /i) s−1 ( R B f (τ ) − c f (0)) and (τ /i) s−1 P B f (τ ). These integrands are defined by the property that all poles "below" τ = i∞ and above the line at Im(τ ) = B are projected out of R B f and are entirely contained in P B f . We call the corresponding integrated expressions, respectively, L * Rf (s|Λ, ǫ|B) reg and L * P f (s|Λ, ǫ|B) reg . The key point here is that regardless of reflecting about the line Im(τ ) = 1/B or Im(τ ) = B, to obtain convergent expressions for the integrals L * Rf (s|Λ, ǫ|B) reg and L * P f (s|Λ, ǫ|B) reg we must project-out all poles above the lowest of the two lines Im(τ ) = B or Im(τ ) = 1/B.
We now obtain L * f (k − s|Λ, ǫ) reg in a form that directly compares with L * f (s) reg in the following way. We preform the same procedure on L * f (k − s|Λ, ǫ) reg , save that we now reflect the contour γ(Λ, ǫ) about the line Re(τ ) = 1/B. Now, we project-out all of the poles of f that lie below the lowest of the two lines Im(τ ) = B or Im(τ ) = 1/B, to obtain the Proof. Special values of L * f (s) reg within the critical strip, s ∈ [0, k], for f ∈ F k are given by Eq. (2.50). When f ∈ F 0 , then by Theorem II.14 the residues of the polar terms B s However when f ∈ F k for k = 0, the two polar terms do not cancel, thus: We now consider k < 0. Here, integers in the critical strip are N ∈ {k, k + 1, . . . , −1, 0}. Importantly, Γ(s) has a pole at each of these critical values of s. Therefore, We now consider k > 0. Here, the integers in the critical strip are N ∈ {0, 1, . . . , k − 1, k}. Importantly, Γ(s) is regular at each of these critical values of s. So, In this section, we briefly comment on several aspects of our definition of L-functions for meromorphic modular forms, before moving on to applications. First, we discuss writing meromorphic modular forms in terms of a Rademacher sum over polar terms and the related issue of numerical evaluations of our L-functions. Second, we discuss ambiguities in our Lfunctions which arise when poles cross the integration contour, and whether this affects the special values of L reg f (s) in Theorem II.14 and Corollary II. 16. The pole-subtraction procedure used to define R B f and P B f is extremely simple when f ∈ F k has negative-definite weight, k < 0. When k < 0, then meromorphic modular forms can be written as a convergent sum over non-positive weight polylogarithms that have poles at all distinct modular images of all of the poles of f . For example, as Hardy and Ramanujan showed [15] that 1/E 6 (τ ) can be written in the following way: where Γ := SL 2 (Z), and Res(i) is the residue of 1/E 6 (τ ) on the pole at τ = i. Crucially, the q-series coefficients for 1/E 6 (τ ) increase exponentially, as where the first sum is over all coprime integers (c, d) = 1, and the second sum is over the related quantity λ = c 2 + d 2 . For a fixed coprime pair (c, d) there is a pair of elements of Γ/Γ ∞ given by γ ± = ±a b c ±d , which yield D λ(c,d) (n): e (Im(γ ± (i))) . (2.63) More generally, when f ∈ F k has a single simple pole within the fundamental domain, τ p ∈ F , then it can be written similarly as a kind of Rademacher sum over Li 0 (e(τ − γτ p )): This has been widely extended to meromorphic modular forms for k < 0 in [16][17][18][19][20] that have multiple poles within F of arbitrary order. In this context, the pole-subtraction procedure used to define R B f and P B f is particularly simple to implement. When the general results of [19,20] apply, we obtain rapidly convergent expressions for c Rf (n|B) by simply deleting terms from poles above the Im(τ ) = B line. Our ability to numerically evaluate L * f (s) reg , as defined in Theorem II.14, appear to depend crucially on having Rademacher-like sums for meromorphic modular forms f ∈ F k .
However, when k = 0, then sums akin to Eq. (2.64) do not converge absolutely and do not reproduce meromorphic modular functions without modification. When k > 0 then these sum no longer converge absolutely; when k ≥ 1 they do not converge at all. In the absence of convergent Rademacher-sums over polar q-series terms, such as Li 1−p (e(τ − γτ p )), for k ≥ 0, then it is less obvious how to write convergent expressions for the c Rf (n|D) in the pole-subtracted sums that we take to define the regularized L-function for generic meromorphic modular forms. We know of one example where this has been explicitly done for k > 0: When k = 2, Bringmann et al [21] build on the work of Bruinier et al [22] and explicitly construct sums of this sort for weight-two quasimodular forms, given by ∂ τ log f for some modular form f . We expect numerically evaluate our expressions for L * f (s) when f ∈ F k to be crucially tied to the existence of convergent sums such as in [21] when f ∈ F k for k ≥ 0.
We now briefly discuss the path-(in)dependence of our result for L Here R B f is regular in the vertical strip V B , and P B f contains every pole of f in V B . We then integrated f = R B f + P B f along the Im(τ )-axis, from τ = i0 + to τ = i∞. Now, precisely because f has poles within F , this definition is path-dependent.
Had we defined the L-function by a different path γ, denoted L * f (s) reg γ , we would have where m i is the pole-order, Res(m i , s, τ i ) is the residue at τ i in Eq, (2.40), and χ(τ i ) = ±1 tracks if the pole τ i , which lies between γ and the Im(τ )-axis, is to the right or left of the axis. In this sense, the L-functions depend explicitly on path.
Similarly, consider a continuous family of meromorphic modular forms indexed by the location of one of their poles, z. When we smoothly send z to γz for γ ∈ SL 2 (Z), the contour integral that defines L * f (s) reg will acquire residues corresponding to how many times images of the pole at z cross the vertical Im(τ )-axis. These additional pieces resemble the discontinuities in theta-decompositions of meromorphic Jacobi forms in [23,24], and could thus also represent interesting wall-crossing phenomena.
Importantly, these additional contributions do not have any poles at finite values of s ∈ C. Hence, they do not change the residues at the poles 1/s and 1/(k − s) in L * f (s) reg . Therefore, the special values discussed in Corollary II.16 are insensitive to this ambiguity. Because the main results in sections III and IV that we derive from our L-functions for f ∈ F k come from these special values, they are not sensitive to this ambiguity.
Finally, we reiterate that L * f (s) reg applies to f ∈ F k that are regular at cusps. This is complimentary to, but does not directly extend or generalize, the results in [4,5] where they define L * f (s) reg for f ∈ M ! k , which have poles exclusively at cusps. However, as commented in the proof of Lemma II.2, up to a finite number of divergent terms in the infinite sum in Eq. (2.13), the results in Eqs. (2.12) and (2.13) exactly their results. It will be important to fuse these two results, and define regularized L-integrals and L-functions for all f ∈ F k .

III. SPECIAL VALUES OF REGULARIZED L-FUNCTIONS
This section has two principal components. The first part is in sections III A, III B, and III C. The second part is in section III D. Throughout this section, Γ := SL 2 (Z).
In sections III A-III C we present one of the main results of this paper, which concerns L-functions for weight-two meromorphic modular forms Λ d ∈ F 2 , where H(d) is the Hurwitz class-number, and that counts the number of distinct quadratics with negative discriminant −d = 4ac − b 2 > 0, and t n (d) is the trace of the unique modular function J n (τ ) = 1/q n +O(q) ∈ M ! 0 over CM-points of a quadratic with negative discriminant −d, and d ≡ 0, 3 (mod 4). Evaluating the L-function of Λ d (τ ) at s = 0 gives a sum-rule that relates the traces over singular moduli and the Hurwitz class numbers, where "reg" denotes regularization via L-function. This relationship between the Hurwitz class numbers and traces of singular moduli is new: it rests upon the L-functions for meromorphic modular forms defined in section II. Later on, in section IV, we relate the special values in sections III A and III B to novel statements about Casimir energies, central charges, and a new reflection-symmetry, for special two-dimensional conformal field theories (CFTs). Finally, in section III C, we show an amusing way to calculate Hurwitz class-numbers via bounding the exponential growth of traces of singular moduli as a function of d.
Then, in section III D, we study a subset of meromorphic modular forms that vanish identically as τ → i0 + . We relate these results to the special values in Corollary II.16 of Theorem II.14, at negative integer values of s and at s = 0. Each result has the interpretation as the sum of exponentially diverging sequences of numbers. Where they overlap, they agree perfectly. We view this as a check on the consistency L reg f (s) when f ∈ F k .

A. Traces of singular moduli and the Hurwitz class numbers
One of our main applications of our L-functions concerns a new relationship between the number of Γ-inequivalent quadratics with a fixed negative discriminant, −d, and the regularized sum of traces of modular functions, when traced over solutions to Γ-inequivalent quadratics with negative discriminant −d.
To begin, we define −d to be the negative discriminant of a quadratic operator, Q(X, Y ) = aX 2 + bXY + cY 2 , where a, b, and c are all integers. Here, −d = 4ac − b 2 . Further, we define the CM point α Q to be a solution of a given quadratic with negative discriminant −d such that In Eq. (3.4), we called the space of quadratics with negative discriminant −d by the name Q d . We now associate the number w Q to each distinct quadratic with discriminant −d, Summing over all w Q for the finite sum Q ∈ Q d /Γ gives Hurwitz class numbers in Eq. (3.1): Weighting each term in the sum by the value of J m (τ ) at the unique root/CM-point α Q ∈ F of the quadratic Q ∈ Q d /Γ, then we obtain the traces of singular moduli t n (d): Finally, we define the d th Hilbert class polynomial as, It is straightforward to show that, In [2], Zagier noted that this follows directly from the definitions and noted that is the generating function of the traces of J n (τ ), from Eq.
Further, in [2], Zagier proved that the Hilbert class polynomials H d (j(τ )) span a subspace of meromorphic modular forms that have Borcherds product expansions [25]. The connection to Borcherds products is a crucial aspect of [2] and of this paper. We will revisit it first in section III B, and then throughout section IV. However, our current focus is a precise statement of the special values of the L-function of the generating function for the traces of singular moduli, L reg Λ d (s). Using Theorem II.14 and Corollary II.16 from section II, we find the L-function of Λ d (τ ): Here "regular(s)" refers to the terms in L * f (s) reg that are finite for finite s ∈ C (see Theorem II.14 for details). Using this L-function, we see immediately Now, if we interpret the special value of L reg Λ d (s) at s = 0 from Eq. (3.13) as a sum of traces of singular moduli of the J m (τ ) at imaginary quadratic points α Q with discriminant −d, then it is a formal relationship between traces of singular moduli and Hurwitz class numbers. This relationship between traces of singular moduli and the H(d) is new. It is important to note that the location of the simple poles in Λ d (τ ) in Eq. (3.11) at the points α Q in Eq. (3.4), combined with Lemmas II.5 and II.6, together imply the refined information that where γ max · α Q is the first modular image of α Q that is distinct from α Q , and ǫ ≥ 0. As Im(γ max · α Q ) < 1/ √ 2, these growth conditions must be less than e πn √ 2 . (See Theorem III.1 in section III C for details.) Without L-functions for meromorphic modular forms, there would be no meaning to this sum of exponentially growing q-series coefficients of Λ d ∈ F 2 . We explore several of consequences of this in section III B and again in section IV.
We now emphasize that the vanishing in Eq. (3.15) happens term-by-term in the sum over Q ∈ Q d /Γ. If we define, then linearity of the integral L * f +g (s) reg = L * f (s) reg + L * g (s) reg , implies Hence, each L reg Λα Q (s) in the sum over Q ∈ Q d /Γ in Eq. (3.18) individually vanishes as s → 0. Before moving on to Borcherds products, we give a mnemonic for understanding the exponential growth of the q-series coefficients of Λ d (τ ). By Eq. (3.11), we have c Λ d (n) = t n (d). Now, t n (d) is "mostly" J n (i √ d/2 + x), where x = 0 or 1/2 depending on d. Further, J n (τ ) is "mostly" e −2πinτ . Thus the dominant contribution to |t n (d)| is "mostly" e +nπ √ d .

B. Traces of singular moduli and Borcherds products
Traces of singular moduli, and the related Hurwitz class numbers, are famously related to Borcherds infinite product formulas for modular forms. Using this, we reinterpret the special values of L reg Λ d (s) in section III A. This reinterpretation will be of considerable interest in physical applications, when the quantum field theoretic path integral has an infinite product expansion (with integer exponents). In [2] Zagier proved that the modular polynomials H d (j(τ )) in Eq. (3.6) are an important class of Borcherds products. Recall that d = 4ac − b 2 ≡ 0, 3 (mod 4). Thus, where A(n 2 , d) can either be interpreted as c(n 2 ) in the q-series of a particular modular form f d (τ ) that is an element of the Kohnen plus-space of M ! 1/2 (Γ 0 (4)), or as c(d) in the q-series of the related particular modular form g n 2 (τ ) ∈ M ! 3/2 (Γ 0 (4)). The modular forms f d (τ ) and g n 2 (τ ) are defined by their modular weight in the subgroup Γ 0 (4) and their q-series growth: (3.20) The g n 2 (τ ) in Eq. A(n, 0)q n = 1 + n∈Z q n 2 = θ(τ ) ∈ M 1/2 (Γ 0 (4)) , where H(0) = −1/12 and H (τ ) is the holomorphic part of Zagier's weight-3/2 mock modular form,Ĝ(τ ). In this context, we see that the q-series coefficients of the Λ d (τ )s can be recast in terms of the q-series coefficients of g D (τ ) or equivalently of f d (τ ): Removing the factor of ζ(s) implies the two formal sum-rules, This is a new relationship between the q-series coefficients of the various g D (τ ). This relationship will be of special use and importance in section IV where we study CFTs whose path integrals have Borcherds products. As we show in section IV B, Casimir energies for these CFTs are directly analogous to the Hurwitz class numbers, H(d), from Eqs. (3.13) and (3.27). In section IV C we comment that both sum-rules in Eq. (3.29) were anticipated by a recently observed symmetry of path integrals in quantum field theory [6,7,14].
Before concluding, we emphasize that de-convolving the formal Dirichlet series in Eq. (3.25) into the product of formal Dirichlet series in Eq. (3.26) is well-motivated on physical grounds. We discuss the physical basis for this factorization in section IV B.
C. Using Borcherds exponents to count quadratics with discriminant −d Poles above the Im(τ ) = √ 3/2-line correspond to q-series growth in excess of e πn √ 3 , as proven in Lemmas II.5 and II.6. Bounding the growth conditions for the Borcherds exponents, A(n 2 , d), is thus directly sensitive to the unique roots of quadratics with negative discriminant −d. Put differently, bounding the A(n 2 , d) directly counts the Hurwitz class numbers, H(d) = A(0, d). We emphasize this fact, even though it lies somewhat outside the main scope of the paper. It is most efficient to state this result as a Theorem: Theorem III.1. Consider the two expressions for Λ d ∈ F 2 defined in Eqs. (3.10) and (3.11). The q-series coefficientst n (d) of the pole-subtracted function,
Proof. To begin the proof, we note that the pole-subtraction procedure removes all poles from Λ d (τ ) that lie either within, or on the boundary of, the fundamental domain. However, Λ d (τ ) has an infinite number of poles within the strip |Re(τ )| ≤ 1/2 that come from the modular orbits of the finite set of points α Q . It is straightforward to show that the maximum imaginary part of any image of any α Q under Γ, which is neither within the fundamental domain nor its boundary, is Im( To prove the Theorem, we overlay the two expressions for Λ d (τ ) in Eqs. (3.10) and (3.11): Now, Λ d (τ ) has a unique simple pole with residue 1/w Q for every point α Q . Because Λ d ∈ F 2 , it has poles with unit residues at all modular images of every distinct α Q . Every α Q is within the fundamental domain, and hence has Im(α Q ) > 1/ √ 2. By Lemmas II.5 and II.6, then the pole-subtracted function, is regular at all points within the fundamental domain, and at all points on its boundary. So any pole inΛ d (τ ) has Im(τ ) < 1/ √ 2. Thus the only source of exponential growth in the q-series coefficients ofΛ d can come from poles below this line. Thus, |J d (n)| < e πn √ 2 .
Theorem III.1 gives an alternative way to determine the Hurwitz class-numbers, directly from the A(n, d): by cancelling the leading growth of the q-series coefficients of Λ d (τ ). Subtract-off exponentially growing contributions to the A(n 2 , d) for any given fixed-d, until growth is bounded by Exp[πn √ 2]. The number of terms which must be subtracted before hitting this bound gives the number of quadratics with negative discriminant −d. L-functions associated with modular forms represent a regularization of the formally divergent sum, n c(n)/n s . In this section, we derive related results for the behavior of meromorphic modular forms, f = n c(n)q n , in the limit where q goes to one. To this end, we prove Lemmas III.2 and III.3, which agree with the results from Corollary II.16. Specifically, Lemma III.2 applies to the case when f ∈ F k vanishes cusps and but has arbitrary weight, while Lemma III.3 applies to the case when f ∈ F k is bounded at cusps and has negative modular weight.
We begin by defining the space of quasimodular meromorphic forms of weight-(k + 2∆) given by the ∆-fold τ -derivative acting on elements of F k . We denote this space byF (k,2∆) : Further, we define the ∆ th descendant of a specific f ∈ F k by f (∆) ∈F (k,2∆) , where g m (τ ) ∈ F k+2(∆−m) and E 2 (τ ) is the weight-two quasimodular holomorphic Eisenstein series. We now state the first Lemma: Lemma III.2. Suppose f ∈ F k , and f (τ ) = n≥1 c(n)q n . Then for any k and for every non-negative integer ∆ ≥ 0, it follows that n ∆ c(n)q n = 0 . Because f ∈ F k , we can relate behavior at τ → i∞ to behavior at τ → i0 + , and have (3.37) A crucial point here is that the decay in (3.36) is exponential rather than power-law. Thus, f (τ ) decays exponentially both when τ → i∞ and when τ → i0 + , regardless of its weight k. Second, we consider ∆ > 0. When ∆ = 1, it is straightforward to see f (1) (τ ) = g 0 + E 2 g 1 , where g 0 ∈ F k+2 and g 1 = k 12 f ∈ F k . Because f vanishes at the cusp, we see that the product E 2 f also vanishes at cusps. Therefore, so too must g 0 . Iterating this implies that when f ∈ F k vanishes at the cusp, then every τ -derivative of f has powers of E 2 (τ ) multiplied by modular forms that vanish as e 2πniτ as τ → i∞. Thus, if f = O(q) and f ∈ F k , then This holds for ∆ ≥ 1. As τ → i0 + , f (∆) vanishes as e −2πi/τ /τ k+2∆ . This follows from the fact that f (∆) decays as e −2π|τ | when τ goes to ∞. As g m ∈ F k+2m vanishes at q = 0, then which vanishes exponentially quickly for any integer ∆ ≥ 0. This completes the proof.
Crucially, each term in the polynomial expansion of f (∆) in terms of E 2 (τ ) m g m (τ ) has negative weight when ∆ < |k|/2. Inspecting examples, however, shows that while the constant term of f (∆) (τ ) vanishes for ∆ > 0, the individual g m ∈ F k+2(∆−m) do not vanish at cusps. And so the only way to ensure that the individual terms, E ∆−m 2 g m , do not diverge as τ → i0 + is to constrain their quasimodular weight to be negative. Now, consider the boundary case ∆ = |k|/2. Here, f (∆) = f (|k|/2) is a meromorphic quasimodular form of weight-zero whose q-series that is O(q). Hence, the leading polynomial divergences from the g m E m 2 terms cancel when τ → i0 + : Because the leading term is O(τ 0 ) and the sub-leading terms decay as τ → i∞ as e −2π|τ | /τ or faster, we may safely focus on the first term. Now, note that this leading term is simply equal to f (|k|/2) itself. Further, recall that f (|k|/2) does not have a constant term. Thus, This fails when |k|/2 ≤ ∆ ≤ |k|. When |k|/2 ≤ ∆ ≤ |k|, then the corresponding f (∆) is a positive weight quasi-modular form. Because the g m (τ ) have constant terms, and f (∆) in this range has positive weight, we conclude that f (∆) diverges as τ → i0 + . The situation changes when ∆ ≥ |k| + 1. Using Bol's identity [26], we conclude that f (|k|+1) is an exactly modular meromorphic form within F |k|+2 = F k+2|k|+2 that lacks a constant term. Thus, by Lemma III.2, we can conclude that whenever ∆ ≥ |k| + 1. This completes the proof.
We now recap. In Corollary II. 16 we derived concrete expressions that stand in for the L-functions of meromorphic modular forms, at s = 0 and at negative integer values of s. These L-functions formally represent the finite part of the nowhere convergent object, n c(n)/n s , at these special values. Independently, in Lemmas III.2 and III.3, we exploited modularity to explicitly compute the finite limit of meromorphic modular forms, and their τ -derivatives, when τ → i0 + . This gives an independent computation of sums of the form, n n ∆ c(n)q n in the limit where q goes to one. Within their regions of overlap, these two completely independent methods to regularize the sum of q-series coefficients are entirely consistent. This overlap occurs in two situations.
First, in Lemma III.2 we study the behavior of meromorphic modular forms f ∈ F k for any k that vanish at cusps, and their meromorphic quasi-modular descendants f (∆) ∈F (k,2∆) . We show that f (∆) (τ ) vanishes exponentially quickly as τ → i0 + for every integer ∆ ≥ 0. When written in terms of the q-series of f (∆) , this amounts to the statement that, lim q→1 ∞ n=1 n ∆ c(n)q n = 0 . when c(0) vanishes. Second, in Lemma III.3 we study meromorphic modular forms f ∈ F k for k ≤ 0, and their quasi-modular derivatives f (∆) ∈F (k,2∆) , that are bounded at cusps. Specifically, we show that when τ is sent to i0 + , f (∆) ∈F (k,2∆) vanish when τ → i0 + for 0 ≤ ∆ ≤ |k|/2. These vanishing results match the special values of the L-function of f ∈ F k for k ≤ 0, spelled-out in Corollary II.16, at integer values of s within the critical strip k ≤ s ≤ 0. Further, by Bol's identity [26], we show that the limit τ → i0 + of f (∆) vanishes when ∆ ≥ |k| + 1. This matches the trivial zeros of L reg f (s) at negative integers s, enforced by the 1/Γ(s) factor built into the L-function of f ∈ F k . Explicitly, Note that Lemma III.3 evaluates the L-function for negative weight meromorphic modular forms at negative integers in the critical strip, which lies between s = 0 and s = −|k| (recall k = −|k| < 0). As the evaluations in this section rest on a direct appeal to modularity, they are limited to τ -derivatives of the original f ∈ F −|k| with non-positive weight. We can appeal to the symmetry L * f (s) reg = i k L * f (k − s) reg in Corollary II. 15, to show that the other |k|/2 special values within the critical strip for negative weight meromorphic modular forms agree with the vanishing results obtained directly from modularity. We leave explicit confirmation of this agreement at s = k, k + 1, . . . , k/2 + 1 (k < 0) to future work.

IV. SUM-RULES IN CONFORMAL FIELD THEORY
In this section, we use the L-functions developed in the previous sections to extract interesting physical data that characterizes conformal field theories in two dimensions (2d CFTs). In particular, we show that the L-function for the logarithmic derivatives of the path integral at s = 0 gives the Casimir energies and thus central charges of certain unitary CFTs. We further show that the special value at s = 1 of the same L-function confirms a sum-rule that was motivated by a recently noticed symmetry of path integrals in quantum field theory (QFT), when applied to 2d CFT path integrals with infinite product expansions.
The structure of this section is as follows. In section IV A, we briefly describe the physical motivation and setting for our study of sum-rules in 2d CFTs. In particular, we emphasize why modularity appears in 2d CFTs [8], discuss holomorphic factorization, and the role of unique ground-states in rewriting path integrals as infinite products. In section IV B, we introduce the stress-energy tensor, and prove a theorem about the special values of its Lfunction when the path integral is a weight-k modular form whose q-series begins with q −∆ . Finally, in sections IV C and IV D we use this L-function to verify the sum-rules suggested by T-reflection, and point-out a class of functions that seem related to traces of singular moduli, but do not seem to be discussed often.

A. Central charges in 2d CFTs and special values of L-functions
By definition, the path integral for a QFT integrates over all allowed configurations of the fields in the theory over all points on the space-time manifold. Carrying-out the path integral for a 2d CFT placed on the two-torus obliterates all information about the structure of the two-torus, save for the lattice of points identified by the toroidal compactification and periodicity conditions along the non-contractible cycles of the two-torus.
The discussion in this section applies to 2d CFT path integrals that factorize holomorphically and have unique ground-states. When the CFT factorizes holomorphically, we have where Z and Z correspond to the path integrals for the decoupled left-and right-movers of the CFT. Holomorphic factorization conventionally means that Z and Z are separately modular invariant, and thus are modular functions: Z, Z ∈ M ! 0 .
However, one of the main motivations of this work is to think of meromorphic modular forms as a possible testing-ground for more general QFT path integrals that behave well under modular transforms and have (Hagedorn) poles. So, in this spirit, we allow ourselves to consider path integrals Z that have nontrivial modular weight: Z ∈ F k for k ∈ 1 2 Z. For convenience, we will call Z a CFT path integral even if it has nontrivial weight or poles.
When the CFT has a unique ground-state, then the lowest power of q in its q-series has unit coefficient, while the other q-series coefficients are all integers: Here −∆ is called the Casimir energy of the CFT. It represents the vacuum-energy of the CFT. The D(n) count the number of states in the full CFT with energy n above the vacuum. As they count states, the D(n) are positive integers. (Twisted indices and path integrals may have D(n) negative.) Because D(0) = 1 and D(n) ∈ Z, we can rewrite Z(τ ) as where the d(n) ∈ Z for every n. (Had D(0) = 1, then the d(n) fail to be integers.) Correlation functions of the stress-energy tensor of the theory, T µν (τ ), carry crucially important information about the theory. We find it very useful to study the expectationvalue of the stress-energy tensor for the full CFT, T (τ ) : 1 Crucially, if Z(τ ) vanishes or has some pole at some finite value of τ , then T (τ ) has a simple pole. For the theories we are interested in, the stress-tensor has the following q-series: This has a number of important consequences. First, by the bounds in section II B, if Z(τ ) has a finite collection of singular points in the fundamental domain, and if the point with largest imaginary part is at τ = z, then there exists a real number C > 0 such that Ce 2πnIm(z) ≤ |md(m)| < e 2πn(Im(z)+ǫ) for every ǫ > 0. (The upper-bound is always satisfied. The lower bound is satisfied for infinitely many positive integers n.) 1 The statement that the one-point function of the stress-energy tensor for the CFT is given by the logderivative of the path integral can be understood in two ways. First, mechanically the ∂ τ log Z gives the expectation value of energy in a statistical ensemble. Second, from general principles the stress-energy tensor is sensitive to variations in length-scales in the geometry of the spacetime manifold. Thus, it is natural that T is given by a functional derivative of the path integral with respect to the shape of the two-torus, τ . Normalizing this derivative by the needed factor of 1/Z yields T ∝ ∂ τ log Z, as in Eq. (4.4).
Second, from the structure of Theorem II.14, we know that the special value of L reg T (s) at s = 0 is exactly the constant term of T (τ ) : the Casimir energy, −∆. This happens despite the growth of the q-series coefficients. From a physics perspective, this is an amusing result.
In unitary 2d CFTs, the Casimir energy ∆ is directly proportional to the central charge c of the CFT. The central charge is commonly referred to as a measure of the degrees of freedom in the theory. Thus, realizing the Casimir energy as the special value of the Lfunction for the stress-energy tensor is logically equivalent to stating that the central charge can be literally interpreted as the sum of the number of excitations in the CFT. This is striking. Even if the CFT is both (a) free CFT and (b) satisfies the rather strong constraint of holomorphic factorization, in general the path integral will vanishes at some value of τ . Hence, by Proposition II.7, the d(n) exhibit exponential/Hagedorn growth. To even write-down the L-function for the stress-energy tensor, we need L-functions for meromorphic modular forms in Theorem II.14. In sections IV B and IV D we develop this physical picture further.

B. The L-function of the one-point function of the stress-energy tensor
In this section, we prove two results about the structure of L reg T (s). (The proofs rest heavily on the beautiful paper [22].) Before doing so, we must introduce a bit of language from statistical mechanics. To begin, we argue that the following two functions, play distinct and important roles when the 2d CFT is (in some sense) free. Here, Z CAN (τ ) is the canonical partition function which counts the number of distinct excitations in the single-particle Fock-space of the theory with a given energy n above the vacuum, while Z GC (τ ) is the grand canonical partition function which counts the number of distinct states in the full multi-particle Fock-space with a given energy n above the vacuum. The physical justification for naming these two functions is as follows. When the quantum field theory is free, then the generating functions for the single-particle Fock-space Z CAN (τ ) and the generating function for the multi-particle Fock-space Z GC (τ ) are related by the following combinatoric map: With Z GC (τ ) and Z CAN (τ ) as defined in Eqs. (4.6) and (4.7), we see that if D(0) = 1 then applying the map in Eq. (4.8), we find: This is exactly the infinite product factorization of Z(τ ) for 2d CFTs in Eqs. (4.2) and (4.3).
Famously, particle number is not fixed in QFT: particle production can happen. As path integrals consider all fluctuations for all possible field configurations of a given QFT or CFT, then, the path integral for a CFT or QFT does not naturally map onto the canonical partition function with fixed particle-number. Yet, because particle number is not fixed in the grand canonical partition function, which accesses the multi-particle Fock-space, it is natural that the path integral and the grand canonical partition function be equal. In what follows, we use Z(τ ) and Z GC (τ ) interchangeably. Now for the main Theorem of this section: where −∆ is the leading power of q in Z GC (τ ), and "regular(s)" refers to the terms in L * f (s) reg that are finite for finite s ∈ C. Proof. Our proof rests on five results. First, that the weighted sum of the orders of zeros and poles of a modular form f ∈ F k in the fundamental domain, F , is exactly given by k: Here ord f (w) is the order of the pole or zero of f at w, and e(w) = 1 unless w = (−1) 1/2 or (−1) 1/3 . At these special values, e((−1) 1/2 ) = 1/2 and e((−1) 1/3 ) = 1/3. Note that ∆ is the order of the zero or pole at the cusp τ = i∞. Second, we use Theorem 5 of [22]. This Theorem states that for f ∈ F k with a q-series q −∆ n c(n)q n with c(1) = 1, and Λ w (τ ) = 1 2πi ∂ τ log(j(τ ) − j(w)), then Third, we use the fact that the L-function of the quasi-modular E 2 (τ ) is, Fourth, we use the fact that the L-function of the exactly modular, but meromorphic, modular form Λ w (τ ) is, where "regular(s)" refers to the terms in L * f (s) reg that are finite for finite s ∈ C. Fifth, we use the fact that the L-function of a finite sum of modular forms is the sum of the finite number of L-functions of each of the individual modular forms. This follows from the fact that L-functions are fundamentally integral transforms of the modular forms and from the fact that an integral is a linear functional of its arguments. Now, we use these five facts for the Borcherds products that define Z GC (τ ). Thus, using Eq. (4.12) we find that Now, using Eqs. (4.13) and (4.14), we have By Eq. (4.11), we see that the common coefficient of the 1/s and 1/(s − 2) poles is simply ∆. This concludes the proof.
Theorem IV.1 concerns special values of L reg T (s). Before stating a very useful Corollary, we make the following observation. When written as a formal Dirichlet series, this L-function would take the form, Recall that in section III B, we promised a physical motivation for the formal factorization of the L-function for Λ d (τ ) into ζ(s) multiplied by n A(n 2 , d)/n s−1 . To give it, we draw an analogy between H d (j(τ )) and the grand canonical partition function of a CFT, n A(n 2 , d)q n and the single-particle partition function of the CFT, and Λ d (τ ) as the stressenergy tensor of the CFT. This analogy provides the physical motivation for the sum-rules in Eq. (3.29).
The sum-rules Eq. (3.29) are special cases of the following quite general Corollary to Theorem IV.1: Corollary IV.2. If a 2d CFT path integral Z(τ ) can be written as q −∆ n (1 −q n ) −d(n) ∈ F k and Z CAN = n d(n)q n , then (4.20) Proof. This immediately follows from Theorem IV.1 and from the formal definition of the L-function for canonical ensemble in Eq. (4.19): L reg Z CAN (s) = −L reg T (s + 1)/ζ(s + 1).

C. The L-function for the stress-tensor and T-reflection sum-rules
The two general sum-rules in Corollary IV.2 were anticipated in Refs. [6,7,14,28], where it was noticed that many finite-temperature path integrals in QFT were invariant under formally reflecting temperatures to negative values (T-reflection). Particularly in Refs. [6,7] it was argued that if a 2d CFT path integral Z(τ ) had well-defined modular weight k ∈ 1 2 Z, and could be written in terms of an infinite product q −∆ n (1 − q n ) −d(n) with d(n) ∈ Z, then T-reflection invariance would imply two the sum-rules: (4.21) We now note that the L-functions defined in section II provide a concrete context in which we may evaluate these sum-rules. Theorem IV.1 and Corollary IV.2, validate these sum-rules. As emphasized above, the stress-energy tensor plays an absolutely crucial role in QFT, and in particular in CFT. By encoding a CFTs central charge, the stress-tensor both counts the number of degrees of freedom in the system, and, more importantly, describes how a CFT responds to deformations of its spacetime manifold (i.e. variations of τ , the shape of the torus). This latter fact is of chief interest to us, here: The response to the CFT to deformations of its spacetime manifold constitutes a sort of anomalous breaking of conformal invariance: If a theory has a nonzero central charge, it has a conformal anomaly.
Similarly, in [7] we argue that 2d CFT path integrals are naturally both invariant under modular transformations and under T-reflection. Throughout [6,7], we argued that if a path integral were invariant under T-reflection up to such an overall phase, this phase would constitute a global gravitational anomaly [29][30][31]. A main role of the stress-energy tensor is to encode the anomalies of a CFT under diffeomorphisms or variations in the spacetime manifold. The new fact in this section, that the L-function of the stress-energy tensor encodes the anomaly of a CFTs path integral under the T-reflection redundancy in how the two-torus is encoded in the path integral, is thus striking.
Not only striking, this new fact is consistent with recent work on T-reflection [6,7] and related work on modular forms [28]. In [7], we argue that the T-reflection phase is tied directly to the modular weight of a 2d CFT path integral: e iγ = (−1) k . The sum-rule in Corollary IV.2 explicitly verifies that this sum-rule is consistent with our L-functions for meromorphic modular forms. Further, in [28] we explicitly construct an extension of SL 2 (Z) modular forms defined on the upper half-plane to GL 2 (Z) modular forms defined on the double half-plane, where the T-reflection phase is again given by (−1) k . This explicitly agrees with the sum-rules in Corollary IV.2, and fits well with the physical and mathematical structure of Refs. [6,7,28]. In this paper, we have shown that the relevant sum-rule that counts the T-reflection phase comes exactly from the L-function for the stress-energy tensor.

D. Borcherds products and the statistical mechanical ensembles
In this brief section, we note that the objects extracted from Borcherds products that we called Z CAN (τ ) may be of independent mathematical interest. We have argued that a Borcherds product resembles the partition function for a 2d CFT on a two-torus in the grand canonical/multi-particle ensemble. The logarithm of a Borcherds product resembles the associated partition function for the same 2d CFT in the canonical/single-particle ensemble.
This physical picture suggests that the canonical partition functions Z CAN (τ ), (4.22) may be of independent mathematical interest. Certainly, for the simplest example, the partition function in the canonical ensemble is a character of SL 2 (R). Juxtaposing the Borcherds products in section III B and in [2] against the physical picture where a modular form is some "grand canonical partition function" for a free statistical system, suggests that the "canonical partition functions" associated to H d (j(τ )), (4.24) are mathematically interesting objects. However, we are unaware of any study of these functions in the literature. It might be pleasing if the other functions for d > 0 also have a clear meaning in terms of characters of non-compact groups, or played a more significant role in the study of e.g. traces of singular moduli. We leave such questions to future work.

V. SUMMARY AND FUTURE DIRECTIONS
Context for the paper: The motivation for this paper was to verify a conjecture of Treflection [14], in the physics-agnostic setting of the mathematics of modular forms [6,7]. To do this, we extended the recent map between weakly holomorphic modular forms [3][4][5] and regularized L-functions/regularized Dirichlet series -itself, an extension of the classic map between holomorphic modular forms and L-functions/Dirichlet series [1] -to a map between meromorphic modular forms and regularized L-functions/regularized Dirichlet series.
Summary and main results of the paper: We have developed a map between meromorphic modular forms and regularized L-functions, presented in Theorem II.14. To do so, we had to explicitly and precisely understand the exponential growth in q-series coefficients, presented in Proposition II.7. We then used these results to understand aspects of traces of singular moduli, when viewed as q-series coefficients of meromorphic modular forms. Many of these aspects hinge critically on L-functions for meromorphic modular forms, and thus are new.
In a seemingly different direction, we viewed meromorphic modular forms as models for CFT path integrals and partition functions in various ensembles in statistical mechanics. Concretely, we mapped meromorphic modular forms onto partition functions for the full multi-particle Fock-space of free 2d CFTs, and we mapped (roughly) the logarithm of meromorphic modular forms onto partition functions for the single-particle Fock-space of a the same free 2d CFTs. We then studied the L-function of the stress-energy tensor, which is the logarithmic derivative of meromorphic modular partition function, of these CFTs, and studied its special values. The stress-tensor is a weight-two meromorphic quasimodular form. Its regularized L-function has special values at s = 0, s = 1 and s = 2.
In Theorem IV.1 and Corollary IV.2, we explicitly showed that the s = 0 special values of the stress-energy tensor L-function allow us to equate the central charge of a free and unitary 2d CFT with a regularized tally of the total number of states in the single-particle Fock-space. Translated back into the mathematical context of traces of singular moduli in, this is equivalent to the statement that the regularized sum of traces of singular moduli for the family of modular functions J n (τ ) := q −m + O(q) ∈ M ! 0 , when evaluated at CM points of discriminant-d quadratics, equals (minus) the Hurwitz class number H(d): reg n t n (d) = −H(d). Though these statements do not have obvious practical utility, explicitly relating the central charge to the (in general) exponentially divergent over states in a CFT and explicitly relating the exponentially divergent sum of the mathematically interesting t n (d) to the Hurwitz class numbers, are both new and pleasing results. Further, the s = 1 special value of the stress-energy tensor L-function verifies the sum-rule for the T-reflection phase, conjectured in [6,7], and hints at important physical aspects of this newly found symmetry.
Bugs, features, and the future: The motivation for this paper is solidly from theoretical physics. The main advances in this paper are largely mathematical in nature. The interplay between these two is somewhere in between. Bugs, features, and future of the results in this paper, and the lines of reasoning that led to it are most easily discussed in the three following brief sections: one on completely mathematical aspects, one on the interface between mathematics and theoretical physics that spurred this note, and one on physical aspects.
A. Mathematical aspects: comments, room for improvement, and exploration In this section, we describe various different routes for improvement in our treatment of L-functions of meromorphic modular forms, make several comments that did not fit into the narrative of sections II and III, and emphasize directions for future exploration.
L-integrals when f ∈ F k has poles at cusps: As emphasized in section II, in order to make the L-integral for meromorphic modular forms well-defined, we made use of contour deformations. It is important that we regulated the integral while keeping the integrand meromorphic. However, as discussed in [4,5], when f ∈ M ! k has poles at cusps, it is natural to regularize the L-integral by deforming the integrand by the non-holomorphic factor e 2πυIm(τ ) (f (τ ) − c f (0)). When Re(υ) is larger than the order of the pole at the cusp, then the υ-deformed L-integral converges, and can be continued to υ = 0. This procedure yields the terms Eq. (2.12). In this sense, our results agree with and generalize the results for L * f (s) reg in the literature. However, when f ∈ F k has a pole at cusps then there are a finite number of terms in the next line (2.13), which diverge.
It is not clear (to me) how to reconcile the tension between non-holomorphic deformations of the integrand from [4,5] to regularize poles at cusps while also being able to use residue theorems and contour-deformations to regularize divergences from poles away from cusps. I will list some possible routes. First, it is possible that there is a sense in which the Λregulated L-integral is given by the constant term in L * f (s|Λ, ǫ) reg , when expanded for large Λ. This would excise the finite number of problematic (divergent) terms in Eq. (2.13). Second, it would be interesting to consider a hybrid regulator, where we look at where J(τ ) = q −1 + n≥1 c(n)q n ∈ M ! 0 , J(τ ) α = e α log J(τ ) and log x is defined to be in the principal branch: −π < Im(log z) ≤ π. With this deformation, the integrand should be meromorphic throughout, allowing straightforward use of residue theorems to ǫ-deform the contour away from poles away from cusps. For sufficiently large Re(α) ≫ 1 the poles at cusps are also regularized. This could in principle allow one to define a regularized integral of a meromorphic integrand that has a smooth limit as Λ → ∞ and ǫ, α → 0. As our chief interest was in meromorphic modular forms that are regular at cusps, we leave this question for future work.
Poles, path-dependence, and wall-crossing: Meromorphic functions have poles. To define their L-functions, we must take a contour integral. Insisting the L-function integrand be meromorphic in τ , rather than e.g. real-analytic function of t = Im(τ ), then it is natural to consider what happens when a pole smoothly moves from z to γz, when this path crosses the contour. If γ ∈ SL 2 (Z), then the locations of set of poles at {γz | γ ∈ SL 2 (Z)} is mapped back to itself. As the integrand has poles at all modular images of z, it is invariant under discrete modular transformations z → γz. However, the integral is not invariant if we continuously deform z to γz if the path crosses the original contour along the imaginary-τ axis. When this happens, then the L-function picks-up residues.
There seems to be tension here, analogous to the tension between modularity of the theta-decomposition of meromorphic Jacobi forms and modularity in Zwegers work [23] and others [24]. There, as here, the non-invariance of the integral seems to be related to pathdependence, which we discussed in section II D. It would be very interesting to understand this phenomena in more detail. (Note: these residues do not alter the residues of 1 Special values in sections III and IV are thus unaffected by this ambiguity.) Polar Rademacher sums and numerical evaluation of L * f (s) reg : It would be very interesting to systematically isolate the contribution of each individual pole to the q-series coefficients for f ∈ F k when k > 0, in a Rademacher-like sum. Hardy and Ramanujan did this for 1/E 6 (τ ) ∈ F −6 [15]. Recent, beautiful, work extends this to wide classes of meromorphic modular forms f ∈ F k with negative weight [16][17][18][19][20].
However, expressions for q-series coefficients for meromorphic forms with k > 0 seem not to be written in terms of the polylogarithms needed to explicitly relate exponential q-series growth due to explicit pole locations in H. If such an extension existed, by simply deleting exponential contributions to c f (n) from poles above the line Im(τ ) = B, we could directly obtain rapidly convergent expressions for c Rf (n|B) used in L * f (s) reg for any f ∈ F k . (Note: q-series coefficients for some f ∈ F k for k ≥ 0 in terms of exponential sums for k > 0 are known, e.g. [32,33]. Yet, the map between terms in the exponential sums and pole locations, which follows naturally from Rademachers sums, seems less direct here.) The constants c f (0) decouple from L-functions for f ∈ F 0 : The integral that defines (0)). So when f ∈ F 0 , then all L-functions are just functions of the non-constant q-series coefficients: the c f (0) drops out. Amusingly, c f (0) is also an element of F 0 . This may suggest a fruitful reformulation of L-functions for f ∈ F k . Recall that when f ∈ M k , we subtract-off the modular function called c f (0) from f (τ ) within L * f (s). This has the effect of subtracting-out the pole in Mellin-space that comes from a residue at τ = i∞. When f ∈ M ! k or f ∈ F k , rather than subtracting-off the modular function called c f (0) which cancels the pole coming from t → 0 and t → ∞, we could subtract-off a modular function that cancels all poles along the integration contour but introduces no new poles. In this scenario, the integrand would be meromorphic in τ throughout, and subtleties of L * f (s) would be contained within L reg f (s) for f ∈ F 0 (and f ∈ M ! 0 ). For future exploration: It would be good to fit L-functions for f ∈ F k (and f ∈ M ! k [4,5]) into the web of conjectures and facts about L-functions for f ∈ M k and S k . It is natural to extend L reg f (s) for f ∈ F k (Γ) to subgroups Γ < SL 2 (Z), or to half-integral k ∈ 1 2 Z. Lfunctions for weakly holomorphic modular forms with Zagier duality particularly stand-out.

B. Field theory, statistical mechanics, Hagedorn, and number theory
There is a strong interrelation between the mathematical and physical perspectives in this paper. This goes in both directions.
Borcherds products as grand canonical partition functions: In section IV, we emphasized that if a modular form has an infinite product of the form F (τ ) := q −∆ n (1−q n ) −d(n) ∈ F k , then it can be interpreted in some sense as the partition function for the multi-particle Fockspace of a free CFT, i.e. a grand canonical partition function. In this guise, then, one can define the partition function for the associated single-particle Fock-space of the free CFT. It is given by f (τ ) := n d(n)q n . (A combinatoric map relates f and F . See section IV B.) Physically, both F and d play very important roles. However, we have been unable to find any discussion of f (τ ) in the mathematics literature, other than for the simple case of F (τ ) = 1/η(τ ), where f (τ ) = 1/(1 − q) is a character of SL 2 (R). It would be very interesting to see if these "single particle partition functions" f (τ ) played a comparably important role in number theory-particularly when F has a Borcherds product expansion, and f (τ ) = − n A(n 2 , d)q n has such a close relationship with traces of singular moduli.
Strings, modularity and Hagedorn: As emphasized in the Introduction, both the lowenergy limit of QCD and the field theory limit of string theories generically exhibit an exponential/Hagedorn rise in the number of states with energy, d(E): When the inverse-temperature β approaches β H , then the one-loop path integral/partition function Z(β) = E d(E)e −Eβ diverges. Additionally, loops of closed strings, either QCDstrings or more fundamental strings, are topologically equivalent to a torus. Thus, we expect the one-loop path integral in these theories, Z(β), to be modular in the β-parameter. Putting these two features together naturally suggests that the theory of meromorphic modular forms and perhaps the related theory of meromorphic Jacobi forms [24] may play a role in understanding the path integral and observables in these physical contexts. In this paper, we studied whether one can define a regularization for the exponentially divergent sums that appear in Casimir energies in these field theories. Ubiquity of modularity and Hagedorn growth/poles suggests L-functions for meromorphic modular forms may directly yield Casimir energies of models of low-energy QCD. It would be very interesting if other aspects of meromorphic modular or Jacobi forms interplay with low-energy QCD.
There is at least one explicit precedent. Namely, recently, it was realized that the path integral of a famously tractable limit of QCD in four-dimensions [9] is given by [10,11] Z(β) = ∞ n=1 (1 − q n ) (1 + q n )(1 − zq n )(1 − q n /z) , z = 2 + √ 3 . This can be easily recognized as the quotient of Jacobi theta-functions and Dedekind etafunctions. It has simple whenever q n → z ±1 . It would be very interesting to understand if, e.g. the wall-crossing in [24] had some role in models of low-energy QCD, for instance in this particular path integral that is meromorphic and involves Jacobi theta functions.
Fermionic symmetries: Finally, we make an amusing standalone observation. Let F (τ ) = q −∆ n (1 − q n ) −d(n) = q −∆ n D(n)q n be a weakly holomorphic modular form with zeros at finite values of τ and poles at cusps. Its q-series coefficients are bounded by |D(n)| < e C √ n for some C > 0. Because F (z) = 0 for some 0 < |z| < ∞, then Proposition II.7 implies |d(n)| > E(e 2πIm(z) ) n /n for some E > 0 for infinitely many positive integers n.
Thus, there are huge cancellations between the q-series coefficients of f (τ ) = n d(n)q n and F (τ ) = q −∆ n D(n)q n . Very similar cancellations were observed in [34][35][36][37] from a fermionic symmetry in non-supersymmetric models of QCD with Hagedorn growth in the number of states at a given energy E. Despite the lack of supersymmetry and the presence of Hagedorn growth, the bosons and fermions are almost exactly paired and cancel in the twisted path integral, which avoids Hagedorn poles at finite temperature.
C. Central charges, T-reflection phases, and the stress-tensor Finally, in sections IV B and IV C we showed that both the Casimir energy and the Treflection sum-rule of a 2d CFT are captured by the L-function of the stress-energy tensor, L reg T (s), when evaluated respectively at s = 0 and s = 1. Central charges give the conformal anomaly of the CFT when it is on curved manifolds. So it is natural that they are captured by the stress-tensor. However, it was completely unexpected that the T-reflection sum-rules suggested in [6,14] would literally appear next to the central charge, in the stress-energy tensor L-function, at s = 0 and 1. Yet, this may be sensible in light of Refs. [6,7], where we noted that the T-reflection phase has a natural interpretation as a global gravitational anomaly. It seems existentially important to better understand the connection between the T-reflection phases and the regularized L-function for the stress-energy tensor of 2d CFTs.
As there is only one scalar, and only one direction in which to move, we find that d(n) = 1 for n = 0, 1, 2 . . . =⇒ Z CAN (τ ) = ∞ n=0 d(n)q n = 1 1 − q . (A1) If we apply the combinatoric map between the single-particle and multi-particle Fock-spaces for this free CFT, we find that Z GC (τ ) evaluates to Now, if we take the logarithmic derivative of Z GC (τ ) to find the one-point correlation function of this 2d CFTs stress-energy tensor, we find it equals the following expression: It is straightforward to compute that the Mellin transform, i.e. the L-function, for the one-point correlation function for the stress-energy tensor of this 2d CFT is Because T (τ ) neither diverges nor vanishes at any finite value of τ = β/(2πiL), the regularization procedure developed in this paper is not necessary to define L f (s). We can also directly compute the L-function of Z CAN (τ ), and find that it equals It is important to note that the formal factorization in Eq. (4.19) would imply that L Z CAN (s) := −L T (s + 1)/ζ(s + 1). Crucially, these computations match. (Note that T (τ ) = −E 2 (τ ) + O(q 0 ), and is tightly related to holomorphic quasimodular forms. See Eq. (4.13) for an explicit expression for the L-function of E 2 (τ ).) Finally, we would like to formally motivate the T-reflection sum-rules in section IV C for this special case. First, note that Z GC (τ ) resembles a product of decoupled products of harmonic oscillator partition functions, 1/(1 − q n ), stripped of their zero-point energies q n/2 : n Z (n) osc. (τ ) := n q n/2 1 − q n . (A6) We may recover the grand canonical partition function by taking a large (infinite) collection of decoupled oscillators that have been properly endowed with their zero-point energies q n/2 . Now, note that if we regulate the sum of exponents n n/2 by the L-function for the canonical ensemble, i.e. ζ(s), then we would find .
Further, note that each decoupled oscillator is invariant under T-reflection [6,14]. (This follows from the fact that q n/2 /(1 − q n ) equals sin(πnτ ) −1 , an odd function of τ .) Thus: n Z (n) osc. (β) = n q n/2 1 − q n → n Z (n) osc. (−τ ) = n q −n/2 1 − q −n = n (−1) q n/2 1 − q n . (A8) Again, we take the limit where the number of oscillators goes to infinity, which is needed in order for the states accessed in the partition function to span the full multi-particle Fockspace, but now at negative temperature. Here, we must regulate the divergent product of zero-point energies and the divergent product of (−1)-factors. Regulating with ζ(s), we find: Thus, we see that the leading power of q in the q-series of this particular path integral is fixed by T-reflection. Further, the eigenvalue under T-reflection is e iγ T = −i, consistent with the sum-rule in Corollary IV.2 and with Refs. [6,7,14,28], where e iγ = (−1) k = (−1) reg n d(n) .