Jordan Algebras of Symmetric Matrices

Abstract: We study linear spaces of symmetric matrices whose reciprocal is also a linear space. These are Jordan algebras. We classify such algebras in low dimensions, and we study the associated Jordan loci in the Grassmannian.

\section{Introduction}

Let $\mathbb{S}^n$ be the space of symmetric $n \times n$ matrices over the complex numbers $\mathbb{C}$.
The $m$-dimensional subspaces $\mathcal{L} \subset \mathbb{S}^n$
are points in the Grassmannian ${\rm Gr}(m,\mathbb{S}^n)$. The terms
{\em pencil of quadrics} \cite{FMS}  resp.~{\em net of quadrics}~\cite{Wall77} are used when $m=2,3$.
We here consider {\em regular}  $\mathcal{L}$, i.e.
$\mathcal{L}_{\rm inv} = \{ X \in \mathcal{L} \mid {\rm det}(X) \not = 0 \}$ is nonempty.
The {\em reciprocal variety} $\mathcal{L}^{-1}$ is the closure in $\mathbb{S}^n$ of the set
$\{X^{-1} \mid X \in \mathcal{L}_{\rm inv}\}$.
We wish to know under which circumstances
the variety $\mathcal{L}^{-1}$ is a linear space in $\mathbb{S}^n$.
The motivation for this question arises in optimization,
namely in the theory of semidefinite programming \cite{Fay, pablo},
and in statistics, namely from  Gaussian models
that are linear in covariance matrices and  concentration matrices \cite{Jen, STZ}.

\smallskip

The following theorem furnishes a complete answer to our question.\vspace{-0.2cm}

\begin{theorem} \label{thm:characterization}
For $\mathcal{L} \in {\rm Gr}(m,\mathbb{S}^n)$ and $U\in\mathcal{L}_{\rm inv}$,
 the following are equivalent: \vspace{-0.17cm}
 \begin{itemize}
 \item[(a)] The reciprocal variety $\mathcal{L}^{-1}$ is
 also a linear space in $\mathbb{S}^n$. \vspace{-0.17cm}
 \item[(b)] $\mathcal{L}$ is a subalgebra
of the Jordan algebra $(\mathbb{S}^n,\bullet_U)$.\vspace{-0.17cm}
  \item[(c)] $\mathcal{L}^{-1}$ equals $\mathcal{L}$ up to congruence; namely
  $\mathcal{L}^{-1}=U^{-1}\mathcal{L} U^{-1}$.
 \end{itemize}
\end{theorem}\vspace{-0.05cm}

Theorem \ref{thm:characterization} is due to Jensen \cite{Jen}
in the special case when $\mathcal{L}$ contains the
 identity matrix ${\bf 1}_n$. He then uses
$U = {\bf 1}_n$ for the Jordan algebra structure on $\mathbb{S}^n$.\vspace{-0.2cm}

\begin{example}[$m=n=4$]
Consider the two linear subspaces of $\mathbb{S}^4$ given by \vspace{-0.2cm}
$$  \mathcal{L}_1 \,\,=\,\,
\begin{small} \begin{pmatrix}  x & y & 0 & 0 \\
 y & z & 0 & 0 \\
 0 &  0 & w & 0 \\
 0 & 0 & 0 & w \end{pmatrix}\end{small} \quad {\rm and} \quad
\mathcal{L}_2 \,\,= \,\, \begin{small}
\begin{pmatrix}  x & y & 0 & w \\
 y & z & -w & 0 \\
 0 & -w & x & y \\
 w & 0 & y & z \end{pmatrix} \end{small} \quad\,
 {\rm with}\,\,\,\, U = {\bf 1}_4.
  \vspace{-0.2cm}
$$
They satisfy (a),(b),(c) in Theorem \ref{thm:characterization}.
Note that $\mathcal{L}_2$ appears in \cite[eqn (39)]{Jen}. It illustrates
case (ii) in \cite[Theorem 1]{Jen}. The Jordan algebra property fails
when $-w$ is replaced by $w$ in $\mathcal{L}_2$.
The reciprocal variety is a threefold of degree five.
\end{example}\vspace{-0.1cm}

The proof of Theorem~\ref{thm:characterization} is
relatively short. It will be presented in
 Section~\ref{sec2}, along with all relevant definitions.
We also review classification results for Jordan algebras \cite{Jordan_1934}.
This lays the foundation for studying
the {\em Jordan locus} ${\rm Jo}(m,\mathbb{S}^n)$.
This is the subvariety of
the Grassmannian ${\rm Gr}(m,\mathbb{S}^n)$ whose
points are the subspaces $\mathcal{L}$
that satisfy the closed conditions (a), (b) and (c)  above.
To be precise, let $X_1,\ldots,X_m$ represent a basis of $\mathcal{L}$ and write $\adj(U)$ for the adjoint of $U$.
From condition (b) it can be seen that $\Jo(m,\mathbb{S}^n)$ is the subvariety set-theoretically defined (in dual Stiefel coordinates) by the linear dependence~of \vspace{-0.24cm}
\[
X_1,X_2,\ldots,X_m, \,\,X_i\adj(U)X_j+X_j\adj(U)X_i\vspace{-0.1cm}
\]
for all matrices $U\in\mathcal{L}$ and all  indices $1\leq i\leq j\leq m$.

The main goal of this article is the classification of
low-dimensional Jordan subalgebras of $\mathbb{S}^n$
and the description of the corresponding Jordan loci for small~$n$.
In Section \ref{sec3} we determine the loci of Jordan pencils $(m=2)$ and
we study their equations.
Section \ref{sec4} establishes a
lower bound on the codimension of any proper subalgebra of $(\mathbb{S}^n ,\bullet_U)$,
and it identifies the Jordan copencils $(m=4,n=3)$. In Section~\ref{sec5} we classify Jordan subalgebras of dimension $m=3$.
This is applied in Section~\ref{sec6} to study the Jordan loci ${\rm Jo}(3,\mathbb{S}^n)$.
We emphasize the explicit computation of their irreducible components and
defining polynomials. One finding of independent interest is the
Chow matrix in Theorem \ref{thm:chow}
which represents the Chow form of the
determinantal variety  $\{ X \in \mathbb{S}^n \mid  \rk X  \leq n-2 \}$.

\section{Jordan algebras}
\label{sec2}

Fix an invertible matrix $U \in \mathbb{S}^n$. We define an algebra structure on $\mathbb{S}^n$ by setting\vspace{-0.25cm}
\begin{equation} \label{eq:bulletdef} X\bullet_UY \,\,\,:= \,\,\, \frac{1}{2}(XU^{-1}Y\,+\,YU^{-1}X)
\qquad \hbox{for all $X,Y \in \mathbb{S}^n$}.\vspace{-0.25cm}
\end{equation}
When $U$ is clear from context, we drop the subscript.
The product $\bullet$ is bilinear and commutative, with unit $U$.
Associativity fails, but, setting $\,X^{\bullet 2} \,:= \, X \bullet X$,
the following weaker version holds.
This is the axiom for a (unital) {\em Jordan algebra}:\vspace{-0.15cm}
\begin{equation}
\label{eq:jordanaxiom}
X^{\bullet2} \bullet (X \bullet Y) \quad = \quad
X \bullet (X^{\bullet2} \bullet Y) .\vspace{-0.15cm}
\end{equation}
The role of the unit can be played by any matrix $U \in \mathbb{S}^n_{\rm inv}$.
 The resulting Jordan algebra  structure
  is denoted by $(\mathbb{S}^n, \,\bullet_U)$.
  Especially important for applications, notably in  statistics, optimization
  and physics, is the case  when $U$ is real and positive definite.
 Such  algebras  over $\mathbb{R}$ are known as
{\em Euclidean Jordan algebras}.

\begin{remark} \label{rmk:Knotclose}
Over the field $\mathbb{R}$, the isomorphism
type of the Jordan algebra $(\mathbb{S}^n,\,\bullet_U)$ depends on the choice of $U$. To see this,
let  $n=2$ and write $E_{ij}$ for the matrix units.
  If $U = {\bf 1}_2 = E_{11} + E_{22}$ then $(\mathbb{S}^2,\,\bullet_U)$
 has no nilpotent elements.  However, if $\,U' = E_{12} + E_{21}\, $ then the matrices
 $E_{11}$ and $E_{22}$ are nilpotent in   $(\mathbb{S}^2,\,\bullet_{U'})$.
\end{remark}

We return to this issue later: over $\mathbb{C}$, the isomorphism type of $(\mathbb{S}^n,\,\bullet_U)$ is in fact independent of $U$.
First, let us prove the result in the
introduction.

\begin{proof}[Proof of Theorem \ref{thm:characterization}]
Clearly, (c) implies (a).
Suppose that (a) holds, so $\mathcal{L}^{-1}$ is a linear space.
We shall prove that (b) holds.
Consider invertible matrices   $U,X\in\mathcal{L}$. For small $t \in \mathbb{R}$,
the inverse of $X-tU$  is given by the Neumann series\vspace{-0.25cm}
$$
(X-tU)^{-1}\,=\,X^{-1}({\bf 1}_n -tUX^{-1})^{-1}\,=\,
X^{-1}\sum_{k=0}^\infty t^k(UX^{-1})^k \,\,\,\,\in\, \,\, \mathcal{L}^{-1}.\vspace{-0.25cm}
$$
 Since $\mathcal{L}^{-1}$ is linear and closed, the following matrix is in $\mathcal{L}^{-1}$ as well:\vspace{-0.25cm}
 $$
 X^{-1}UX^{-1} \,\, = \,\,\,
\lim_{t\to 0} \, \frac{1}{t}\bigl((X-tU)^{-1}-X^{-1}\bigr) \,\,\,\,\in\, \,\, \mathcal{L}^{-1}.\vspace{-0.25cm}
$$
Its inverse $X\bullet_{U}X$ is in $\mathcal{L}$. Since invertible $X$ are dense in
$\mathcal{L}$, this conclusion holds for all $X\in\mathcal{L}$. Condition (b) states that $X,Y \in \mathcal{L}$
implies $X \bullet_U Y \in \mathcal{L}$. Let $X,Y\in\mathcal{L}$. Then the Jordan algebra
squares $ (X+Y)^{\bullet2}, X^{\bullet2},  Y^{\bullet2}$ are in $\mathcal{L}$,
and we conclude that $ \,X \bullet_U Y  = \frac{1}{2} ((X+Y)^{\bullet2} - X^{\bullet2} -  Y^{\bullet2} )\,$
is also in $\mathcal{L}$. So (b) holds.

Assume (b) and consider any invertible $X\in\mathcal{L} U^{-1}$.
Then $XU\bullet_{U}XU = X^2 U$ is in $\mathcal{L}$.
Next, we see that $X^3 U =X^2U\bullet_{U}XU\in \mathcal{L}$.
By iterating this, every power $X^k$ lies in $\mathcal{L} U^{-1}$.  Since
the characteristic polynomial of $X$ has nonzero constant term,
the Cayley-Hamilton Theorem implies $X^{-1} \in\mathcal{L} U^{-1}$.
Since invertible matrices are dense in $\mathcal{L} U^{-1}$, we get
$\mathcal{L} U^{-1} = (\mathcal{L} U^{-1})^{-1} = U\mathcal{L}^{-1}$, i.e.~(c) holds.
\end{proof}

We next review material from the theory of abstract Jordan algebras.
Already in 1934, Jordan, von Neumann and Wigner~\cite{Jordan_1934}
 established a structure theorem for Euclidean Jordan algebras.
We state a variant that holds over $\mathbb{C}$.
An {\em ideal} $\mathcal{I}$ in a Jordan algebra $\mathcal{A}$ is a subspace such that
$X\bullet Y\in\mathcal{I}$ for all $X \in\mathcal{I}, \,Y \in\mathcal{A}$. The {\em radical} of $\mathcal{A}$ is
the ideal $\rad \mathcal{A}:=\{X\in \mathcal{A}\mid X\bullet Y \textit{ is nilpotent for  all $Y \in \mathcal{A}$}\}$.

\begin{theorem}\label{thm: strucure Jordan Algebras}
For a finite-dimensional complex Jordan algebra $\mathcal{A}$,
 the quotient  $\mathcal{A}/\rad \mathcal{A}$ is a direct sum of
simple Jordan algebras. Every simple Jordan algebra of dimension $\leq 8$ is  isomorphic to either $(\mathbb{C}, \, \cdot \,)$, or
to $(\mathbb{S}^3, \bullet_{{\bf 1}_3})$, or to a {\em spin factor} $\mathbb{C} \times \mathbb{C}^d$,
 $d\geq 2$, with product $(\lambda,a)\bullet (\mu,b)=(\lambda\mu+a^\top b,\lambda b+\mu a)$
 and unit $(1,0)$.
\end{theorem}

\begin{proof}
We find this by applying known structure theorems to semi-simple Jordan algebras of dimension $\leq 8$.  See
 \cite[Ch.~X, \S 3.3]{BraunKoecher1966} or \cite[Part~I, \S 2.13]{McCrimmon2004}. \end{proof}

Theorem~\ref{thm: strucure Jordan Algebras} will be used to classify low-dimensional Jordan algebras.
First, however, we highlight the difference between $\mathbb{R}$ and $\mathbb{C}$
concerning Remark~\ref{rmk:Knotclose}.

\begin{lemm}\label{lemma: squareroot}
Square roots exist in any  complex Jordan subalgebra $\mathcal{L}\subset\mathbb{S}^n$ with unit $U$.
To be precise, for every matrix $X\in\mathcal{L}_{\rm inv}$ there exists $Y\in\mathcal{L}$ with $Y^{\bullet2}=X$.
\end{lemm}

\begin{proof}
The subalgebra $\mathcal{B}=\spa\{X^{\bullet\ell}\mid \ell\geq0\}$ of $\mathcal{L}$ is associative as $Y\bullet Z= YU^{-1}Z$ $=ZU^{-1}Y$ for all $Y,Z\in\mathcal B$. By \cite[Part~I, \S 2.13]{McCrimmon2004}, the only associative simple Jordan algebra in $\mathbb{C}$. So by Theorem~\ref{thm: strucure Jordan Algebras}, the quotient $\mathcal{B}/\rad \mathcal{B}$ is a direct sum of copies of $\mathbb{C}$. In particular, every element of $\mathcal{B}/\rad \mathcal{B}$ is a square. So there exist  $Z \in \mathcal{B}$ and $W\!\in \rad \mathcal{B}$ with $Z^{\bullet2}+W=X$. We have $X^{-1}\in U^{-1}\mathcal{B}U^{-1}$ and so, since $W\in\rad\mathcal{B}$, the matrix $WX^{-1}$ is nilpotent. Hence the matrix $X-W = Z^{\bullet2}$ is invertible and using Theorem~\ref{thm:characterization} we can take $Z^{\bullet-2}:= U\bullet_{Z^{\bullet 2}}U\in\mathcal{B}$ which satisfies $Z^{\bullet2}\bullet Z^{\bullet -2}=U$.
 Since $W\in\rad \mathcal{B}$, the element $Z^{\bullet-2}\bullet W$ is nilpotent, say $(Z^{\bullet-2}\bullet W)^{\bullet p}=0$. Using the Taylor series
 $\sqrt{1+\lambda}= \sum_{\ell=0}^\infty a_\ell \lambda^\ell $, we~see that\vspace{-0.35cm}
$$
\begin{matrix}
X \,\,\,=\,\,\,Z^{\bullet2}\bullet\left(U+Z^{\bullet-2}\bullet W\right)\,\,\,=\,\,\,
Z^{\bullet2}\bullet \left(\,\sum_{\ell=0}^{p-1} a_\ell \left(Z^{\bullet-2}\bullet W\right)^{\bullet\ell}\right)^{\bullet2} . \end{matrix}\vspace{-0.25cm}
$$
Hence $X$ has a square root $Y$ in the subalgebra $\mathcal{B}$ of $\mathcal{L}$.
\end{proof}

Having access to a square root within the algebra, we can find isomorphisms
between $(\mathcal{L},\bullet_U)$ and $(\mathcal{L},\bullet_V)$ for different choices of  matrices $U,V\in\mathcal{L}_{\rm inv}$.

\begin{proposition} \label{prop:compU} Consider a complex Jordan algebra $(\mathcal{L},\bullet_U)$
as in Theorem~\ref{thm:characterization}. There is an isomorphism
    $(\mathcal{L},\bullet_U)\cong(\mathcal{L},\bullet_V)$ of Jordan algebras for all $V \in \mathcal{L}_{\rm inv}$.
\end{proposition}

\begin{proof}
By Lemma \ref{lemma: squareroot}, there exists $W\in\mathcal{L}_{\rm inv}$ with $WU^{-1}W=V$.
Consider the~map\vspace{-0.15cm}
\[
\phi : \mathcal{L} \rightarrow \mathcal{L} \,, \,\,\,
X\,\,\mapsto \,\,VW^{-1}XW^{-1}V\,\,=\,\,2V\bullet_W(X\bullet_WV)-(V\bullet_WV)\bullet_W X. \vspace{-0.15cm}
\]
This linear automorphism of $\mathcal{L}$
is  invertible since $VW^{-1}$ is invertible.  We have\vspace{-0.15cm}
\[
\phi(XU^{-1}X)
\,=\,\phi(X W^{-1} V W^{-1} X)
\,=\,VW^{-1}XW^{-1}VW^{-1}XW^{-1}V
\,=\,\phi(X)V^{-1}\phi(X). \vspace{-0.15cm}
\]
This shows that $\phi$ is a Jordan algebra isomorphism from $(\mathcal{L},\bullet_U)$ to $(\mathcal{L},\bullet_V)$.
\end{proof}

We now come to the classification of low-dimensional Jordan algebras.

\begin{lemm} \label{lem:k+1}
If $k=\dim(\rad\mathcal{A})$, then $X^{\bullet k+1}=0$ for all $X\in\rad\mathcal{A}$.
\end{lemm}

\begin{proof}
Let $X\in\rad\mathcal{A}$. As $X,X^{\bullet2},\ldots,X^{\bullet k+1}\in\rad\mathcal{A}$ are linearly dependent,  the minimum polynomial of $X$ has degree\,$\leq k+1$. As $X$ is nilpotent, $X^{\bullet k+1}=0$.
\end{proof}

\begin{proposition}\label{thm:2d_classification}
Every Jordan algebra $\mathcal A$ of dimension two over $\mathbb{C}$ is isomorphic to
the Jordan algebra $\mathbb{C}\{U,X\}$ with unit $U$, where the  product is given by
$
\begin{array}{llllll}
1\colon\!\!\!\!&X^{\bullet2} =X,&\mbox{\!\!when }\mathcal{A}\cong\mathbb{C}\times\mathbb{C}; \mbox{ or}\\
2\colon\!\!\!\!&X^{\bullet2} =0,&\mbox{\!\!when }\mathcal{A}/\rad\mathcal{A}\cong\mathbb{C}.
\end{array}
$
\end{proposition}

\begin{proof}
By Theorem~\ref{thm: strucure Jordan Algebras}, either
$\mathcal{A}\cong\mathbb{C}\times\mathbb{C}$ or $\mathcal{A}/\rad\mathcal{A}\cong\mathbb{C}$.
We seek $X$ such that $\{U,X\}$ is a basis of $\mathcal{A}$. In the first case, take
the idempotent  $X=(1,0)$  in $\mathbb{C}\times\mathbb{C}$. In the second case, we choose $X$ to span $\rad\mathcal{A}$.
 Then $X^{\bullet 2}=0$ by Lemma \ref{lem:k+1}.
\end{proof}

\begin{theorem}\label{thm:3d_classification}
Every Jordan algebra $\mathcal A$ of dimension three over
$\mathbb{C}$ is isomorphic to the Jordan algebra $\mathbb{C}\{U,X,Y\}$
with unit $U$, where the product is given~by
$
\begin{array}{llllll}
1\!\!\!\!&{\rm (a)}\colon\!\!\!\!&X^{\bullet2} =X,&\!\!\!\! Y^{\bullet2}=Y& \mbox{\!\!\!\!and } X\bullet Y=0,&\mbox{\!\!when }\mathcal{A}\cong\mathbb{C}\times\mathbb{C}\times\mathbb{C}; \\
&{\rm (b)}\colon\!\!\!\!&X^{\bullet2} =U,&\!\!\!\! Y^{\bullet2}=U& \mbox{\!\!\!\!and } X\bullet Y=0,&\mbox{\!\!when }\mathcal{A}\cong\mathbb{C}\times\mathbb{C}^2; \\
2\!\!\!\!&{\rm (a)}\colon\!\!\!\!&X^{\bullet2} =X,&\!\!\!\! Y^{\bullet2}=0& \mbox{\!\!\!\!and } X\bullet Y=0\mbox{ or}&\\
&{\rm (b)}\colon\!\!\!\!&X^{\bullet2} =X,&\!\!\!\! Y^{\bullet2}=0& \mbox{\!\!\!\!and } X\bullet Y=Y/2,&\mbox{\!\!when }\mathcal{A}/\rad\mathcal{A}\cong\mathbb{C}\times\mathbb{C}; \mbox{ or}\\
3\!\!\!\!&{\rm (a)}\colon\!\!\!\!&X^{\bullet2} =Y,&\!\!\!\! Y^{\bullet2}=0& \mbox{\!\!\!\!and } X\bullet Y=0\mbox{ or} \\
&{\rm (b)}\colon\!\!\!\!&X^{\bullet2} =0,&\!\!\!\! Y^{\bullet2}=0& \mbox{\!\!\!\!and } X\bullet Y=0,&\mbox{\!\!when }\mathcal{A}/\rad\mathcal{A}\cong\mathbb{C}.
\end{array}
$
\end{theorem}

\begin{proof}
By Theorem~\ref{thm: strucure Jordan Algebras},  there are four possible cases to consider:
 \vspace{-0.25cm}
\[
\mathcal{A}\,\cong\,\mathbb{C}\times\mathbb{C}\times\mathbb{C},~ \mathcal{A}\,\cong\,\mathbb{C}\times\mathbb{C}^2,
~\mathcal{A}/\rad\mathcal{A}\,\cong\,\mathbb{C}\times\mathbb{C}
\,\, \mbox{ or } \,\,\mathcal{A}/\rad\mathcal{A}\,\cong\,\mathbb{C}.\vspace{-0.25cm}
\]
We seek $X,Y$  such that $\{U,X,Y\}$ is a basis of $\mathcal{A}$.
 In the first case, the choices $X=(1,0,0)$ and $Y=(0,1,0)$ in $\mathbb{C}\times\mathbb{C}\times\mathbb{C}$
 have the desired properties. In the second case, we can choose $X,Y$ to be elements of order two with product zero.
Two such elements in $\mathbb{C}\times\mathbb{C}^2$ are $(0,e_1),(0,e_2)$. So, case 1(a) or case 1(b) holds.

Now suppose $\mathcal{A}/\rad\mathcal{A}\cong\mathbb{C}\times\mathbb{C}$.
Let $X \in \mathcal{A}$ be a preimage of $(1,0)$, so
$X^{\bullet2}\equiv X\!\!\mod\rad\mathcal{A}$. Choose $Y$ to span $\rad\mathcal{A}$. Then $Y^{\bullet2}=0$ and $X\bullet Y\in\rad\mathcal{A}$. Here we use that $\rad\mathcal{A}$ is an ideal. So $X^{\bullet2}=X+\lambda Y$ and $X\bullet Y=\mu Y$ for some $\lambda,\mu\in\mathbb{C}$.

We determine all $\lambda,\mu \in \mathbb{C}$ such that the commutative bilinear form $\bullet$ satisfies \vspace{-0.20cm}
$$
Z^{\bullet2}\bullet(Z\bullet W)=Z\bullet(Z^{\bullet2}\bullet W)\vspace{-0.10cm}
$$
 with $Z=c_1U+c_2X+c_3Y$ and $W=d_1U+d_2X+d_3Y$ for all $c_1,c_2,c_3,d_1,d_2,d_3\in\mathbb{C}$.
 Using computer algebra,
  we find that either $\mu=0$, $\mu=1$ or $(\lambda,\mu)=(0,1/2)$. The first two cases are isomorphic via the base change that replaces $X$ by $U-X$ or $X\pm\lambda Y$. From this we can conclude that case 2(a) holds or case 2(b) holds.

If $\mathcal{A}/\rad\mathcal{A}\cong\mathbb{C}$, we choose $X,Y$ to span $\rad\mathcal{A}$. In this case,  either the square of every element in $\rad\mathcal{A}$ is $0$, or we can choose $X,Y$ such that $X^{\bullet2}=Y$, since $X$ is nilpotent. In the first case, we find that $X^{\bullet2}=Y^{\bullet2}=X\bullet Y=0$. In the second case, we find that $Y^{\bullet2}=X\bullet Y=0$ as $X^{\bullet 3}=0$. So case 3(a) or case 3(b) holds.
\end{proof}

\section{Jordan Pencils}
\label{sec3}

We now embark on the study of the Jordan locus ${\rm Jo}(m,\mathbb{S}^n)$
in ${\rm Gr}(m,\mathbb{S}^n)$. It is convenient to work in
the sub-Grassmannian
${\rm Gr}_{\bf 1}(m,\mathbb{S}^n)$ of subspaces $\mathcal{L}$
that contain the identity matrix ${\bf 1}_n$. The
{\em restricted Jordan locus} is the intersection
\begin{equation}
\label{eq:restrictedJordan} {\rm Jo}_{\bf 1}(m,\mathbb{S}^n) \,\,\, =\,\,\,
{\rm Gr}_{\bf 1}(m,\mathbb{S}^n) \, \cap \, {\rm Jo}(m,\mathbb{S}^n).
\end{equation}

We next examine the variety (\ref{eq:restrictedJordan}) for $m=2$.
Pencils in ${\rm Gr}_{\bf 1}(2,\mathbb{S}^n)$ have the form
$\mathcal{L} = \mathbb{C} \{ {\bf 1}_n, X \}$, where ${\rm trace}(X) = 0$.
This identifies  ${\rm Gr}_{\bf 1}(2,\mathbb{S}^n)$ with the projective space
$ \mathbb{P}^{\binom{n+1}{2}-2}$ of traceless symmetric matrices.
With these conventions,
$\mathcal{L}$ is a Jordan algebra if and only if
the minimal polynomial of $X$ has degree two.

\begin{proposition}
The restricted Jordan locus ${\rm Jo}_{\bf 1}(2,\mathbb{S}^n)$ is a variety with
 $\lfloor n/2 \rfloor$ irreducible components $V_i$ for $1 \leq i \leq n/2$.
The component $V_i$ has codimension $\binom{i+1}{2} + \binom{n-i+1}{2} - 2$.
It  is parametrized by diagonalizable matrices
that have two distinct eigenvalues of multiplicities $i$ and $n-i$.
This is case 1 in Proposition \ref{thm:2d_classification}.
\end{proposition}

\begin{proof}
Following \cite{FMS}, regular pencils are classified by their Segre symbols $\sigma$.
By \cite[Theorem 3.2]{FMS}, a pencil $\mathcal{L}$ has
${\rm deg}(\mathcal{L}^{-1}) = 1$ if and only if
its Segre symbol is $\sigma = [(1,\ldots, 1), (1 ,\ldots, 1)]\,$ or
$\,\sigma = [(2,\ldots, 2 ,1 ,\ldots, 1)]$. The varieties $V_i$ are the
closures of the former strata.
These contain the latter in their closure, by \cite[Theorem 5.1]{FMS}.
The formula for ${\rm codim}(V_i)$ appears in  \cite[Proposition~5.4]{FMS}.
\end{proof}

\begin{remark}
Nondiagonalizable pencils in ${\rm Jo}_{\bf 1}(2,\mathbb{S}^n)$
have only one eigenvalue and Jordan blocks of sizes one or two.
This is case 2 in Proposition \ref{thm:2d_classification}.
\end{remark}

Equations defining  $V_i$ are computed as follows.
Fix a symmetric $n \times n$ matrix of unknowns, $X  = (x_{ij}) $.
Consider the entries of the product $(X-a {\bf 1}_n)(X-b {\bf 1}_n)$,
the $(i+1)$-minors of $X-a {\bf 1}_n$, and the
$(n-i+1)$-minors of $X-b {\bf 1}_n$. By eliminating $a$ and $b$
from these equations, we obtain homogeneous polynomials
in the entries of $X$ that cut out $V_i$ set-theoretically. We
discuss these varieties for $n \leq 5$.
 All pencils of binary quadrics $(n=2)$ are Jordan: ${\rm Jo}_{\bf 1}(2,\mathbb{S}^2) = V_1 ={\rm Gr}_{\bf 1}(2,\mathbb{S}^2)$.

\bigskip

We thus start with Jordan pencils of plane conics ($n=3$).

\begin{example}[$m=2,n=3$]
The variety ${\rm Jo}_{\bf 1}(2,\mathbb{S}^3) = V_1$ is irreducible of codimension $2$
and degree $4$. Its ideal is generated by the following seven cubics: \vspace{-0.17cm}
$$ \begin{small} \begin{matrix}
     x_{11} x_{13} x_{22}-x_{11} x_{13} x_{33}-x_{12}^2 x_{13}+x_{12} x_{22} x_{23}
    -x_{12} x_{23} x_{33}+x_{13}^3-x_{13} x_{22}^2+x_{13} x_{22} x_{33}, \\
     x_{11} x_{12} x_{22}-x_{11} x_{12} x_{33}-x_{12}^3+x_{12} x_{13}^2
    -x_{12} x_{22} x_{33}+x_{12} x_{33}^2+x_{13} x_{22} x_{23}-x_{13} x_{23} x_{33}, \\
     x_{11}^2 x_{23}-x_{11} x_{12} x_{13}-x_{11} x_{22} x_{23}-x_{11} x_{23} x_{33}
    +x_{12} x_{13} x_{22}+x_{13}^2 x_{23}+x_{22} x_{23} x_{33}-x_{23}^3, \\
    x_{12}^2 x_{23}-x_{12} x_{13} x_{22}+x_{12} x_{13} x_{33}-x_{13}^2 x_{23}, \quad
x_{11} x_{13} x_{23}-x_{12} x_{13}^2+x_{12} x_{23}^2-x_{13} x_{22} x_{23}, \\ \!\!\!\!
x_{11} x_{12} x_{23}-x_{12}^2 x_{13}-x_{12} x_{23} x_{33}+x_{13} x_{23}^2\,, \quad
 x_{11}^2 x_{22}-x_{11}^2 x_{33}-x_{11} x_{12}^2  +x_{11} x_{13}^2 \quad \\ \qquad -x_{11} x_{22}^2
+x_{11} x_{33}^2+x_{12}^2 x_{22}-x_{13}^2 x_{33}+x_{22}^2 x_{33}
-x_{22} x_{23}^2-x_{22} x_{33}^2+x_{23}^2 x_{33}.
\end{matrix}
\end{small} \vspace{-0.13cm}
$$
These are the expressions in the sum of squares representation of the
discriminant of the characteristic polynomial of $X$, seen in \cite[page 97]{CBMS}.
Indeed, $V_1$ is also the Zariski closure of all \underbar{real}
 matrices in $\mathbb{S}^3$ with a double eigenvalue.
 \end{example}

\begin{example}[$m=2,n=4$]
We compute  ${\rm Jo}_{\bf 1}(2,\mathbb{S}^4) $ in {\tt Macaulay2} as follows: \vspace{-0.13cm}
\begin{small}
\begin{verbatim}
  R = QQ[a,b,x11,x12,x13,x14,x22,x23,x24,x33,x34,x44];
  X = matrix {{x11,x12,x13,x14},{x12,x22,x23,x24},
              {x13,x23,x33,x34},{x14,x24,x34,x44}};
  I = eliminate({a,b},minors(1,(X-a)*(X-b)))
\end{verbatim}
\end{small} \vspace{-0.13cm}
The ideal ${\tt I}$ is generated by $30$ cubics. It is the intersection
of two prime ideals, corresponding to
${\rm Jo}_{\bf 1}(2,\mathbb{S}^4) = V_1 \cup V_2$.
The component $V_1$ has codimension $5$ and degree $8$. Its
prime ideal is generated by $10$ quadrics, including
$\,x_{13} x_{24}-x_{12} x_{34},\, x_{14} x_{23}-x_{12} x_{34}$ and
$x_{13} x_{23}-x_{14} x_{24}-x_{12} x_{33}+x_{12} x_{44}$.
The component $V_2$ has codimension $4$ and degree  $6$. Its prime ideal
is generated by $9$ quadrics.
\end{example}

\begin{example}[$m=2,n=5$] The ideal of
${\rm Jo}_{\bf 1}(2,\mathbb{S}^5) = V_1 \cup V_2$ is generated by $81$ cubics.
The variety $V_1$ has codimension $9$ and degree $16$.
Its ideal is generated by $35$ quadrics, including ten
$2 \times 2$-minors,
like $x_{12} x_{34}-x_{13} x_{24}$.
 The variety $V_2$ has codimension $7$ and degree $40$.
Its  ideal is generated by $95$ cubics, e.g.~$ \begin{small}
x_{12} x_{13} x_{35}
- x_{13} x_{15} x_{23}
+ x_{22} x_{23} x_{35}
- x_{23}^2 x_{25}
- x_{23} x_{34} x_{45}
- x_{23} x_{35} x_{55}
+ x_{24} x_{34} x_{35}
+ x_{25} x_{35}^2.
\end{small}
$
\end{example}

We now observe that ${\rm Jo}(m,\mathbb{S}^n)$ is the orbit of (\ref{eq:restrictedJordan}) under the
congruence action by ${\rm GL}(n)$. This implies the following result in the
Grassmannian ${\rm Gr}(2,\mathbb{S}^n)$.

\begin{corollary} \label{cor:n2floor}
The Jordan locus ${\rm Jo}(2,\mathbb{S}^n)$ has $\lfloor n/2\rfloor$ irreducible components.
Using the notation from  \cite[Section 5]{FMS}, these components are
the Grassmann strata $\overline{{\rm Gr}_\sigma}$  that are
associated with the Segre symbols
$\sigma = [(1,\ldots, 1), (1 ,\ldots, 1)]$.
\end{corollary}

\begin{remark}
Jordan algebras can be viewed as a nonabelian generalization of {\em partition matroids}.
These are the matroids that are direct sums of uniform matroids. Indeed,
the diagonalizable types in Corollary \ref{cor:n2floor} correspond to
the rank-$2$ partition matroids, and those in
Theorem \ref{thm:classification m=3,n=4}
to the rank-$3$ partition matroids.
\end{remark}

\section{Jordan Copencils}
\label{sec4}

This section concerns
Jordan algebras of low codimension.
These are quite rare:

\begin{theorem} \label{Thm: codimension}
Let $\mathcal{L}\subset\mathbb{S}^n$ be a proper Jordan subalgebra. Then $\codim(\mathcal{L})\geq n-1$.
\end{theorem}

To prove this theorem, we need the following result
and some  lemmas.\vspace{-0.08cm}

\begin{theorem}[Peirce decomposition, {\cite[Part II, Chapter 8]{McCrimmon2004}}]
\label{thm:peirce}
Let $X_1,\dots X_d$ be orthogonal idempotents in $\mathcal{L}$ with $U=X_1+\dots+X_d$. Then $\mathcal{L}=\bigoplus_{1\leq i\leq j\leq d} \mathcal{X}_{i,j}$ where
 $\,\mathcal X_{i,i}=\{\,Y\in \mathcal{L}\,\mid \,Y= X_i\bullet Y\}\,$ and
$\,\mathcal X_{i,j}=\{\,Y\in \mathcal{L}\,\mid \, Y = 2 X_i\bullet Y= 2 X_j\bullet Y\}\,$ for $\,i<j$.
\end{theorem}

Here, being  {\em orthogonal} means that $X_i \bullet X_j =0$ whenever $i \not= j$. An idempotent $X$ is called
{\em primitive} if there are no nonzero orthogonal idempotents $X_1$ and $ X_2$ with $X=X_1+X_2$.
We denote by $J_n$ the $n\times n$ matrix with ones on its anti-diagonal.

\begin{lemm}[{\cite[Lemma 1]{bukovsek-omladic}}]\label{lemma:conjugation0}
If  $X,Y\in\mathbb{S}^n$ are similar, then they are orthogonally congruent.
\end{lemm}

\begin{lemm}\label{lemm:antitriagle}
If the unit $U$ of a Jordan subalgebra $\mathcal{L} \subset \mathbb{S}^n$ is primitive, then $\mathcal{L}$ is congruent to $\,\mathbb{C} J_n\oplus \mathcal{L}'$ where $\mathcal{L}'$ consists of upper anti-triangular matrices.
\end{lemm}

\begin{proof} After congruence, we may assume $U = {\bf 1}_n$.
Let $X\in\mathcal{L}$.  By considering the Jordan canonical form of $X$, adding multiples of ${\bf 1}_n$ to $X$, and taking powers, we see that ${\bf 1}_n$ cannot be primitive if $X$ has distinct eigenvalues. So $\mathcal{L}=\mathbb{C} {\bf 1}_n\oplus \mathcal{L}'$ where $\mathcal{L}'\subset\mathcal{L}$ is the subset of nilpotent matrices. By Lemma~\ref{lemma:conjugation0}, we see that $XY+YX\in\mathcal{L}'$ for  $X,Y\in\mathcal{L}'$. So $\mathcal{L}'$ is anti-triangularizable by \cite[Theorem 9]{bukovsek-omladic}. The transformation given there takes the unit matrix ${\bf 1}_n$ to the matrix $J_n$.
\end{proof}

Two spaces $\mathcal{L},\mathcal{L}'\subset\mathbb{S}^n$ are congruent if $\mathcal{L}'=P\mathcal{L} P^\top$ for some $P\in\GL(n)$.
We say that $\mathcal{L},\mathcal{L}'$ are {\em orthogonally congruent} if $P$ lies in
the orthogonal group $O(n)$.\vspace{-0.08cm}

\begin{lemm}\label{lemma:conjugation}
If  $\mathcal{L},\mathcal{L}'\subset\mathbb{S}^n$ are similar, then they are orthogonally congruent.
\end{lemm}

\begin{proof}
We have $\mathcal{L} P=P\mathcal{L}'$ for some  $P\in\GL(n)$.
Consider  $A\in\mathcal{L}$ and $B\in\mathcal{L}'$ with $AP=PB$. By the proof of \cite[Lemma 1]{bukovsek-omladic}, we get $AU=UB$ for $U=(PP^\top)^{-1/2}P$.
The matrix $U$ is orthogonal, and it maps $\mathcal{L}$ into $\mathcal{L'}$
under congruence.
\end{proof}

\begin{lemm}\label{lemma:idempotent}
Let $X\in\mathbb{S}^n$ be an idempotent matrix of rank $r$. Then $X$ is orthogonally congruent to $\Diag({\bf 1}_r,{\bf 0}_{n-r})$. If $X=\Diag({\bf 1}_r,{\bf 0}_{n-r})$ and $Y\in\mathbb{S}^n$ is a matrix with $XY+YX=0$, then $Y=\Diag({\bf 0}_r,Z)$ for some $Z\in\mathbb{S}^{n-r}$.
\end{lemm}

\begin{proof} This follows from   Lemma~\ref{lemma:conjugation0}.
\end{proof}

\begin{proof}[Proof of Theorem \ref{Thm: codimension}]
Either $\mathcal{L}$ contains a primitive idempotent of rank $>1$, or
the unit $U$ of $\mathcal{L}$ is a sum of orthogonal rank-$1$ idempotents. We shall examine the associated Peirce decompositions
(Theorem \ref{thm:peirce}). These have $d=2$ and $d=n$.

Suppose that $\mathcal{L}$ contains a primitive idempotent $X_1$ of rank $r>1$ and take $X_2=U-X_1$. Then $\mathcal{L}=\mathcal{X}_{1,1}\oplus\mathcal{X}_{1,2}\oplus\mathcal{X}_{2,2}$. By Lemmas \ref{lemm:antitriagle} and \ref{lemma:idempotent}, we assume that $X_1=\Diag(J_r,{\bf 0}_{n-r})$, $X_2=\Diag({\bf 0}_r,{\bf 1}_{n-r})$ and $\mathcal{X}_{1,1}$ consists of upper anti-triangular $r \times r$ matrices. So $\mathcal{X}_{1,1}$ has codimension $\geq r-1$ in $\mathbb{S}^r$. Elements in $\mathcal X_{1,2}$ are of the form $Z={\footnotesize
 \begin{pmatrix}
 {\bf 0}_r & V^\top\\
 V &  {\bf 0}_{n-r}
 \end{pmatrix}
}$ with $Z^{\bullet 2}=\Diag(V^\top V,VJ_rV^\top)$. It follows that $V^\top V\in\mathcal{X}_{1,1}$.
So, the last column $v$ of $V$ satisfies $v^\top v=0$. Hence $v$ is contained in a subspace $\mathcal Y\subset\mathbb{C}^{n-r}$ of codimension $c\geq (n-r)/2$. We have $\mathcal{X}_{2,2}=\Diag({\bf 0}_r,\mathcal{Z})$ where $\mathcal{Z}\subset\mathbb{S}^{n-r}$ consists of matrices that map $\mathcal Y$ into itself. Therefore $\mathcal{Z}$ has codimension $\geq c(n-r-c)$. So $\mathcal{L}$ has codimension $\geq r-1+c+c(n-r-c)\geq n-1$.

If all primitive idempotents of $\mathcal{L}$ have rank $1$, then $U=X_1+\dots+X_n$ for rank-$1$ orthogonal idempotents $X_1,\ldots,X_n$. By Lemma~\ref{lemma:idempotent}, we may assume $X_i=E_{ii}$. We
note that the spaces of the Peirce decomposition are either $\mathcal{X}_{i,j} =\mathbb{C}\{E_{i,j}+E_{j,i}\}$ or $\mathcal{X}_{i,j}=\{0\}$.  Since $\mathcal{L}$ is a proper subalgebra, at least one of these spaces is $\{0\}$, say $\mathcal X_{1,2}=\{0\}$. Since $2(E_{i,j}+E_{j,i})\bullet (E_{i,k}+E_{k,i})=E_{k,j}+E_{j,k}$, this implies either $\mathcal X_{1,i}=\{0\}$ or $\mathcal X_{2,i}=\{0\}$ for each $i \in \{3,\dots, n\}$. Therefore $\codim (\mathcal{L})\geq n-1$.
\end{proof}

A linear subspace $\mathcal{L}\subset\mathbb{S}^n$ is said to be a
 {\em copencil} if $\codim (\mathcal{L}) = 2$, i.e.~the orthogonal complement $ \mathcal{L}^\perp =  \{\,
X \in \mathbb{S}^n \mid {\rm trace}(XZ) = 0 \,\,\hbox{for} \,\, Z \in \mathcal{L} \} $ is a pencil.
We say that $\mathcal{L}$ is a {\em Jordan copencil} if $\mathcal{L}$ is a copencil that satisfies the equivalent
conditions (a), (b) and (c) from Theorem \ref{thm:characterization}.
Theorem \ref{Thm: codimension} implies:

\begin{corollary} There are no Jordan copencils unless $n \leq 3$.
Every copencil with $n=2$ is a Jordan copencil since it is the span of
one invertible matrix $U $ in $\mathbb{S}^2$.
\end{corollary}

It remains to study $n=3$. We represent  Jordan  copencils  by the variety
\begin{equation}
\label{eq:copencils3}
 \left\{\, (X,Y) \in \mathbb{S}^3 \times \mathbb{S}^3 \,\,\middle|\,\,
(\mathbb{C} \{X,Y\} )^\perp \,\,\hbox{is a point in} \,\, \,{\rm Jo}(4,\mathbb{S}^3) \, \right\} .
\end{equation}

\vspace{5pt}

\begin{proposition}\label{thm:classification m=4,n=3}
There are two congruence orbits of Jordan copencils in $\mathbb{S}^3$:\vspace{-0.15cm}
$$ \mathcal{L}_1 \,\,=\,\, \begin{small} \begin{pmatrix}
\, x & 0 & 0 \, \\ \,0 &y&w \, \\ \, 0 &w&z\,
\end{pmatrix} \end{small} \quad \mbox{ and } \quad
\mathcal{L}_2\, \,= \,\,
\begin{small} \begin{pmatrix}
\, x&y&w \, \\ \, y&z & 0 \,\\ \, w & 0 & 0 \,
\end{pmatrix}. \end{small}\vspace{-0.15cm}
$$
The orbit of $\mathcal{L}_2$ is in the closure of
the orbit of  $\mathcal{L}_1$. These orbits have codimensions $5$ and $4$, so the Jordan
locus ${\rm Jo}(4, \mathbb{S}^3)$ is irreducible of codimension~$4$.
The prime ideal of (\ref{eq:copencils3}) has degree $21$ and is
 generated by $4$ cubics and $15$ quartics.
\end{proposition}

\begin{proof}
Theorem~\ref{thm: strucure Jordan Algebras} gives a finite list of cases to examine.
By multiple applications of Lemma~\ref{lemma:idempotent}, we see that $\mathbb{C}\times\mathbb{C}\times\mathbb{C}\times\mathbb{C}$ cannot be embedded into~$\mathbb{S}^3$. By Theorem~\ref{thm:classification m=3,n=n}, $\mathbb{C}\times\mathbb{C}^2$ cannot be embedded into $\mathbb{S}^3$.
And, a computation shows that the Jordan algebra $\mathbb{C}\{X,Y\}$ with zero bilinear form $\bullet$ cannot be embedded into $\mathbb{S}^3$. Using the help of a computer as in the proof of Theorem~\ref{thm:3d_classification}, we generate a finite list of
Jordan algebras that do not contain any of these Jordan algebras as a subalgebra. Using Lemma~\ref{lemma:idempotent}, we find that all embeddings of these Jordan algebras are
congruent  to $\mathcal{L}_1$ or $\mathcal{L}_2$.

The last two statements are proved by computer algebra.
Note that
$\mathcal{L}_1^\perp $ and $\mathcal{L}_2^\perp$ are singular pencils.
 The four cubics in the ideal for
$\mathcal{L}_1$ are the four coefficients of  ${\rm det}(x X + y Y)$. The
$15$ quartics can also be derived using linear algebra. For the containment of orbits, we
check that $\mathcal{L}_2$ satisfies these $19$ equations.
\end{proof}

\section{Three-Dimensional Jordan Algebras}
\label{sec5}

The conditions in Theorem~\ref{thm:characterization} are invariant under congruence.
Ideally, one would classify Jordan subalgebras of $\mathbb{S}^n$ up to congruence and
describe the poset of orbit closures.
We tackle this problem for $m=3$.
The following two theorems are our main results in Section \ref{sec5}.
The first gives a complete solution for $4 \times 4$ matrices.
The second gives a solution for Jordan subalgebras without a radical.

\begin{theorem}\label{thm:classification m=3,n=4}
Every Jordan net in $\mathbb{S}^4$ is congruent to one of the following nets:

\noindent$
\begin{array}{llllll}
1\!\!\!\!&{\rm (a)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&\Diag(x{\bf 1}_2,y,z);\\
&{\rm (b)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&\mbox{the Kronecker product of  ${\bf 1}_2$
with $\mathbb{S}^2$;
\qquad (see $\mathcal{L}_2$ in Example   \ref{ex:drei})}
\\
2\!\!\!\!&{\rm (a1)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&\Diag(xJ_2+yE_{11},z{\bf 1}_2);\\
&{\rm (a2)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&\Diag(xJ_3+yE_{11},z);\\
&{\rm (b)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&\Diag(xJ_2,yJ_2)+z(E_{13}+E_{31});\\
3\!\!\!\!&{\rm (a)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&x\Diag(J_3,1)+y(E_{12}+E_{21})+zE_{11};\\
&{\rm (b1)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&xJ_4+yE_{11}+zE_{22};\mbox{ or}\\
&{\rm (b2)}\!\!\!\!\!\!\!\!\!\!&\colon\!\!\!\!&xJ_4+yE_{11}+z(E_{13}+E_{31}).
\end{array}
$

\smallskip

\noindent The numberings of the congruence classes are the same as in
   Theorem \ref{thm:3d_classification}.
Their orbit closures are subvarieties in the Grassmannian ${\rm Gr}(3,\mathbb{S}^4)$ of dimension~$21$.

\begin{samepage}

The codimensions of these Jordan loci
and their containments are as follows:

\vspace*{-0.07cm}

\begin{center}
\begin{tikzpicture}

	\node () at (-3,2.4) {$\codim$};
\node () at (-3,1.6) {$\codim$};
\node () at (-3,.8) {$\codim$};
\node () at (-3,0) {$\codim$};
	\node () at (-2,2.4) {$9$};
	\node () at (-2,1.6) {$10$};
	\node () at (-2,.8) {$11$};
	\node () at (-2,0) {$12$};

    \node (1a) at (1,2.4) {1(a)};
    \node (1b) at (4,1.6) {1(b)};
    \node (2a1) at (0,1.6) {2(a1)};
    \node (2a2) at (2,1.6) {2(a2)};
    \node (2b) at (5.,.8) {2(b)};
    \node (3a) at (1.,.8) {3(a)};
    \node (3b1) at (3.,.8) {3(b1)};
    \node (3b2) at (3.,0) {3(b2)};

    \draw [thick] (1a) -- (2a1);
    \draw [thick] (1a) -- (2a2);
    \draw [thick] (1b) -- (2b);
    \draw [thick] (1b) -- (3b1);
    \draw [thick] (2a1) -- (3a);
    \draw [thick] (2a2) -- (3b1);
    \draw [thick] (2b) -- (3b2);
    \draw [thick] (2a2) -- (3a);
    \draw [thick] (3a) -- (3b2);
    \draw [thick] (3b1) -- (3b2);
\end{tikzpicture}
\end{center}
\end{samepage}
\end{theorem}

\begin{corollary} \label{cor:Jo34}
The Jordan locus ${\rm Jo}(3,\mathbb{S}^4)$ has two irreducible components.
\end{corollary}

Jordan algebras of types 1(a) and 1(b) are of primary interest also for $n \geq 5$.

\begin{theorem}\label{thm:classification m=3,n=n}
Every Jordan net in $\mathbb{S}^n$ of type 1(a) is congruent to a diagonal~net
$ \Diag(x{\bf 1}_{k_1},y{\bf 1}_{k_2},z{\bf 1}_{k_3})$
for $k_1\geq k_2\geq k_3\geq1$ with sum $n$.
Jordan nets of type 1(b)~exist only for even $n$.
They are congruent to the Kronecker product of ${\bf 1}_{n/2}$ with~$\mathbb{S}^2$.
\end{theorem}

We note that
Theorem \ref{thm:classification m=3,n=n}
does not account for all irreducible components of
${\rm Jo}(3,\mathbb{S}^n)$ when $n \geq 5$. At present we do not
know these decompositions.

\begin{proposition} \label{prop:tau}
The Jordan net $\mathcal{L}_* =uJ_5+x(E_{15}+E_{51})+y(E_{14}+E_{41}) $ in $\mathbb{S}^5$ has type 2(b).
It  is not in the closure of the diagonalizable locus 1(a) in ${\rm Jo}(3,\mathbb{S}^5)$.
\end{proposition}

\begin{proof}
For any linear space $\mathcal{L} \subset \mathbb{S}^n$ we consider
 $\tau(\mathcal{L}):=\min\{\rk X\mid X\in\mathcal{L} \setminus \{0\}\}$.
 This invariant equals $1$ for all Jordan algebras $\mathcal{L}$ of type 1(a), since
 $k_3 = 1$ for $n=5$.
 For the special net above, we see that $\tau(\mathcal{L}_*) = 2$.
 It now suffices to note that $\{\mathcal{L}\in\Gr(m,\mathbb{S}^n)\mid \tau(\mathcal{L})\leq k\}$  is a
 closed set because it is the projection along $\mathbb{P}(\mathbb{S}^n)$ of the incidence variety
$\{(\mathcal{L},X)\in\Gr(m,\mathbb{S}^n)\times\mathbb{P}(\mathbb{S}^n)\mid X\in\mathcal{L},\rk X\leq k\}$.
\end{proof}

The situation is easier if we restrict to the case of interest in applications.
Let $\mathcal{L}$ be a real Jordan subalgebra of symmetric matrices whose unit $U$ is positive definite. Then $\mathcal{L}$ is {\em formally real}, i.e. if $X_1^{\bullet2}+\cdots+X_k^{\bullet 2}=0$, then $X_1,\ldots,X_k=0$.  Formally real Jordan algebras have zero radical \cite{Jordan_1934} and this persists if we view
them as Jordan algebras over $\mathbb{C}$.
We define the {\em Euclidean Jordan locus} ${\rm EJo}(3,\mathbb{S}^n)$ to be the closure
in ${\Gr}(3,\mathbb{S}^n) $ of the set of such Jordan algebras.
 For our study with $m=3$, we get that formally real Jordan nets
  have type 1(a) or~1(b).

\begin{corollary}
The number of irreducible components of
${\rm EJo}(3,\mathbb{S}^n)$ is the coefficient of the $t^n$ in the generating function
\[
t^3/((1-t)(1-t^2)(1-t^3)) + t^2/(1-t^2).
\]
\end{corollary}
\begin{proof}
The first summand is the generating function of triples of $a,b,c\geq0$ with $a+2b+3c=n-3$. These correspond to triples of $k_1\geq k_2\geq k_3\geq 1$ with sum $n$ via  \vspace{-0.16cm}
\[
(a,b,c)\mapsto(a+b+c+1,b+c+1,c+1). \vspace{-0.16cm}
\]
The second generating function accounts for the Jordan nets of type 1(b).
\end{proof}

In the remainder of this section we present the proofs of our two theorems.

\begin{proof}[Proof of Theorem \ref{thm:classification m=3,n=n}]
Let $\mathcal{L}$ be a Jordan net of type 1(a). Then $\mathcal{L}$ is congruent to
$\,\Diag(x{\bf 1}_{k_1},y{\bf 1}_{k_2},z{\bf 1}_{k_3}) $.
This is seen by two applications of Lemma~\ref{lemma:idempotent}.
Suppose $\mathcal{L}$ has unit $U = {\bf 1}_n$ and type 1(b). After applying an orthogonal congruence, we may assume that $X=\Diag({\bf 1}_r,-{\bf 1}_{n-r})$ for some $1\leq r\leq n-1$. The condition $XY+YX=0$ now implies
$ Y= \begin{small} \begin{pmatrix} 0&P\\
P^\top&0\end{pmatrix} \end{small} $
for some $P\in\mathbb{C}^{r\times(n-r)}$. The condition $Y^2={\bf 1}_n$ implies  $PP^\top={\bf 1}_r$ and $P^\top P={\bf 1}_{n-r}$. This is only possible when $2r=n$. By applying the orthogonal
 congruence $\Diag(P^\top, {\bf 1}_r)$, we get the desired form.
\end{proof}

\begin{proof}[Proof of Theorem \ref{thm:classification m=3,n=4}]
We proceed case-by-case, starting from the classification of
Jordan algebras in Theorem~\ref{thm:3d_classification}.
Let $U,X,Y$ be a basis of $\mathcal{L}$, where we
 assume $U={\bf 1}_4$. The types 1(a) and 1(b) are already covered by the previous proof.

Suppose that $\mathcal{L}$ has type 2(a). We can assume $X=\Diag({\bf 1}_{4-k},{\bf 0}_k)$ and $Y=\Diag({\bf 0}_{4-k},Z)$ with $Z^2={\bf 0}_k$ for some $k\in\{2,3\}$.
Up to orthogonal congruence, there is only one
$Z \in \mathbb{S}^k$ with $Z^2 = {\bf 0}_k\neq Z$,
as there is only one option for its Jordan canonical form.
This uses Lemma~\ref{lemma:conjugation0}.
 So, we get two classes of this type.

Suppose that $\mathcal{L}$ has type 2(b). Then we can assume that $X=\Diag({\bf 1}_k,{\bf 0}_{4-k})$. The condition $XY+YX=Y$ implies
$ Y= \begin{small} \begin{pmatrix} 0&P\\
P^\top&0\end{pmatrix} \end{small} $
for some matrix $P\in\mathbb{C}^{k\times(4-k)}$. The condition $Y^2={\bf 0}_4$ implies $PP^\top={\bf 0}_k$ and $P^\top P={\bf 0}_{4-k}$. As $P$ is nonzero, this is only possible when $k=2$. So $P=ab^\top$ for some nonzero vectors $a,b\in\mathbb{C}^2$ with $a^\top a=b^\top b=0$. Such matrices $Y$ form a single orbit of $\Diag({\rm O}(2),{\rm O}(2))$. So we get one class of this type.
Suppose that  $\mathcal{L}$ has type 3(a). Then $X^2=Y$ and $X^3=0$. By Lemma~\ref{lemma:conjugation0}
and appealing to Jordan canonical forms as above, there is a unique $X$ up to
    orthogonal congruence. So we get one class of this type.

Consider $\mathcal{L}$ of type 3(b). We have $X^2=Y^2=XY+YX=0$. Since $X$, $Y$ or~$X+Y$ has rank~$2$, we may assume $\rk X=2$. One can show, using Lemma~\ref{lemma:conjugation0},
 that we may assume $\rk Y=1$.  Using Lemma~\ref{lemma:conjugation}
 and  conjugating $\mathcal{L}$ in $\mathbb{C}^{n \times n}$, we~get\vspace{-0.15cm}
$$
X=\begin{pmatrix}
0&{\bf 1}_2\\
0&0
\end{pmatrix}\,\,\,\mbox{ and } \,\,\, Y=\begin{pmatrix}
0&Z\\0&0
\end{pmatrix}\vspace{-0.15cm}
\,\,\quad \hbox{where $\,\rk Z =1$.}
$$
 By conjugating with a matrix $\Diag(P,P)$ and scaling, we find $Z=E_{22}$ or $Z=E_{12}$.
 So we get two classes. In one of them,  $\rad\mathcal{L}$ is spanned by two rank-$1$ matrices.

To identify the poset of inclusions,
 consider the associated nets of quadrics
 $ (a,b,c,d)\mathcal{L}(a,b,c,d)^T$.
In this notation, type 1(a) is written as $x(a^2+b^2)+yc^2+zd^2$.
Replacing $(a,b,c,d)$ by $(a,b,c,c+td)$ in 1(a),
where $t \in \mathbb{R} \setminus \{0\}$, we get \vspace{-0.15cm}
$$
x(a^2+b^2)+y'c^2+z'(2cd+td^2)\vspace{-0.15cm}
$$
for $y'=y+z$ and $z'=tz$. As $t\to0$, we obtain type 2(a1).  Similarly:

$
\begin{array}{llclcl}
\bullet\mbox{ We get ~1(a)}&\!\!\!\to\!\!\!&\mbox{2(a2)}&\mbox{using }(c+t^2a,tb,c,d).\\
\bullet\mbox{ We get ~1(b)}&\!\!\!\to\!\!\!&\mbox{2(b)}&\mbox{using }(-a,b,i(a+tc),i(b+td)).
\end{array}
$

$
\begin{array}{llclcl}
\bullet\mbox{ We get ~1(b)}&\!\!\!\to\!\!\!&\mbox{3(b1)}&\mbox{using }(a,tb,tc,d).\\
\bullet\mbox{ We get ~2(a1)}&\!\!\!\to\!\!\!&\mbox{3(a)}&\mbox{using }(a,b,a+tb+t^2c,td).\\
\bullet\mbox{ We get ~2(a2)}&\!\!\!\to\!\!\!&\mbox{3(a)}&\mbox{using }(a,td,t(tc-b),a+tb).\\
\bullet\mbox{ We get ~2(a2)}&\!\!\!\to\!\!\!&\mbox{3(b1)}&\mbox{using }(a,d+tb,c,d).\\
\bullet\mbox{ We get ~2(b)}&\!\!\!\to\!\!\!&\mbox{3(b2)}&\mbox{using }(c+ta,d+tb,c,d).\\
\bullet\mbox{ We get ~3(a)}&\!\!\!\to\!\!\!&\mbox{3(b2)}&\mbox{using }(a,b,tc,i(b-td)).\\
\bullet\mbox{ We get ~3(b1)}&\!\!\!\to\!\!\!&\mbox{3(b2)}&\mbox{using }(a,a+tb,c,-c+td).
\end{array}
$

\noindent The codimensions were found by direct computation.
It remains to show that there is no inclusion of orbit closures when
there is no edge in our Hasse diagram.
The arrow 1(a) $\rightarrow$ 2(b) is missing because $\tau$ cannot go up in limits;
see Proposition \ref{prop:tau}.
  The absence of other inclusions is certified by explicit
  polynomials that vanish on one orbit, but not on another. For example, the polynomial for
  the orbit 1(b)    displayed in (\ref{eq:onerep})
  does not vanish on the  orbits 1(a), 2(a1), 2(a2) and 3(a).
  Similarly, $4 p_{012} p_{146} - 4 p_{013} p_{145} -
   p_{014} p_{056} - p_{014} p_{126} + p_{014} p_{135} - p_{014} p_{234}
       + p_{024} p_{046} - p_{034}  p_{045}    \,   $ vanishes on
        2(a1) but not on  3(b1).
  \end{proof}

\section{Nets of Quadrics}
\label{sec6}

In this section we offer a more detailed study of Jordan loci in the case $m=3$.

\begin{proposition} \label{prop:33}
The Jordan locus ${\rm Jo}(3,\mathbb{S}^3) $ consists
of nets that are generated by three rank-$1$ conics in $\mathbb{P}^2$.
This  is an irreducible variety of codimension $3$
in the Grassmannian ${\rm Gr}(3,\mathbb{S}^3) $, and it has
degree $57$ in the Pl\"ucker embedding in $ \mathbb{P}^{19}$.
Its prime ideal is generated by $62$ quadrics
in the $20$ Pl\"ucker coordinates.
\end{proposition}

\begin{proof}
Here we use the Wall's classification \cite{Wall77} of nets of conics.
For each net $\mathcal{L}$ in his list, we computed
$\mathcal{L}^{-1}$. The surface $\mathcal{L}^{-1}$ is  linear
only for Wall's types E, G and H. According to
\cite[Figure 5]{Wall77}, the strata G and H are contained in the closure of
stratum E, which has codimension $3$, by \cite[Table 2]{Wall77}.
 We see in \cite[Table 1]{Wall77} that E corresponds to nets
spanned by rank-$1$ conics. Hence the closure of stratum E is
an irreducible variety. It equals ${\rm Jo}(3,\mathbb{S}^3) $.
The results about the ideal and degree of this Jordan locus are  in
\cite[Proposition 3.1]{BS}.
\end{proof}

The ideal in Proposition \ref{prop:33} simplifies greatly if
we restrict to  $ {\rm Jo}_{\bf 1}(3,\mathbb{S}^3)$. Namely, for the net
$\mathcal{L} = {\rm span} \{{\bf 1}_3,\, X,\, Y \}$, with
$X = (x_{ij})$ and $Y= (y_{ij})$, we obtain \vspace{-0.14cm} \begin{equation}
\label{eq:buchloe} \begin{small}
\begin{matrix}
\!\! \bigl\langle \,
 (x_{11} y_{12}-x_{12} y_{11})+(x_{12} y_{22}-x_{22} y_{12})+(x_{13} y_{23}-x_{23} y_{13}), \\ \quad
 (x_{11} y_{23} - x_{23}  y_{11} )+(x_{12} y_{13}-x_{13} y_{12})+2 (x_{22} y_{23}-x_{23} y_{22}), \\ \qquad
  2 (x_{11} y_{13}-x_{13} y_{11})+(x_{12} y_{23}-x_{23} y_{12})-(x_{13} y_{22}-x_{22} y_{13}) \bigr\rangle .
\end{matrix} \end{small} \vspace{-0.14cm}
\end{equation}
This prime ideal of codimension $3$ defines the Jordan locus
in  ${\rm Gr}_{\bf 1}(3, \mathbb{S}^3)$.

We now assume $n \geq 3$.
For any net $\mathcal{L} \in {\rm Gr}(3,\mathbb{S}^n)$,
the reciprocal  $\mathcal{L}^{-1}$
 is a surface in $\mathbb{P} (\mathbb{S}^n)
= \mathbb{P}^{\binom{n+1}{2}-1}$.
This is a linear plane precisely when $\mathcal{L}$ is a Jordan algebra.
We now examine the other extreme case, when
$\mathcal{L}$ is generic in~${\rm Gr}(3,\mathbb{S}^n)$.

\begin{theorem} \label{thm:chow}
 The following conditions are equivalent for a net of quadrics $\mathcal{L}$:  \vspace{-0.15cm}
\begin{enumerate}
\item Every nonzero complex matrix in $\mathcal{L}$ has rank $\geq n-1$.  \vspace{-0.15cm}
\item The linear span of the surface $\mathcal{L}^{-1}$ is the ambient projective space $\mathbb{P}(\mathbb{S}^n)$.\vspace{-0.15cm}
\item $\mathcal{L}^{-1}$ is linearly isomorphic to the $(n-1)$st Veronese embedding of $\mathbb{P}^2$.
 \vspace{-0.15cm}
\item The Chow matrix of $\mathcal{L}$ is invertible.
\end{enumerate}
\end{theorem}

To explain this result, we must first define the {\em Chow matrix}.
Fix a basis $X,Y,Z$  of $\mathcal{L}$, and let $x,y,z$ be unknowns.
The Chow matrix is
square of size $\binom{n+1}{2} \times \binom{n+1}{2}$, and its
entries are homogeneous polynomials of degree $n-1$ in
the entries of $X,Y,Z$. The columns are labeled by the monomials
of degree $n-1$ in $x,y,z$, and the rows are labeled by
pairs $(i,j)$ where $1 \leq i \leq j \leq n$. Consider the entry in position
$(i,j)$  of the adjoint matrix of $xX+yY+zZ$. This is a polynomial of
degree $n-1$ in $x,y,z$. Its coefficients form the
row indexed $(i,j)$ of the Chow matrix. The entry in a given column is
the coefficient of the monomial label.

\begin{example}[$n=3$]
The Chow matrix of three plane conics has format $6 \times 6$: \vspace{-0.15cm}
$$
\begin{tiny}
\begin{pmatrix}
x_{22} x_{33}-x_{23}^2 &
 x_{22} y_{33}-2 x_{23} y_{23}+x_{33} y_{22} &
  x_{22} z_{33}-2 x_{23} z_{23}+x_{33} z_{22} &
      y_{22} y_{33}-y_{23}^2 &
  \cdots & \cdots  \phantom{spa} \\
x_{13} x_{23} -x_{12} x_{33} &x_{13} y_{23}-x_{12} y_{33}+x_{23} y_{13}-x_{33} y_{12} &
 x_{13} z_{23}-x_{12} z_{33}+x_{23} z_{13}-x_{33} z_{12} &
  y_{13} y_{23} -y_{12} y_{33}&
  \cdots & \cdots \phantom{spa} \\

 x_{12} x_{23}-x_{13} x_{22} & x_{12} y_{23}-x_{13} y_{22}-x_{22} y_{13}+x_{23} y_{12} &
  x_{12} z_{23}-x_{13} z_{22}-x_{22} z_{13}+x_{23} z_{12} & y_{12} y_{23}-y_{13} y_{22} &
  \cdots & \cdots \phantom{spa} \\
x_{11} x_{33}-x_{13}^2 & x_{11} y_{33}-2 x_{13} y_{13}+x_{33} y_{11} &
x_{11} z_{33}-2 x_{13} z_{13}+x_{33} z_{11} &
 y_{11} y_{33}-y_{13}^2 &
 \cdots & \cdots \phantom{spa} \\
x_{12} x_{13} -x_{11} x_{23} &
x_{12} y_{13} -x_{11} y_{23}+x_{13} y_{12}-x_{23} y_{11} &
x_{12} z_{13} -x_{11}  z_{23}+x_{13} z_{12}-x_{23} z_{11} &
  y_{12} y_{13} -y_{11} y_{23} &
\cdots & \cdots \phantom{spa} \\
x_{11} x_{22}-x_{12}^2 & x_{11} y_{22}-2 x_{12} y_{12}+x_{22} y_{11} &
 x_{11} z_{22}-2 x_{12} z_{12}+x_{22} z_{11} &
  y_{11} y_{22}-y_{12}^2 &
 \cdots & \cdots \phantom{spa} \\
\end{pmatrix}.
\end{tiny} \vspace{-0.13cm}
$$
The row labels $(1,1),(1,2),(1,3),(2,2),(2,3),(3,3)$ refer
to the entries of the adjoint of $xX+yY+zZ$. The
column labels are the monomials $x^2, xy, xz, y^2, yz, z^2$.
The determinant of the Chow matrix is an irreducible homogeneous
polynomial of degree $12$ with $22659$ terms.
It vanishes if and only if
the net contains a rank-$1$ matrix.
Thus, it equals the Chow
form \cite{FKO}  of the Veronese surface in~$\mathbb{P}^5$.
\end{example}

The {\em Chow form} of a projective variety is
a hypersurface in the Grassmannian. Its points are
linear spaces of complementary dimension that
unexpectedly intersect the variety. If the linear
space is expressed as the row space of a matrix
then the entries of the matrix are the
{\em dual Stiefel coordinates} of that linear space.

\begin{lemm} \label{lem:stiefel} Consider the net $\mathcal{L} = \{xX + yY+zZ\}$ spanned by
$X,Y,Z \in \mathbb{S}^n$.
The Chow matrix is singular if and only if $\mathcal{L}\setminus\{{\bf 0}_n\}$
contains a matrix of rank $\leq n-2$.
So the determinant of the Chow matrix is the Chow form,
written in dual Stiefel coordinates, of the subvariety of $\mathbb{P}(\mathbb{S}^n)$ defined by the
$(n-1) \times (n-1)$-minors.
\end{lemm}

\begin{proof}
If $X$ has rank $\leq n-2$ then the coefficient of $x^{n-1}$ is zero
in every entry of the adjoint of $xX + yY+zZ$, so the first row
of the Chow matrix is zero. Such a rank drop happens whenever
some matrix in $\mathcal{L}$
has rank $\leq n-2$ since the determinant of the Chow matrix is invariant of the
choice of basis in $\mathcal{L}$. That determinant is a nonzero
homogeneous polynomial
in the entries of $X,Y,Z$ of degree
$(n-1)\binom{n+1}{2}$. We already argued that it vanishes on the hypersurface defined
by the Chow form of the rank-$( n-2)$ variety.
That variety has degree $\binom{n+2}{3}$, so its Chow form
 in dual Stiefel coordinates is an irreducible
polynomial of degree $3 \binom{n+2}{3}$. This irreducible polynomial
divides our determinant. Since $(n-1)\binom{n+1}{2} = 3 \binom{n+2}{3}$,
the two polynomials agree up to a nonzero constant. \end{proof}

\begin{remark}
The Chow matrix is new and interesting even for $n=4$.
It gives the Chow form for
symmetric $4 \times 4$ matrices of rank~$\leq~2$.
A formula for that Chow form was the
main result in \cite[Section 3]{FKO}. Our construction is much simpler.
 \end{remark}

\begin{proof}[Proof of Theorem \ref{thm:chow}]
Lemma \ref{lem:stiefel} says that  conditions (1) and (4) are equivalent.
The surface $\mathcal{L}^{-1}$ is parametrized by the entries
of the adjoint of $xX+yY+zZ$.
By inverting the Chow matrix, we can express every monomial
of degree $n-1$ in $x,y,z$ as a linear combination of those entries.
Hence $\mathcal{L}^{-1}$ is linearly isomorphic to the Veronese surface,
so (4) implies  (3). Clearly (3) implies (2). Finally,
 if (4) fails then $\mathcal{L}^{-1}$ lies in a hyperplane in $\mathbb{P}(\mathbb{S}^n)$, so (2) fails.
\end{proof}

\begin{corollary} \label{cor:chowrank}
Given any regular net $\mathcal{L}$ of symmetric $n \times n$ matrices,
the rank of its Chow matrix
equals the dimension of the linear span of $\mathcal{L}^{-1}$ inside $\mathbb{S}^n$.
\end{corollary}

\begin{proof}
A linear form vanishes on $\mathcal{L}^{-1}$ if and only if it vanishes
on the adjoint of $xX + yY+zZ$. These are the
 vectors in the left kernel of the Chow matrix.
\end{proof}

Every regular subspace $\mathcal{L}\subset\mathbb{S}^n$ generates a unique Jordan subalgebra
${\rm Jor}(\mathcal{L})$ inside $\mathbb{S}^n$.
The dimension of ${\rm Jor}(\mathcal{L})$ is bounded below by the
rank of the Chow matrix of $\mathcal{L}$.
However, this bound is not tight. To see this, consider the~net
\begin{equation}
\label{eq:netrank8}
\mathcal{L} \,\,=\,\,\, \begin{small} \begin{pmatrix}
\, x +y \, & z & z & 0 \,\, \\
z & \!\! x-y \!\! & 0 & z \,\, \\
z & 0 & x & z \,\,\\
0 & z & z & x \,\, \\
\end{pmatrix}
\end{small}.
\end{equation}
Here, ${\rm Jor}(\mathcal{L})= \mathbb{S}^4$ but
 the $10 \times 10$ Chow matrix has rank~$8$.  Its left kernel
 gives the  two linear forms
 $\, z_{14}-z_{23}-z_{33}+z_{44} \,$ and $\, 2 z_{12}-z_{13}-z_{24}  \,$ which
 vanish on $\mathcal{L}^{-1}$.

\smallskip

We now return to the polynomial equations that vanish on the Jordan locus
in the Grassmannian.
The following result is immediate from Corollary \ref{cor:chowrank}.

\begin{proposition} \label{prop:choweqns} Using dual Stiefel coordinates
on ${\rm Gr}(m,\mathbb{S}^n)$, the Jordan locus  $\,{\rm Jo}(m,\mathbb{S}^n)\,$
is cut out by the $(m+1) \times (m+1)$-minors of the Chow matrix.
\end{proposition}
\begin{proof}
For any regular $\mathcal{L}\in\Gr(m,\mathbb{S}^n)$, the reciprocal variety $\mathcal{L}^{-1}$ is irreducible of dimension $m$ and contained in its span, which is itself an irreducible variety. These two varieties coincide if and only if the latter has dimension~$\leq m$.
\end{proof}

\begin{remark}
These $4 \times 4$-minors do not generate
the prime ideal of ${\rm Jo}(3,\mathbb{S}^n)$. For instance, if
$n=3$ and we fix ${\bf 1}_3 \in \mathcal{L}$ then
Proposition \ref{prop:choweqns} yields $30$ quartics in $x_{ij},y_{ij}$.
These cut out the same variety
${\rm Jo}_{\bf 1}(3,\mathbb{S}^3)$ as the three quadrics in (\ref{eq:buchloe}).
\end{remark}

We know from Theorem \ref{thm:classification m=3,n=4}
that  ${\rm Jo}(3,\mathbb{S}^4)$ has two irreducible components,
of types 1(a) and 1(b).
They are given by the nets $\mathcal{L}_1$ and $\mathcal{L}_2$ in the next example.

\begin{example} \label{ex:drei}
We consider the congruence orbits of the following nets in $\mathbb{S}^4$:  \vspace{-0.1cm}
$$\mathcal{L}_1 \,\,= \,\,
\begin{small} \begin{pmatrix}
x & 0 & 0 & 0 \\
0 & x & 0 & 0 \\
0 & 0 & y & 0 \\
0 & 0 & 0 & z
\end{pmatrix} \end{small}\, , \,\,\,
\mathcal{L}_2 \,\,= \,\,
\begin{small} \begin{pmatrix}
x & y & 0 & 0 \\
y & z & 0 & 0 \\
0 & 0 & x & y \\
0 & 0 & y & z
\end{pmatrix} \end{small}
\, , \,\,\,
\mathcal{L}_3 \,\,= \,\,
\begin{small} \begin{pmatrix}
0 & 0 & x & y \\
0 & -2x & -y & z \\
x & -y & -2z & 0 \\
y & z & 0 & 0
\end{pmatrix} \end{small}.  \vspace{-0.1cm}
$$
They  have codimensions $9$, $10$ and $9$ in ${\rm Gr}(3,\mathbb{S}^4)$.
The first two fill out ${\rm Jo}(3,\mathbb{S}^4)$.

It is interesting to  compare the non-diagonalizable nets
$\mathcal{L}_2$ and $\mathcal{L}_3$.
Both represent double conics in $\mathbb{P}^2$, since
 ${\rm det}(\mathcal{L}_2) =  {\rm det}(\mathcal{L}_3) = (xz-y^2)^2$.
 But, while $\mathcal{L}_2^{-1}$ is linear,
  the Chow matrix of $\mathcal{L}_3$ is invertible, so
 $\mathcal{L}_3^{-1}$ is a Veronese surface.
 The congruence orbit of $\mathcal{L}_3$ is the
 smooth variety of codimension $9$
  studied~in \cite{EPS}. It is essentially the
 Hilbert scheme of
  twisted cubic curves in $\mathbb{P}^3$.
  \end{example}

We end this section by reporting a challenging computation.
It is predicated on the idea that
Pl\"ucker coordinates are best when working with Grassmannians.

\begin{example}\label{example: quadrics in plucker}
The $21$-dimensional Grassmannian ${\rm Gr}(3,\mathbb{S}^4)$ lives in $\mathbb{P}^{119}$.
This space has
 $\binom{10}{3} = 120$ dual Pl\"ucker coordinates $p_{ijk}$,
 one for each  $3 \times 3$-minor~of  \vspace{-0.14cm}
$$ \begin{small} \begin{pmatrix}
\, u_{11} & u_{12} & u_{13} & u_{14} & u_{22} & u_{23} & u_{24} & u_{33} & u_{34} & u_{44} \, \\
\, x_{11} & x_{12} & x_{13} & x_{14} & x_{22} & x_{23} & x_{24} & x_{33} & x_{34} & x_{44} \, \\
\, y_{11} & y_{12} & y_{13} & y_{14} & y_{22} & y_{23} & y_{24} & y_{33} & y_{34} & y_{44}\  \\
\end{pmatrix}. \end{small}  \vspace{-0.14cm}
$$
We computed the quadrics that vanish on subvarieties of ${\rm Gr}(3,\mathbb{S}^4)$.
Modulo Pl\"ucker relations, there are $4950$ linearly independent quadrics.
Each quadric has the form $\sum \lambda^{ijk}_{lmn} p_{ijk} p_{lmn}$
where the $\lambda^{ijk}_{lmn}$ are unknown coefficients and
$i \leq l, j \leq m, k \leq n$.
For each of the three orbits in
Example \ref{ex:drei},
we set up a linear system of equations in the $\lambda^{ijk}_{lmn}$
 whose solutions  are quadrics that vanish on that variety.
We then solved these equations using {\tt Maple}. Here are the results:

We found a basis of $85$ quadrics for the orbit of $\mathcal{L}_1$. One representative is
\vspace{-0.16cm}
$$ \begin{small}
\begin{matrix} p_{124} p_{789}-p_{125} p_{589}+p_{125} p_{679}
   -p_{126} p_{678} + p_{127} p_{569} - p_{128} p_{568} - p_{234} p_{579} \\
      + p_{234} p_{678} + p_{235}  p_{479} + p_{235} p_{568}
    - p_{236} p_{478} - p_{236} p_{567} - p_{237} p_{459} + p_{238} p_{458}. \end{matrix}
    \end{small} \vspace{-0.16cm} $$
We found a basis for $189$ quadrics for the orbit of $\mathcal{L}_2$. One representative is \vspace{-0.16cm}
\begin{equation} \label{eq:onerep} \begin{small} \begin{matrix}
 p_{012} p_{018} - p_{012} p_{026} - p_{012} p_{035} - 2 p_{012} p_{123}
     - p_{013} p_{017} + 2 p_{013} p_{025} - p_{023} p_{024}.  \end{matrix}
     \end{small} \vspace{-0.16cm}
\end{equation}
Each of these $189$ quadrics also vanishes on $\mathcal{L}_3$. But, there are
$918$ quadrics vanishing on the orbit of $\mathcal{L}_3$. One representative is
$\, 8 p_{012} p_{457} \,-\, 4 p_{045} p_{057} \,+\, p_{047}^2 $.
To see that $\mathcal{L}_3$ is not contained in $\mathcal{L}_2$,
we can use the quartics in  Proposition~\ref{prop:choweqns}.
\end{example}

\begin{small}

Rings and Algebras (math.RA); Algebraic Geometry (math.AG)

Jordan Algebras of Symmetric Matrices


Optimal supply against fluctuating demand

Abstract: Sornette et al. claimed that the optimal supply does not agree with the average demand, by analyzing a bakery model where a daily demand fluctuates with a uniform distribution. In this note, we extend the model to general probability distributions, and obtain the formula of the optimal supply for Gaussian distribution, which is more realistic. Our result is useful in a real market to earn the largest income on average.

\section{Introduction}

Sornette et al.\ (1999) claimed that the optimal supply does
not agree with the average demand, contrary to the common sense in
economy. They considered a bakery model where a daily demand
fluctuates with a uniform distribution, and derived the formula of the
optimal supply. Although their result is reasonable and
meaningful, it is not clear how the result will change if we consider
a different distribution.
In this note, we extend the model to general probability
distributions, and calculate the optimal supply for the
Gaussian distribution, which is more realistic in a market.

\section{Model and analysis of Sornette, Stauffer and Takayasu}

In this section we review the model and analysis of Sornette
et al.\  Let us consider a bakery shop where baked croissants are sold
every day. A question is how many croissants should be baked a day to
make the maximal profit.
We define the variables as follows.
\begin{itemize}
\item $x$: the selling price of a croissant.
\item $y$: its production cost.
\item $s$: the production number of croissants per day (supply).
\item $n$: the number of croissants requested by customers per day
(demand).
\item $D$: the average demand, i.e., $D\equiv<n>$.
\end{itemize}
The expectation of the total profit $L(s)$ is given by
\begin{equation}\label{L}
L(s)\equiv<x~{\rm min}(n,s)-ys>
=x\int^s_0nP(n)dn+xs\int^{\infty}_sP(n)dn-ys,
\end{equation}
where $P(n)$ is the probability distribution of $n$.

Sornette et al.\ assumed, for simplicity, a uniform distribution,
\begin{equation}\label{uniform}
P_u(n)\equiv
\left\{\begin{array}{ll}
1/2\delta ~& {\rm for} ~~~D-\delta\le n\le D+\delta\\
0         ~& {\rm for} ~~~n<D-\delta,~D+\delta<n.
\end{array}\right.
\end{equation}
Then one can integrate (\ref{L}) as
\begin{equation}
L(s)=-{x\over4\delta}\left\{s-D-\delta\left(1-{2y\over
x}\right)\right\}^2
+(x-y)\left(D-{\delta y\over x}\right).
\end{equation}
$L(s)$ takes the maximam value when $s$ takes
\begin{equation}\label{su}
s_{{\rm max}}\equiv D+\delta\left(1-{2y\over x}\right).
\end{equation}
This shows that, if the cost per price, $y/x$, is larger (smaller)
than half, the optimal demand, $s_{{\rm max}}$, is smaller (larger)
than the average demand, $D$.

\section{Optimal supply for Gaussian distribution}

Let us re-analyze (\ref{L}) for general probability distributions. We
do not have to integrate (\ref{L}) directly, because what we want to
know is the optimal supply $s_{{\rm max}}$, which is given by
\begin{equation}\label{dL}
{dL\over ds}(s_{{\rm max}})=x\int^{\infty}_{s_{{\rm max}}}P(n)dn-y=0.
\end{equation}
This simple equation gives the optimal supply for general probability
distributions.

If we assume Gaussian distribution,
\begin{equation}
P_G(n)\equiv{1\over\sqrt{2\pi}\sigma}
\exp\left[-{(n-D)^2\over2\sigma^2}\right],
\end{equation}
the above integration is expressed as
\begin{equation}
\int^{\infty}_{s_{{\rm max}}}P_G(n)dn=\frac12-\frac12
{\rm Erf}\left({s_{{\rm max}}-D\over\sqrt{2}\sigma}\right),
\end{equation}
where Erf is the error function, which is defined as
\begin{equation}
{\rm Erf}~z\equiv{2\over\sqrt{\pi}}\int^z_0e^{-t^2}dt.
\end{equation}
Then we arrive at the formula of the optimal supply for the Gaussian
distribution,
\begin{eqnarray}\label{sG}
{s_{{\rm max}}-D\over\sigma}
&=&\sqrt{2}{\rm Erf}^{-1}\left(1-\frac{2y}x\right) \nonumber\\
&=&\sqrt{{\pi\over2}}\left(1-\frac{2y}x\right)
+{\sqrt{2}\pi^{\frac32}\over24}\left(1-\frac{2y}x\right)^3
+O\left[\left(1-\frac{2y}x\right)^5\right]
~~~({\rm Gaussian}).
\end{eqnarray}
Because the cost is usually in the range $0.3x<y<0.7x$, which reads
$|1-2y/x|<0.4$, the first-order approximation in
(\ref{sG}) is sufficient in most cases.

For reference, we rewrite the result for the uniform distribution
(\ref{su}). Because the variance of the uniform distribution (\ref{uniform}) is
evaluated as $\sigma^2=\delta^2/3$, (\ref{su}) is rewritten as
\begin{equation}\label{su2}
{s_{{\rm max}}-D\over\sigma}=\sqrt{3}\left(1-\frac{2y}x\right)
~~~({\rm uniform}).
\end{equation}

We see that the difference between $\sqrt{\pi/2}\approx1.25$ in
(\ref{sG}) and $\sqrt{3}\approx1.73$ in (\ref{su2}) is not negligible.
Contrary to the speculation of
Sornette {\it et al.}, however, the critical value of the
cost-to-price ratio, $y/x=1/2$, is unchanged.
Because Gaussian distribution is more realistic,
our simple formula (\ref{sG}) is useful
in a real market to earn the largest income on average.

Physics and Society (physics.soc-ph); General Finance (q-fin.GN)


Minimal Area of a Voronoi Cell in a Packing of Unit Circles

Abstract: We present a new self-contained proof of the well-known fact that the minimal area of a Voronoi cell in a unit circle packing is equal to $2\sqrt{3}$, and the minimum is achieved only on a perfect hexagon. The proof is short and, in our opinion, instructive.

\section{Introduction}

This work originated from attempts to rely on the proof in \cite{H} of the fact that the minimal area of a Voronoi cell in a unit circle packing in $\bbR^2$ is equal to $2\sqrt{3}$. Unfortunately, this proof contains some gaps (see Appendix), and our efforts to recover the proof resulted in an alternative approach presented below.

The result itself immediately implies the strong/local version of the theorem of Thue and Fejes T\'oth \cite{Th, To} on dense unit circle packings in $\bbR^2$. Additional applications include the derivation of so-called Peierls estimates in two-dimensional lattice hard-core models of statistical mechanics \cite{MSS1}.

\section{Definitions and basic properties}

For any $\bx \in \bbR^2$ and $r > 0$ denote by $B_r(\bx)$ an open disk of radius $r$ centered at $\bx$. A shorthand notation $B(\bx)$ is used for the {\it unit} disk having $r=1$. A collection of non-overlapping open unit disks is called {\it admissible} and is denoted by $\{B(\bx_i)\}$. An admissible collection $\{B(\bx_i)\}$ represents a {\it unit circle packing}. (It is clear that an admissible collection is finite or countable, as inside each disk there exists a point with rational coordinates.) An admissible collection $\{B(\bx_i)\}$ can be identified with the corresponding collection of centers $\{\bx_i\}$ which is also called admissible. Each element of $\{\bx_i\}$ is called an {\it occupied} point in $\bbR^2$; respectively, we speak about admissible collections of occupied points. Clearly, $|\bx_{i'}-\bx_{i''}| \ge 2$ for any distinct $\bx_{i'}, \bx_{i''} \in \{\bx_i\}$, where $|\bx_{i'}-\bx_{i''}|$ denotes the Euclidean distance between $\bx_{i'},\bx_{i''} \in \bbR^2$.

For each $\bx_{i'} \in \{\bx_i\}$ the corresponding {\it Voronoi cell} is defined as
$$V(\bx_{i'}) = \left\{\bz \in \bbR^2:\;\;|\bx_{i'}-\bz| \le \inf_{i'' \not = i'}|\bx_{i''}-\bz|\right\}.$$
If for $\bx_{i''}$ the intersection $V(\bx_{i'}) \cap  V(\bx_{i''})$ contains more than one point of $\bbR^2$ then we say that $\bx_{i''}$ is a {\it
neighboring occupied point} for $\bx_{i'}$ or simply $\bx_{i''}$ is a {\it neighbor} of $\bx_{i'}$.

Observe that $V(\bx)$ can be unbounded. If $V(\bx)$ is bounded, i.e., $V(\bx) \subset B_{r}(\bx)$ for some $0 < r <\infty$, then $B(\by) \subset B_{4r}(\bx)$ for any neighboring occupied point $\by$. Due to the admissibility requirement the number of such $\by$ cannot exceed $|B_{4r}(\cdot)|/|B(\cdot)| < \infty$, where $|\cdot|$ denotes the area of the corresponding disk.

Our aim is to find the minimal possible area of a Voronoi cell among all Voronoi cells in all admissible collections of occupied points. Suppose that a Voronoi cell with minimal or close to minimal area can be found in an admissible collection containing unbounded Voronoi cells. Then this collection can be completed without breaking the admissibility by a finite or countable set of additional occupied points such that the resulting admissible collection has bounded Voronoi cells only. For that reason from now on we mainly consider admissible collections without unbounded Voronoi cells.

The rest of this section verifies some standard properties of Voronoi cells which makes the entire argument self-contained.

\medskip
{\bf Lemma 1.} {\sl For an occupied point $\bx$ in an admissible collection the corresponding Voronoi cell $V(\bx)$ contains the closure of $B(\bx)$.}

\medskip
{\bf Proof.} If for a point $\bz \in B(\bx)$ there exists an occupied point $\by \not = \bx$ with $|\by-\bz| < |\bx-\bz| < 1$ then by the triangle inequality $|\by-\bx| < 2$ which contradicts the admissibility of the collection. \qed

\medskip
{\bf Lemma 2.} {\sl For an occupied point $\bx$ in an admissible collection the corresponding Voronoi cell $V(\bx)$ is a convex subset of $\bbR^2$.}

\medskip
{\bf Proof.} For an occupied point $\by \not = \bx$ and any $\bz \in V(\bx)$ one has $|\bz-\bx|^2 \le |\bz-\by|^2$ and consequently
$$|\bx|^2 - 2\bz\cdot \bx \le |\by|^2 - 2\bz\cdot \by.$$ Consider two distinct points $\bz', \bz'' \in V(\bx)$ and suppose that $\bz = \lambda \bz' + (1-\lambda) \bz''$, where $0 \le \lambda \le 1$. Then
$$\beacl
|\bz -\bx|^2 &= |\bz|^2 + \lambda(|\bx|^2-2 \bz'\cdot \bx) + (1-\lambda)(|\bx|^2-2 \bz''\cdot \bx)\cr \\
&\le |\bz|^2 + \lambda(|\by|^2-2 \bz'\cdot \by) + (1-\lambda)(|\by|^2-2 \bz''\cdot \by) \cr
&= |\bz -\by|^2,\ena
$$
which establishes the lemma. \qed

\medskip
{\bf Lemma 3.} {\sl For an occupied point $\bx$ in an admissible collection the boundary of the corresponding Voronoi cell $V(\bx)$ is piecewise linear.}

\medskip
{\bf Proof.} Consider two different occupied points $\bx, \by$ with $V(\bx)\cap V(\by)$ containing more than one $\bbR^2$ point. Let $\bz', \bz'' \in V(\bx)\cap V(\by)$ and $\bz' \not= \bz''$. Now suppose that $\bz = \lambda \bz' + (1-\lambda) \bz''$, $0 \le \lambda \le 1$. Then
$$\beacl
|\bz -\bx|^2 &= |\bz|^2 + \lambda(|\bx|^2-2 \bz'\cdot \bx) + (1-\lambda)(|\bx|^2-2 \bz''\cdot \bx)\cr
&= |\bz|^2 + \lambda(|\by|^2-2 \bz'\cdot \by) + (1-\lambda)(|\by|^2-2 \bz''\cdot \by) \cr
&= |\bz -\by|^2,\ena
$$
i.e., $\bz \in V(\bx)\cap V(\by)$. Also note that for any two distinct occupied points $\bx, \by$ the set
$$\left\{ \bz \in \bbR^2:\;\; |\bz -\bx| \le |\bz -\by| \right\}$$
is a closed half-plane. Consequently, $V(\bx)$ is an intersection of a finite or countable number of closed half-planes. \qed

\medskip
Lemmas~1-3 are valid for both bounded and unbounded Voronoi cells in an admissible collection of occupied points. If $V(\bx)$ is bounded then these lemmas imply that $V(\bx)$ is a convex polygon containing the unit disk $B(\bx)$. Furthermore, each neighboring occupied point $\by$ is the reflection of $\bx$ with respect to the common side of $V(\bx)$ and $V(\by)$.

The clockwise circular order of sides of a bounded polygon $V(\bx)$ generates the clockwise circular order of the corresponding neighboring occupied points. Connecting the neighboring points in this circular order we obtain a so-called {\it polygon of neighbors}. Together with $\bx$, each side of this polygon uniquely defines a {\it constituting triangle} such that the entire polygon of neighbors is partitioned into constituting triangles. By construction, the vertices of $V(\bx)$ are the centers of the circumcircles of these constituting triangles. The angle of the constituting triangle at vertex $\bx$ is called the {\it constituting angle}.

\medskip
{\bf Lemma 4.} {\sl The circumradius of a constituting triangle is not shorter than $2/\sqrt{3}$.}

\medskip
{\bf Proof.} At least one angle of a constituting triangle, say angle $\alpha$, is not larger than $\pi / 3$. By the admissibility requirement the length $a$ of the opposite triangle side is not shorter than $2$. According to the sine theorem for triangles the triangle circumradius
$$r = {a \over 2 \sin \alpha} \le {2 \over 2 \sin {\pi \over 3}} = {2\over\sqrt{3}},$$
which establishes the lemma. \qed

\section{Results}

Consider an admissible collection $\{\bx_i\}$ which forms a triangular lattice with the shortest distance between sites equal to  $2$. Then for each $\bx_i$ the corresponding polygon of neighbors is a perfect hexagon with the side length equal to $2$. Correspondingly, $V(\bx_i)$ is a perfect hexagon with the side length equal to $2 / \sqrt{3}$ and the area $|V(\bx_i)|= 2 \sqrt{3}$.

\medskip
{\bf Theorem.} {\sl The minimal area of a Voronoi cell in an admissible configuration of occupied points is equal to $2 \sqrt{3}$. The minimum is achieved only on occupied points $\bx$ having a perfect hexagon with side length $2/\sqrt{3}$ as its Voronoi cell, or equivalently, a perfect hexagon with side length} $2$ as the corresponding polygon of neighbors.

\medskip
{\bf Proof.} In view of properties of a Voronoi cell presented in Lemmas~1-4 the problem is reduced to finding the polygon of minimal area among all polygons $P$ having the following properties:
\begin{description}
\item{(i)} $P$ is convex.
\item{(ii)} $P$ contains a unit disk.
\item{(iii)} The distances from the vertices of $P$ to the center of this disk are not shorter than~${2 \over \sqrt{3}}$.\end{description}

\noindent
The last property is a weaker replacement of the requirement for a Voronoi cell to have the corresponding neighboring occupied points at distances not shorter than  $2$ from each other. It turns out that this weaker requirement is enough to establish the theorem.

For a convex polygon of area $a$ and one of its angles of measure $\alpha$ define the {\it angular density} (with respect to this angle) as the ratio $E={a \over \alpha}$. Now take any polygon satisfying (i)-(iii) and consider several rays originating at the center $\bo$ of the contained unit disk.
Let the angles between the clockwise consecutive rays be
smaller than
$\pi$. These rays partition the polygon into several convex polygons each located inside
the angle $\alpha_j$ between the corresponding two clockwise consecutive rays. Clearly, the area $a$ of the original polygon $P$ can be calculated as
$a= \sum_j E_j \alpha_j$, where $E_j$ is the corresponding angular density (with respect to angle $\alpha_j$).

For any vertex $\bv$ of $P$
the convex hull of this vertex and the unit disk centered at $\bo$ belongs to $P$.
If $|\bo-\bv|=r$ then this convex hull is bounded by two straight segments $[\bv\ba]$ and $[\bv\bc]$ of length $\sqrt{r^2-1}$ and the arc of the unit circle connecting $\ba$ with $\bc$ and having length $2\pi - 2 \arctan\sqrt{r^2-1}$. The area of the quadrilateral $[\bo\ba\bv\bc]$ is equal to $\sqrt{r^2-1}$ and $|\angle \ba\bo\bc|=2\arctan\sqrt{r^2-1}$. Therefore, the corresponding angular density is $E_{\bv} = {\sqrt{r^2-1}\over 2 \arctan\sqrt{r^2-1}}$. Take any point $\bb$ inside the segment $[\bv\bc]$ and consider the angular density $\overline E_{\bv}$ of the smaller quadrilateral $[\bo\ba\bb\bv]$ with respect to the angle $\angle \ba\bo\bb$. Let $|\bb-\bc| = x < \sqrt{r^2-1}$. Then
$$\overline E_{\bv} = {\sqrt{r^2-1} - {x \over 2} \over 2\arctan \sqrt{r^2-1} - \arctan x} > {\sqrt{r^2-1} \over 2\arctan \sqrt{r^2-1}} = E_{\bv}.$$

For two clockwise consecutive polygon vertices $\bv_i$ and $\bv_{i+1}$ consider the corresponding convex hulls and two corresponding quadrilaterals $[\bo\ba_i\bv_i\bc_i]$ and $[\bo\ba_{i+1}\bv_{i+1}\bc_{i+1}]$. Observe that consecutive open quadrilaterals are either adjacent (have a common side) or intersecting.  (In the case when they are separated by
a non-zero angle $\angle \bc_i\bo\ba_{i+1}$ the segment $[\bv_i \bv_{i+1}]$ intersects the interior of the unit disk which contradicts property (ii).) With this observation at hand, denote by $\bb_i$ the intersection point of segments $[\bv_i\bc_i]$ and $[\bv_{i+1}\ba_{i+1}]$. According to the displayed equation above, the angular density of $[\bo\ba_i\bv_i\bb_i]$ is larger than the angular density of $[\bo\ba_i\bv_i\bc_i]$ and the angular density of $[\bo\bb_i\bv_{i+1}\bc_{i+1}]$ is larger than the angular density of $[\bo\ba_{i+1}\bv_{i+1}\bc_{i+1}]$.

\def\cQ{\mathcal Q}

Consider now the union of the quadrilaterals $[\bo\ba_i\bv_i\bc_i]$ over all vertices $\bv_i$. It is a polygon $Q$ (generally, non-convex) contained in $P$. Between each two vertices $\bv_i$ and $\bv_{i+1}$ the polygon $Q$ may contain an additional vertex $\bb_i$ that is the intersection point introduced above. The angular density of $[\bo\bb_{i-1}\bv_i\bb_i]$ is not smaller than the angular density of $[\bo\ba_i\bv_i\bc_i]$, and the two angular densities are equal only if $\bb_{i-1}=\ba_i$ and $\bb_i=\bc_i$. For $r \ge {2 \over \sqrt{3}}$ the minimum of $E_{\bv} = {\sqrt{r^2-1}\over 2\arctan\sqrt{r^2-1}}$ equals $\sqrt{3} \over \pi$ and is achieved only at $r = {2 \over \sqrt{3}}$. Thus, the total area of the union of the quadrilaterals $[\bo\ba_i\bv_i\bc_i]$ (or equivalently the union of mutually disjoint open quadrilaterals $[\bo\bb_{i-1}\bv_i\bb_{i}]$) is not smaller than ${\sqrt{3} \over \pi} 2\pi = 2\sqrt{3}$. Obviously, this minimum is achieved only when the interiors of the quadrilaterals $[\bo\ba_i\bv_i\bc_i]$ are disjoint (i.e., $[\bo\ba_i\bv_i\bc_i] = [\bo\bb_{i-1}\bv_i\bb_{i}]$) and $|\bo-\bv_i| = {2 \over \sqrt{3}}$ for all $i$. In this case the corresponding polygon is the perfect hexagon with side length ${2 \over \sqrt{3}}$.

Indeed, if $V(\bx)$ is a hexagon with $r_i > {2 \over \sqrt{3}}$ for some $i$ then $|\angle \ba_i \bo \bc_i| > {2 \pi \over 6}$. Therefore, such a hexagon has the angular density larger than minimal for some non-zero angle. Consequently, its area is larger than $2\sqrt{3}$.

Any polygon with $n > 6$ vertices contains at least two overlapping quadrilaterals $[\bo \ba_i\bv_i\bc_i]$ even if $r_i = {2 \over \sqrt{3}}$ for all $i$ because $n {2 \pi \over 6} > 2\pi$. Therefore, the angular density becomes larger than minimal for some non-zero angle. Hence, any polygon satisfying (i)-(iii) with more than 6 vertices has a non-minimal area.

A polygon with $n = 3,4$ or 5 vertices necessarily has at least one $|\angle \ba_i \bo \bc_i| \ge {2\pi \over n}$ and therefore $r_i > {1 \over \cos{\pi \over n}} > {2 \over \sqrt{3}}$. Thus, the corresponding angular density $E_{\bv_i}$ is again greater than $\sqrt{3} \over \pi$ for some non-zero angle; consequently, the area of the entire polygon is larger than the minimal area $2\sqrt{3}$.~\qed

\section{Appendix}

The problem with the argument in \cite{H} is that at some point the proof considers the ``most critical case'', but it is not specified what quantity is optimized at this ``most critical case'' and why this quantity is important. The desired quantity seems to be the area excess in $V(\bx)$ over the area of the contained unit disk which can be attributed to a single non-close neighbor. (In the terminology of \cite{H} a neighbor $\by$ is non-close to the center $\bx$ of $V(\bx)$ iff $|\bx-\by| > 2.3$.) By considering this ``most critical case'' for a single non-close neighbor the proof in \cite{H} concludes that the corresponding area excess is at least $0.21$. After that the proof claims that the presence of two non-close neighbors implies that the area excess is at least $0.42$. This additivity assumption is actually wrong as individual excesses can overlap, and one can give a counterexample showing two non-close neighbors with the total attributed area excess smaller than $0.42$.