\section{Introduction}

The $F$-pure threshold was introduced by Musta\c t\u a, Takagi, and Watanabe~\cite{TW, MTW}; it is a positive characteristic invariant, analogous to log canonical thresholds in characteristic zero. We calculate the possible values of the $F$-pure threshold of the affine cone over a Calabi-Yau hypersurface, and relate the threshold to the order of vanishing of the Hasse invariant, and to a numerical invariant introduced by van~der~Geer and Katsura in~\cite{GK}.

\begin{theorem}
\label{thm:intro}
Suppose $R=K[x_0,\dots,x_n]$ is a polynomial ring over a field $K$ of characteristic $p>n+1$, and $f$ is a homogeneous polynomial in $R$ defining a smooth Calabi-Yau hypersurface $X=\operatorname{Proj} R/fR$. Then the $F$-pure threshold of $f$ has the form
\[
\mathrm{fpt}(f)\ =\ 1-\frac{h}{p}\,,
\]
where $h$ is an integer with $0\leqslant h\leqslant\dim X$. If $p\geqslant n^2-n-1$, then the integer $h$ equals the order of vanishing of the Hasse invariant on the versal deformation space of $X \subset \mathbf{P}^n$.
\end{theorem}

Hern\'andez has computed $F$-pure thresholds for binomial hypersurfaces~\cite{Hernandez1} and for diagonal hypersurfaces~\cite{Hernandez2}. The $F$-pure threshold is computed for a number of examples in~\cite[Section~4]{MTW}. Example~4.6 of that paper computes the $F$-pure threshold in the case of an ordinary elliptic curve, and raises the question for supersingular elliptic curves; this is answered by the above theorem.

The theory of $F$-pure thresholds is motivated by connections to log canonical thresholds; for simplicity, let $f$ be a homogeneous polynomial with rational coefficients. Using~$f_p$ for the corresponding prime characteristic model, one has
\[
\mathrm{fpt}(f_p)\leqslant\mathrm{lct}(f)\qquad\text{ for all }p\gg0\,,
\]
where $\mathrm{lct}(f)$ denotes the log canonical threshold of $f$, and
\[
\lim_{p\longrightarrow\infty}\mathrm{fpt}(f_p)=\mathrm{lct}(f)\,,
\]
see~\cite[Theorems~3.3,~3.4]{MTW}; this builds on the work of a number of authors, primarily Hara and Yoshida~\cite{Hara-Yoshida}. It is conjectured that $\mathrm{fpt}(f_p)$ and $\mathrm{lct}(f)$ are equal for infinitely many primes; see~\cite{MS} for more in this direction.

The $F$-pure threshold is known to be rational in a number of cases, including for principal ideals in an excellent regular local ring of prime characteristic~\cite{KLZ:JALG}. Other results on rationality include~\cite{BMS:MMJ, BMS:TAMS, BSTZ, Hara, ST}. For more on $F$-pure thresholds, we mention~\cite{BHMM, HMTW, MSV, Mustata, MY}.

\section{The $F$-pure threshold}

In~\cite{TW} the $F$-pure threshold is defined for a pair $(R,\mathfrak{a})$, where $\mathfrak{a}$ is an ideal in an $F$-pure ring of prime characteristic. The following special case is adequate for this paper:

\begin{definition}
Let $(R,\mathfrak{m})$ be a regular local ring of characteristic $p>0$. For an element $f$ in $\mathfrak{m}$, and integer $q=p^e$, we define
\[
\mu_f(q)\ :=\ \min\big\{k\in\mathbf{N}\mid f^k\in\mathfrak{m}^{[q]}\big\}\,,
\]
where $\mathfrak{m}^{[q]}$ denotes the ideal generated by the elements $r^q$ for $r\in\mathfrak{m}$. Note that $\mu_f(1)=1$, and that $1\leqslant\mu_f(q)\leqslant q$. Moreover, $f^{\mu_f(q)}\in\mathfrak{m}^{[q]}$ implies that $f^{p\mu_f(q)}\in\mathfrak{m}^{[pq]}$, and it follows that $\mu_f(pq)\leqslant p\mu_f(q)$. Thus,
\[
\Big\{\frac{\mu_f(p^e)}{p^e}\Big\}_{e\ge0}
\]
is a non-increasing sequence of positive rational numbers; its limit is the \emph{$F$-pure threshold} of $f$, denoted $\mathrm{fpt}(f)$.
\end{definition}

By definition, $f^{\mu_f(q)-1}\notin\mathfrak{m}^{[q]}$. Taking $p$-th powers, and using that $R$ is $F$-pure, we get
\[
f^{p\mu_f(q)-p}\ \notin\ \mathfrak{m}^{[pq]}\,.
\]
Combining with the observation above, one has
\begin{equation}
\label{eqn:mu}
p\mu_f(q)-p+1\ \leqslant\ \mu_f(pq)\ \leqslant\ p\mu_f(q)\,.
\end{equation}
Note that this implies
\[
\mu_f(q)\ =\ \left\lceil{\frac{\mu_f(pq)}{p}}\right\rceil\qquad\text{ for each }q=p^e\,.
\]

The definition is readily adapted to the graded case where $R$ is a polynomial ring with homogeneous maximal ideal $\mathfrak{m}$, and $f$ is a homogeneous polynomial.

\begin{remark}
\label{rmk:fptthickenings}
The numbers $\mu_f(p^e)$ may be interpreted in terms of thickenings of the hypersurface $f$ as follows. Let $K$ be a field of characteristic $p>0$, and $f$ a homogeneous polynomial of degree $d$ in $R=K[x_0,\dots,x_n]$. Fix integers $q=p^e$ and $t\leqslant q$. The Frobenius iterate $F^e\colon R/fR\longrightarrow R/fR$ lifts to a map $R/fR\longrightarrow R/f^qR$; composing this with the canonical surjection $R/f^qR\longrightarrow R/f^tR$, we obtain a map
\[
\widetilde{F^e_t}\colon R/fR\longrightarrow R/f^tR\,.
\]

Consider the commutative diagram with exact rows
\[
\CD
0@>>>R(-d)@>f>>R@>>>R/fR@>>>0\phantom{\,,}\\
@. @VVf^{q-t}F^eV @VVF^eV @VV\widetilde{F^e_t}V \\
0@>>>R(-dt)@>f^t>>R@>>>R/f^tR@>>>0\,,
\endCD
\]
and the induced diagram of local cohomology modules
\begin{equation}
\label{equation:cd}
\CD
0@>>>H^n_\mathfrak{m}(R/fR)@>>>H^{n+1}_\mathfrak{m}(R)(-d)@>f>>H^{n+1}_\mathfrak{m}(R)@>>>0\phantom{\,.}\\
@. @VV\widetilde{F^e_t}V @VVf^{q-t}F^eV @VVF^eV \\
0@>>>H^n_\mathfrak{m}(R/f^tR)@>>>H^{n+1}_\mathfrak{m}(R)(-dt)@>f^t>>H^{n+1}_\mathfrak{m}(R)@>>>0\,.
\endCD
\end{equation}
Since the vertical map on the right is injective, it follows that $\widetilde{F^e_t}$ is injective if and only if the middle map is injective, i.e., if and only if the element
\[
f^{q-t}F^e\left(\left[\frac{1}{x_0\cdots x_n}\right]\right)\ =\ \left[\frac{f^{q-t}}{x_0^q\cdots x_n^q}\right]
\]
is nonzero, equivalently, $f^{q-t}\notin\mathfrak{m}^{[q]}$. Hence $\widetilde{F^e_t}\colon H^n_\mathfrak{m}(R/fR)\longrightarrow H^n_\mathfrak{m}(R/f^tR)$ is injective if and only if $\mu_f(q)>q-t$.
\end{remark}

We note that the generating function of the sequence $\{\mu_f(p^e)\}_{e\ge1}$ is a rational function:

\begin{theorem}
Let $(R,\mathfrak{m})$ be a regular local ring of characteristic $p>0$, and let $f$ be an element of $\mathfrak{m}$. Then the generating function
\[
G_f(z)\ :=\ \sum_{e\geqslant 0}\mu_f(p^e)z^e
\]
is a rational function of $z$ with a simple pole at $z=1/p$; the $F$-pure threshold of $f$ is
\[
\mathrm{fpt}(f)\ =\ \lim_{z\longrightarrow 1/p}(1-pz)\,G_f(z)\,.
\]
\end{theorem}

\begin{proof}
Since the numbers $\mu_f(p^e)$ are unchanged when $R$ is replaced by its $\mathfrak{m}$-adic completion, there is no loss of generality in assuming that $R$ is a complete regular local ring; the rationality of $\mathrm{fpt}(f)$ now follows from~\cite[Theorems~3.1,~4.1]{KLZ:JALG}. Let $\mathrm{fpt}(f)=a/b$ for integers $a$ and $b$. By~\cite[Proposition~1.9]{MTW}, one has
\[
\mu_f(p^e)\ =\ \lceil p^e\,\mathrm{fpt}(f)\rceil\ =\ \left\lceil{\frac{ap^e}{b}}\right\rceil
\qquad\text{for each }q=p^e\,.
\]
Suppose $ap^{e_0}\equiv ap^{e_0+e_1}\mod b$ for integers $e_0$ and $e_1$. Then $ap^{e_0}\equiv ap^{e_0+ke_1}\mod b$ for each integer~$k\ge0$. Hence there exists an integer $c$ such that
\begin{align*}
H(z)\ := \ \sum_{k\ge0}\mu_f(p^{e_0+ke_1})z^{e_0+ke_1}\ &=\ \sum_{k\ge0}\left\lceil{\frac{ap^{e_0+ke_1}}{b}}\right\rceil z^{e_0+ke_1}\\
&=\ \sum_{k\ge0}\frac{ap^{e_0+ke_1}+c}{b}\ z^{e_0+ke_1}\,,
\end{align*}
is a rational function of $z$ with a simple pole at $z=1/p$. Moreover,
\[
\lim_{z\longrightarrow 1/p}(1-pz)\,H(z)\ =\ \frac{a}{be_1}\,.
\]
Partitioning the integers $e\geqslant e_0$ into the congruence classes module $e_1$, it follows that $G_f(z)$ is the sum of a polynomial in $z$ and $e_1$ rational functions of the form
\[
\sum_{k\geqslant 0}\mu_f(p^{\ell+ke})z^{\ell+ke}\,.
\]
The assertions regarding the pole and the limit now follow.
\end{proof}

The theorem holds as well in the graded setting.

\section{Preliminary results}

We record some elementary calculations that will be used later. Here, and in the following sections, $R$ will denote a polynomial ring $K[x_0,\dots,x_n]$ over a field $K$ of characteristic~$p>0$, and $\mathfrak{m}$ will denote its homogeneous maximal ideal. By the \emph{Jacobian ideal} of a polynomial~$f$, we mean the ideal generated by the partial derivatives
\[
f_{x_i}\ :=\ \partial\!f/\partial x_i\quad\text{ for }\ 0\leqslant i\leqslant n\,.
\]
If $f$ is homogeneous of degree coprime to $p$, then the Euler identity ensures that $f$ is an element of the Jacobian ideal; this is then the defining ideal of the singular locus of $R/fR$.

\begin{lemma}
\label{lemma:jacobian}
Let $f$ be a homogeneous polynomial of degree $d$ in $K[x_0,\dots,x_n]$ such that the Jacobian ideal $J$ of $f$ is $\mathfrak{m}$-primary. Then
\[
\mathfrak{m}^{(n+1)(d-2)+1}\ \subseteq\ J\,.
\]
\end{lemma}

\begin{proof}
Since $J$ is $\mathfrak{m}$-primary, it is a complete intersection ideal. As it is generated by forms of degree $d-1$, the Hilbert-Poincar\'e series of $R/J$ is
\[
P(R/J,t)\ =\ \frac{(1-t^{d-1})^{n+1}}{(1-t)^{n+1}}\ =\ (1+t+t^2+\dots+t^{d-2})^{n+1}\,.
\]
It follows that $R/J$ has no nonzero elements of degree greater than $(n+1)(d-2)$.
\end{proof}

\begin{lemma}
\label{lemma:colon}
Let $R=K[x_0,\dots,x_n]$ and $\mathfrak{m}^{[q]}=(x_0^q,\dots,x_n^q)$. Then
\[
\mathfrak{m}^{[q]}:_R\mathfrak{m}^{(n+1)(d-2)+1}\ \subseteq\ \mathfrak{m}^{[q]}+\mathfrak{m}^{(n+1)(q-d+1)}\,,
\]
where $\mathfrak{m}^i=R$ for $i\le0$.
\end{lemma}

\begin{proof}
We prove, more generally, that
\[
\mathfrak{m}^{[q]}:_R\mathfrak{m}^k\ =\
\begin{cases}
\mathfrak{m}^{[q]}+\mathfrak{m}^{nq+q-n-k}&\text{ if }\ 0\leqslant k\leqslant nq+q-n\,,\\
R&\text{ if }\ k\geqslant nq+q-n\,.
\end{cases}
\]
Suppose $r$ is a homogeneous element of $\mathfrak{m}^{[q]}:_R\mathfrak{m}^k$. Computing the local cohomology module~$H^{n+1}_\mathfrak{m}(R)$ via a \v Cech complex on $x_0,\dots,x_n$, the element
\[
\left[\frac{r}{x_0^q\cdots x_n^q}\right]\ \in\ H^{n+1}_\mathfrak{m}(R)
\]
is annihilated by $\mathfrak{m}^k$, and hence lies in ${\left[H^{n+1}_\mathfrak{m}(R)\right]}_{\geqslant -n-k}$. If $r\notin\mathfrak{m}^{[q]}$, then
\[
\deg r-(n+1)q\ \geqslant\ -n-k\,,
\]
i.e., $r\in\mathfrak{m}^{nq+q-n-k}$. The pigeonhole principle implies that $\mathfrak{m}^{nq+q-n}$ is contained in $\mathfrak{m}^{[q]}$, which gives the rest.
\end{proof}

\begin{lemma}
\label{lemma:bound}
Let $f$ be a homogeneous polynomial of degree $d$ in $K[x_0,\dots,x_n]$, such that the Jacobian ideal of $f$ is $\mathfrak{m}$-primary. If $\mu_f(q)$ is not a multiple of $p$, then
\[
\mu_f(q)\ \geqslant\ \frac{(n+1)(q+1)-nd}{d}\,.
\]
\end{lemma}

\begin{proof}
Set $k:=\mu_f(q)$, i.e., $k$ is the least integer such that
\[
f^k\ \in\ \mathfrak{m}^{[q]}\,.
\]
Applying the differential operators $\partial/\partial x_i$ to the above, we see that
\[
kf^{k-1}f_{x_i}\ \in\ \mathfrak{m}^{[q]}\qquad\text{ for each }i\,,
\]
since $\partial/\partial x_i$ maps elements of $\mathfrak{m}^{[q]}$ to elements of $\mathfrak{m}^{[q]}$. As $k$ is nonzero in $K$, one has
\[
f^{k-1}J\ \subseteq\ \mathfrak{m}^{[q]}\,,
\]
where $J$ is the Jacobian ideal of $f$. Lemma~\ref{lemma:jacobian} now implies that
\[
f^{k-1}\mathfrak{m}^{(n+1)(d-2)+1}\ \subseteq\ \mathfrak{m}^{[q]}\,.
\]
By Lemma~\ref{lemma:colon}, we then have
\[
f^{k-1}\ \in\ \mathfrak{m}^{[q]}+\mathfrak{m}^{(n+1)(q-d+1)}\,.
\]
But $f^{k-1}\notin\mathfrak{m}^{[q]}$ by the minimality of $k$, so $f^{k-1}$ has degree at least $(n+1)(q-d+1)$, i.e.,
\[
d(k-1)\ \geqslant\ (n+1)(q-d+1)\,;
\]
rearranging the terms, one obtains the desired inequality
\[
k\ \geqslant\ \frac{(n+1)(q+1)-nd}{d}\,.
\qedhere
\]
\end{proof}

\begin{lemma}
\label{lemma:patterns}
Let $f$ be a homogeneous polynomial of degree $d$ in $K[x_0,\dots,x_n]$, such that the Jacobian ideal of $f$ is $\mathfrak{m}$-primary.
\begin{enumerate}[\quad\rm(1)]
\item If \ $\displaystyle{\frac{\mu_f(q)-1}{q-1}=\frac{n+1}{d}}$ \ for some $q=p^e$, then \ $\displaystyle{\frac{\mu_f(pq)-1}{pq-1}=\frac{n+1}{d}}$.

\item Suppose $p\geqslant nd-d-n$. If \ $\displaystyle{\frac{\mu_f(q)}{q}<\frac{n+1}{d}}$ \ for some $q=p^e$, then $\mu_f(pq)=p\mu_f(q)$.
\end{enumerate}
\end{lemma}

\begin{proof}
(1) Since $f^{\mu_f(q)-1}$ has degree $(q-1)(n+1)$ and is not an element of $\mathfrak{m}^{[q]}$, it must generate the socle in $R/\mathfrak{m}^{[q]}$. But then
\[
\left(f^{\mu_f(q)-1}\right)^{\frac{pq-1}{q-1}}
\]
generates the socle in $R/\mathfrak{m}^{[pq]}$, so
\[
\mu_f(pq)-1\ =\ \big(\mu_f(q)-1\big)\left(\frac{pq-1}{q-1}\right)\,.
\]

For (2), suppose that $\mu_f(pq)<p\mu_f(q)$. Then $\mu_f(pq)$ is not a multiple of $p$ by~\eqref{eqn:mu}. Lemma~\ref{lemma:bound} thus implies that
\[
(n+1)(pq+1)-nd\ \leqslant\ d\mu_f(pq)\,.
\]
Combining with $\mu_f(pq)\leqslant p\mu_f(q)-1$ and $d\mu_f(q)\leqslant q(n+1)-1$, we obtain
\[
p\ \leqslant\ nd-d-n-1\,,
\]
which contradicts the assumption on $p$.
\end{proof}

We next prove a result on the injectivity of the Frobenius action on negatively graded components of local cohomology modules:

\begin{theorem}
\label{thm:highdeg}
Let $K$ be a field of characteristic $p>0$. Let $f$ be a homogeneous polynomial of degree $d$ in $R=K[x_0,\dots,x_n]$, such that the Jacobian ideal of $f$ is primary to the homogeneous maximal ideal $\mathfrak{m}$ of $R$.

If $p\geqslant nd-d-n$, then the Frobenius action below is injective:
\[
F\colon{\left[H^n_\mathfrak{m}(R/fR)\right]}_{<0}\longrightarrow{\left[H^n_\mathfrak{m}(R/fR)\right]}_{<0}\,.
\]
\end{theorem}

\begin{proof}
Using~\eqref{equation:cd} in the case $t=1$ and $e=1$, and restricting to the relevant graded components, we have the diagram with exact rows
\[
\CD
0@>>>{\left[H^n_\mathfrak{m}(R/fR)\right]}_{\leqslant-1}@>>>{\left[H^{n+1}_\mathfrak{m}(R)\right]}_{\leqslant-d-1}@>>>\cdots\phantom{\,.}\\
@. @VVFV @VVf^{p-1}FV\\
0@>>>{\left[H^n_\mathfrak{m}(R/fR)\right]}_{\leqslant-p}@>>>{\left[H^{n+1}_\mathfrak{m}(R)\right]}_{\leqslant-d-p}@>>>\cdots\,.
\endCD
\]
Thus, it suffices to prove the injectivity of the map
\[
f^{p-1}F\colon{\left[H^{n+1}_\mathfrak{m}(R)\right]}_{\leqslant-d-1}\longrightarrow{\left[H^{n+1}_\mathfrak{m}(R)\right]}_{\leqslant-d-p}\,.
\]
A homogeneous element of ${[H^{n+1}_\mathfrak{m}(R)]}_{\leqslant-d-1}$ may be written as
\[
\left[\frac{g}{(x_0\cdots x_n)^{q/p}}\right]
\]
for some $q$, where $g\in R$ is homogeneous of degree at most $(n+1)q/p-d-1$. Suppose
\[
f^{p-1}F\left( \left[\frac{g}{(x_0\cdots x_n)^{q/p}}\right] \right)\ =\ 0\,,
\]
then it follows that $f^{p-1}g^p\in\mathfrak{m}^{[q]}$. Let $k$ be the least integer with
\[
f^kg^p\ \in\ \mathfrak{m}^{[q]}\,,
\]
and note that $0\leqslant k\leqslant p-1$. If $k$ is nonzero, then applying $\partial/\partial x_i$ we see that
\[
f^{k-1}g^pJ\ \subseteq\ \mathfrak{m}^{[q]}\,.
\]
Lemma~\ref{lemma:jacobian} and Lemma~\ref{lemma:colon} show that
\[
f^{k-1}g^p\ \in\ \mathfrak{m}^{[q]}+\mathfrak{m}^{(n+1)(q-d+1)}\,.
\]
Since $f^{k-1}g^p\notin\mathfrak{m}^{[q]}$, we must have
\[
\deg f^{k-1}g^p\ \geqslant\ (n+1)(q-d+1)\,.
\]
Using $k\leqslant p-1$ and $\deg g^p\leqslant q(n+1)-pd-p$, this gives
\[
nd-d-n-1\ \geqslant\ p\,,
\]
contradicting the assumption on $p$. It follows that $k=0$, i.e., that $g^p\in\mathfrak{m}^{[q]}$. But then
\[
\left[\frac{g}{(x_0\cdots x_n)^{q/p}}\right]\ =\ 0
\]
in $H^{n+1}_\mathfrak{m}(R)$, which proves the desired injectivity.
\end{proof}

\begin{remark}
Theorem~\ref{thm:highdeg} is equivalent to the following geometric statement: if $X$ is a smooth hypersurface of degree $d$ in $\mathbf{P}^n$, then the map
\[
H^{n-1}(X,{\mathcal{O}}_X(j))\longrightarrow H^{n-1}(X,{\mathcal{O}}_X(jp))\,,
\]
induced by Frobenius map on $X$, is injective for $j<0$ and $p\geqslant nd-d-n$. This statement indeed admits a geometric proof based on the Deligne-Illusie method \cite{DI}. One views the de Rham complex $\Omega^*_{X/K}$ as an ${\mathcal{O}}_{X^{(1)}}$-complex, where $X^{(1)}$ is the Frobenius twist of $X$ over~$K$, and twists it over the latter with ${\mathcal{O}}_{X^{(1)}}(j)$. For $p>n-1$, the Deligne-Illusie decomposition $\Omega^*_{X/K} \simeq \oplus_i \Omega^i_{X^{(1)}/K}[-i]$, which is available as $X$ clearly lifts to $W_2(K)$, reduces the above injectivity statement to proving $H^{n-1-i}(X,\Omega^i_{X/K}(jp)) = 0$ for $i > 0$ and $j < 0$. If $p \geqslant nd-d-n$, this vanishing can be proven using standard sequences (details omitted).
\end{remark}

\section{Calabi-Yau hypersurfaces}

We get to the main theorem; see below for the definition of the Hasse invariant.

\begin{theorem}
\label{thm:CalabiYau}
Let $K$ be a field of characteristic $p>0$, and $n$ a positive integer. Let $f$ be a homogeneous polynomial of degree $n+1$ in $R=K[x_0,\dots,x_n]$, such that the Jacobian ideal of~$f$, i.e., the ideal $(f_{x_0},\dots,f_{x_n})$, is primary to the homogeneous maximal ideal of $R$. Then:
\begin{enumerate}[\quad\rm(1)]
\item $\mu_f(p)=p-h$, where $h$ is an integer with $0\leqslant h\leqslant n-1$,
\item $\mu_f(pq)=p\mu_f(q)$ for all $q=p^e$ with $q\geqslant n-1$.
\item If $p\geqslant n-1$, then $\displaystyle{G_f(z)=\frac{1-hz}{1-pz}}$ \ and \ $\displaystyle{\mathrm{fpt}(f)=1-\frac{h}{p}}$\,, where $0\leqslant h\leqslant n-1$.
\item Set $X=\operatorname{Proj} R/fR$. If $p\geqslant n^2-n-1$, then the integer $h$ in $(1)$ is the order of vanishing of the Hasse invariant on the versal deformation space of $X \subset \mathbf{P}^n$.
\end{enumerate}
\end{theorem}

The deformation space in (4) refers to embedded deformations of $X\subset\mathbf{P}^n$; if $n\geqslant 5$, this coincides with the versal deformation space of $X$ as an abstract variety (see Remark~\ref{rmk:VersalDefHyp}). The following example from~\cite{Hernandez2} shows that all possible values of $h$ from (1) above are indeed attained:

\begin{example}
\label{ex:CalabiYau}
Consider $f=x_0^{n+1}+\dots+x_{n+1}^{n+1}$ over a field of prime characteristic $p$ not dividing $n+1$. Let $h$ be an integer such that $p\equiv h+1\mod n+1$ and $0\leqslant h\leqslant n-1$. Then
\[
\mathrm{fpt}(f)\ =\ 1-h/p\,,
\]
for a proof, see~\cite[Theorem~3.1]{Hernandez2}.
\end{example}

\begin{proof}[Proof of Theorem~\ref{thm:CalabiYau}]
If $\mu_f(p)=p$, then Lemma~\ref{lemma:patterns}\,(1) shows that $\mu_f(q)=q$ for all $q$, and assertions (1--3) follow. Assume that $\mu_f(p)<p$. Lemma~\ref{lemma:bound} gives $\mu_f(p)\geqslant p-n+1$, which proves (1). As $\mu_f(p)\leqslant p-1$, it follows~that
\[
\mu_f(q)\ \leqslant\ q-q/p\qquad\text{ for each }q=p^e\,.
\]
If $\mu_f(pq)<p\mu_f(q)$, then $\mu_f(pq)$ is not a multiple of $p$ by~\eqref{eqn:mu}. Lemma~\ref{lemma:bound} now implies
\[
\mu_f(pq)\ \geqslant\ pq-n+1\,,
\]
and combining with $\mu_f(pq)\leqslant p\mu_f(q)-1\leqslant pq-q-1$, we see that $pq-n+1\leqslant pq-q-1$, i.e., that $q\leqslant n-2$. This completes the proof of (2), and then (3) follows immediately.

The proof of (4) and the surrounding material occupy the rest of this section.
\end{proof}

\subsection*{The Hasse invariant}

We briefly review the construction of the Hasse invariant for suitable families of varieties in characteristic $p$. Fix a proper flat morphism $\pi\colon{\mathcal{X}}\longrightarrow S$ of relative dimension $N$ between noetherian $\mathbf{F}_p$-schemes. Assume that the formation of~$\mathrm{R}^i\pi_*{\mathcal{O}}_{\mathcal{X}}$ is compatible with base change, and that $\omega := \omega_{{\mathcal{X}}/S} := \mathrm{R}^N\pi_*{\mathcal{O}}_{\mathcal{X}}$ is a line bundle; the key example is a family of degree $(n+1)$ hypersurfaces in $\mathbf{P}^n$. The standard diagram of Frobenius twists of $\pi$ takes the shape
\[
\xymatrix{
{\mathcal{X}} \ar@/_/[ddr]_{\pi} \ar@/^/[drr]^{\mathrm{Frob}_{\mathcal{X}}} \ar[dr]^-{\mathrm{Frob}_\pi} & & \\
& {\mathcal{X}}^{(1)} \ar[d]^-{\pi^{(1)}} \ar[r]^-{\mathrm{Frob}_S} & {\mathcal{X}} \ar[d]^-{\pi} \\
& S \ar[r]^-{\mathrm{Frob}_S} & S\,,}
\]
where the square is Cartesian. Our assumption on $\pi$ shows that
\[
\omega_{{\mathcal{X}}^{(1)}/S}\ :=\ \mathrm{R}^N \pi_*{\mathcal{O}}_{{\mathcal{X}}^{(1)}}\ \simeq\ \mathrm{Frob}_S^*\omega\ \simeq\ \omega^p\,.
\]
Using this isomorphism, we define:

\begin{definition}
The \emph{Hasse invariant} $H$ of the family $\pi$ is the element in
\[
\mathrm{Hom}(\omega_{{\mathcal{X}}^{(1)}/S},\omega)\ \simeq\ \mathrm{Hom}(\mathrm{Frob}_S^*\omega,\omega)\ \simeq\ \mathrm{Hom}(\omega^p,\omega)\ \simeq\ H^0(S,\omega^{1-p})\,,
\]
defined by pullback along the relative Frobenius map $\mathrm{Frob}_\pi\colon{\mathcal{X}}\longrightarrow{\mathcal{X}}^{(1)}$.
\end{definition}

\begin{remark}
\label{rmk:hasseflatbasechange}
The formation of the relative Frobenius map $\mathrm{Frob}_\pi\colon{\mathcal{X}}\longrightarrow{\mathcal{X}}^{(1)}$ is compatible with base change on $S$. It follows by our assumption on $\pi$ that the formation of $H$ is also compatible with base change. In particular, given a flat morphism $g\colon S' \longrightarrow S$ and a point $s'\in S'$, the order of vanishing of $H$ at $s'$ coincides with that at $g(s')$. Thus, in proving Theorem~\ref{thm:CalabiYau}\,(4), we may assume that $K$ is perfect.
\end{remark}

To analyze $H$, fix a point $s$ in $S$ and an integer $t\ge0$. Write $t[s]\subset S$ for the order $t$ neighbourhood of $s$, and let $t\!{\mathcal{X}}_s\subset{\mathcal{X}}$ and $t\!{\mathcal{X}}^{(1)}_s\subset{\mathcal{X}}^{(1)}$ be the corresponding neighbourhoods of the fibres of $\pi$ and $\pi^{(1)}$. The map $\mathrm{Frob}_\pi$ induces maps $t\!{\mathcal{X}}_s\longrightarrow t\!{\mathcal{X}}^{(1)}_s$, and hence maps
\[
\phi_t\colon H^N\big(t\!{\mathcal{X}}_{s}^{(1)},\ {\mathcal{O}}_{t\!{\mathcal{X}}_s^{(1)}}\big)\longrightarrow H^N\big(t\!{\mathcal{X}}_s,\ {\mathcal{O}}_{t\!{\mathcal{X}}_s}\big)\,.
\]
The order of vanishing of $H$ at $s$ is, by definition, the maximal $t$ such that this map is zero. In favourable situations, one can give a slightly better description of this integer:

\begin{lemma}
\label{lem:identifyhasse}
If the map
\[
\psi_t\colon H^N\big({\mathcal{X}}_s,\ {\mathcal{O}}_{X_s}\big) \longrightarrow H^N\big(t\!{\mathcal{X}}_s,\ {\mathcal{O}}_{t\!{\mathcal{X}}_s}\big)
\]
induced by $\mathrm{Frob}_{{\mathcal{X}}}$ is nonzero for some $t\leqslant p$, then the minimal such $t$ is $\mathrm{ord}_s H+1$.
\end{lemma}

\begin{proof}
For $t \leqslant p$, by the base change assumption on $\mathrm{R}^N\pi_*{\mathcal{O}}_{{\mathcal{X}}}$, one has
\[
H^N\big({\mathcal{X}}_s,\ {\mathcal{O}}_{X_s}\big)\otimes_{\kappa(s)}{\mathcal{O}}_{t[s]}\ \simeq\ H^N\big(t\!{\mathcal{X}}^{(1)}_s,\ {\mathcal{O}}_{t\!{\mathcal{X}}^{(1)}_s}\big)\,,
\]
where ${\mathcal{O}}_{t[s]}$ is viewed as $\kappa(s)$-algebra via the composite
\[
\CD
\kappa(s)@>{\mathrm{Frob}_S}>>{\mathcal{O}}_{\mathrm{Frob}_S^{-1}[s]}@>{\mathrm{can}}>>{\mathcal{O}}_{t[s]}\,,
\endCD
\]
and the isomorphism is induced by the base change ${\mathcal{X}}^{(1)} \longrightarrow {\mathcal{X}}$ of $\mathrm{Frob}_S\colon S \longrightarrow S$. Hence, for such $t$, by adjunction, the map $\phi_t$ induced by $\mathrm{Frob}_\pi$ is nonzero if and only the map $\psi_t$ induced by $\mathrm{Frob}_{\mathcal{X}}$ is nonzero. But $\mathrm{ord}_s H$ is the maximal integer $t$ with $\phi_t = 0$.
\end{proof}

It is typically hard to calculate $H$, or even bound its order of vanishing. However, for families of Calabi-Yau hypersurfaces, we have the following remarkable theorem due to Deuring and Igusa; see~\cite{Igusa} for $n=2$, and Ogus~\cite[Corollary~3]{Ogus} in general:

\begin{theorem}
\label{thm:ogusigusadeuring}
Let $\pi\colon{\mathcal{X}}\longrightarrow\mathrm{Hyp}_{n+1}$ be the universal family of Calabi-Yau hypersurfaces in $\mathbf{P}^n$. For any point $[Y]\in\mathrm{Hyp}_{n+1}(K)$ corresponding to a smooth hypersurface $Y\subset\mathbf{P}^n$, we have $\mathrm{ord}_{[Y]}(H)\leqslant n-1$ if $n\leqslant p$.
\end{theorem}

Ogus's proof relies on crystalline techniques: he relates $\mathrm{ord}_{[Y]}(H)$ to the relative position of the conjugate and Hodge filtrations on a crystalline cohomology group of $Y$ (following an idea of Katz), and then exploits the natural relation between the Hodge filtration and deformation theory of $Y$. His result will not be used in proving Theorem~\ref{thm:CalabiYau}; in fact, our methods yield an alternative proof of Theorem~\ref{thm:ogusigusadeuring} avoiding crystalline methods under a mild additional constraint on the prime $p$, see Remark~\ref{remark:Ogus}.

\begin{remark}
\label{rmk:VersalDefHyp}
The universal family $\pi\colon{\mathcal{X}}\longrightarrow\mathrm{Hyp}_{n+1}$ is, in fact, versal at $[X]$ if $n \geqslant 5$ so $\mathrm{ord}_{[X]}(H)$, i.e., the order of vanishing of $H$ at $[X]\in\mathrm{Hyp}_{n+1}(K)$, is completely intrinsic to $X$. To see versality, it suffices to show that the map $\mathrm{Hom}(I_X/I_X^2, {\mathcal{O}}_X)\longrightarrow H^1(X,T_X)$ obtained from the adjunction sequence
\[
\CD
0@>>>I_X/I_X^2@>>>\Omega^1_{\mathbf{P}^n}|_X@>>>\Omega^1_X@>>>0
\endCD
\]
is surjective, and that $H^2(X,T_X)=0$. By the long exact sequence, it suffices to show the vanishing of $H^1(X,T_{\mathbf{P}^n}|_X)$ and $H^2(X,T_X)$. The Euler sequence
\[
\CD
0@>>>{\mathcal{O}}_{\mathbf{P}^n}@>>>{\mathcal{O}}_{\mathbf{P}^n}(1)^{\oplus n+1}@>>>T_{\mathbf{P}^n}@>>>0
\endCD
\]
restricted to $X$ immediately shows that $H^i(X,T_{\mathbf{P}^n}|_X) = 0$ for $i = 1,2$ if $n \geqslant 5$; here we use that $H^i(X,{\mathcal{O}}_X(j)) = 0$ for $0 < i < n-1$ and all $j$. The cohomology sequence for the dual of the adjunction sequence then shows that $H^2(X,T_X)=0$.
\end{remark}

\subsection*{The universal family}

Fix notation as in Theorem~\ref{thm:CalabiYau}, with $K$ a perfect field. Let
\[
\mathrm{Hyp}_{n+1}\ := \ \mathbf{P}\big(H^0(\mathbf{P}^n,{\mathcal{O}}_{\mathbf{P}^n}(n+1))^\vee\big)
\]
be the space of hypersurfaces of degree $(n+1)$ in $\mathbf{P}^n$; we follow Grothendieck's conventions regarding projective bundles. Let $\pi\colon{\mathcal{X}}\longrightarrow\mathrm{Hyp}_{n+1}$ be the universal family, and let $\mathrm{ev}\colon{\mathcal{X}}\longrightarrow\mathbf{P}^n$ be the evaluation map. Informally, ${\mathcal{X}}$ parametrizes pairs $(x,Y)$ where $x\in\mathbf{P}^n$ and $Y \in \mathrm{Hyp}_{n+1}$ is a degree $(n+1)$ hypersurface containing $x$. This description shows that $\mathrm{ev}\colon{\mathcal{X}}\longrightarrow\mathbf{P}^n$ is a projective bundle, and we can formally write it as $\mathbf{P}({\mathcal{K}}^\vee)\longrightarrow\mathbf{P}^n$, where ${\mathcal{K}}\in\mathrm{Vect}(\mathbf{P}^n)$ is defined as
\[
{\mathcal{K}}\ :=\ \ker\big(H^0(\mathbf{P}^n,{\mathcal{O}}_{\mathbf{P}^n}(n+1))\otimes{\mathcal{O}}_{\mathbf{P}^n}\ \longrightarrow\ {\mathcal{O}}_{\mathbf{P}^n}(n+1)\big)\,,
\]
and the map is the evident one. The resulting map
\[
{\mathcal{X}}\longrightarrow\mathbf{P}\big(H^0(\mathbf{P}^n,{\mathcal{O}}_{\mathbf{P}^n}(n+1))^\vee\big)
\]
is identified with $\pi$. Our chosen hypersurface $X \in \mathbf{P}^n$ gives a point $[X] \in \mathrm{Hyp}_{n+1}(K)$ with $X' := \pi^{-1}([X])$ mapping isomorphically to $X$ via $\mathrm{ev}$. For an integer $t \geqslant 1$, let $tX \subset\mathbf{P}^n$ and $tX' \subset {\mathcal{X}}$ be the order $t$ neighbourhoods of $X \subset\mathbf{P}^n$ and $X' \subset {\mathcal{X}}$ respectively.

\begin{lemma}
\label{lem:pullbackevinjective}
The map $H^{n-1}(tX,{\mathcal{O}}_{tX}) \longrightarrow H^{n-1}(tX',{\mathcal{O}}_{tX'})$ is injective for all $t$.
\end{lemma}

\begin{proof}
Let $V = H^0(\mathbf{P}^n,{\mathcal{O}}_{\mathbf{P}^n}(1))$, so $\mathbf{P}^n = \mathbf{P}(V)$ and $\mathrm{Hyp}_{n+1} = \mathbf{P}(\mathrm{Sym}^{n+1}(V)^\vee)$. For each~$t$, the sheaf ${\mathcal{O}}_{tX}$ admits a filtration defined by powers of the ideal defining $X \subset tX$, and similarly for $tX'$. The map $tX' \longrightarrow tX$ is compatible with this filtration as it sends $X' \subset tX'$ to~$X \subset tX$. Hence, it suffices to check that the induced map
\[
\phi_j\colon H^{n-1}(X,I_X^j/I_X^{j+1})\longrightarrow H^{n-1}(X',I_{X'}^j/I_{X'}^{j+1})
\]
is injective for each $j \geqslant 0$. Fix an isomorphism $\det(V)\simeq K$ and $f \in R_{n+1}$ defining $X$. These choices determine isomorphisms $I_X \simeq {\mathcal{O}}_{\mathbf{P}^n}(-n-1) \simeq K_{\mathbf{P}^n}$, and hence an isomorphism
\[
H^{n-1}(X,O_X)\ \simeq\ H^n(\mathbf{P}^n,I_X)\ \simeq\ H^n(\mathbf{P}^n,K_{\mathbf{P}^n})\ \simeq\ K\,.
\]
Tensoring the exact sequence $0\longrightarrow I_X\longrightarrow{\mathcal{O}}_{\mathbf{P}^n}\longrightarrow{\mathcal{O}}_X\to0$ with $I_X^j$ and using Serre duality shows that
\begin{align*}
H^{n-1}\big(X,I_X^j/I_X^{j+1}\big)\ &=\ \ker\Big(H^n\big(\mathbf{P}^n,I_X^{j+1}\big)\longrightarrow H^n\big(\mathbf{P}^n,I_X^j\big)\Big)\\
&=\ \mathrm{coker}\Big(H^0\big(\mathbf{P}^n,{\mathcal{O}}_{\mathbf{P}^n}((j-1)(n+1))\big)\longrightarrow H^0\big(\mathbf{P}^n,{\mathcal{O}}_{\mathbf{P}^n}(j(n+1))\big)\Big)^\vee\\
&=\ \mathrm{coker}\Big(\mathrm{Sym}^{(j-1)(n+1)}(V) \stackrel{f}{\longrightarrow} \mathrm{Sym}^{j(n+1)}(V)\Big)^\vee.
\end{align*}

As $X' \subset {\mathcal{X}}$ is a fibre of $\pi$, one has $tX' = \pi^{-1}(t[X])$, where $t[X] \subset \mathrm{Hyp}_{n+1}$ is the order $t$ neighbourhood of $[X] \in \mathrm{Hyp}_{n+1}(K)$. Using flatness of $\pi$ and the aforementioned isomorphism $H^{n-1}(X,O_X) \simeq K$, we get
\[
H^{n-1}(X',\ I_{X'}^j/I_{X'}^{j+1})\ =\ \big(\mathrm{Sym}^j(\mathrm{Sym}^{n+1}(V)/(f))\big)^\vee.
\]
One can check that the pullback $\phi_j$ above is dual to the map induced by the composition map $\mathrm{Sym}^j(\mathrm{Sym}^{n+1}(V)) \longrightarrow \mathrm{Sym}^{j(n+1)}(V)$ by passage to the appropriate quotients. In particular, the dual map is surjective, so $\phi_j$ is injective.
\end{proof}

\begin{proof}[Proof of Theorem~\ref{thm:CalabiYau}\,(4)]
By Remark~\ref{rmk:hasseflatbasechange}, we may assume that the field $K$ is perfect. Fix some $1\leqslant t\leqslant p$. We get a commutative diagram
\[
\CD
X' @>i>> tX' @>j>> pX' @>a>> \mathrm{Frob}_{\mathcal{X}}^*X' @>c>> X'\\
@VVV @VVV @VVV @VVV @VVV\\
X @>k>> tX @>\ell>> pX @>b>> \mathrm{Frob}_{\mathbf{P}^n}^*X @>d>> X
\endCD
\]
Here all vertical maps are induced by $\mathrm{ev}\colon{\mathcal{X}}\longrightarrow\mathbf{P}^n$, the maps $c$ and $d$ are induced by the Frobenius maps on ${\mathcal{X}}$ and $\mathbf{P}^n$ respectively, and $i$, $j$, $k$, $\ell$, $a$ and $b$ are the evident closed immersions; the map $b$ is an isomorphism as $X \subset \mathbf{P}^n$ is a Cartier divisor. In particular, the composite map $X' \longrightarrow X'$ and $X \longrightarrow X$ obtained from each row are the Frobenius maps on~$X'$ and $X$ respectively. Passing to cohomology gives a commutative diagram
\[
\CD
H^{n-1}(X,{\mathcal{O}}_X) @= H^{n-1}(X',{\mathcal{O}}_{X'})\\
@VV{a_t}V @VV{b_t}V \\
H^{n-1}(tX,{\mathcal{O}}_{tX}) @>{c_t}>> H^{n-1}(tX',{\mathcal{O}}_{tX'})
\endCD
\]
where $a_t$ and $b_t$ are induced by the Frobenius maps on $\mathbf{P}^n$ and ${\mathcal{X}}$ respectively, while $c_t$ is injective by Lemma~\ref{lem:pullbackevinjective}. To finish the proof, observe that Lemma~\ref{lem:identifyhasse} shows that $h+1$ is the minimal value of $t$ for which $b_t$ is injective. Since $c_t$ is injective as well, this is also the minimal value of $t$ for which $a_t$ is injective. It follows by Remark~\ref{rmk:fptthickenings} and Theorem~\ref{thm:highdeg} that $\mu_f(p) = p - (h+1) + 1 = p-h$.
\end{proof}

\begin{remark}
\label{remark:Ogus}
We recover Ogus's result, Theorem~\ref{thm:ogusigusadeuring}, for primes $p\geqslant n^2-n-1$; this is immediate from Theorem~\ref{thm:CalabiYau}\,(1) and (4).
\end{remark}

We conclude by giving a de~Rham interpretation for $\mathrm{fpt}(f)$; write $\mathrm{Fil}^\mathrm{conj}_\bullet$ and $\mathrm{Fil}^\bullet_H$ for the increasing conjugate and decreasing Hodge filtrations on $H^{n-1}_\mathrm{dR}(X)$ respectively. Then:

\begin{corollary}
One has $\mathrm{fpt}(f)=1-a/p$, where $a$ is the largest $i$ such that
\[
\mathrm{Frob}_X^* H^{n-1}_\mathrm{dR}(X)\ \simeq\ \mathrm{Fil}^\mathrm{conj}_0(H^{n-1}_\mathrm{dR}(X))\ \subset\ \mathrm{Fil}^i_H (H^{n-1}_\mathrm{dR}(X))\,.
\]
\end{corollary}

\begin{proof}
This follows from Theorem~\ref{thm:CalabiYau} and Ogus's result~\cite[Theorem~1]{Ogus}.
\end{proof}

Note that the integer $a$ appearing above is the $a$ number defined by van~der~Geer and Katsura~\cite{GK} for the special case of a Calabi-Yau family.

\section{Quartic hypersurfaces in $\mathbf{P}^2$}

Our techniques also yield substantive information for hypersurfaces other than Calabi-Yau hypersurfaces; as an example, we include here the case of quartic hypersurfaces in~$\mathbf{P}^2$.

When $f$ defines a Calabi-Yau hypersurface $X$, it is readily seen that the Frobenius action on the vector space $H^{\dim X}(X,{\mathcal{O}}_X)$ is injective if and only if $\mathrm{fpt}(f)=1$, i.e., if and only
\[
\mathrm{fpt}(f)=\mathrm{lct}(f)\,.
\]
For hypersurfaces $X$ of general type, the injectivity of the Frobenius on $H^{\dim X}(X,{\mathcal{O}}_X)$, or even the ordinarity of $X$ in the sense of Bloch and Kato~\cite[Definition~7.2]{Bloch-Kato}---a stronger condition---does not imply the equality of the $F$-pure threshold and the log canonical threshold: for example, for each $f$ defining a quartic hypersurface in $\mathbf{P}^2$ over a field of characteristic $p\equiv 3\mod 4$, we shall see that $\mathrm{fpt}(f)<\mathrm{lct}(f)$; we emphasize that this includes the case of generic hypersurfaces, and that these are ordinary in the sense of Bloch and Kato by a result of Deligne; see~\cite{Illusie}. More generally:

\begin{lemma}
\label{lemma:pigeon}
Let $f$ be a homogeneous polynomial of degree $d$ in $K[x_0,\dots,x_n]$. Then:
\begin{enumerate}[\quad\rm(1)]
\item For each $q=p^e$, one has $\displaystyle{\mu_f(q)\ \leqslant\ \left\lceil\frac{nq+q-n}{d}\right\rceil}\,,$

\item If $nq+q$ is congruent to any of $1,2,\dots,n\mod d$ \ for some $q$, then $\displaystyle{\mathrm{fpt}(f)<\frac{n+1}{d}}$.
\end{enumerate}
\end{lemma}

For quartics in $\mathbf{P}^2$, one has $n=2$ and $d=4$; thus, if $p\equiv 3\mod 4$, then $np+p\equiv1\mod d$, so the $F$-pure threshold is strictly smaller than the log canonical threshold by (2).

\begin{proof}
The pigeonhole principle implies that $f^k\in\mathfrak{m}^{[q]}$ whenever $dk\geqslant (n+1)(q-1)+1$, which proves (1). For (2), suppose $nq+q$ is congruent to any of $1,2,\dots,n\mod d$. Then
\[
\left\lceil\frac{nq+q-n}{d}\right\rceil\ \leqslant \ \frac{nq+q-n+(n-1)}{d}\ =\ \frac{nq+q-1}{d}\,,
\]
and it follows using (1) that $\mu_f(q)<(nq+q)/d$. Thus,
\[
\mathrm{fpt}(f)\ \leqslant\ \mu_f(q)/q\ <\ (n+1)/d\,.\qedhere
\]
\end{proof}

\begin{theorem}
\label{theorem:quartic}
Let $K$ be a field of characteristic $p> 2$. Let $f$ be a homogeneous polynomial of degree $4$ in $K[x_0,x_1,x_2]$, such that the Jacobian ideal of $f$ is $\mathfrak{m}$-primary. Then the possible values for $\mu_f(q)$ and the $F$-pure threshold are:
\begin{align*}
p&\equiv 1\mod 4:\qquad \mu_f(q)=
\begin{cases}
\frac{q(3p-3)}{4p}&\text{ for all }q,\quad\mathrm{fpt}(f)=\frac{3p-3}{4p}\,,\\
\frac{3q+1}{4}&\text{ for all }q,\quad\mathrm{fpt}(f)=\frac{3}{4}\,,\\
\end{cases}\\
p&\equiv 3\mod 4:\qquad \mu_f(q)=
\begin{cases}
\frac{q(3p-5)}{4p}&\text{ for all }q,\quad\mathrm{fpt}(f)=\frac{3p-5}{4p}\,,\\
\frac{q(3p-1)}{4p}&\text{ for all }q,\quad\mathrm{fpt}(f)=\frac{3p-1}{4p}\,.\\
\end{cases}
\end{align*}
\end{theorem}

\begin{proof}
Lemma~\ref{lemma:bound} and Lemma~\ref{lemma:pigeon}\,(1) provide the respective inequalities
\[
\left\lceil\frac{3p-5}{4}\right\rceil\ \leqslant\ \mu_f(p)\ \leqslant\ \left\lceil\frac{3p-2}{4}\right\rceil\,.
\]

If $p\equiv 3\mod 4$, this reads
\[
\frac{3p-5}{4}\ \leqslant\ \mu_f(p)\ \leqslant\ \frac{3p-1}{4}\,,
\]
so there are two possible values for the integer $\mu_f(p)$. The sequence $\{\mu_f(q)/q\}_q$ is constant by Lemma~\ref{lemma:patterns}\,(2), which completes the proof in this case.

If $p\equiv 1\mod 4$, the inequalities read
\[
\frac{3p-3}{4}\ \leqslant\ \mu_f(p)\ \leqslant\ \frac{3p+1}{4}\,.
\]
Again, there are two choices for $\mu_f(p)$. If $\mu_f(p)=(3p-3)/4$, then $\{\mu_f(q)/q\}_q$ is a constant sequence by Lemma~\ref{lemma:patterns}\,(2), whereas if $\mu_f(p)=(3p+1)/4$, then Lemma~\ref{lemma:patterns}\,(1) implies that $\mu_f(q)=(3q+1)/4$.
\end{proof}

\begin{remark}
Similarly, for a quintic $f$ in $K[x_0,x_1,x_2]$ with an $\mathfrak{m}$-primary Jacobian ideal, the possibilities for $\mathrm{fpt}(f)$ are easily determined; we do not list the corresponding $\mu_f(q)$ as these are dictated by $\mathrm{fpt}(f)$. We assume below that $p>5$.
\begin{alignat*}7
p&\equiv1\mod 5:\quad && \mathrm{fpt}(f):\quad &&(3p-3)/5p\,,\ && \text{ or }\ &&4/5\,,\\
p&\equiv2\mod 5:\quad && \mathrm{fpt}(f):\quad &&(3p-6)/5p\,,\ && \text{ or }\ &&(3p-1)/5p\,,\\
p&\equiv3\mod 5:\quad && \mathrm{fpt}(f):\quad &&(3p-4)/5p\,,\ && \text{ or }\ &&(3p^2-7)/5p^2\,,\ && \text{ or }\ &&(3p^2-2)/5p^2\,,\\
p&\equiv4\mod 5:\quad && \mathrm{fpt}(f):\quad &&(3p-7)/5p\,,\ && \text{ or }\ &&(3p-2)/5p\,.
\end{alignat*}
\end{remark}

The log canonical threshold of a smooth quartic in $\mathbf{P}^2$ is $3/4$; except for the case where it equals $3/4$, the denominator of $\mathrm{fpt}(f)$ in Theorem~\ref{theorem:quartic} is $p$. For a quintic as above, if $\mathrm{fpt}(f)\neq\mathrm{lct}(f)$, then the denominator of $\mathrm{fpt}(f)$ is a power of $p$. More generally:

\begin{proposition}
Let $K$ be a field of characteristic $p>0$. Let $f$ be a homogeneous polynomial of degree $d$ in $K[x_0,\dots,x_n]$ with an $\mathfrak{m}$-primary Jacobian ideal. If $p\geqslant nd-d-n$, then either $\mathrm{fpt}(f)=(n+1)/d$, or else the denominator of $\mathrm{fpt}(f)$ is a power of $p$.
\end{proposition}

\begin{proof}
If $\mathrm{fpt}(f)<(n+1)/d$, then there exists an integer $q$ such that $\mu_f(q)/q<(n+1)/d$. But then $\mathrm{fpt}(f)=\mu_f(q)/q$ by Lemma~\ref{lemma:patterns}\,(2).
\end{proof}

\section*{Acknowledgments}

We thank Mircea Musta\c t\u a for raising the question for elliptic curves at the AMS-MRC workshop on Commutative Algebra, Snowbird, 2010, the workshop participants, and the American Mathematical Society. We also thank Johan de Jong for a useful conversation.


\section{Introduction}
\label{Introduction}
The cubic nonlinear Schr\"odinger (NLS) equation,
which describes wave propagation in the presence of
second-order dispersion and cubic (Kerr) nonlinearity,
is one of the most widely researched nonlinear wave models
in physics. It was successfully employed to describe
water wave dynamics \cite{Zakharov84},
nonlinear waves in plasma \cite{Horton96},
Bose-Einstein condensates (BECs) \cite{Dalfovo99},
and pulse propagation in optical waveguides \cite{Hasegawa95}.
The fundamental NLS solitons are
the most ubiquitous solutions of the cubic NLS equation
due to their stability and to the fact that a generic
wave pattern typically evolves into a sequence
of fundamental solitons in the presence of
anomalous dispersion and Kerr nonlinearity.
In the absence of additional physical
processes (perturbations), the fundamental solitons propagate
without any change in their amplitude, group velocity, and shape.
However, the presence of perturbations usually breaks
this ideal picture, by inducing changes in the solitons'
amplitude, group velocity, and shape.
Dissipative perturbations  due to linear and nonlinear
gain or loss are of particular interest, since they are very
common in soliton systems. In optical waveguides, for example,
nonlinear loss or gain arise due to multiphoton absorption or
emission, respectively \cite{Boyd2008}.
Moreover, nonlinear loss and gain play an important role
in many phenomena described by the complex
Ginzburg-Landau (GL) equation \cite{Kramer2002},
such as convection and pattern formation in fluids
and mode-locked lasers.

Despite its success, the single NLS equation is limited to
describing a scalar physical field and is also unsuitable for
handling generic broadband nonlinear wave systems.
If the physical field is a vector, or if the spectra of the
waves in a broadband wave system are concentrated about multiple
widely separated wavelengths, the single NLS equation
should be replaced by a system of coupled-NLS equations.
Indeed, in recent years, coupled-NLS models have been
employed in studies of a wide range of phenomena in fluid
dynamics \cite{Fluids,Onorato2006},
nonlinear optics \cite{Hasegawa95,NLoptics},
multiphase BECs \cite{BECs}, and plasma physics \cite{plasma}.
Coupled-NLS models are particularly useful in describing
broadband wave systems, in which the waves are organized in
sequences (trains) moving with very different group velocities,
such as in crossing seas \cite{Onorato2006},
or in broadband optical waveguide transmission \cite{Hasegawa95}.
Due to the large group velocity differences
in these systems, collisions between pulses from different
sequences are very frequent, and thus play a major role in
the dynamics.

In the current paper, we study the dynamics
of two colliding sequences of fundamental NLS solitons
in the presence of dissipative perturbations due to
linear loss, cubic gain, and quintic loss
(i.e., a GL gain-loss profile). In the absence of perturbations,
the solitons' amplitudes and group velocities do not
change in the collisions. However, the presence of
dissipative interaction due to cubic gain and quintic loss
induces additional amplitude shifts during the collisions,
whose magnitude depends on the initial amplitudes of the two
colliding solitons. Consequently, amplitude dynamics of
solitons in the two sequences become nonlinearly coupled,
and important questions arise regarding the nature of
this dynamics. In three recent studies we showed that
amplitude dynamics of $N$ sequences of colliding solitons
in the presence of dissipative perturbations might be
described by Lotka-Volterra (LV) models for $N$ species,
where the exact form of the LV model depends
on the nature of the perturbations \cite{NP2010a,PNC2010,PC2012b}.
This potential relation between coupled-NLS equations and LV
systems is of great interest, due to the central role
of these seemingly unrelated models in the natural sciences.
While the importance of coupled-NLS equations was explained
in the previous paragraph, LV models are widely used in
environmental sciences to describe population dynamics
\cite{Lotka25,Volterra28,LV_ecology}, in chemistry, to study
dynamics of chemical reactions \cite{LV_chem}, in economics,
to explain interaction between different technologies
\cite{LV_economics}, and in neural networks, to compute neuron
firing rates \cite{LV_nueral}. Due to the importance
of coupled-NLS and LV models, it is essential to further enhance
and extend the understanding of the relation between them.
Moreover, a better understanding of this relation can be used
for stabilizing the propagation, for controlling and tuning
of soliton amplitudes and group velocities, and even for
broadband switching.

In order to give further motivation for the current study
and to clarify the importance of its results, we provide a
brief critical description of our three earlier studies of
the problem. In Ref. \cite{NP2010a}, we studied the propagation
of colliding soliton sequences in the presence of delayed Raman
response. We showed that amplitude dynamics of $N$ soliton
sequences is described by an $N$-dimensional
predator-prey model, but the results were not checked
by numerical simulations with the corresponding coupled-NLS
model. Furthermore, the equilibrium states of the predator-prey
model were found to be centers, i.e., the equilibria were not
asymptotically stable. In order to resolve the problem of
asymptotic stability in {\it the LV model},
we turned to investigate collision dynamics
of soliton sequences in the presence of weak linear gain
and cubic loss \cite{PNC2010}. We showed that in this
case dynamics of soliton amplitudes might be approximately
described by a LV model for competing species with quadratic
interaction terms. Stability analysis of the equilibrium
states of the model yielded conditions on the
physical parameters for stable propagation of two sequences
with equal soliton amplitudes. Numerical simulations with the full
coupled-NLS model confirmed stability at short-to-intermediate
distances, but uncovered an instability at longer
distances \cite{PNC2010}. The latter finding was attributed
to weak (second-order) radiative instability due to the impact
of linear gain on the two sequences. In order to overcome this
destabilizing effect, we turned to study a different setup,
where the soliton sequences propagate in the presence of weak
linear loss, cubic gain, and quintic loss \cite{PC2012b}.
We showed that in this setup, the LV model for amplitude dynamics
contains quadratic and quartic interaction terms. The presence of linear
loss was expected to suppress the weak radiative instability
and by this to enable stable long-distance propagation.
However, numerical solution of the corresponding
perturbed coupled-NLS model still showed instability
at intermediate-to-long propagation distances \cite{PC2012b}.
Thus, the important problem of further stabilization
of the propagation remained unresolved.
Moreover, due to the instability at intermediate
distances, all the previous studies were unable to explore
the possibility of achieving stable on-off and off-on
switching, i.e., the turning on and off of transmission
of one sequence by a fast change of one or more of the
physical parameters. In addition, the earlier studies were
limited in scope to small values of the dissipative coefficients
and to uniform initial distribution of soliton amplitudes
and phases in each sequence.

In the current paper we address these important deficiencies
of the earlier studies. Considering propagation of two colliding
soliton sequences in the presence of a GL gain-loss profile,
we observe stable long-distance transmission over a wide range
of the gain-loss coefficients, including values that are outside
of the perturbative regime. The stable transmission is enabled
by the presence of linear loss, a sufficiently large initial
inter-soliton separation, and the absence of significant
initial position shift between the two soliton sequences.
Moreover, we show that robust on-off and off-on switching
of one soliton sequence can be realized by an abrupt change
in the ratio of the cubic gain and quintic loss coefficients.
In addition, we extend the results to soliton sequences with
periodically alternating phases. Our findings significantly
enhance the surprising relation between propagation of
colliding NLS soliton sequences and population dynamics
in LV models. The results also indicate that the relation
might be further extended to sequences of solitons of
the cubic-quintic GL equation.

The rest of the paper is organized as follows.
In Sec. \ref{models}, we present the coupled-NLS model
for the propagation along with the reduced LV model
for amplitude dynamics. We also discuss the predictions
of the latter model for stability of propagation.
In Sec. \ref{simu}, we present the results of numerical
simulations with the coupled-NLS model for long distance
transmission, as well as for on-off and off-on switching.
In addition, we compare the simulation results
with the predictions of the LV model.
Our conclusions are presented in
Sec. \ref{conclusions}.

\section{Coupled-NLS and Lotka-Volterra models}
\label{models}

We consider propagation of sequences of soliton pulses in
the presence of second-order dispersion, Kerr nonlinearity,
and a GL gain-loss profile, consisting of
linear loss or gain, cubic gain, and quintic loss.
We denote the physical field of the $j$th sequence by $\psi_{j}$.
In the context of propagation of light through optical waveguides,
for example,  $\psi_{j}$ is proportional to
the envelope of the electric field of the $j$th sequence.
In the absence of gain and loss, the propagation
of the $j$th soliton-sequence is described by the
cubic NLS equation
\begin{eqnarray} &&
i\partial_z\psi_{j}+\partial_{t}^2\psi_{j}+2|\psi_{j}|^2\psi_{j}=0,
\label{crosstalk1}
\end{eqnarray}
where we adopt the waveguide optics notation, in which
$z$ is propagation distance, and $t$ is time.
The fundamental soliton solution of the NLS equation
with group velocity $2\beta_{j}$
is  $\psi_{sj}(t,z)=\eta_{j}\exp(i\chi_{j})\mbox{sech}(x_{j})$,
where $x_{j}=\eta_{j}\left(t-y_{j}-2\beta_{j} z\right)$,
$\chi_{j}=\alpha_{j}+\beta_{j}(t-y_{j})+
\left(\eta_{j}^2-\beta_{j}^{2}\right)z$,
and $\eta_{j}$, $y_{j}$, and $\alpha_{j}$
are the soliton amplitude, position, and phase, respectively.

We focus attention on the dynamics of two sequences of fundamental
NLS solitons propagating with group velocities $2\beta_{j}$,
where $j=1,2$. Assuming a large group velocity difference
$|\beta_{1}-\beta_{2}|\gg 1$, the solitons undergo a large number
of fast inter-sequence collisions. We now take into  account
the effects of linear gain-loss, cubic gain, and quintic loss
as well as inter-sequence interaction due to Kerr nonlinearity.
Thus, the propagation is described by the following system
of perturbed coupled-NLS equations:
\begin{eqnarray} &&
i\partial_z\psi_{j}+\partial_{t}^2\psi_{j}+2|\psi_{j}|^2\psi_{j}
+4|\psi_{k}|^2\psi_{j}=
\nonumber \\&&
ig_{j}\psi_{j}/2
+i\epsilon_{3}|\psi_{j}|^2\psi_{j}+2i\epsilon_{3}|\psi_{k}|^2\psi_{j}
\nonumber \\&&
-i\epsilon_{5}|\psi_{j}|^4\psi_{j}-3i\epsilon_{5}|\psi_{k}|^4\psi_{j}
-6i\epsilon_{5}|\psi_{k}|^2|\psi_{j}|^2\psi_{j},
\label{crosstalk2}
\end{eqnarray}
where $j=1,2$, $k=1,2$, $g_{j}$ is the linear
gain-loss coefficient for the $j$th sequence,
and $\epsilon_{3}$ and $\epsilon_{5}$ are the cubic gain
and quintic loss coefficients, respectively.
The term $4|\psi_{k}|^2\psi_{j}$ in Eq. (\ref{crosstalk2})
describes inter-sequence interaction due to Kerr nonlinearity,
while $ig_{j}\psi_{j}/2$, $i\epsilon_{3}|\psi_{j}|^2\psi_{j}$, and
$-i\epsilon_{5}|\psi_{j}|^4\psi_{j}$ correspond to
intra-sequence effects due to linear gain-loss, cubic gain and quintic loss,
respectively. In addition, the term $2i\epsilon_{3}|\psi_{k}|^2\psi_{j}$
represents inter-sequence effects due to cubic gain, while
$-3i\epsilon_{5}|\psi_{k}|^4\psi_{j}$ and
$-6i\epsilon_{5}|\psi_{k}|^2|\psi_{j}|^2\psi_{j}$
describe dissipative inter-sequence interaction due to quintic loss.

In Ref. \cite{PC2012b}, we showed that under certain assumptions,
amplitude dynamics of solitons in the two sequences can be approximately
described by a LV model for two species with quadratic
and quartic interaction terms. The derivation of the LV model was based
on the following assumptions.
(1) The temporal separation $T$ between
adjacent solitons in each sequence is a large constant: $T \gg 1$.
In addition, the amplitudes are equal for all solitons from the
same sequence, but are not necessarily equal for solitons from
different sequences. (2) The pulses circulate in a closed
loop, e.g., in an optical waveguide ring.
(3) As $T\gg 1$, the pulses in each sequence are
temporally well-separated. As a result, intra-sequence interaction is
exponentially small and is neglected. (4) High-order effects due to
collision-induced frequency shift and emission of radiation
are also neglected.

Since the soliton sequences are periodic,
the amplitudes of all pulses in a given sequence undergo the same
dynamic evolution. Taking into account collision-induced
and single-pulse amplitude changes, we obtain the following equation
for the rate of change of the amplitude of the solitons in the $j$th
sequence $\eta_{j}$ \cite{PC2012b}:
\begin{eqnarray} &&
\frac{d \eta_{j}}{dz}=
\eta_{j}\left[g_{j}+\frac{4}{3}\epsilon_{3}\eta_{j}^{2}
-\frac{16}{15}\epsilon_{5}\eta_{j}^{4}
+\frac{8}{T}\epsilon_{3}\eta_{k}
-\frac{8}{T}\epsilon_{5}\eta_{k}\left(2\eta_{j}^{2}+\eta_{k}^{2}\right)
\right],
\label{crosstalk3}
\end{eqnarray}
where $j=1,2$. Note that Eq. (\ref{crosstalk3}) can be described as
a LV model for two species with quadratic cooperation terms
and quartic competition terms.

The results described in the previous paragraph indicate that
there is a relation between collision-induced dynamics of soliton
sequences in nonlinear waveguides with a GL gain-loss profile
and dynamics of population size of species with nonlinear
competition and cooperation. Alternatively, one may speak
about a connection between the
perturbed coupled-NLS model (\ref{crosstalk2})
and the LV model (\ref{crosstalk3}). This relation
along with knowledge about properties of LV models
can be used to develop ways for controlling the dynamics
of the colliding solitons. A straightforward way for
achieving this goal is by tuning of the
linear gain-loss coefficients $g_{j}$. In optical waveguide
systems, for example, this can be realized  by adjusting
the linear amplifier gain.

As a particular example, we consider the important case where the
values of the $g_{j}$ coefficients are chosen such that the LV model
(\ref{crosstalk3}) has a steady state $(\eta,\eta)$ with equal
amplitudes for both sequences. This case is of special interest,
since in optical waveguide systems it is usually easier to control
the transmission when the amplitudes of all pulses are equal.
Requiring that $(\eta,\eta)$ is a steady state, we obtain
$g_{j}=4\epsilon_{5}\eta\left(-\kappa\eta/3+4\eta^{3}/15
-2\kappa/T+6\eta^{2}/T\right)$, where
$\kappa=\epsilon_{3}/\epsilon_{5}$ and $\epsilon_{5}\ne 0$.
Substituting this relation into Eq. (\ref{crosstalk3}),
we arrive at the following LV model for amplitude dynamics:
\begin{eqnarray} &&
\frac{d \eta_{j}}{dz}=
\epsilon_{5}\eta_{j}
\left\{\frac{4\kappa}{3}(\eta_{j}^{2}-\eta^{2})
-\frac{16}{15}(\eta_{j}^{4}-\eta^{4})
+\frac{8\kappa}{T}(\eta_{k}-\eta)
-\frac{8}{T}\left[\eta_{k}\left(2\eta_{j}^{2}+\eta_{k}^{2}\right)
-3\eta^{3}\right]\right\}.
\label{crosstalk4}
\end{eqnarray}
We note that  $(\eta,\eta)$ and $(0,0)$ are equilibrium points
of Eq. (\ref{crosstalk4}) for any positive values of $\eta$,
$\kappa$, and $T$.

Let us discuss the predictions of the LV model (\ref{crosstalk4})
regarding stability of propagation of the two soliton sequences.
For concreteness, we choose $\eta=1$, and require that the
equilibrium points at $(1,1)$ and $(0,0)$ are stable nodes (sinks).
The stability requirement on $(0,0)$, which is important
for suppression of weak instability due to radiation emission,
leads to $g_{j}<0$ for $j=1,2$, i.e., the solitons are subject to
net linear loss. The above stability requirements yield the
following conditions on $T$ and $\kappa$ \cite{PC2012b}:
$(4T+90)/(5T+30)<\kappa<(8T+135)/(5T+15)$ for $3 < T < 60/17$,
and $(4T+90)/(5T+30)<\kappa<(8T-15)/(5T-15)$ for $T \ge 60/17$.

\section{Numerical simulations with the coupled-NLS model}
\label{simu}
As explained above, the LV model (\ref{crosstalk4}) provides
an approximate description for collision-induced
dynamics of circulating soliton sequences.
The propagation is fully described by the coupled-NLS system
(\ref{crosstalk2}) with periodic boundary conditions.
A key question about the collision-induced dynamics
concerns the possibility to suppress the instability
observed in numerical simulations in Ref. \cite{PC2012b}.
Here we address this important question.
More specifically, we show that the instability can be effectively
mitigated by using a larger temporal separation value T, compared
with the one used in \cite{PC2012b}, and by eliminating the
initial inter-sequence position-shift, which was introduced in
\cite{PC2012b} to avoid incomplete collisions.
Thus, the numerical simulations setup in the current paper
consists of the coupled-NLS model (\ref{crosstalk2}) with
periodic boundary conditions and an initial condition in the form of
two periodic sequences of $2J+1$ {\it overlapping} solitons
with amplitudes $\eta_{j}(0)$ and zero phase:
\begin{eqnarray} &&
\psi_{j}(t,0)\!=\!\sum_{k=-J}^{J}
\frac{\eta_{j}(0)\exp[i\beta_{j}(t-kT)]}
{\cosh[\eta_{j}(0)(t-kT)]},
\label{crosstalk5}
\end{eqnarray}
where $j=1,2$. The temporal separation $T$ is taken as $T=20$,
compared with $T=10$ that was used in Ref. \cite{PC2012b}.
Four different values of $\epsilon_{5}$
are used: $\epsilon_{5}=0.01$, $\epsilon_{5}=0.06$,
$\epsilon_{5}=0.1$, and $\epsilon_{5}=0.5$,
compared with $\epsilon_{5}=0.01$ that was
considered in Ref. \cite{PC2012b}.
Note that a value of $\epsilon_{5}=0.5$ is outside
of the parameter range, where the perturbative approach
leading to the LV model (\ref{crosstalk3}) is expected to be valid.
The other physical parameters in the simulations
are taken as  $\kappa=1.5$, $\eta=1$, $\beta_{1}=0$,
$\beta_{2}=40$, and $J=2$.

The $z$-dependence of $\eta_{1}$ obtained by numerical
solution of Eq. (\ref{crosstalk2}) with initial amplitudes
$\eta_{1}(0)=1.25$ and $\eta_{2}(0)=0.7$ is shown in
Fig. \ref{fig1}(a) for  $\epsilon_{5}=0.01$, $\epsilon_{5}=0.06$,
and $\epsilon_{5}=0.1$, and in Fig. \ref{fig1}(b) for $\epsilon_{5}=0.5$.
The predictions of the LV model (\ref{crosstalk4}) are also presented.
The agreement between the coupled-NLS simulations and the predictions
of the LV model are excellent. Moreover, in all
four cases $\eta_{1}(z)$ tends to
the equilibrium value of $\eta=1$ and the approach to
$\eta=1$ is faster as $\epsilon_{5}$ increases.
Similar results are obtained for $\eta_{2}(z)$ and
for other initial amplitude values.
Note that the distances over which stable propagation is
observed are larger by a factor of 6.67 compared with the
the distances in Ref. \cite{PNC2010}, and by a factor of 2
or more compared with the distances in Ref. \cite{PC2012b}.
The fact that the predictions of the LV model hold even
for $\epsilon_{5}$ values as large as 0.5, i.e., outside
of the perturbative regime, is quite surprising.
Furthermore, as shown in the inset of Fig. \ref{fig1}(b),
the shape of the two soliton trains is retained during
the propagation despite of the magnitude of $\epsilon_{5}$
and the large number of collisions.

\begin{figure}[ptb]
\begin{tabular}{cc}
\epsfxsize=6.0cm  \epsffile{fig_1a.eps} \\
\epsfxsize=6.0cm  \epsffile{fig_1b.eps}
\end{tabular}
\caption{(Color online)
(a) The $z$-dependence of $\eta_{1}$ for various $\epsilon_{5}$ values
and $T=20$, $\kappa=1.5$, $\eta_{1}(0)=1.25$, and $\eta_{2}(0)=0.7$.
The solid black, dotted blue, and dashed red lines correspond to
$\eta_{1}(z)$ values obtained by numerical solution of the
coupled-NLS model (\ref{crosstalk2}) with $\epsilon_{5}=0.01$,
$\epsilon_{5}=0.06$, and $\epsilon_{5}=0.1$, respectively.
The dashed-dotted green, short dashed purple, and short dashed-dotted
orange lines represent $\eta_{1}(z)$ values predicted by the LV model
(\ref{crosstalk4}) with  $\epsilon_{5}=0.01$, $\epsilon_{5}=0.06$,
and $\epsilon_{5}=0.1$.
(b) The z-dependence of $\eta_{1}$ and $\eta_{2}$ for $\epsilon_{5}=0.5$
and the same values of $T$, $\kappa$, $\eta_{1}(0)$, and $\eta_{2}(0)$
as in (a). The solid red and dashed blue lines represent
$\eta_{1}(z)$ and $\eta_{2}(z)$ as obtained by numerical solution of
Eq. (\ref{crosstalk2}), while the dashed and short dashed black lines
correspond to $\eta_{1}(z)$ and $\eta_{2}(z)$ values as obtained by
the LV model (\ref{crosstalk4}). The inset shows the final pulse
patterns $|\psi_{1}(t,z_{f})|$ (solid red) and $|\psi_{2}(t,z_{f})|$
(dashed blue) obtained in the coupled-NLS simulation.}
\label{fig1}
\end{figure}

The observation of long-distance propagation of the two
soliton sequences opens the way for studying broadband
on-off and off-on switching, i.e., the turning on and
off of one of the propagating pulse sequences. The two
switching scenarios are based on bifurcations of the
steady state $(1,1)$ of the LV model. More specifically,
in on-off switching, the value of the parameter $\kappa$
is abruptly increased at the switching distance $z_{s}$
from $\kappa_{i}<\kappa_{c}$ to $\kappa_{f}>\kappa_{c}$,
where $\kappa_{c}=(8T-15)/(5T-15)$,
such that the steady state $(1,1)$ becomes unstable,
while another steady state at $(\eta_{s},0)$ is stable.
As a result, soliton amplitudes in sequences 2 and 1
tend to 1 for $z<z_{s}$, but tend to 0 and $\eta_{s}$,
respectively, for $z>z_{s}$. This means that transmission
of sequence 2 is effectively turned off at $z_{s}$.
Off-on switching is realized in a similar manner,
by decreasing the value of $\kappa$ at the switching distance
$z_{s}$ from $\kappa_{i}>\kappa_{c}$ to $\kappa_{f}<\kappa_{c}$
such that the steady state $(1,1)$ becomes stable.
Consequently, for $z<z_{s}$, the amplitudes of solitons in sequences
1 and 2 tend to $\eta_{s}$ and 0, respectively,
but for $z>z_{s}$, both amplitudes tend to 1.
Thus, in this case transmission of sequence 2 is turned
on at $z_{s}$.

In order to check if these switching scenarios can be realized
with sequences of colliding solitons, we numerically solve
Eq. (\ref{crosstalk2}) with initial pulse patterns of the
form (\ref{crosstalk5}), where $T=20$, $\eta_{1}(0)=1.05$,
and $\eta_{2}(0)=0.9$. Since $T=20$, the critical
$\kappa$-value is $\kappa_{c}=29/17$.
As an example, the switching distance is taken as $z_{s}=175$.
In on-off switching, $\epsilon_{5}=0.1$, and $\kappa$ is
increased from $1.5$ to $2.0$ at $z_{s}$.
In off-on switching, we use $\epsilon_{5}=0.01$
and $\kappa=2.0$ for $z \le 175$, and $\epsilon_{5}=0.1$
and $\kappa=1.5$ for $z > 175$. The results of the
simulations for $\eta_{j}(z)$ are shown in
Fig. \ref{fig2} (a) and (b) for on-off and off-on switching,
respectively. Also shown is the prediction of the LV
model (\ref{crosstalk4}). The overall agreement
between the coupled-NLS simulations and
the predictions of the LV model is very good for both
switching scenarios. The only exception is for small values
of $\eta_{2}(z)$ in on-off switching, where the LV model
underestimates the value obtained with the coupled-NLS model.
This is due to the fact that in this regime, the linear loss
term is comparable to or larger than the Kerr nonlinearity term,
leading to the breakdown of the perturbative description.
Despite of this fact, the amplitudes of the solitons
in sequence 2 do tend to 0, and thus,
the on-off switching is fully realized.
Based on the good agreement between coupled-NLS simulations
and predictions of the LV model we conclude that on-off and
off-on switching can be realized for the two sequences
of colliding solitons. It should be noted that stable
off-on switching can be achieved only when
$\eta_{2}(z_{s})$ is larger than some threshold
value $\eta_{th}$. For the parameter values used in our
simulations, $\eta_{th}=0.65$. Thus, in order to implement
off-on switching with the current setup of two soliton sequences,
the value of the decision level distinguishing between 0 and 1
states should be set larger than $\eta_{th}$.

\begin{figure}[ptb]
\begin{tabular}{cc}
\epsfxsize=6.0cm  \epsffile{fig_2a.eps} \\
\epsfxsize=6.0cm  \epsffile{fig_2b.eps}
\end{tabular}
\caption{(Color online)
The $z$-dependence of soliton amplitudes in on-off switching (a),
and in off-on switching (b) with $z_{s}=175$ and parameter values
as described in the text. The red circles and blue squares
represent $\eta_{1}(z)$ and $\eta_{2}(z)$ values obtained
with the coupled-NLS model (\ref{crosstalk2}),
while the solid and dashed black lines correspond to
$\eta_{1}(z)$ and $\eta_{2}(z)$ values obtained with the
LV model (\ref{crosstalk4}).}
\label{fig2}
\end{figure}

All the propagation setups discussed so far are limited,
in the sense that all amplitudes and phases within each soliton
sequence are equal. It is therefore important to extend the results
to more general setups. Possible extension can be achieved by launching
soliton sequences with periodically alternating phases.
Indeed, the LV model (\ref{crosstalk4}) is independent
of the soliton phases, and as a result,
amplitude dynamics for sequences with periodically
alternating phases is expected to be the same as the
dynamics for soliton sequences with a uniform phase pattern.
On the other hand, it is known that intra-sequence soliton
interaction strongly depends on the relative phase
difference between the solitons. The interplay between this
interaction and other effects, such as loss or radiation
emission might induce instabilities, which would lead to
the breakdown of the LV model description.
For this reason, it is important to check whether the predictions
of the LV model  (\ref{crosstalk4}) do apply for soliton
sequences with periodically alternating phases.
For this purpose, we numerically solve
Eq. (\ref{crosstalk2}) with pulse sequences consisting of two
solitons each with initial amplitudes $\eta_{1}(0)=1.05$ (sequence 1)
and $\eta_{2}(0)=0.9$ (sequence 2).
Two sets of initial phases are used:
$\alpha_{j1}(0)=0$ and $\alpha_{j2}(0)=\pi$ [set(a)];
and $\alpha_{j1}(0)=0$ and $\alpha_{j2}(0)=\pi/2$ [set(b)].
The other physical parameter values are taken as
$T=20$, $\kappa=1.5$, and $\epsilon_{5}=0.01$.
The results of the numerical simulations for
$\eta_{j}(z)$ are shown in Fig. \ref{fig3}
along with the predictions of the LV model.
As can be seen, for both sets of initial phases
the agreement between full-scale simulations
and the predictions of the LV model is excellent
over the entire propagation distance.
Furthermore, the soliton amplitudes tend to the
equilibrium value $\eta=1$, i.e., the propagation
is indeed linearly stable. Thus, our numerical
simulations demonstrate that the predictions of
the LV model (\ref{crosstalk4}) can be applied to
generic setups with periodic phase patterns.

\begin{figure}[ptb]
\begin{tabular}{cc}
\epsfxsize=6.0cm  \epsffile{fig_3a.eps}  \\
\epsfxsize=6.0cm  \epsffile{fig_3b.eps}
\end{tabular}
\caption{(Color online)
The $z$-dependence of pulse amplitudes for soliton sequences
with alternating initial phases and parameter values $T=20$,
$\kappa=1.5$, $\epsilon_{5}=0.01$, $\eta_{1}(0)=1.05$,
and $\eta_{2}(0)=0.9$. (a) $\eta_{j}$ vs $z$ for
$\alpha_{jk}(0)=0$, $\alpha_{jk+1}(0)=\pi$.
(b) $\eta_{j}$ vs $z$ for $\alpha_{jk}(0)=0$, $\alpha_{jk+1}(0)=\pi/2$.
The solid red and dashed blue lines represent
$\eta_{1}(z)$ and $\eta_{2}(z)$ values obtained
with the coupled-NLS model (\ref{crosstalk2}),
while the solid and dashed black lines correspond to
$\eta_{1}(z)$ and $\eta_{2}(z)$ values obtained with the
LV model (\ref{crosstalk4}).}
\label{fig3}
\end{figure}

\section{Conclusions}
\label{conclusions}
In summary, we investigated long-distance propagation
as well as on-off and off-on switching
of two colliding soliton sequences in the presence
of second-order dispersion, Kerr nonlinearity, and a GL gain-loss
profile. Using coupled-NLS simulations and a reduced LV model,
we showed that stable long-distance propagation
is possible for a wide range of the gain-loss coefficients,
including values that are outside of the perturbative regime.
The stable propagation is enabled by the presence of linear loss,
the choice of a sufficiently large initial inter-pulse
separation within each sequence, and the absence
of significant initial time-delay between the two sequences.
Furthermore, we found that robust on-off and off-on switching
of one of the propagating sequences can be achieved by an abrupt
change in the ratio of cubic gain and quintic loss coefficients.
Additionally, we demonstrated that the results can be extended
to pulse sequences with periodically alternating phases.
Our study significantly enhances the recently found relation
between collision-induced dynamics of sequences of NLS solitons
and evolution of population sizes in LV models. Moreover, it
indicates that the relation might be further extended to sequences
of solitary waves of the cubic-quintic GL equation, at least
in some regions of parameter space. Finally, our results might
find useful applications in nonlinear waveguides with saturable
absorption for stable energy equalization and broadband switching.


\section{Introduction}

One of the principal goals of modern cosmology is to characterize the statistical properties of the primordial density fluctuations, i.e. the seeds of the present large-scale structure. As widely presumed, the initial perturbations are associated with quantum properties of an inflationary field \citep{guth1981}. If this model is correct, the primordial fluctuations should be overwhelmingly Gaussian \citep{bardeen,bond} along with the small temperature
fluctuations of the Cosmic Microwave Background (CMB) sky.

Gaussianity is the most fundamental prediction of inflation.
Randomness of the complex phases of the harmonic coefficients
of small CMB temperature fluctuations provides natural constraints, since
departures from Gaussian behavior typically cause deviations from
randomness \citep{colesNat}. There are several methods
constraining non-Gaussianity from phase information: phase mapping and uniformity tests \citep{chiang2002,chiang2004}, Shannon entropy of phases \citep{chiang2000}, surrogates \citep{raeth2010}, random walks \citep{stannard2005,hansen2011}, etc. These have been applied to WMAP all-sky maps, and most
recently to Planck \citep{planck_isotropy}. In some cases, non-Gaussian residuals have been detected \citep{chiang2003,naselsky2005},
although no primordial non-Gaussianity has been found with any
certainty.

Other studies, such as \cite{land2005,land2005b,land2006,copi2004b,copi2006} and \cite{bielewicz2005} defined directions on a sphere at each $\ell$ to construct estimators constraining unusual alignments and correlations in the harmonic series representing the CMB maps. Several ``anomalies'' and alignments were
identified, and several tests have been performed to explore their origin
\citep{francis2010,frommert2010,rassat2013}.
These marginally significant anomalies were originally detected in WMAP,
and recently confirmed in Planck \citep{planck_isotropy}.

Complex phases correspond to a unit vector in the complex plane, where the U(1) group acts as a rotation.
Based on this observation we generalize the usual U(1) phases for the group SO(3), relevant to the CMB or any full-sky map, as unit vectors in
$(2\ell+1)$ dimensional representation spaces. These Generalized Phases in the sense of SO(3) respond to SO(3) rotations analogously to complex phases responding to U(1) rotations.
In the rest of this paper we only deal with the SO(3) group, therefore without ambiguity we can call them Generalized Phases, or GPs, hereafter.
For an isotropic Gaussian field, they correspond to a random direction by symmetry, represent
most of the information beyond the measured power spectrum, and they are independent from it. Nevertheless, two observations
of the same CMB realization should have exactly the same phases.
The principal aim of this work is to use this simple property
to construct a rigorous and concise $\ell$-by$\ell$ comparison of WMAP and Planck maps that emphasizes information beyond the
power spectrum. In particular, we will characterize coherence of two maps by the
angle of the unit vectors corresponding to their GPs, that also
corresponds to a correlation coefficient in harmonic space.

We organize this paper as follows. In Section 2 we describe the data we used, and introduce our methods including theoretical expectations, simulations and measurements. In Section 3 our results, and statistical significances of our findings are presented. Finally, we briefly summarize our results in Section 4.
The appendix contains derivations of formulae used in the main text.

\section{Data and methods}

To quantify the coherence of WMAP and Planck we first prepare maps
of the same resolution.
The WMAP team provides $N_{side}=512$ CMB temperature maps, therefore we choose this as our base resolution.
The Planck CMB products have higher resolution, $N_{side}=2048$, thus
we downgraded Planck maps using HEALPIX \citep{healpix} for $N_{side}=512$.
We also used the $N_{side}=512$ WMAP9 Temperature Analysis Mask that
leaves 78\% of the sky for our analysis.

For WMAP, we used the Q,V,W frequency bands downloaded from the LAMBDA website \footnote{\texttt{http://lambda.gsfc.nasa.gov/}}, using both original and foreground reduced versions \citep{wmap,bennett2012}. For Planck, we downloaded the
NILC and Smica CMB maps \citep{pla2013} from the Planck Legacy Archive \footnote{\texttt{http://www.sciops.esa.int/index.php?project=PLANCK}}. They
already have galactic foregrounds and known point sources removed.

\subsection{Generalized Phases}

The CMB temperature fluctuations can be expanded into
spherical harmonics:
\begin{equation}
\frac{\Delta T}{\bar T} (\vartheta,\varphi) = \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{\ell} a_{\ell m} Y_{\ell m}(\vartheta,\varphi)
\label{eq_dT}
\end{equation}
Phases are defined by complex $a_{\ell m}$ coefficients of CMB multipoles as follows
\begin{equation}
a_{\ell m} = |a_{\ell m}| \cdot \exp(i\phi_{\ell m})
\label{eq35}
\end{equation}
These Fourier phases generate rotations around the $z$-axis,
corresponding to to the U(1) subgroup of the full SO(3) symmetry of the
harmonic coefficients. For Gaussian random fields (GRF),
Fourier phases are random and uniformly distributed between 0 and $2\pi$.
Testing the randomness of these phases therefore provides an
interesting diagnostic of the Gaussianity of the fluctuation field \citep{colesNat}. Note that the power at each $\ell$ and these phases do not fully determine the random field.

To generalize complex Fourier phases, we first build
$(2\ell+1)$ dimensional vectors using real and imaginary parts of $a_{lm}$ coefficients:
\begin{equation}
\varepsilon_{\ell}=(a_{\ell 0}/\sqrt{2}, Re[a_{\ell 1}], .... Re[a_{\ell \ell}], Im[a_{\ell 1}], .... Im[a_{\ell \ell}])
\label{xl}
\end{equation}
These vectors contain all the information due to the reality of the underlying
random field. For a Gaussian field, this is a random vector, with
each elements of $\varepsilon_{\ell}$ having a variance of $C_\ell/2$.
Generalized Phases are now defined as $(2\ell+1)$ dimensional unit vectors
\begin{equation}
\hat{\varepsilon}_{\ell} = \frac{\varepsilon_{\ell}}{\sqrt{\sum_{k}{ {\varepsilon^{2}_{\ell,k}} }}}.
\label{xl2}
\end{equation}
As $a_{\ell m}$ coefficients of different multipoles are independent, GPs are uncorrelated for a Gaussian distribution. Moreover, they follow uniform distributions over the sphere $S^{2\ell}$ for each $\ell$ \citep{cai2013}.
The statistics of GPs contain information complementary to the
power spectrum, and for mildly non-Gaussian distributions, they
should contain most of the non-Gaussian information. If the power
and the GPs are given, the realization of a random field is fully constrained.

In this work, we compare Generalized Phases to quantify the (generalized)
phase coherence of WMAP and Planck maps, i.e. the $\ell$-by$\ell$
coherence of the maps beyond and independently of the match of their
power spectra.

To quantify the coherence of the two maps, we calculated dot products of unit vectors defined by individual datasets at each $\ell$ multipole as
\begin{equation}
\cos\Theta_\ell = \sum_{k}\hat{\varepsilon}^{\rm Planck}_{\ell,k} \cdot \hat{\varepsilon}^{\rm WMAP}_{\ell,k}.
\label{scal}
\end{equation}
\begin{figure}
\begin{center}
\includegraphics[width=90mm]{fig1.eps}
\caption{Distributions of angles between random unit vectors in $(2\ell+1)$ dimensions. We compare analytic expectations and simulations, and test properties of galactic masks, as well. We present $\ell=2,5,25,75$ cases. These illustrate that angles are concentrated around $\pi /2$. This concentration becomes stronger as the dimension $n$ grows. These results are insensitive to the galactic
mask as long as the unit vectors are truly random.}
\label{sim}
\end{center}
\end{figure}
\subsection{Random angle statistics in $n$ dimensions}

Angles between Generalized Phases of two uncorrelated datasets - e.g.  CMB realizations - fluctuate around  $\pi/2$, their distributions are characterized by analytic formulae \citep{cai2013}. When the dimension $n=2\ell+1$ is fixed,
the distribution of angles has a density function given by
\begin{equation}
h(\Theta) = \frac{1}{\sqrt{\pi}} \frac{\Gamma(\frac{n}{2})}{\Gamma(\frac{n-1}{2})} \cdot \sin^{n-2} \Theta.
\label{htheta}
\end{equation}
Note that if $n$ = 2, $h(\Theta$) is the uniform density on [0,$\pi$]. When $n$ $\geq$ 3, h($\Theta$) is a unimodal  distribution with peak position of $\Theta$ = $\pi/2$. The concentration around $\pi/2$ becomes stronger as $n$ grows, since $\sin^{n-2}\Theta$ is driven to zero quickly for $\Theta \neq \pi /2$ \citep{cai2013}. This means that uncorrelated vectors in high dimensions tend to be perpendicular. As expected, in large dimensions, the distribution tends to a Gaussian distribution centered on  $\pi/2$.
In Fig. \ref{sim} we show estimates of distributions of angles between unit vectors in higher dimensions. We simulated 500 CMB skies to test Eq. \eqref{htheta}. Simulations were made by using WMAP9 cosmological parameters, and WMAP9 noise.
We randomly choose 10,000 pairs of CMB simulations, and calculate Generalized Phases. Four examples of $\ell=2,5,25,75$ illustrate that individual distributions of angles between random unit vectors in $(2\ell+1)$ dimensions follow Eq. \eqref{htheta} closely. We checked that these results hold up to $\ell=1535$, the
maximum we can measure with our maps.

\begin{figure}
\begin{center}
\includegraphics[width=80mm]{fig0a.eps}
\includegraphics[width=80mm]{fig0b.eps}
\caption{We show original (amplitude-shuffled) version of a simulated CMB map on the top (bottom). These CMB maps have the same phases and pseudo power spectra, but different GPs.}
\label{GPs}
\end{center}
\end{figure}

We repeated our measurements on masked CMB skies using WMAP9 Temperature Analysis Mask. According to Fig. \ref{sim}, and perhaps somewhat surprisingly,
no difference was found. While galactic mask strongly affects statistical analysis of normal phases \citep{chiang2007}, the distribution of $\Theta$ is insensitive to the mask. The CMB mask is centered on $\vartheta=\pi/2$ in the spherical coordinate, which causes strong phase correlation only among phases of $\ell \approx m$.

It is useful to consider the closely related
correlation coefficient $\cos \Theta = C_\ell^{AB}/\sqrt{C_\ell^A C_\ell^B}$ in
addition to $\Theta$. In this case, $x = (\cos \Theta + 1)/2$ follows the Beta distribution on $[0,1]$, i.e. $ x^{\alpha -1}(1-x)^{\beta - 1} / B(\alpha,\beta)$ with parameters $\alpha = \beta = (n-1)/2 = \ell$. Thus the exact first two moments of $\cos \Theta$ are
\begin{equation}
\left \langle {\cos \Theta}\right \rangle  = 0,\quad \left\langle\cos^2\Theta\right \rangle = \frac{1}{n}
\end{equation}

\newcommand{\textrm{CMB}}{\textrm{CMB}}

Finally, we quantified the resolving power of GPs by the following procedure. We shuffled $|a_{\ell m}|$ amplitudes of a simulation for a given $\ell$, keeping both pseudo power spectrum and phases unchanged. Fig. \ref{GPs} shows the original and the "shuffled" CMB maps. We measured $\Theta_\ell$ angles between GPs of the maps (Fig. \ref{GPs2}), finding values fluctuating around $\Theta_\ell \approx 38^{\circ}$. We  integrated the Gaussian distributions of the $a_{\ell m}$'s to find the
average value $\avg{\cos \Theta_\ell} = \pi / 4$. This corresponds to $\Theta_\ell = 38.24^{\circ}$, i.e. $78.5\%$ correlation.

\begin{figure}
\begin{center}
\includegraphics[width=90mm]{shuffled.eps}
\caption{Measured $\Theta_\ell$ angles of GPs of the original and $|a_{\ell m}|$-shuffled maps are illustrated. See text for details. }
\label{GPs2}
\end{center}
\end{figure}

\subsection{CMB and noise}

WMAP and Planck measurements of the CMB sky contain noise. This noise induces a rotation of the unit vectors $\hat{\varepsilon}^{\rm CMB}_{l}$ on the $2\ell$ dimensional sphere. Assuming full sky coverage and isotropic Gaussian noise, these rotations will only depend on the respective spectra of the $\textrm{CMB}$ and that of the noise. The angles obey
\newcommand{\begin{equation}}{\begin{equation}}
\newcommand{\end{equation}}{\end{equation}}
\newcommand{\textrm{noise}}{\textrm{noise}}
\newcommand{\left (}{\left (}
\newcommand{\right )}{\right )}
\begin{equation}
\cos \Theta_\ell = \frac{\epsilon_\ell^{\textrm{CMB}}  \cdot  \left ( \epsilon_\ell^{\textrm{CMB}}  +  \epsilon_\ell^{\textrm{noise}} \right )}{\left | \epsilon_\ell^{\textrm{CMB}} \right |  \left| \epsilon_\ell^{\textrm{CMB}}+ \epsilon_\ell^{\textrm{noise}} \right | }.
\end{equation}
In the case of Gaussian noise, it is possible to obtain an explicit form for the distribution of the angle, generalizing \eqref{htheta}. Introducing the signal to noise $SN$ as the ratio of the norms of the two vectors,
\begin{equation}
SN = \frac{\left | \epsilon_\ell^\textrm{CMB} \right |  }{ \left| \epsilon_\ell^\textrm{noise} \right | } = \sqrt \frac{C_\ell^\textrm{CMB}}{C_\ell^\textrm{noise}},
\end{equation}
one finds (see Appendix~\ref{appendixa} for details)
\begin{equation}\label{pdfsn}
\begin{split}
&h_N\left ( \Theta \right ) = \frac{\Gamma(n)}{\Gamma\left ( \frac{n-1}{2}\right )} \sin^{n-2} \Theta \\
& \cdot  \exp\left (-\frac n2 SN^2 \sin^2\Theta \right )
\textrm{i}^{n-1} \textrm{erfc}\left (- \sqrt{\frac n 2 }  SN\cos \Theta\right ),
\end{split}
\end{equation}
where the special functions
\begin{equation}
\textrm{i}^{n} \textrm{erfc}\left ( z\right ) = \frac 2 {\sqrt \pi} \int_z^\infty dt\: \frac{\left ( t-z \right )^n}{n!} e^{-t^2}
\end{equation}
are the iterated integrals of the complementary error function \citep{abramowitz70a}. These functions satisfy convenient recursion relations allowing easy generation of $h_N(\Theta)$. With the help of $\textrm{i}^{n} \textrm{erfc}\left ( 0 \right ) = 2^{-n}/\Gamma(n/2 + 1)$ we can check that we recover the corresponding distribution \eqref{htheta} for $\Theta$ in the limit of vanishing signal to noise, as expected.

Again, the density function is very close to a Gaussian. Useful simple approximations for its mean and variance are
\begin{equation}
\left \langle{\cos \Theta} \right \rangle \approx \frac{SN}{\sqrt{1 +  SN^2}} = \sqrt{ \frac{C^{\textrm{CMB}}_l}{C_l^{\textrm{CMB}} + C_l^{\textrm{noise}}}}
\label{mean}
\end{equation}
and
\begin{equation}
\left \langle{\cos^2 \Theta} \right \rangle - \left\langle \cos \Theta \right \rangle^2 \approx  \frac{1}{n} \frac{1 + SN^2/2}{\left ( 1 + SN^2 \right )^3}.
\label{var}
\end{equation}
Both of these approximations are already at least $5\%$ accurate for any value of SN at $\ell = 5$.
We evaluate Eqs. \eqref{mean} and \eqref{var} using WMAP Q, V, and W noise realizations, that are white noise to a good approximation, and represent different variances. We compared our model with simulations on Fig. \ref{noises}, and found that higher variance causes decoherence at lower $\ell$. Besides, different realizations of WMAP noise produced almost identical curves, in agreement with our model.

\begin{figure}
\begin{center}
\includegraphics[width=90mm]{noises3.eps}
\caption{We measured $\Theta_{\ell}$ angles between a Gaussian simulation, and the same simulation with WMAP noise added. We show two noise realizations of WMAP's Q, V, and W measurements, and compare them with 2$\sigma$ limits of our noise model. Dashed red line illustrates the expected value of $\pi / 2$ for no correlation, while the solid black curve shows 5$\sigma$ difference from this at each $\ell$.}
\label{noises}
\end{center}
\end{figure}

\section{Results}

\begin{figure}
\begin{center}
\includegraphics[width=90mm]{fig2.eps}
\caption{The measured and theoretically predicted angles between Generalized Phases of Planck Smica map, and WMAP products. A black solid line defines the 5$\sigma$ alignment confidence level, while dashed black curves correspond to 4$\sigma$, 3$\sigma$, 2$\sigma$, and 1$\sigma$ values. The top inset zooms on $\ell<50$, while bottom inset shows the same without foreground cleaning.}
\label{res1}
\end{center}
\end{figure}

We obtained Generalized Phases of WMAP and Planck datasets by applying Equations \eqref{xl} and \eqref{xl2}. We present our results for the Planck Smica map,
but repeating all our analysis with the NILC map produced
virtually identical results.
We used Eq. \eqref{scal} to characterize the coherence of the maps.
While this angle does not contain all information, indeed there are many
ways of constructing a unit vector that is at angle $\Theta$ with respect
to another one, it corresponds to a concise way of expressing coherence,
and we can additionally interpret $\cos\Theta_{\ell}$ in terms of $C_\ell$'s is a
of correlation coefficient  $C^{\rm WMAP, Planck}_\ell/\sqrt{C_\ell^{\rm WMAP}C_\ell^{\rm Planck}}$,
i.e. $60^\circ$  means 50\% correlation between the two maps.

To quantify the coherence, we choose as our null hypothesis that the
two maps are {\em not} correlated. In that case the distribution
$\Theta_\ell$ follows analytic distributions of Eq. \eqref{htheta}, and $p$-values can be calculated by integrating Eq. \eqref{htheta} to the measured $\Theta_\ell$. We define the two maps
as significantly correlated if the null hypothesis can be rejected at the
$5\sigma$ level.

Figure \ref{res1} shows our results, where we compare Planck Smica map to WMAP Q, V and W band measurements. In general, the correlation between the maps
decreases with $\ell$, as qualitatively expected in the presence of uncorrelated
noise. For the lowest $\ell$'s the null hypothesis cannot be rejected at the
$5\sigma$ level, especially for the Q band, but using foreground reduced maps
improves the correlation to the point that maybe only the dipole
is incoherent. This, however, only reflects the different cleaning procedures used by WMAP and Planck. In particular, the Smica algorithm sets $\ell=0,1$ exactly to zero, therefore it contains no information on the CMB (Jean-Francois Cardoso, private communication). The pattern illustrated on Figure \ref{diff} was also detected by  \cite{Fre2013}.

For higher $\ell$'s, the monotonically increasing $p$-values reach the limit confidence levels corresponding to $5\sigma$. We define these $\ell$'s corresponding to decoherence at $\ell\approx700$, $\ell\approx900$ and $\ell\approx1100$ for Q, V and W maps, respectively. This result is robust whether we use foreground removed WMAP maps or not, or Planck Smica/NILC maps. The observed decoherence can be fully explained based on a WMAP noise model, as illustrated
in our Figure~\ref{res1}, and explained in more detail next. Our interpretation is that WMAP GPs are dominated by noise above these $\ell$'s.

Our theory of Eq.~\ref{pdfsn} using simple Gaussian assumption for
both the CMB and noise provides a
prediction for  the expected coherence angle between the
maps. The agreement is excellent with both
simulations and measurements at all $\ell$'s,
although there appears to be small but significant bias
in the measurements at low-intermediate $\ell$'s.
Figure~\ref{bias} displays the residual $\Theta_\ell$, i.e.
the difference between our theoretical predictions for the decoherence
based on our noise model, and the measured angle. For each Q,V and W,
there appears to be an excess angle, i.e. more decoherence than predicted,
for $\ell \lower.5ex\hbox{\ltsima} 500,400$ and
$300$, respectively.

\begin{figure}
\begin{center}
\includegraphics[width=80mm]{fig3.eps}
\caption{Difference map of Planck Smica and WMAP9 Internal Linear Combination (ILC) maps in $\mu K$, smoothed at $2\deg$. See text for details.}
\label{diff}
\end{center}
\end{figure}

At face value in the framework
of our simple assumptions, this would be a sign of
excess noise not taken into account in our noise model. It needs to
be emphasized though that this is a small, (although) significant
effect, and therefore should be interpreted cautiously, given the assumption of uncorrelated
Gaussian noise; noise correlations, foregrounds, and/or leakage from
the dipole (e.g. \cite{PrunetEtal2005,DasSouradeep}) could all influence the
coherence angle in subtle ways.

For completeness, we measured power spectrum of the Planck Smica map, cross-power spectra of WMAP9 Q1-Q2, V1-V2, and an average cross-spectrum of six combinations of W1-W4 differential assemblies with SpICE \citep{SzapudiB}: the power spectrum is complementary to the GPs, corresponds to the
amplitude of the vector we defined in Eq.~\eqref{xl}, and might give
additional insight into the decoherence at low-intermediate $\ell$'s.
We used WMAP9 beam transfer function products for
Q1, Q2, V1, V2, W1, W2, W3, and W4 maps,
and a 5' Gaussian smoothing for the Smica map.
We emphasize that we used again the same resolution maps, with the same mask,
and the same method to measure the power spectrum for all maps, thus our
comparison is more immediate than taking final products from the WMAP
and Planck team, respectively.

The power spectra are consistent with each other for the most part,
but curiously, in approximately the same range of $\ell$'s, where we found
less coherence than predicted by our theory, we find that $\tilde C_{\ell}^{\rm WMAP}$ is on average $2.6\%$ higher than $\tilde C_{\ell}^{\rm Planck}$ in the three Q,V,W maps.
For the sake of consistency, we consider multipoles between $10$ and $300$
for each band, and find that
the WMAP spectra are $2.7\%,2.6\% $ and $2.5\%$ higher than Smica, respectively. While visual inspection confirms the significance of this bias, we estimated it quantitatively in Appendix~\ref{appendixb} to be in the range of $10$'s of
$\sigma$'s. This bias is confined to these scales, the inclusion of higher multipoles result in a non-detection of significant bias.
While it would be difficult to assess quantitatively whether the bias
persists on larger $\ell$'s, at least qualitatively, it appears from
Figure~\ref{bias} that the bias is not significant above the
the same $\ell \lower.5ex\hbox{\gtsima} 500,400$ and $300$ for Q,V, and W, respectively,
where our theory predicts the decoherence based on the simple
Gaussian WMAP noise model. This might be a tantalizing hint, but
more investigations are needed to establish whether the two small,
but significant effects are related.

We repeated our measurements with WMAP 7 year foreground cleaned data, and found similar trends in terms of $\Theta_\ell$ angles. The agreement with WMAP9 results is less accurate, when we analyze maps without cleaning of foregrounds, but the difference is only significant at low $\ell$'s.
The most important observation, however, is that the estimated
$5\sigma$ decoherence is at
slightly lower $\ell$ if we use WMAP 7 year products. This is consistent
with WMAP7 having more noise than WMAP9 further supporting the thesis
that all experiments observe the same underlying CMB, and that instrumental
noise causes the observed decoherence.

\begin{figure}
\begin{center}
\includegraphics[width=90mm]{fig4.eps}
\caption{We show measured biases of power spectra for Q, V and W bands, while estimated $2\sigma$ deviations are shown by solid lines. In addition, discrepancies between modeled and measured $\Theta_{\ell}$ are illustrated for Q, V and W, where dashed lines correspond to $2\sigma$ differences in our model.}
\label{bias}
\end{center}
\end{figure}
\newcommand{\textrm{masked}}{\textrm{masked}}
\newcommand{\textrm{unmasked}}{\textrm{unmasked}}
\newcommand{f_{\textrm{sky}}}{f_{\textrm{sky}}}

\subsection{Decoherence from WMAP noise and impact of mask at high $\ell$}

So far we established that the decoherence observed on Figure \ref{res1} is expected to originate primarily from the noise in WMAP. Assuming that the level of noise in the Smica map is negligible with respect to that of the Q,V,W maps on these scales, we can test this hypothesis using our density functions in Eq.~\eqref{pdfsn}. We proceeded as follows. Assuming white noise $\sigma^2_N$ in each WMAP maps, the signal to noise is $SN = C_{\ell}^\textrm{CMB} / \sigma^2_N$, where $C_{\ell}^\textrm{CMB}$ was generated with the CAMB package \footnote{\texttt{http://camb.info/}} with Planck's best fit parameters \citep{Planck_15}, multiplied by the respective beam window function of the Q,V or W  maps. The solid lines  in Figure \ref{res1} show the mean of the density function, and the dashed ones correspond to $2\sigma$ deviations. The decoherence is in excellent quantitative agreement with this simple model. It makes no difference to use the exact $h_N(\Theta)$ in  Eq.~\eqref{pdfsn} or the approximations in Eqs.~\eqref{mean}~and~\eqref{var}.

The noise dominates by orders of magnitude at the highest $\ell$'s, therefore
an angle of $90$ degrees is expected naively. The observed angles, however,
deviate slightly from this theoretical prejudice,
indicating a few percent residual correlation. As we show
in more detail in Appendix~\ref{appendixc},
this correlation is due to leakage of low $\ell$ power
into higher $\ell$'s, and essentially white noise.
We can obtain accurate analytic approximations assuming
an azimuthally symmetric mask centered on the equator and white noise. The mask is an equatorial band sustaining an angle $b$ with the equator, so that $f_{\textrm{sky}} = 1-\sin b$. Using the explicit formula relating the spectrum $\tilde C_{\ell}$ of the masked field to that of the unmasked field $C_{\ell}$ \citep{2002ApJ...567....2H}, we derive in Appendix~\ref{appendixc} the asymptotic behavior of the spectrum,
\begin{equation}
\tilde C_\ell \rightarrow \quad \frac{16 \sigma^2_T \cos b}{\left ( 2\ell +1 \right )^3},\quad \ell \rightarrow \infty
\end{equation}
where $\sigma^2_T$ is the variance of the unmasked map
\begin{equation}
\sigma^2_T = \sum_{\ell} \frac{2\ell+1}{4\pi}C_{\ell} = \left \langle{ \left ( \frac{\Delta T}{\bar T}\right )^2} \right\rangle.
\end{equation}
On the other hand the white noise spectra are simply multiplied by $f_{\textrm{sky}}$. Since $\left\langle \cos \Theta_l \right\rangle \approx \sqrt{\tilde C_l / (\tilde C_l + \sigma^2_N)}$, we obtain in the very low signal to noise regime
\begin{equation}
\begin{split}
\cos \Theta_{\ell}  \rightarrow   \frac{4}{\left ( 2\ell + 1 \right )^{3/2}}  \left ( \frac{\sigma_T}{\sigma_N} \right ) \left ( \frac{2-
f_{\textrm{sky}}}{f_{\textrm{sky}}} \right )^{1/4}
\end{split}
\end{equation}
Despite the above approximations, these ideas explain the shape of measured $\Theta_{\ell}$ curves extremely well, and
predict asymptotic properties at high $\ell$ in virtually perfect agreement with simulations and measurements. Note that these considerations do not affect
our 5$\sigma$ decoherence limits, as our null hypothesis of no correlations
(corresponding to infinite noise) has no bias.

We used our well calibrated decoherence model to forecast GP angles of Planck and a hypothetical perfect CMB experiment without noise (Fig. \ref{estim}). Decoherence is predicted at $\ell\approx 2900$, beyond which any non-Gaussian information should be dominated by noise.

\begin{figure}
\begin{center}
\includegraphics[width=90mm]{fig5.eps}
\caption{Estimated decoherence properties of a Planck-like simulation ($N_{side}=2048$, assuming Smica noise) with a hypothetical perfect CMB experiment with no noise, and our WMAP results for comparison. We show simulated decoherence using single realizations, together with theory and $2\sigma$ errors. Black solid line illustrates $5\sigma$ significance level for our null hypothesis, while dashed red line shows $90^{\circ}$ decoherence level.}
\label{estim}
\end{center}
\end{figure}

\section{Conclusions}

We quantified the $\ell$-by-$\ell$ coherence of
latest  WMAP, and Planck CMB maps. We introduced a new set of statistics, Generalized Phases, that are complementary to the (pseudo-)power spectrum, and can
be used to characterize the phase-coherence of two CMB maps.
We compared GP's of the two maps by simply calculating the
angles between the corresponding unit vectors. These angles, while
do not contain all non-Gaussian information, concisely summarize the
coherence properties of two maps at each $\ell$.
Using the statistics of random vectors in $(2\ell+1)$ dimensions,
we defined the $\ell$ of decoherence where the null hypothesis of
no correlation between the maps could not be rejected at the $5\sigma$
level. We controlled any effect of the masks, typically a problem
with statistics based on phases, with careful simulations and analytical
models that, albeit based on simplifying assumptions, appear to provide
an excellent quantitative framework. To check for systematics, we
repeated all our measurements of the Planck Smica map with the NILC maps
finding virtually identical results. According to our definition,
decoherence from Planck was found above $\ell\approx 700$, $\ell\approx 900$ and $\ell\approx 1100$ for WMAP9 Q, V and W.
Our theoretical description is in excellent agreement with the measured
coherence angles, with a slight bias for low-intermediate $\ell$'s
We also find
a small bias of the WMAP pseudo-$\tilde C_{\ell}$ at  $10 \le \ell \le 300$
at an average $2.6$\% level with very high significance.
It appears that for high $\ell$'s, where our theoretical prediction
for the coherence angle is accurate based on a simple Gaussian WMAP noise
model, there is no significant bias in the power spectra either.
Qualitatively, there is a slight color dependency as well based on
Figure~\ref{bias}. From the excess decoherence we can calculate the
amount of excess noise it corresponds to. We found that,
with the exception of the Q map in the range of
$250 \lower.5ex\hbox{\ltsima} \ell \lower.5ex\hbox{\ltsima}  500$ the noise
corresponding to the excess decoherence is below what is needed
to fully explain the bias in the power spectra. Nevertheless,
the qualitative behaviour of the noise is similar to the observed
one, and it is different than our simulations.
In conclusion, there are tantalizing coincidences hinting
that the excess decoherence and power spectrum bias are related,
but no consistent picture emerged. Note that our simulations
do not contain correlated noise, we did not check for any effect
of foregrounds or low-$\ell$ leakage, especially from the dipole,
into higher $\ell$'s; such investigations
are left for future research.

Our analytical and simulation framework can be used to
forecast the coherence of Planck with a noise-free experiment (the true
CMB). We find that below $\ell\le 2900$ Planck is coherent with the
CMB according to our $5\sigma$ criterion,
thus non-Gaussian information can be best gleaned from below these
$\ell$'s.

\section*{Acknowledgments}
We acknowledge support from NASA grants NNX12AF83G and NNX10AD53G. In addition, AK acknowledges support from Campus Hungary fellowship program, and OTKA through grant no. 101666.
\bibliographystyle{mn2e}