Specify the positive spectrum
For integer d≥2, consider one eigenvalue per positive integer index,
with no additional multiplicity:
Pd(n)=S(n,d)=(dn+d−1)=d!n(n+1)⋯(n+d−1),Zd(s)=n=1∑∞Pd(n)−s.
Comparison with n−ds proves absolute convergence for
Res>1/d. All eigenvalues in this definition are
dimensionless. The integer d labels the polynomial sequence; interpreting
an actual spatial operator requires its eigenvalues and multiplicities.
The associated positive heat sum is developed in
Simplex heat trace.
Continuation by subtracting the large-index terms
Set
Ad=2d(d−1),bn=j=0∏d−1(1+nj).
Uniformly for s in a compact set,
bn−s=1−sAd/n+O(n−2). Therefore
Rd(s)=n=1∑∞n−ds[bn−s−1+nsAd]
converges locally uniformly and defines a holomorphic function for
Res>−1/d. Initially in the convergent half-plane,
then by continuation,
Zd(s)=(d!)s[ζ(ds)−sAdζ(ds+1)+Rd(s)].
Here ζ is the Riemann zeta function. The remainder bound follows
by expanding logbn=Ad/n+O(n−2) and its exponential; the tail
is O(n−dRes−2). This also justifies local termwise
differentiation.
In this half-plane, the pole at s=1/d has residue
(d!)1/d/d. At s=0, the apparent pole in the second term cancels
against its factor s.
Value and derivative at zero
Use ζ(0)=−1/2 and
ζ(1+ds)=1/(ds)+γ+O(s), where γ is Euler's
constant. Since Rd(0)=0,
Zd(0)=−21−dAd=−2d.
The derivative of the remainder is an ordinarily convergent series:
Rd′(0)=n≥1∑[−logbn+nAd]=j=0∑d−1logΓ(1+j)+γAd.
To obtain the last equality, take the logarithm of the
gamma-function product:
logΓ(1+j)=−γj+n≥1∑[nj−log(1+nj)].
Only finitely many j are summed. Differentiating the continuation formula
then cancels the γAd terms and gives
Zd′(0)=−2dlog(d!)+dζ′(0)+j=0∑d−1log(j!)=−2dlog(2πd!)+j=0∑d−1log(j!).
The Riemann values used here, including
ζ′(0)=−21log(2π), are recorded in
DLMF §25.6.
Regularized product and scale
Define the zeta determinant by
detζPd=exp[−Zd′(0)]. The result is
detζPd=∏j=0d−1j!(2πd!)d/2.
This regularized product has the following simple benchmarks:
The first row follows directly from Z1=ζ. If the physical
eigenvalue is μPd(n), choose a fixed comparison unit E∗ and set
r=μ/E∗>0. Then ZrPd(s)=r−sZd(s) and
logdetζ(rPd)=logdetζPd−2dlogr.
Thus the normalization convention remains explicit. Zd(0) is an
analytically continued spectral quantity, not a count of spatial
directions. A gravitational or vacuum-energy use additionally needs a
specified physical operator, action, regularization prescription, and
scale matching.
spectral_analysis_checks.py
checks convergent sums, remainder bounds, and independent centered-polynomial
continuations for d=2,3, including numerical derivatives at zero and scale
changes. The general formulas follow from the analytic derivation above.