We consider the linear system $|2{\Theta }_0|$ of second order theta functions over the Jacobian $JC$ of a non-hyperelliptic curve $C$. A result by J.Fay says that a divisor $D\in |2{\Theta }_0|$ contains the origin $𝒪\in JC$ with multiplicity $4$ if and only if $D$ contains the surface $C-C=\left\{𝒪\right(p-q)\mid p,q\in C\}\subset JC$. In this paper we generalize Fay’s result and some previous work by R.C.Gunning. More precisely, we describe the relationship between divisors containing $𝒪$ with multiplicity $6$, divisors containing the fourfold $C_2-C_2=\{𝒪(p+q-r-s)\mid p,q,r,s\in C\}$, and divisors singular along $C-C$, using the third exterior...