1 Introduction
Let
$\mathbb {S}^{2n-1}$
denote the unit sphere in
$\mathbb {C}^n$
, where
$n\geq 2$
. The hypersurface
$\mathbb {S}^{2n-1}$
is an embedded CR manifold, and the tangential Cauchy–Riemann operators
$\overline {\partial }_b$
and
$\overline {\partial }_b^*$
are defined on the corresponding Hilbert spaces. Furthermore, the Kohn Laplacian on
$L^2(\mathbb {S}^{2n-1})$

is a linear self-adjoint densely defined closed operator. We refer to [Reference BoggessBog91, Reference Chen and ShawCS01] for the detailed definitions.
Our inspiration is the celebrated Weyl’s law for Riemannian manifolds. In that setting, the eigenvalue counting function of the Laplace–Beltrami operator on a manifold M has leading coefficient proportional to the volume of M with a constant depending only on the dimension of M. There is also a corresponding result for the extension of this operator (also called the Hodge Laplacian, or the Laplace–de Rham operator) to differential forms.
Motivated by the spectral theory for the Laplacian on Riemannian manifolds, one can investigate the spectrum of
$\square _b$
on CR manifolds and its relation to the complex geometry of the underlying manifold. For example, in [Reference FuFu05, Reference FuFu08], Fu studied the spectrum of the
$\overline {\partial }$
-Neumann Laplacian
$\square $
on smooth pseudoconvex domains
$\Omega $
and related the distribution of the eigenvalues of
$\square $
to the D’Angelo type of
$b\Omega $
. In [Reference FollandFol72], Folland computed the spectrum of
$\square _b$
on
$\mathbb {S}^{2n-1}$
on all differential form levels. Recently, in [Reference Ahn, Bansil, Brown, Cardin and ZeytuncuABB+19, Reference Abbas, Brown, Ramasami and ZeytuncuABRZ19], the authors used the spectrum of the Kohn Laplacian to prove the nonembeddability of the Rossi sphere.
For
$\lambda>0$
, let
$N(\lambda )$
denote the number of positive eigenvalues (counting multiplicity) of
$\square _b$
on
$L^2(\mathbb {S}^{2n-1})$
that are less than or equal to
$\lambda $
. It was noted in [Reference Ahn, Bansil, Brown, Cardin and ZeytuncuABB+19] that
$N(\lambda )$
grows on the order of
$\lambda ^n$
and later, in [Reference Bansil and ZeytuncuBZ20] (see the erratum), the leading coefficient in the asymptotic expansion was calculated as an infinite sum. In particular, [Reference Bansil and ZeytuncuBZ20] obtained

by using careful counting arguments.
1.1 Main result
In this paper, we obtain the same leading coefficient by another argument, namely by Karamata’s Tauberian Theorem. We highlight that this technique is different than the one in [Reference Bansil and ZeytuncuBZ20]. Furthermore, we formulate the leading coefficient as a product of the volume of
$\mathbb {S}^{2n-1}$
and an integral that only depends on n. This representation resonates better with Weyl’s law and is more amenable to generalization to other CR manifolds.
Theorem 1.1 Let
$N(\lambda )$
be the eigenvalue counting function for
$\square _b$
on
$L^2(\mathbb {S}^{2n-1})$
as above. Then,

The integral representation is strikingly similar to (and indeed inspired by) a similar formula for the leading coefficient of the eigenvalue counting function for
$\square _b$
acting on
$(0,q)$
-forms (where
$q\geq 1$
) in [Reference Stanton and TartakoffST84, Reference StantonSta84]. However, we note that the expression in Theorem 1.1 is not a special case of the expression in [Reference Stanton and TartakoffST84]. Stanton and Tartakoff’s formula does not cover the case of functions; the improper integral in their representation diverges for
$q=0$
. Furthermore, in the case of functions, the kernel of
$\square _b$
is infinite-dimensional (unlike
$(0,q)$
-forms where
$q\geq 1$
). Therefore, a careful regularization idea and some nontrivial adjustments are necessary to formulate the correct expression for the case
$q=0$
.
We anticipate that, for a general pseudoconvex CR manifold M of hypersurface type embedded in
$\mathbb {C}^n$
, the leading coefficient of the Kohn Laplacian on functions should be given by

This conjecture is modeled after Theorem 1.1 and is further supported by a connection with Stanton and Tartakoff’s formula. In the last section, we recover formula (1.1) by extending Stanton and Tartakoff’s formula to a meromorphic function and evaluating the nonsingular part at
$q=0$
. Although in Section 3 we present only the computation for the case of spheres, the same computation holds for general hypersurfaces in
$\mathbb {C}^n$
.
1.2 Ingredients
Before we present a proof of Theorem 1.1 in the next section, we state some known facts about the eigenvalues of
$\square _b$
on
$\mathbb {S}^{2n-1}$
and Karamata’s Tauberian Theorem.
The Kohn Laplacian acts on the space of
$L^2$
differential forms on the sphere. In [Reference FollandFol72], Folland uses unitary representations to explicitly compute the eigenvalues and corresponding eigenspaces for the Kohn Laplacian on
$(0,j)$
-forms. In particular, he shows that eigenforms

have corresponding eigenvalues
$2 (q+j) (p+n-j-1)$
, where the hat over a form indicates its exclusion from the wedge product.
We are interested in the case where
$j=0$
, corresponding to functions on the sphere. Folland (see also [Reference Ahn, Bansil, Brown, Cardin and ZeytuncuABB+19, Reference Abbas, Brown, Ramasami and ZeytuncuABRZ19]) explicitly shows that, in this case,

where
$\mathcal {H}_{p,q}(\mathbb {S}^{2n-1})$
is the space of spherical harmonics of bidegree
$p,q$
. Furthermore, each space
$\mathcal {H}_{p,q}(\mathbb {S}^{2n-1})$
is an eigenspace of
$\square _b$
with dimension

as computed by an inclusion–exclusion argument (see [Reference KlimaKli04]). These spaces have corresponding eigenvalues
$2q (p+n-1)$
.
Karamata’s Tauberian Theorem has been used in [Reference KacKac66, Reference Stanton and TartakoffST84] to understand the distribution of eigenvalues. We follow the statement in [Reference Arendt, Nittka, Peter, Steiner, Arendt and SchleichANPS09, Theorem 1.1, p. 57], where the reader can find further references.
Theorem 1.2 (Karamata)
Let
$\{\lambda _j\}_{j\in \mathbb {N}}$
be a sequence of positive real numbers such that
$\sum _{j\in \mathbb {N}}e^{-\lambda _j t}$
converges for every
$t>0$
. Then, for
$n>0$
and
$a\in \mathbb {R}$
, the following are equivalent.
-
(1)
$\lim _{t\rightarrow {0}^+}t^n\sum _{j\in \mathbb {N}}e^{-\lambda _jt}=a$ ;
-
(2)
$\lim _{\lambda \rightarrow {\infty }}\frac {N(\lambda )}{\lambda ^n}=\frac {a}{\Gamma (n+1)}$ ,
where
$N(\lambda )=\#\{\lambda _j:\lambda _j\leq \lambda \}$
is the counting function.
The next section is dedicated to the proof of Theorem 1.1. First, we prove the expression for the leading coefficient as a series (recovering the result in [Reference Bansil and ZeytuncuBZ20]). Next, we express the leading coefficient as the volume of
$\mathbb {S}^{2n-1}$
times an integral.
2 Tonelli’s tactic to tie together two Tauberian terms
Instead of considering
$N(\lambda )$
directly, we define the function
$G(t)=\sum _{j\in \mathbb {N}} e^{-\lambda _j t}$
and consider

where the sequence
$\{\lambda _j\}_{j\in \mathbb {N}}$
is the sequence of all positive eigenvalues (with multiplicity) of
$\square _b$
on the sphere. Once we compute this limit, we can invoke Karamata’s Theorem. Using the preliminaries from Section 1.2, we have

Applying the standard recursive formula for binomial coefficients to the first product of binomial coefficients gives usFootnote 1

This allows us to rewrite
$G(t)$
as a sum of two positive pieces by noticing that

We label the parts as

and

so that
$G(t)= G_1(t)+ G_2(t)$
.
Our goal is to calculate
$\lim _{t \to 0^+} t^n G(t)$
, which we do by calculating
$ \lim _{t \to 0^+} t^n G_1(t)$
and
$\lim _{t \to 0^+} t^n G_2(t)$
separately.
Theorem 1.1 claims that
$\lim _{\lambda \to \infty } \frac {N(\lambda )}{\lambda ^n}$
can be written either as an infinite series or as an improper integral. The key to obtaining these two distinct forms is Tonelli’s Theorem (see [Reference RudinRud87, Theorem 8.8]). Because the summands in
$G_1(t)$
and
$G_2(t)$
are positive, we can exchange the order of the infinite sums in each. When the outer sum is over q, we can show that
$\lim _{t\to 0^+} G_1(t)$
is an infinite series and
$\lim _{t\to 0^+} G_2(t)$
is an improper integral; exchanging the order of summation allows us to express
$\lim _{t\to 0^+} G_1(t)$
as an integral and
$\lim _{t\to 0^+} G_2(t)$
as an infinite series.
2.1 Serious series for spherical spectra
In this part, we prove the infinite series formula for
$\lim _{\lambda \to \infty } \frac {N(\lambda )}{\lambda ^n}$
. The computations for
$G_1(t)$
and
$G_2(t)$
are quite similar. In each instance, we will apply the Dominated Convergence Theorem by leveraging Lemma 2.2.
In each computation, we eventually exchange the limit as
$t\to 0^+$
with a sum, and in each case, we require the following limit calculation, which follows from L’Hopital’s rule.
Lemma 2.1 Let
$\alpha> 0$
. Then,

In order to exchange limits with sums, we will apply the Dominated Convergence Theorem. The next technical lemma will come in handy for showing the conditions of the Dominated Convergence Theorem are satisfied in each case.
Lemma 2.2 For
$n\in \mathbb {N}$
, there exists
$M>0$
such that

for all
$x\in (0,\infty )$
.
Proof The function f is continuous on
$(0,\infty )$
. If we show that
$\lim _{x\to \infty } f(x)<\infty $
and
$\lim _{x\to 0^+} f(x)<\infty $
, it will follow that there exists some M such that
$f(x)<M$
on
$(0,\infty )$
. Indeed,
$\lim _{x\to \infty } f(x) = 0$
. Furthermore, L’Hopital’s rule gives
$\lim _{x\to 0^+}\frac {x}{e^{2x}-1}=\frac {1}{2}$
and so
$\lim _{x\to 0^+}\frac {x^n e^{2x(n-1)}}{(e^{2x}-1)^n} =\frac {1}{2^n}$
. ▪
Now, we are ready to compute
$\lim _{t\to 0^+}t^n G_1(t)$
.
Proposition 2.3

Proof We start by coming up with an expression for
$G_1(t)$
containing only a single sum. For
$|z| < 1$
, there is the Taylor series expansion

Let
$z = e^{-2tq}$
; then,

For positive t,
$\frac {t^n(e^{2 t q})}{(e^{2 t q}-1)^n}\leq \frac {1}{q^n}\frac {(qt)^n(e^{2 t q (n-1)})}{(e^{2 t q}-1)^n} =\frac {f(qt)}{q^n}\leq \frac {M}{q^n}$
, where f is as in Lemma 2.2. Hence,
$t^n \binom {n+q-2}{q} \frac {(e^{2 t q})}{(e^{2 t q}-1)^n}\leq M \binom {n+q-2}{q} \frac {1}{q^n}$
, for all t, and because
$\sum _{q=1}^\infty M \binom {n+q-2}{q} \frac {1}{q^n}$
converges, we are free to apply the Dominated Convergence Theorem to
$\lim _{t\to 0^+}t^n G_1(t)$
.
This technical justification allows us to exchange the order of the limit and summation in
$\lim _{t\to 0^+}t^n G_1(t)$
and conclude that

where the final equality follows from Lemma 2.1. ▪
Next, we move to the second piece and compute
$\lim _{t\to 0^+} t^n G_2(t)$
.
Proposition 2.4

Proof We start out by manipulating the form of
$G_2(t)$
as we did with
$G_1(t)$
, but this time, we apply Tonelli’s Theorem to switch the order of summation. In our calculation, we will make use of the substitutions
$z=e^{-2t(p+n-1)}$
and
$w=p+n-1$
, and we will apply the power series expansion of
$\frac {1}{(1-z)^{n-1}}$
.

By Lemma 2.2,
$t^n \binom {w-1}{n-2} \frac {e^{2tw(n-1)}}{(e^{2tw}-1)^n}=\binom {w-1}{n-2}\frac {f(tw)}{w^n} \leq M\binom {w-1}{n-2}\frac {1}{w^n}.$
As
$\sum_{w=0}^\infty M \binom {w-1}{n-2}\frac {1}{w^n}.$
converges, we apply the Dominated Convergence Theorem to
$\lim _{t\to 0^+} t^nG_2(t)$
in order to exchange the limit and the sum. Thus,

▪
Applying the Tauberian Theorem 1.2 in combination with Propositions 2.3 and 2.4 proves that

This concludes the proof of the first identity in Theorem 1.1.
2.2 Is it not interesting: intense integral is identical in immensity
In this subsection, we will prove the second part of Theorem 1.1. We start off by reanalyzing
$G_1(t)$
:

where the switching of the order of summation is justified by Tonelli’s Theorem and we have taken
$w = n + p - 1$
. To calculate the limit as
$t\to 0^+$
, we need the following lemma.
Lemma 2.5 Fix an integer
$r \geq 1$
and consider the function

defined for
$x> 0$
. Then,
-
(1)
$f_{r}(x)> 0$ .
-
(2)
$f_{r}^{\prime }(x) < 0$ for sufficiently large x.
-
(3)
$f_{r}$ is bounded and
$\int _{0}^\infty f_{r}(x) dx < \infty .$
Proof We can tell
$f_{r}(x)>0$
because
$e^{-2x}<1$
, so
$0<1-e^{-2x}<1$
, and hence
$\frac {1}{(1-e^{-2x})^r}>1$
, proving (1).
To establish (2), a derivative calculation shows that

Define
$l(x) = e^{2x} - 1 - (1 - e^{-2x})^{r+1}e^{2x} - 2x.$
To prove (2), it suffices to show that
$l(x) < 0$
for sufficiently large x. For this, let
$y = e^{-2x}$
and note that

Looking at the definition of l, the first three terms converge to r and the last goes to
$-\infty $
, so
$l(x) < 0$
for sufficiently large x.
To prove (3), we will show that
$f_r(x)$
is bounded by a constant near 0 and bounded by an exponentially decaying function for large x. The function
$f_{r}(x)$
is not defined at
$x=0$
, but Lemma 2.1 implies that
$\lim _{x\to 0} f_{r}(x)<\infty $
, so the discontinuity at 0 is removable. This implies that
$f_r(x)$
is bounded on
$[0,\frac {\log (2)}{2}]$
. To bound
$f_r(x)$
on
$[\frac {\log 2}{2},\infty )$
, we note

Using the formula for a difference of rth powers, we obtain the following expression for the numerator of the above fraction:

Hence,

Note
$e^{2x}-1 \geq \frac {e^{2x}}{2}$
for
$x \geq \frac {\log (2)}{2}$
, so

Applying this to
$f_r(x)$
gives that when
$x \geq \frac {\log (2)}{2},$

Thus,
$f_r(x)$
is dominated by an exponentially decaying function for sufficiently large x, so
$\int _0^\infty f_r(x)\,dx<\infty .$
▪
It will be helpful to make the following definition.
Definition 2.6 For real
$\alpha \neq 0$
, define the scaled ceiling function
$\lceil \cdot \rceil _{\alpha }: \mathbb {R} \to \mathbb {R}$
by

For example,
$\lceil 7 \rceil _{3} = 3 \lceil \frac {7}{3} \rceil = 3 \cdot 3 = 9$
and
$\lceil 6 \rceil _{2} = 2 \lceil \frac {6}{2} \rceil = 2 \cdot 3 = 6$
. Then, we have the following proposition.
Proposition 2.7
$\lceil x \rceil _{\alpha }$
is x rounded up to the nearest integral multiple of
$\alpha $
, i.e.,

Proof

and the result follows. ▪
The scaled ceiling function has the following properties:
-
(1) Fix
$x \in \mathbb {R}$ ,
$\alpha> 0$ . Then,
$0 \leq \lceil x \rceil _{\alpha } - x < \alpha $ .
-
(2) Fix
$x \in \mathbb {R}$ . Then,
$\lim _{\alpha \to 0^+} \lceil x \rceil _{\alpha } = x$ .
-
(3) Let
$f: [a,b] \to \mathbb {R}$ be monotonically decreasing. Fix
$0 < \alpha \leq b - a$ . Then, for all
$x \in [a,b-\alpha ]$ ,
$f(\lceil x \rceil _{\alpha }) \leq f(x)$ .
The next few propositions will allow us to compute
$\lim _{t\to 0^+} t^n G_1(t)$
.
Proposition 2.8

Proof Manipulating the sum, we have

Fix C such that
$f_{n-1}^{\prime }(x) < 0$
for
$x \geq C$
and fix M such that
$f_{n-1}(x)<M$
. Then,

for all
$x> 0$
. Because the integral of the right-hand side is finite, we may apply dominated convergence to see that

which completes the proof. ▪
The following proposition will help us compute the rest of the limits we need before we can tackle
$\lim _{t\to 0^+}t^nG(t)$
in its entirety.
Proposition 2.9 Suppose that
$a_t(w)$
is a positive function of real t and integer
$w\geq 1$
such that
$\lim _{t\to 0^+}a_t(w)=0$
for each w and
$\lim _{t\to 0^+}\sum _{w=1}^\infty a_t(w) =M<\infty $
. Then,

Proof Let
$\epsilon>0$
be arbitrary and let
$k>2\frac {M}{\epsilon }-1$
be an integer. Because
$\lim _{t\to 0^+} a_t(w)=0$
for each w, there exists some T such that for all
$0<t<T$
and all
$w\leq k$
,
$a_t(w)\leq \frac {\epsilon }{2k}$
. Thus, for
$t<T$
, we have

Therefore,
$\lim _{t\to 0^+} \sum _{w=1}^\infty \frac {a_t(w)}{w} = 0.$
▪
Proposition 2.10 Fix
$0 \leq k < n -1$
. Then,

Proof The observation that
$0\leq \sum _{w=n-1}^{\infty }t^n w^k ( \frac {1}{(1 - e^{-2tw})^{n-1}} - 1 )\leq \sum _{w=n-1}^{\infty }\frac {1}{w} t^n w^{n-1} ( \frac {1}{(1 - e^{-2tw})^{n-1}} - 1 )$
allows us to apply Propositions 2.8 and 2.9 and arrive at the desired conclusion.▪
With Propositions 2.8 and 2.10 in hand, we can move forward. We break
$\binom {w}{n-1}$
into a polynomial in w as follows:

We write


The final equality was obtained by applying Propositions 2.8 and 2.10.
Now, we move on to
$G_2(t)$
. We have


The next lemma plays a role in our analysis of
$G_2(t)$
, which is analogous to the role of Lemma 2.5 in our analysis of
$G_1(t)$
.
Lemma 2.11 Fix integers
$r \geq 0$
and
$s \geq 1$
, and consider the function

defined for
$x> 0$
. Then:
-
(1)
$g_{r,s}(x)> 0.$
-
(2)
$g^{\prime }_{r,r} (x)< 0.$
-
(3) If
$r \geq s$ ,
$g_{r,s}(x)$ is bounded and
$\int _{0}^\infty g_{r,s}(x) dx < \infty $ .
-
(4) If
$r \geq s + 1$ ,
$\lim _{x \to 0} g_{r,s}(x) = 0$ .
Proof Suppose
$r\geq 0$
and
$s\geq 1$
are integers and that
$x>0$
. As
$x^r>0$
,
$e^{-2xs}>0$
, and
$1 - e^{-2x}> 0$
, we have
$g_{r,s}(x)> 0$
, which shows (1). To prove (2), notice that we may write

Because the function
$x/(e^{2x}-1)$
is decreasing, it follows that
$g^{\prime }_{r,r}(x)<0$
.
For part (3), fix
$M \in \mathbb {R}$
such that for all
$x \geq M$
,
$e^{2x} - 1 \geq \frac {e^{2x}}{2}$
and
$(e^{2x})^{1/2} \geq x^r$
. Then, for
$x \geq M$
,

As
$s \geq 1$
, the right-hand side has finite integral over
$[M,\infty )$
and hence
$\int _M^\infty g_{r,s}(x)\ dx < \infty $
. The integral of
$g_{r,s}(x)$
over
$[0,M]$
is also finite, because we can extend
$g_{r,s}(x)$
to a continuous, bounded function on this compact interval. Adding these two parts together shows that

It remains to show
$\lim _{x\rightarrow 0}g_{r,s}(x)=0$
whenever
$r\geq s+1$
. First, notice that Lemma 2.1 can be used to show that

Therefore, when
$r\geq s+1$
, we have

as we needed to show, this proves (4). ▪
Now, we move on to propositions which, similarly to our strategy for analyzing
$G_1(t)$
, allow us to break up the binomial coefficient in
$G_2(t)$
into a polynomial and analyze the resulting sums separately.
Proposition 2.12 Fix
$0 \leq k < n -1$
. Then,

Proof We recognize that this expression can be rewritten in terms of
$g_{n,n-1}(tq)$
:

Because
$g_{n,n-1}$
is bounded and
$n-k \geq 2$
, the summand is dominated by some constant times
$\frac {1}{q^{2}}$
. Therefore, we can apply dominated convergence to conclude that the limit is 0, because the summand converges pointwise to 0 by the previous lemma. ▪
Proposition 2.13

Proof Write

The integrand is exactly
$g_{n-1,n-1}(\lceil x \rceil _{t})$
, and is thus dominated by
$g_{n-1,n-1}$
, because
$g_{n-1,n-1}$
is a decreasing function. Therefore, we can apply dominated convergence to find that

▪
We can expand
$\binom {n+q-2}{n-1}$
as

Therefore,

Next, we put the two parts of G back together in the limit. This gives us

We now have an expression for
$\lim _{t\to 0^+} G(t)$
in terms of an integral, so we could have chosen to stop here. Instead, we will press on and apply a few more tricks in order to arrive at the expression in Theorem 1.1. One reason we prefer the expression in Theorem 1.1 is its similarity to a closely related result in [Reference Stanton and TartakoffST84].
To continue on our path of manipulating the combined integral, we apply integration by parts with

and
$v=n^{-1}x^n$
. This gives us

The boundary term clearly vanishes at
$\infty $
, because
$x^n \left ( \frac {1}{(1 - e^{-2x})^{n-1}} - 1 \right )$
and
$\frac {x^n}{(e^{2x} - 1)^{n-1}}$
each vanish at
$\infty $
. To see that it vanishes at 0, apply L’Hopital’s rule to the quotients
$\frac {x}{1-e^{-2x}}$
and
$\frac {x}{e^{2x}-1}$
. Because the integrand is even, we further have

As
$\frac {x}{\sinh (x)}$
is an even function, the value of the first integral is unchanged if we change the
$e^{x(n-2)}$
in the integrand to
$e^{-x(n-2)}$
. Thus,

Therefore, by the Tauberian Theorem due to Karamata, we obtain the following limit:

and complete the proof of Theorem 1.1.
3 The formula for functions versus the formula for forms
In [Reference Stanton and TartakoffST84], Stanton and Tartakoff prove the following formula, reminiscent of Weyl’s law, for the Kohn Laplacian on CR manifolds of hypersurface type.
Theorem 3.1 [Reference Stanton and TartakoffST84, Theorem 6.1]
Let M be a CR submanifold of
$\mathbb {C}^n$
,
$n \geq 3$
. Let
$N(\lambda )$
be the eigenvalue counting function of
$\square _b$
on M acting on
$(p,q)$
-forms,
$0 \leq p < n$
,
$0 < q < n-1$
. Then, we have the asymptotic equivalence

where

and the Kohn Laplacian and the volume of M are defined with respect to a Levi metric.
For spheres embedded in
$\mathbb {C}^n$
, the induced metric is a Levi metric (see Definition 1.5 and the following remark in [Reference Stanton and TartakoffST84]). However, as stated, this only applies to
$(p,q)$
-forms with
$q\geq 1$
. We analyze how this expression relates to Theorem 1.1. Toward this end, we define the function

which, for
$q=1,\dots ,n-2$
, is the leading coefficient on
$\lambda ^{n}$
for the asymptotic growth of
$N(\lambda )$
, the eigenvalue counting function of
$\square _b$
on M acting on
$(0,q)$
-forms.
The following statement shows that this function is closely related to our formula.
Theorem 3.2 The definition of f given above is convergent for complex q satisfying
$0 < \Re (q) < n-1$
. Furthermore, f is holomorphic on this strip, and has an analytic continuation to a meromorphic function on the strip
$-1 < \Re (q) < n-1$
whose only pole is at
$q=0$
. Finally, the Laurent expansion of f about 0 is

where
$a_n \neq 0$
.
In other words, the constant term in the Laurent expansion of f about 0 is the expression from Theorem 1.1.
Proof For notational convenience, let
$m = n - 1$
. Then, we may write f as

We first prove that the integrand is integrable (
$L^1$
) whenever
$0 < \Re (q) < m$
. For this, choose
$C> 0$
and
$\alpha> 0$
, so that whenever
$|x| \geq C$
,
$|\sinh (x)| \geq \alpha e^{|x|}$
. Therefore, for
$|x| \geq C$
, we have

Because
$x/\sinh (x)$
has a removable singularity at 0, it is continuous and thus bounded on
$[-C,C].$
Hence,
$x^m/\sinh (x)^m$
is bounded on
$[-C,C]$
as well, so we can choose D, so that
$x^m/\sinh (x)^m \leq D$
whenever
$|x| \leq C$
. Thus, we have the bound

for all
$x \in \mathbb {R}$
.
Therefore, we may write

A repeated integration by parts shows that these integrals are finite if
$0 < \Re (q) < m.$
To define the binomial coefficient, use the gamma function, i.e.,

which is defined because
$\Gamma $
has no zeros. Thus, f is well defined on its domain of definition.
Now, rewrite f as follows:

As
$\frac {x}{\sinh (x)}$
is an even function, the value of the second integral is unchanged if we change the
$e^{-(m-2q)\tau }$
in the integrand to
$e^{(m-2q)\tau }$
. Therefore,

where the last step follows because the integrand is even. Now, define the function g as follows:

We claim that g is holomorphic on the strip
$-1 < \Re (q) < m.$
Assuming this for now, consider

defined for
$0 < \Re (q) < m$
. This is easy to evaluate explicitly with the substitution
$u = 2q\tau$
:

so

This is meromorphic in q on the whole complex plane. Thus, the function
$g + (f - g)$
is meromorphic in q for
$-1 < \Re (q) < n-1$
, and equal to f if
$0 < \Re (q) < n-1$
. Because f converges on its domain of definition, f is holomorphic in q, and
$g + (f -g)$
is the desired continuation. To complete the proof, we need to show that
$g(0)$
is our formula. We have

where we have used the expression appearing in the discussion after Proposition 2.13, and Karamata’s Theorem. Thus, we have our theorem, modulo showing that g is holomorphic for
$-1 < \Re (q) < m$
. We move on to this now. It clearly suffices to show that the expression

is holomorphic for
$-1 < \Re (q) < m$
.
As a stepping stone toward proving that
$h(q)$
is holomorphic, we will first prove that, for
$\beta>0$
,

is convergent and continuous. By differentiating with respect to
$\beta $
, we see that the integrand is motonically increasing, for each
$\tau $
, as
$\beta \to 0.$
Fix some
$\beta _0> 0$
. Now, for
$\tau $
near 0, the integrand is bounded because it has a limit at 0. For large
$\tau $
, it is bounded by a constant times
$e^{-2 \beta \tau } \tau ^m$
, which is integrable on
$[0,\infty ).$
To see continuity at
$\beta _0$
, fix some
$\beta _1$
with
$0 < \beta _1 < \beta _0$
, and note that the integrand at
$\beta _1$
dominates the integrand at
$\beta $
for all
$\beta> \beta _1$
. Continuity of
$\phi $
at
$\beta _0$
follows then by dominated convergence. Thus, we have shown that
$\phi $
is convergent and continuous for
$\beta>0$
.
Next, we compute



To analyze the first term, we have

For the second term, we have

We bound the latter expression in the integrand by

Thus, for the second term, we have

In total, we have shown that

Using this bound, we now show that h is holomorphic via Morera’s Theorem. Fix some triangle
$\Delta \subset \{q \in \mathbb {C}: -1 < \Re (q) < m\}.$
Parameterize
$\Delta $
by arc length with the piecewise differentiable curve
$\gamma (t)$
,
$a \leq t \leq b$
(so
$|\gamma '(t)| = 1 )$
for all t). Then,

The estimate above,

is uniformly bounded in t by the compactness of
$\Delta $
and the continuity of the upper bound in t, so

Therefore, we may apply Fubini’s Theorem to see that


by Cauchy’s Theorem. This shows (by Morera’s Theorem) that

is a holomorphic function of q for
$-1 < \Re (q) < m$
, and thus
$g(q)$
is holomorphic for
$-1 < \Re (q) < m$
, completing the proof. ▪
Acknowledgment
We would like to thank the referees for constructive feedback. We would also like to thank Mohit Bansil and Michael Dabkowski for careful comments on an earlier version of this paper. This research was conducted at the NSF REU Site (DMS-1950102 and DMS-1659203) in Mathematical Analysis and Applications at the University of Michigan–Dearborn. We would also like to thank the National Science Foundation, National Security Agency, and University of Michigan–Dearborn for their support.