\section{Introduction}

The celebrated Eckart-Young Theorem says that, for a real $m \times
n$-matrix $A$ with $m \leq n$ and for an integer $k \leq m$, a matrix
$B$ of rank at most $k$ nearest to $A$ is obtained from $A$ as follows:
Compute the singular value decomposition $A=U \Sigma V^T$, where $U,V$
are orthogonal matrices and where $\Sigma=\mathrm{diag}(\sigma_1,\ldots,\sigma_m)$
is the ``diagonal'' $m \times n$-matrix with the singular values $\sigma_1
\geq \cdots \geq \sigma_m \geq 0$ on its main diagonal, and set $B:=U
\mathrm{diag}(\sigma_1,\ldots,\sigma_k,0,\ldots,0) V^T$. Such a best rank-$k$
approximation is unique if and only if $\sigma_k > \sigma_{k+1}$,
and for us ``nearest'' refers to the Frobenius norm (but in fact,
the result holds for arbitrary $\O_m \times \O_n$-invariant norms
\cite{Mirsky}).

For higher-order tensors, an analogous approach for finding best rank-$k$
approximations fails in general \cite{VNVM}. It succeeds, with respect
to the Frobenius norm, for {\em orthogonally decomposable tensors}
\cite{VNVM,BDHR}, but this is a very low-dimensional real-algebraic
variety in the space of all tensors.  In this paper, we will establish
versions of the Eckart-Young Theorem and the Spectral Theorem that {\em
do} hold for general tensors.

To motivate this theorem, consider matrices once again, and
assume that the $\sigma_i$ are distinct and positive. A statement
generalizing the Eckart-Young Theorem says that we obtain all critical
points of the distance function $d_A(B):=||A-B||^2$ on the manifold of
rank-$k$ matrices by setting any $m-k$ of the singular values equal to
zero \cite{DHOST}, so as to obtain a matrix
\[ B_{i_1,\ldots,i_k}:=U
\mathrm{diag}(0,\ldots,0,\sigma_{i_1},0,\ldots,0,\sigma_{i_2},0,
\ldots,0, \sigma_{i_k},0,\ldots) V^T \]
for any ordered $k$-tuple $i_1<\ldots<i_k$ in $\{1,\ldots,m\}$. We will
call these critical points {\em critical rank-$k$ matrices for $A$}.
In particular, the critical rank-one matrices are $B_1,\ldots,B_m$, and
we draw attention to the fact that for each $k \geq 1$ and each $k$-tuple
$i_1<\ldots<i_k$ the critical rank-$k$ matrix $B_{i_1,\ldots,i_k}$ lies
in the linear span of $B_1,\ldots,B_m$. Moreover, this linear span has
a direct description in terms of $A$: it consists of all matrices $B$
such that both $A^TB$ and $AB^T$ are symmetric matrices.

Taking cue from this observation, we will associate a {\em critical space}
$H_f$ to a tensor $f$,
show that $H_f$ contains the critical \red{rank-at-most-$k$}
tensors for $f$ for each value of $k$ (see below for a definition),
and that $H_f$ is spanned by the critical rank-one tensors for $f$. We
will establish these results for sufficiently general partially symmetric tensors, and we work over the base field $\mathbb{C}$ rather than $\mathbb{R}$.

\begin{thm} \label{thm:Main}
Let $f$ be a sufficiently general tensor in $S^{d_1} \mathbb{C}^{n_1+1} \otimes
\cdots \otimes S^{d_p} \mathbb{C}^{n_p+1}$. Then for each natural number $k$,
the critical \red{rank-at-most-$k$} tensors for $f$ lie in the critical
space $H_f$.  Moreover, if for each $\ell$ with $d_\ell=1$ the {\em
triangle inequality} $n_\ell\leq\sum_{i \neq \ell} n_i $ holds, then
$\operatorname{codim} H_f=\sum_\ell \binom{n_\ell + 1}{2}$ and $H_f$ is spanned by
the critical rank-one tensors for $f$.
\red{In particular, $f$ itself lies in the linear
span of the critical rank-one tensors for $f$.}
\end{thm}

We record the following two corollaries over the real numbers.

\begin{cor} \label{cor:tensor}
If $n_1,\ldots,n_p$ satisfy the triangle inequality $n_\ell \leq \sum_{i
\neq \ell} n_i$ for each $\ell=1,\ldots,p$, then for a sufficiently
general tensor $f \in \bigotimes_{i=1}^p \mathbb{R}^{n_i+1}$ and any natural
number $k$, any real tensor of real rank at most $k$ closest to $f$
in the Frobenius norm lies in the linear span of the complex critical
rank-one tensors for $f$. In particular, $f$ itself lies in the linear
span of the complex critical rank-one tensors for $f$.
\end{cor}

\begin{cor}\label{cor:sym}
For a sufficiently general symmetric tensor $f \in S^d \mathbb{R}^{n+1}$ and any
natural number $k$, any real symmetric tensor of real symmetric rank at
most $k$ closest to $f$ in the Frobenius norm lies in the linear span of
the complex critical symmetric rank-one tensors for $f$.  In particular,
$f$ itself lies in the linear span of the complex critical rank-one
tensors for $f$.
\end{cor}

In the case of symmetric tensors, these critical rank-one tensors
correspond to the so-called {\em eigenvectors} of $f$ \cite{Mac},
while in the case of ordinary tensors, they correspond to {\em singular
vector tuples} \cite{Lim}. In the case $n=1$ of binary forms, Corollary
\ref{cor:sym} was proved in \cite{OT}. The two corollaries above can be
regarded as generalizations of the Eckart-Young Theorem and the Spectral
Theorem from matrices to tensors.

Several remarks are in order here. First, we complexify $d_f$ to
the quadratic polynomial $d_f(u):=(u-f|u-f)$, where $(.|.)$ is the
standard complex bilinear form on the space of tensors (and not a
Hermitian form). The point of doing this is that, unlike for matrices,
the critical points of this function on low-rank tensors are in general
not real anymore, even if $f$ is real. Accordingly, the critical space
$H_f$, while defined by linear equations over $\mathbb{R}$ if $f$ is real, is
taken to be the space of complex solutions to those equations. Second,
we denote the dimensions by $n+1$ rather than $n$ since we will be using
methods from projective algebraic geometry where the formulas look
more appealing in terms of the projective dimension $n$ than in the
affine dimension $n+1$. An example of this phenomenon is the triangle
inequalities in the theorem, which hold if and only if the variety dual
to the Segre-Veronese embedding of the product $\mathbb{P}^{n_1} \times
\cdots \times \mathbb{P}^{n_p}$ via degrees $d_1,\ldots,d_p$
is a hypersurface \cite[Corollary 5.11]{GKZ}.

The remainder of this paper is organized as follows. In
Section~\ref{sec:CritSpace} we define the critical space
$H_f$ for a partially symmetric tensor $f$ and prove that the
critical \red{rank-at-most-$k$} tensors for $f$ lie in $H_f$. Then, in
Section~\ref{sec:Span} we use vector bundle techniques to compute the
dimension of the space spanned by the critical rank-one tensors for $f$
and to show that this space equals $H_f$. In Section~\ref{sec:Final}
we combine these ingredients to establish the results above.

\section{The critical space of a tensor}
\label{sec:CritSpace}

\subsection*{Partially symmetric tensors and their ranks}

Let $p \in \mathbb{Z}_{\geq 1}$, let $V_1,\ldots,V_p$ be complex vector spaces,
and let $d_1,\ldots,d_p \in \mathbb{Z}_{\geq 1}$. Let $S^dV$ be the $d$-th symmetric power of $V$. We will study
tensors in the space
\[ T:=S^{d_1} V_1 \otimes \cdots \otimes S^{d_p} V_p. \]
So for $p=1$, $T$ is a symmetric power of $V_1$, which is canonically
isomorphic with the space of symmetric, $d_p$-way $n_1 \times \cdots
\times n_1$-tensors. On the other hand, when all $d_i$ are equal to $1$,
then $T$ is a space of $p$-way ordinary tensors. We will
write $[p]:=\{1,\ldots,p\}$.

Inside $T$, let $X$ be the set of all tensors of the form
\[ v_1^{d_1} \otimes \cdots \otimes v_p^{d_p}\ \quad
(v_\ell \in V_\ell, \ell \in [p]).\]
Then $X$ is a closed subvariety of $T$ known as the affine cone over the
Segre-Veronese embedding of $\mathbb{P}^{n_1} \times \cdots \times \mathbb{P}^{n_p}$
of degrees $(d_1,\ldots,d_p)$. Let $kX$ denote the set of sums of $k$
elements of $X$.  An arbitrary element $t$ of $T$ lies in $kX$ for some $k$,
and the minimal such $k$ is called the rank of $t$ \cite[Definition
5.2.1.1]{Landsberg12}. For $p=1$
this is the symmetric or Waring rank, and if all $d_q$ are $1$, this notion is
ordinary tensor rank. We write $\operatorname{Sec}_k(X)$ for the Zariski
(or Euclidean) closure of $kX$ in $T$.

For real tensors a few modifications are needed. The real rank of a real
tensor $t$ is the minimum $k$ such that $t=\sum_{i=1}^k\lambda_ix_i$
with $\lambda_i\in{\mathbb R}$ and $x_i\in X_{\mathbb{R}}$ (it is enough to allow $\lambda_i=\pm 1$). For example {$(e_1+ie_2)^3+(e_1-ie_2)^3$} has complex rank $2$ and real rank $3$. Real rank is subtle for low-rank approximation of tensors.
A classical example of de Silva and Lim \cite{dSL} shows that for almost all $2\times 2\times 2$-tensors of real rank $3$ (like the above one) does not exist a closest tensor
of real rank $2$, while such phenomena may happen only for measure zero subsets of
the set of complex tensors of given rank.

\subsection*{Symmetric bilinear forms and pairings}

If $V,W$ are complex vector spaces with symmetric bilinear forms $(.|.)$,
and if $d \in \mathbb{Z}_{\geq 0}$, then $S^d V$ and $V \otimes W$ carry unique
symmetric bilinear forms, also denoted $(.|.)$, that satisfy

\begin{align*}
(v^d|v'^d)&:=(v|v')^d\text{ and} \\
(v \otimes w | v' \otimes w') &:= (v|v') (w|w').
\end{align*}

{The first of these equalities implies

\begin{align*}
(v_1\ldots v_d|v'^d)&=\prod_{i=1}^d(v_i|v')
\end{align*}
}
and more in general
\begin{align*}
(v_1 \cdots v_d | v'_1 \cdots v'_d)&= \frac{1}{d!} \sum_{\pi \in S_d}
\prod_{i=1}^d (v_i|v'_{\pi(i)})
\end{align*}

We now fix nondegenerate symmetric bilinear forms on each $V_\ell,\
\ell \in [p]$. Then, iterating
these constructions, we obtain a canonical bilinear form on $T$.

Using the bilinear forms on $V$ and $W$, we can also build more
general bilinear maps whose outputs are vectors or tensors rather than
scalars. We will call these bilinear maps {\em pairings} and denote them
by $[.|.]$. Of particular relevance to us is the skew-symmetric pairing $S^d V \times
S^d V \to \bigwedge^2 V$
determined by

\[
[v^d | w^d]:=
 (v|w)^{d-1}
v \wedge w
\]

which implies

{
\[
[v_1\ldots v_d | w^d]=\frac{1}{d}\sum_{i'\in[d]}
\left(\prod_{i\neq i'}(v_i|w)\right)v_{i'}\wedge w
\]
}

and more in general

\[
[v_1 \cdots v_d | w_1 \cdots w_d]=
\frac{1}{d\cdot d!}\sum_{i',j' \in [d]} \sum_{\pi:[d]\setminus i' \to
[d]\setminus j'}
\left( \prod_{i \neq i'} (v_i|w_{\pi(i)}) \right)
v_{i'} \wedge w_{j'}
\]
where $\pi$ runs over all bijections $[d]\setminus i' \to [d]\setminus
j'$.

\begin{remark}
In the case of binary forms ($\dim V=2$ and arbitrary $d$), the pairing $[f|g]$
coincides (up to scalar multiples) with $\left( f| D(g)\right)$, where
	$D(g)=g_xy-g_yx$\red{;} see \cite{OT}. Note the skew-symmetry property $\left(
f| D(g)\right) = -\left( g| D(f)\right)$. On the other hand, in the case
of symmetric matrices ($d=2$ and arbitrary $V$), the pairing $[f|g]$
coincides (up to scalar multiples) with the bracket $fg-gf$.
\end{remark}

Building on this construction, for each $\ell \in [p]$ we define
a pairing $[.|.]_\ell:T \times T \to \bigwedge^2 V_\ell$ by
\begin{equation} \label{eq:Pairing}
[f_1 \otimes \cdots \otimes f_p|g_1 \otimes \cdots \otimes
g_p]_\ell:=\left(\prod_{i \neq \ell} (f_i|g_i)\right) [f_\ell|g_\ell],\
f_i,g_i \in S^{d_i} V_i.
\end{equation}
which we will use to define the critical space.

\begin{remark}\label{rem:bracket2}In the case of matrices $T=V_1\otimes V_2$, the pairing $[f,g]_1$
coincides (up to scalar multiples) with $fg^t-gf^t$,
while  $[f,g]_2$
is (up to scalar multiples) $f^tg-g^tf$.
\end{remark}

{
\subsection*{Basis, orthogonal basis and monomials}

If $B$ is a basis of $V$, then the degree-$d$ monomials in
the elements of $B$ form a basis of $S^d V$.
Such a basis is orthogonal if $B$ is orthogonal.
Hence if we fix bases (respectively, orthogonal bases) of $V_1,\ldots,V_p$, then by taking
tensor products we obtain a basis (respectively, orthogonal basis) of $T$, whose elements we
will call monomials of degree $D:=\sum_{\ell=1}^p d_\ell$.} We will use
the word {\em gcd} of two such monomials $x,y$ for the highest-degree
monomial $z$ such that both $x$ and $y$ can be obtained from $z$ by
multiplying $z$ with suitable monomials.

\begin{example}
If $p=3$ and $V_1=V_2=V_3=\mathbb{C}^3$ with the standard bilinear
form, and $d_1=d_2=3$ and $d_3=2$, then the gcd of
$(e_1^2 e_2) \otimes (e_1e_2e_3) \otimes (e_1^2)$ and
$(e_1 e_2 e_3) \otimes (e_3^3) \otimes (e_2e_3)$ equals
$(e_1 e_2) \otimes (e_3) \otimes (1)$.
\end{example}

\begin{lem} \label{lem:Pairing}
For two monomials $f=f_1 \otimes \cdots \otimes f_p,g=g_1 \otimes
\cdots \otimes g_p$ in $T$ relative to the same orthogonal bases of
$V_1,\ldots,V_p$ and for $\ell \in [p]$
we have $[f|g]_\ell=0$ unless $f_i=g_i$ for all $i \neq \ell$ and
$h:=\gcd(f_\ell,g_\ell)$ has degree $d_\ell-1$; in this case $\gcd(f,g)$
has degree $D-1$ and $[u|v]_\ell \in \mathbb{C}^* (f_\ell/h) \wedge (g_\ell/h)$.
\end{lem}

This is immediate from the definition of the pairing in
\eqref{eq:Pairing}.

\subsection*{Critical rank-one tensors}

Let $f \in T$. Then the critical points of the distance function $d_f:
x \mapsto (f-x|f-x)$ on $X$ are by definition those $x \in X \setminus
\{0\}$ for which $f-x$ is perpendicular to the tangent space $T_x X$ to
$X$ at $x$; we write this as $f-x \perp T_x X$. We call these tensors the
{\em critical rank-one tensors} for $f$. For sufficiently general $f$,
each of these critical rank-one tensors is {\em non-isotropic}, i.e.,
satisfies $(x|x) \neq 0$ (see \cite[Lemma 4.2]{DLOT}, in next Proposition
\ref{prop:RankOne} we will prove a slightly more general fact).

\comment{
Slightly
more generally, we call a non-isotropic $x$ a {\em singular rank-one
tensor} for $f$ if $f-cx \perp T_x X$ for some $c \in \mathbb{C}$.  In that
case, since $x$ is non-isotropic and since $x \in T_xX$, the $c$ with
this property is the unique $c$ with $(f-cx|x)=0$; and if $c \neq 0$,
then $cx \in X$ is a critical point of $d_f$.
}

We will establish a bilinear characterization of these critical rank-one
tensors for $f$. First, we describe the tangent space of $X$ {at a point $x$} in more
detail.
For this, write
\begin{equation} \label{eq:RankOne}
x=v_1^{d_1} \otimes \cdots \otimes v_p^{d_p}.
\end{equation}
  {Hence we may
extend each $v_\ell$ to a basis of $V_\ell$. We then obtain
an $x$-{\em adapted} basis of $T$ consisting of monomials.
If moreover $x$ is non-isotropic, we have $(v_\ell|v_\ell) \neq 0$ and we may extend each $v_\ell$ to an orthogonal basis.
 We then obtain
an $x$-{\em adapted} orthogonal basis of $T$.
}

{
\begin{lem} \label{lem:Tangent}
Let $x\in X$ as in \eqref{eq:RankOne}.
\begin{enumerate}
\item   Then, relative
to any $x$-adapted basis, $T_x X$ is spanned by all degree-$D$
monomials whose gcd with $x$ has degree at least $D-1$.

\item  Assume moreover that $x$ is non-isotropic. Then, relative
to any $x$-adapted orthogonal basis, $(T_x X)^\perp$ is spanned by all degree-$D$
monomials whose gcd with $x$ has degree at most $D-2$.
\end{enumerate}
\end{lem}

\begin{proof}
Part (1) follows by applying the Leibniz rule to the
parameterisation \eqref{eq:RankOne}
of $X$; part (2) is a straightforward consequence.
\end{proof}
}

\begin{prop} \label{prop:RankOne}
Let $f \in T$ and let $x \in X$ be
non-isotropic. Then the following two statements are equivalent:
\begin{enumerate}
\item some (nonzero) scalar multiple of $x$ is a critical rank-one tensor for
$f$ and
\item a unique (nonzero) scalar multiple of $x$ is a critical rank-one tensor
for $f$;
\end{enumerate}
and they imply the following statement:
\begin{enumerate}
\addtocounter{enumi}{2}
\item for each $\ell \in [p]$, $[f|x]_\ell \in \bigwedge^2 V_\ell$ is
zero.
\end{enumerate}
Moreover, if $f$ is sufficiently general, then {\em every} nonzero $x
\in X$ satisfying (3) is non-isotropic and satisfies (1) and (2).
\end{prop}

The pairing in item (3) is the pairing from \eqref{eq:Pairing}.

\begin{proof}
For the equivalence of the first two statements, we note that if $cx,c'x$
with $c,c' \neq 0$ are critical rank-one tensors for $f$, then $T_{cx}
X = T_{c'x} X=T_x X$ and $f-cx \perp T_x X$ and $f-c'x \perp T_x X$.
Since $x \in T_x X$, we find that $(c-c')x \perp x$, and using that $x$
is non-isotropic we find that $c=c'$.

To prove that (1) implies (3), write $x$ as in \eqref{eq:RankOne} and
extend each $v_\ell$ to an orthogonal basis of $V_\ell$, so as to obtain
an $x$-adapted orthogonal basis of $T$. Now assume that $f-cx \perp T_x
X$. Then by Lemma~\ref{lem:Tangent}, $f-cx$ is a linear combination of
degree-$D$ monomials whose gcds with $x$ have degrees at most $D-2$.
Hence by Lemma~\ref{lem:Pairing}, $[x|f-cx]_\ell=0$. By the skew-symmetry,
$[x|x]_\ell=0$, so $[x|f]_\ell=0$.

For the last statement, consider an $x=v_1^{d_1} \otimes \cdots
\otimes v_p^{d_p} \in X$ where, say, $v_1,\ldots,v_a$ with $a>0$ are isotropic
but the remaining factors are not. Extend each $v_\ell$, $\ell > a$
to an orthogonal basis of $V_\ell$, and for $v_\ell$ with $\ell \leq a$
find an isotropic $w_\ell \in V_i$ with $(v_\ell|w_\ell) =1$ and extend
$v_\ell,w_\ell$ with an orthogonal basis of the orthogonal complement
of $\langle v_\ell, w_\ell \rangle^\perp$ to a basis of $V_\ell$.  In the
corresponding (non-orthogonal) monomial basis of $T$, the monomials $y$
with $[y|x]_\ell \neq 0$ for $\ell \leq a$ are those of the form
\[ w_1^{d_1} \otimes \cdots \otimes w_\ell^{d_\ell-1} u_\ell \otimes
\cdots \otimes w_a^{d_a} \otimes v_{a+1}^{d_{a+1}} \otimes \cdots
\otimes v_p^{d_p} \]
where $u_\ell$ is a basis vector of $V_\ell$ that is distinct from
$v_\ell$ but possibly equal to $w_\ell$. These monomials all satisfy
$[y|x]_i=0$ for $i \neq \ell$. Similarly, the monomials $y$ with
$[y|x]_\ell \neq 0$ for $\ell > a$ are those of the form
\[ w_1^{d_1} \otimes \cdots \otimes w_a^{d_a} \otimes
v_{a+1}^{d_{a+1}} \otimes \cdots \otimes v_l^{d_l-1} u_l \otimes
\cdots \otimes v_p^{d_p} \]
with $u_\ell$ a basis vector of $V_\ell$ distinct from $v_\ell$;
they, too, satisfy $[y|x]_i=0$ for $i \neq \ell$.
The remaining monomials span the space of $f$s with $[x|f]_\ell=0$
for all $\ell$; this space therefore has dimension
\[ \dim T - (n_1+\ldots+n_p), \]
and it does not change when we scale $x$. Since the isotropic projective
points $\langle x \rangle \in \mathbb{P} T$ form a subvariety of positive
codimension in the $(n_1+\ldots+n_p)$-dimensional projective variety
$\mathbb{P} X$, the locus of all $f$ for which there is a nonzero isotropic $x
\in X$ with $[f|x]_\ell=0$ for all $\ell$ has dimension less than
$\dim T$.

Now assume that $f$ is sufficiently general and let $x \in X \setminus
\{0\}$ satisfy $[x|f]_\ell=0$ for all $\ell$. By the above, $x$ is
non-isotropic. Suppose that $f$, expanded on the $x$-adapted orthogonal
basis, contains a monomial $y$ whose gcd with $x$ has degree exactly
$D-1$. If $y$ agrees with $x$ except in the factor $S^{d_\ell} V_\ell$
where it equals $v_\ell^{d_\ell-1} u_\ell$, then in $[x|f]_\ell$,
expanded on the standard basis of $\bigwedge^2 V_\ell$ relative to the
chosen basis of $V_\ell$, the term $v_\ell \wedge u_\ell$ has a nonzero
coefficient. Hence $[x|f]_\ell$ is nonzero, a contradiction.

Therefore, $f$ contains only monomials whose gcds with $x$ have degrees
at most $D-2$, and possibly the monomial $x$ itself. Then $f-c x \perp
T_x X$ for a unique constant $c$. By generality of $f$, it does not lie
in $(T_x X)^\perp$ for any $x \in X \setminus \{0\}$ (the union of these
orthogonal complements is the cone over the variety dual to the projective
variety defined by $X$, and of positive codimension). Hence $c \neq 0$,
and $cx$ is a critical rank-one tensor for $f$.
\end{proof}

{
\begin{remark}
The implication (1) $\Longrightarrow$ (3) in Proposition \ref{prop:RankOne} holds without the assumption of non-isotropy of $x$.
This follows from the fact that the ED correspondence $$\{(x,f)\in X\times V|\ x\textrm{\ is critical for\ }f\}$$
is a irreducible variety (see \cite[\S 4 and Lemma 2.1]{DHOST}) and the nonempty open part in it where $x$ is non-isotropic
lies in the variety defined by $[f|x]_\ell=0\ \forall\ell \in [p]$ by Proposition \ref{prop:RankOne}.
\end{remark}
}

\subsection*{The critical space}

In view of Proposition~\ref{prop:RankOne}, we introduce the following
notion.

\begin{defn}
For a tensor $f \in T$, {\em the critical space $H_f \subseteq T$ of $f$}
is defined as
\[ H_f:=\{g \in T \mid [f|g]_\ell=0 \text{ for all } \ell \in [p]\}.
\]
\end{defn}

\begin{remark}\label{rem:fincritical} By the skew-symmetry, it follows immediately that
$f\in H_f$.
\end{remark}

\begin{remark}
In the case of binary forms ($\dim V=2$), $H_f$ is the hyperplane orthogonal to $D(f)$ \cite{OT}. In the case of ordinary tensors,
$H_f$ was first defined in \cite{OP} where it was called singular space, but
in view of the results in this paper we feel that {\em critical space}
is a better name.
\end{remark}

Proposition~\ref{prop:RankOne} establishes that the non-isotropical
critical rank-one tensors all lie inside $H_f$; hence for a sufficiently
general $f$, all critical rank-one tensors lie in $H_f$.  In the next
subsection we will establish an analogous statement for higher ranks.

Note that the number of linear conditions for $g$ to lie in $H_f$
is at most $\sum_{\ell=1}^p \dim \bigwedge^2 V_\ell = \sum_{\ell=1}^p
\binom{n_\ell+1}{2}$---the linear conditions in the definition may not
all be linearly independent. In Proposition~\ref{prop:dim} we will see
that, assuming the triangle inequalities from Theorem~\ref{thm:Main} and
assuming that $f$ is sufficiently general, equality holds.

\subsection*{Higher rank}

We will now establish a generalization of Proposition~\ref{prop:RankOne}
to higher-rank tensors.

\red{
\begin{defn}
Let $f \in T$ and let $k$ be any nonnegative integer.
A {\em critical rank-at-most-$k$ tensor} for $f$
is a tensor $g \in kX$ such that $f-g \perp T_g \operatorname{Sec}_k(X)$.
\end{defn}
}

\red{Note that by \cite[Lemma 4.2]{DLOT}, all the critical
rank-at-most-$k$ tensors for a sufficiently general $f \in T$ are smooth
points of $\operatorname{Sec}_k(X)$ and can be written as a sum of $k$ non-isotropic
rank-one tensors. Moreover, if we assume that $k$ is at most the generic
rank of tensors in $T$, then these critical tensors to a sufficiently
general $f$ have rank equal to $k$. If $k$ is at least the generic
rank of tensors in $T$, then the only critical
rank-at-most-$k$ tensor for a sufficiently general $f$ is $f$ itself.}

\begin{prop} \label{prop:HigherRank}
Let $f \in T$ be sufficiently general \red{and let $k$ be a nonnegative
integer}. Then all the critical \red{rank-at-most-$k$} tensors for $f$ lie in
the critical space $H_f$.
\end{prop}

\begin{proof}
Let $g$ be a critical \red{rank-at-most-$k$} tensor. By generality of $f$, $g$
can be written as $x_1+\cdots+x_k$ with each $x_i \in X$ non-isotropic.
Then $T_g \operatorname{Sec}_k X \supseteq \sum_{i=1}^k T_{x_i} X$, so that for each $i \in
[k]$ we have $f-g \perp T_{x_i} X$. By
Lemma~\ref{lem:Tangent} this means that, in the $x_i$-adapted
orthogonal basis, $f-g$ is a linear combination of monomials whose gcds
with $x_i$ have degrees at most $D-2$. Hence, by Lemma~\ref{lem:Pairing},
$[f-g|x_i]_\ell=0$ for all $\ell=1,\ldots,p$. We conclude that, for
each $\ell$,
\[ [f-g|g]_\ell=\sum_{i=1}^k [f-g|x_i]_\ell=0, \]
and therefore
\[ [f|g]_\ell=[f-g|g]_\ell + [g|g]_\ell=0 + 0, \]
where in the last step we used that $[.|.]_\ell$ is
skew-symmetric. Hence $g \in H_f$.
\end{proof}

In the next section we compute the dimension of the space spanned by
the critical rank-one tensors for a general tensor, and show that this
space equals $H_f$.

\section{The scheme of critical rank-one tensors}
\label{sec:Span}

\subsection*{Critical rank-one tensors as the zero locus of a vector
bundle section}

Let $f \in T=\bigotimes_{\ell=1}^p S^{d_\ell} V_\ell$ be a tensor. We
assume that $p \geq 2$, $d_\ell \geq 1$, and $\dim V_\ell=n_\ell+1 \geq 1$
for all $\ell$. We adapt the notation of \cite[Section 5.1]{OP} to our
current setting.

Consider the Segre-Veronese variety $\mathbb{P} X=\mathbb{P} V_1\times\ldots\times\mathbb{P}
V_p$ embedded with $\mathcal{O}(d_1,\ldots,d_p)$ in $\mathbb{P} T$; so $\mathbb{P} X$
is the projective variety associated to the affine cone $X \subseteq T$.
Let $\pi_\ell:\mathbb{P} X \rightarrow \mathbb{P} V_\ell$ be the projection on the
$\ell$-th factor and set $N:=\dim \mathbb{P} X = n_1+\ldots+n_p$. For each $\ell
\in [p]$ let $\mathcal{Q}_\ell$ be the quotient bundle on $\mathbb{P} V_\ell$, whose
fibre over a point $\langle v \rangle$ is $V_\ell/\langle v \rangle$. From
these quotient bundles, we construct the following vector bundles on
$\mathbb{P} X$:
\[ \mathcal{E}:=\bigoplus_{\ell=1}^{p} \mathcal{E}_l \quad \text{where}
\quad \mathcal{E}_l:= (\pi_\ell^*\mathcal{Q_\ell})\otimes
\mathcal{O}(d_1,\ldots,d_{\ell-1},d_\ell-1,d_{\ell+1},\ldots,d_p). \]
Note that $\mathcal{E}$ has rank $N$.  The fibre of $\mathcal{E}_\ell$ over a point
$v:=(\langle v_1 \rangle,\ldots,\langle v_p \rangle) \in \mathbb{P} X$
consists of polynomial maps $\prod_{i=1}^p \langle v_i \rangle
\to V_\ell/\langle v_\ell \rangle$ that are multi-homogeneous of
multi-degree $(d_1,\ldots,d_\ell-1,\ldots,d_p)$.  The tensor $f$
yields a global section of $\mathcal{E}_\ell$ which over the point $v$ is
the map sending $(c_1 v_1,\ldots,c_p v_p)$ to the natural pairing of
$f$ with $(c_1 v_1)^{d_1} \cdots (c_\ell v_\ell)^{d_\ell-1} \cdots
(c_p v_p)^{d_p}$---a vector in $V_\ell$---taken modulo $\langle v_\ell
\rangle$. Combining these $p$ sections, $f$ yields a global section $s_f$
of $\mathcal{E}$. By Proposition~\ref{prop:RankOne}, for $f$ sufficiently general,
the tensor $x:=v_1^{d_1} \otimes \cdots
\otimes v_p^{d_p}$ is a nonzero scalar multiple of a critical rank-one
tensor for $f$ if and only if the point $(\langle v_1 \rangle, \ldots,
\langle v_p \rangle)$ is in the zero locus $Z_f$ of the section $s_f$. In
\cite{FriOtt}, this is used to compute the cardinality of $Z_f$ for $f$
sufficiently general as the top Chern class of $\mathcal{E}$. Our current task is
different: we want to show that, if we assume the triangle inequalities of
Theorem~\ref{thm:Main} and that $f$ is sufficiently general, the linear
span $\langle Z_f \rangle$ equals the projectivised critical space $\mathbb{P}
H_f$; this is the second part of Theorem~\ref{thm:Main}.

\subsection*{Bott's formulas and a consequence}

Our central tool will be the following formulas for the cohomology
of vector bundles over projective spaces \cite{OSSBott}. Recall that
$\Omega_{\mathbb{P}^n}^r (k)$ is the $\mathcal{O}(k)$-twisted sheaf of differential
$r$-forms on $\mathbb{P}^n$.

\begin{lem}[Bott's formulas]
For $q,n,r \in \mathbb{Z}_{\geq 0}$ and $k \in \mathbb{Z}$ we have
\[ h^q(\mathbb{P}^n,\Omega_{\mathbb{P}^n}^r(k))
=
\begin{cases}
\binom{k+n-r}{k} \binom{k-1}{r} & \text{if $q=0 \leq r \leq n$ and
$k>r$,} \\
1 & \text{if $0 \leq q=r \leq n$ and $k=0$,}\\
\binom{-k+r}{-k} \binom{-k-1}{n-r} & \text{if $q=n \geq r \geq 0$ and
$k<r-n$, and}\\
0 & \text{otherwise.}
\end{cases}
\]
\end{lem}

A consequence featuring the triangle inequalities of
Theorem~\ref{thm:Main} is the following.

\begin{lem} \label{lem:ConsBott}
Suppose that $n_\ell \leq \sum_{i \neq \ell} n_{i}$ holds for all $\ell$ with
$d_\ell=1$. Let $k \geq 2$ be an integer, $q_1,\ldots,q_p$ be nonnegative
integers with $\sum_{\ell=1}^p q_\ell < k$ and $r_1,\ldots,r_p$ be
nonnegative integers with $\sum_{\ell=1}^p r_\ell =k$. Then
\[ \bigotimes_{\ell=1}^p H^{q_\ell}(\mathbb{P} V_\ell,
\Omega_{\mathbb{P} V_\ell}^{r_\ell}(-d_\ell(k-1)+2r_\ell)) =0. \]
\end{lem}

\begin{proof}
Assume that all factors in the tensor product are nonzero.

First, if all of the factors were nonzero by virtue of the second and
third line in Bott's formulas, then we would have $q_\ell \geq r_\ell$ for
all
$\ell$, and hence $k>\sum_\ell q_\ell \geq \sum_\ell r_\ell=k$, a contradiction.

So some factor is nonzero by virtue of the first line in Bott's
formulas;
without loss of generality this is the first factor. Hence we have
$q_1=0
\leq r_1 \leq n_1$ and $-d_1(k-1)+2r_1>r_1$. This last inequality
reads
$r_1>d_1(k-1)$. Combining this with $\sum_\ell r_\ell=k$ and the fact that
$d_1$
is a positive integer, we find that
$r_1=k$, $d_1=1$, and $r_\ell=0$ for $\ell>1$.
In particular, there are no $\ell>1$ for which the first line in Bott's
formulas applies.

For any $\ell > 1$, if the second line applies, then
$0=r_\ell=q_\ell=-d_\ell(k-1)+2r_\ell$, which contradicts that both
$d_\ell$ and
$k-1$
are positive. Hence the third line applies for all $\ell>1$, and in
particular we have $q_\ell=n_\ell$. But then
\[ n_1 \geq r_1 = k > \sum_{l=1}^p q_l =
\sum_{l=2}^p
n_l,\]
which together with $d_1=1$ contradicts the triangle inequality in the
lemma.
\end{proof}

\subsection*{Vanishing cohomology}

The vanishing result in this subsection uses Lemma~\ref{lem:ConsBott}
and the following version of K\"unneth's formula.

\begin{lem}[K\"unneth's formula]
For vector bundles $\mathcal{G}_\ell$ on $\mathbb{P} V_\ell$ for
$\ell=1,\ldots,p$ and a nonnegative integer $q$ we have
\[ H^q(\mathbb{P} X,\bigotimes_\ell \pi_\ell^* \mathcal{G}_\ell)
\cong \bigoplus_{q_1+\ldots+q_p=q} \bigotimes_\ell H^{q_\ell}(\mathbb{P}
V_\ell,\mathcal{G}_\ell), \]
where the sum is over all $p$-tuples of nonnegative integers summing
to $q$.
\end{lem}

\begin{lem}\label{lema3-SV}
Suppose that $n_\ell \leq \sum_{i\neq \ell} n_i$ holds for all $\ell$ such
that $d_\ell=1$. Let $k \geq 2$ be an integer and $q \in \{0,\ldots,k-1\}$. Then
we have
\[
H^q(\mathbb{P} X,(\bigwedge^{k}\mathcal{E}^*) \otimes \mathcal{O}(d_1,\ldots,d_p))=0.
\]
\end{lem}

\begin{proof}
First,
\[ \mathcal{E}^*=\bigoplus_{\ell=1}^{p}(\pi_\ell^*\mathcal{Q_\ell}^*)\otimes
\mathcal{O}(-d_1,\ldots,-d_{\ell-1},-(d_\ell-1),-d_{\ell+1},\ldots,-d_p).
\]
A well-known formula for $k$-th wedge power of a direct sum yields
\[
\bigwedge^k\mathcal{E}^*=\bigoplus_{r_1+\ldots+r_p=k}
\bigotimes_\ell
\bigwedge^{r_\ell}(\pi_\ell^*\mathcal{Q_\ell}^*\otimes \mathcal{O}(-d_1,
\ldots,-(d_\ell-1),\ldots,-d_p)). \]
Using
$\bigwedge^r(\mathcal{F}\otimes
\mathcal{O}(\omega))=(\bigwedge^r\mathcal{F})(\red{r}\omega)$,
$\mathcal{Q}^*=\Omega^1(1)$, and
$\bigwedge^r (\Omega^1(1))=\Omega^{r}(r)$,
we obtain
\[
\bigwedge^k\mathcal{E}^*=\bigoplus_{r_1+\ldots+r_p=k}
\bigotimes_{\ell} (\pi_\ell^*\Omega_{\mathbb{P} V_\ell}^{r_\ell}(r_\ell)\otimes
\mathcal{O}(-r_\ell d_1,\ldots,-r_\ell(d_\ell-1),\ldots,-r_\ell d_p).
\]
Twisting by $\mathcal{O}(d_1,\ldots,d_p)$, regrouping in
each projection, and using $\sum_\ell r_\ell=k$ we find:
\[
(\bigwedge^k\mathcal{E}^*) \otimes \mathcal{O}(d_1,\ldots,d_p)=\bigoplus_{r_1+\ldots+r_p=k}
\bigotimes_{\ell} \pi_\ell^*\Omega_{\mathbb{P} V_\ell}^{r_\ell}(-d_\ell(k-1)+2 r_\ell).
\]
To compute $H^q$ of each summand we apply K\"unneth's formula, and
obtain subsummands which are exactly of the form in
Lemma~\ref{lem:ConsBott}, hence zero.
\end{proof}

\subsection{Comparing $\mathbb{P} H_f$ and $\langle Z_f \rangle$}

Assume that $f$ is sufficiently general in $T$.  By the first subsection
of this section and by Proposition~\ref{prop:RankOne}, $Z_f$ is contained
in the projectivised critical space $\mathbb{P} H_f$, hence so is $\langle
Z_f \rangle$.  Our goal now is to show that $\langle Z_f \rangle$ is
{\em equal} to $\mathbb{P} H_f$ and to compute its dimension. Both of these
goals are achieved through the following lemma. The section $s_f$ of
$\mathcal{E}$ yields a homomorphism $\mathcal{E}^* \to \mathcal{O}$ of sheaves whose image is
contained in the ideal sheaf $\mathcal{I}_{Z_f}$ of the zero locus of $s_f$.

\begin{lem} \label{lem:linforms}
Assume that for each $\ell \in [p]$ we have $n_\ell \leq \sum_{i \neq
\ell} n_i$ and that $f$ is sufficiently general. Then the induced
homomorphism $\mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p) \to \mathcal{I}_{Z_f} \otimes
\mathcal{O}(d_1,\ldots,d_p)$ induces an isomorphism at the level of global
sections.
\end{lem}

The following proof can be shortened considerably using spectral
sequences, but we found it more informative in its current form. To make
the formulas more transparent, we write $H^q(.)$ instead of $H^q(\mathbb{P}
X,.)$ everywhere.

\begin{proof}
To establish the desired isomorphism
\[ H^0(\mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p))
\cong H^0(\mathcal{I}_{Z_f} \otimes \mathcal{O}(d_1,\ldots,d_p)) \]
we use the following Koszul complex (see, e.g., \cite[Chapter III,
Proposition 7.10A]{Hartshorne}):
\[ 0=\bigwedge^{N+1} \mathcal{E}^* \rightarrow\bigwedge^N\mathcal{E}^*\rightarrow\cdots\rightarrow\bigwedge^2
\mathcal{E}^*\rightarrow \mathcal{E}^*\rightarrow \mathcal{I}_Z\rightarrow 0. \]
Letting $\mathcal{F}_k$ be the quotient of $\bigwedge^k \mathcal{E}^*$ by the image of
$\bigwedge^{k+1} \mathcal{E}^*$, this yields the short exact sequence
\[ 0 \to \mathcal{F}_2 \to \mathcal{E}^* \to \mathcal{I}_Z \to 0. \]
Tensoring with $\mathcal{O}(d_1,\ldots,d_p)$ yields the short exact sequence
\[ 0 \to \mathcal{F}_2 \otimes \mathcal{O}(d_1,\ldots,d_p) \to
\mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p) \to
\mathcal{I}_Z \otimes \mathcal{O}(d_1,\ldots,d_p) \to 0, \]
and this gives a long exact sequence in cohomology beginning with
\begin{align*}
0 &\to H^0(\mathcal{F}_2 \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
H^0 (\mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
H^0 (\mathcal{I}_Z \otimes \mathcal{O}(d_1,\ldots,d_p))\\ & \to
H^1 (\mathcal{F}_2 \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
\end{align*}
So to obtain the desired isomorphism we want that
\[ H^q(\mathcal{F}_2 \otimes \mathcal{O}(d_1,\ldots,d_p))=0 \quad \text{for $q=0,1$.}\]
For each $k=2,\ldots,N$ we have the short exact sequence
\[ 0 \to \mathcal{F}_{k+1} \to \bigwedge^k \mathcal{E}^* \to \mathcal{F}_k \to 0 \]
which yields the long exact sequence
\begin{align*}
& \to
H^{k-2} (\bigwedge^k \mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
H^{k-2} (\mathcal{F}_k \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
H^{k-1} (\mathcal{F}_{k+1} \otimes \mathcal{O}(d_1,\ldots,d_p)) \\ & \to
H^{k-1} (\bigwedge^k \mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
H^{k-1} (\mathcal{F}_k \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
H^k (\mathcal{F}_{k+1} \otimes \mathcal{O}(d_1,\ldots,d_p)) \to
\end{align*}
Using Lemma~\ref{lema3-SV} the two leftmost spaces are zero, so that
\begin{align*}
H^{k-2} (\mathcal{F}_k \otimes \mathcal{O}(d_1,\ldots,d_p))
&\cong H^{k-1} (\mathcal{F}_{k+1} \otimes \mathcal{O}(d_1,\ldots,d_p)) \text{ and}\\
H^{k-1} (\mathcal{F}_k \otimes \mathcal{O}(d_1,\ldots,d_p))
&\subseteq H^{k} (\mathcal{F}_{k+1} \otimes \mathcal{O}(d_1,\ldots,d_p)) \text{ and}\\
\end{align*}
Hence using that $\mathcal{F}_{N+1}=0$ we find that
\begin{align*}
H^0 (\mathcal{F}_2 \otimes \mathcal{O}(d_1,\ldots,d_p)) &\cong \cdots \cong H^{N-1}(\mathcal{F}_{N+1} \otimes
\mathcal{O}(d_1,\ldots,d_p)) = 0 \text{ and}\\
H^1 (\mathcal{F}_2 \otimes \mathcal{O}(d_1,\ldots,d_p)) &\subseteq \cdots \subseteq H^{N}(\mathcal{F}_{N+1} \otimes
\mathcal{O}(d_1,\ldots,d_p)) = 0,
\end{align*}
as desired.
\end{proof}

\begin{prop} \label{prop:dim}
Suppose that for each $\ell \in [p]$ we have $n_\ell \leq \sum_{i \neq
\ell} n_i$ and that $f$ is sufficiently general. Then $\langle Z_f
\rangle=\mathbb{P} H_f$ and $\operatorname{codim}_T H_f=\sum_\ell \binom{n_\ell+1}{2}$.
\end{prop}

\begin{proof}
Since $\mathbb{P} X$ is embedded by $\mathcal{O}(d_1,\ldots,d_p)$, the
space of linear forms on $T$ vanishing on $Z_f$ is $H^0(\mathcal{I}_{Z_f} \otimes
\mathcal{O}(d_1,\ldots,d_p))$. By Lemma~\ref{lem:linforms} this space is
isomorphic to
\begin{align*}
&H^0(\mathcal{E}^* \otimes \mathcal{O}(d_1,\ldots,d_p)) =
\bigoplus_\ell H^0(\pi_\ell^* \mathcal{Q}_l^* \otimes
\mathcal{O}(0,\ldots,1,\ldots,0))\\
&=\bigoplus_\ell H^0(\pi_{\ell^*}(\Omega^1_{\mathbb{P}
V_\ell}(2)))
=\bigoplus_\ell H^0(\mathbb{P} V_\ell, \Omega^1_{\mathbb{P} V_\ell}(2)),
\end{align*}
which by the first line in Bott's formulas has dimension $\sum_{\ell}
\binom{n_\ell+1}{2}$. This means that $\operatorname{codim}_{\mathbb{P} T} \langle Z_f
\rangle=\sum_\ell \binom{n_\ell+1}{2}$, so the second statement in the
proposition follows from the first statement.

To establish the first statement, we spell out the map
\[ H^0(\mathbb{P} V_\ell, \mathcal{Q}_\ell^* \otimes \mathcal{O}(1))=H^0(\mathbb{P} V_\ell, \Omega_{\mathbb{P} V_\ell}^1 (2)) \to H^0(\mathcal{I}_{Z_f} \otimes
\mathcal{O}(d_1,\ldots,d_p)) \]
in greater detail. The space on the left is canonically $(\bigwedge^2
V_\ell)^*$, and an element $\xi$ in this space is mapped to the linear form
$T \to \mathbb{C}, g \mapsto \xi([f|g]_\ell)$. As $\ell$ varies, these
are precisely the linear forms that cut out $H_f$. This
proves that $\mathbb{P} H_f = \langle Z_f \rangle$.
\end{proof}

\begin{remark}
In general, for the equality $\langle Z_f \rangle=\mathbb{P} H_f$ we only need
that the linear equations cutting out $\mathbb{P} H_f$ also cut out $Z_f$,
i.e., we only need that the linear map in Lemma~\ref{lem:linforms}
is surjective.  One might wonder whether this surjectivity remains true
when the triangle inequalities fail. In the case of $(n_1+1) \times
(n_2+1)$-matrices, it does indeed---there we already knew the critical
rank-one approximations span the critical space---but for $p=3$ and
$2 \times 2 \times 4$-tensors (so that $n_3=3>1+1=n_1+n_2$)
the space $\langle Z_f \rangle$ has dimension $6$ while computer
experiments suggest that the space $\mathbb{P} H_f$ has dimension $7$
, hence the surjectivity fails.  Still, in these experiments,
$f$ itself seems to lie in the span of $Z_f$. This leads to the open
problem whether our analogue of the Spectral Theorem and the Eckart-Young
Theorem persists when the triangle inequalities fail.
\end{remark}

\section{Proofs of the main results}
\label{sec:Final}

\begin{proof}[Proof of Theorem~\ref{thm:Main}]
The first statement is Proposition~\ref{prop:HigherRank}; the second
and \red{third} statement are Proposition~\ref{prop:dim}. \red{The last statement follows from Remark \ref{rem:fincritical}.}
\end{proof}

\begin{proof}[Proof of Corollaries~\ref{cor:tensor} and~\ref{cor:sym}.]
If $g$ is a real tensor of real rank at most $k$ closest to $f$, then one
can write it as $x_1+\ldots+x_k$ with $x_1,\ldots,x_k$ real points of $X$.
In particular, all of these points are non-isotropic, and the argument
of Proposition~\ref{prop:HigherRank} applies. Hence $g$ lies in $H_f$.
Now the result follows from Proposition~\ref{prop:dim}. The argument
applies, in particular, to $k$ equal to the rank of $f$, which gives
the last statement of the corollaries.
\end{proof}

Note that, if $f$ is {\em any} real tensor, then any real tensor of real
rank at most $k$ closest to $f$ lies in $H_f$ by the argument above. Only
for the conclusion that it lies in the span of the complex critical
rank-one tensors of $f$ do we use that $f$ is sufficiently general. We
do not know whether this generality is really needed.

Also note that we do not shed new light on the question of when for
sufficiently general $f$ there {\em exists} a closest real tensor of
rank at most $k$. For an update on the complex case see \cite{Qi}.


\section{Introduction}
\vskip12pt

 Analyses of the flare activity and of variations in the UV emission of the sun are important for research on
solar and terrestrial physics. The accelerated particles generated by powerful flares have a fundamental influence on
space weather and have many practical applications relating to the study of physical processes in the heliosphere
and in the earth's atmosphere.

\vskip12pt
\section{Features of activity cycle 24}
\vskip12pt

The current solar activity cycle (No. 24) is the weakest solar cycle over
the last 100 years or more. At present the transition to the following activity minimum between cycles 24 and 25, which is expected in roughly 2018-2019, is taking place.

Solar cycle 24 continues the trend of recent years of a
reduction in the number of sunspots since cycle 21, which reached its maximum in roughly 1980 (see Fig. 1). Many scientists at NASA who are involved in the prediction of solar activity assume that cycle 25 will be roughly the same or weaker than 24. The variations in the solar activity have been modelled numerically [1] using observational data
on the number of sunspots in the period 1750-2050.

The dynamo theory was used to solve the evolution equation
for a number of relative sunspot number (SSN) with the mechanism for magnetic field formation in sunspots taken into account. The
model results [1] correlate well with observations and predict a prolonged period of low activity up to 2050 similar to the Dalton minimum.

Fig. 1 shows the variations in the number of sunspots over 300 years. The current cycle 24 is clearly one of the lowest and is comparable to the activity during the epoch of the Dalton minimum. It is expected that cycles
25 and 26 will be similar to cycle 24.

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=120mm]{FIGURE1.eps}}
 \caption{Yearly average number of sunspots (SSN) from 1700 to the present.}
{\label{Fi:Fig1}}
\end{figure}

\vskip12pt
\subsection{Solar activity in the ultraviolet range and flare activity during cycle 24}
\vskip12pt

Studies in the ultraviolet (UV) are an important part of research on the development of solar flares. The UV part of the electromagnetic spectrum covers the range from 5 to 400 nm. UV photons are absorbed in the upper layers of the  atmosphere; short-wavelength
fluxes of solar radiation cause ionization and dissociation of atmospheric species and lead to formation of the ionosphere. The solar UV emission is formed in the upper chromosphere and the transition region, while X-rays are
formed in the corona; the fluxes in these ranges make up a comparatively small fraction of the overall radiative flux.

As opposed to most stars of the sun's type, the sun is characterized by comparatively low spottedness, moderate flare
activity, and a low level of corona emission flux [2,3].

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE2.eps}}
 \caption{Cyclical variations in the fluxes $F_{10.7}$, $F_{Ly -\alpha}$, and MgII 280 nm
(core/wings) and the flare index (FI) observed during cycles 21-24.}
{\label{Fi:Fig2}}
\end{figure}

Against this background, even very small variations in the fluxes associated with the appearance and disappearance of groups of spots in the active regions or variations in
activity during the solar cycle and large flares can cause significant changes in the UV and X-ray activity indices.

The UV fluxes analysed in this paper vary by tens of percent on time scales ranging from a few minutes to several hours (during solar flares); by a few percent over times ranging from a few days to several months (owing to the movement of active regions over the disk because of the sun's rotation); and by as much as a factor of two over times
ranging from a few years to decades (variations in the solar cycle). This variability produces variations in the earth's ionosphere and upper atmosphere that are similar in amplitude and time. Since the solar UV emission is mainly absorbed in the upper layers of the atmosphere, it is measured using instruments mounted on rockets and satellites.

Because of the related difficulties, for many years different predictions of the UV fluxes were used that employed correlations with the numbers of sunspots and 10.7 cm radio wave emission ($F_{10.7}$). The $F_{10.7}$ flux is now often used
as an objective index for the current level of solar activity; this flux is measured in sfu (solar flux units), with 1 sfu corresponding to a flux of $10^{-22}$ $W/m^2/Hz$. Thus, a comprehensive analysis of observational data on solar UV emission
is needed, since the emission in this region of the sun's spectrum is the most important input parameter for modelling the state of the earth's ionosphere and determines the space weather in the earth's vicinity.

It has been shown [4] that the relationship between the UV solar activity index and the general activity of the sun based on the value of $F_{10.7}$ differed for cycles 21-23, with changes in the magnitudes of the regression coefficients. The correlation coefficient for a linear regression varies within a single cycle, with two maxima during
the rising and decreasing phases and two minima during the cycle minimum and maximum [4].

Fig. 2 shows time series of the ultraviolet activity indices and the flare index compared with a series of $F_{10.7}$ as indicators of the overall activity of the sun's atmosphere.
Fig. 2 shows that the current cycle 24 has
amplitudes of the UV indices and the flare index that are roughly half those for cycles 21-23. This corresponds to the lower level of solar activity as a whole at the present time.

The goals of this paper are:

(1) to analyse the variations in the UV Lyman alpha 121.7 nm and MgII 280 nm lines, as well as of the flare activity
of the sun in cycle 24;

(2) to compare the fluxes at wavelengths ranging from soft X-ray (data from GOES satellites for 0.1-0.8 nm) to
ultraviolet (data on fluxes in five UV lines from the Solar Dynamics Observatory - SDO) in order to track the
propagation of the largest flare in cycle 24 in the solar atmosphere;

(3) to estimate the delay in the onset of proton events relative to the onset of flares in the 0.1-0.8 nm range associated with these events based on data from the GOES 13-15 satellites.

\vskip12pt
\section{UV fluxes and flare activity in cycle 24}
\vskip12pt

\vskip12pt
\subsection{UV radiation flux in the MgII 280 nm line}
\vskip12pt

The index characterizing the ratio of the fluxes in the center
and wings of the atmospheric MgII (280 nm) line has been observed continuously using instruments on satellites since 1978. The core of the doublet MgII 279.56 and 280.27 nm lines is formed in the chromosphere, but the wings are
mainly of photospheric origin. This index is a good indicator of the state of the solar atmosphere.

A close coupling is observed [5,6] between it and the flux in different lines and ranges of UV emission from the sun, which are important
for modelling the state of the earth's atmosphere. The index MgII (core/wings) is measured in relative units. The
observations from Ref. 7 are used.

Fig. 3 shows the MgII 280 nm (c/w) UV radiation flux as a function of the radio emission flux at a wavelength of 10.7 cm, an indicator of the level of solar activity. The solid points show the dependence for cycles
21-23 and the hollow circles, the dependence for cycle 24. Separate quadratic regression fits are shown for cycles
21-23 and for cycle 24. These fits differ significantly: in cycle 24 the maximum fluxes were roughly half the
magnitude of the maxima in cycles 21-23, the spread in the values is larger in cycle 24, and the regression coefficients
differ substantially.

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE3.eps}}
 \caption{MgII c/w  as a function of
the radio emission flux at 10.7 cm}
{\label{Fi:Fig}}
\end{figure}

\vskip12pt
\subsection{The UV flux in the Lyman $ \alpha $ 121.6 nm line}
\vskip12pt

The flux in the Lyman-alpha hydrogen line ($F_{Ly -\alpha}$) from the entire
solar disk is an important indicator that characterizes the state of the chromosphere and the lower part of the transition region [8]. The index $F_{Ly -\alpha}$ is measured in units of $10^{11} photons/cm^2/s$. The observations from Ref. 9 are used.

Fig. 4 shows the dependence of the ultraviolet emission $F_{Ly -\alpha}$
at 121.6 nm on the 10.7 cm radio emission.
The solid points show this dependence for cycles 21-23 and the hollow triangles, the dependence of cycle 24.

Quadratic regression fits are shown here: the regression curves and fits for these two dependences differ more distinctly
than for the case of the MgII 280 nm c/w index.

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE4.eps}}
 \caption{The $F_{Ly -\alpha}$ 121.6 nm UV radiation flux as a function of
the radio emission flux at 10.7 cm}
{\label{Fi:Fig}}
\end{figure}

\vskip12pt
\subsection{The flare index in cycle 24}
\vskip12pt

The distinctive feature of the time variations in the flare index at the end
of the 20th and beginning of the 21st centuries is that by cycle 23 the maximum of FI had fallen by a factor of two compared to cycles 21 and 22. Archived data on FI from 1975 to 2008 are available at
$$http://www.ngdc.noaa.gov/
stp/space-weather/solar-data/solar-indices$$
For cycle 24 we calculated FI from Ref. 10 using data on flare activity that are available at
$$http://www.wdcb.ru/
stp/data/Solar_Flare_Events/Fl_XXIV.pdf$$
(see Fig. 2). It was found that the amplitude of FI during cycle 24 decreased
by a further factor of 2 compared to cycle 23.

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE5.eps}}
 \caption{The flare index FI as a function of the radio emission
flux at 10.7 cm}
{\label{Fi:Fig}}
\end{figure}

Fig. 5 shows the dependence of FI on the radio emission flux at 10.7 cm. The solid points are for cycles 21-23 and the solid squares are for cycle 24. Quadratic fit curves are shown: the differences in the regression curves for the two dependences are not as marked as in the case of the UV indices, but the spread in the values is greater.

\vskip12pt
\section{ Large flares in cycle 24}
\vskip12pt

Fig. 6 shows the distribution of the MI-X7 (740) flares in cycle 24. The points indicate the months for which the number of flares per month was N<6; the asterisks, the months for which 5 < N < 10; and the large hollow circles, the months for which N > 20.

The data for the flares are plotted on a curve showing the monthly-mean numbers of solar spots. It is clear that the largest number of flares took place in the falling branch after the first maximum,
and during the second maximum and the falling branch of cycle 24. The largest number of large flares occurred
in 2014. Before 2011 and since mid-2015, almost no large flares have been observed. As compared to
cycles 21-23, where the most powerful flares were observed during the decay phase, in cycle 24 the greatest number of
powerful flares (>X2.7) occurred in the rising phase and at the maximum with the exception of the two largest flares of cycle 24, which occurred almost at the minimum between 24 and 25 cycle in September 2017.

\vskip12pt
\section{North-south (N-S) asymmetry in the localization of large flares in cycle 24}
\vskip12pt

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE6.eps}}
 \caption{Monthly averaged distribution of large flares
during cycle 24}
{\label{Fi:Fig}}
\end{figure}

The basic features of the N-S asymmetry have been studied and described in recent years. These show up in all the indices of solar activity, but the reason for this asymmetry is not yet understood. Probably the N-S asymmetry in the indices is a fundamental property of solar activity that is indicative of significant differences in processes taking place in the northern and southern hemispheres. In a study of the N-S asymmetry [11] it was found that many features of the cyclical activity of the sun show up more strongly in the asymmetry of the activity indices
than in the values of these indices as such.

An analysis of solar activity in cycles 21-23 [12] showed that the observed large flares in the soft X-ray range (GOES, 0.1-0.8 nm) in different phases of solar cycles 21-23 had an N-S asymmetry. It was also found that over the entire cycle 21, a small excess of flare activity was observed in the northern hemisphere, while during cycles 22 and 23, a southern excess predominated. It is important to note that an N-S asymmetry of solar flares shows
up both in the number of flares and in the intensities of the flares; this should be taken into account in models of the solar dynamo [12].

We have studied the N-S distribution of flares during cycle 24.

Table 1 lists data from the analysis of the N-S asymmetry for large flares of X-ray class > M5 (119) in different phases of cycle 24. It is clear that a strong asymmetry in the flare activity was observed during 2011 (15 flares in the northern hemisphere and 4 in the southern) and 2014 (11 in the northern hemisphere versus 28 in the southern). During the rest of cycle 24, the number of flares in the two hemispheres were roughly equal.

Fig. 7 illustrates the N-S distribution of large flares (monthly average data for X-ray class > M1) in different phases of cycle 24.
On the plot of the sliding monthly average numbers of sunspots, different symbols are used to denote the months with distinct N-S asymmetries. In this figure, the solid circles indicate the years when the number of flares in the northern hemisphere is substantially higher than in the southern hemisphere, the asterisks indicate predominant flares in the southern hemisphere, and the hollow circles indicate roughly equal number of flares in the two hemispheres.

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE7.eps}}
 \caption{The N-S asymmetry in the number of large
(>M1) flares during cycle 24}
{\label{Fi:Fig}}
\end{figure}

It can be seen from Fig. 7 that in the rising phase of cycle 24 the numbers of large flares are initially roughly equal in the northern and southern hemispheres. Then in 2011 and 2012 the number of flares is seen to be larger in the northern hemisphere; this corresponds to an asymmetry in the number of spots.

During the time of the first maximum, a roughly equal number of flares is observed in the two hemispheres. During the second maximum, in the decreasing phase of the cycle from 2013 to 2016, we see that flares predominate in the southern hemisphere;
this corresponds to the asymmetry of sunspots during cycles 21-23. In 2016, we again see more flares in the northern hemisphere.

An analysis of the rotation velocities of sunspots and the evolution of the active longitudes [13] shows that the northern and southern hemispheres rotate with different velocities and that the rotation velocities in both hemispheres have increased since 1990. It is assumed that the observed north-south asymmetries of groups of sunspots and flare activity are consequences of the temporal evolution of the large-scale magnetic field and the solar dynamo.

It is clear that the asymmetry in the flare activity correlates well with the asymmetry in the appearance of groups
of sunspots according to data in the archives of the Solar Influences Data Analysis Center (SIDC) of the Royal Observatory of Belgium, which includes the World Data Center for the study of sunspots.

According to these data, for cycles 21-23 and the rising phase of cycle 24, there is an N-S asymmetry in the number of spots that differs slightly from the N-S asymmetry in the flare activity; for all of cycles 21-24 the number of spots in the northern hemisphere
predominates in the rising phases of the cycles, while the number of spots in the southern hemisphere predominates in the decreasing phases; see
$$http://www.sidc.be/silso/datafiles$$
The observations of N-S asymmetry in the various
indices are of great significance, since they provide valuable information for modelling the solar dynamo. The appearance of a substantial southern asymmetry at the start of the 21st century indicates the onset of a period of small
activity cycles [14], which we are observing at the present time.

\vskip12pt
\section {Analysis of satellite observations of the X6.9 flare of August 9, 2011}
\vskip12pt

\begin{center}
\begin{table}
\caption{The Asymmetry in Flare Activity During Cycle 24}
\vskip12pt
\begin{tabular}{||c|c|c|c|c||}

\hline \hline

  Year &Number flares&Number flares &Number flares&Number flares\\

     &     > M5   &   > X1     &   in N hemisph. &  in S hemisph. \\  \hline

2010 &  3   &   0 & 2 &  1 \\ \hline

2011 &   19   &  8   &  15 & 4  \\ \hline

 2012 &    21  &  6   & 12  &  9 \\ \hline

2013   &   21   &  12   &  10 & 11 \\ \hline

 2014 &  39    &  16   & 11  &  28 \\ \hline

2015 &  12    &   2  &   5  &  7  \\ \hline

2016 &  4    &  -   &  4   & - \\ \hline

\hline \hline
\end{tabular}
\end{table}
\end{center}

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE8.eps}}
\caption{The evolution of the flare of August 9, 2011, in
different UV lines and in the X-ray range (the 30.4, 33.5,
28.4, and 9.4 nm lines and the 0.1-0.8 nm range)}
{\label{Fi:Fig}}
\end{figure}

UV fluxes have been measured with high accuracy by satellites since 1999. This has made it possible for us to analyse the onset times in different lines of the largest flare during cycle 24, X 6.9, which took place on August
9, 2011 using observational data from SDO/EVE and GOES-15.

This is important for refining the models of flares.
We note the times associated with flare intensification of the fluxes in coronal lines and in lines formed in the chromosphere and transition region. The regions of the coronal flare plasma in which UV and X-ray emission is formed have different heights and temperatures. The gradual increase in the fluxes in these lines during the development of a flare can be used to determine the temperatures within the regions where the flare begins, as well as the in the regions where it propagates.

Fig. 8 shows the evolution of the flare of August 9, 2011, in the coronal lines in the ranges of 9.4, 33.5, and 0.1-0.8 nm, as well as in lines (30.4 and 28.4 nm) of the chromosphere and transition region. The wavelengths are indicated on the right hand side of the figure.

The horizontal smooth lines indicate the level of the background radiation flux; this makes it possible to estimate more accurately the time when the fluxes begin to increase. Note
that for shorter wavelengths of a line or spectral range, the ratio of the maximum flare amplitude to the background level is higher. It can be seen that the rise in the fluxes begins gradually (pre-flare enhancement).

In the 0.1-0.8 nm range and the 9.4 nm line, the pre-flare begins at 7h 48m and 7h 54m, while the pre-flare (which is less distinct) in the other lines begins at 8h 00m.

Thus, the pre-flare enhancement propagates from the higher levels of the corona into
the lower corona and chromosphere. We see the onset of the flare as a sharp, almost vertical increase in the flux
for 0.1-0.8 nm and a steep rise for 9.4 nm at 8h 02m, while the onset of the flare for the other lines occurs a few minutes
later at 8h 04m-8h 05m.

\vskip12pt
\section{Characteristics of the proton events in cycle 24}
\vskip12pt

Solar proton events (SPE) and the flares accompanying them during cycle 24 were analysed on the basis of
data from GOES 13-15. The most powerful SPE with proton energies E>10 MeV and E>100 MeV were identified
and studied.

SPE are observed in interplanetary space as a result of flares on the sun, as a result of ejection of coronal
mass, and sometimes after the disappearance of filaments.

The duration of a proton event depends on the energy
of the protons. For energies ~10 MeV the duration is a few hours; for higher energies they can last up to several days.

At present there are no unique answers to the question of where, when, and how the protons are accelerated and how they propagate to the observer. The power of proton events is measured in

pfu ($particle flux units, protons/
cm^2/s/sr$).

\begin{center}
\begin{table}
\caption{Flares in Cycle 24 Accompanied by the Largest-scale Proton Events}
\vskip12pt
\begin{tabular}{||c|c|c|c|c||}

\hline \hline

  Date &X-ray/optic scale&Position &Proton flux
& Proton flux  \\

     &       &   on the disc     &  (>10 MeV), pfu&  (>100 MeV), pfu\\  \hline

04.08.11 & X2.3/2B  & N19W76 & 99 &  2\\ \hline

23.01.12&  M8.3/2B   &  N28W21   & 4500 & 3 \\ \hline

 27.01.12 & X1.7/2F   &  N27W71  & 800  &  11 \\ \hline

07.03.12  &  X5.4/2B   & N17E27   & 200 & 70 \\ \hline

07.03.12 &  X1.3/SF   &  N23E12   & 35 &  3 \\ \hline

13.03.12 &  M7.9/1B    &   N19W59 &   500  &  2  \\ \hline

17.05.12 &  M5.1/1F   &  N11W76  &  250  & 20 \\ \hline

12.07.12 &  X1.4/2B   &  N13W15  & 100 & 0.2\\ \hline

17.07.12 & M1.7/1F    &   S28W90  &  130  &  0.1  \\ \hline

19.07.12 &  M7.7/SF   &  S16W90   &  90  & 0.7 \\ \hline

11.04.13 &  M6.5/3B    &   N09E12  &   100  &  2  \\ \hline

22.05.13 & M5.0/3N   & N15W70  &  2000   & 4 \\ \hline

07.01.14 &  X1.2/2N    &   S15W11  &   1000 &  4 \\ \hline

10.09.14 &  X1.6/2B   &  N12E02  &  130   & 1\\ \hline

21.06.15 &  M3.8/2B    &   S18W57  &   1000  &  0.1  \\ \hline

06.09.17 &  X9.3/3B    &  S15W45   &  50   & 0.7 \\ \hline

10.09.17 &  X8.2/3B    &  S12W85   &  650   & 70 \\ \hline

\hline \hline
\end{tabular}
\end{table}
\end{center}
\bigskip

\begin{figure}[tbh!]
\centerline{
\includegraphics[width=100mm]{FIGURE9.eps}}
\caption{Histogram of the time delays of the start of proton events relative to the onset
of the flares that produce the events in the soft X-ray range}
{\label{Fi:Fig}}
\end{figure}
\bigskip

Proton events produced by flares and not by other causes (ejection of coronal mass, eruptive prominences) were selected on the basis of the observational data. Proton events in cycles 24 (black) and 23 (grey) were examined separately. The distributions of the delays between the flare onset in the X-ray range and the start of a proton event are similar in cycles 24 and 23.

The maxima of the distributions occur for a delay of 2 hours, both for >10 MeV
protons and >100 MeV protons. Note that the most powerful flares are accompanied by proton events with the shortest delays. The largest delays are observed for proton flares lying in the eastern hemisphere because of a longer trajectory for propagation of the protons from the flare region.

Cycle 24 was distinguished by the fact that no proton fluxes with energies >100 MeV and characterized by a maximum amplitude of more than 100 pfu were observed.
During activity cycle 24, the number of flares accompanied by proton fluxes was lower by roughly a factor of 2.

The characteristics of the energetic solar protons in the first 64 months of cycle 24 (the rising phase and maximum) have been compared with those of the previous cycles (21-23) [15]. It was found that, despite the low solar activity during the rising phase and maximum of cycle 24, the number of SPE with proton energies E>10 MeV and E>100 Mev in the current cycle differs little from the number of these events in cycles 21-23. During the period from 2010-2016, we identified 62 SPE with proton energies E>10 MeV and 24 with E>100 MeV (the most powerful
of these are listed in Table 2).

Some of the proton events (including the event of March 7, 2012) have a complicated
time profile with several maxima and have a high power (see Table 2).
It can be seen that a large fraction of these events originated in the northern hemisphere during 2011-2013, while two out of three of the large flares with proton fluxes during 2014-2015 occurred in the southern
hemisphere.

It is clear from Table 2 that the largest proton events were associated with flares of class M5 or higher. For the large flares in classes M1-X7 accompanied by fluxes of protons with energies >10 MeV and >100 MeV, we have analysed the time delay between the onset of a flare in the X-ray range and the onset of the proton event caused by this flare (see Fig. 9). An estimate of the delay in the arrival of the protons in the earth's atmosphere is
important for predicting space weather and the state of the upper atmosphere.

\vskip12pt
\section{Conclusion}
\vskip12pt

It follows from this paper that the flare activity during cycle 24 was considerably lower than in cycles 21-23. As compared to cycles 21-23, where the most powerful flares were observed during the decay phase, in cycle 24 the greatest number of
powerful flares (>X2.7) occurred in the rising phase and at the maximum with the exception of the two largest flares of cycle 24, which occurred almost at the minimum between 24 and 25 cycles in September 2017.

An analysis of the ultraviolet fluxes in the MgII 280 nm and La 121.6 nm lines showed that the regression fits as a function of the overall emission level of the sun differ substantially in cycle 24 and in cycles 21-23.

The evolution of the flare of August 9, 2011, as observed in coronal lines (9.4, 33.5, and 0.1-0.8 nm) and in lines of the chromosphere and transition region (30.4 and 28.4 nm) proceeds in a direction from the upper corona toward the lower corona and chromosphere.

A study of the N-S asymmetry in the flare distribution during cycle 24 showed that a strong N-asymmetry of the flares was observed in 2011 and an S-asymmetry in 2014. In the other years of cycle 24, the numbers of flares in the two hemispheres were roughly the same.

During cycle 24, proton events with E>100 MeV are less intense; this confirms the conclusions of Ref. 15. The largest-scale proton events were associated with X-ray flares of class >M5.

Histograms of the time delays of the onset of proton events relative to the onset in the soft X-ray range of the flares causing these events during cycles 24 and 23 show that the maximum of the distribution occurs for a delay of 2 hours for protons with energies >10 MeV, as well as for those with energies >100 MeV.