Multiple zeta values for classical special functions

We compute multiple zeta values (MZVs) built from the zeros of various entire functions, usually special functions with physical relevance. In the usual case, MZVs and their linear combinations are evaluated using a morphism between symmetric functions and multiple zeta values. We show that this technique can be extended to the zeros of any entire function, and as an illustration, we explicitly compute some MZVs based on the zeros of Bessel, Airy, and Kummer hypergeometric functions. We highlight several approaches to the theory of MZVs, such as exploiting the orthogonality of various polynomials and fully utilizing the Weierstrass representation of an entire function. On the way, an identity for Bernoulli numbers by Gessel and Viennot is revisited and generalized to Bessel–Bernoulli polynomials, and the classical Euler identity between the Bernoulli numbers and Riemann zeta function at even argument is extended to this same class.