  The heat kernel (or propagator) of the quantum harmonic oscillator (qHO) is
given by the Mehler formula, and the partition function is obtained by taking
its trace. In general, the spectral zeta function of the given system is
obtained by the Mellin transform of its partition function. In the case of
non-commutative harmonic oscillators (NCHO), however, the heat kernel and
partition functions are still unknown, although meromorphic continuation of the
corresponding spectral zeta function and special values at positive integer
points have been studied. On the other hand, explicit formulas for the heat
kernel and partition function have been obtained for the quantum Rabi model
(QRM), which is the simplest and most fundamental model for light and matter
interaction in addition to having the NCHO as a covering model. In this paper,
we propose a notion of the quasi-partition function for a quantum interaction
model if the corresponding spectral zeta function can be meromorphically
continued to the whole complex plane. The quasi-partition function for qHO and
QRM actually gives the partition function. Assuming that this holds for the
NCHO (currently a conjecture), we can find various interesting properties for
the spectrum of the NCHO. Moreover, although we cannot expect any functional
equation of the spectral zeta function for the quantum interaction models, we
try to seek if there is some relation between the special values at positive
and negative points. Attempting to seek this, we encounter certain divergent
series expressing formally the Hurwitz zeta function by calculating integrals
of the partition functions. We then give two interpretations of these divergent
series by the Borel summation and $p$-adically convergent series defined by the
$p$-adic Hurwitz zeta function.
