---
abstract: 'In this paper, we prove a conjecture of Broadhurst and Zudilin concerning a divisibility property of the Fourier coefficients of a meromorphic modular form using the generalization of the Shimura lift by Borcherds and Hecke operators on vector-valued modular forms developed by Bruinier and Stein. Furthermore, we construct a family of meromorphic modular forms with this property.'
address: 'Fachbereich Mathematik, Technische Universität Darmstadt, Schlossgartenstrasse 7, D–64289 Darmstadt, Germany'
author:
- Yingkun Li and Michael Neururer
bibliography:
- 'magnetic.bib'
title: A Magnetic Modular Form
---

[^1]

Introduction.
=============

Let $f$ be a weakly holomorphic modular form with the Fourier expansion $\sum_{n \gg -\infty} a(n) q^n_z$, where $q_z := {{\mathbf{e}}}(z) := e^{2\pi i z}$ for $z$ in the upper-half complex plane ${{\mathcal{H}}}$. Suppose $f$ has weight 2, level 1 and integral Fourier coefficients. Then it possesses the following divisibility property $$\label{eq:divprop}
  n \mid a(n) \text{ for all } n \in {\mathbb{N}}.$$ Indeed, it is the image of a polynomial in the Klein $j$-invariant with integral coefficients under the differential operator $q\frac{d}{dq}$. More generally, one can apply this derivative $k-1$ times to weakly holomorphic modular forms of weight $2-k$ to produce forms of weight $k$ with this property (see e.g. [@KazalickiScholl]). Analogues in higher weights have been observed in [@DukeJenkins08] as a consequence of the action of Hecke operators on weakly holomorphic modular forms. The same phenomenon in the half-integral weight case has also been studied [@Edixhoven05; @Jenkins05; @ZhangYC16].

In [@Ausserlechner2016], Ausserlechner studied how the output voltage of a Hall plate is affected by the shape of the plates and sizes of the contacts, where he encountered a double integral generalizing the classical elliptic integral used to evaluate the arithmetic geometric mean. $$\label{eq:I2}
  I_2(f) := (1 + f) \int \int_{\pi/2 > \alpha > \beta > 0} \frac{d\alpha d\beta}{\sqrt{(( 1 + f)^2 - 4 f \cos^2 \alpha)(( 1 + f)^2 - 4f \sin^2 \beta)}}.$$In [@BZ17], Broadhurst and Zudilin studied this integral $I_2(f)$ in detail and showed that it satisfies $$I_2(f) = I_2 \left( \frac{1-f}{1+f} \right)\text{ for } f \in [0, 1],$$ which was conjectured by Ausserlechner and also proved by Glasser and Zhou [@GlasserZhou18]. After applying a modular parametrization, they also showed that $I_2(f)$ satisfies an inhomogeneous differential equation, whose constant term is the following modular form $$\label{eq:phi}
  \begin{split}
  {\phi}(z) &:=  (\eta(2z) \eta(4z))^4 \frac{1 - 96 \psi(2z) + 256\psi(2z)^2}{(1 + 16\psi(2z))^2} \\
&= \sum_{n\in 2{\mathbb{N}}-1} a(n) q^{n}_z = q_z - 132q_z^{3} + 5630q_z^{5} - 189672q_z^{7} + 5768181q_z^{9} + O(q^{10})
.    
  \end{split}
$$ Here, $$\label{eq:psi}
\psi(z) := \frac{\eta(z)^8 \eta(4z)^{16} }{\eta(2z)^{24}}$$ is a Hauptmodul on the modular curve $X_0(4)$. The function ${\phi}$ is a meromorphic modular form of weight 4 on $\Gamma_0(8)$ and has double poles at $z_\pm := \frac{\pm 1 + i}{4}$. The numbers $a(n)$ are a priori integers. From numerical computations, Broadhurst and Zudilin made the following conjecture.

\[conj:main\] The meromorphic modular form $\phi(z)$ satisfies .

On the one hand, the conjecture above seems like a natural and expected extension of for weakly holomorphic forms. On the other hand, $p^3 \nmid a(p)$ for every odd prime $p$ (see ). Therefore, unlike elements of the bases in [@DukeJenkins08], $\phi$ does not come from applying the operator $(q \tfrac{d}{dq})^3$ to a modular form of weight $-2$, which makes this conjecture surprising.

In this note, we will prove this conjecture by first realizing ${\phi}(z)$ as the regularized theta lift of a half-integral weight, vector-valued modular form. In [@Shimura73], Shimura initiated the study of modular forms of half-integral weight, and showed they correspond to modular forms of integral weight. This is called the Shimura lift. Using techniques introduced by Shintani [@Shintani75], Niwa expressed the Shimura lift of a holomorphic, half-integral weight modular form as its integral against a suitable theta kernel [@Niwa75]. In [@Borcherds98], Borcherds expanded the input space to include vector-valued modular forms on the metaplectic cover ${{\mathrm{Mp}}}_2({\mathbb{Z}})$ of ${{\mathrm{SL}}}_2({\mathbb{Z}})$ with singularities at the cusps. This will be the setting that we work in.

To describe the vector-valued input, we need the following 3-dimensional representation $\varrho: {{\mathrm{Mp}}}_2({\mathbb{Z}}) \to {{\mathrm{GL}}}_3({\mathbb{C}})$ $$\label{eq:varrho}
  \varrho(T) := 
\begin{pmatrix}
\zeta_8  & & \\ & \zeta_8^7 & \\ && \zeta_8^5 
\end{pmatrix}, ~
  \varrho(S) := \frac{\overline{\zeta_8}}{2} \begin{pmatrix}
1 & 2 & 1 \\ {1} & 0 & -{1} \\ 1 & -2 & 1
\end{pmatrix},$$ where $\zeta_8 := e^{2\pi i /8 }$ is an $8$[^th^ ]{}root of unity, and $T, S$ are the generators of ${{\mathrm{Mp}}}_2({\mathbb{Z}})$ (see ). It is isomorphic to a subrepresentation of the 64 dimensional Weil representation $\rho_L$ associated to a particular lattice $L$ (see section \[subsec:lattice\]). Let $\{{\mathfrak{e}}_\ell: \ell = 0, 1, 2\}$ be the standard basis of ${\mathbb{C}}^3$. Note that if we scale ${\mathfrak{e}}_1$ by $\sqrt{2}$, then the representation $\varrho$ becomes unitary with respect to the standard inner product on ${\mathbb{C}}^3$. If $q:= {{\mathbf{e}}}(\tau)$ with $\tau \in {{\mathcal{H}}}$, then the action of $\varrho(T)$ implies that any modular form $G$ in $M^!_{5/2, \varrho}$ (see section \[subsec:vvmf\] for notation) has a Fourier expansion in $q^{1/8}$, and the coefficient of $q^{m/8}$ is zero if $m \not\equiv 1, 5$ or $7 \bmod{8}$. For any $m \in {\mathbb{Z}}$ congruent to $1, 7$ or $5$ modulo $8$, let $\ell(m) \in \{0, 1, 2\}$ be the unique index satisfying $1 - m \equiv 2\ell(m) \bmod{8}$. For any modular form $G \in M^!_{5/2, \varrho}$, we can write out its Fourier expansion as $$\label{eq:GFE}
G(\tau) = \sum_{m \in {\mathbb{Z}}} c(G,  m )\mathfrak{q}^m,~ \mathfrak{q}^m :=  q^{m/8} {\mathfrak{e}}_{\ell(m)}$$ and define the formal power series $$\label{eq:Psi}
  {{\Phi}}(z, G) := \sum_{n \in 2{\mathbb{N}}- 1} (-1)^{(n-1)/2} q^n_z \sum_{r \mid n} r \cdot c(G, (n/r)^2),$$ which converges absolutely for $y := {\mathrm{Im}}(z)$ sufficiently large, and analytically continues to a meromorphic function in $z \in {{\mathcal{H}}}$, which we also denote by ${{\Phi}}(z, G)$. Before plunging into heavier notations, we will first state a result in terms of scalar-valued modular forms. For a weakly holomorphic modular form $F \in M^!_1(\Gamma_1(4))$ with the Fourier expansion $F(\tau)\eta(\tau)^3 = \sum_{n \gg -\infty} c(n) q^{n/8}$, we can define a meromorphic function on ${{\mathcal{H}}}$ by Our first result is as follows.

\[thm:ShimLift\] In the notation above, the function ${{\Phi}}(z, G)$ is in $M^\mathrm{mero}_{4, \chi}(\Gamma^+_0(8))$, where $\Gamma_0^+(8)\subset {{\mathrm{SL}}}_2({\mathbb{R}})$ is the extension of $\Gamma_0(8)$ by the matrices $R, U, W$ defined in and $\chi$ is the character in In particular, there exists a unique $G_1
\in M^!_{5/2, \varrho}$ with $ -4 {\phi}(z) =  {{\Phi}}(z, G_1)$ (see ).

\[rmk:vanish\] Using SageMath [@sagemath], it is easy to check that $M_4(\Gamma_0(8))$ is 5-dimensional, with a $W$-eigenbasis consisting of 4 Eisenstein series and 1 cusp form. The Eisenstein series are $E_4(z) \pm 64 E_4(8z)$ and $E_4(2z) \pm 4 E_4(4z)$ with $E_4 \in M_4$ the Eisenstein series of weight 4. Among the eigenbasis, only $E_4(z) - 64 E_4(8z)$ and $E_4(2z) - 4E_4(4z)$ have $W$-eigenvalue $-1$. From their Fourier expansions, it is clear that the only linear combination with $R$-eigenvalue $-1$ is $0$. Therefore, $M_{4, \chi}(\Gamma^+_0(8))$ is trivial.

$$\label{eq:G11}
  G_{1, 1}(\tau)  = 2^7 \frac{\eta(2\tau)^{10} \eta(4\tau)^4}{\eta(\tau)^{9}} + 4 \eta(\tau)^3 E_4(\tau) \in M^!_1(\Gamma_1(4)),$$ where $E(\tau, \chi_4) := \tfrac{1}{4} + \sum_{n \ge 1} q^n \sum_{r \mid n} \chi_4(r) \in M_1(\Gamma_1(4))$ is the Eisenstein series with $\chi_4(r) := (\tfrac{-4}{r})$, such that

To prove the conjecture, it is necessary to study the integral structure of $M^!_{5/2, \varrho}$. In section \[subsec:orthogp\], we will see that $\varrho$ can be realized as a rational subrepresentation of a Weil representation $\rho_L$ attached to the lattice $L$ in section \[subsec:lattice\]. Therefore, the ${\mathbb{Z}}$-module ${\mathbb{M}}^!_{5/2, \varrho}$ of modular forms in $M^!_{5/2, \varrho}$ with integral Fourier coefficients is free and a complete lattice in $M^!_{5/2, \varrho}$ by a theorem of McGraw [@McGraw03] (see section \[subsec:vvmf\]). Through explicit construction in section \[subsec:family\] and the theory of Hecke operators on vector-valued modular forms developed in [@BS10], we will prove the following result.

\[thm:HeckeSc\] The free ${\mathbb{Z}}$-module ${\mathbb{M}}^!_{5/2, \varrho}$ has a canonical basis $
\{G_{d}: d \in {\mathbb{N}}, d \equiv 1, 3, 7 \bmod{8}\}
$ characterized by the property $$\label{eq:polecond}
  G_{d}(\tau) =  q^{-d/8} {\mathfrak{e}}_{\ell(-d)} + O(q^{1/8}).$$ Suppose $G_d$ has Fourier coefficients $c(G_{d}, m) \in {\mathbb{Z}}$ as in . Then $$\label{eq:thm2}
  c(G_d, m^2) \in m {\mathbb{Z}}$$ for all $m \in {\mathbb{N}}$ and square-free $d \in {\mathbb{N}}$.

From the definition of the map ${\Phi}$, it is clear that Theorems \[thm:ShimLift\] and \[thm:HeckeSc\] together imply Conjecture \[conj:main\], i.e., $n|a(n)$. In fact the Hecke theory allows us to study the coefficients $a(n)$ modulo $n^3$. For instance Corollary \[cor:p\^3 congruence\] implies that $$\label{eq:pdiv}
  p^3 \mid (a(p)-p)$$ for every odd prime $p$. Furthermore, this puts ${\phi}_1(z) = -4 {{\phi}(z)}$ into a family of modular forms $\{{\phi}_d(z) := {\Phi}(z, G_d): d \equiv 1, 3, 7 \bmod{8} \text{ square-free}\}$ in ${\mathbb{M}}^\mathrm{mero}_{4, \chi}(\Gamma^+_0(8))$ that all satisfy the divisibility property . It is worth mentioning that the phenomenon of Zagier duality [@Za98; @ZhangYC16] is also present between this basis and the canonical basis of the free ${\mathbb{Z}}$-module ${\mathbb{M}}^!_{-1/2, {\varrho^*}}$, where $\varrho^*$ is the unitary dual of $\varrho$ with respect to the standard inner product on ${\mathbb{C}}^3$ (see Prop. \[prop:family\]). Finally, we can apply the same idea to Example 14.4 in [@Borcherds98] to prove the following result.

\[thm:scalar\] Let $\Delta(z)$ be Ramanujan’s delta function, and $E_4(z) = 1 + 240q_z + O(q_z^2)$ the Eisenstein series of weight 4. Then the meromorphic modular form $64\frac{\Delta(z)}{E_4(z)^2}$ satisfies .

In fact numerically, also $\Delta(z)/E_4(z)^2$ satisfies . It is ongoing work of the authors to remove the factor $64$ from the above theorem.

Can a non-constant holomorphic modular form have the divisibility property ? It seems likely that even the weaker condition, that every sufficiently large prime divides the corresponding Fourier coefficient, cannot be satisfied by a holomorphic modular form.

This note is organized as follows. In section \[sec:prelim\], we give the preliminaries concerning vector-valued modular forms as input to Borcherds’ lift and prove Theorem \[thm:ShimLift\] as a special case of Borcherds’ result. In section \[sec:3\], we construct the family $\{G_d\}$, prove Conjecture \[conj:main\], and give some numerical data of the Fourier expansions of the bases.

**Acknowledgement.** We thank Paul Jenkins for helpful conversations.

Preliminaries. {#sec:prelim}
==============

In this section we will introduce the vector-valued theta lift à la Borcherds [@Borcherds98].

Vector-Valued Modular Forms. {#subsec:vvmf}
----------------------------

Denote ${{\mathcal{H}}}^* := {{\mathcal{H}}}\cap \mathbb{P}^1({\mathbb{Q}})$ the extended upper half plane, which is acted on by ${{\mathrm{SL}}}_2({\mathbb{R}})$ via linear fractional transformation. Let ${{\mathrm{Mp}}}_2({\mathbb{R}})$ be the metaplectic cover of ${{\mathrm{SL}}}_2({\mathbb{R}})$ consisting of pairs $(\gamma,\varphi)$, where $\gamma = {\left(\begin{smallmatrix}
                 a & b\\
                 c & d
\end{smallmatrix}\right)}\in{{\mathrm{SL}}}_2({\mathbb{R}})$ and $\varphi$ is a holomorphic function on ${{\mathcal{H}}}$ satisfying $
\varphi(\tau)^2 = (c\tau+d).
$ We denote the preimage of ${{\mathrm{SL}}}_2({\mathbb{Z}})$ in ${{\mathrm{Mp}}}_2({\mathbb{R}})$ under the covering map by ${{\mathrm{Mp}}}_2({\mathbb{Z}})$, which is generated by $$\label{eq:ST}
  T:=\left({\begin{pmatrix}
                 1 & 1\\
                 0 & 1
\end{pmatrix}},1\right),~
  S:=\left({\begin{pmatrix}
                 0 & -1\\
                 1 & 0
\end{pmatrix}},\sqrt{\tau}\right).$$ Here we take the principal branch of the square root.

Let $\Gamma \subset {{\mathrm{Mp}}}_2({\mathbb{R}})$ be a subgroup commensurable with ${{\mathrm{Mp}}}_2({\mathbb{Z}})$. A meromorphic function $f : {{\mathcal{H}}}^* \to {\mathbb{C}}$ is called a meromorphic modular form of weight $k \in \tfrac{1}{2} {\mathbb{Z}}$ with respect to a unitary representation $\rho: \Gamma \to {{\mathrm{GL}}}(W)$ on a finite dimensional ${\mathbb{C}}$-vector space $W$ if it satisfies $$\label{eq:modular}
  (f \mid_k (\gamma, \varphi))(\tau) := \varphi(\tau)^{-2k} f (\gamma \cdot \tau) = \rho((\gamma, \varphi)) f$$ for all $(\gamma, \varphi) \in \Gamma$. We use $M^{\mathrm{mero}}_{k, \rho}(\Gamma)$ to denote the vector space of these meromorphic modular forms. It contains the subspaces $M^!_{k, \rho}(\Gamma), M_{k, \rho}(\Gamma), S_{k, \rho}(\Gamma)$ of weakly holomorphic, holomorphic, and cuspidal modular forms.

When $k \in {\mathbb{Z}}$, it is necessary for $\rho$ to factor through the image of $\Gamma$ in ${{\mathrm{SL}}}_2({\mathbb{Z}})$, and we replace $\Gamma$ with its image in ${{\mathrm{SL}}}_2({\mathbb{Z}})$ in the above notation. Also, we drop the subscript $\rho$, resp. $\Gamma$, if it is 1-dimensional and trivial, resp. ${{\mathrm{Mp}}}_2({\mathbb{Z}})$. For $N \in {\mathbb{N}}$, we denote by $\Gamma^*_0(N)$ the extension of the congruence subgroup $\Gamma_0(N)$ by Atkin-Lehner operators [@AtkinLehner] A common type of representation comes from arithmetic. Let $L$ be an even, integral lattice with quadratic form $Q$. Let $(b^+,b^-)$ be the signature of $L_{\mathbb{R}}:= L \otimes_{\mathbb{Z}}{\mathbb{R}}$ and $L^\vee := \mathrm{Hom}_{{\mathbb{Z}}}(L, {\mathbb{Z}})$ be the dual lattice. Through the bilinear form $(\cdot, \cdot)$ induced by $Q$ on $L$, we identify $L^\vee$ with a sublattice of $L_{\mathbb{R}}$ containing $L$. The quotient ${{\mathcal{A}}}_L:= L^\vee/L$ is then a finite abelian group, on which $Q$ becomes a quadratic form valued in ${\mathbb{Q}}/{\mathbb{Z}}$. On the vector space ${\mathbb{C}}[{{\mathcal{A}}}_L]$ with basis $\{{\mathfrak{e}}_h: h \in {{\mathcal{A}}}_L\}$, there is a hermitian inner product $\langle\cdot, \cdot \rangle$ defined by $$\label{eq:lrangle}
  \langle v, w \rangle := \sum_{h \in {{\mathcal{A}}}_L} v_h \overline{w_h}$$ for $v = \sum_{h \in {{\mathcal{A}}}} v_h {\mathfrak{e}}_h$ and $w = \sum_{h \in {{\mathcal{A}}}} w_h {\mathfrak{e}}_h$ in ${\mathbb{C}}[{{\mathcal{A}}}_L]$, which induces the norm $$\label{eq:norm}
 \| v \|:= \sqrt{\langle v, v \rangle }$$ on ${\mathbb{C}}[{{\mathcal{A}}}]$. The group ${{\mathrm{Mp}}}_2({\mathbb{Z}})$ acts through the Weil representation $\rho_L$ defined by $$\label{eq:weil}
\rho_L(T){\mathfrak{e}}_h={{\mathbf{e}}}(Q(h)){\mathfrak{e}}_h,~\rho_L(S){\mathfrak{e}}_h=
\frac{{{\mathbf{e}}}((b^+-b^-)/{8})}{|{{\mathcal{A}}}_L|}\sum_{\ell\in {{\mathcal{A}}}_L} {{\mathbf{e}}}(-(h,\ell)){\mathfrak{e}}_\ell,$$ which is unitary with respect to the hermitian inner product in . By a theorem of McGraw [@McGraw03], the space $M^!_{k, \rho_L}$ has a basis with Fourier coefficients in ${\mathbb{Q}}$. We will use ${\mathbb{M}}^!_{k, \rho_L}$ to denote the ${\mathbb{Z}}$-module of modular forms in $M^!_{k, \rho_L}$ with integral Fourier coefficients.

One can define the orthogonal group of ${{\mathcal{A}}}_L$ as the following finite group $$\label{eq:OAL}
   \mathrm{O}({{\mathcal{A}}}_L) := \{ \sigma: {{\mathcal{A}}}_L \to {{\mathcal{A}}}_L \text{ group automorphism } \mid Q(\sigma(h)) = Q(h) \text{ for all } h \in {{\mathcal{A}}}_L\}.$$ Every element in $\mathrm{O}({{\mathcal{A}}}_L)$ induces a $\rho$-linear automorphism on ${\mathbb{C}}[{{\mathcal{A}}}_L]$, hence also acts on $M^{\mathrm{mero}}_{k, \rho_L}$, which decomposes according to the irreducible representations of $\mathrm{O}({{\mathcal{A}}}_L)$. For each $h \in {{\mathcal{A}}}_L$, we have the normal subgroup $\mathrm{O}({{\mathcal{A}}}_L)_h \subset \mathrm{O}({{\mathcal{A}}}_L)$ consisting of the stabilizers of $h$ in $\mathrm{O}({{\mathcal{A}}}_L)$.

Symmetric Space.
----------------

Let $V_{\mathbb{R}}:= M_2({\mathbb{R}})^0$ be the real vector space of 2 by 2 matrices with trace 0. It becomes a real quadratic space of signature $(2, 1)$ with respect to the quadratic form $Q := -N \cdot \det$ for any natural number $N$. The group ${{\mathrm{GL}}}_2({\mathbb{R}})$ acts isometrically on $V_{\mathbb{R}}$ via conjugation, which is explicitly given by $$\label{eq:SOaction}
\begin{split}
  \gamma \cdot {\begin{pmatrix}
                 B & C\\
                 -A & -B
\end{pmatrix}} &= \gamma {\begin{pmatrix}
                 B & C\\
                 -A & -B
\end{pmatrix}} \gamma^{-1} \\
&= {\begin{pmatrix}
                 ad B - bdA - acC + bcB & -2abB + b^2 A + a^2C\\
                 2cdB - d^2 A - c^2C & -bcB - adB + bdA + acC
\end{pmatrix}}
\end{split}$$ for $\gamma = {\left(\begin{smallmatrix}
                 a & b\\
                 c & d
\end{smallmatrix}\right)} \in {{\mathrm{GL}}}_2({\mathbb{R}})$. This identifies ${{\mathrm{SL}}}_2({\mathbb{R}})$ with ${{\mathrm{Spin}}}(V_{\mathbb{R}})$ and ${{\mathrm{PSL}}}_2({\mathbb{R}})$ with ${{\mathrm{SO}}}^+(V_{\mathbb{R}})$, the connected component of the special orthogonal group ${{\mathrm{SO}}}(V_{\mathbb{R}})$ containing the identity. For any $\gamma \in {{\mathrm{SL}}}_2({\mathbb{R}})$, we also use $\gamma$ to represent its image in ${{\mathrm{SO}}}^+(V_{\mathbb{R}})$.

Let ${{\mathcal{D}}}$ be the symmetric space of oriented negative lines in $V_{\mathbb{R}}$ and ${{\mathcal{D}}}^0 \subset {{\mathcal{D}}}$ the connected component containing ${\mathbb{R}}\cdot {\begin{pmatrix}
                 0 & 1\\
                 -1 & 0
\end{pmatrix}}$. As usual, we can use ${{\mathcal{H}}}$ to parametrize ${{\mathcal{D}}}^0$ by $$\label{eq:parametrize}
  \begin{split}
    {{\mathcal{H}}}& \to {{\mathcal{D}}}^0 \\
z &\mapsto {\mathbb{R}}\cdot {\begin{pmatrix}
                 -x & |z|^2\\
                 -1 & x
\end{pmatrix}}.
  \end{split}$$ As usual, we can use ${{\mathcal{H}}}$ to parametrize ${{\mathcal{D}}}^0$ by defining $$\label{eq:parametrization}
  Z(z) :=  {\begin{pmatrix}
                 -z & z^2\\
                 -1 & z
\end{pmatrix}}
$$ for each $z \in {{\mathcal{H}}}$. Then $\{{\mathrm{Re}}(Z(z)), {\mathrm{Im}}(Z(z))\}$ always span a positive definite 2-plane in $V_{\mathbb{R}}$ and its orthogonal complement is an element of ${{\mathcal{D}}}^0$. Furthermore, $$\label{eq:invariance}
  \gamma \cdot Z(z) = (cz + d)^{-2} Z(\gamma z)$$ for all $\gamma = {\left(\begin{smallmatrix}
                 a & b\\
                 c & d
\end{smallmatrix}\right)} \in {{\mathrm{GL}}}^+_2({\mathbb{R}})$.

Lattice. {#subsec:lattice}
--------

Let $L \subset M_2({\mathbb{Q}})$ be the following lattice $$\label{eq:L}
  L := 
\left\{
{\begin{pmatrix}
                 B & C/2\\
                 -4A & -B
\end{pmatrix}}: A, B, C \in {\mathbb{Z}}\right\},$$ which is even integral with respect to the quadratic form $Q = -2\cdot {\det}{}$ with dual lattice $$\label{eq:Lvee}
L^\vee := 
\left\{
{\begin{pmatrix}
                 b/4 & c/8\\
                 -a & -b/4
\end{pmatrix}}: a, b, c \in {\mathbb{Z}}\right\}.$$ The real quadratic space $L_{\mathbb{R}}$ has signature $(2, 1)$. It is not hard to see that $L$ is the direct sum of the following two sublattices $L_1$ and $L_2$ of signatures $(1,1)$ and $(1,0)$: $$\label{eq:L12}
  L_1 := 
\left\{
  {\begin{pmatrix}
                 0 & C/2\\
                 -4A & 0
\end{pmatrix}} \in L
\right\}, \quad
L_2 := \left\{
{\begin{pmatrix}
                 B & 0\\
                 0 & -B
\end{pmatrix}} \in L
\right\}.$$ Furthermore, the dual lattice $L_j^\vee \subset L_{j, {\mathbb{R}}} \subset L_{\mathbb{R}}$ is contained in $L^\vee$ for $j = 1, 2$.

The finite quadratic module ${{\mathcal{A}}}:= L^\vee/L$ is isometric to $({\mathbb{Z}}/4{\mathbb{Z}})^3$ via the map $$\label{eq:Z4isom}
  \begin{split}
  {{\mathcal{A}}}&\stackrel{\cong}{\to} ({\mathbb{Z}}/4{\mathbb{Z}})^3 \\
{\begin{pmatrix}
                 b/4 & c/8\\
                 -a & -b/4
\end{pmatrix}} & \mapsto (a, b, c),
\end{split}$$ where the quadratic form on $({\mathbb{Z}}/4{\mathbb{Z}})^3$ is $Q((a, b, c)) := \frac{b^2 - 2ac}{8} \in {\mathbb{Q}}/{\mathbb{Z}}$. The isotropic elements are $$\label{eq:isoA}
{{\mathrm{Iso}}}({{\mathcal{A}}}) = \left\{
\begin{aligned}
&( 0 , 0 , 0 ) ,
( 0 , 0 , 1 ) ,
( 0 , 0 , 2 ) ,
( 0 , 0 , 3 ) ,
( 1 , 0 , 0 ) ,
( 1 , 2 , 2 ) ,\\
&( 2 , 0 , 0 ) ,
( 2 , 0 , 2 ) ,
( 2 , 2 , 1 ) ,
( 2 , 2 , 3 ) ,
( 3 , 0 , 0 ) ,
( 3 , 2 , 2 )
\end{aligned}
\right\}.$$ For $j = 1, 2$, let ${{\mathcal{A}}}_j$ be the finite quadratic module $L_j^\vee/L_j$ and $\rho, \rho_j$ be the Weil representations associated to $L, L_j$ respectively. Then ${{\mathcal{A}}}_j \subset {{\mathcal{A}}}$ and ${{\mathcal{A}}}= {{\mathcal{A}}}_1 \oplus {{\mathcal{A}}}_2$, which induces $$\label{eq:rhoequal}
 {\mathbb{C}}[{{\mathcal{A}}}] = {\mathbb{C}}[{{\mathcal{A}}}_1] \otimes {\mathbb{C}}[{{\mathcal{A}}}_2],~\rho = \rho_1 \otimes \rho_2.$$ We also identify ${{\mathcal{A}}}_1$, resp. ${{\mathcal{A}}}_2$, with $({\mathbb{Z}}/4{\mathbb{Z}})^2$, resp. ${\mathbb{Z}}/4{\mathbb{Z}}$, by sending $(a, 0, c)$, resp. $(0, b, 0)$, to $(a, c) \in ({\mathbb{Z}}/4{\mathbb{Z}})^2$, resp. $b \in {\mathbb{Z}}/4{\mathbb{Z}}$.

Let ${{\mathrm{SO}}}(L) \subset {{\mathrm{SO}}}(V)$ be the special orthogonal group stabilizing the lattice. It also fixes the dual lattice, hence induces an action on the finite quadratic module ${{\mathcal{A}}}$. Denote $\Gamma_L \subset {{\mathrm{SO}}}(L)$ the kernel of this action. In fact, $\Gamma_L$ is contained in ${{\mathrm{SO}}}^+(L) := {{\mathrm{SO}}}(L) \cap {{\mathrm{SO}}}^+(L_{\mathbb{R}})$ After scaling by $\sqrt{2}$, the lattice $L$ is the same as the lattice $L(8, 4)$ in Section 4 of [@Zemel17a], where ${{\mathrm{SO}}}^+(L)$ and $\Gamma_L$ were given explicitly in Theorem 4.2 loc. cit. as $$\label{eq:SO}
  {{\mathrm{SO}}}^+(L) := {\begin{pmatrix}
                 1/2 & 0\\
                 0 & 1
\end{pmatrix}} \Gamma^*_0(2) {\begin{pmatrix}
                 2 & 0\\
                 0 & 1
\end{pmatrix}}, ~ \Gamma_L := \Gamma_0(8).$$ Note that $\Gamma^*_0(2) \subset {{\mathrm{SL}}}_2({\mathbb{R}})$ is obtained from $\Gamma_0(2)$ by adjoining the Fricke-involution $W_2 := {\left(\begin{smallmatrix}
                 0 & 1/\sqrt{2}\\
                 -\sqrt{2} & 0
\end{smallmatrix}\right)}$. Consider the following elements in ${{\mathrm{SL}}}_2({\mathbb{R}})$ $$\label{eq:gens}
R:= {\begin{pmatrix}
                 1 & 1/2\\
                 0 & 1
\end{pmatrix}}, ~ U:= {\begin{pmatrix}
                 1 & 0\\
                 4 & 1
\end{pmatrix}}, ~ W := {\begin{pmatrix}
                 0 & 1/(2\sqrt{2})\\
                 -2\sqrt{2} & 0
\end{pmatrix}}.$$ Then their images in ${{\mathrm{PSL}}}_2({\mathbb{R}})$, along with $\Gamma_L$, generate ${{\mathrm{SO}}}^+(L)$, which is called the group of Atkin-Lehner operators. This is a reasonable name since if $L$ is the lattice studied in [@BO10], then $\Gamma_L$, resp. ${{\mathrm{SO}}}^+(L)$, is isomorphic to the image of the congruence subgroup $\Gamma_0(N) \subset {{\mathrm{SL}}}_2({\mathbb{Z}})$, resp. $\Gamma_0^*(N)$, in ${{\mathrm{PSL}}}_2({\mathbb{R}})$. We denote the preimage of ${{\mathrm{SO}}}^+(L)$ in ${{\mathrm{SL}}}_2({\mathbb{R}})$ by $\Gamma_0^+(8)$, which contains $\Gamma_0(8)$.

Orthogonal Groups. {#subsec:orthogp}
------------------

It is important to understand the finite group $\mathrm{O}(A)$, since ${\mathbb{C}}[{{\mathcal{A}}}]$ decomposes according to its irreducible representations. By considering ${{\mathcal{A}}}$ as a free ${\mathbb{Z}}/4{\mathbb{Z}}$-module of rank 3 through the map , we see that the group $\mathrm{O}({{\mathcal{A}}})$ defined in has the following form $$\label{eq:SOAc}
  \mathrm{O}({{\mathcal{A}}}) = \left\{g \in {{\mathrm{GL}}}_3({\mathbb{Z}}/4{\mathbb{Z}}): Q(g\cdot h) = Q(h) \text{ for all } h \in ({\mathbb{Z}}/4{\mathbb{Z}})^3 \cong {{\mathcal{A}}}\right\}.$$ The same holds for the free ${\mathbb{Z}}/4{\mathbb{Z}}$-modules ${{\mathcal{A}}}_1, {{\mathcal{A}}}_2$, and we can then define ${{\mathrm{SO}}}({{\mathcal{A}}})$, resp. ${{\mathrm{SO}}}({{\mathcal{A}}}_j)$, as the subgroup of $\mathrm{O}({{\mathcal{A}}})$, resp. $\mathrm{O}({{\mathcal{A}}}_j)$, consisting of the elements with determinant $1 \bmod 4$. Let $\nu$ be the negative identity matrix, which is contained in all three orthogonal groups above. Note that for the finite quadratic modules ${{\mathcal{A}}}, {{\mathcal{A}}}_1$ and ${{\mathcal{A}}}_2$ above, we always have the following isometry $$\label{eq:nu}
  \nu(h) := -h$$ in the orthogonal group. It is even in ${{\mathrm{SO}}}({{\mathcal{A}}}_1)$ since ${{\mathcal{A}}}_1$ has even rank.

There is a natural map ${{\mathrm{SO}}}(L)/\Gamma_L \hookrightarrow {{\mathrm{SO}}}({{\mathcal{A}}})$. By checking every element in the finite group ${{\mathrm{SO}}}({{\mathcal{A}}})$, we know that this map is an isomorphism. We abuse the notation slightly by using $R, U, W$ to represent the images of the elements in under this map. It is now easy to check that $R, U, W$ all have order 2 and $
U = WRW.
$ Therefore, the group ${{\mathrm{SO}}}^+(L)/\Gamma_L$ is the wreath product of ${\mathbb{Z}}/2{\mathbb{Z}}$ by ${\mathbb{Z}}/2{\mathbb{Z}}$, which is isomorphic to the dihedral group $D_8$ of order 8, and ${{\mathrm{SO}}}({{\mathcal{A}}}) \cong D_8 \times {\mathbb{Z}}/2{\mathbb{Z}}$. Using , their actions on ${{\mathcal{A}}}$ are given by $$\label{eq:AL}
  \begin{split}
R(a, b, c) &= (a, b + 2a, 2a + 2b + c),  \;
W(a, b, c) = (c, -b, a).
  \end{split}$$ To make the picture complete, we also consider the following central element in $\mathrm{O}({{\mathcal{A}}})$ $$\label{eq:mu}
  \begin{split}
  \mu: {{\mathcal{A}}}& \to {{\mathcal{A}}}\\
(a, b, c) &\mapsto (a, -b, c)    ,
  \end{split}$$ which generates $\mathrm{O}({{\mathcal{A}}}) \cong D_8 \times ({\mathbb{Z}}/2{\mathbb{Z}})^2$ along with $R, U, W$ and $\nu$.

By acting on the basis $\{{\mathfrak{e}}_h: h \in {{\mathcal{A}}}\}$, the group $\mathrm{O}({{\mathcal{A}}})$ naturally acts on ${\mathbb{C}}[{{\mathcal{A}}}]$, which commutes with the action of ${{\mathrm{Mp}}}_2({\mathbb{Z}})$ under $\rho$. Therefore, ${{\mathrm{Mp}}}_2({\mathbb{Z}})$ acts on the $\chi$-isotypic subspace ${\mathbb{C}}[{{\mathcal{A}}}]^\chi \subset {\mathbb{C}}[{{\mathcal{A}}}]$ as well, where $\chi$ is a character of order 2 on $\mathrm{O}({{\mathcal{A}}})$ defined by $$\label{eq:chi}
  \chi(R) = \chi(U) = \chi(W) = \chi(\mu) = -\chi(\nu) = -1.$$ After composing with the map ${{\mathrm{SO}}}^+(L) \to {{\mathrm{SO}}}^+(L)/\Gamma_L \hookrightarrow \mathrm{O}({{\mathcal{A}}})$, we can also view $\chi$ as a character of ${{\mathrm{SO}}}^+(L) = \Gamma^+_0(8)$. Consider the projection map $\varpi: {\mathbb{C}}[{{\mathcal{A}}}] \to {\mathbb{C}}[{{\mathcal{A}}}]^\chi$ defined by $$\label{eq:varpi}
\varpi(  v) :=  \frac{1}{|\mathrm{O}({{\mathcal{A}}})|}\sum_{s \in \mathrm{O}({{\mathcal{A}}})} \chi(s)^{-1} s(v),~ v \in {\mathbb{C}}[{{\mathcal{A}}}],
  $$ which restricts to the identity map on ${\mathbb{C}}[{{\mathcal{A}}}]^\chi \subset {\mathbb{C}}[{{\mathcal{A}}}]$. For any $h \in {{\mathcal{A}}}$, the vector $\varpi({\mathfrak{e}}_h)$ does not vanish if and only if the stabilizer $\mathrm{O}({{\mathcal{A}}})_h$ is contained in the kernel of $\chi$. Let $A \subset {{\mathcal{A}}}$ denote the subset of such elements, which is $\mathrm{O}({{\mathcal{A}}})$-invariant and decomposes into the following $\mathrm{O}({{\mathcal{A}}})$-orbits $\subset A_0$, which can be decomposed into the following five ${{\mathrm{SO}}}^+({{\mathcal{A}}})$-orbits $$  \begin{split}
      {{\mathrm{Iso}}}({{\mathcal{A}}}) =& \{(0, 0, 0)\} \sqcup \{(2, 0, 2)\} \sqcup 
\{(1, 0, 0), (0, 0, 1), (1, 2, 2), (2, 2, 1)\} \\ 
&\sqcup \{(2, 0, 0), (0, 0, 2)\} \sqcup \{(3, 0, 0), (0, 0, 3), (3, 2, 2), (2, 2, 3)\}.
  \end{split}$$ $$\label{eq:nonfixorbit}
  \begin{split}
      A &= A_0 \sqcup A_1 \sqcup A_2, \\
A_0 &:= \{(0, 1, 1), (3, 3, 0), (0, 3, 3), (1, 1, 0), (0, 3, 1), (1, 3, 0), (0, 1, 3), (3, 1, 0)\}, \\
A_1 &:= \{(1, 1, 1), (1, 3, 1),(3, 3, 3), (3, 1, 3)\}, \\
A_2 &:= \{(1, 1, 2), (2, 1, 1), (3, 3, 2),  (2, 3, 3), (3, 1, 2), (2, 3, 1),(1, 3, 2), (2, 1, 3)\}.
  \end{split}$$ Therefore, the vector space ${\mathbb{C}}[{{\mathcal{A}}}]^\chi$ is three dimensional and we identify it with ${\mathbb{C}}^3$ using the basis of orthogonal vectors $\{{\mathfrak{e}}_\ell: \ell = 0, 1, 2\}$ given by $$\label{eq:efell}
  {\mathfrak{e}}_\ell := {|\mathrm{O}({{\mathcal{A}}})|} \cdot {\varpi({\mathfrak{e}}_{(1, 1, \ell)})}.$$ After a straightforward calculation, we see that $\rho$ becomes the representation $\varrho$ in when restricted to ${\mathbb{C}}[{{\mathcal{A}}}]^\chi$.

\[lemma:subrep\] In the notations above, the map $\varpi: {\mathbb{C}}[{{\mathcal{A}}}] \to {\mathbb{C}}[{{\mathcal{A}}}]^\chi \cong {\mathbb{C}}^3$ $$\label{eq:subrep}
    \begin{split}
      {\mathbb{C}}[{{\mathcal{A}}}] &\to {\mathbb{C}}[{{\mathcal{A}}}]^\chi\\
(a_0, a_1, a_2) &\mapsto \sum_{\ell = 0, 1, 2} a_\ell {\mathfrak{e}}_{\ell}
    \end{split}$$ intertwines the representations $\rho$ and $\varrho$ from .

\[rmk:subspace\] For any $k \in \frac{1}{2} {\mathbb{Z}}$, we can view $M^!_{k, \varrho}$ as a subspace of $M^!_{k, \rho}$.

For any $v\in{\mathbb{C}}[{{\mathcal{A}}}]^\chi$ we have $v=v^\chi$ and for any $t\in O({{\mathcal{A}}})$, $(t(v))^\chi = \chi(t^{-1})v^\chi$. Hence ${\mathbb{C}}[{{\mathcal{A}}}]$ is generated by $e_{(1,1,\ell)}$ for $\ell=0,1,2$ and the above map is an isomorphism of vector spaces. The fact that it is an intertwining map follows from expressing $\rho(T)e_{(1,1,\ell)}$ and $\rho(T)e_{(1,1,\ell)}$ in terms of $e_{(1,1,\ell)}$.

Similarly, we can analyze $\mathrm{O}({{\mathcal{A}}}_j)$. For $j = 2$, this is generated by $\nu$. For $j = 1$, one needs the additional generator $$\label{eq:sigma}
  \begin{split}
    \sigma: {{\mathcal{A}}}_1 &\to {{\mathcal{A}}}_1 \\
(a, c) &\mapsto (c, a).
  \end{split}$$ The product $\mathrm{O}({{\mathcal{A}}}_1) \times \mathrm{O}({{\mathcal{A}}}_2)$ canonically embeds into $\mathrm{O}({{\mathcal{A}}})$ via , and we use $\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2) \subset \mathrm{O}({{\mathcal{A}}})$ to denote this image. When we restrict $\chi$ to $\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)$, it decomposes as $\chi_1 \otimes \chi_2$, where $$\label{eq:chij}
  \chi_1(\sigma) = - \chi_1(\nu) = - \chi_2(\nu) = 1.$$ Let $\varrho_j$ denote the restriction of $\rho_j$ to the subspace ${\mathbb{C}}[{{\mathcal{A}}}_j]^{\chi_j} \subset {\mathbb{C}}[{{\mathcal{A}}}_j]$ and $\varpi_j: {\mathbb{C}}[{{\mathcal{A}}}_j] \to {\mathbb{C}}[{{\mathcal{A}}}_j]^{\chi_j}$ the projection defined as in for $j = 1, 2$. The same analysis as before shows that ${\mathbb{C}}[{{\mathcal{A}}}_j]^{\chi_j}$ has dimension 3 and 1 for $j = 1, 2$ respectively. Therefore, $\varrho_2$ is a character of ${{\mathrm{Mp}}}_2({\mathbb{Z}})$ given by $$\label{eq:varrho2}
  \varrho_2(T) = \zeta_8,~ \varrho_2(S) = \overline{\zeta_8}.$$ Now, the identification in means that $
{\mathbb{C}}[{{\mathcal{A}}}_1]^{\chi_1} \otimes {\mathbb{C}}[{{\mathcal{A}}}_2]^{\chi_2} \subset {\mathbb{C}}[{{\mathcal{A}}}]^\chi.
$ By considering the dimensions, we see that this is in fact an equality, and intertwines the representations $\varrho_1 \otimes \varrho_2$ and $\varrho$. We can now identify ${\mathbb{C}}[{{\mathcal{A}}}_1]^{\chi_1}$ and ${\mathbb{C}}[{{\mathcal{A}}}_2]^{\chi_2}$ with ${\mathbb{C}}^3$ and ${\mathbb{C}}$ respectively via the bases $$\label{eq:ef'}
  \begin{split}
      {\mathfrak{b}}_\ell &:= \sum_{s \in \mathrm{O}({{\mathcal{A}}}_1)} \chi(s) s({\mathfrak{e}}_{(1, 0, \ell)}) \in {\mathbb{C}}[{{\mathcal{A}}}_1]^{\chi_1}, \ell = 0, 1, 2 \\
 {\mathfrak{e}}&:= {\mathfrak{e}}_{(0, 1, 0)} - {\mathfrak{e}}_{(0, 3, 0)} \in {\mathbb{C}}[{{\mathcal{A}}}_2]^{\chi_2}.
  \end{split}$$ Then ${\mathfrak{b}}_\ell \otimes {\mathfrak{e}}= {\mathfrak{e}}_\ell$ for $\ell = 0, 1, 2$ and the equality ${\mathbb{C}}[{{\mathcal{A}}}_1]^{\chi_1} \otimes {\mathbb{C}}[{{\mathcal{A}}}_2]^{\chi_2} = {\mathbb{C}}[{{\mathcal{A}}}]^\chi$ becomes ${\mathbb{C}}^3 \otimes {\mathbb{C}}= {\mathbb{C}}^3$. With respect to the basis $\{{\mathfrak{b}}_\ell: \ell = 0, 1, 2\}$, the representation $\varrho_1$ has the matrix representation $$\label{eq:varrho1}
    {\varrho_1}(T) := \begin{pmatrix}
1 & & \\ & -i & \\ && -1
\end{pmatrix}, \quad
  {\varrho_1}(S) := \frac{i}{2} \begin{pmatrix}
1 & 2 & 1 \\ 1 & 0 & -1 \\ 1 & -2 & 1
\end{pmatrix}.$$ As in Remark \[rmk:subspace\], we have $M_{k, \varrho_j} = M_{k, \rho_j}^{\chi_j} \subset M_{k, \rho_j}$ for $j = 1, 2$. We use ${\mathfrak{e}}_\ell^*, {\mathfrak{b}}_\ell^*$ and ${\mathfrak{e}}^*$ to represent basis vectors of the unitary dual representations $\varrho^*, \varrho_1^*$ and $\varrho_2^*$. Furthermore, $\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)$ acts transitively on the orbit $A_\ell$ for $\ell = 0, 1, 2$, i.e.$$\label{eq:Ai}
  A_\ell = \{s \cdot (1, 1, \ell): s \in \mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)\}$$ for $\ell = 0, 1, 2$.

This important observation enables us to prove the following result.

Heegner Divisor.
----------------

For $m \in {\mathbb{Q}}$ and $h \in {{\mathcal{A}}}$, we define the $\Gamma_L$-invariant subset $$L_{m, h} := \{\lambda \in L+h: Q(\lambda) = m\}\subseteq L^\vee.$$ When $m < 0$, the following analytic subset of ${{\mathcal{H}}}$ $$\label{eq:Zmh}
  Z_{m, h} := \{ z \in {{\mathcal{H}}}: (Z(z) , \lambda) = 0 \text{ for some } \lambda \in L_{m, h}\}$$ is $\Gamma_L$-invariant and descends to an algebraic divisor on $\Gamma_L \backslash {{\mathcal{H}}}$. Note that $Z_{m, h}$ is empty when $m > 0$. The singularities of Borcherds lifts are supported on these divisors, which are called *Heegner divisors*. For a function to be a regularized theta lift, it is then necessary for it to have singularity along Heegner divisors of a certain lattice. In our case, the singularity of ${\phi}$ appears at $z = z_\pm$. It is straightforward to check that $$(Z(z_\pm), \lambda_\pm) = 0\text{ for } \lambda_\pm := {\begin{pmatrix}
                 \pm 1/4 & 1/8\\
                 -1 & \mp 1/4
\end{pmatrix}} \in L_{-1/8, (1, \pm 1, 1)}.$$ We now have the following lemma.

\[lemma:orbit\] The set $\Gamma_L \backslash L_{-1/8, (1, \pm 1, 1)}$ has size one.

We associate to $\lambda = {\left(\begin{smallmatrix}
                 b/4 & c/8\\
                 -a & -b/4
\end{smallmatrix}\right)}\in L_{-1/8,(1,\pm 1,1)}$ the binary quadratic form $[\lambda](x,y) = [8a,4b,c](x,y)=8ax^2+4bxy+cy^2$. This identifies $L_{-1/8,(1,\pm 1,1)}$ with the set $$\mathcal{Q}^{0}_{8,-16,\pm 4} = 
\{[8a,4b,c]|\, b\equiv \pm 4\bmod {16},~\gcd(8a,b,c)=1\}$$ If $[8a,b,c]\in \mathcal{Q}^{0}_{8,-16,\pm 4}$, we must have $a\equiv c\equiv \pm 1\bmod 8$. So $\mathcal{Q}^{0}_{8,-16,\pm 4}$ is the union of the images of $L_{-1/8,(1,\pm 1,1)}$ and $L_{-1/8,(-1,\pm 1,-1)}$ under the map $\lambda\mapsto [\lambda]$. If $\gamma\in\Gamma_L$ we have $f_{\gamma\lambda}(x,y) = f((x,y)(\gamma^{-1})^T)$. Hence $\Gamma_L\setminus (L_{-1/8,(1,\pm 1,1)}\cup  L_{-1/8,(-1,\pm 1,-1)})$ is in bijection with $\Gamma_L\setminus \mathcal{Q}^0_{8,-16,\pm 4}$. We also note that while, $\Gamma_L\setminus L_{-1/8,(1,\pm 1,1)}$ and $\Gamma_L\setminus L_{-1/8,(-1,\pm 1,-1)}$ are disjoint, they are in bijection to each other via the map that sends ${\left(\begin{smallmatrix}
                 b/4 & c/8\\
                 -a & -b/4
\end{smallmatrix}\right)}$ to ${\left(\begin{smallmatrix}
                 b/4 & -c/8\\
                 a & -b/4
\end{smallmatrix}\right)}$. In [@GrossKohnenZagier §I.1] the classes of $\mathcal{Q}^0_{8,-16,\pm 4}$ modulo $\Gamma_L = \Gamma_0(8)$ are classified. In our particular case these $\Gamma_0(8)$-classes correspond to ${{\mathrm{SL}}}_2({\mathbb{Z}})$-classes of primitive binary quadratic forms of discriminant $-16$, of which there are $2$. Hence $|\Gamma_L\setminus L_{-1/8,(1,\pm 1,1)}|=|\Gamma_L\setminus \mathcal{Q}^0_{8,-16,\pm 4}|/2= 1$.

Additive Borcherds’ Lift.
-------------------------

In [@Borcherds98], Borcherds extended the input space of theta lift from ${{\mathrm{SL}}}_2$ to $\mathrm{O}(p, q)$ to include weakly holomorphic, vector-valued modular forms. The outputs are then automorphic forms on orthogonal Shimura varieties with singularities along Heegner divisors. For $(p, q) = (2, 1)$, the orthogonal Shimura variety is again the modular curve, and one can obtain a generalization of the Shimura lift to include weakly holomorphic modular forms (see [@ChoieLim16] and [@LiZemel17]). For the lattice $L$ in section \[subsec:lattice\], Borcherds’ result implies Theorem \[thm:ShimLift\].

Let $G  \in M^!_{5/2, \varrho} = M^{!, \chi}_{5/2, \rho} \subset M^!_{5/2, \rho}$ with Fourier expansion $$G(\tau) = \sum_{h \in {{\mathcal{A}}}} {\mathfrak{e}}_h \sum_{m \in {\mathbb{Z}}}  c_h(m) q^{m/8}$$ and $\rho = \rho_L$ with $L$ as in section \[subsec:lattice\]. We choose the isotropic vector ${\left(\begin{smallmatrix}
                 0 & 0\\
                 -4 & 0
\end{smallmatrix}\right)} \in L$ with $z' := {\left(\begin{smallmatrix}
                 0 & 1/8\\
                 0 & 0
\end{smallmatrix}\right)} \in L^\vee$ in the notation of Theorem 14.3 of [@Borcherds98]. The lattice $K$ is then our $L_2$, $m^+ = 2$ and $c_{h}(0) = 0$ whenever $h$ is not in the set $A$ defined in . Therefore, the Fourier expansion of the regularized theta lift in part 5 of Theorem 14.3 loc. cit. becomes $$\sum_{n > 0} \sum_{b \in {\mathbb{Z}}, b \text{ odd}, b < 0}q_z^{-(nb)/8} n \sum_{a \in {\mathbb{Z}}/4{\mathbb{Z}}} {{\mathbf{e}}}(-na) c_{(a, b, 0)}\left( {b^2}{} \right).$$ Since $G$ is in the $\chi$-isotypic component and $$\mu((1, 1, 0)) = (1, 3, 0), \nu((1, 1, 0)) = (3, 3, 0), (\mu \circ \nu)((1, 1, 0)) = (3, 1, 0),$$ we know that $c_{(1, 1, 0)}(m) = c_{(3, 3, 0)}(m) = -c_{(1, 3, 0)}(m) = - c_{(3, 1, 0)}(m)$ for all $m \in {\mathbb{Z}}$. By Lemma \[lemma:subrep\], we can write $G(\tau) =  \sum_{m \in {\mathbb{Z}}} c(G, m) q^{m/8} {\mathfrak{e}}_{\ell(m)}$ as in in the introduction with ${\mathfrak{e}}_\ell$ defined in . Then $
c(G, m) ={c_{(1, 1, 0)}(m) }
$ for all $m \in 8{\mathbb{Z}}+ 1$ Substitute this into the expression above gives us $2i {\Phi}(z, G)$.

The action of the orthogonal group on the theta kernel and imply that ${\Phi}(z, G) \mid_4 \gamma = {\Phi}(z, \gamma \cdot G)$ for any $\gamma \in {{\mathrm{PSL}}}_2({\mathbb{R}})$, hence ${\Phi}(z, G) \in M^\mathrm{mero}_{4, \chi}({{\mathrm{SO}}}^+(L))$ with ${{\mathrm{SO}}}^+(L) = \Gamma^+_0(8)$. By Theorem 6.2 of [@Borcherds98], the function $${\Phi}(z, G) - \frac{i}{16\pi^2} \sum_{\lambda \in L^\vee, \, (\lambda, Z(z_0)) = 0}\frac{ c_\lambda(Q(\lambda)) }{(\lambda, Z(z))^2}$$ is holomorphic when $z \in {{\mathcal{H}}}$ is near $z_0\in {{\mathcal{H}}}$. When $G = -G_1/4 = -\frac{q^{-1/8}}{4} {\mathfrak{e}}_1 + O(q^{1/8})$, the polar part of the expansion of ${\Phi}(z, G)$ near $z = z_\pm$ matches exactly with that of ${\phi}(z)$. Therefore, the difference ${\Phi}(z, G_1) - {\phi}(z) \in M_{4, \chi}(\Gamma^+_0(8))$ is zero by Remark \[rmk:vanish\]. The existence follows from Serre duality [@BF04], as there is no non-trivial cusp form in $S_{-1/2, \overline{\varrho}}$. In fact, we will explicitly construct it in the following section. It is unique since $M_{5/2, \varrho} = \{0\}$ by Prop. \[prop:family\] below.

Bases, Dualities and Hecke Operators. {#sec:3}
=====================================

Let ${\phi}(z) \in M^{\mathrm{mero}}_{4, \chi}(\Gamma^+_0(8))$ be the meromorphic modular form defined in , which has singularity along CM points. From Prop. \[prop:Shimura\], we saw that the Shimura lift ${\Phi}(z, g)$ of a form $g \in M^{!}_{5/2, \varrho} \subset M^{!}_{5/2, \rho}$ is in $M^{\mathrm{mero}}_{4, \chi}(\Gamma^+_0(8))$ with singularity on Heegner divisors.

.... In this section, we will construct the ${\mathbb{Z}}$-basis $\{G_d: d \equiv 1, 3 \text{ or } 7 \bmod{8}\}$ of ${\mathbb{M}}^{!}_{5/2, \varrho}$, and show that it has duality with respect to a ${\mathbb{Z}}$-basis of ${\mathbb{M}}^!_{-1/2, {\varrho^*}}$, where $\varrho^*$ is the unitary dual of $\varrho$, i.e., $\varrho^*(T)$ and $\varrho^*(S)$ are the conjugate transpose of $\varrho(T)$ and $\varrho(S)$ respectively. Using Hecke operators on vector-valued modular forms [@BS10], we will prove Conjecture \[conj:main\].

Two Bases of Modular Forms. {#subsec:family}
---------------------------

In the notations of the previous section, we will prove the following result. Let $L, \rho, {{\mathcal{A}}}, \chi$ and $M^{!, \chi}_{5/2, \rho}$ be as in section \[subsec:lattice\]. In this section, we will study a *unique* family of vector-valued, weakly holomorphic modular form $g_d  \in M^{!, \chi}_{5/2, \rho}$ for odd $d \in (8{\mathbb{N}}+ 1)$. The main result is as follows.

\[prop:family\] For every $d \in {\mathbb{N}}$ congruent to 1, 3 or 7 modulo 8, there exists a unique $G_d$ in ${\mathbb{M}}^{!}_{5/2, \varrho}$ with the Fourier expansion $$\label{eq:Gdqexp}
    G_d(\tau) = \mathfrak{q}^{-d}    + O(1) = \sum_{D \ge -d,~ D \equiv 1, 5, 7 \bmod{8}} B(D, d) \mathfrak{q}^{D},~\mathfrak{q}^D := q^{D/8} {\mathfrak{e}}_{\ell(D)}.$$ For every $D \in {\mathbb{N}}$ congruent to 1, 5 or 7 modulo 8, there exists a unique $F_D$ in ${\mathbb{M}}^{!}_{-1/2, {\varrho^*}}$ with the Fourier expansion $$\label{eq:FDqexp}
    F_D(\tau) = \tilde{\mathfrak{q}}^{-D}    + O(1) = \sum_{d \ge -D,~ d \equiv 1, 3, 7 \bmod{8}} A(D, d) \tilde{\mathfrak{q}}^{d},~ \tilde{\mathfrak{q}}^d := q^{d/8} {\mathfrak{e}}^*_{\ell(-d)}.$$ Furthermore, $A(D, d) + B(D, d) = 0$ for all $D, d \in {\mathbb{N}}$.

To prove this proposition, we will first reduce it to the case of weight one using the presence of the unary theta function $$\label{eq:theta}
\vartheta(\tau) := \sum_{b \in {\mathbb{Z}}/4{\mathbb{Z}}} {\mathfrak{e}}_{(0, b, 0)} \sum_{n \in 4{\mathbb{Z}}+ b} n q^{n^2} = \eta^3(\tau) {\mathfrak{e}}\in S_{3/2, \varrho_2} \cap {\mathbb{M}}^!_{3/2, \varrho_2},$$ which only vanishes at the cusp infinity. Since tensoring with ${\mathfrak{e}}$ gives an isomorphism between ${\mathbb{C}}[{{\mathcal{A}}}_1]^{\chi_1}$ and ${\mathbb{C}}[{{\mathcal{A}}}]^{\chi}$ as ${{\mathrm{Mp}}}_2({\mathbb{Z}})$-modules under the representations $\varrho_1$ and $\varrho$, we obtain the following lemma.

\[lemma:wt1\] In the notations above, the following map is an isomorphism of ${\mathbb{Z}}$-modules $$\label{eq:varthetaisom}
    \begin{split}
      \iota_\vartheta: {\mathbb{M}}^{!}_{1, \varrho_1} & \to {\mathbb{M}}^{!}_{5/2, \varrho} \\
      f &\mapsto f   \otimes \vartheta .
    \end{split}$$ Furthermore, the preimage of ${\mathbb{M}}_{5/2, \varrho}$ under this isomorphism is the trivial subspace ${\mathbb{M}}^{}_{1, \varrho_1}$ and hence ${\mathbb{M}}_{5/2, \varrho}$ is also trivial.

\[rmk:-1/2\] By taking the conjugate transpose of the representation $\varrho_1$ and dividing by $\eta^3(\tau)$, we also obtain an explicit isomorphism ${\mathbb{M}}^!_{1, {\varrho^*_1}}$ and ${\mathbb{M}}^!_{-1/2, {\varrho^*}}$. However, we will see later that ${\mathbb{M}}_{1, {\varrho^*_1}}$ is one dimensional.

The first part follows from the argument above. The inverse map is simply dividing by $\eta^3(\tau)$ on each component. For the second part, notice that ${\mathbb{M}}_{5/2, \varrho} \subset S_{5/2, \varrho}$ and the order of vanishing at each component of any $G \in {\mathbb{M}}_{5/2, \varrho}$ is at least $q^{1/8}$. Therefore, the result is still holomorphic after dividing $G$ by $\eta^3(\tau)$.

For any $f = \sum_{\ell = 0,1, 2} f_\ell{\mathfrak{b}}_\ell \in M_{1, \varrho_1}$, the function $f_1(4\tau)$ is in the 7-dimensional space $M_{1}(\Gamma_1(16))$ with the Fourier expansion $$f_1(4\tau) = \sum_{n \ge 0} c(n) q^{4n + 3}$$ at infinity. A quick calculation with SageMath [@sagemath] shows that $f_1$ is identically zero. Then $$f \mid_1 S = (f_0 \mid S) {\mathfrak{e}}_0 + (f_2 \mid S) {\mathfrak{e}}_2 = \varrho_1(S) \cdot f = \frac{i}{2} ((f_0 + f_2)({\mathfrak{e}}_0 + {\mathfrak{e}}_2) + 2(f_0 - f_2) {\mathfrak{e}}_1)$$ implies that $f_0 = f_2$. Since $f_\ell(\tau + 1) = (-i)^\ell f_\ell(\tau)$, we conclude that $f_0$ and $f_2$ are both identically zero. Therefore, $M^{}_{1, \varrho_1}$ is trivial and so is ${\mathbb{M}}_{1, \varrho_1}\subset M^{}_{1, \varrho_1}$. The same procedure shows that $S_{1, {\varrho^*_1}}$ is trivial and $M_{1, {\varrho^*_1}}$ is at most one dimensional. We will construct this non-trivial element later.

The proof of this result follows from a sequence of lemmas, and will be given at the end of the section. First, it is helpful to consider the decomposition ${{\mathcal{A}}}= {{\mathcal{A}}}_1 \oplus {{\mathcal{A}}}_2$ as in section \[subsec:lattice\]. It induces an isomorphism $\rho \cong \rho_1 \otimes \rho_2$ and embedding $\mathrm{O}({{\mathcal{A}}}_1) \times \mathrm{O}({{\mathcal{A}}}_2) \hookrightarrow \mathrm{O}({{\mathcal{A}}})$, whose image was denoted by $\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)$. When we restrict $\chi$ to $\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)$, it decomposes as $\chi_1 \otimes \chi_2$, where $\chi_j(\nu) = -1$ for $j = 1, 2$ and $\chi_1(\sigma) = 1$. We would hope to write modular forms in $M^{!, \chi}_{5/2, \rho}$ as sums of tensor products of modular forms in $M^{!, \chi_1}_{k, \rho_1}$ and $M^{!, \chi_2}_{5/2 - k, \rho_2}$, where character in the superscript denotes the corresponding isotypic subspace. In this case, we are very lucky due to the presence of the unary theta function

The first lemma identifies $M^{!, \chi}_{5/2, \rho}$ and $M^{!, \chi_1}_{1, \rho_1}$ through tensoring with $\vartheta$.

\[lemma:isom\] In the notations above, the following map is an isomorphism of ${\mathbb{C}}$-vector spaces $$\label{eq:varthetaisom}
    \begin{split}
      \iota_\vartheta: M^{!, \chi_1}_{1, \rho_1} & \to M^{!, \chi}_{5/2, \rho} \\
      f &\mapsto f   \otimes \vartheta .
    \end{split}$$ Furthermore, the preimage of $M^\chi_{5/2, \rho}$ under this isomorphism is the trivial subspace $M^{\chi_1}_{1, \rho_1}$.

\[rmk:Qstr\] From the main result of [@McGraw03], we can deduce that $M^{!, \chi_1}_{1, \rho_1}$ and $M^{!, \chi}_{5/2, \rho}$ both have a basis of modular forms with rational Fourier coefficients. From the product expansion of $\eta(\tau)$, it is immediate that $\iota_\vartheta$ preserves this rational structure. Furthermore, it even preserves integrality, i.e. $f \in M^{!}_{1, \rho_1}$ has integral Fourier coefficients if and only if $\iota_\vartheta(f)$ does.

It is not difficult to see that $\iota_\vartheta$ is an injective map from $M^{!}_{1, \rho_1}$ to $M^!_{5/2, \rho}$. If $f = \sum_{h \in {{\mathcal{A}}}_1} f_h {\mathfrak{e}}_h \in M^{!, \chi_1}_{1, \rho_1}$, then it is straightforward to check that $f_h$ vanishes identically unless $h$ is in the $\mathrm{O}({{\mathcal{A}}}_1)$-orbit of $(1, 0), (1, 1)$ or $(1, 2)$. Using the character $\chi_1$, we can then write $$\label{eq:f2}
f = \sum_{\ell = 0, 1, 2} f_{(1, \ell)} 
\sum_{s \in \mathrm{O}({{\mathcal{A}}}_1)/\mathrm{O}({{\mathcal{A}}}_1)_{(1, \ell)}} \chi_1(s)  {\mathfrak{e}}_{s \cdot (1, \ell)},$$ where for a group $G$ acting on a set $X$, we denote by $G_x \subset G$ the stabilizer of $x \in X$. Let $W := \oplus_{\ell = 0, 1, 2}{\mathbb{C}}{\mathfrak{e}}_\ell$ be a ${\mathbb{C}}$-vector space of dimension 3 and define a function $\iota(f) : {{\mathcal{H}}}\to W$ by $$\iota(f) := \sum_{\ell = 0, 1, 2} f_{(1, \ell)} {\mathfrak{e}}_\ell.$$ It is easy to check that $\iota$ gives an isomorphism between $M^{!, \chi_1}_{1, \rho_1}$ and $M^!_{1, \widetilde{\rho_1}}$ as ${\mathbb{C}}$-vector spaces, where $  \widetilde{\rho_1}$ is a representation of ${{\mathrm{SL}}}_2({\mathbb{Z}})$ on $W$ given by $$\label{eq:trho1}
  \widetilde{\rho_1}(T) := \begin{pmatrix}
1 & & \\ & -i & \\ && -1
\end{pmatrix}, \quad
  \widetilde{\rho_1}(S) := \frac{i}{2} \begin{pmatrix}
1 & 2 & 1 \\ 1 & 0 & -1 \\ 1 & -2 & 1
\end{pmatrix}.$$ $$\widetilde{\rho_1}(T) := \begin{pmatrix}
1 & & \\ & -i & \\ && -1
\end{pmatrix}, \quad
  \widetilde{\rho_1}(S) := \frac{i}{2} \begin{pmatrix}
1 & 1 & 1 \\ 2 & 0 & -2 \\ 1 & -1 & 1
\end{pmatrix}.$$

It is clear that $\iota$ preserves the holomorphicity condition.

From the expression for $f$ in , it follows that $$\label{eq:iotatheta}
  \begin{split}
\iota_\vartheta(f)(\tau) &=  \eta^3(\tau)  
\sum_{\ell = 0, 1, 2}
f_{(1, \ell)}(\tau) {\mathfrak{e}}^\chi_{(1, 1, \ell)}, \\
{\mathfrak{e}}^\chi_h &:= \sum_{s \in \mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)/\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)_h} \chi(s) {\mathfrak{e}}_{s \cdot h} \in {\mathbb{C}}[{{\mathcal{A}}}].    
  \end{split}$$ Suppose from now on $h \in  \{ (1, 1, 0), (1, 1, 1),(1, 1, 2)\}$, as $\mathrm{O}({{\mathcal{A}}})_h \subset \ker(\chi)$ in these cases and the vector ${\mathfrak{e}}^\chi_h$ is well-defined. By examining , we see that $\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)$ acts transitively on $A_\ell$ for $\ell = 0, 1, 2$. Therefore, the following inclusion of left cosets is a bijection $$\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)/\mathrm{O}({{\mathcal{A}}}_1, {{\mathcal{A}}}_2)_h \stackrel{\cong}{\hookrightarrow} \mathrm{O}({{\mathcal{A}}})/\mathrm{O}({{\mathcal{A}}})_h.$$ Furthermore, the character $\chi$ is a well-defined function on these sets above. Therefore, we can rewrite $$\label{eq:efchi}
{\mathfrak{e}}^\chi_h =  \sum_{s \in \mathrm{O}({{\mathcal{A}}})/\mathrm{O}({{\mathcal{A}}})_h} \chi(s) {\mathfrak{e}}_{s \cdot h}.$$ Given any $r \in \mathrm{O}({{\mathcal{A}}})$, it is easy to check $$\begin{aligned}
  r \cdot {\mathfrak{e}}^\chi_h
&=  r \cdot \sum_{s \in \mathrm{O}({{\mathcal{A}}})/\mathrm{O}({{\mathcal{A}}})_h} \chi(s) {\mathfrak{e}}_{s \cdot h}
= \chi(r) \sum_{s \in \mathrm{O}({{\mathcal{A}}})/\mathrm{O}({{\mathcal{A}}})_h} \chi(rs) {\mathfrak{e}}_{(rs) \cdot h}
= \chi(r) {\mathfrak{e}}_h^\chi,\end{aligned}$$ which implies that $\iota_\vartheta(f) \in M^{!, \chi}_{5/2, \rho}$.

Now, we want to prove that $\iota_\vartheta$ is surjective. Suppose $g = \sum_{h \in {{\mathcal{A}}}} g_h {\mathfrak{e}}_h \in M^{!, \chi}_{5/2, \rho}$ is a vector-valued modular form. If $h$ is not in the set $A$ in , then there exists $s \in \mathrm{O}({{\mathcal{A}}})_h$ with $\chi(s) = -1$, and $s \cdot g = -g$ implies that $g_h$ vanishes identically. For $\ell = 0, 1, 2$, define $$f(\tau) := \eta^{-3}(\tau) \sum_{\ell = 0, 1, 2}  g_{(1, 1, \ell)}(\tau) {\mathfrak{e}}_\ell: {{\mathcal{H}}}\to W.$$ Then one can use the modularity of $g$ and the eta function to verify that $f \in M^!_{1, \widetilde{\rho_1}}$. From the first half of the proof, we have $\iota(f) \in M^{!, \chi_1}_{1, \rho_1}$ and $\iota_\vartheta(\iota(f)) = g$. If $g \in M^\chi_{1, \rho}$, then $g_{(1, 1, \ell)}(8\tau) = O(q)$ for every $\ell = 0, 1, 2$, and the function $f$ constructed above is in $M_{1, \widetilde{\rho_1}}$. This finishes the proof.

With the lemma above, it suffices to study the space $M^{!}_{1, \varrho_1} = M^{!, \chi_1}_{1, \rho_1}$, which has the following property.

\[lemma:wt1family\] For every $m$ in ${\mathbb{N}}\backslash (4{\mathbb{N}}-1)$, there exists a unique $g_m \in {\mathbb{M}}^{!}_{1, \varrho_1}$ with the Fourier expansion $$g_m(\tau) =  q^{-m/4} {\mathfrak{b}}_{m \bmod 4} + O(1).$$ For every $m$ in ${\mathbb{N}}\cup \{0\} \backslash (4{\mathbb{N}}-3)$, there exists a unique $f_m \in {\mathbb{M}}^{!}_{1, {\varrho^*_1}}$ with the Fourier expansion $$f_m(\tau) = q^{-m/4} {\mathfrak{b}}^*_{-m \bmod 4} + O(q^{1/4}).$$

The uniqueness easily follows from the previous lemma, where we saw that $M^{}_{1, \varrho_1}$ and $S^{}_{1, {\varrho^*_1}}$ are trivial. For the existence, we will give an explicit construction, which can be implemented numerically.

Since $\rho_1$ is the representation attached to a scaled hyperbolic plane, we will use Lemma 2.6 in [@Borcherds98] to construct $\tilde{g}_m$ from scalar-valued modular forms in $M^!_1(\Gamma_1(4))$. For $\tilde{g} \in M^!_1(\Gamma_1(4))$ and $(a, c) \in {{\mathcal{A}}}_1$, define $$\label{eq:sc2vv}
{{\mathscr{B}}}(\tilde{g})_{(a,c)} := \sum_{d \in {\mathbb{Z}}/4{\mathbb{Z}}, (c, d, 4) = 1} i^{ad} \tilde{g} \mid_1 {\begin{pmatrix}
                 * & *\\
                 c & d
\end{pmatrix}} + 
\sum_{d \in {\mathbb{Z}}/4{\mathbb{Z}}, (a, d, 4) = 1} i^{cd} \tilde{g} \mid_1 {\begin{pmatrix}
                 * & *\\
                 a & d
\end{pmatrix}}.$$ Then by the same proof of Lemma 2.6 in [@Borcherds98], we know that $$\label{eq:Bsc}
 {{\mathscr{B}}}(\tilde{g}) = \sum_{h \in {{\mathcal{A}}}_1} {{\mathscr{B}}}(\tilde{g})_h {\mathfrak{e}}_h = {{\mathscr{B}}}(\tilde{g})_{(1,0)} {\mathfrak{e}}_0 + \frac{{{\mathscr{B}}}(\tilde{g})_{(1,1)}}{2} {\mathfrak{e}}_1 + {{\mathscr{B}}}(\tilde{g})_{(1,0)} {\mathfrak{e}}_0$$ is in $M^{!}_{1, \varrho_1} = M^{!, \chi_1}_{1, \rho_1}$. For $(a, c) = (1,0), (1, 1)$ and $(1, 2)$, we can explicitly write ${{\mathscr{B}}}(\tilde{g})_{(a, c)}$ as $$\begin{aligned}
  {{\mathscr{B}}}(\tilde{g})_{(1, 0)} &= 2i \cdot \tilde{g} + \sum_{d = 0}^3 \tilde{g} \mid_1 S T^d, \quad
  {{\mathscr{B}}}(\tilde{g})_{(1, 1)} = 2 \sum_{d = 0}^3 i^d (\tilde{g} \mid_1 ST^d), \\
  {{\mathscr{B}}}(\tilde{g})_{(1, 2)} &= 2i (\tilde{g} \mid_1 ST^{2}S^{-1}) + \sum_{d = 0}^3 (-1)^d (\tilde{g} \mid_1 S T^d).\end{aligned}$$ So the Fourier expansion of ${{\mathscr{B}}}(\tilde{g})$ is directly related to those of $\tilde{g}$ at the three cusps of $\Gamma_1(4)$.

To construct the family $g_m$, we start with the Eisenstein series $$\tilde{g}_0(\tau)=iE_1^{\varphi,\textbf{1}}=iL(0,\varphi) + 2i\sum_{n\geq 1}\sum_{m|n}\varphi(m)q^n=\frac{i}{2} (1 + 4q + O(q^2)) \in \frac{i}{2} {\mathbb{Z}}\llbracket q \rrbracket$$ that generates ${M}_1(\Gamma_1(4))$. Here $\varphi$ is the Kronecker symbol modulo $4$. The Fourier expansions of $E_1^{\varphi,\textbf{1}}$ at the cusps of $\Gamma_1(4)$ are given in [@Weisinger1977 Chapter 2], from which one calculates $$\begin{aligned}
\tilde{g}_0|_1 S &= -\frac{i}{2} \tilde{g}_0\left(\frac{\tau}{4}\right) = \frac{1}{4}(1+4q^{1/4}+O(q^{2/4}))\in  \frac{1}{4} {\mathbb{Z}}\llbracket q^{1/4} \rrbracket,\\
\tilde{g}_0|_1 ST^2 S^{-1} &=
2i\sum_{n\geq 1~ \mathrm{odd}}\sum_{m|n}\varphi(m) q^{n/2}
 = 2iq^{1/2}(1+2q^2+O(q^4)) \in 2i q^{1/2} {\mathbb{Z}}\llbracket q \rrbracket.\end{aligned}$$ Since $M_{1,\varrho_1}=\{0\}$, $\tilde{g}_0$ must lift to $0$ under ${{\mathscr{B}}}$ and using the Fourier expansions we can check that this is indeed the case. Let $\psi$ be the Hauptmodul from . Then $$\label{eq:tphi}
  {{\tilde{\psi}}}(\tau) := \psi(-1/(4\tau)) = 2^4 \eta(\tau)^{-16} \eta(2\tau)^{24} \eta(4\tau)^{-8} = 2^4(1 + O(q)) \in 2^4{\mathbb{Z}}\llbracket q \rrbracket$$ is also a Hauptmodul of $X_1(4)$. At the other cusps, it has the expansions $$\begin{aligned}
  ({{\tilde{\psi}}}\mid S)(\tau) &= \eta(\tau)^{-16} \eta(\tau/2)^{24} \eta(\tau/4)^{-8} = q^{-1/4} + 8 + O(q^{1/4}) \in q^{-1/4} {\mathbb{Z}}\llbracket q^{1/4} \rrbracket, \\
({{\tilde{\psi}}}\mid ST^2S^{-1})(\tau) &= -2^8 \eta(\tau)^{-8} \eta(4\tau)^{8} = -2^8q ( 1 + O(q)) \in 2^8 q {\mathbb{Z}}\llbracket q \rrbracket.\end{aligned}$$ Multiplying $\tilde{g}_0$ with monic polynomials in ${{\tilde{\psi}}}$ with integral coefficients, we can recursively construct $\tilde{g}_m \in M^!_1(\Gamma_1(4))$ such that $    2i \cdot \tilde{g}_m(\tau), iq^{-1/2} (\tilde{g}_m \mid_1 ST^2S^{-1})(\tau)
 \in {\mathbb{Z}}\llbracket q \rrbracket$ and $$\label{eq:cuspexp}
(\tilde{g}_m \mid S)(\tau) = \frac{1}{4} ( q^{-m/4} + O(q^{1/4})) \in \frac{1}{4} {\mathbb{Z}}\llbracket q^{1/4} \rrbracket
$$ for any $m \in {\mathbb{N}}$. For example $$\label{eq:tg1}
 \tilde{g}_1 = \tilde{g}_0 ({{\tilde{\psi}}}- 12).$$ Let $a_m(n) \in {\mathbb{Z}}$ be the $(n/4)$[^th^ ]{}Fourier coefficient of $4\tilde{g}_m \mid S$. Then the ${\mathfrak{e}}_0, {\mathfrak{e}}_1$ and ${\mathfrak{e}}_2$ components of $g_m := {{\mathscr{B}}}(\tilde{g}_m)$ are given by $$\begin{aligned}
 & 2i \cdot \tilde{g}_m + \sum_{n \in 4 {\mathbb{Z}}} a_m(n) q^{n/4},~
\sum_{n \in 4{\mathbb{Z}}-1} a_m(n) q^{n/4}, ~
 2i (\tilde{g}_m \mid_1 ST^{2}S^{-1}) + \sum_{n \in 4{\mathbb{Z}}+ 2} a_m(n) q^{n/4}
$$ respectively. Note that $g_m$ is identically zero if $m \equiv 3 \bmod{4}$. Therefore the $g_m$’s satisfy the condition in the lemma. The family $\{f_m: m \in {\mathbb{N}}\cup \{0\} \backslash (4{\mathbb{N}}-3)\} \subset {\mathbb{M}}^!_{1, {\varrho^*_1}}$ can be constructed similarly.

For $d \equiv 1, 3, 7 \bmod{8}$, there is a unique $G_d \in M^!_{5/2, \varrho}$ satisfying . By Lemma \[lemma:wt1\], there exists $g \in M^!_{1, \varrho_1}$ such that $\iota_\vartheta(g) = G_d$. Since the principal part of $g$ has integral Fourier coefficients, we can express it as an integral linear combination of the $g_m$’s from Lemma \[lemma:wt1family\]. Therefore $g$ is contained in ${\mathbb{M}}^!_{1, \varrho_1}$ and $G_d = \iota_\vartheta(g)$ is contained in ${\mathbb{M}}^!_{5/2, \varrho}$. The same proof works for the $F_D$’s.

The last statement is a simple consequence that every weakly holomorphic modular form of weight 2 on ${{\mathrm{SL}}}_2({\mathbb{Z}})$ has vanishing constant term, as it is the derivative of a modular function. If we view ${\mathbb{C}}[{{\mathcal{A}}}_L] \times {\mathbb{C}}[{{\mathcal{A}}}_L]$ as an ${{\mathrm{SL}}}_2({\mathbb{Z}})$-module with respect to $\rho_L \otimes {\rho^*_L}$, then the hermitian inner product $\langle \cdot, \cdot \rangle: {\mathbb{C}}[{{\mathcal{A}}}_L] \times {\mathbb{C}}[{{\mathcal{A}}}_L] \to {\mathbb{C}}$ in is ${{\mathrm{SL}}}_2({\mathbb{Z}})$-linear, whose action on ${\mathbb{C}}$ is trivial. Therefore, $\langle F_D, G_d \rangle$ is in $M_2^!$ and its constant term is given by $\sum_{D', d' \in {\mathbb{Z}}, D' + d' = 0} A(D, d') B(D', d)$. We are then done by and .

Since the proof is constructive, we can explicitly give the Fourier expansions of $G_d$ for any particular $d$. For example when $d = 1$, the ${\mathfrak{e}}_0$-component of the modular form $$\label{eq:G1}
 G_{1} = \iota_\vartheta(g_1) = \iota_\vartheta({{\mathscr{B}}}(\tilde{g}_1))$$ has the Fourier expansion $$\begin{aligned}
\eta(\tau)^3 {{\mathscr{B}}}(\tilde{g}_1)_{(1, 0)} &= 
 -4(q^{1/8} + 129 q^{9/8} + 1144 q^{17/8} + 5625 q^{25/8} + O(q^{33/8})),
$$ and one can check both numerically and from the proof of Theorem \[thm:ShimLift\] that ${\Phi}(z, G_1) = -4 {\phi}(z)$.

Hecke operators
---------------

In [@BS10 Theorem 4.10] Hecke operators on vector-valued modular forms are defined. They act on $G\in M_{5/2,\rho}^!$ as follows. If $$G(\tau)=\sum_{h\in\mathcal{A}} {\mathfrak{e}}_{h} \sum_{n\in{\mathbb{Z}}}c_h(n)q^{n/8}$$ and $p$ is an odd prime, then multiplying by $p$ on ${{\mathcal{A}}}$ is an isometry and $$(G|T_{p^2})(\tau) = \sum_{h\in\mathcal{A}}{\mathfrak{e}}_h \sum_{n\in {\mathbb{Z}}}b_h(n) q^{n/8},~
b_h(n) := 
c_{ph}(p^2n)+
p\left(\frac{n}{p}\right)c_h(n)+p^{3}c_{p^{-1} h}(n/p^2),$$ with $\left(\frac{n}{p}\right)$ the Kronecker symbol. It is easily checked that multiplying by $p$ acts as either the identity or $\nu$ on ${{\mathcal{A}}}$. Therefore, we have $\varpi(p \cdot v) = \varpi(v)$ for any $v \in {\mathbb{C}}[{{\mathcal{A}}}]$ and $T_{p^2}$ preserves the subspace $M^{!}_{5/2, \varrho} \subset M^!_{5/2, \rho}$ and the lattice ${\mathbb{M}}^{!}_{5/2, \varrho} \subset M^!_{5/2, \varrho}$. Given $G(\tau) = \sum_{n \in {\mathbb{Z}}} c(G, n)\mathfrak{q}^n$ in $M^!_{5/2, \varrho}$, the Hecke operator $T_{p^2}$ then acts as $$\label{eq:Heckep}
(G|T_{p^2})(\tau) = \sum_{n \in {\mathbb{Z}}}
\left(c(G, p^2n)+
p\left(\frac{n}{p}\right)c(G, n)+p^{3}c(G, n/p^2) \right) \mathfrak{q}^n.$$ Using Prop. \[prop:family\] and comparing the principal parts, we can deduce the equality $$\label{eq:Heckeeq}
  G_d \mid T_{p^2} = p^3 G_{p^2d} + p \left(\frac{-d}{p}\right) G_d + G_{d/p^2}$$ for any $G_d \in {\mathbb{M}}^!_{5/2, \varrho}$, where $G_{d/p^2}$ is zero if $p^2 \nmid d$. This leads to the following result, which implies the second half of Theorem \[thm:HeckeSc\]

\[prop:divisibility\] Let $G_d(\tau) = \sum_{n \in {\mathbb{Z}}} B(n, d) \mathfrak{q}^n \in {\mathbb{M}}^!_{5/2, \varrho}$ be as in Prop. \[prop:family\]. Suppose $d \equiv 1, 3, 7\bmod{8}$ is square-free. Then $$\label{eq:divisibility}
B(p^2 n, d) \equiv p \left( \left( \frac{-d}{p} \right) - \left( \frac{n}{p} \right) \right) B(n, d) \bmod{p^3}$$ for any $n \in {\mathbb{N}}$ and odd prime $p$. In particular, we have $B(m^2, d) \in m{\mathbb{Z}}$ for all $m \in {\mathbb{N}}$.

Since $d$ is square-free, becomes $  G_d \mid T_{p^2} = p^3 G_{p^2d} + p \left(\frac{-d}{p}\right) G_d$. Comparing the $n$[^th^ ]{}Fourier coefficients of both sides then gives us the congruence . The last claim then follows since $B(m, d) = 0$ whenever $2\mid m$.

Proposition \[prop:divisibility\] now leads to a proof of Conjecture \[conj:main\]. A more detailed analysis reveals higher power congruences:

\[cor:p\^3 congruence\] Let $d$ be square-free and ${{\Phi}}(z, G_d) = \sum_{n \in 2{\mathbb{N}}- 1} a_d(n)q_z^n$. Then for every odd prime $p$ we have $$\label{eq:p3congruence}
a_d(p)\equiv p\left( \frac{d}{p} \right) B(1, d)\bmod p^3.$$

The coefficient $a_d(p)$ is given by and equals $(-1)^{(p-1)/2}(p\cdot B(1, d) + B(p^2, d))$. By Prop. \[prop:divisibility\] we have $B(p^2, d) \equiv p\left( \left( \frac{-d}{p} \right) - 1 \right) B(1, d) \bmod{p^3}$ and follows from $(-1)^{(p-1)/2}=\left(\frac{-1}{p}\right)$ and the multiplicativity of the Kronecker symbol.

Finally, we record some examples of the bases. Two instances of Zagier duality are marked in color. $$\begin{aligned}
  G_1(\tau) &= \mathfrak{q}^{-1} {{\color{blue}-4}}\mathfrak{q}  {{\color{blue}+112}}\mathfrak{q}^{5} {{\color{blue}+19}}\mathfrak{q}^{7} {{\color{blue}-516}}\mathfrak{q}^{9} {{\color{blue}+1712}}\mathfrak{q}^{13}  {{\color{blue}-87}} \mathfrak{q}^{15}+ O(q^{16}),\\
  G_3(\tau) &= \mathfrak{q}^{-3} {{\color{red}- 4}}\mathfrak{q} {{\color{red}- 267}}\mathfrak{q}^{5} {{\color{red}+ 1024}}\mathfrak{q}^{7} {{\color{red}- 3012}}\mathfrak{q}^{9} {{\color{red}- 19666}}\mathfrak{q}^{13} {{\color{red}+ 44032}} \mathfrak{q}^{15} + O(q^{16}),\\
    G_7(\tau) &= \mathfrak{q}^{-7} - 7\mathfrak{q} - 3136\mathfrak{q}^{5} - 20480\mathfrak{q}^{7} - 102396\mathfrak{q}^{9} - 1546048\mathfrak{q}^{13}  - 5074944 \mathfrak{q}^{15} + O(q^{16}), \\                          
  G_9(\tau) &= \mathfrak{q}^{-9}- 20\mathfrak{q} + 16944\mathfrak{q}^{5} - 172\mathfrak{q}^{7} - 854548\mathfrak{q}^{9} + 18047344\mathfrak{q}^{13}  + 5031 \mathfrak{q}^{15} + O(q^{16}),\\
  G_{11}(\tau) &= \mathfrak{q}^{-11}- 12\mathfrak{q} - 21303\mathfrak{q}^{5} + 216064\mathfrak{q}^{7} - 1566540\mathfrak{q}^{9} - 44627503\mathfrak{q}^{13} + 193840128 \mathfrak{q}^{15} + O(q^{16}),\\
  G_{15}(\tau) &= \mathfrak{q}^{-15}- 25\mathfrak{q} - 111552\mathfrak{q}^{5} - 1617920\mathfrak{q}^{7} - 15953955\mathfrak{q}^{9} - 770664640\mathfrak{q}^{13} - 4226125824 \mathfrak{q}^{15} + O(q^{16}).
            $$

$$\begin{aligned}
  F_{1}(\tau) &=   \tilde{\mathfrak{q}}^{-1}  {{\color{blue}+4}}\tilde{\mathfrak{q}} {{\color{red}+4}}\tilde{\mathfrak{q}}^{3} + 7\tilde{\mathfrak{q}}^{7} + 20\tilde{\mathfrak{q}}^{9} + 12\tilde{\mathfrak{q}}^{11} + 25\tilde{\mathfrak{q}}^{15} + O(q^{16}), \\
  F_{5}(\tau) &=   \tilde{\mathfrak{q}}^{-5} {{\color{blue}-112}}\tilde{\mathfrak{q}} {{\color{red}+267}}\tilde{\mathfrak{q}}^{3} + 3136\tilde{\mathfrak{q}}^{7} - 16944\tilde{\mathfrak{q}}^{9} + 21303\tilde{\mathfrak{q}}^{11} + 111552\tilde{\mathfrak{q}}^{15} + O(q^{16}),\\
  F_{7}(\tau) &=   \tilde{\mathfrak{q}}^{-7} {{\color{blue}-19}}\tilde{\mathfrak{q}} {{\color{red}- 1024}}\tilde{\mathfrak{q}}^{3} + 20480\tilde{\mathfrak{q}}^{7} + 172\tilde{\mathfrak{q}}^{9} - 216064\tilde{\mathfrak{q}}^{11} + 1617920\tilde{\mathfrak{q}}^{15} + O(q^{16})
 \\
  F_{9}(\tau) &=   \tilde{\mathfrak{q}}^{-9} {{\color{blue}+ 516}}\tilde{\mathfrak{q}} {{\color{red}+ 3012}}\tilde{\mathfrak{q}}^{3} + 102396\tilde{\mathfrak{q}}^{7} + 854548\tilde{\mathfrak{q}}^{9} + 1566540\tilde{\mathfrak{q}}^{11} + 15953955\tilde{\mathfrak{q}}^{15} + O(q^{16}),
\\
  F_{13}(\tau) &=   \tilde{\mathfrak{q}}^{-13}{{\color{blue}- 1712}}\tilde{\mathfrak{q}} {{\color{red}+ 19666}}\tilde{\mathfrak{q}}^{3} + 1546048\tilde{\mathfrak{q}}^{7} - 18047344\tilde{\mathfrak{q}}^{9} + 44627503\tilde{\mathfrak{q}}^{11} + 770664640\tilde{\mathfrak{q}}^{15} + O(q^{16}),
\\
  F_{15}(\tau) &=   \tilde{\mathfrak{q}}^{-15} {{\color{blue}+ 87}}\tilde{\mathfrak{q}} {{\color{red}- 44032}}\tilde{\mathfrak{q}}^{3} + 5074944\tilde{\mathfrak{q}}^{7} - 5031\tilde{\mathfrak{q}}^{9} - 193840128\tilde{\mathfrak{q}}^{11} + 4226125824\tilde{\mathfrak{q}}^{15} + O(q^{16}).
            $$

[^1]: Both authors are partially funded by the LOEWE research unit USAG and the DFG grant BR-2163/4-2. The first author is also partially supported by an NSF postdoctoral fellowship.
