Mixed Parity Variants of Apéry-Type Binomial Series and Level Four Colored Multiple Zeta Values

In this paper, we study an Apéry-type series involving the central binomial coefficients ∑n1>…>nd>014n12n1n1n1−s1…nd−sd and its variations where the summation indices may have mixed parities and some or all “>” are replaced by “≥”, as long as the series are defined. These sums have naturally appeared in the calculation of massive Feynman integrals by the work of Jegerlehner, Kalmykov, and Veretin. We show that all these sums can be expressed as Q-linear combinations of the real and/or imaginary parts of the colored multiple zeta values at level four, i.e., special values of multiple polylogarithms at fourth roots of unity. For example, our main theorem shows that when n1−s1 is replaced by (2n1)−s1 and other nj−sj’s are replaced by either (2nj)−sj or (2nj+1)−sj, then all the colored multiple zeta values can be chosen to have the same weight s1+…+sd, but the weights of these values are only bounded by s1+…+sd for general variant Apéry-type series of mixed parities. We also show that the corresponding series where 2n1n1/4n1 is replaced by 2n1n12/16n1 can be expressed in a similar way except for a possible extra factor of 1/π, with the weight of the colored multiple zeta values similarly bounded.