The author was partially supported by NSF grant DMS-0140172.
<ph f="cmbx">On Wronskians of weight one Eisenstein series</ph>

### Lev A. Borisov

Department of Mathematics, University of Wisconsin, 417 Van Vleck Hall, 480 Lincoln Dr, Madison, WI 53706, USA. email:borisov@math.wisc.edu.
• Abstract. We describe the span of Hecke eigenforms of weight four with nonzero central value of $L$  -function in terms of Wronskians of certain weight one Eisenstein series.

1 Introduction

For any positive integer $l$  we consider the congruence subgroup ${\Gamma }_{1}\left(l\right)\subseteq S{l}_{2}\left(\mathbb{Z}\right)$  . The space of cusp forms for ${\Gamma }_{1}\left(l\right)$  of a given weight $k$  splits according to the eigenvalues of Hecke operators. We say that a Hecke eigenform has analytic rank zero, if the central value of the corresponding $L$  -function is nonzero.
It has been shown in [BG1that the span of Hecke eigenforms of weight two coincides with the span of the cuspidal parts of products of certain weight $1$  Eisenstein series for the group ${\Gamma }_{1}\left(l\right)$  . These series are the logarithmic derivatives in the $z$  direction of the standard $\theta$  -function, evaluated at $\frac{a}{l}$  for $a=1,\dots ,l-1$  . It is convenient to look at the Fricke involutions of these Eisenstein series. These are linear combinations of the original series and are given by ${s}_{a}\left(q\right)=\left(\frac{1}{2}-\left\{\frac{a}{l}\right\}\right)+{\sum }_{n>0}{q}^{n}{\sum }_{d|n}\left({\delta }_{d}^{amodl}-{\delta }_{d}^{-amodl}\right),$  where $q=exp\left(2\pi i\tau \right)$  and $\delta$  is a version of Kronecker symbol. In this paper we look at the Wronskians $W\left({s}_{a},{s}_{b}\right)$  defined as usual by $W\left({s}_{a},{s}_{b}\right)=\left(\frac{d}{d\tau }{s}_{a}\right){s}_{b}-\left(\frac{d}{d\tau }{s}_{b}\right){s}_{a}$  . It is easy to see that $W\left({s}_{a},{s}_{b}\right)$  is always a cusp form of weight $4$  , and the main result of this paper relates the span of such forms with the span of Hecke eigenforms of analytic rank zero.
Theorem  6.5 . For arbitrary $l>1$  the span of Hecke eigenforms of weight four and analytic rank zero is equal to the span of the Wronskians $W\left({s}_{a}\left(\tau \right),{s}_{b}\left(\tau \right)\right)$  for all $a,b\in \mathbb{Z}/l\mathbb{Z}$  .
Before we explain the idea of the proof of this paper, we remark that it should be possible to prove Theorem  6.5 using Rankin-Selberg method, by combining the formulas [Z,4.3,equation(4)and [Sc,Theorem4.6.3. However, we chose to use the technique of [BG1and [BG2that emphasizes the map from modular symbols to modular forms.
The space ${M}_{4}\left(l\right)$  of modular symbols of weight four can be thought of as a combinatorial counterpart to the space of modular forms. It is a vector space of roughly twice the dimension, and it contains subspaces ${S}_{4}\left(l{\right)}_{+}$  and ${S}_{4}\left(l{\right)}_{-}$  which are naturally dual to the space ${\mathcal{S}}_{4}\left(l\right)$  of cusp forms of weight four.
Moreover, the action of Hecke operators on the space of modular symbols is given explicitly, see [M1. Ignoring minor complications due to old forms, the span of Hecke eigenforms of weight four and analytic rank zero can be seen as the image of the endomorphism $\rho :{\mathcal{S}}_{4}\left(l\right)\to {\mathcal{S}}_{4}\left(l\right)$  given by $\rho \left(f\right)={\sum }_{n>0}L\left({T}_{n}f,2\right){q}^{n}$  where ${T}_{n}$  denote the Hecke operator. We observe that $L\left(f,2\right){q}^{n}$  is the result of the pairing $〈f,xy\left(0,1{\right)}_{-}〉$  of $f$  a certain element of ${S}_{4}\left(l{\right)}_{-}$  to calculate $\rho$  in terms of modular symbols as a composition of maps ${\mathcal{S}}_{4}\left(l\right)\stackrel{Int}{\to }\left({S}_{4}\left(l{\right)}_{-}{\right)}^{*}\stackrel{PD}{\to }{S}_{4}\left(l{\right)}_{+}\stackrel{\mu }{\to }{\mathcal{S}}_{4}\left(l\right)$  where $Int$  is induced by the integration pairing of ${\mathcal{S}}_{4}\left(l\right)$  and ${S}_{4}\left(l{\right)}_{-}$  , the $PD$  is the Poincaré duality map which we define in Section  3 , and $\mu$  is the Wronskian map, defined in Section  4 .
The map $PD$  is a weight four analog of the intersection pairing on weight $2$  symbols considered in [BG1. It is shown to be nondegenerate in Section  3 as a consequence of a modular symbol formula for Petersson inner product.
The map $\mu$  is the main novelty of this paper. It is a map from the space of modular symbols to the space of modular forms, which in particular maps $xy\left(a,b\right)$  to the Wronskian $W\left({s}_{a},{s}_{b}\right)$  . Our calculations are purely elementary and rely on properties of the Euclid algorithm and some explicit calculations with modular symbols.
There are several directions in which one can try to extend the results of this paper. For example, one can look at the subspaces in the spaces of modular forms of higher weight that are spanned by Wronskians of Eisenstein series of higher weight. Intuition derived from [BG2and [Zsuggests that these would be related to values of the $L$  -function at $2$  . Consequently, we expect the Wronskians to span the whole space in the higher weight setting.
It is worth mentioning that the product and the Wronskian are the first two cases of Cohen operators (see [Z). One can wonder if the forms of higher weight of analytic rank zero can be described in terms of higher Cohen operators of ${s}_{a}$  . Clearly, for a high enough weight this seems impossible for dimension reasons. On the other hand, one could perhaps apply Cohen operators to the theta function itself, rather than its logarithmic derivatives, similar to the definition of $\mu$  on the noncuspidal symbols of weight four. But this is all but a speculation at this point.
One might hope to use the construction of this paper to give upper bounds on the number of Hecke eigenforms of higher analytic rank. However, analogous statements for weight two, at least so far, has not lead to such results. It can also be argued that there may be some deeper reason behind the results of this paper and [BG1which is yet to be uncovered. From this point of view, it would be tempting to try to see the sums along the runs of Euclid algorithm as a calculation of an Euler characteristics of some complex, whose cohomology is located at top and bottom location only. But at the moment we do not have a suitable candidate for it. Finally, one can wonder whether derivatives of $L$  -function at the central value can be somehow seen in terms of Eisenstein series and Cohen operators.
Notations. We denote by $\mathcal{ℋ}$  the upper half-plane and denote by $\tau$  , $\Im \left(\tau \right)>0$  the complex coordinate on it. We use the notation $q=exp\left(2\pi i\tau \right)$  when writing Fourier expansion of modular forms. Throughout the paper $l$  denotes the level, and it is generally fixed, except for the proof of Theorem  6.5 that requires induction on the level. We use a slightly modified Kronecker $\delta$  notation ${\delta }_{u}^{vmodw}$  which gives $1$  when $u=vmodw$  and $0$  otherwise.
Acknowledgments. This paper grew out of a search of a (weight two) skew-symmetric analog of [BG1which the author talked about on and off for a few years with Paul Gunnells. The author also thanks Loïc Merel for helpful remarks regarding the Poincaré duality map.

2 Modular symbols of weight four

Our main reference for modular symbols is the paper [M1by Merel, which in turn builds on the work of Manin and Shokurov. In this section we recall the purely combinatorial description of modular (Manin, in the terminology of [M1) symbols of weight four for the group ${\Gamma }_{1}\left(l\right)$  .
The modular symbols of weight four and level $l$  is a quotient of the vector space with basis ${x}^{2}\left(u,v\right)$  , $xy\left(u,v\right)$  , ${y}^{2}\left(u,v\right)$  , with $\left(u,v\right)\in \left(\mathbb{Z}/l\mathbb{Z}{\right)}^{2}$  , $gcd\left(u,v,l\right)=1$  by the span of the relations
 $\begin{array}{c}\begin{array}{c}{x}^{2}\left(u,v\right)+{y}^{2}\left(v,-u\right),xy\left(u,v\right)-xy\left(v,-u\right),{y}^{2}\left(u,v\right)+{x}^{2}\left(v,-u\right)\\ xy\left(v,-u-v\right)-xy\left(-u-v,u\right)+{y}^{2}\left(-u-v,u\right)+{x}^{2}\left(u,v\right)-xy\left(u,v\right)\end{array}\end{array}$ (2.1)
for all $u,v$  with $gcd\left(u,v,l\right)=1$  .
Remark 2.1. Our set of relations looks somewhat smaller than that of [M1, where the relations are
 $\begin{array}{c}P\left(x,y\right)\left(u,v\right)+P\left(y,-x\right)\left(u,v\right)\end{array}$
 $\begin{array}{c}P\left(x,y\right)\left(u,v\right)+P\left(y-x,-x\right)\left(-u-v,u\right)+P\left(-y,x-y\right)\left(v,-u-v\right)\end{array}$
for an arbitrary degree two homogeneous polynomial $P\left(x,y\right)$  . The ”missing” relations are obtained by cyclic permutations of $\left(u,v,-u-v\right)$  in the last line of  2.1 , so the two definitions of modular symbols are equivalent.
Recall that the subspace ${S}_{4}\left(l\right)\subset {M}_{4}\left(l\right)$  of cuspidal modular symbols is characterized as follows. The cusps of the modular curve ${X}_{1}\left(l\right)=\overline{\mathcal{ℋ}/{\Gamma }_{1}\left(l\right)}$  are in one-to-one correspondence with elements of the set $I=\left\{\left(a,b\right),a\in \mathbb{Z}/l\mathbb{Z},b\in \left(\mathbb{Z}/\left(a,l\right)\mathbb{Z}{\right)}^{*}\right\}/±$  . This correspondence maps $\left(a,b\right)$  to $\frac{{b}^{*}}{a}\in \mathbb{Q}\cup i\infty$  where ${b}^{*}$  is the inverse of $bmod\left(a,l\right)$  . For every element of $I$  there is a map ${M}_{4}\left(l\right)\to \mathbb{C}$  defined by
 $\begin{array}{c}\begin{array}{ccc}{x}^{2}\left(u,v\right)& ↦& {\delta }_{u}^{amodl}{\delta }_{v}^{bmod\left(a,l\right)}+{\delta }_{u}^{-amodl}{\delta }_{v}^{-bmod\left(a,l\right)}\\ {y}^{2}\left(u,v\right)& ↦& -{\delta }_{v}^{amodl}{\delta }_{u}^{-bmod\left(a,l\right)}-{\delta }_{v}^{-amodl}{\delta }_{u}^{bmod\left(a,l\right)}\\ xy\left(u,v\right)& ↦& 0\end{array}\end{array}$ (2.2)
Then the space of cuspidal symbols ${S}_{4}\left(l\right)$  is defined as the intersection of the kernels of all these maps.
We are now ready to formulate the main result of this section.
Proposition 2.2. The space of cuspidal symbols ${S}_{4}\left(l\right)$  is spanned by the modular symbols of the form $xy\left(u,v\right)$  .
• Proof. We can use the first set of equations to solve for ${y}^{2}\left(u,v\right)$  . Then we can think of modular symbols of weight four as being spanned by ${x}^{2}\left(u,v\right)$  and $xy\left(u,v\right)$  , subject to conditions  $\begin{array}{c}\begin{array}{c}{x}^{2}\left(u,v\right)-{x}^{2}\left(-u,-v\right)=0,xy\left(u,v\right)-xy\left(v,-u\right)=0,\\ {x}^{2}\left(u,u+v\right)-{x}^{2}\left(u,v\right)=xy\left(v,-u-v\right)-xy\left(-u-v,u\right)-xy\left(u,v\right).\end{array}\end{array}$ (2.3)
Clearly, $xy\left(u,v\right)$  are cuspidal. On the other hand, if a linear combination of $w={\sum }_{u,v}{\alpha }_{u,v}{x}^{2}\left(u,v\right)$  is cuspidal, then for each $amodl$  and each $bmod\left(a,l\right)$  ${\sum }_{v=bmod\left(a,l\right)}\left({\alpha }_{u,v}+{\alpha }_{-u,-v}\right)=0$  By using relations ${x}^{2}\left(u,v\right)={x}^{2}\left(-u,-v\right)$  we can write $w$  as a linear combination of ${x}^{2}\left(u,v+ku\right)-{x}^{2}\left(u,v+\left(k-1\right)u\right)$  which is then written as a linear combination of $xy\left({u}^{\prime },{v}^{\prime }\right)$  .
Remark 2.3. In addition to the obvious symmetry relations $xy\left(u,v\right)=xy\left(v,-u\right)$  there are still some other linear relations among the symbols $xy\left(u,v\right)$  in ${S}_{4}\left(l\right)$  . In fact, one can show that for $l\ge 5$  the linear relations on $xy\left(u,v\right)$  in ${M}_{4}\left(l\right)$  (or ${S}_{4}\left(l\right)$  ) are spanned by the symmetry relations and ${\sum }_{k=0}^{l-1}\left(xy\left(v+ku,-\left(k+1\right)u-v\right)-xy\left(-\left(k+1\right)u-v,u\right)-xy\left(u,v+ku\right)\right)=0$  for all $u$  and $v$  with $gcd\left(u,v,l\right)=1$  . We leave the proof of this claim to the reader, as it will not be used elsewhere in the paper.
We now recall that ${M}_{4}\left(l\right)$  and ${S}_{4}\left(l\right)$  naturally split according to the eigenvalue of the involution $i$  given by ${x}^{2}\left(u,v\right)↦{x}^{2}\left(-u,v\right),xy\left(u,v\right)↦-xy\left(-u,v\right),{y}^{2}\left(u,v\right)↦{y}^{2}\left(-u,v\right).$  We define the corresponding eigenspaces by ${M}_{4}\left(l{\right)}_{+}$  , ${M}_{4}\left(l{\right)}_{-}$  , ${S}_{4}\left(l{\right)}_{+}$  and ${S}_{4}\left(l{\right)}_{-}$  . There are symmetrization maps ${M}_{4}\left(l\right)\to {M}_{4}\left(l{\right)}_{±}$  given by $t\to \frac{1}{2}\left(t±i\left(t\right)\right)$  , and similarly for ${S}_{4}\left(l\right)\to {S}_{4}\left(l{\right)}_{±}$  . We use a subscript to indicate the symmetrization of a symbol. We can now apply Proposition  2.2 to ${S}_{4}\left(l{\right)}_{±}$  .
Corollary 2.4. The space ${S}_{4}\left(l{\right)}_{±}$  is a linear span of the symbols $xy\left(u,v{\right)}_{±}$  with $\left(u,v\right)\in \left(\mathbb{Z}/l\mathbb{Z}{\right)}^{2}$  and $gcd\left(u,v,l\right)=1$  .
Remark 2.5. It is amusing to observe that for prime $l\ge 3$  the space ${S}_{4}\left(l{\right)}_{+}$  is a linear span of symbols $xy\left(u,v{\right)}_{+}$  with $\left(u,v\right)\in \left(\mathbb{Z}/l\mathbb{Z}{\right)}^{2}-\left(0,0\right)$  with linear relations among these symbols generated by $xy\left(u,v{\right)}_{+}=-xy\left(-u,v{\right)}_{+}=-xy\left(v,u{\right)}_{+}.$  Clearly, these relations hold in ${S}_{4}\left(l{\right)}_{+}$  and, by themselves, they cut its dimension down to at most $\frac{1}{8}\left(l-1\right)\left(l-3\right)$  . On the other hand, by [M1, ${S}_{4}\left(l{\right)}_{+}$  is dual to the space ${\mathcal{S}}_{4}\left(l\right)$  of cusp forms of weight four, which has dimension $\frac{1}{8}\left(l-1\right)\left(l-3\right)$  by the usual Riemann-Roch calculation. This shows that all other relations on $xy\left(u,v{\right)}_{+}$  follow from the above symmetry relations (which can also be checked directly along the lines of Remark  2.3 ). One can thus identify ${S}_{4}\left(l{\right)}_{+}$  with the second exterior power of the vector space of dimension $\left(l-1\right)/2$  which is generated by the symbols ${r}_{a}$  for $amodl$  with ${r}_{-a}=-{r}_{a}$  .

3 Poincare Duality for Modular Symbols

The goal of this section is to explicitly describe a certain map $PD:{M}_{4}\left(l{\right)}^{*}\to {M}_{4}\left(l\right)$  which is a weight four analog of the Poincaré duality for the weight two cuspidal symbols. It is rather easy to show that $PD\left({M}_{4}\left(l\right)\right)\subseteq {S}_{4}\left(l\right)$  . In fact, we will see that $PD\left({M}_{4}\left(l\right)\right)={S}_{4}\left(l\right)$  , which is crucial for the argument of this paper. This is proved by comparison of $PD$  and the expression of the Petersson inner product of cusp forms of weight four in terms of their period integrals.
Definition 3.1. The linear map $PD:{M}_{4}\left(l{\right)}^{*}\to {M}_{4}\left(l\right)$  is defined by sending any linear function $\phi :{M}_{4}\left(l\right)\to \mathbb{C}$  to the element of ${M}_{4}\left(l\right)$  given by
 $\begin{array}{c}\frac{1}{24}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v\right)-\left(y+x{\right)}^{2}\left(v,v-u\right)\right){x}^{2}\left(u,v\right)\end{array}$
 $\begin{array}{c}-2\phi \left(y\left(y-x\right)\left(-v,u+v\right)-\left(-y\right)\left(y+x\right)\left(v,v-u\right)\right)xy\left(u,v\right)\end{array}$
 $\begin{array}{c}+\phi \left({y}^{2}\left(-v,u+v\right)-\left(-y{\right)}^{2}\left(v,v-u\right)\right){y}^{2}\left(u,v\right)\right)\end{array}$
where we adopt the convention $\phi \left(P\left(x,y\right)\left(u,v\right)\right)=0$  for any $P\left(x,y\right)$  if $gcd\left(u,v,l\right)>1$  .
Proposition 3.2. The bilinear form on ${M}_{4}\left(l\right)$  induced by $PD$  is skew-symmetric. Namely, for any $\phi ,\lambda \in {M}_{4}\left(l{\right)}^{*}$  one has $\lambda \left(PD\left(\phi \right)\right)=-\phi \left(PD\left(\lambda \right)\right).$
• Proof. We can express $\lambda \left(PD\left(\phi \right)\right)$  as  $\begin{array}{c}\frac{1}{24}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v\right)-\left(y+x{\right)}^{2}\left(v,v-u\right)\right)\lambda \left({x}^{2}\left(u,v\right)\right)\end{array}$
 $\begin{array}{c}-2\phi \left(y\left(y-x\right)\left(-v,u+v\right)-\left(-y\right)\left(y+x\right)\left(v,v-u\right)\right)\lambda \left(xy\left(u,v\right)\right)\end{array}$
 $\begin{array}{c}+\phi \left({y}^{2}\left(-v,u+v\right)-\left(-y{\right)}^{2}\left(v,v-u\right)\right)\lambda \left({y}^{2}\left(u,v\right)\right)\right).\end{array}$
We then use the relations $P\left(x,y\right)\left(c,d\right)=-P\left(y,-x\right)\left(-d,c\right)$  to rewrite the terms with $\left(v,v-u\right)$  in terms of $\left(u-v,v\right)$  . Afterwards, we use the relation $P\left(x,y\right)\left(c,d\right)=-P\left(y-x,-x\right)\left(-c-d,c\right)-P\left(-y,x-y\right)\left(d,-c-d\right)$  to further rewrite them in terms of $\left(v,-u\right)$  and $\left(-u,u-v\right)$  . Finally, we use the relations $P\left(x,y\right)\left(c,d\right)=-P\left(y,-x\right)\left(d,-c\right)$  to write the result in terms of $\left(u,v\right)$  and $\left(u,u-v\right)$  . We also rewrite the terms with $\left(-v,u+v\right)$  in terms of $\left(u,v\right)$  and $\left(-u-v,u\right)$  . After simplifications, this gives  $\begin{array}{cc}\left(y-x{\right)}^{2}\left(-v,u+v\right)-\left(y+x{\right)}^{2}\left(v,v-u\right)=& {x}^{2}\left(u-v,u\right)\end{array}$
 $\begin{array}{cc}& -{x}^{2}\left(-u-v,u\right)\end{array}$
 $\begin{array}{cc}y\left(y-x\right)\left(-v,u+v\right)-\left(-y\right)\left(y+x\right)\left(v,v-u\right)=& -x\left(y+x\right)\left(u-v,u\right)\end{array}$
 $\begin{array}{cc}& -x\left(x-y\right)\left(-u-v,u\right)\end{array}$
 $\begin{array}{cc}{y}^{2}\left(-v,u+v\right)-\left(-y{\right)}^{2}\left(v,v-u\right)=& \left(y+x{\right)}^{2}\left(u-v,u\right)\end{array}$
 $\begin{array}{cc}& -\left(x-y{\right)}^{2}\left(-u-v,u\right)\end{array}$
which allows us to write $\lambda \left(PD\left(\phi \right)\right)$  as  $\begin{array}{c}\frac{1}{24}{\sum }_{u,v}\left(\phi \left({x}^{2}\left(u-v,u\right)-{x}^{2}\left(-u-v,u\right)\right)\lambda \left({x}^{2}\left(u,v\right)\right)\end{array}$
 $\begin{array}{c}-2\phi \left(-x\left(y+x\right)\left(u-v,u\right)-x\left(x-y\right)\left(-u-v,u\right)\right)\lambda \left(xy\left(u,v\right)\right)\end{array}$
 $\begin{array}{c}+\phi \left(\left(y+x{\right)}^{2}\left(u-v,u\right)-\left(x-y{\right)}^{2}\left(-u-v,u\right)\right)\lambda \left({y}^{2}\left(u,v\right)\right)\right).\end{array}$
It remains to switch the indexing in ${\sum }_{u,v}$  so that $\phi \left(...\right)$  becomes $\phi \left(...\left(u,v\right)\right)$  and simplify to get $-\phi \left(PD\left(\lambda \right)\right)$  .
Proposition 3.3. The map PD passes through the spaces of cusp symbols, namely there is a commutative diagram
 $\begin{array}{ccc}{M}_{4}\left(l{\right)}^{*}& \stackrel{PD}{\to }& {M}_{4}\left(l\right)\end{array}$
 $\begin{array}{ccc}↓& & ↑\end{array}$
 $\begin{array}{ccc}{S}_{4}\left(l{\right)}^{*}& \to & {S}_{4}\left(l\right)\end{array}$
with the side maps coming from the natural inclusions ${S}_{4}\left(l\right)\to {M}_{4}\left(l\right)$  .
• Proof. First, let ${\psi }_{\left(a,b\right)}:{M}_{4}\left(l\right)\to \mathbb{C}$  be the evaluation at the cusp $±\left(a,b\right)$  with $a\in \mathbb{Z}/l\mathbb{Z}$  and $b\in \mathbb{Z}/\left(a,l\right)\mathbb{Z}$  as in  2.2 . Let us show that ${\psi }_{\left(a,b\right)}\left(PD\left(\phi \right)\right)=0$  for any $\phi$  . Using  2.2 , we get $24{\psi }_{\left(a,b\right)}\left(PD\left(\phi \right)\right)=2\phi \left({\sum }_{v=bmod\left(a,l\right)}\left(\left(y-x{\right)}^{2}\left(-v,a+v\right)-\left(y+x{\right)}^{2}\left(v,v-a\right)\right)-{\sum }_{u=-bmod\left(a,l\right)}\left({y}^{2}\left(-a,u+a\right)-{y}^{2}\left(a,a-u\right)\right)\right)$  By switching from $u$  to $u\mp a$  in the two terms of the last sum, we see that it vanishes. For the first sum, we switch from $v$  to $v-a$  in the first term. Then we use $\left(x-y{\right)}^{2}\left(c,d\right)=-\left(x+y{\right)}^{2}\left(-d,c\right)$  to reduce it (up to a constant) to the value of $\phi$  on ${\sum }_{v=bmod\left(a,l\right)}\left(y-x{\right)}^{2}\left(a-v,v\right)$  . We now apply the relations on modular symbols to rewrite this as $-{\sum }_{v=bmod\left(a,l\right)}\left({x}^{2}\left(v,-a\right)+{y}^{2}\left(-a,a-v\right)\right)$  $=-{\sum }_{v=bmod\left(a,l\right)}\left({x}^{2}\left(v,-a\right)+{y}^{2}\left(a,v\right)\right)=0.$  This shows that the image of $PD$  sits inside ${S}_{4}\left(l\right)$  .
By Proposition  3.2 , we now see that $\phi \left(PD\left({\psi }_{a,b}\right)\right)=0$  for any $\phi$  . This shows that $PD$  passes through ${S}_{4}\left(l{\right)}^{*}$  which finishes the proof.
The key result of this section hinges on a formula for the Petersson inner product of cusp forms in terms of period integrals. Recall that the Petersson inner product of two holomorphic cusp forms of weight four with respect to ${\Gamma }_{1}\left(l\right)$  is defined as $\left(f,g{\right)}_{Petersson}=\int {\int }_{\mathcal{ℋ}/{\Gamma }_{1}\left(l\right)}f\left(\tau \right)\overline{g\left(\tau \right)}\Im \left(\tau {\right)}^{2}d\tau d\overline{\tau }.$  Period integrals define a pairing between the space ${\mathcal{S}}_{4}\left(l\right)$  of cusp forms of weight four and ${M}_{4}\left(l\right)$  , which we denote by $〈,〉$  . This pairing is a crucial feature of the theory of modular symbols, and we refer the reader to [M1for its definitions and properties.
Theorem 3.4. For any two holomorphic weight four forms $f$  and $g$  there holds
 $\begin{array}{cc}\left(1\right)& \left(f,g{\right)}_{Petersson}=-\frac{1}{24}{\sum }_{c,d\in \mathbb{Z}/l\mathbb{Z},gcd\left(c,d,l\right)=1}\end{array}$
 $\begin{array}{cc}& \left(\overline{\left(〈g,\left(y-x{\right)}^{2}\left(-d,c+d\right)〉-〈g,\left(y+x{\right)}^{2}\left(d,d-c\right)〉\right)}〈f,{x}^{2}\left(c,d\right)〉\end{array}$
 $\begin{array}{cc}& -2\overline{\left(〈g,y\left(y-x\right)\left(-d,c+d\right)〉-〈g,-y\left(y+x\right)\left(d,d-c\right)〉\right)}〈f,xy\left(c,d\right)〉\end{array}$
 $\begin{array}{cc}& +\overline{\left(〈g,{y}^{2}\left(-d,c+d\right)〉-〈g,\left(-y{\right)}^{2}\left(d,d-c\right)〉\right)}〈f,{y}^{2}\left(c,d\right)〉\right)\end{array}$
 $\begin{array}{cc}\left(2\right)& 0=-\frac{1}{24}{\sum }_{c,d\in \mathbb{Z}/l\mathbb{Z},gcd\left(c,d,l\right)=1}\end{array}$
 $\begin{array}{cc}& \left(\left(〈g,\left(y-x{\right)}^{2}\left(-d,c+d\right)〉-〈g,\left(y+x{\right)}^{2}\left(d,d-c\right)〉\right)〈f,{x}^{2}\left(c,d\right)〉\end{array}$
 $\begin{array}{cc}& -2\left(〈g,y\left(y-x\right)\left(-d,c+d\right)〉-〈g,-y\left(y+x\right)\left(d,d-c\right)〉\right)〈f,xy\left(c,d\right)〉\end{array}$
 $\begin{array}{cc}& +\left(〈g,{y}^{2}\left(-d,c+d\right)〉-〈g,\left(-y{\right)}^{2}\left(d,d-c\right)〉\right)〈f,{y}^{2}\left(c,d\right)〉\right).\end{array}$
• Proof. Consider the cosets ${\Gamma }_{1}\left(l\right)\lambda$  for $\lambda \in S{l}_{2}\left(\mathbb{Z}\right)$  . Then the fundamental domain $\mathcal{ℋ}/{\Gamma }_{1}\left(l\right)$  can be chosen as ${\cup }_{{\Gamma }_{1}\left(l\right)\lambda }\lambda \left({D}_{0}\right)$  for any fundamental domain ${D}_{0}$  of ${\Gamma }_{1}\left(l\right)$  . Moreover, we can use a union of three different such ${D}_{0}$  to write the Petersson pairing as $\left(f,g{\right)}_{Petersson}=\frac{1}{3}{\sum }_{{\Gamma }_{1}\left(l\right)\lambda }\int {\int }_{\lambda \left(D\right)}f\left(\tau \right)\overline{g\left(\tau \right)}\Im \left(\tau {\right)}^{2}d\tau d\overline{\tau }$  where $D$  is the geodesic triangle in $\mathcal{ℋ}$  with vertices $i\infty$  , $-1$  and $0$  . The boundary of $D$  consists of the vertical lines $\Re \left(\tau \right)=0$  and $\Re \left(\tau \right)=-1$  , as well as the upper half of the circle of radius $\frac{1}{2}$  centered at $-\frac{1}{2}$  .
We use $\Im \left(\lambda \left(\tau \right){\right)}^{2}=\Im \left(\tau {\right)}^{2}|c\tau +d{|}^{-4}$  where $\lambda \left(\tau \right)=\frac{a\tau +b}{c\tau +d}$  to rewrite each term of the above sum as $\int {\int }_{D}f\left(\lambda \left(\tau \right)\right)\overline{g\left(\lambda \left(\tau \right)\right)}\Im \left(\tau {\right)}^{2}|c\tau +d{|}^{-8}d\tau d\overline{\tau }.$  For each such $\lambda$  we introduce for $i=0,1,2$  ${G}_{i,\lambda }\left(\tau \right)=\overline{{\int }_{-1}^{\tau }g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}{s}^{i}ds}$  and ${f}_{i,\lambda }\left(\tau \right)\right)=f\left(\lambda \left(\tau \right)\right)\left(c\tau +d{\right)}^{-4}{\tau }^{i}$  . Then we write $\Im \left(\tau {\right)}^{2}=-\frac{1}{4}\left(\tau -\overline{\tau }{\right)}^{2}$  and use Stokes's Theorem to derive $\left(f,g{\right)}_{Petersson}=-\frac{1}{12}{\sum }_{{\Gamma }_{1}\left(l\right)\lambda }{\int }_{\partial D}\left({G}_{0,\lambda }\left(\tau \right){f}_{2,\lambda }\left(\tau \right)-2{G}_{1,\lambda }\left(\tau \right){f}_{1,\lambda }\left(\tau \right)+{G}_{2,\lambda }\left(\tau \right){f}_{0,\lambda }\left(\tau \right)\right)d\tau .$  The boundary of $D$  splits into three geodesics, and our first claim is that the terms of the integration for the ${\int }_{i\infty }^{-1}$  and ${\int }_{-1}^{0}$  of $\partial D$  cancel each other.
Consider the map $\sigma \left(\tau \right):=-\frac{1}{\tau +2}$  . Element $\sigma \in S{l}_{2}\left(\mathbb{Z}\right)$  acts on the set of cosets ${\Gamma }_{1}\left(l\right)\lambda$  by right multiplication. We observe that ${\int }_{-1}^{0}\left({G}_{0,\lambda \sigma }\left(\tau \right){f}_{2,\lambda \sigma }\left(\tau \right)-2{G}_{1,\lambda \sigma }\left(\tau \right){f}_{1,\lambda \sigma }\left(\tau \right)+{G}_{2,\lambda \sigma }\left(\tau \right){f}_{0,\lambda \sigma }\left(\tau \right)\right)d\tau$  $={\int }_{-1}^{0}\left({G}_{0,\lambda }\left(\sigma \left(\tau \right)\right){f}_{2,\lambda }\left(\sigma \left(\tau \right)\right)-2{G}_{1,\lambda }\left(\sigma \left(\tau \right)\right){f}_{1,\lambda }\left(\sigma \left(\tau \right)\right)$  $+{G}_{2,\lambda }\left(\sigma \left(\tau \right)\right){f}_{0,\lambda }\left(\sigma \left(\tau \right)\right)\right)d\sigma \left(\tau \right)$  $={\int }_{-1}^{i\infty }\left({G}_{0,\lambda }\left(\tau \right){f}_{2,\lambda }\left(\tau \right)-2{G}_{1,\lambda }\left(\tau \right){f}_{1,\lambda }\left(\tau \right)+{G}_{2,\lambda }\left(\tau \right){f}_{0,\lambda }\left(\tau \right)\right)d\tau .$  The first equality is verified by a lengthy but straightforward calculation, which is left to the reader, since we will perform a similar calculation below.
It is crucial that we chose $\left(-1\right)$  as the lower limit of integration in the definition of ${G}_{i,\lambda }$  and that $\sigma$  preserves $\left(-1\right)$  .
So now we are left with $\left(f,g{\right)}_{Petersson}=-\frac{1}{12}{\sum }_{{\Gamma }_{1}\left(l\right)\lambda }{\int }_{0}^{i\infty }\left({G}_{0,\lambda }\left(\tau \right){f}_{2,\lambda }\left(\tau \right)-2{G}_{1,\lambda }\left(\tau \right){f}_{1,\lambda }\left(\tau \right)+{G}_{2,\lambda }\left(\tau \right){f}_{0,\lambda }\left(\tau \right)\right)d\tau .$  We will do a similar trick, this time with $\nu \left(\tau \right)=-\frac{1}{\tau }$  instead of $\sigma \left(\tau \right)$  .
It has an effect of switching the direction of integration, but since it does not preserve $\left(-1\right)$  , the functions ${G}_{i,\lambda }$  acquire extra additive terms. More specifically, one has ${G}_{i,\lambda \nu }\left(\tau \right)=\overline{{\int }_{-1}^{\tau }g\left(\lambda \left(-\frac{1}{s}\right)\right)\left(ds-c{\right)}^{-4}{s}^{i}ds}=\overline{{\int }_{1}^{\nu \left(\tau \right)}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}{s}^{2-i}\left(-1{\right)}^{i}ds}$  $=\left(-1{\right)}^{i}{G}_{2-i,\lambda }\left(\nu \left(\tau \right)\right)+\overline{{\int }_{1}^{-1}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}{s}^{2-i}\left(-1{\right)}^{i}ds}$  ${f}_{i,\lambda \nu }\left(\tau \right)d\tau =\left(-1{\right)}^{i}{f}_{2-i,\lambda }\left(\nu \left(\tau \right)\right)d\nu \left(\tau \right).$  We rewrite $\left(f,g{\right)}_{Petersson}$  as $-\frac{1}{24}\left({\sum }_{{\Gamma }_{1}\left(l\right)\lambda }{\int }_{0}^{i\infty }\left({G}_{0,\lambda }\left(\tau \right){f}_{2,\lambda }\left(\tau \right)-2{G}_{1,\lambda }\left(\tau \right){f}_{1,\lambda }\left(\tau \right)+{G}_{2,\lambda }\left(\tau \right){f}_{0,\lambda }\left(\tau \right)\right)d\tau +{\sum }_{{\Gamma }_{1}\left(l\right)\lambda }{\int }_{0}^{i\infty }\left({G}_{0,\lambda \nu }\left(\tau \right){f}_{2,\lambda \nu }\left(\tau \right)-2{G}_{1,\lambda \nu }\left(\tau \right){f}_{1,\lambda \nu }\left(\tau \right)+{G}_{2,\lambda \nu }\left(\tau \right){f}_{0,\lambda \nu }\left(\tau \right)\right)d\tau \right)$  which together with transformation formulas for $G$  and $f$  implies, after cancelling $±{\int }_{0}^{i\infty }$  $\left(f,g{\right)}_{Petersson}=-\frac{1}{24}{\sum }_{{\Gamma }_{1}\left(l\right)\lambda }\left(\overline{{\int }_{1}^{-1}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}{s}^{2}ds}{\int }_{0}^{i\infty }{f}_{0,\lambda }\left(\tau \right)d\tau -2\overline{{\int }_{1}^{-1}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}sds}{\int }_{0}^{i\infty }{f}_{1,\lambda }\left(\tau \right)d\tau +\overline{{\int }_{1}^{-1}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}ds}{\int }_{0}^{i\infty }{f}_{2,\lambda }\left(\tau \right)d\tau \right).$  It remains to write ${\int }_{1}^{-1}$  in terms of pairings with modular symbols by writing the arc from $1$  to $\left(-1\right)$  in terms of the unimodular arcs from $1$  to $i\infty$  and from $i\infty$  to $\left(-1\right)$  . Namely, for a homogeneous degree two polynomial $P\left(x,y\right)$  one has ${\int }_{1}^{-1}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}P\left(s,1\right)ds$  $={\int }_{1}^{0}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}P\left(s,1\right)ds-{\int }_{-1}^{0}g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}P\left(s,1\right)ds$  $={\int }_{0}^{i\infty }g\left(\lambda \left(\frac{1}{1-s}\right)\right)\left(-ds+\left(c+d\right){\right)}^{-4}P\left(1,1-s\right)ds$  $-{\int }_{0}^{i\infty }g\left(\lambda \left(-\frac{1}{1+s}\right)\right)\left(ds+\left(d-c\right){\right)}^{-4}P\left(-1,1+s\right)dt$  $=〈g,P\left(y,y-x\right)\left(-d,c+d\right)〉-〈g,P\left(-y,y+x\right)\left(d,d-c\right)〉.$  Finally, one observes that cosets ${\Gamma }_{1}\left(l\right)\lambda$  are in one-to-one correspondence with pairs $\left(c,d\right)$  with $gcd\left(c,d,l\right)=1$  , and the first claim of the theorem follows.
The second claim of the theorem is proved by the same technique. This time we define ${G}_{i,\lambda }\left(\tau \right)$  as ${\int }_{-1}^{\tau }g\left(\lambda \left(s\right)\right)\left(cs+d{\right)}^{-4}{s}^{i}ds$  . Consequently, it is holomorphic, and the Stokes's Theorem gives $0$  instead of the Petersson product. The rest of the calculations are unchanged.
Remark 3.5. Similar formulas for Petersson product are already present in the literature. They seem to go back to at least as far as [Hand [KZ. We learned the argument (in weight two case) from [M2. In addition to extending it to weight four, we streamlined it just slightly by looking at the union of three fundamental domains for $S{l}_{2}\left(\mathbb{Z}\right)$  , rather than one. This allowed us to avoid integration between elliptic points.
Corollary 3.6. The pairing on ${S}_{4}\left(l\right)$  induced by $PD$  is nondegenerate.
• Proof. By [Sh, the integration pairing is a perfect pairing between ${S}_{4}\left(l\right)$  and the direct sum ${V}_{hol}\oplus {\overline{V}}_{hol}$  of the spaces of holomorphic and anti-holomorphic cusp forms. Every element $\phi \in {S}_{4}\left(l{\right)}^{*}$  can therefore be written as $〈f,\cdot 〉+〈\overline{g},\cdot 〉$  . Suppose $PD\left(\phi \right)=0$  . Denote by $\overline{}$  the anti-isomorphism of ${M}_{4}\left(l\right)$  that sends $\alpha {x}^{i}{y}^{2-i}\left(u,v\right)$  to $\overline{\alpha }{x}^{i}{y}^{2-i}\left(u,v\right)$  . Then Theorem  3.4 shows that $0=\phi \left(\overline{PD\left(\phi \right)}\right)=-〈f,f{〉}_{Petersson}-〈g,g{〉}_{Petersson}\le 0$  with equality only if $f=g=0$  .
Remark 3.7. The arguments of our proof of Theorem  3.4 extend naturally to arbitrary integer weights $k\ge 2$  and arbitrary subgroups $\Gamma$  of finite index in $S{l}_{2}\left(\mathbb{Z}\right)$  . We expect Corollary  3.6 to extend to arbitrary weight and to arbitrary group $\Gamma$  , after an appropriate definition of $PD$  . One would need to interpret the arguments of Propositions  3.2 and  3.3 to extend them to this more general setting. For instance, we expect that the bilinear form on ${S}_{k}\left(\Gamma \right)$  induced by $PD$  is $\left(-1{\right)}^{k+1}$  -symmetric. We believe that maps $PD$  can be interpreted as an intersection pairings in the middle cohomology of the Kuga varieties, although we do not need this for the purposes of this paper. Nevertheless, this is why we refer to $PD$  as the Poincaré duality map.
The map $PD$  behaves well with respect to the involution $i$  . We denote by ${M}_{4}\left(l{\right)}_{±}^{*}$  the eigenspaces of the dual involution ${i}^{*}$  on ${M}_{4}\left(l{\right)}^{*}$  .
Proposition 3.8. $PD\left({M}_{4}\left(l{\right)}_{±}^{*}\right)\subseteq {M}_{4}\left(l{\right)}_{\mp }$  .
• Proof. If $\phi \in {M}_{4}\left(l{\right)}_{±}^{*}$  , then for any modular symbol $P\left(x,y\right)\left(u.v\right)$  there holds $\phi \left(P\left(x,y\right)\left(u,v\right)\right)=\phi \left(P\left(x,y\right)\left(u,v{\right)}_{±}\right)$  . Consequently, $PD\left(\phi \right)=\frac{1}{24}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v{\right)}_{±}-\left(y+x{\right)}^{2}\left(v,v-u{\right)}_{±}\right){x}^{2}\left(u,v\right)$  $-2\phi \left(y\left(y-x\right)\left(-v,u+v{\right)}_{±}-\left(-y\right)\left(y+x\right)\left(v,v-u{\right)}_{±}\right)xy\left(u,v\right)$  $+\phi \left({y}^{2}\left(-v,u+v{\right)}_{±}-\left(-y{\right)}^{2}\left(v,v-u{\right)}_{±}\right){y}^{2}\left(u,v\right)\right)$  $=\frac{1}{24}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v{\right)}_{±}\mp \left(y-x{\right)}^{2}\left(-v,v-u{\right)}_{±}\right){x}^{2}\left(u,v\right)$  $-2\phi \left(y\left(y-x\right)\left(-v,u+v{\right)}_{±}\mp \left(-y\right)\left(y-x\right)\left(-v,v-u{\right)}_{±}\right)xy\left(u,v\right)$  $+\phi \left({y}^{2}\left(-v,u+v{\right)}_{±}\mp \left(-y{\right)}^{2}\left(-v,v-u{\right)}_{±}\right){y}^{2}\left(u,v\right)\right)$  $=\frac{1}{24}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v{\right)}_{±}\right)\left({x}^{2}\left(u,v\right)\mp {x}^{2}\left(-u,v\right)\right)$  $-2\phi \left(y\left(y-x\right)\left(-v,u+v{\right)}_{±}\right)\left(xy\left(u,v\right)±xy\left(-u,v\right)\right)$  $+\phi \left({y}^{2}\left(-v,u+v{\right)}_{±}\right)\left({y}^{2}\left(u,v\right)\mp {y}^{2}\left(-u,v\right)\right)\right)$  $=\frac{1}{12}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v{\right)}_{±}\right){x}^{2}\left(u,v{\right)}_{\mp }$  $-2\phi \left(y\left(y-x\right)\left(-v,u+v{\right)}_{±}\right)xy\left(u,v{\right)}_{\mp }+\phi \left({y}^{2}\left(-v,u+v{\right)}_{±}\right){y}^{2}\left(u,v{\right)}_{\mp }\right).$
Remark 3.9. In what follows, we will abuse the notations somewhat to denote the induced map ${M}_{4}\left(l{\right)}_{-}^{*}\to {M}_{4}\left(l{\right)}_{+}$  by $PD$  as well. By Proposition  3.3 , this map comes from a map ${S}_{4}\left(l{\right)}_{-}^{*}\to {S}_{4}\left(l{\right)}_{+}$  .
Corollary 3.10. The induced map $PD:{S}_{4}\left(l{\right)}_{-}^{*}\to {S}_{4}\left(l{\right)}_{+}$  is an isomorphism.
• Proof. Combine Corollary  3.6 and Proposition  3.8 .

4 The Wronskian map

In this section we define the map from the modular symbols of weight four to cusp forms of weight four. First, we need to define the Eisenstein series ${s}_{a}\left(q\right)$  , ${t}_{a}\left(q\right)$  and ${r}_{a}\left(q\right)$  for $a\in \mathbb{Z}/l\mathbb{Z}$  . Our notation for ${s}_{a}$  differs from that of [BG1by a Fricke involution. We recall that quasimodular forms of weight two are linear combinations of the usual modular forms of weight two and the (slightly non-modular) Eisenstein series ${E}_{2}\left(q\right)$  of weight $2$  . In weight one, quasimodular forms are modular.
Proposition 4.1. For each $amodl$  there exist ${\Gamma }_{1}\left(l\right)$  -quasimodular forms ${s}_{a}\left(q\right)$  , ${t}_{a}\left(q\right)$  and ${r}_{a}\left(q\right)$  of weights $1$  , $2$  and $2$  respectively given by ${s}_{a}\left(q\right)=\left(\frac{1}{2}-\left\{\frac{a}{l}\right\}\right)+{\sum }_{n>0}{q}^{n}{\sum }_{d|n}\left({\delta }_{d}^{amodl}-{\delta }_{d}^{-amodl}\right),$  $ifa\ne 0modl,{s}_{0}\left(q\right)=0,$  ${t}_{a}\left(q\right)=constant+{\sum }_{n}{q}^{n}{\sum }_{d|n}\frac{n}{k}\left({\delta }_{d}^{amodl}+{\delta }_{d}^{-amodl}\right)$  ${r}_{a}\left(q\right)=constant+{\sum }_{n}{q}^{n}{\sum }_{d|n}d\left({\delta }_{d}^{amodl}+{\delta }_{d}^{-amodl}\right),$  where the exact values of the constants depend on $a$  and $l$  and are determined uniquely by the quasimodularity.
• Proof. These series are obtained as linear combinations of the weight one and two Eisenstein series considered in [BG1. Details are left to the reader.
Definition 4.2. We define the map $\mu :{S}_{4}\left(l\right)\to {\mathcal{S}}_{4}\left(l\right)$  by the formula ${x}^{2}\left(u,v\right)↦-2{t}_{u}{r}_{v}-\frac{1}{l}q\frac{\partial }{\partial q}{r}_{v}-{\delta }_{v}^{0modl}q\frac{\partial }{\partial q}{t}_{u},xy\left(u,v\right)↦\frac{1}{2\pi i}W\left({s}_{u},{s}_{v}\right),$  ${y}^{2}\left(u,v\right)↦2{r}_{u}{t}_{v}+\frac{1}{l}q\frac{\partial }{\partial q}{r}_{u}+{\delta }_{u}^{0modl}q\frac{\partial }{\partial q}{t}_{v}.$  Note that the ring of quasimodular forms is closed under $q\frac{\partial }{\partial q}$  . We will show in Theorem  4.3 below that $\mu$  is well-defined, i.e. it is compatible with the relations on modular symbols.
Theorem 4.3. The map $\mu$  of definition  4.2 is well-defined.
• Proof. We need to check that $\mu$  maps the relations  2.1 to zero. We have $\mu \left({x}^{2}\left(u,v\right)+{y}^{2}\left(v,-u\right)\right)=0$  and $\mu \left(xy\left(u,v\right)-xy\left(v,-u\right)\right)=\frac{1}{2\pi i}\left(W\left({s}_{u},{s}_{v}\right)-W\left({s}_{v},{s}_{-u}\right)\right)=0$  by the symmetry properties ${s}_{-a}=-{s}_{a}$  , ${r}_{a}={r}_{-a}$  and ${t}_{a}={t}_{-a}$  .
The difficult part is to show that $\mu$  maps $R=xy\left(v,-u-v\right)-xy\left(-u-v,u\right)+{y}^{2}\left(-u-v,u\right)+{x}^{2}\left(u,v\right)-xy\left(u,v\right)$  to zero. For each positive integer $n$  let us denote by $I\left(n\right)$  the set of fourtuples $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)\in {\mathbb{Z}}_{>0}^{4}$  that satisfy ${m}_{1}{k}_{1}+{m}_{2}{k}_{2}=n$  . Let us denote by $\sim$  the equality of power series in $q$  up to linear combinations of quasimodular forms of weights $0$  , $1$  , $2$  , and the derivatives of ${s}_{a}\left(q\right)$  with respect to $\tau$  . This allows us to avoid looking at the specific values of the constant terms in Proposition  4.1 .
We have $\mu \left(R\right)\sim -\frac{1}{l}q\frac{\partial }{\partial q}{r}_{v}-{\delta }_{v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+{\delta }_{u+v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+\frac{1}{l}q\frac{\partial }{\partial q}{r}_{u+v}$  $+{\sum }_{n>0}{q}^{n}{\sum }_{{I}_{n}}\left(\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right)\left({\delta }_{{k}_{1}}^{vmodl}-{\delta }_{{k}_{1}}^{-vmodl}\right)\left({\delta }_{{k}_{2}}^{-u-vmodl}-{\delta }_{{k}_{2}}^{u+vmodl}\right)-\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right)\left({\delta }_{{k}_{1}}^{-u-vmodl}-{\delta }_{{k}_{1}}^{u+vmodl}\right)\left({\delta }_{{k}_{2}}^{umodl}-{\delta }_{{k}_{2}}^{-umodl}\right)+2{k}_{1}{m}_{2}\left({\delta }_{{k}_{1}}^{-u-vmodl}+{\delta }_{{k}_{1}}^{u+vmodl}\right)\left({\delta }_{{k}_{2}}^{umodl}+{\delta }_{{k}_{2}}^{-umodl}\right)-2{m}_{1}{k}_{2}\left({\delta }_{{k}_{1}}^{umodl}+{\delta }_{{k}_{1}}^{-umodl}\right)\left({\delta }_{{k}_{2}}^{vmodl}+{\delta }_{{k}_{2}}^{-vmodl}\right)-\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right)\left({\delta }_{{k}_{1}}^{umodl}-{\delta }_{{k}_{1}}^{-umodl}\right)\left({\delta }_{{k}_{2}}^{vmodl}-{\delta }_{{k}_{2}}^{-vmodl}\right)\right).$  We introduce the notation ${A}_{{k}_{1},{k}_{2}}={\delta }_{{k}_{1}}^{umodl}\delta {{k}_{2}}^{vmodl}+{\delta }_{{k}_{1}}^{-umodl}{\delta }_{{k}_{2}}^{-vmodl}$  to rewrite the above as $\mu \left(R\right)\sim -\frac{1}{l}q\frac{\partial }{\partial q}{r}_{v}-{\delta }_{v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+{\delta }_{u+v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+\frac{1}{l}q\frac{\partial }{\partial q}{r}_{u+v}$  $+{\sum }_{n>0}{q}^{n}{\sum }_{I\left(n\right)}\left(\left({m}_{2}{k}_{2}-{m}_{1}{k}_{1}-2{m}_{1}{k}_{2}\right){A}_{{k}_{1},{k}_{2}}+\left({m}_{2}{k}_{2}-{m}_{1}{k}_{1}+2{k}_{1}{m}_{2}\right){A}_{{k}_{2},-{k}_{1}-{k}_{2}}+\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right){A}_{-{k}_{1}-{k}_{2},{k}_{1}}-\left({m}_{2}{k}_{2}-{m}_{1}{k}_{1}+2{m}_{1}{k}_{2}\right){A}_{-{k}_{1},{k}_{2}}-\left({m}_{2}{k}_{2}-{m}_{1}{k}_{1}-2{k}_{1}{m}_{2}\right){A}_{{k}_{2},{k}_{1}-{k}_{2}}-\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right){A}_{{k}_{1}-{k}_{2},-{k}_{1}}\right)$  We recall (see [BG2) that $I\left(n\right)$  is a disjoint union of the runs of Euclid algorithm. The algorithm is given by the partially defined map $up:I\left(n\right)\to I\left(n\right)$  $up:\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)↦\left\{\begin{array}{cc}\left({m}_{2},{k}_{1}+{k}_{2},{m}_{1}-{m}_{2},{k}_{1}\right),& {m}_{1}>{m}_{2}\\ \left({m}_{2}-{m}_{1},{k}_{2},{m}_{1},{k}_{1}+{k}_{2}\right),& {m}_{1}<{m}_{2}.\end{array}$  Repeated applications of this map go from the subset of $I\left(n\right)$  with ${k}_{1}={k}_{2}$  to the subset of $I\left(n\right)$  with ${m}_{1}={m}_{2}$  , where $up$  is not defined. As in [BG2, we will show that for each run of the algorithm the above sum is telescoping. Namely, the ”plus” terms with ${A}_{{k}_{1},{k}_{2}}$  , ${A}_{{k}_{2},-{k}_{1}-{k}_{2}}$  , ${A}_{-{k}_{1}-{k}_{2},{k}_{1}}$  for $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)$  cancel the ”minus” terms with ${A}_{-{k}_{1},{k}_{2}}$  , ${A}_{{k}_{2},{k}_{1}-{k}_{2}}$  , ${A}_{{k}_{1}-{k}_{2},-{k}_{1}}$  for $up\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)$  . There are two cases to check, depending on whether ${m}_{1}>{m}_{2}$  or ${m}_{1}<{m}_{2}$  . In the case of ${m}_{1}>{m}_{2}$  the ”minus” terms for $up\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)=\left({m}_{2},{k}_{1}+{k}_{2},{m}_{1}-{m}_{2},{k}_{1}\right)$  are $-\left(\left({m}_{1}-{m}_{2}\right){k}_{1}-{m}_{2}\left({k}_{1}+{k}_{2}\right)+2{m}_{2}{k}_{1}\right){A}_{-{k}_{1}-{k}_{2},{k}_{1}}$  $-\left(\left({m}_{1}-{m}_{2}\right){k}_{1}-{m}_{2}\left({k}_{1}+{k}_{2}\right)-2\left({k}_{1}+{k}_{2}\right)\left({m}_{1}-{m}_{2}\right)\right){A}_{{k}_{1},{k}_{2}}$  $-\left({m}_{2}{k}_{1}+{m}_{2}{k}_{2}\right)-\left({m}_{1}-{m}_{2}\right){k}_{1}\right){A}_{{k}_{2},-{k}_{1}-{k}_{2}}$  $=-\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right){A}_{-{k}_{1}-{k}_{2},{k}_{1}}-\left(-{k}_{1}{m}_{1}+{k}_{2}{m}_{2}-2{k}_{2}{m}_{1}\right){A}_{{k}_{1},{k}_{2}}$  $-\left({m}_{2}{k}_{2}-{m}_{1}{k}_{1}+2{m}_{2}{k}_{1}\right){A}_{{k}_{2},-{k}_{1}-{k}_{2}}$  which is seen to equal the ”plus” terms for $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)$  . The case of ${m}_{1}<{m}_{2}$  is similar and left to the reader. One needs to use there the symmetry ${A}_{{k}_{1},{k}_{2}}={A}_{-{k}_{1},-{k}_{2}}$  . Consequently, the only terms that will not be cancelled are the ”plus” terms for the subset of $I\left(n\right)$  with ${m}_{1}={m}_{2}$  and the ”minus” terms for the subset of $I\left(n\right)$  with ${k}_{1}={k}_{2}$  . This gives $\mu \left(R\right)\sim -\frac{1}{l}q\frac{\partial }{\partial q}{r}_{v}-{\delta }_{v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+{\delta }_{u+v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+\frac{1}{l}q\frac{\partial }{\partial q}{r}_{u+v}$  $+{\sum }_{n>0}{q}^{n}\left({\sum }_{{m}_{1},{k}_{1},{k}_{2}>0{m}_{1}\left({k}_{1}+{k}_{2}\right)=n}\left(-n{A}_{{k}_{1},{k}_{2}}+n{A}_{{k}_{2},-{k}_{1}-{k}_{2}}+{m}_{1}\left({k}_{1}-{k}_{2}\right){A}_{-{k}_{1}-{k}_{2},{k}_{1}}\right)-{\sum }_{{m}_{1},{m}_{2},{k}_{1}>0\left({m}_{1}+{m}_{2}\right){k}_{1}=n}\left(n{A}_{-{k}_{1},{k}_{1}}-n{A}_{{k}_{1},0}+\left({m}_{1}-{m}_{2}\right){k}_{1}{A}_{0,-{k}_{1}}\right)\right)$  $=-\frac{1}{l}q\frac{\partial }{\partial q}{r}_{v}-{\delta }_{v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+{\delta }_{u+v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+\frac{1}{l}q\frac{\partial }{\partial q}{r}_{u+v}$  $+{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}\left({\sum }_{0  We observe that ${\sum }_{0  and ${\sum }_{0  $+\left(1-\left\{\frac{d-u}{l}\right\}-\left\{-\frac{u}{l}\right\}\right).$  Consequently, $\mu \left(R\right)\sim -\frac{1}{l}q\frac{\partial }{\partial q}{r}_{v}-{\delta }_{v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+{\delta }_{u+v}^{0modl}q\frac{\partial }{\partial q}{t}_{u}+\frac{1}{l}q\frac{\partial }{\partial q}{r}_{u+v}$  $+{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}\left(-{\delta }_{u}^{0modl}\left({\delta }_{d}^{vmodl}+{\delta }_{d}^{-vmodl}\right)-{\delta }_{d}^{u+vmodl}\left(\frac{d}{l}-\left\{\frac{v}{l}\right\}+\left\{-\frac{u}{l}\right\}\right)-{\delta }_{d}^{-u-vmodl}\left(\frac{d}{l}-\left\{-\frac{v}{l}\right\}+\left\{\frac{u}{l}\right\}\right)+{\delta }_{d}^{vmodl}\left(\frac{d}{l}-\left\{-\frac{u}{l}\right\}+\left\{-\frac{u+v}{l}\right\}\right)+{\delta }_{d}^{-vmodl}\left(\frac{d}{l}-\left\{\frac{u}{l}\right\}+\left\{\frac{u+v}{l}\right\}\right)+\frac{n}{d}{\delta }_{v}^{0modl}\left({\delta }_{d}^{umodl}+{\delta }_{d}^{-umodl}\right)-\frac{n}{d}{\delta }_{u+v}^{0modl}\left({\delta }_{d}^{vmodl}+{\delta }_{d}^{-vmodl}\right)+{\delta }_{d}^{-umodl}\cdot \cdot \left(\frac{l}{d}\left(\left\{\frac{-u-v}{l}{\right\}}^{2}-\left\{\frac{-u-v}{l}\right\}-\left\{-\frac{v}{l}{\right\}}^{2}+\left\{-\frac{v}{l}\right\}\right)+\left(1-\left\{\frac{-u-v}{l}\right\}-\left\{-\frac{v}{l}\right\}\right)\right)+{\delta }_{d}^{umodl}\left(\frac{l}{d}\left(\left\{\frac{u+v}{l}{\right\}}^{2}-\left\{\frac{u+v}{l}\right\}-\left\{\frac{v}{l}{\right\}}^{2}+\left\{\frac{v}{l}\right\}\right)+\left(1-\left\{\frac{u+v}{l}\right\}-\left\{\frac{v}{l}\right\}\right)\right)\right).$  We use $\left\{t\right\}+\left\{-t\right\}=1-{\delta }_{t}^{0mod1}$  and, after tedious but straightforward calculations, get $\mu \left(R\right)\sim 0.$  Since $\mu \left(R\right)$  is quasimodular of weight four and $\sim$  is equality modulo forms of weight less than four, we get $\mu \left(R\right)=0$  .
Remark 4.4. There is a natural projection map from the space of quasimodular forms of weight four to the space of modular forms of weight four, which sends all forms divisible by ${E}_{2}$  to zero. So one can compose $\mu$  with this projection and have a map ${\mu }_{1}$  to the space of modular forms of weight four.
Proposition 4.5. Map $\mu$  sends ${M}_{4}\left(l{\right)}_{-}$  to zero.
• Proof. The statement immediately follows from the symmetry properties of $r$  , $s$  and $t$  .
Proposition 4.6. The image of ${S}_{4}\left(l{\right)}_{+}$  under $\mu$  is a subspace of ${\mathcal{S}}_{4}\left(l\right)$  which is the linear span of $W\left({s}_{a},{s}_{b}\right)$  with $gcd\left(a,b,l\right)=1$  .
• Proof. By Proposition  2.2 , $\mu \left({S}_{4}\left(l\right)\right)$  is spanned by $\mu \left(xy\left(a,b\right)\right)=W\left({s}_{a},{s}_{b}\right)$  for all $gcd\left(a,b,l\right)=1$  . By Proposition  4.5 , $\mu \left({S}_{4}\left(l\right){\right)}_{+}=\mu \left({S}_{4}\left(l\right)\right)$  .

5 The composition map

In this section we calculate the composition of the duality map $PD:{S}_{4}\left(l{\right)}_{-}^{*}\to {S}_{4}\left(l{\right)}_{+}$  of Section  3 and the Wronskian map $\mu$  of Section  4 . Our arguments are purely elementary. The result of this calculation will be used in the next section.
We need to introduce some additional notation. For any $\phi \in {S}_{4}\left(l{\right)}_{-}^{*}$  we will set $\phi \left(P\left(x,y\right)\left(u,v{\right)}_{-}\right)=0$  if $gcd\left(u,v,l\right)>1$  . We will also use the notation $\sim$  , namely, $f\sim g$  would mean that $f-g$  is a linear combination of quasimodular forms of weight at most two and the quasimodular forms of weight three that are derivatives of ${s}_{u}\left(\tau \right)$  . Finally, for every $n>0$  we introduce the set $H\left(n\right)$  of fourtuples of integers $\left(a,b,c,d\right)$  that satisfy $ad-bc=n,a>b\ge 0,d>c\ge 0$  .
Proposition 5.1. For any $\phi \in {S}_{4}\left(l{\right)}_{-}^{*}$  there holds $\mu \circ PD\left(\phi \right)\sim {\sum }_{n>0}{q}^{n}{\sum }_{H\left(n\right)}\phi \left(\left(ax+by\right)\left(cx+dy\right)\left(c,d{\right)}_{-}\right).$
• Proof. We will use the notations ${A}_{{k}_{1},{k}_{2}}$  and $I\left(n\right)$  from Section  4 . We use the last identity in the proof of Proposition  3.8 to get $12\mu \circ PD\left(\phi \right)=\mu \left({\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\left(\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v{\right)}_{-}\right){x}^{2}\left(u,v{\right)}_{+}$  $-2\phi \left(y\left(y-x\right)\left(-v,u+v{\right)}_{-}\right)xy\left(u,v{\right)}_{+}+\phi \left({y}^{2}\left(-v,u+v{\right)}_{-}\right){y}^{2}\left(u,v{\right)}_{+}\right)\right)$  $\sim 2{\sum }_{n>0}{q}^{n}\phi \left({\sum }_{I\left(n\right)}\left(\left(-2{m}_{1}{k}_{2}\left(y-x{\right)}^{2}-2\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right)y\left(y-x\right)+2{m}_{2}{k}_{1}{y}^{2}\right)\left(-{k}_{2},{k}_{1}+{k}_{2}{\right)}_{-}+\left(-2{m}_{1}{k}_{2}\left(y-x{\right)}^{2}+2\left({m}_{1}{k}_{1}-{m}_{2}{k}_{2}\right)y\left(y-x\right)+2{m}_{2}{k}_{1}{y}^{2}\right)\left(-{k}_{2},-{k}_{1}+{k}_{2}{\right)}_{-}\right)\right)\right)$  $-\frac{1}{l}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\phi \left(\left(y-x{\right)}^{2}\left(-v,u+v{\right)}_{-}-{y}^{2}\left(-u,v+u{\right)}_{-}\right)q\frac{\partial {r}_{v}}{\partial q}$  $-{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}\phi \left(\left(y-x{\right)}^{2}\left(0,u{\right)}_{-}-{y}^{2}\left(-u,u{\right)}_{-}\right)q\frac{\partial {t}_{u}}{\partial q}.$  Let us simplify the last two lines of the above equation. We use ${y}^{2}\left(-u,v+u{\right)}_{-}=-{x}^{2}\left(-v,-u{\right)}_{-}-\left(x-y{\right)}^{2}\left(v+u,-v{\right)}_{-}$  , ${y}^{2}\left(-u,u{\right)}_{-}=-{x}^{2}\left(0,-u{\right)}_{-}-\left(x-y{\right)}^{2}\left(u,0{\right)}_{-}$  , ${x}^{2}\left(v,u{\right)}_{-}+{x}^{2}\left(v,-u{\right)}_{-}={y}^{2}\left(v,u{\right)}_{-}+{y}^{2}\left(v,-u{\right)}_{-}=0$  and ${x}^{2}\left(0,u{\right)}_{-}={y}^{2}\left(0,u{\right)}_{-}=0$  and $xy\left(0,u{\right)}_{-}=xy\left(u,0{\right)}_{-}$  to rewrite them as $\frac{4}{l}{\sum }_{u,v\in \mathbb{Z}/l\mathbb{Z}}\phi \left(xy\left(v,u{\right)}_{-}\right)q\frac{\partial {r}_{v}}{\partial q}+4{\sum }_{u}\phi \left(xy\left(u,0{\right)}_{-}\right)q\frac{\partial {t}_{u}}{\partial q}.$  To handle the sum over $I\left(n\right)$  , for each $n$  we observe that $I\left(n\right)$  can be embedded into the disjoint union of two copies of $H\left(n\right)$  in two different ways as follows. The subset of $I\left(n\right)$  with ${m}_{1}\ge {m}_{2}$  can be identified with the subset of $H\left(n\right)$  with $c>0$  via $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)=\left(a,d-c,a-b,c\right)$  . The subset of $I\left(n\right)$  with ${m}_{1}<{m}_{2}$  can be identified with the subset of $H\left(n\right)$  with $bc>0$  via $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)=\left(a-b,c,a,d-c\right)$  . This describes the first embedding of $I\left(n\right)$  into the disjoint union of two copies of $H\left(n\right)$  . The second embedding is obtained by comparing ${k}_{i}$  . Namely, the subset of $I\left(n\right)$  with ${k}_{1}>{k}_{2}$  can be identified with the subset of $H\left(n\right)$  with $bc>0$  via $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)=\left(a-b,d,b,d-c\right)$  , and the subset of $I\left(n\right)$  with ${k}_{1}\le {k}_{2}$  can be identified with the subset of $H\left(n\right)$  with $b>0$  via $\left({m}_{1},{k}_{1},{m}_{2},{k}_{2}\right)=\left(b,d-c,a-b,d\right)$  . We will use these embeddings in order to rewrite the above sum over $I\left(n\right)$  in terms of $H\left(n\right)$  as follows. For the terms with $\left(-{k}_{2},{k}_{1}+{k}_{2}{\right)}_{-}$  we will use the first embedding, and for the terms with $\left(-{k}_{2},-{k}_{1}+{k}_{2}{\right)}_{-}$  we will use the second one. After some straightforward simplifications, we get $12\mu \circ PD\left(\phi \right)\sim 4{\sum }_{n>0}{q}^{n}\phi \left({\sum }_{H\left(n\right),bc>0}\left(\left(-ac{x}^{2}+\left(ad+bc\right)xy-bd{y}^{2}\right)\left(-c,d{\right)}_{-}+\left(\left(-ad-bc+ac+bd\right){x}^{2}+\left(ad+bc-2bd\right)xy+bd{y}^{2}\right)\left(c-d,d{\right)}_{-}+\left(\left(-ad-bc+ac+bd\right){x}^{2}+\left(ad+bc-2ac\right)xy+ac{y}^{2}\right)\left(c-d,-c{\right)}_{-}+\left(-bd{x}^{2}+\left(ad+bc\right)xy-ac{y}^{2}\right)\left(-d,c{\right)}_{-}\right)+{\sum }_{H\left(n\right),b=0,c>0}\left(-ac{x}^{2}+\left(ad+bc\right)xy-bd{y}^{2}\right)\left(-c,d{\right)}_{-}+{\sum }_{H\left(n\right),b>0,c=0}\left(\left(-bd\right){x}^{2}+\left(ad+bc\right)xy+\left(-ac\right){y}^{2}\right)\left(-d,c{\right)}_{-}+{\sum }_{d|n}\left(\frac{2{n}^{2}}{d}xy\left(0,d{\right)}_{-}+\frac{2nd}{l}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}xy\left(d,u{\right)}_{-}\right)\right)$  $=4{\sum }_{n>0}{q}^{n}\phi \left({\sum }_{H\left(n\right),bc>0}\left(\left(2ac{x}^{2}+2\left(ad+bc\right)xy+2bd{y}^{2}\right)\left(c,d{\right)}_{-}+\left(\left(-ad-bc+ac+bd\right){x}^{2}+\left(ad+bc-2bd\right)xy+bd{y}^{2}\right)\left(c-d,d{\right)}_{-}+\left(\left(-ad-bc+ac+bd\right){x}^{2}+\left(ad+bc-2ac\right)xy+ac{y}^{2}\right)\right)\left(c-d,-c{\right)}_{-}\right)+{\sum }_{H\left(n\right),bc=0}\left(ac{x}^{2}+\left(ad+bc\right)xy+bd{y}^{2}\right)\left(c,d{\right)}_{-}-{\sum }_{d|n}nxy\left(0,d{\right)}_{-}+{\sum }_{d|n}\left(\frac{2{n}^{2}}{d}xy\left(0,d{\right)}_{-}+\frac{2nd}{l}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}xy\left(d,u{\right)}_{-}\right)\right).$  We used various symmetries of $\left(u,v{\right)}_{-}$  to derive the last identity. A fortunate observation allows one to simplify the second and third lines of the last formula. Indeed, the relations on modular symbols imply $\left(\left(-ad-bc+ac+bd\right){x}^{2}+\left(ad+bc-2bd\right)xy+bd{y}^{2}\right)\left(c-d,d{\right)}_{-}$  $+\left(ac{x}^{2}+\left(ad+bc-2ac\right)xy+\left(-ad-bc+ac+bd\right){y}^{2}\right)\left(-c,c-d{\right)}_{-}$  $=\left(ac{x}^{2}+\left(ad+bc\right)xy+bd{y}^{2}\right)\left(c,d{\right)}_{-}.$  Then one gets $12\mu \circ PD\left(\phi \right)\sim 4{\sum }_{n>0}{q}^{n}\phi \left({\sum }_{H\left(n\right)}\left(3ac{x}^{2}+3\left(ad+bc\right)xy+3bd{y}^{2}\right)\left(c,d{\right)}_{-}-2{\sum }_{H\left(n\right),bc=0}\left(ac{x}^{2}+\left(ad+bc\right)xy+bd{y}^{2}\right)\left(c,d{\right)}_{-}-{\sum }_{d|n}nxy\left(0,d{\right)}_{-}+{\sum }_{d|n}\left(\frac{2{n}^{2}}{d}xy\left(0,d{\right)}_{-}+\frac{2nd}{l}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}xy\left(d,u{\right)}_{-}\right)\right)=4{\sum }_{n>0}{q}^{n}\phi \left({\sum }_{H\left(n\right)}\left(3ac{x}^{2}+3\left(ad+bc\right)xy+3bd{y}^{2}\right)\left(c,d{\right)}_{-}-2{\sum }_{d|n,d>c>0}\frac{n}{d}\left(c{x}^{2}+dxy\right)\left(c,d{\right)}_{-}-{\sum }_{d|n}nxy\left(0,d{\right)}_{-}+{\sum }_{d|n}\frac{2nd}{l}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}xy\left(u,d{\right)}_{-}\right).$  Using calculations similar to that of Section  4 , we can rewrite the last two lines in terms of fractional parts as $S=\phi \left(4{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}\left(\left(\frac{l}{d}\left(\left\{\frac{d-u}{l}\right\}-\left\{\frac{d-u}{l}{\right\}}^{2}\right)-\frac{l}{d}\left(\left\{-\frac{u}{l}\right\}-\left\{-\frac{u}{l}{\right\}}^{2}\right)+\frac{l-d}{l}-2\left\{\frac{u-d}{l}\right\}\right){x}^{2}+\left(-2\left\{\frac{u-d}{l}\right\}+1-\left\{-\frac{u}{l}\right\}+\left\{\frac{u}{l}\right\}\right)xy\right)\left(u,d{\right)}_{-}\right).$  We observe that for any $t$  there holds $\left\{t\right\}-\left\{t{\right\}}^{2}=\left\{-t\right\}-\left\{-t{\right\}}^{2}$  , and then use symmetries of $P\left(x,y\right)\left(±u,d{\right)}_{-}$  to see that $S\sim \phi \left(4{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}\left(-2\left\{\frac{u-d}{l}\right\}{x}^{2}+\left(-2\left\{\frac{u-d}{l}\right\}+1\right)xy\right)\left(u,d{\right)}_{-}\right)$  $=\phi \left(4{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}\left(\left(\left\{\frac{-u-d}{l}\right\}-\left\{\frac{u-d}{l}\right\}\right){x}^{2}+\left(-\left\{\frac{u-d}{l}\right\}-\left\{\frac{-u-d}{l}\right\}+1\right)xy\right)\left(u,d{\right)}_{-}\right)\sim \phi \left(4{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}{\sum }_{u\in \mathbb{Z}/l\mathbb{Z}}\left(\left(\left\{1-{\delta }_{d}^{-umodl}\right){x}^{2}+{\delta }_{d}^{umodl}xy\right)\left(u,d{\right)}_{-}\right)=\phi \left(4{\sum }_{n>0}n{q}^{n}{\sum }_{d|n}\left({x}^{2}+xy\right)\left(d,d{\right)}_{-}\right)\right)\sim 0.$  This finishes the proof.

6 Relation to Hecke eigenforms of rank zero

In this section we prove our main result that relates Wronskians of weight one Eisenstein series and the Hecke eigenforms of weight four with nonzero central value of $L$  -function.
Let ${T}_{n}$  denote the Hecke operators for ${\Gamma }_{1}\left(l\right)$  and let $L\left(f,s\right)$  denote the Hecke $L$  -function. We will normalize it so the central value is $L\left(f,2\right)={\int }_{0}^{i\infty }f\left(\tau \right)\tau d\tau$  . We say that a weight four Hecke eigenform $f$  has analytic rank zero if $L\left(f,2\right)\ne 0$  .
Definition 6.1. Let $f\in {\mathcal{S}}_{4}\left(l\right)$  be a weight four cusp form for ${\Gamma }_{1}\left(l\right)$  .
Define $\rho \left(f\right)={\sum }_{n>0}L\left({T}_{n}f,2\right){q}^{n}$  .
The following statements are analogous to the weight two calculation of [BG1.
Proposition 6.2. Definition  6.1 gives a linear map $\rho :{\mathcal{S}}_{4}\left(l\right)\to {\mathcal{S}}_{4}\left(l\right)$  , which commutes with ${\Gamma }_{0}\left(l\right)/{\Gamma }_{1}\left(l\right)$  -action. The image of $\rho$  contains all newforms $f$  with $L\left(f,2\right)\ne 0$  , and is contained in the span of all Atkin-Lehner lifts of all Hecke eigenforms $f$  of analytic rank zero.
Similar to [BG1, the key idea of this paper is to relate the map $\rho$  to the map $\mu$  of Section  4 .
Proposition 6.3. The map $\rho$  is the composition of the maps ${\mathcal{S}}_{4}\left(l\right)\stackrel{Int}{\to }\left({S}_{4}\left(l{\right)}_{-}{\right)}^{*}\stackrel{PD}{\to }{S}_{4}\left(l{\right)}_{+}\stackrel{\mu }{\to }{\mathcal{S}}_{4}\left(l\right)$  where $Int$  is induced by the integration pairing of ${\mathcal{S}}_{4}\left(l\right)$  and ${S}_{4}\left(l{\right)}_{-}$  , the $PD$  is the Poincaré duality map of Section  3 , and $\mu$  is the Wronskian map of Section  4 .
• Proof. We denote by $〈,〉$  the integration pairing between ${\mathcal{S}}_{4}\left(l\right)$  and ${S}_{4}\left(l{\right)}_{-}$  .
For a given $f\in {\mathcal{S}}_{4}\left(l\right)$  we calculate $\rho \left(f\right)={\sum }_{n>0}L\left({T}_{n}f,2\right){q}^{n}={\sum }_{n>0}〈{T}_{n}f,xy\left(0,1{\right)}_{-}〉{q}^{n}.$  By [M1,Theorem2andProposition10, $〈{T}_{n}f,xy\left(0,1{\right)}_{-}〉=〈f,{T}_{n}xy\left(0,1{\right)}_{-}〉=〈f,{\sum }_{H\left(n\right)}\left(ax+by\right)\left(cx+dy\right)\left(c,d{\right)}_{-}〉$  which leads to $\rho \left(f\right)={\sum }_{n>0}{q}^{n}〈f,{\sum }_{H\left(n\right)}\left(ax+by\right)\left(cx+dy\right)\left(c,d{\right)}_{-}〉$  Proposition  5.1 now shows $\mu \circ PD\circ Int\left(f\right)\sim \rho \left(f\right)$  and since both sides are quasimodular forms of weight four, the claim follows.
Corollary 6.4. The image of $\rho$  equals the linear span of $W\left({s}_{a},{s}_{b}\right)$  for $gcd\left(a,b,l\right)=1$  .
• Proof. Recall that $Int$  and $PD$  are isomorphisms, by [M1and Corollary  3.10 respectively. Then Proposition  4.6 finishes the proof.
We are now ready to formulate our main result.
Theorem 6.5. For arbitrary $l>1$  the span of Hecke eigenforms of weight four and analytic rank zero is equal to the span of the Wronskians $W\left({s}_{a},{s}_{b}\right)$  for all $a,b\in \mathbb{Z}/l\mathbb{Z}$  .
• Proof. In one direction, consider $f=W\left({s}_{a},{s}_{b}\right)$  . If $gcd\left(a,b,l\right)=d$  , then Corollary  6.4 applied to $\frac{l}{d}$  shows that $f$  is in $\rho \left({\mathcal{S}}_{4}\left(\frac{l}{d}\right)\right)$  . Indeed, $f$  is, up to a nonzero multiple, the $d$  -lift of $W\left({s}_{\frac{a}{d},\frac{l}{d}},{s}_{\frac{b}{d},\frac{l}{d}}\right)$  where the second subscript in $s$  is used to indicate the level. By Proposition  6.2 , $W\left({s}_{\frac{a}{d},\frac{l}{d}},{s}_{\frac{b}{d},\frac{l}{d}}\right)$  lies in the linear span of eigenforms of analytic rank zero, hence $f$  does as well.
To prove the opposite inclusion, it is enough to show that for any $d|l$  and any newform $g\left(\tau \right)\in {\mathcal{S}}_{4}\left(\frac{l}{d}\right)$  of analytic rank zero, its lift $g\left(k\tau \right)\in {\mathcal{S}}_{4}\left(l\right)$  lies in the span of $W\left({s}_{a},{s}_{b}\right)$  for any $k|d$  . By Proposition  6.2 , $g\in \rho \left({\mathcal{S}}_{4}\left(\frac{l}{d}\right)\right)$  . Then by Corollary  6.4 , $g$  is a linear combination of Wronskians of Eisenstein series ${s}_{i,\frac{l}{d}}$  of level $\frac{l}{d}$  . Then $g\left(k\tau \right)$  is a linear combination of Wronskians of $s$  -series of level $\frac{kl}{d}$  , since ${s}_{i,\frac{l}{d}}\left(k\tau \right)={s}_{ki,\frac{kl}{d}}\left(\tau \right)$  . Finally, $s$  -series of level $\frac{kl}{d}$  are sums of ${s}_{a}$  of level $l$  , which shows that $g\left(k\tau \right)$  lies in the span of the Wronskians $W\left({s}_{a},{s}_{b}\right)$  , as claimed.
Corollary 6.6. The span of Hecke eigenforms of weight four and analytic rank zero for the group ${\Gamma }_{0}\left(l\right)$  coincides with the span of ${\sum }_{j\in \left(\mathbb{Z}/l\mathbb{Z}{\right)}^{*}}W\left({s}_{aj},{s}_{bj}\right)$  for all $a,b\in \mathbb{Z}/l\mathbb{Z}$  .
• Proof. Use the formulas for the action of ${\Gamma }_{0}\left(l\right)$  on ${s}_{a}$  from [BG2.
References

1. L. Borisov, P. Gunnells, Toric modular forms and nonvanishing of $L$  -functions, J. Reine Angew. Math. 539 (2001), 149–165.
2. L. Borisov, P. Gunnells, Toric modular forms of higher weight, J. Reine Angew. Math. 560 (2003), 43–64.
3. W. Kohnen, D. Zagier, Modular forms with rational periods. Modular forms (Durham, 1983), 197–249, Ellis Horwood Ser. Math. Appl.: Statist. Oper. Res., Horwood, Chichester, 1984.
4. L. Merel, Universal Fourier expansions of modular forms. On Artin's conjecture for odd $2$  -dimensional representations, 59–94, Lecture Notes in Math., 1585, Springer, Berlin, 1994.
5. L. Merel, Modular forms and $L$  -functions, lectures in Taiwan, August 2001.
6. A. Scholl, An introduction to Kato's Euler systems. Galois representations in arithmetic algebraic geometry (Durham, 1996), 379–460, London Math. Soc. Lecture Note Ser., 254, Cambridge Univ. Press, Cambridge, 1998.
7. V. Shokurov, Shimura integrals of cusp forms Izv. Akad. Nauk SSSR Ser. Mat. 44 (1980), no. 3, 670–718.
8. K. Haberland, Periods of modular forms of one variable and group cohomology. I, II, III, Math. Nachr. 112 (1983), 245–315.
9. D. Zagier, Modular forms of one variable, notes of a course given in Utrecht, Spring 1991.

Department of Mathematics, University of Wisconsin, 417 Van Vleck Hall, 480 Lincoln Dr, Madison, WI 53706, USA. email:borisov@math.wisc.edu.