First published in Asian J Math, 8 (2004) no. 1, 161-172, special issue dedicated to Yum-Tong Siu on his 60th birthday. Partially supported by NSF Mathematical Sciences Postdoctoral Research Fellowship.
<ph f="cmbx">Cubic equations for the hyperelliptic locus</ph>

### Samuel Grushevsky

Mathematics Department, Princeton University, Fine Hall, Washington Road, Princeton, NJ 08544, USA E-mail address : sam@math.princeton.edu
• Abstract. We discuss the conjecture of Buchstaber and Krichever from [BK2that the multi-dimensional vector addition formula for Baker-Akhiezer functions obtained there characterizes Jacobians among principally polarized abelian varieties, and prove that it is indeed a weak characterization, i.e. that it is true up to additional components. We also show that this addition formula is equivalent to Gunning's multisecant formula for the Kummer variety obtained in [Gu2.
We then use the computation of the coefficients in the addition formula from [BK2to obtain cubic relations among theta functions that (weakly) characterize the locus of hyperelliptic Jacobians among irreducible abelian varieties. In genus 3 our equations are equivalent to the vanishing of one theta-null, and thus are classical (see [M, [P), but already for genus 4 they appear to be new.

1 Definitions and notations

We work over $\mathbb{C}$  , and fix the dimension/genus $g>1$  . Let ${\mathcal{ℋ}}_{g}$  be the Siegel upper half-space — the set of all $g×g$  period matrices $\tau$  , i.e.
symmetric complex $g×g$  matrices with positive definite imaginary part.
Each such $\tau$  corresponds to an abelian variety ${X}_{\tau }:={\mathbb{C}}^{g}/\tau {\mathbb{Z}}^{g}+{\mathbb{Z}}^{g}$  , and the moduli space ${\mathcal{A}}_{g}$  of principally polarized abelian varieties (ppavs) is then the quotient of ${\mathcal{ℋ}}_{g}$  by a certain action of the symplectic group $Sp\left(2g,\mathbb{Z}\right)$  .
A ppav is called irreducible if it is not isomorphic to a product of two lower-dimensional ppavs (with polarization). For convenience we denote by $\mathcal{A}irrg$  the moduli space of irreducible ppavs of genus $g$  . When in the following we say “abelian variety”, we actually mean a ppav.
Denoting $\mathbf{e}\left(x\right):=exp\left(\pi ix\right)$  , for a period matrix $\tau$  and a vector $z\in {\mathbb{C}}^{g}$  we define the theta function with characteristics $\varepsilon ,\delta \in \left(\mathbb{Z}/2\mathbb{Z}{\right)}^{g}$  , thought of as vectors consisting of zeros and ones, to be $\theta \left[\begin{array}{c}\varepsilon \\ \delta \end{array}\right]\left(\tau ,z\right):={\sum }_{m\in {\mathbb{Z}}^{g}}\mathbf{e}\left[\left(m+\frac{\varepsilon }{2},\tau \left(m+\frac{\varepsilon }{2}\right)\right)+2\left(m+\frac{\varepsilon }{2},z+\frac{\delta }{2}\right)\right],$  where $\left(\cdot ,\cdot \right)$  denotes the scalar product. A theta function with characteristics is even or odd as a function of $z$  depending on whether the scalar product $\left(\varepsilon ,\delta \right)$  is even or odd.
We denote by $\theta \left(\tau ,z\right):=\theta \left[\begin{array}{c}0\\ 0\end{array}\right]\left(\tau ,z\right)$  the classical Riemann's theta function. Theta functions with characteristics are, up to a constant factor, just the values of Riemann's theta function of a shifted argument:
 $\begin{array}{c}\begin{array}{c}\theta \left(\tau ,z+\frac{\tau \varepsilon +\delta }{2}\right)={\sum }_{m\in {\mathbb{Z}}^{g}}\mathbf{e}\left[\left(m,\tau m\right)+2\left(m,z+\frac{\tau \varepsilon +\delta }{2}\right)\right]\\ ={\sum }_{m\in {\mathbb{Z}}^{g}}\mathbf{e}\left[\left(m+\frac{\varepsilon }{2},\tau \left(m+\frac{\varepsilon }{2}\right)\right)+2\left(m+\frac{\varepsilon }{2},z+\frac{\delta }{2}\right)-\left(\frac{\varepsilon }{2},\tau \frac{\varepsilon }{2}\right)-\left(\varepsilon ,z+\frac{\delta }{2}\right)\right]\\ =\left(-1{\right)}^{\left(\varepsilon ,\delta \right)}\mathbf{e}\left[-\frac{1}{4}\left(\varepsilon ,\tau \varepsilon \right)-\left(\varepsilon ,z\right)\right]\theta \left[\begin{array}{c}\varepsilon \\ \delta \end{array}\right]\left(\tau ,z\right).\end{array}\end{array}$ (1)
Thus instead of thinking of a characteristic $\left[\begin{array}{c}\varepsilon \\ \delta \end{array}\right]$  as two integer vectors it sometimes is better to think of it as the point $\frac{\tau \varepsilon +\delta }{2}$  of order two on the abelian variety ${X}_{\tau }$  .
We further define theta functions of the second order to be $\Theta \left[\varepsilon \right]\left(\tau ,z\right):=\theta \left[\begin{array}{c}\varepsilon \\ 0\end{array}\right]\left(2\tau ,2z\right).$  For a fixed $\tau$  the theta functions, as functions of the variable $z$  , are sections of certain bundles on the abelian variety ${X}_{\tau }$  , which is to say that if the variable $z$  is translated by a vector of the lattice $\tau {\mathbb{Z}}^{g}+{\mathbb{Z}}^{g}$  , theta functions multiply by a certain number. In fact it is known that all theta functions of the second order are sections of the same bundle, denoted $2\Theta$  , and transform as follows:
 $\begin{array}{c}\Theta \left[\varepsilon \right]\left(\tau ,z+{e}_{j}+\tau {e}_{k}\right)=\mathbf{e}\left(-2\left({e}_{k},\tau {e}_{k}\right)-4\left({e}_{k},z\right)\right)\Theta \left[\varepsilon \right]\left(\tau ,z\right),\end{array}$ (2)
where we denote by ${e}_{k}$  the basis vector for the $k$  'th direction in ${\mathbb{C}}^{g}$  .
Theta functions of the second order form a basis for the sections of $2\Theta$  over ${X}_{\tau }$  . The square of any theta function with characteristics is also a section of $2\Theta$  , and thus is expressible as a linear combination of theta functions of the second order. In fact a slightly more general formula, Riemann's bilinear addition theorem, holds:
 $\begin{array}{c}\theta \left[\begin{array}{c}\varepsilon \\ \delta \end{array}\right]\left(z\right)\theta \left[\begin{array}{c}\varepsilon \\ \delta \end{array}\right]\left(w\right)={\sum }_{\sigma \in \left(\mathbb{Z}/2\mathbb{Z}{\right)}^{g}}\left(-1{\right)}^{\left(\delta ,\sigma \right)}\Theta \left[\sigma +\varepsilon \right]\left(\frac{z+w}{2}\right)\Theta \left[\sigma \right]\left(\frac{z-w}{2}\right).\end{array}$ (3)
For a fixed $\tau$  the map $z\to \left\{\Theta \left[\varepsilon \right]\left(\tau ,z\right){\right\}}_{all\varepsilon }$  defines the Kummer embedding $K:{X}_{\tau }/±1\to {\mathbb{P}}^{{2}^{g}-1}$  . This map is well-defined since all theta functions of the second order are sections of the same line bundle, and are even in $z$  .
The values of theta functions at $z=0$  are called the associated theta constants. Theta constants of the second order are modular forms of weight one half with respect to a certain finite index normal subgroup $\Gamma \left(2,4\right)\subset Sp\left(2g,\mathbb{Z}\right)$  , which is to say that if we act upon $\tau$  by some $\gamma \in \Gamma \left(2,4\right)$  , then $\Theta \left[\varepsilon \right]\left(\gamma \tau \right)=k\left(\gamma ,\tau \right)\Theta \left[\varepsilon \right]\left(\tau \right)$  , where $k$  is some multiplier depending on $\gamma$  and $\tau$  , but independent of $\varepsilon$  . Thus letting ${\mathcal{A}}_{g}\left(2,4\right):={\mathcal{ℋ}}_{g}/\Gamma \left(2,4\right)$  , we see that theta constants of the second order define a map $Th:{\mathcal{A}}_{g}\left(2,4\right)\to {\mathbb{P}}^{{2}^{g}-1}$  , which is known to be generically injective for all genera, and injective for $g\le 3$  . The level moduli space ${\mathcal{A}}_{g}\left(2,4\right)$  is a finite cover of ${\mathcal{A}}_{g}$  .
Let us denote by ${\mathcal{J}}_{g}\subset {\mathcal{A}}_{g}$  the locus of Jacobians of Riemann surfaces of genus $g$  , and by ${\mathcal{ℐ}}_{g}\subset {\mathcal{J}}_{g}$  the locus of Jacobians of hyperelliptic Riemann surfaces. The question of characterizing ${\mathcal{J}}_{g}$  within ${\mathcal{A}}_{g}$  is called the Schottky problem, and that of characterizing ${\mathcal{ℐ}}_{g}$  — the Schottky problem for the hyperelliptics. More precisely, one takes the preimages ${\mathcal{J}}_{g}\left(2,4\right)$  and ${\mathcal{ℐ}}_{g}\left(2,4\right)$  of ${\mathcal{J}}_{g}$  and ${\mathcal{ℐ}}_{g}$  , respectively, under the covering map $\pi :{\mathcal{A}}_{g}\left(2,4\right)\to {\mathcal{A}}_{g}$  , and asks to describe $Th\left({\mathcal{J}}_{g}\left(2,4\right)\right)$  and $Th\left({\mathcal{ℐ}}_{g}\left(2,4\right)\right)$  inside $Th\left({\mathcal{A}}_{g}\left(2,4\right)\right)$  . The question of describing $Th\left({\mathcal{A}}_{g}\left(2,4\right)\right)\subset {\mathbb{P}}^{{2}^{g}-1}$  , i.e. determining all the relations in the subring of the ring of modular forms generated by theta constants is also of interest, but we will not discuss it here. Notice that ${\mathcal{J}}_{g}$  and ${\mathcal{ℐ}}_{g}$  are irreducible, while ${\mathcal{J}}_{g}\left(2,4\right)$  and ${\mathcal{ℐ}}_{g}\left(2,4\right)$  have many irreducible components. We refer the reader to [Ifor more details on theta functions, and to [Grfor more details on the Schottky problem.
We will always think of a curve $C$  embedded in its Jacobian by the Abel-Jacobi map $A:C\to J\left(C\right)$  with some choice of the basis for the space of holomorphic differentials and of the starting point $P$  made.
This choice will be made explicitly when necessary. To avoid technical difficulties in the following sections, it will often be easier to work with the universal cover $\stackrel{~}{C}$  of a curve $C$  and the universal cover ${\mathbb{C}}^{g}$  of the abelian variety ${X}_{\tau }$  , and later take the automorphy properties of theta functions into account. The abelian variety will be fixed throughout, and thus we will often omit $\tau$  from the notations for theta functions and constants.
In this work we first prove the conjecture of Buchstaber and Krichever stated in [BK2that the validity of a certain $g$  -dimensional addition formula developed there and in [BK1characterizes Jacobians, but only up to additional components (i.e. that the locus of Jacobians is an irreducible component of the locus where the addition formula is satisfied), and then proceed to obtain from this some explicit identities for theta functions of hyperelliptic curves, using the explicit coefficients for the addition formula from [BK2.
In [M, theorem 9.1, and references therein Mumford showed that the hyperelliptic locus is characterized by a certain set of vanishing and non-vanishing conditions for theta constants with characteristics (the idea goes back at least to Thomae,see [T). In [PPoor showed that on $\mathcal{A}irrg$  Mumford's vanishing conditions by themselves (without the non-vanishing) define precisely the hyperelliptic locus, i.e. that there are no extra components. However, it is still not known how to obtain an ideal-theoretic description of the closure of $Th\left({\mathcal{ℐ}}_{g}\right)$  inside $Th\left({\mathcal{A}}_{g}\right)$  . It is known that if the vanishing holds and we also have some vanishing instead of non-vanishing, the abelian variety must be reducible, but then it does not necessarily have to be a limit of hyperelliptic Jacobians. Thus it would be interesting to study our equations on the reducible locus. It would also be very interesting to compare Mumford's equations to ours, but we have not been able to achieve this yet. We would also like to refer to a recent work [SMfor a further discussion of these issues as well as a description of components of ${\mathcal{ℐ}}_{g}^{2,4}$  as locally complete intersections.

Acknowledgements

We would like to thank Igor Krichever for bringing to our attention and explaining to us [BK1and [BK2and the conjectural characterization of Jacobians by the addition formula, which got this work started, Robert Gunning for helpful discussions on addition formulas for theta functions and multisecants, and Emma Previato for valuable discussions of the equations for Kummer and modular varieties.
We would like to especially thank Giuseppe Pareschi and Mihnea Popa, the authors of related recent preprint [PP(see remark  4 for details), for pointing out that the assumption of theorem 3 in the published version of this paper was too weak, and thus that some general position assumption was needed. The current text is the corrected version of the paper (the mistake was that lemma 4 was false), and the appropriate erratum is to appear in print shortly.
I owe a debt of gratitude to Professor Yum-Tong Siu for teaching me and sharing with me his insights and ideas in moduli theory. It is an honor for me to dedicate this paper to Professor Yum-Tong Siu on the occasion of his 60th birthday.

From [BK2we know the following addition formula (called “formula” to distinguish it from Riemann's addition “theorem”):
Fact 1 ([BK2, theorem 1). Let $P:={A}_{0},{A}_{1},\dots ,{A}_{g},Q:={A}_{g+1}\in \stackrel{~}{C}\subset {\mathbb{C}}^{g}$  , and denote by $R:=K-{\sum }_{i=1}^{g}{A}_{i}$  the vector of Riemann constants shifted by $-\sum {A}_{i}$  . Then for all $x,y\in {\mathbb{C}}^{g}$  the following identity is satisfied:
$\frac{\theta \left(Q+x+y+R\right)\theta \left(P+R\right)}{\theta \left(Q+R\right)\theta \left(P+x+y+R\right)}=$  $=\frac{\theta \left(Q+x+R\right)\theta \left(P+R\right)}{\theta \left(Q+R\right)\theta \left(P+x+R\right)}\cdot \frac{\theta \left(Q+y+R\right)\theta \left(P+R\right)}{\theta \left(Q+R\right)\theta \left(P+y+R\right)}$  $-{\sum }_{k=1}^{g}\frac{\theta \left(Q+R+{A}_{k}-P\right)\theta \left(Q+R-{A}_{k}+P+x+y\right)}{{\theta }^{2}\left(Q+R\right)\theta \left(P+x+y+R\right)\theta \left(2{A}_{k}-P+R\right)}$  $\cdot \frac{\theta \left(R+{A}_{k}+x\right)\theta \left(P+R\right)}{\theta \left(P+x+R\right)}\cdot \frac{\theta \left(R+{A}_{k}+y\right)\theta \left(P+R\right)}{\theta \left(P+y+R\right)}.$
Though this formula may look formidable, it is very explicit and is written entirely in terms of theta functions. In the following, we take the Abel-Jacobi map to start at $P$  , so that $P=0\in {\mathbb{C}}^{g}$  . Upon cancellations and multiplication by the common denominators, the above formula becomes simply $\theta \left(Q+x+y+R\right)\theta \left(Q+R\right)\theta \left(x+R\right)\theta \left(y+R\right)$  $=\theta \left(x+y+R\right)\theta \left(R\right)\theta \left(Q+x+R\right)\theta \left(Q+y+R\right)$  $-{\sum }_{k=1}^{g}\frac{\theta \left(R\right)\theta \left(Q+{A}_{k}+R\right)\theta \left(Q-{A}_{k}+x+y+R\right)\theta \left({A}_{k}+x+R\right)\theta \left({A}_{k}+y+R\right)}{\theta \left(2{A}_{k}+R\right)}.$  To see that this is in fact equivalent to Gunning's general multisecant formula from [Gu2(see Poor's work [Pfor an in-depth discussion) we use Riemann's bilinear addition theorem for the last two factors of each term. Denoting $z:=\frac{x+y}{2}$  and $w:=\frac{x-y}{2}$  , notice that the half-difference is always simply $\frac{x-y}{2}$  , so we will have a common factor of $\Theta \left[\sigma \right]\left(w\right)$  , and the resulting equation will be $0={\sum }_{\sigma }\left[\theta \left(Q+2z+R\right)\theta \left(Q+R\right)\Theta \left[\sigma \right]\left(z+R\right)-\theta \left(2z+R\right)\theta \left(R\right)\Theta \left[\sigma \right]\left(Q+z+R\right)$
 $\begin{array}{c}+{\sum }_{k=1}^{g}\frac{\theta \left(R\right)\theta \left(Q+{A}_{k}+R\right)\theta \left(Q-{A}_{k}+2z+R\right)\Theta \left[\sigma \right]\left({A}_{k}+z+R\right)}{\theta \left(2{A}_{k}+R\right)}\right]\Theta \left[\sigma \right]\left(w\right)\end{array}$ (4)
In the above the coefficient in the square brackets does not depend on $w$  , so we have an equation ${\sum }_{\sigma }{b}_{\sigma }\Theta \left[\sigma \right]\left(w\right)=0\forall w$  , where the coefficients ${b}_{\sigma }$  do not depend on $w$  . Since theta functions of the second order are a basis for sections of $2\Theta$  and thus linearly independent, it means that all coefficients ${b}_{\sigma }$  must be zero. Then since ${b}_{\sigma }$  actually are some functions of $z$  multiplied by $\Theta \left[\sigma \right]\left({A}_{i}+z+R\right)$  ( $i$  here ranges from $0$  to $g+1$  , i.e.
includes $P$  and $Q$  ), we have $\forall \sigma ,\forall z{\sum }_{i=0}^{g+1}{c}_{i}\left(z,{A}_{i}\right)\Theta \left[\sigma \right]\left({A}_{i}+z+R\right)=0$  for appropriate ${c}_{i}$  's.
Since $R$  does not depend on $z$  , in the above we can shift $z$  by $R$  and redefine the ${c}_{i}$  to see that the addition formula of [BK2implies the existence for any $z$  of some complex numbers ${c}_{i}$  , not all simultaneously zero, such that for some fixed $A$  's lying on the image $A\left(C\right)\subset J\left(C\right)$  we have
 $\begin{array}{c}{\sum }_{i=0}^{g+1}{c}_{i}\left(z\right)K\left({A}_{i}+z\right)=0\forall z\in {\mathbb{C}}^{g},\end{array}$ (5)
which is equivalent to saying that the $g+2$  points $K\left({A}_{i}+z\right)$  lie on a $g$  -plane in ${\mathbb{P}}^{{2}^{g}-1}$  . In fact Gunning proves a more general theorem (with a different and seemingly much more complicated expression for ${c}_{i}$  's in terms of the prime form):
Fact 2 ([Gu2, theorem 2). For any curve $C$  of genus $g$  , for any $1\le m\le g$  and for any points ${x}_{1},\dots ,{x}_{m},{A}_{0},\dots ,{A}_{m+1}\in A\left(C\right)$  the $m+2$  points $K\left({A}_{i}+\sum x-\sum A\right)$  are collinear, i.e. lie on the intersection of the Kummer variety of $C$  with an $m$  -plane in the projective space ${\mathbb{P}}^{{2}^{g}-1}$  .
In particular since the $g$  'th symmetric power of the curve is its Jacobian, ${S}^{g}\left(A\left(C\right)\right)=J\left(C\right)$  , the case $m=g$  of this theorem is formula ( 5 ), while the case $m=1$  is the case of a family of trisecant lines. It is shown in [Gu1that the existence of a family of trisecants characterizes Jacobians among irreducible ppavs.

3 Characterizing Jacobians by families of multisecants

We will now show that Buchstaber-Krichever's addition formula and Gunning's multisecant formula weakly characterize Jacobians.
Theorem 3. Let $X$  be an irreducible principally polarized abelian variety of dimension $g$  , and let ${A}_{0},\dots ,{A}_{g+1}$  be distinct points of $X$  . Suppose that $\forall z\in X$  the $g+2$  points $K\left({A}_{i}+z\right)$  are linearly dependent. Assume moreover the following general position condition: that there exist some $k$  and $l$  such that for $y:=-\frac{{A}_{k}+{A}_{l}}{2}$  the linear span of the points $K\left({A}_{i}+y\right)$  for $i=0\dots g+1$  is of dimension precisely $g+1$  , and not less.
Then $X$  is the Jacobian of some curve $C$  , and all ${A}_{i}\in A\left(C\right)$  .
• Proof. Assume $g\ge 4$  , otherwise the theorem is trivial. Working in the spirit of Gunning's work [Gu1, we reduce the theorem to the case of the trisecant.
Since the rank of the $\left(g+2\right)×{2}^{g}$  matrix $K\left({A}_{i}+y\right)$  is equal to exactly $g+1$  , it means that there exists unique $\stackrel{⃗}{c}\left(y\right)\in {\mathbb{P}}^{g+1}$  such that $\sum {c}_{i}\left(y\right)K\left({A}_{i}+y\right)=0$  . Moreover, since we have ${A}_{k}+y=-\left({A}_{l}+y\right)$  and thus $K\left({A}_{k}+y\right)=K\left({A}_{l}+y\right)$  , we must then have $p:=\stackrel{⃗}{c}\left(y\right)={1}_{k}-{1}_{l}$  , where ${1}_{i}$  denotes the basis vector of ${\mathbb{C}}^{g}$  in the $i$  'th direction. Since the rank of $K\left({A}_{i}+y\right)$  is equal to $g+1$  , the rank of $K\left({A}_{i}+z\right)$  must be equal to $g+1$  identically in a neighborhood of $y$  , and thus for all $z$  sufficiently close to $y$  there is a unique projective solution $\stackrel{⃗}{c}\left(z\right)$  to $\sum {c}_{i}\left(z\right)K\left({A}_{i}+z\right)=0$  .
We will show that the differential $d\stackrel{⃗}{c}:{T}_{y}{\mathbb{C}}^{g}\to {T}_{p}{\mathbb{P}}^{g+1}$  is of maximal rank. This will then imply that locally near $p$  the image $\stackrel{⃗}{c}\left({\mathbb{C}}^{g}\right)$  is of dimension $g$  and thus locally the preimage of any coordinate plane ${\mathbb{P}}^{2}\subset {\mathbb{P}}^{g+1}$  containing $p$  is at least one-dimensional. Then we have a one-dimensional family of trisecants of the Kummer variety, and by Welters' [Weinfinitesimal version of Gunning's trisecant criterion from [Gu1, $X$  is a Jacobian of some curve $C$  with the points ${A}_{i},{A}_{k},{A}_{l}$  lying on $A\left(C\right)$  . Thus we see that all ${A}_{i}$  lie on $A\left(C\right)$  .
For contradiction, suppose that the rank of $d\stackrel{⃗}{c}{|}_{y}$  is not maximal, i.e.
that there is some vector $v\in {\mathbb{C}}^{g}$  such that we have $\frac{\partial }{\partial v}{c}_{i}\left(z\right){|}_{y}=\lambda {c}_{i}\left(y\right)$  for some constant $\lambda$  independent of $i$  — this means that the derivative $\frac{\partial }{\partial v}$  of the projective point $\stackrel{⃗}{c}$  is zero.
Let us take the derivative $\frac{\partial }{\partial v}$  of ( 5 ) at $y$  :  $\begin{array}{c}{\sum }_{i=0}^{g+1}{\frac{\partial {c}_{i}\left(z\right)}{\partial v}K\left({A}_{i}+z\right)+{c}_{i}\left(z\right)\frac{\partial K}{\partial v}\left({A}_{i}+z\right)|}_{y}\end{array}$
 $\begin{array}{c}=\lambda \left(K\left({A}_{k}+y\right)-K\left({A}_{l}+y\right)\right)+\left(\frac{\partial K}{\partial v}\left({A}_{k}+y\right)-\frac{\partial K}{\partial v}\left({A}_{l}+y\right)\right)\end{array}$
 $\begin{array}{c}=2\frac{\partial K}{\partial v}\left({A}_{k}+y\right)=0,\end{array}$
because ${A}_{k}+y=-\left({A}_{l}+y\right)$  , theta functions are even and their derivatives are odd. But this then implies that $\frac{\partial K}{\partial v}\left({A}_{k}+y\right)=0$  , which is impossible by Wirtinger's theorem unless ${A}_{k}+y=\frac{{A}_{k}-{A}_{l}}{2}$  is a point of order two in $X$  . If this is the case, though, it would mean that ${A}_{k}-{A}_{l}=0\in X$  , which contradicts the assumption that all ${A}_{i}$  are distinct. Thus we have arrived at a contradiction and showed that the differential $d\stackrel{⃗}{c}:{T}_{y}{\mathbb{C}}^{g}\to {T}_{P}{\mathbb{P}}^{g+1}$  is injective.
Remark 4. In their recent preprint [PPPareschi and Popa prove the following statement.
Castelnuovo-Schottky lemma ([PP). Let ${A}_{0}\dots {A}_{n}\in X$  be a set of points in $\theta$  -general position imposing only $g+1$  conditions on $|2\Theta {|}_{a}$  for $a$  general. Then if $n\ge g+1$  , $X$  is the Jacobian of some curve $C$  , and moreover ${A}_{i}\in A\left(C\right)$  .
All sections of $|2\Theta {|}_{a}$  are scalar products $v\cdot K\left(z+a\right)$  for some ${2}^{g}$  -dimensional vector $v$  , and the condition for a section to vanish at ${A}_{i}$  is simply that $v\cdot K\left({A}_{i}+a\right)=0$  . Thus the points imposing only $g+1$  conditions on $|2\Theta {|}_{a}$  means that for general $a$  the rank of the $\left(n+1\right)×{2}^{g}$  matrix $K\left({A}_{i}+a\right)$  is equal to $g+1$  , which is the collinearity condition of our theorem 3.
The $\theta$  -general position means that for any $g+1$  points ${p}_{1}\dots {p}_{g},{p}_{g+1}$  among $\left\{{A}_{i}\right\}$  there exists some translate $a\in X$  such that $\theta \left({p}_{i}+a\right)=0\ne \theta \left(q+a\right)$  — this is different from our rank condition.
Since both our theorem 3 and the Castelnuovo-Schottky lemma characterize Jacobians, these assumptions must be equivalent, but we do not know a way to derive one from the other directly.
With more work it can also be shown that Gunning's addition formula for all other values of $m$  also serves to characterize Jacobians:
Proposition 5. Let $X$  be an irreducible principally polarized abelian variety of dimension $g$  , and let ${A}_{0},\dots ,{A}_{m+1}$  be different points of $X$  . Suppose that the $m+2$  points $K\left({A}_{i}+z\right)$  are linearly dependent for any $z\in M$  , where the set $M\subset X$  is at least $m$  -dimensional at the point $y=-\frac{{A}_{k}+{A}_{k}}{2}$  for some $k$  and $jk$  (i.e. has a non-degenerate $m$  -jet at $y$  ). Assume moreover that the rank of the $\left(m+2\right)×{2}^{g}$  matrix $K\left({A}_{i}+y\right)$  is equal to precisely $m+1$  . Then $X$  is the Jacobian of some curve.
• Proof. We again study the local situation near $y$  and work with the germ of $M$  at this point. Then we want to show that $d\stackrel{⃗}{c}$  is injective on ${T}_{y}\stackrel{~}{M}$  . But since $d\stackrel{⃗}{c}$  is non-degenerate on ${T}_{y}{\mathbb{C}}^{g}$  , it is also non-degenerate on a subspace, and we are done.
Remark 6. By imitating Gunning's proof of the trisecant theorem in [Gu1it seems to be also possible to show that if $M$  is $m$  -dimensional at some point, not necessarily $-\frac{{A}_{i}+{A}_{j}}{2}$  and under some general position assumption, then $X$  a Jacobian. However, as we do not need this result in what follows, and the proof would be technically rather complicated, we will not give it.
In [Gu2Gunning shows that for a Jacobian $X$  with fixed $A$  's the coefficients ${c}_{i}$  in ( 5 ) are unique up to scaling and expressible in terms of the Klein-Gunning prime form. The Klein-Gunning prime form expression is rather hard to deal with (see [Pfor a detailed discussion and computations); however, the expression for ${c}_{i}$  obtained in formula ( 4 ) seems more amenable.
Equation ( 4 ) includes, however, both theta functions of the second order and the classical Riemann's theta function with zero characteristic.
Let us use the addition theorem once again for the last two factors of the type $\theta \left(\dots \right)\theta \left(\dots \right)$  in each term. We then get for all $\sigma$
 $\begin{array}{c}\begin{array}{c}{\sum }_{\varepsilon }\Theta \left[\varepsilon \right]\left(Q+z+R\right)\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\sigma \right]\left(z+R\right)-\Theta \left[\varepsilon \right]\left(z+R\right)\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\sigma \right]\left(Q+z+R\right)\\ +{\sum }_{k,\varepsilon }\frac{\theta \left(R\right)}{\theta \left(2{A}_{k}+R\right)}\Theta \left[\varepsilon \right]\left(Q+z+R\right)\Theta \left[\varepsilon \right]\left(z-{A}_{k}\right)\Theta \left[\sigma \right]\left({A}_{k}+z+R\right)=0.\end{array}\end{array}$ (6)
Since all the transformations that we have done so far are equivalencies, we do not lose any information, unless all the coefficients ${c}_{i}$  in the formula above are identically zero — in which case the projective point $\stackrel{⃗}{c}\left(z\right)$  is never defined.
Since the coefficients in ( 5 ) are unique, they must be exactly the ones given by formula ( 6 ), and thus formula ( 6 ) is satisfied if and only if the collinearity condition ( 5 ) is satisfied.
Proposition 7. With the same assumptions and notations as above, identity ( 6 ) characterizes Jacobians among all irreducible abelian varieties.

4 Addition formula for the hyperelliptic case

Formula ( 6 ) we obtained is in terms of theta functions evaluated at different points. In the classical approach to the Schottky problem (see [Sfor the origins and [F, [Grfor a review) one wants to characterize the Jacobian locus by some algebraic relations among theta constants.
Thus it would be nice if for some special values of $z,P,Q$  and $A$  's the addition formula yielded such equations.
However, from the transformation rule ( 2 ) for theta functions it is easy to see that if a vector $v$  is not a point of order two, then $\Theta \left[\varepsilon \right]\left(z+v\right)$  is not a section of the bundle $2\Theta$  and thus cannot be expressed as a linear combination of theta functions of the second order. Thus we only have a reasonable hope of getting from ( 6 ) some equations for theta constants if we are so lucky that all the “shifts” of $z$  that appear there are points of order two. In particular, this means that the points ${A}_{i}+z+R$  must be of order two for all $i$  .
Now suppose that indeed both ${A}_{i}+z+R$  and ${A}_{j}+z+R$  for $i\ne j$  are points of order two on $J\left(C\right)$  . Then their difference, ${A}_{i}-{A}_{j}$  , is also of order two, so we have $2{A}_{i}-2{A}_{j}=0\in J\left(C\right)$  . By Abel's theorem this then means that there is a function $f$  on $C$  whose divisor is equal to $2{A}_{i}-2{A}_{j}$  , i.e. with a double pole at ${A}_{j}$  and holomorphic on $C-\left\{{A}_{j}\right\}$  . Since the existence of such a function characterizes hyperelliptic curves, it means that $C$  then has to be hyperelliptic. For the hyperelliptic curves it is known (see, for example, [M) that if we take $P=0$  to be the image of one of the Weierstrass points, then the other $2g+1$  Weierstrass points will also map to points of order two on the Jacobian. Thus let us assume that all $A$  's and $Q$  in formula  6 are chosen to be points of order two, i.e. that we are dealing with a hyperelliptic curve, and let us rewrite the addition formula in this case. Denote $Q=\frac{\tau {\alpha }_{0}+{\beta }_{0}}{2}$  , ${A}_{k}=\frac{\tau {\alpha }_{k}+{\beta }_{k}}{2}$  and $R=\frac{\tau \alpha +\beta }{2}$  — in fact $R$  is expressible in terms of $A$  's and Riemann constants, but we will deal with this later.
Now we rewrite formula ( 6 ) for these $A$  's and $Q$  . In doing this, we need to be extra careful to remember that we are actually working on $\stackrel{~}{C}$  and ${\mathbb{C}}^{g}$  , as not to omit any important automorphy factors. Indeed, from the automorphy properties of $\theta$  it follows that $\frac{\theta \left(R\right)}{\theta \left(2{A}_{k}+R\right)}=\frac{\theta \left(\frac{\tau \alpha +\beta }{2}\right)}{\theta \left(\tau {\alpha }_{k}+{\beta }_{k}+\frac{\tau \alpha +\beta }{2}\right)}=\mathbf{e}\left[\left({\alpha }_{k},\tau {\alpha }_{k}\right)+\left(\alpha ,\tau {\alpha }_{k}\right)\right].$  Also for any integers $a$  and $b$  it follows from ( 1 ) and ( 2 ) that
 $\begin{array}{c}\begin{array}{c}\Theta \left[\delta \right]\left(\tau ,z+\frac{\tau a+b}{2}\right)=\theta \left[\begin{array}{c}\delta \\ 0\end{array}\right]\left(2\tau ,2z+\tau a+b\right)\\ =\theta \left(2\tau ,2z+\tau \delta +\tau a\right)\mathbf{e}\left[\frac{1}{2}\left(\delta ,\tau \delta \right)+2\left(\delta ,z\right)+\left(\delta ,\tau a\right)+\left(\delta ,b\right)\right]\\ =\theta \left[\begin{array}{c}\delta +a\\ 0\end{array}\right]\left(2\tau ,2z\right)\mathbf{e}\left[-\frac{1}{2}\left(\delta +a,\tau \left(\delta +a\right)\right)-2\left(\delta +a,z\right)\right]\\ \cdot \mathbf{e}\left[\frac{1}{2}\left(\delta ,\tau \delta \right)+2\left(\delta ,z\right)+\left(\delta ,\tau a\right)+\left(\delta ,b\right)\right]\\ =\left(-1{\right)}^{\left(\delta ,b\right)}\mathbf{e}\left[-\frac{1}{2}\left(a,\tau a\right)-2\left(a,z\right)\right]\Theta \left[\delta +a\right]\left(\tau ,z\right).\end{array}\end{array}$ (7)
When we substitute this into ( 6 ) notice that as functions of $z$  all terms are actually sections of the same bundle, $6\Theta$  , as each is cubic in theta functions of the second order. Thus the $\mathbf{e}\left[\left(*,z\right)\right]$  factors must cancel everywhere. Also evaluating at $z=0$  and noticing that all terms are modular forms in $\tau$  with respect to $\Gamma \left(2,4\right)$  of the same weight, we expect the factors $\mathbf{e}\left[\left(*,\tau *\right)\right]$  to cancel as well. A trivial but tedious computation confirms this, and we arrive at
Proposition 8. An irreducible abelian variety ${X}_{\tau }$  with some points $P=0,Q,{A}_{1},\dots ,{A}_{g}$  with ${A}_{i}=\frac{\tau {\alpha }_{k}+{\beta }_{k}}{2}\in {\mathbb{C}}^{g}$  is the Jacobian of a hyperelliptic curve $C$  , and $P,Q,{A}_{i}\in A\left(\stackrel{~}{C}\right)\subset {\mathbb{C}}^{g}$  if and only if the following is satisfied for all $\sigma \in \left(\mathbb{Z}/2\mathbb{Z}{\right)}^{g}$  and for all $z\in {\mathbb{C}}^{g}$  :
 $\begin{array}{c}\begin{array}{cc}& {\sum }_{\varepsilon }\left(-1{\right)}^{\left(\varepsilon ,\beta +{\beta }_{0}\right)}\Theta \left[\varepsilon +\alpha +{\alpha }_{0}\right]\left(z\right)\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\sigma +\alpha \right]\left(z\right)\\ =& {\sum }_{e,k}\left(-1{\right)}^{\left(\varepsilon ,\beta +{\beta }_{0}+{\beta }_{k}\right)+\left(\sigma ,{\beta }_{k}\right)}\Theta \left[\varepsilon +\alpha +{\alpha }_{0}\right]\left(z\right)\Theta \left[\varepsilon +{\alpha }_{k}\right]\left(z\right)\Theta \left[\sigma +\alpha +{\alpha }_{k}\right]\left(z\right)\\ +& {\sum }_{\varepsilon }\left(-1{\right)}^{\left(\varepsilon ,\beta \right)+\left(\sigma ,{\beta }_{0}\right)}\Theta \left[\varepsilon +\alpha \right]\left(z\right)\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\sigma +\alpha +{\alpha }_{0}\right]\left(z\right).\end{array}\end{array}$ (8)
and not all the coefficients in front of the Kummer images (i.e. in front of the last theta function factor, for $\sigma$  varying) appearing here are identically zero in $z$  .

5 Cubic equations for the hyperelliptic locus

To make formula ( 8 ) entirely explicit, we now need to pick some specific way to map a hyperelliptic curve into its Jacobian, and pick some $g+2$  Weierstrass points on it in a certain way. This is indeed a very classical construction.
Let us think of a hyperelliptic curve sitting on a skewer that intersects it in precisely the $2g+2$  Weierstrass points. Label them ${p}_{1},\dots ,{p}_{2g+2}$  going from left to right along the skewer. Then pick for the basis of the cycles ${a}_{i}$  to be the loop around the $i$  'th handle, passing through points ${p}_{2i-1}$  and ${p}_{2i}$  , and ${b}_{i}$  to be the loop around the $i$  'th hole, passing through ${p}_{2i}$  and ${p}_{2i+1}$  . Then thinking of the skewer as being the $x$  axis and the whole picture being that of ${y}^{2}=\prod \left(x-{p}_{i}\right)$  , we can compute the images of ${p}_{i}$  in the Jacobian. Indeed, let us use ${p}_{1}$  as the starting point, so that $A\left({p}_{1}\right)=0$  . Then we see that $A\left({p}_{2}\right)=\frac{{e}_{1}}{2}$  , $A\left({p}_{3}\right)=\frac{\tau {e}_{1}+{e}_{1}}{2}$  , $A\left({p}_{4}\right)=\frac{\tau {e}_{1}+{e}_{1}}{2}+\frac{{e}_{1}+{e}_{2}}{2}=\frac{\tau {e}_{1}+{e}_{2}}{2}+{e}_{1}$  , $A\left({p}_{5}\right)=\frac{\tau \left({e}_{1}+{e}_{2}\right)+{e}_{2}}{2}+{e}_{1}$  , and in general we have $A\left({p}_{2i}\right)=\frac{\tau {s}_{i-1}+{e}_{i}}{2}+{s}_{i-1}$  and $A\left({p}_{2i+1}\right)=\frac{\tau {s}_{i}+{e}_{i}}{2}+{s}_{i-1}$  for $1\le i\le g$  , while $A\left({p}_{2g+2}\right)=\frac{\tau {s}_{g}}{2}$  (where for convenience we have denoted ${s}_{k}:={\sum }_{i=1}^{k}{e}_{i}$  ).
For our purposes let us choose $P:=0=A\left({p}_{1}\right)$  , $Q:=A\left({p}_{2}\right)$  and ${A}_{k}:=A\left({p}_{2k+2}\right)$  for $1\le k\le g$  . We now need to compute the vector $R$  , i.e to compute the vector of Riemann constants and subtract from it the sum of $A$  's. The result is certainly classical: $R=Q=A\left({p}_{2}\right)$  .
To prove this one can note that by definition $R$  is the unique vector such that $\theta \left(A\left(p\right)+R\right)$  as a function of $p\in C$  has precisely $g$  zeroes at ${A}^{-1}\left({A}_{k}\right)$  , i.e. at ${p}_{2k+2}$  . To check that this is the case we note that $A\left({p}_{2}\right)+A\left({p}_{2i}\right)$  is always odd, as a theta characteristic, so that even Riemann's theta function will vanish at the point $A\left({p}_{2}\right)+A\left({p}_{2i}\right)$  , while $A\left({p}_{2}\right)+A\left({P}_{2i+1}\right)$  is even; thus $R=A\left({p}_{2}\right)$  .
Let us now substitute all this into formula ( 8 ). We have $\alpha ={\alpha }_{0}=0$  , $\beta ={\beta }_{0}={e}_{1}$  , ${\alpha }_{k}={s}_{k}$  and ${\beta }_{k}={e}_{k+1}$  , where we understand ${e}_{g+1}$  to be zero. Since all theta functions of the second order are periodic with respect to $z\to z+{e}_{i}$  , the additional integer shifts by ${s}_{i}$  do not matter, and finally ( 8 ) yields
Theorem 9. An irreducible period matrix $\tau \in {\mathcal{ℋ}}_{g}$  is the period matrix of a hyperelliptic Jacobian with the basis of cycles chosen as above if and only if the following cubic identity for theta functions of the second order is satisfied for all $\sigma \in \left(\mathbb{Z}/2\mathbb{Z}{\right)}^{g}$  and for all $z\in {X}_{\tau }$  (and thus for all $z\in {\mathbb{C}}^{g}$  ):
 $\begin{array}{c}\begin{array}{cc}& {\sum }_{\varepsilon }\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\sigma \right]\left(z\right)\\ =& {\sum }_{\varepsilon }{\sum }_{k=0}^{g}\left(-1{\right)}^{\left(\varepsilon +\sigma ,{e}_{k+1}\right)}\Theta \left[\varepsilon \right]\left(z\right)\Theta \left[\varepsilon +{s}_{k}\right]\Theta \left[\sigma +{s}_{k}\right]\left(z\right),\end{array}\end{array}$ (9)
where we understand ${e}_{g+1}$  to be zero, and asumme moreover that not all the coefficients appearing in front of $\Theta \left[\sigma +{s}_{k}\right]\left(z\right)$  and $\Theta \left[\sigma \right]\left(z\right)$  are identically zero in $z$  .
To check that this makes sense let us do the computations in low genus and see what we get. To simplify formulas, we write $\left[\varepsilon \right]$  for $\Theta \left[\varepsilon \right]\left(z\right)$  . We order the $\varepsilon$  for summation of the terms of ( 9 ) lexicographically to keep track of where we are.
Genus 2: We do not expect to get any meaningful equations, as any irreducible abelian variety of dimension two is a hyperelliptic Jacobian, so our characterization should be vacuous. We verify this; here is what formula ( 9 ) yields for $\sigma =00$  :
 $\begin{array}{cc}& \left[00\right]\left[00\right]\left[00\right]+\left[01\right]\left[01\right]\left[00\right]+\left[10\right]\left[10\right]\left[00\right]+\left[11\right]\left[11\right]\left[00\right]\end{array}$
 $\begin{array}{cc}=& \left[00\right]\left[00\right]\left[00\right]+\left[01\right]\left[01\right]\left[00\right]-\left[10\right]\left[10\right]\left[00\right]-\left[11\right]\left[11\right]\left[00\right]\end{array}$
 $\begin{array}{cc}+& \left[00\right]\left[10\right]\left[10\right]-\left[01\right]\left[11\right]\left[10\right]+\left[10\right]\left[00\right]\left[10\right]-\left[11\right]\left[01\right]\left[10\right]\end{array}$
 $\begin{array}{cc}+& \left[00\right]\left[11\right]\left[11\right]+\left[01\right]\left[10\right]\left[11\right]+\left[10\right]\left[01\right]\left[11\right]+\left[11\right]\left[00\right]\left[11\right]\end{array}$
and all the terms cancel. A similar computation shows that the identity is also trivial for all other choices of $\sigma$  .
Genus 3: Here we have $dim{\mathcal{J}}_{3}=dim{\mathcal{A}}_{3}=6$  , while $dim{\mathcal{ℐ}}_{3}=5$  , so we should have a non-trivial identity. Indeed let us choose $\sigma =000$  and write down ( 9 ) in this case; after multiple cancellations and dividing by two it becomes simply
 $\begin{array}{cc}& \Theta \left[000\right]\left(z\right)\Theta \left[101\right]\left(z\right)\Theta \left[101\right]\left(z\right)+\Theta \left[011\right]\left(z\right)\Theta \left[101\right]\left(z\right)\Theta \left[110\right]\left(z\right)\end{array}$
 $\begin{array}{cc}=& \Theta \left[010\right]\left(z\right)\Theta \left[101\right]\left(z\right)\Theta \left[111\right]\left(z\right)+\Theta \left[001\right]\left(z\right)\Theta \left[100\right]\left(z\right)\Theta \left[101\right]\left(z\right),\end{array}$
which using formula ( 3 ) is equivalent to
 $\begin{array}{c}\Theta \left[101\right]\left(z\right)\cdot \theta \left[\begin{array}{c}101\\ 111\end{array}\right]\left(2z\right)\cdot \theta \left[\begin{array}{c}101\\ 111\end{array}\right]\left(0\right)=0.\end{array}$ (10)
Choosing a different $\sigma$  yields a different equation: in general we get $\Theta \left[101+\sigma \right]\left(z\right)\cdot \theta \left[\begin{array}{c}101\\ 111\end{array}\right]\left(2z\right)\cdot \theta \left[\begin{array}{c}101\\ 111\end{array}\right]\left(0\right)$  . All of these equations together are equivalent to $\theta \left[\begin{array}{c}101\\ 111\end{array}\right]\left(0\right)=0$  , since theta functions of the second order never all vanish simultaneously and $\theta \left[\begin{array}{c}101\\ 111\end{array}\right]\left(2z\right)$  cannot be identically zero in $z$  . Now to actually characterize ${\mathcal{ℐ}}_{3}\subset {\mathcal{A}}_{3}$  we need to get rid of the condition “that the basis of cycles is chosen as above” in theorem  9 . But choosing a different basis of cycles means acting on the period matrix by a symplectic transformation. Since theta constants are modular with respect to $\Gamma \left(2,4\right)$  , conjugating equation ( 10 ) by any $\gamma \in \Gamma \left(2,4\right)$  would not change it, so we only need to act by the finite group $Sp\left(2g,\mathbb{Z}\right)/\Gamma \left(2,4\right)$  . It is well known that the action of this group is transitive on the set of even theta characteristics (see [I), so we can get the vanishing of a theta constant with any even characteristic.
Proposition 10. An irreducible abelian variety of genus 3 is a hyperelliptic Jacobian if it has a vanishing theta constant with even characteristic. This is known classically, see [M.
Genus 4: here the situation is more interesting: Mumford's conditions include some non-vanishing, so getting the explicit equations for the closure of $Th\left({\mathcal{ℐ}}_{4}\right)\subset Th\left({\mathcal{A}}_{4}\right)$  , without any inequalities serving to cut off the extra components inside the reducible locus, would be interesting.
For the case of $\sigma =0000$  the equation we get from ( 9 ) after cancellations becomes the following cubic (we have rearranged the terms lexicographically and omitted square brackets):
 $\begin{array}{c}0000\cdot 1001\cdot 1001+0000\cdot 1010\cdot 1010+0000\cdot 1011\cdot 1011+0000\cdot 1101\cdot 1101\end{array}$
 $\begin{array}{c}0011\cdot 1101\cdot 1110+0101\cdot 1000\cdot 1101+0101\cdot 1011\cdot 1110+0110\cdot 1010\cdot 1100\end{array}$
 $\begin{array}{c}0111\cdot 1001\cdot 1110+0111\cdot 1011\cdot 1100=0001\cdot 1000\cdot 1001+0001\cdot 1100\cdot 1101\end{array}$
 $\begin{array}{c}0010\cdot 1000\cdot 1010+0010\cdot 1101\cdot 1111+0011\cdot 1000\cdot 1011+0100\cdot 1010\cdot 1110\end{array}$
 $\begin{array}{c}0100\cdot 1011\cdot 1111+0101\cdot 1001\cdot 1100+0101\cdot 1010\cdot 1111+0110\cdot 1001\cdot 1111,\end{array}$
while for example for $\sigma =0001$  we get
 $\begin{array}{c}0000\cdot 1000\cdot 1001+0000\cdot 1100\cdot 1101+0010\cdot 1001\cdot 1010+0011\cdot 1001\cdot 1011\end{array}$
 $\begin{array}{c}0011\cdot 1100\cdot 1110+0100\cdot 1000\cdot 1101+0100\cdot 1011\cdot 1110+0101\cdot 1010\cdot 1110\end{array}$
 $\begin{array}{c}0101\cdot 1011\cdot 1111+0111\cdot 1000\cdot 1110=0001\cdot 1000\cdot 1000+0001\cdot 1010\cdot 1010\end{array}$
 $\begin{array}{c}0001\cdot 1011\cdot 1011+0001\cdot 1100\cdot 1100+0010\cdot 1100\cdot 1111+0100\cdot 1001\cdot 1100\end{array}$
 $\begin{array}{c}0100\cdot 1010\cdot 1111+0110\cdot 1000\cdot 1111+0110\cdot 1010\cdot 1101+0111\cdot 1011\cdot 1101.\end{array}$
Neither of these cubics is equal to $\Theta \left[\delta \right]\left(z\right)\theta \left[\begin{array}{c}\alpha \\ \beta \end{array}\right]\left(2z\right)\theta \left[\begin{array}{c}\alpha \\ \beta \end{array}\right]\left(0\right)$  for any $\alpha ,\beta ,\delta$  . However, from theorem  9 we see that
Proposition 11. The vanishing of the full set of $16$  cubics similar to the ones above, for all $\sigma$  , identically in $z$  characterizes a component of $Th\left({\mathcal{ℐ}}_{4}\left(2,4\right)\right)\subset Th\left({\mathcal{A}}_{4}^{irr}\left(2,4\right)\right)$  .
By Mumford's and Poor's results such a component is also determined by identical vanishing of some set of theta constants with characteristics.
Thus the vanishing of our 16 cubics should imply, for irreducible abelian varieties, the vanishing of some theta constants with characteristics and vice versa, but we are now unable to see this directly.
The difficulty in doing so is not only due to the fact that the cubic equations are very complicated. Indeed, thinking of each cubic ${f}_{\sigma }$  , evaluated at $z=0$  , as a polynomial on ${\mathbb{P}}^{15}$  with zero locus $Z\left({f}_{\sigma }\right)$  , we can only say that ${\cap }_{\sigma }Z\left({f}_{\sigma }\right)\cap Th\left({\mathcal{A}}_{4}^{irr}\left(2,4\right)\right)$  is contained in the common zero locus of some quadrics (which are by ( 3 ) the expressions for theta constants with characteristics in terms of theta constants of the second order) on $Th\left({\mathcal{A}}_{4}^{irr}\left(2,4\right)\right)$  . It may in fact not be the case that the whole $\cap Z\left({f}_{\sigma }\right)$  is contained in the zero locus of these quadrics in ${\mathbb{P}}^{15}$  . Thus to be able to see the relation of the vanishing of our cubics to the vanishing of theta constants with characteristics, we may need to know the equations for the closure of $Th\left({\mathcal{A}}_{4}^{irr}\left(2,4\right)\right)\subset {\mathbb{P}}^{15}$  , which are not known.
The above discussion was for just one component of $Th\left({\mathcal{ℐ}}_{4}\left(2,4\right)\right)$  , which projects to just one component of ${\mathcal{ℐ}}_{4}\subset {\mathcal{A}}_{4}$  . The equations for the other components corresponding to different choices of the basis of cycles on the curve are of course obtained by acting on the set of 16 cubics by elements of $Sp\left(8,\mathbb{Z}\right)/\Gamma \left(2,4\right)$  .
Final remark. It seems likely that in any genus evaluating equations ( 9 ) only at $z=0$  for all $\sigma$  should yield the defining set of equations for a component of $Th\left({\mathcal{ℐ}}_{g}\left(2,4\right)\right)$  . Indeed using ( 7 ) it can be easily shown that if these are satisfied, then ( 9 ) is satisfied for $z$  being any point of order two. Thus both sides of ( 9 ) are sections of $6\Theta$  that agree at all points of order two, and one would hope that then they agree everywhere and give the same function of $z$  , so that ( 9 ) is true identically and we can apply theorem  9 . References

1. Buchstaber, V.M., Krichever, I.M.: Vector addition theorems and Baker-Akhiezer functions. Teor. i Mat. Fizika 94 (1993), no. 2, 200–212.
2. Buchstaber, V.M., Krichever, I.M.: Multidimensional vector addition theorems and the Riemann theta functions. Internat. Math. Res. Notices 1996, no. 10, 505–513.
3. Igusa, J.-I.: Theta functions. Die Grundlehren der mathematischen Wissenschaften, Band 194. Springer-Verlag, New York-Heidelberg, 1972.
4. Farkas, H.: Schottky-Jung theory, Theta functions — Bowdoin 1987, 459–483
5. Grushevsky, S.: Effective Schottky problem, Ph.D. Dissertation, Harvard University, 2002
6. Gunning, R. C.: Some curves in abelian varieties. Invent. Math. 66 (1982), no. 3, 377–389.
7. Gunning, R. C.: Some identities for abelian integrals. Amer. J. Math. 108 (1986), no. 1, 39–74.
8. Mumford, D.: Tata lectures on Theta II. Progress in Math. vol. 43, Birkhäuser, Boston, 1984
9. Poor, C.: The hyperelliptic locus, Duke Math. J. 76 (1994) 3, 809–884
10. Pareschi, G., Popa, M. Castelnuovo theory and the geometric Schottky problem, preprint math.AG/0407370
11. Salvati Manni, R.: Modular forms vanishing on the hyperelliptic locus, Japan. J. Math. (N.S.) 29 (2003), no. 1, 135–142
12. Schottky, F.: Zur Theorie der Abelschen Functionen vor vier Variablen. J. Reine Angew. Math. 102 (1888), 304–352
13. Thomae, J.: Beitrag zur Bestimmung von $\theta \left(0,0,\dots ,0\right)$  durch die Klassenmoduln algebraischer Funktionen, J. Reine Angew. Math. 71 (1870), 201–222
14. Welters, G.: A criterion for Jacobi varieties. Ann. of Math. (2) 120 (1984), no. 3, 497–504
15. Wirtinger, W.: Untersuchungen über Thetafunktionen. Leipzig: Teubner 1895

Mathematics Department, Princeton University, Fine Hall, Washington Road, Princeton, NJ 08544, USA E-mail address : sam@math.princeton.edu