Symmetric Multiple Eisenstein Series in Modular Forms

• Symmetric multiple Eisenstein series are holomorphic q-series that generalize classic Eisenstein series and symmetric multiple zeta values.
• Their construction employs shuffle-regularization and Hopf algebra techniques to establish essential algebra homomorphism properties.
• They reveal deep connections with modular forms, period polynomials, and elliptic zeta values, offering a modular framework for multiple zeta relations.

Symmetric multiple Eisenstein series are a class of holomorphic q-series on the modular group, motivated by the theory of symmetric multiple zeta values and possessing deep connections to modular forms, period polynomials, and elliptic zeta values. They generalize both classical Eisenstein series and the symmetric multiple zeta values introduced by Kaneko and Zagier, and are constructed so as to satisfy distinguished shuffle-type relations within a noncommutative algebraic framework. Recent work by Hara, Sakugawa, and Tasaka details their definition, algebraic properties, structural results, and interplay with modular forms and elliptic zeta relations (Hara et al., 20 Jan 2026).

1. Construction and Regularization of Symmetric Multiple Eisenstein Series

Let $\tau$ lie in the complex upper half-plane, and consider the lattice $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$. For strictly positive integers $k_1, \dots, k_d \geq 2$, the (depth-$d$) multiple Eisenstein series is defined as the iterated, conditionally convergent sum

$$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$

with respect to a total order $<$ on $\Lambda$ determined by a positive cone $P$. This generalizes the classical Eisenstein series ($d=1$).

Multiple Eisenstein series $G_{k_1,\ldots,k_d}(\tau)$ do not a priori satisfy the shuffle product relations characteristic of multiple zeta values. To remedy this, a shuffle-regularized version $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$0 is defined, providing an algebra homomorphism from the shuffle algebra on noncommutative generators $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$1. This construction utilizes generating functions $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$2 satisfying (quasi-)shuffle or stuffle relations, and their Taylor expansions give rise to $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$3, which do satisfy shuffle product rules. The regularization is made precise using Hopf algebra techniques, notably Goncharov's coproduct.

Motivated by the notion of symmetric multiple zeta values, the symmetric shuffle-regularized multiple Eisenstein series is defined as

$\Lambda = \mathbb{Z}\tau + \mathbb{Z}$4

Alternatively, this can be expressed as a sum over strictly ordered indices with respect to the total Kontsevich order on $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$5, which extends the lattice order to a cyclic, infinite total order identifying $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$6 with $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$7.

2. Linear Shuffle Relations

Symmetric multiple Eisenstein series satisfy a fundamental class of algebraic identities: linear shuffle relations. For a decomposition $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$8 and the shuffle product $\Lambda = \mathbb{Z}\tau + \mathbb{Z}$9 in the noncommutative algebra $k_1, \dots, k_d \geq 2$0, the following holds: $k_1, \dots, k_d \geq 2$1 This relation is most naturally phrased via generating series in commuting variables $k_1, \dots, k_d \geq 2$2: $k_1, \dots, k_d \geq 2$3 with corresponding identities involving application of shuffle operators $k_1, \dots, k_d \geq 2$4. The only obstruction to $k_1, \dots, k_d \geq 2$5 being an algebra homomorphism arises from these specific “twisted” linear shuffle relations, which mirror the algebraic constraints found for symmetric multiple zeta values.

3. Structure and Dimensions of Spaces Spanned by Symmetric Double Eisenstein Series

Define, for weight $k_1, \dots, k_d \geq 2$6,

$k_1, \dots, k_d \geq 2$7

The structure of $k_1, \dots, k_d \geq 2$8 depends on the parity of $k_1, \dots, k_d \geq 2$9:

• Even weight $d$0:

  $d$1

  with explicit dimension

  $d$2

  The proof relies on shuffle and stuffle arguments, explicit formulas for $d$3, and known decompositions for classical modular forms.

• Odd weight $d$4:

  $d$5

  is a $d$6-basis, so

  $d$7

  The argument uses the linear shuffle relations, an $d$8-representation argument for dimension, and linear independence established via a rank argument on a matrix of binomial coefficients extracted from $d$9-expansions.

4. Fay–Shuffle Correspondence and Elliptic Double Zeta Values

Elliptic double zeta values, as introduced by Matthes, satisfy two classes of linear relations: classical shuffle, and Fay–shuffle stemming from Fay’s identity for elliptic integrals. The space of polynomial Fay–shuffle relations in weight $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$0,

$$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$1

has dimension $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$2. There is an explicit correspondence between these Fay–shuffle spaces and the linear shuffle space $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$3 satisfied by the symmetric Eisenstein series:

• For $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$4 even, the odd period polynomial subspace $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$5.
• For $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$6 odd, there is a $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$7-linear bijection:

  $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$8

achieved by transformations involving generators $$G_{k_1,\ldots,k_d}(\tau) = \lim_{M\to\infty} \lim_{N\to\infty} \sum_{\substack{0<\lambda_1<\cdots<\lambda_d \ \lambda_i\in(-M\tau-N,\dotsc, M\tau+N)} } \lambda_1^{-k_1}\cdots \lambda_d^{-k_d}$$9 of $<$0-type symmetries, with linear conditions checked via explicit substitution.

This correspondence elucidates the parallel between the algebraic relations among symmetric double Eisenstein series and the structure of the elliptic double zeta value relations.

5. Modular Forms as Sums of Symmetric Triple Eisenstein Series

Every modular form for $<$1 with rational coefficients and weight $<$2 can be written as a finite $<$3-linear combination of symmetric triple Eisenstein series: $<$4 Defining $<$5 as the span of all depth-3 symmetric Eisenstein series in weight $<$6, the argument leverages the fact that products $<$7, $<$8, and $<$9 belong to $\Lambda$0 and that these are known to span all modular forms. Thus, symmetric triple Eisenstein series suffice to generate the entire space of modular forms of a given weight.

Taking the limit as $\Lambda$1, the constant terms give rise to linear relations among symmetric multiple zeta values $\Lambda$2, in accordance with the modular relations originally posited by Kaneko–Zagier.

6. Implications and Connections

The introduction and systematic study of symmetric multiple Eisenstein series reveals several significant structures:

• Linear shuffle relations satisfied by these series directly parallel those found in the theory of symmetric multiple zeta values.
• Their expansion space in depth two and three aligns with the structure and dimension formulas of modular forms and the corresponding period polynomial spaces.
• They establish the first concrete modular framework to systematically derive Kaneko–Zagier type relationships for multiple zeta values, particularly the “modular phenomena” conjectured in higher weight and depth (Hara et al., 20 Jan 2026).
• The correspondence with Fay–shuffle relations for elliptic zeta values indicates a unity between the algebraic structures of modular forms, period polynomials, and elliptic and classical multiple zeta values.

A plausible implication is that symmetric multiple Eisenstein series act as a bridge between classical modular form theory and the arithmetic of multiple zeta values, providing a unifying perspective for further study of their special values, $\Lambda$3-expansions, and motivic or cohomological interpretations.

7. Summary Table

Structure	Object/Invariants	Dimension/Decomposition
Symmetric double Eisenstein	$\Lambda$4	Even $\Lambda$5: $\Lambda$6; Odd $\Lambda$7: $\Lambda$8
Fay–shuffle polynomials	$\Lambda$9	$P$0
Modular forms in triple series	$P$1	All cusp forms admit decomposition by symmetric triple series

Further technical details, Fourier expansions, and coproduct computations are elaborated in Hara–Sakugawa–Tasaka’s paper (Hara et al., 20 Jan 2026).