Function Repository Resource:

 MultipleZetaStar

Source Notebook

Compute the multiple zeta–star special function

Contributed by: Jayanta Phadikar

ResourceFunction["MultipleZetaStar"][{r1,r2,…,rk}]

gives the multiple zeta star function ζ*(r1,r2,…,rk) defined by the nested sum .

Details

Series definition: .
ResourceFunction["MultipleZetaStar"] is also known as "Multiple Zeta Star". It is similar to MultipleZeta function. The difference is that MultipleZeta uses strictly ordered indices in the series definition, whereas MultipleZetaStar uses weakly ordered indices.
ResourceFunction["MultipleZetaStar"] gives nested sums of reciprocal powers, extending Zeta to multiple weakly ordered indices.
ResourceFunction["MultipleZetaStar"] appears as periods in number theory, as coefficients in ε-expansions of Feynman integrals, and in knot invariants.
MultipleZetaStar is a mathematical function, suitable for both symbolic and numerical manipulation.
The parameters can be arbitrary complex numbers. The defining series converges whenever the real part of each initial partial sum of the arguments exceeds its length, i.e. Re[r1+⋯+rm]>m for m=1,…,k.
In ζ*(r1,r2,…,rk), r1, r2,… are the weights, and k is the depth.
ResourceFunction["MultipleZetaStar"][{r}] is same as MultipleZeta[{r}] and Zeta[r].

Examples

Basic Examples (2) 

Exact values are generated automatically for small weights:

In[1]:=
ResourceFunction["MultipleZetaStar"][{2, 3}]
Out[1]=

Evaluate numerically:

In[2]:=
N[ResourceFunction["MultipleZetaStar"][{12, 2, 4}]]
Out[2]=

Scope (14) 

Numerical Evaluation (3) 

Evaluate with large weight:

In[3]:=
N[ResourceFunction["MultipleZetaStar"][{6, 12, 18}]]
Out[3]=

Evaluate to arbitrary precision:

In[4]:=
N[ResourceFunction["MultipleZetaStar"][{5, 3}], 50]
Out[4]=

Evaluate with rational and complex weights:

In[5]:=
N[ResourceFunction["MultipleZetaStar"][{6, 1/3, 2/3}]]
Out[5]=
In[6]:=
N[ResourceFunction["MultipleZetaStar"][{6, 1/3, 2 + I}]]
Out[6]=

Specific Values (5) 

For total weight ≤ 7, all possible convergent MultipleZetaStar will autoevaluate in terms of Zeta:

In[7]:=
Table[ResourceFunction["MultipleZetaStar"][{n, 1}], {n, 2, 6}]
Out[7]=

For total weight > 7, certain MultipleZetaStar usages may not be reducible in terms of Zeta:

In[8]:=
ResourceFunction["MultipleZetaStar"][{5, 3}]
Out[8]=

Reduction of a MultipleZetaStar of weight 10:

In[9]:=
ResourceFunction["MultipleZetaStar"][{9, 1}]
Out[9]=

MultipleZetaStar starting with 1 in the parameter list is divergent:

In[10]:=
ResourceFunction["MultipleZetaStar"][{1, 6, 2}]
Out[10]=

Failure of the convergence criterion can make a MultipleZetaStar divergent:

In[11]:=
ResourceFunction["MultipleZetaStar"][{1/2, 2}]
Out[11]=

Visualization (5) 

Log-scale discrete plot of the diagonal depth-2 values MultipleZetaStar[n, n] across n, highlighting rapid decay with increasing weight:

In[12]:=
DiscretePlot[ResourceFunction["MultipleZetaStar"][{n, n}], {n, 2, 12},
  ScalingFunctions -> {None, "Log"}]
Out[12]=

Side-by-side log-scale traces of MultipleZetaStar[{n,2}] and MultipleZetaStar[{2,n}] that emphasize order asymmetry under index swap:

In[13]:=
DiscretePlot[
 Evaluate@{ResourceFunction["MultipleZetaStar"][{n, 2}], ResourceFunction["MultipleZetaStar"][{2, n}]}, {n, 2, 16}, ScalingFunctions -> {None, "Log"}, PlotStyle -> {
Directive[Thick, 
RGBColor[0.1, 0.47, 0.71]], 
Directive[Thick, 
RGBColor[0.89, 0.1, 0.11]]}]
Out[13]=

Discrete 3-D bars of Log10[MultipleZetaStar[{s1,s2}]] over admissible indices, showing how magnitude varies across the grid:

In[14]:=
DiscretePlot3D[
 Abs@Log10@ResourceFunction["MultipleZetaStar"][{s1, s2}], {s1, 2, 10}, {s2, 1, 10}, Sequence[
 ExtentSize -> Full, ColorFunction -> "NeonColors", AxesLabel -> Automatic]]
Out[14]=

Log-scale discrete plot of the diagonal depth-3 values MultipleZetaStar[{n,n,n}] versus n, showing depth-driven suppression with weight:

In[15]:=
DiscretePlot[
 ResourceFunction["MultipleZetaStar"][{n, n, n}], {n, 2, 12}, ScalingFunctions -> {None, "Log"}]
Out[15]=

Dependence of MultipleZetaStar[{4,4,n}] on the tail index n at fixed head (4, 4):

In[16]:=
DiscretePlot[ResourceFunction["MultipleZetaStar"][{4, 4, n}], {n, 10}]
Out[16]=

Function Representation (1) 

Obtain MultipleZetaStar from its series representation:

In[17]:=
Sum[1/(n1^2 n2^4 n3^1), {n1, 1, Infinity}, {n2, 1, n1}, {n3, 1, n2}]
Out[17]=
In[18]:=
% == ResourceFunction["MultipleZetaStar"][{2, 4, 1}]
Out[18]=

Properties and Relations (2) 

At depth-1, MultipleZetaStar, MultipleZeta and Zeta are the same:

In[19]:=
ResourceFunction["MultipleZetaStar"][{10}] === MultipleZeta[{10}] === Zeta[10]
Out[19]=

A convergent MultipleZetaStar can be written as a linear combination of MultipleZeta's:

In[20]:=
Inactive[
ResourceFunction["MultipleZetaStar"]][{2, 6, 8}] == Sum[1/(n1^2 n2^6 n3^8), {n1, 1, Infinity}, {n2, 1, n1}, {n3, 1, n2}]
Out[20]=

Verify the identity numerically:

In[21]:=
N[Activate[%]]
Out[21]=

Possible Issues (1) 

The argument of MultipleZetaStar must be a list:

In[22]:=
ResourceFunction["MultipleZetaStar"][2]
Out[22]=

Publisher

Jayanta Kumar Phadikar

Requirements

Wolfram Language 13.0 (December 2021) or above

Version History

  • 1.0.0 – 22 August 2026

License Information