---
title: Multiple Zeta Star Functions
url: https://www.emergentmind.com/topics/multiple-zeta-star-functions-mzsfs
type: topic
---

# Multiple Zeta Star Functions

Multiple zeta star functions (MZSFs) are star-variants of Euler–Zagier multiple zeta functions in which weak inequalities replace strict ones. In depth \(r\), they are given by
\[
\zeta_r^\star(s_1,\dots,s_r)=\sum_{0<n_1\le n_2\le \cdots \le n_r}\frac{1}{n_1^{s_1}\cdots n_r^{s_r}},
\]
and the corresponding multiple zeta-star values arise by specializing to positive integer indices. Recent work places MZSFs at the intersection of analytic continuation, iterated integrals, generating functions, Schur-theoretic determinant identities, Bell- and Stirling-polynomial expansions, positive-characteristic zeta theory, and Diophantine approximation. The subject is now broad enough that the same object is routinely studied as a multivariable meromorphic function, a period-like integral, a quasi-symmetric specialization, and a source of explicit special-value formulae [2508.12577].

## 1. Definitions, domains, and basic variants

For an index \(\mathbf{k}=(k_1,\ldots,k_n)\) of positive integers with \(k_1\ge 2\), the multiple zeta-star value is
\[
\zeta^\star(\mathbf{k})=\sum_{m_1\ge \cdots \ge m_n\ge 1}\frac{1}{m_1^{k_1}\cdots m_n^{k_n}}.
\]
The weak inequalities distinguish the star theory from the ordinary multiple zeta function, where the summation is strict. In the function-theoretic setting, \(\zeta_r^\star(s_1,\dots,s_r)\) converges in a domain \(D_r\), admits meromorphic continuation, and has known singularities. For identical arguments, the standard shorthand is
\[
\zeta_r^\star(s):=\zeta_r^\star(s,\ldots,s),
\]
with the convention \(\zeta_0^\star(\varnothing)=1\) [2508.12577, 2603.27266].

At non-positive integers, two limiting procedures are used. For \(l=(l_1,\dots,l_r)\in \mathbb Z_{\ge 0}^r\),
\[
\zeta_r^{\star,\mathrm{reg}}(-l):=\lim_{s_1\to -l_1}\cdots\lim_{s_r\to -l_r}\zeta_r^\star(s_1,\dots,s_r),
\]
while
\[
\zeta_r^{\star,\mathrm{rev}}(-l):=\lim_{s_r\to -l_r}\cdots\lim_{s_1\to -l_1}\zeta_r^\star(s_1,\dots,s_r).
\]
These regular and reverse values play a central role in explicit evaluation theorems at non-positive integers [2508.12577].

The same star principle extends to a substantial family of variants. Finite and symmetric multiple zeta-star values support Ohno-type generating functions [1905.04875]. Multiple Hurwitz zeta-star functions admit analogous sum formulae [1704.05636]. Multiple star \(t\)-values
\[
t^\star(\alpha_1,\ldots,\alpha_k)=\sum_{1\le j_1\le \cdots \le j_k}\prod_{i=1}^k (2j_i-1)^{-\alpha_i}
\]
form the odd-denominator analogue of MZSFs [1609.01362]. In positive characteristic, one has \(\infty\)-adic and \(v\)-adic multiple zeta star functions defined by weakly decreasing degree indices rather than integer summation variables [2201.12953].

## 2. Relations with ordinary multiple zeta functions and Schur-type generalizations

A foundational fact is that MZSFs and ordinary multiple zeta functions (MZFs) are mutually convertible by summing over all ways of merging adjacent arguments:
\[
\zeta^\star(\mathbf{s})=\sum_{\mathbf{t}\preceq \mathbf{s}}\zeta(\mathbf{t}),\qquad
\zeta(\mathbf{s})=\sum_{\mathbf{t}\preceq \mathbf{s}}(-1)^{n-\ell(\mathbf{t})}\zeta^\star(\mathbf{t}),
\]
where \(\ell(\mathbf{t})\) is the length of \(\mathbf{t}\). This already shows that the star operation is not merely a notational variant: it is tied to a specific composition-refinement combinatorics [1704.08511].

That relation becomes especially sharp at non-positive integers. If \(l=(l_1,\dots,l_r)\in\mathbb Z_{\ge0}^r\) with \(l_1\ge1\), then
\[
\zeta_r^{\mathrm{reg}}(-l)=(-1)^{r+|l|}\zeta_r^{\star,\mathrm{reg}}(-l),\qquad
\zeta_r^{\mathrm{rev}}(-l)=(-1)^{r+|l|}\zeta_r^{\star,\mathrm{rev}}(-l),
\]
where \(|l|=l_1+\cdots+l_r\). A common misconception is that star and non-star values are always substantially different; the non-positive integer theory shows that, under this hypothesis, the difference collapses to a sign, even though other parts of the theory retain distinct combinatorial content [2508.12577].

Schur multiple zeta functions provide a combinatorial interpolation between MZFs and MZSFs. For a partition \(\lambda\), the Schur multiple zeta function is a sum over semi-standard Young tableaux. The one-column case \(\lambda=(1^n)\) recovers the ordinary multiple zeta function, while the one-row case \(\lambda=(n)\) recovers the MZSF. Under content parametrization, hook-type Schur multiple zeta functions decompose explicitly into products of ordinary and star zeta functions:
\[
\zeta_{(p+1,1^q)}(\mathbf{s})
=
\sum_{j=0}^q (-1)^j
\zeta^\star(z_{-j},\dots,z_{-1},z_0,\ldots,z_p)\,
\zeta(z_{-j+1},\dots,z_{-q}),
\]
and arbitrary shapes are then reconstructed via the Giambelli determinant. This places MZSFs inside a larger Schur-theoretic framework governed by Frobenius coordinates, determinant identities, and root-system analogues [2301.05801, 1704.08511].

## 3. Explicit special values: Stirling polynomials, Bell polynomials, and singularities

A major recent development is the explicit evaluation of MZSFs at non-positive integers by Stirling polynomials. Let \(l_1,\dots,l_r\in\mathbb Z_{\ge0}\), \(L_j=l_1+\cdots+l_j\), \(K_j=k_1+\cdots+k_j\), and \(L_0=K_0=0\). If \(S(n,m,X)\) denotes the Stirling polynomial of the second kind, defined by
\[
(X+Y)^n=\sum_{m=0}^n S(n,m,Y)(X)_m,
\]
then the reverse value of the depth-\(r\) MZSF is given explicitly by
\[
\zeta_r^{\star,\mathrm{rev}}(-l)
=
\sum_{0\le k_i\le l_i}
\left[
\prod_{j=1}^r
(-1)^{l_j-k_j}
S(l_j,k_j,L_{j-1}+j)
\frac{(K_j+j-1)!}{(K_{j-1}+j-1)!}
\right]
\zeta_r^{\star,\mathrm{rev}}(\mathbf{0}).
\]
The same paper gives recurrences and a Gregory-coefficient expansion for reverse values, extending a line of work associated with Matsusaka, Murahara, and Onozuka. The generalized Gregory coefficients \(G_{m,n}\) are defined by
\[
\mathcal G(x,y)=\sum_{m,n\ge0}G_{m,n}x^m y^n
=
\frac{y\log^2(1+x)-x\log^2(1+y)}{\log(1+x)-\log(1+y)},
\]
and reverse values are expressed as finite sums built from these coefficients and Stirling-polynomial weights [2508.12577].

For identical arguments, harmonic-product recurrences lead to Bell-polynomial formulae. For \(r\ge1\) and \(s>1\),
\[
r\,\zeta_r^\star(s)=\sum_{j=1}^r \zeta_{r-j}^\star(s)\,\zeta(js),
\]
and more explicitly
\[
\zeta_r^\star(s)=
\sum_{\substack{c_1,\ldots,c_r\ge0\\ c_1+2c_2+\cdots+rc_r=r}}
\frac{1}{1^{c_1}c_1!\,2^{c_2}c_2!\cdots r^{c_r}c_r!}
\zeta(s)^{c_1}\zeta(2s)^{c_2}\cdots \zeta(rs)^{c_r}.
\]
Equivalently,
\[
\zeta_r^\star(s)=\frac1{r!}\mathbf{Y}_r\bigl(0!\zeta(s),1!\zeta(2s),\ldots,(r-1)!\zeta(rs)\bigr),
\]
where \(\mathbf Y_r\) is the complete Bell polynomial. This makes the singularity structure accessible through the poles of \(\zeta(ks)\). The same work notes that the star function has the same possible singularities as the non-star function at the points \(s=1/k\), and that the Bell-polynomial formula gives a direct route to residues and pole orders. It also records special values such as
\[
\zeta_r^\star(2)=\frac{(-1)^{r+1}(2^{2r}-2)B_{2r}\pi^{2r}}{(2r)!},
\qquad
\zeta_r^\star(0)=-\frac{1}{r\,2^{2r-1}}\binom{2r-2}{r-1}.
\]
These formulae complement the non-positive integer theory by showing that the identical-argument specialization has its own closed algebraic model [2603.27266].

## 4. Generating functions, block structures, and explicit families

Generating functions are one of the most effective ways to organize MZSF families. For repeated blocks of length \(1\),
\[
\sum_{j\ge0}\zeta^\star(\{a\}^j)x^j
=
\prod_{n\ge1}\left(1-\frac{x}{n^a}\right)^{-1}.
\]
For more complicated patterns, generating series can often be converted into hypergeometric expressions and then attacked by creative telescoping [2404.16199].

One influential line of work develops generating functions for arbitrary patterns of blocks of \(2\)s and intervening indices, expressing MZSFs in terms of multiple alternating Euler sums and reducing the lengths of blocks of twos. In particular, explicit series are obtained for
\[
\zeta^\star(\{2\}^{a_0},c_1,\{2\}^{a_1},\ldots,c_d,\{2\}^{a_d}),
\]
with the coefficients controlled by multiple sharp sums \(S_{k,m}^\sharp(\{1\}^c)\). This yields a systematic block-reduction principle for star values that had previously been available only in special configurations [1806.09142].

Specific index families have been worked out in closed form. For the 3–2–1 sector, generating-function methods and Bell-polynomial expansions yield explicit formulae for \(\zeta^\star(\{3,1\}^d)\), \(\zeta^\star(\{3,1\}^d,2)\), and the generalized families
\[
\zeta^\star(\{\{2\}^m,3,\{2\}^m,1\}^d),\qquad
\zeta^\star(\{\{2\}^m,3,\{2\}^m,1\}^d,\{2\}^{m+1}),
\]
including product formulae for the corresponding generating functions in terms of trigonometric factors and coefficient formulae in Bernoulli numbers [1806.10510]. A different strand, based on creative telescoping, establishes new evaluations for block patterns such as \(\zeta^\star(\{1,3\}^n,1,2)\), reducing them to polynomials in Riemann zeta values after first passing through interpolated half-values:
\[
\zeta^\star(\{1,3\}^n,1,2)=2^{2n+1}\zeta^{1/2}(\{2\}^{2n},3).
\]
The paper also treats \(\zeta^\star(2,\{1,3\}^n,1,2)\), \(\zeta^\star(1,2,\{1,3\}^n)\), and \(\zeta^\star(2,3,\{1,3\}^n)\) by the same pipeline [2404.16199].

Finite and symmetric star values support their own generating-function theory. For an index \(\mathbf{k}\), Kaneko’s conjecture for \(\mathcal O_F(2,1,2)\) was confirmed and generalized to arbitrary depth-three indices. If
\[
\mathcal O_F(\mathbf{k})
=
\sum_{\mathbf e\ge0}\zeta_F^\star(\mathbf{k}\oplus \mathbf e)X^{|\mathbf{k}|+|\mathbf e|}
+
\sum_{\mathbf e\ge0}(-1)^{|\mathbf e|}\zeta_F^\star((\mathbf{k}\oplus \mathbf e)^\vee)X^{|\mathbf{k}|+|\mathbf e|},
\]
then for \(k_2=1\),
\[
\mathcal O_F(k_1,1,k_3)=-F_{k_1,1}(X)F_{k_3,1}(X),
\]
and for \(k_2\ge2\),
\[
\mathcal O_F(k_1,k_2,k_3)
=
-\sum_{\substack{i,j\ge2\\ i+j=k_2+1}}F_{k_1,i}(X)F_{k_3,j}(X).
\]
These formulas encode duality through the Hoffman dual and show that, in the finite and symmetric settings, depth-three Ohno-type sums reduce to quadratic forms in depth-one series [1905.04875].

## 5. Integral, cyclic, and combinatorial frameworks

An integral theory for MZSFs parallel to the classical iterated-integral theory of ordinary multiple zeta values was established by Yamamoto. For an index \(\mathbf{k}=(k_1,\ldots,k_n)\), let \(K=|\mathbf{k}|\) and
\[
J(\mathbf{k})=\{0,k_1,k_1+k_2,\ldots,k_1+\cdots+k_{n-1}\}.
\]
With \(\omega_0(t)=dt/t\) and \(\omega_1(t)=dt/(1-t)\), one has
\[
\zeta^\star(\mathbf{k})
=
\int_{\Delta(\mathbf{k})}\omega_{\delta(1)}(t_1)\cdots \omega_{\delta(K)}(t_K),
\]
where \(\Delta(\mathbf{k})\subset[0,1]^K\) is determined by alternating inequalities \(t_j<t_{j+1}\) or \(t_j>t_{j+1}\) according to membership in \(J(\mathbf{k})\). For \(\mathbf{k}=(2,1)\),
\[
\zeta^\star(2,1)=\int_{0<t_1<t_2>t_3<1}\frac{dt_1}{1-t_1}\frac{dt_2}{t_2}\frac{dt_3}{1-t_3}.
\]
This integral formalism is subsumed by the more general theory of admissible 2-labeled posets. If \(X\) is such a poset, then
\[
I(X)=\int_{A(X)}\prod_{x\in X}\omega_{\delta_X(x)}(t_x),
\]
with multiplicativity, decomposition under refinement of incomparable elements, and duality \(I(X^\ast)=I(X)\). The same framework covers MZVs, MZSFs, Arakawa–Kaneko zeta values, Mordell–Tornheim zeta values, and root-system zeta functions of type \(A\) [1405.6499].

Cyclic relations extend the algebraic structure beyond integral representations. The cyclic relation of Hirose, Murakami, and Murahara was generalized to complex variables, and the specialization to positive integers recovers the cyclic sum formula for multiple zeta-star values:
\[
\sum_{i=1}^d\sum_{m=1}^{k_i-1}
\zeta^\star(k_i-m,k_{i+1},\dots,k_d,k_1,\dots,k_{i-1},m+1)
=
\zeta(k_1+\cdots+k_d+1).
\]
The complex-variable version places star and non-star cyclic identities inside a unified analytic framework that also includes derivation relations [2003.08068].

Several explicitly evaluable star families were obtained by integral and iterated-integral methods. Examples include polynomial expressions for
\[
\zeta^\star(\bar1,\{1\}_m,\bar1),\qquad
\zeta^\star(2,\{1\}_m,\bar1),
\]
in terms of classical zeta values, polylogarithms, and powers of \(\ln 2\) [1702.03868], as well as explicit relations between MZSVs, Kaneko–Yamamoto type MZVs, alternating MZVs, and the Apéry-type multiple zeta \(B\)-star values
\[
\zeta_B^\star(k_1,\ldots,k_r)=
\sum_{n=1}^\infty \frac{S_n^\star(k_2,\ldots,k_r)}{n^{k_1}4^n\binom{2n}{n}},
\]
which are shown in many cases to be \(\mathbb Q\)-linear combinations of alternating MZVs [2012.02588].

## 6. Arithmetic extensions, positive characteristic, and limiting-set phenomena

In positive characteristic, MZSFs admit \(\infty\)-adic and \(v\)-adic analogues. With \(A=\mathbb F_q[\theta]\), \(k=\mathbb F_q(\theta)\), and \(s_j\in X_\infty\), the \(\infty\)-adic star function is
\[
\zeta_\infty^\ast(s_1,\ldots,s_r)
=
\sum_{i_1\ge \cdots \ge i_r\ge0} S_{i_1}(s_1)\cdots S_{i_r}(s_r),
\]
while the \(v\)-adic version is
\[
\zeta_v^\ast(s_1,\ldots,s_r)
=
\sum_{i_1\ge \cdots \ge i_r\ge0}\tilde S_{i_1}(s_1)\cdots \tilde S_{i_r}(s_r).
\]
Both are rigid analytic; for \(r=1\) they recover the single-variable zeta functions of Carlitz–Goss and Goss. A central relation is the orthogonal property
\[
\sum_{l=0}^r(-1)^l\,
\zeta_\infty(s_r,\ldots,s_{r-l+1})\,
\zeta_\infty^\ast(s_1,\ldots,s_{r-l})=0,
\]
with a parallel identity in the \(v\)-adic setting. The same theory gives integral expressions for negative integer values, Kummer-type congruences such as
\[
\zeta_v^\ast(m_1,\ldots,m_r)\equiv \zeta_v^\ast(l_1,\ldots,l_r)\pmod{\mathfrak m_v^e},
\]
and recursive links between \(v\)-adic and \(\infty\)-adic star values [2201.12953].

A separate arithmetic direction studies the set of multiple zeta-star values itself as a subset of the real line. The set \(\mathbb Z^\ast\) of all multiple zeta-star values is a countable dense subset of \((1,+\infty)\). A Dirichlet-type criterion characterizes when \(a>1\) belongs to \(\mathbb Z^\ast\), and a Khinchin-type zero-one law states that for
\[
F_f=\left\{a\in(1,\infty):\ |a-\zeta^\ast(k_1,\ldots,k_r)|<\sum_{n_1\ge\cdots\ge n_{r+1}\ge2}\frac{1}{n_1^{k_1}\cdots n_r^{k_r}n_{r+1}f(r)}
\ \text{infinitely often}\right\},
\]
one has \(m(F_f)=0\) if \(\sum_r \frac1{2f(r)}\) converges and \(m((1,\infty)\setminus F_f)=0\) if it diverges [2503.23286].

Finite multiple star harmonic sums admit an analytic interpolation built from divided differences. For fixed \(n\),
\[
H_n^\star(k_1,\ldots,k_r)=\sum_{n\ge m_1\ge \cdots \ge m_r\ge1}\frac1{m_1^{k_1}\cdots m_r^{k_r}},
\]
and the interpolated function \(H_s^\star(k_1,\ldots,k_r)\) has the integral representation
\[
H_s^\star(k_1,\ldots,k_r)
=
\int_{(0,1)^{k_1+\cdots+k_r}}
[1,X_1,\ldots,X_r]\, t^{s+r-1}\,dx_1\cdots dx_{k_1+\cdots+k_r},
\]
where \([1,X_1,\ldots,X_r]\) is a divided difference and the Hermite–Genocchi formula is used to prove analyticity. For real \(s>1\), the resulting map \(\mathfrak h_s\) is a bijection onto \((1,s]\), giving a finite zeta-star correspondence; for complex \(s\) with \(\Re s>1\), injectivity is conjectured [2606.16627].

Taken together, these developments show that MZSFs are no longer studied only as a star-version of Euler–Zagier sums. They form a nexus joining explicit special-value theory, determinant and tableau formalisms, iterated-integral and generating-function methods, positive-characteristic interpolation, and metric questions about the distribution and coding of zeta-star values themselves.

Source: https://www.emergentmind.com/topics/multiple-zeta-star-functions-mzsfs