\section{Introduction}
Although the mathematical study of the golden ratio dates to antiquity, one of its unusual properties appears to have passed the millennia unnoticed: its minimal polynomial $f(x) = x^2 - x - 1$ is irreducible over $\Q$, as is the second iterate $f(f(x))$ of this polynomial, but
\begin{equation} \label{goldenfact}
f(f(f(x)))
= (x^4 - 3x^3 + 4x - 1)(x^4 - x^3 - 3x^2 + x + 1),
\end{equation}
where the two quartics are irreducible. The factorization in \eqref{goldenfact} together with the irreducibility over $\Q$ of $f(f(x))$ is a rare phenomenon: among polynomials $x^2 + ax + b$ with $a, b$ integers satisfying $|a| \leq 100,000$ and $|b| \leq 1,000,000,000$, only 8 others have this property. One might wonder whether there are infinitely many such polynomials. Indeed there are, as we show
in Theorem \ref{summary} below.

More generally, let $K$ be a field and $f \in K[x]$ a polynomial of degree $d \geq 2$. Write $f^n(x)$ for the $n$th iterate of $f$. We say that $f$ has a \textit{newly reducible $n$th iterate} (over $K$) for some $n \geq 2$ if $f^{n-1}$ is irreducible over $K$ but $f^n$ is not irreducible over $K$. We note that the irreducibility of $f^{n-1}$ implies that each of $f, f^2, \ldots, f^{n-2}$ is also irreducible over $K$. Several recent papers have given conditions ensuring that all iterates of a polynomial remain irreducible (see e.g. \cite{ayad}, \cite{ayadcor}, \cite{danielson}, \cite{goksel}, \cite{quaddiv}, \cite{itconst}). Many fewer have studied newly reducible iterates; a few examples are \cite{HCsummer}, \cite{fein}, and \cite{preszler}. In this paper our main objects of study are the following.

\begin{definition}
Let $d,n$ be integers that are at least two. Define
$\mathcal{N}_{d,n}$ (resp. $\mathcal{N}_{d,n}^{\infty})$ to be the class\footnote{Because the collection of all fields is not a set, neither are $\mathcal{N}_{d,n}$ and $\mathcal{N}_{d,n}^\infty$. We use them in this article not as objects, but as notational devices. In referring to them we use the usual notation and language of set theory.} of all fields $K$ such that there exists at least one (resp. infinitely many) $f \in K[x]$ of degree $d$ with newly reducible $n$th iterate over $K$.
\end{definition}

When $K$ is a perfect field, the question of whether $K \in \mathcal{N}_{d,n}$ can be phrased as an inverse Galois-type problem: does there exist $f \in K[x]$ such that the absolute Galois group of $K$ acts transitively on the roots of $f^{n-1}$ but not on the roots of $f^n$? Fein and Schacher \cite{fein} appear to be the first to have studied newly reducible iterates, and they used results of Odoni to obtain the fundamental result in this area: for all $d,n \geq 2$, $\mathcal{N}_{d,n}$ is non-empty (\cite[Corollary 1.3]{fein}). Their method relies on knowledge of the Galois groups of iterates of generic polynomials, and then an appeal to the Hilbert irreducibility theorem; the resulting field $K$ depends on the choice of specialization. Our aim is to address at least some cases of the following questions, each of which is unresolved by the methods of \cite{fein}.

\begin{question} \label{mainquest}
Let $d,n$ be integers that are at least $2$.
\begin{enumerate}
\item Precisely which fields belong to $\mathcal{N}_{d,n}$ and $\mathcal{N}_{d,n}^\infty$?
\item For which $d,n$ is $\Q \in \mathcal{N}_{d,n}$? For which $d,n$ is $\Q \in \mathcal{N}_{d,n}^\infty$?
\item If $K$ is a number field with ring of integers $\mathcal{O}_K$ and $K \in \mathcal{N}_{d,n}^\infty$, do there exist infinitely many monic $f \in \mathcal{O}_K[x]$ of degree $d$ with newly reducible $n$th iterate over $K$?
\end{enumerate}
\end{question}

The case $K = \mathbb{Q}$ has received significant attention recently. The recent paper \cite{goksel} of Goksel describes an infinite family of monic quadratic $f \in \Z[x]$ with newly reducible third iterate \cite[Lemma 3.9]{goksel}, thereby showing that $\Q \in \mathcal{N}_{2,3}$.
The paper \cite{preszler} gives an infinite family of non-monic cubic $f \in \Z[x]$ with newly reducible second iterate, thus proving $\Q \in \mathcal{N}_{3,2}$,
and conjectures that no such family exists among monic cubic $f \in \Z[x]$ \cite[Conjecture 4.1]{preszler}.

We summarize our main results. Throughout the article, we sometimes use $\text{char}(K)$ to denote the characteristic of a field $K$.

\begin{theorem} \label{summary}
Let $K$ be an infinite field with $\text{char}(K) \neq 2$, and for $n \geq 2$ set $K^n = \{k^n : k \in K\}$.
\begin{enumerate}
\item We have $K \in \mathcal{N}_{2,2}$ if and only if $K \neq K^{2}$, and $K \in \mathcal{N}_{2,2}^\infty$ if and only if $K \setminus K^{2}$ is infinite (Proposition \ref{n22}).
\item We have $K \in \mathcal{N}_{2,3}^\infty$ if
\begin{itemize}
\item $-1 \not\in K^2$ and $K$ has a discrete valuation $v$ with $v(5)$ odd (Proposition \ref{m-1}); or
\item $K$ is a totally real number field and there is a prime of residue degree 1 lying over $(3)$ (Corollary \ref{newfamily}).
\end{itemize}
\item We have $K \in \mathcal{N}_{3,2}^\infty$ if $\text{char}(K) \neq 3$ and $2 \not\in K^3$ (Theorem \ref{cubic2}).
\item We have $K \in \mathcal{N}_{4,2}^\infty$ if $\text{char}(K) \neq 3$, $3 \not\in K^2$, and $-3 \not\in K^4$ (Corollary \ref{deg-4-infinitely-many-corollary}).
\item We have $K \in \mathcal{N}_{d,2}^\infty$ for all $d \equiv 2 \pmod{4}$ if $-1 \not\in K^2$ and $K$ has a non-trivial discrete valuation (Corollary \ref{genbigd}).
\end{enumerate}
\end{theorem}
In each case, we exhibit an explicit infinite family of polynomials with the desired properties, and if $K$ is a number field the polynomials may be taken to be monic with coefficients in $\mathcal{O}_K$ without loss of generality. For instance, a consequence of Theorem \ref{cubic2} is that if
\begin{equation} \label{cubicfamily}
f(x) = \left(x + 93312 t^9+36 t^3\right)^3-93312 t^9
\end{equation}
for $t \in \Z \setminus \{0\}$, then $f$ has a newly reducible second iterate over $\Q$. This disproves Conjecture 4.1 of \cite{preszler}, which is based on a search of all polynomials of the form $x^3 + ax^2 + bx + c,$ where $a,b,c$ are integers with absolute value at most 500; it is noteworthy that none of the polynomials in \eqref{cubicfamily} are of this form.

When $K$ is a finite field of characteristic not equal to $2$, all the statements of Theorem \ref{summary} hold with $\mathcal{N}_{d,n}^\infty$ replaced by $\mathcal{N}_{d,n}$, although (5) and the first statement of (2) become vacuous.  We have $\mathbb{F}_2 \not\in \mathcal{N}_{2,2}$, since $f(x) = x^2 + x + 1$ is the only irreducible quadratic over $\mathbb{F}_2$, and $f^2(x)$ is irreducible over $\mathbb{F}_2$. We prove in Proposition \ref{fintwo} that $\mathbb{F}_{2^n} \in \mathcal{N}_{2,2}$ for all $n > 1$ and $\mathbb{F}_{2^n} \in \mathcal{N}_{2,3}$ for all $n \geq 1$. Surprisingly, $\mathbb{F}_{2^n} \not\in \mathcal{N}_{2,k}$ for any $k \geq 4$ and any $n \geq 1$, due to a result of Ahmadi et al \cite[Theorem 10]{shparostafe} (see also \cite{ahmadi}) that every $f \in \mathbb{F}_{2^n}[x]$ has $f^3(x)$ reducible.

When $K$ is a general field of characteristic 2, we know of no results addressing whether $K \in \mathcal{N}_{d,n}$ for any $d,n$, though we give a condition that applies to the case $d = n = 2$ in Proposition \ref{secondchar2}. This leads us to pose the following question.

\begin{question} \label{charpquest}
Let $p$ be a prime number and $d \geq p$. Which fields of characteristic $p$ belong to $\mathcal{N}_{d,n}$ for various $n \geq 2$?
\end{question}

The results of Theorem \ref{summary} prompt the following two conjectures.
\begin{conjecture}
For each $d, n \geq 2$, $\mathcal{N}_{d,n}^\infty$ is non-empty.
\end{conjecture}

\begin{conjecture}
Let $K$ be a finite extension of $\Q$ or $\Fp(t)$. Then $K \in \mathcal{N}_{2,3}^\infty$ and $K \in \mathcal{N}_{d,2}^\infty$ for every $d \geq 2$.
\end{conjecture}

The proofs of the statements in Theorem \ref{summary} proceed by enumerating all $f$ of a specified form such that $f^n(x)$ is reducible, and then giving conditions under which infinitely many $f$ remain after discarding those with $f^m$ reducible for $m<n$. In the case where $f$ is quadratic, the specified form is $f(x) = (x-\gamma)^2 + \gamma + m$, so that $\gamma$ is the critical point of $f$, and this includes all $f$. We thus obtain the following result.

\begin{theorem} \label{characterization}
Let $K$ be a field of characteristic not equal to $2$, and let $f\left(x\right)=\left(x-\gamma\right)^2+\gamma+m\in K[x]$. Then $f$ has a newly reducible third iterate over $K$ if and only if the following both hold:
\begin{enumerate}
\item There exist $r,s\in K$ such that
\begin{align*}
\gamma & = \frac{1}{256r^2}(-2 r^5 s^2+9 r^4 s^4-4 r^4 s^2-16 r^3 s^6+32 r^3 s^4+16 r^3 s^2+32 r^3+14 r^2 s^8 \\ &\phantom{=\frac{1}{256r^2}(}-64 r^2 s^6-8 r^2 s^4+96 r^2 s^2+128 r^2-6 r s^{10}+48 r s^8-64 r s^6 \nonumber \\ &\phantom{=\frac{1}{256r^2}(} -160 r s^4+96 r s^2+128 r+s^{12}-12 s^{10}+40 s^8-112 s^4-64 s^2), \nonumber \\
m & = \frac{-r^2-2 r s^2-4 r+s^4-4 s^2-4}{8 r}.
\end{align*}
\item None of $\sqrt{-m-\gamma}$, $\sqrt{-2m + 2\sqrt{m^2+m+\gamma}}$, $\sqrt{-2m - 2\sqrt{m^2+m+\gamma}}$ is in $K$.
\end{enumerate}
\end{theorem}
We note that taking $K = \Q$ and $r = s = 1$ in Theorem \ref{characterization} gives $\gamma = 1/2$ and $m = -7/4$. This yields $f(x) = (x-1/2)^2 + 1/2 - 7/4 = x^2 - x - 1$, and thus we recover the minimal polynomial of the golden ratio, mentioned in the first paragraph. Observe that for this polynomial, $-m - \gamma = 5/4$ and $m^2 + m + \gamma = 29/16$, ensuring that condition (2) in Theorem \ref{characterization} holds.

In general, a convenient way to verify the conditions in part (2) of Theorem \ref{characterization} is to show that neither $-m - \gamma$ nor $m^2 + m + \gamma$ is a square in $K$. With $\gamma, m$ as in part (1) of Theorem \ref{characterization}, up to squares in $K(r,s)$ we have $-m-\gamma = 2r - s^2 + 4$ and
$$
m^2 +m + \gamma = -2s^2r^3 + (5s^4 + 4s^2 + 4)r^2 + (-4s^6 + 12s^4 + 32s^2 + 16)r + s^8 - 8s^6 + 8s^4 + 32s^2 + 16
$$

In \cite[Lemma 3.9]{goksel}, Goksel studies the case $m = -1$ and gives algebraic characterizations of all $\gamma$ making $f^2(x)$ and $f^3(x)$ reducible. We recover his results as an outcome of our proof of Theorem \ref{characterization} (see the discussion following the proof of Theorem \ref{symmchar}). We use Theorem \ref{characterization} to give another family with non-constant $m$ in Theorem \ref{newfamily}.

We close this introduction with three additional questions, and an outline of the paper.

\begin{question} \label{fourquest}
Is $\Q \in \mathcal{N}_{2,4}$?
\end{question}

We have not been able to find any quadratic $f \in \Q[x]$ with newly reducible fourth iterate; we briefly discuss some work on this question in Section \ref{fourth}. We remark that  \cite[Theorem 1.2]{goksel} shows that for $m \in \{0,-1,-2\}$ (i.e. when $f$ is post-critically finite), there is no $\gamma \in \Z$ that works.

\begin{question} \label{cubicquest}
Does there exist $f \in \Q[x]$ such that $\deg f = 3$ and $f^2$ is newly reducible with three distinct factors, each of degree $3$?
\end{question}

In the course of our study of cubic polynomials in Section \ref{cubics} we find infinite families with newly reducible second iterate, but all such families with $f \in \Q[x]$ have $f^2(x)$ that factors as a product of an irreducible sextic and cubic. It would be interesting to find an $f \in \Q[x]$ with second iterate factoring as in Question \ref{cubicquest}. More generally, it would be interesting to study polynomials with newly reducible iterates that have more than two irreducible factors.

To motivate our last question, we remark that the fields in $\mathcal{N}_{d,n}$ constructed in \cite{fein} are number fields whose degrees grow very rapidly with $n$ and $d$.

\begin{question}
Fix $d \geq 2$, and let $m_n$ be the minimal degree of a number field $K$ over which there is $f \in K[x]$ of degree $d$ with newly reducible $n$th iterate. Is $(m_n)_{n \geq 2}$ unbounded? Does $\lim_{n \to \infty} m_n = \infty$?
\end{question}

In Section \ref{symmetric}, we prove under mild hypotheses that if $f$ has a newly reducible $n$th iterate, then the factorization of $f^n$ must have a certain form. This generalizes some of the results from \cite{settled}. These results are particularly useful in the case where $f$ is quadratic, which we study in Section \ref{surface}, proving Theorem \ref{characterization} and parts (1) and (2) of Theorem \ref{summary}, as well as some results for fields of characteristic 2, such as Proposition \ref{fintwo}. Section \ref{fourth} contains brief remarks on Question \ref{fourquest}. In Section \ref{cubics}, we study cubics with newly reducible second iterate, and prove part (3) of Theorem \ref{summary}. Sections \ref{quartics} and \ref{largerdegrees} study polynomials of higher degree with newly reducible second iterate, and contain the proofs of parts (4) and (5) of Theorem \ref{summary}.

\medskip

\noindent \textbf{Acknowledgements.} We thank the anonymous referee for helpful comments. We are also grateful to Carleton College's Towsley Endowment for the Sciences, which partially supported the research of the third author.

\section{The form of factors} \label{symmetric}

We begin with some very general results on factorization of iterates, culminating in Theorem \ref{symmetric thm}, which gives a generalization of \cite[Proposition 2.6]{settled} to arbitrary characteristic.
The key lemma, Lemma \ref{surjective}, can be proven using Capelli's Lemma, as in \cite[Proposition 2.6]{settled}, but we give here a self-contained proof.

Fix an algebraic closure $\overline{K}$ and separable closure $K^{\text{sep}}$ of $K$. Recall that $f \in K[x]$ is separable over $K$ if it has $\deg f$ distinct roots in $\overline{K}$, or equivalently, all roots of $f$ in $\overline{K}$ lie in $K^{\text{sep}}$.

\begin{definition}
Let $K$ be a field and $f \in K[x]$ have degree $2n$ for some $n \geq 1$. We say that $f$ is \textit{symmetrically reducible} over $K$ if there exists a monic $g \in K[x]$ with $f(x) = Cg(x)g(-x)$ for some $C \in K$.
\end{definition}

\begin{lemma} \label{deglem}
Let $K$ be a field, and let $f(x) \in K[x]$ have degree $d \geq 2$. If $f^n(x)$ is newly reducible for some $n \geq 2$, then each irreducible factor of $f^n(x)$ has degree divisible by $d^{n-1}$.
\end{lemma}

\begin{proof}
Let $\alpha \in \overline{K}$ be a root of $f^n(x)$. Then $f(\alpha) \in K(\alpha)$, and so $K(f(\alpha)) \subseteq K(\alpha)$. Hence
\begin{equation} \label{tower}
[K(\alpha) : K] = [K(\alpha) : K(f(\alpha))][K(f(\alpha)) : K] = d^{n-1}[K(\alpha) : K(f(\alpha))],
\end{equation}
where the last equality follows because $f(\alpha)$ is a root of $f^{n-1}(x)$ and $f^{n-1}$ is assumed to be irreducible over $K$. Now let $h$ be any irreducible factor of $f^n$, and $\alpha$ a root of $h$. Then $\deg h = [K(\alpha) : K]$ is divisible by $d^{n-1}$ from \eqref{tower}.
\end{proof}

\begin{lemma} \label{surjective}
Let $K$ be a field and let $g,f \in K[x]$ be separable over $K$. Assume $g$ is irreducible over $K$, and let $h \in K[x]$ be any non-constant factor of $g \circ f$. Consider the map
\begin{equation}
\Phi : \{\text{roots of $h$ in $K^{\text{sep}}$}\} \to \{\text{roots of $g$ in $K^{\text{sep}}$}\} \label{fmap}
\end{equation}
defined by $\Phi(\alpha) = f(\alpha)$. Then $\Phi$ is a surjective $k$-to-1 map for some $k$ with $1 \leq k \leq \deg f$.
\end{lemma}

\begin{proof}
First note that a root $\alpha$ of $h$ satisfies $g(f(\alpha)) = 0$, and so $f(\alpha)$ is indeed a root of $g$, showing that $\Phi$ is well-defined. Because $h$ is non-constant, the image of $\Phi$ is non-empty.

Let $G_K := \Gal(K^{\text{sep}}/K)$ be the absolute Galois group of $K$. The key observation is that since $f$ is defined over $K$, it must commute with the action of $G_K$ on $K^{\text{sep}}$. Hence if $\alpha'$ is a root of $g$ and $f(\alpha) = \alpha'$, then for any $\sigma \in G_K$ we have
\begin{equation} \label{galois}
\sigma(\alpha') = \sigma(f(\alpha)) = f(\sigma(\alpha)).
\end{equation}
Now let $\alpha', \beta'$ be roots of $g$. Because $g$ is irreducible over $K$, $G_K$ acts transitively on the roots of $g$, and hence there is some $\sigma \in G_K$ with $\sigma(\alpha') = \beta'$. From \eqref{galois} we have $\beta' = f(\sigma(\alpha))$. If we further assume that $\alpha$ is a root of $h$, then so must be $\sigma(\alpha)$, since $h$ is defined over $K$ and hence the set of its roots is preserved by $G_K$. Therefore $\sigma$ induces a map $\Phi^{-1}(\alpha') \to \Phi^{-1}(\beta')$. This map is injective since $\sigma$ is an injection from the set of all roots of $g \circ f$ into itself. Similarly, $\sigma^{-1}$ gives an injection $\Phi^{-1}(\beta') \to \Phi^{-1}(\alpha')$. It follows that $\#\Phi^{-1}(\alpha') = \#\Phi^{-1}(\beta')$, and thus all fibers of $\Phi$ have equal cardinality. But $\Phi$ has non-empty image, and the lemma is proved.
\end{proof}

\begin{lemma} \label{seplem}
Let $K$ be a field and $f \in K[x]$. If $f$ is irreducible over $K$ and $f'$ is not identically zero, then all iterates of $f$ are separable over $K$
\end{lemma}

\begin{proof}
Fix an algebraic closure $\overline{K}$ of $K$, and recall the well-known fact that $g \in K[x]$ is separable over $K$ if and only if $g$ and $g'$ have no common roots in $\overline{K}$. We proceed by induction on $n$, the number of iterations of $f$. Let $n = 1$, and suppose there is $\alpha \in \overline{K}$ with $f(\alpha) = f'(\alpha) = 0$. Because $f$ is irreducible over $K$, we have $[K(\alpha) : K] = \deg f$. But $\deg f' < \deg f$, and so $f'(\alpha) = 0$ forces $f'$ to be identically zero, which is a contradiction.

Assume now that $f^{n}$ is separable over $K$ for some $n \geq 1$. The chain rule gives
\begin{equation} \label{chain}
(f^{n+1})'(x) = (f^{n})'(f(x)) \cdot f'(x).
\end{equation}
Suppose that $\alpha \in \overline{K}$ satisfies $f^{n+1}(\alpha) = (f^{n+1})'(\alpha) = 0$. From \eqref{chain}, we have either $(f^{n})'(f(\alpha)) = 0$ or $f'(\alpha) = 0$. In the former case, we also have $0 = f^{n+1}(\alpha) = f^{n}(f(\alpha))$, contradicting the separability of $f^n$. In the latter case, we also have $0 = f^{n+1}(\alpha) = f(f^n(\alpha))$, implying that $f^n(\alpha)$ is a root of $f$ with $[K(f^n(\alpha)) : K] \leq [K(\alpha) : K] \leq \deg f' < \deg f$. This contradicts the irreducibility of $f$.
\end{proof}

\begin{theorem} \label{symmetric thm}
Let $K$ be a field, let $f(x) = ax^2 + bx + c \in K[x]$ with $a \neq 0$, and assume that $f'$ is non-zero. If $f^n(x)$ is newly reducible over $K$ for some $n \geq 2$, then there is a monic irreducible $h \in K[x]$ of degree $2^{n-1}$ such that
$$f^n(x) = a^{2^n-1}h(x)h(-(x+(b/a))).$$
In particular, if $K$ has characteristic not equal to $2$, then $f^n(x + \gamma)$ is symmetrically reducible over $K$, where $\gamma = -b/2a$ is the critical point of $f$.
\end{theorem}

\begin{proof}
Assume that $f^n$ is newly reducible over $K$ for some $n \geq 2$, and observe that by Lemma \ref{deglem} we have $f^n(x) = \ell(f^n)h_1(x)h_2(x)$, where $\ell(f^n)$ is the leading coefficient of $f^n$, and $h_1, h_2 \in K[x]$ are monic irreducibles of degree $2^{n-1}$. An easy induction shows that $\ell(f^n) = a^{2^n-1}$. By Lemma \ref{seplem}, both $f^{n-1}$ and $f^n$ are separable over $K$, and it follows that $h_1$ and $h_2$ are separable over $K$ as well. By Lemma \ref{surjective}, we have that for $i=1,2$ the maps
  \begin{equation*}
    \Phi_i : \{\text{roots of $h_i$ in $K^{\text{sep}}$}\} \to \{\text{roots of $f^{n-1}$ in $K^{\text{sep}}$}\}
  \end{equation*}
  defined by $\Phi_i(\alpha)=f(\alpha)$, are bijections. Therefore $\Phi_2^{-1} \circ \Phi_1$ is a bijection from the roots of $h_1$ to the roots of $h_2$. Letting $\alpha$ be a root of $h_1$, we have that $\Phi_2^{-1}(\Phi_1(\alpha))$ is a root $\beta$ of $h_2$ satisfying $f(\alpha) = f(\beta)$. This gives
  $a\alpha^2 + b \alpha = a \beta^2 + b \beta$, and because $\alpha \neq \beta$ we further deduce $\beta = - \alpha - \frac{b}{a}$.

  Now let $\alpha_1, \ldots, \alpha_{2^{n-1}}$ be the roots of $h_1$ in $K^{\text{sep}}$ and $\beta_1, \ldots, \beta_{2^{n-1}}$ be the roots of $h_2$ in $K^{\text{sep}}$. We have
  $$
  h_2(x) = \prod_{i=1}^{2^{n-1}} (x - \beta_i) = \prod_{i=1}^{2^{n-1}} (x - (-\alpha_i - (b/a))) = \prod_{i=1}^{2^{n-1}} ((x + (b/a)) + \alpha_i).
  $$
  Because $n \geq 1$, this last expression is the same as $\prod_{i=1}^{2^{n-1}} (-(x + \frac{b}{a}) - \alpha_i)$, which is $h_1(-(x+\frac{b}{a}))$.

  If $K$ has characteristic not equal to $2$, then from $f^n(x) = \ell(f^n)h_1(x)h_1(-(x+\frac{b}{a}))$ we have $f^n(x - \frac{b}{2a}) = \ell(f^n)h_1(x-\frac{b}{2a})h_1(-(x+\frac{b}{2a})),$ whence
  $f^n(x - \frac{b}{2a}) = \ell(f^n)g(x)g(-x)$ for $g(x) = h_1(x-\frac{b}{2a})$. Thus $f^n(x - \frac{b}{2a})$ is symmetrically reducible over $K$.
  \end{proof}

\section{Quadratic polynomials with newly reducible second and third iterate} \label{surface}

A monic quadratic polynomial in $K[x]$ has the form $x^2+ax+b$, but when $K$ has characteristic not equal to 2 we may write it as $(x-\gamma)^2+m+\gamma$ with $\gamma = -a/2$ and $m = b + a/2 - a^2/4$. This latter form emphasizes the critical point $\gamma$ of $f$, and behaves more simply with respect to iteration.
In this section, we prove Theorem \ref{characterization} and give some results in the case where $K$ has characteristic 2.

To prove Theorem \ref{characterization}, it is enough, in light of Theorem \ref{symmetric thm}, to characterize quadratic polynomials with symmetrically reducible third iterate, and then discard those with a reducible first or second iterate. Clearly $f(x)$ is reducible over $K$ if and only if $\sqrt{-m-\gamma}$ is in $K$.

We begin by giving a criterion for when $f^2(x+\gamma)$ has symmetrically reducible second iterate.

\begin{lemma}\label{second-iterate-criterion}
Let $K$ be a field of characteristic not equal to $2$, and let $f\left(x\right)=\left(x-\gamma\right)^2+\gamma+m$ for $\gamma, m\in K$. Then $f^2\left(x+\gamma\right)$ is symmetrically reducible over $K$ if and only if at least one of
\begin{align*}
\sqrt{-2m\pm2\sqrt{m^2+m+\gamma}
\end{align*}
is in $K$.
\end{lemma}

\begin{proof}
We have
$$f^2\left(x\right) = \left(x-\gamma\right)^4 +2 m \left(x-\gamma \right)^2 + m^2 + m + \gamma.$$
On the other hand, by definition $f^2\left(x+\gamma\right)$ is symmetrically reducible over $K$ if and only if $f^2\left(x+\gamma\right) = (x^2 + cx + d)(x^2 - cx + d)$ for some $c,d \in K$.
Equating coefficients gives
\begin{align}
m^2+m+\gamma &= d^2 \label{one}  \\
2m &= 2d-c^2 \label{two}
\end{align}
Clearly \eqref{one} has a solution $d \in K$ if and only if $\sqrt{m^2+m+\gamma} \in K$. Substituting $d = \pm \sqrt{m^2+m+\gamma}$ into \eqref{two}, we see that \begin{align*}
c &= \pm \sqrt{-2m\pm 2\sqrt{m^2+m+\gamma}}.
\end{align*}
If one of these choices of $c$ lies in $K$, then so does the corresponding choice of $d$, and hence $f^2\left(x+\gamma\right)$ is symmetrically reducible over $K$.
On the other hand, if neither of $\sqrt{-2m\pm2\sqrt{m^2+m+\gamma}}$ is in $K$, then there is no $c \in K$ satisfying \eqref{two}, and hence $f^2(x + \gamma)$ is not symmetrically reducible over $K$.
\end{proof}

\begin{proposition} \label{n22}
Let $K$ be a field of characteristic not equal to $2$, and let $K^2 = \{k^2 : k \in K\}$. Then $K \in \mathcal{N}_{2,2}$ (resp. $K \in \mathcal{N}_{2,2}^\infty$) if and only if $K \neq K^{2}$ (resp. $K \setminus K^{2}$ is infinite).
\end{proposition}

\begin{proof}
Equations \eqref{one} and \eqref{two} imply that for any $a,b \in K$ we can construct $f$ with symmetrically reducible second iterate, simply by taking $m = a - \frac{1}{2}b^2$ and $\gamma = a^2 - m^2 - m$. Taking $b = 2$ yields the family $f(x) = (x - (3a-2))^2 + 4a - 4$ for $a \in K$, for which we have
\begin{equation*}
f^2(x) = \left(x^2 + (-6a + 2)x + 9a^2 - 5a \right)\left(x^2 + (-6a + 6)x + 9a^2 - 17a + 8\right).
\end{equation*}
Thus $f$ has newly reducible second iterate provided $f(x)$ is irreducible over $K$, and because $K$ has characteristic not equal to $2$, this is equivalent to $1-a \not\in K^2$. Because $x \mapsto 1-x$ is a bijection on $K$, it follows that $K \in \mathcal{N}_{2,2}$ (resp. $K \in \mathcal{N}_{2,2}^\infty$) if $K \neq K^{2}$ (resp. $K \setminus K^{2}$ is infinite). The converse follows from the observation that $K = K^2$ implies that no quadratic polynomial over $K$ is irreducible, and $K \setminus K^{2}$ finite implies that only finitely many quadratic polynomials over $K$ are irreducible.
\end{proof}

When $K$ has characteristic 2, the corresponding result to Lemma \ref{second-iterate-criterion} is more complicated.

\begin{lemma} \label{secondchar2}
Let $K$ be a field of characteristic 2, and let $f(x) = x^2 + ax + b$. Then $f^2(x) = h(x)h(x+a)$ for some $h \in K[x]$ if and only if $b = 0$, $b = a+1$, or
\begin{equation} \label{poly}
x^4 + a^2x^3 + a^3x^2 + a^2(b^2 + ab + b)x + (b^2 + ab + b)^2
\end{equation}
has a root in $K$.
\end{lemma}

\begin{proof}
We have $f^2(x) = x^4 + (a^2 + a)x^2 + a^2x + (b^2 + ab + b)$. Letting $h(x) = x^2 + cx + d$ and equating coefficients in $f^2(x) = h(x)h(x+a)$ gives
\begin{align*}
a^2 + ac + c^2 & = a^2 + a \\
a^2c + ac^2 &= a^2 \\
a^2d + acd + d^2 &= b^2 + ab + b
\end{align*}
The first equation gives $ac + c^2 + a = 0$, which renders the second equation redundant. Multiplying through by $a^2d^2$ gives $a^2d(acd) + (acd)^2 + a^3d^2 = 0$. From the third equation we have
$acd = a^2d + d^2 + b^2 + ab + b$, and substitution yields
$$a^2d(a^2d + d^2 + b^2 + ab + b) + (a^2d + d^2 + b^2 + ab + b)^2 + a^3d^2 = 0.$$
Expanding, simplifying, and writing this as a polynomial in $d$ yields
$$d^4 + a^2d^3 + a^3d^2 + (a^3b + a^2b^2 + a^2b)d + a^2b^2 + b^4 + b^2 = 0,$$
which is equivalent to $x^4 + a^2x^3 + a^3x^2 + a^2(b^2 + ab + b)x + (b^2 + ab + b)^2$ having a root in $K$. When such a root exists, we may solve for $c$ provided that $a \neq 0$ and $d \neq 0$. But $a=0$ forces $c = 0$, and so we obtain $c$ anyway. If $d = 0$, then we get $b^2 + ab + b = 0$, and so $b = 0$ or $b = a+1$.
\end{proof}

\begin{remark}
It follows from Theorem \ref{symmetric thm} and Lemma \ref{secondchar2} that $f(x)$ has a newly reducible second iterate over $K$ if and only if $f$ is irreducible and the polynomial in \eqref{poly} has a root in $K$. This is because if $b = 0$ then $f(0) = 0$, and if $b = a+1$ then $f(1) = 0$, and thus $f$ is reducible in both cases.
\end{remark}

We now turn to finite fields of characteristic 2. In particular, one might ask which of the fields $\mathbb{F}_{2^n}$ belong to $\mathcal{N}_{2,k}$ for various $k \geq 2$? When $k \geq 4$, this question has been settled by Ahmadi et al \cite[Theorem 10]{shparostafe} (see also \cite{ahmadi}), and the answer is surprising: none of them. The result \cite[Theorem 10]{shparostafe} gives that if $f \in \mathbb{F}_{2^n}[x]$ is quadratic, then $f^3(x)$ is reducible over $\mathbb{F}_{2^n}$.

As for $\mathcal{N}_{2,2}$ and $\mathcal{N}_{2,3}$, observe that $\mathbb{F}_2 \not\in \mathcal{N}_{2,2}$, since the unique irreducible quadratic polynomial $f(x) = x^2 + x + 1 \in \mathbb{F}_2[x]$ has $f^2(x)$ irreducible. By adapting the methods of \cite{ahmadi}, we show the following.
\begin{proposition} \label{fintwo}
Let $K = \mathbb{F}_{2^n}$. Then $K \in \mathcal{N}_{2,2}$ for all $n \geq 2$ and $K \in \mathcal{N}_{2,3}$ for all $n \geq 1$.
\end{proposition}

\begin{proof}
Let $K = \mathbb{F}_{2^n}$, and $f(x) = x^2 + ax + b \in K[x]$, and denote the trace map $K \to \mathbb{F}_2$ by $\text{Tr}_{K/\mathbb{F}_2}$. By a standard result in field theory \cite[Corollary 3.6]{appfin}, we have that $f$ is irreducible over $K$ if and only if $a \neq 0$ and $\text{Tr}_{K/\mathbb{F}_2}(b/a^2) = 1$. Assume that this holds. Then by Capelli's Lemma (\cite{ahmadi, cohen2}), $f^2(x)$ is reducible over $K$ if and only if $f(x) - \alpha$ is reducible over $K(\alpha)$, where $f(\alpha) = 0$. Applying \cite[Corollary 3.6]{appfin} again, this is equivalent to $\text{Tr}_{K(\alpha)/\mathbb{F}_2}((b-\alpha)/a^2) = 0$. Properties of the trace now give
\begin{align*}
\text{Tr}_{K(\alpha)/\mathbb{F}_2}((b-\alpha)/a^2) &= \text{Tr}_{K(\alpha)/\mathbb{F}_2}(b/a^2) - \text{Tr}_{K(\alpha)/\mathbb{F}_2}(\alpha/a^2) \\
&= 2\text{Tr}_{K/\mathbb{F}_2}(b/a^2) + \text{Tr}_{K/\mathbb{F}_2}(\text{Tr}_{K(\alpha)/K}(\alpha/a^2)) \\
&= \text{Tr}_{K/\mathbb{F}_2}(a/a^2).
\end{align*}
To show that $K \in \mathcal{N}_{2,2}$, we thus seek $a,b \in K$ with $a \neq 0$, $\text{Tr}_{K/\mathbb{F}_2}(b/a^2) = 1$, and $\text{Tr}_{K/\mathbb{F}_2}(1/a) = 0$. If there are $r,s \in K$ with $r \neq 0$, $\text{Tr}_{K/\mathbb{F}_2}(r) = 0$, and $\text{Tr}_{K/\mathbb{F}_2}(s) = 1$, then we can take $a = 1/r$ and $b = s/r^2$. Now because $K$ is a separable extension of $\mathbb{F}_2$, the bilinear form $(x,y) = \text{Tr}_{K/\mathbb{F}_2}(xy)$ is non-degenerate, and taking $y=1$ gives that $\text{Tr}_{K/\mathbb{F}_2}$ is a surjective homomorphism from the additive group of $K$ to the additive group of $\mathbb{F}_2$. Therefore we can find the desired $r,s$ provided that $|K| > 2$.

To show that $K \in \mathcal{N}_{2,3}$, we seek $a,b \in K$ with $a \neq 0$ and $\text{Tr}_{K/\mathbb{F}_2}(b/a^2) = \text{Tr}_{K/\mathbb{F}_2}(1/a) = 1$. This ensures that $f^2(x)$ is irreducible over $K$, and then by \cite[Theorem 10]{shparostafe} we have that $f^3(x)$ is reducible over $K$. If we find $r \in K$ with $r \neq 0$ and $\text{Tr}_{K/\mathbb{F}_2}(r) = 1$, then we may take $a = b = 1/r$. Because $\text{Tr}_{K/\mathbb{F}_2} : K^+ \to \mathbb{F}_2^+$ is a surjective homomorphism, the desired $r$ must exist.
\end{proof}

The techniques of Proposition \ref{fintwo} and \cite[Theorem 10]{shparostafe} do not apply to general fields of characteristic 2, and it remains an open question which of them belong to $\mathcal{N}_{2,2}$ and $\mathcal{N}_{2,3}$. See Question \ref{charpquest}.

We now give a characterization of quadratic polynomials $f(x)$ such that $f^3\left(x+\gamma\right)$ is symmetrically reducible.
\begin{theorem} \label{symmchar}
Let $K$ be a field of characteristic not equal to 2, and let $f\left(x\right)=\left(x-\gamma\right)^2+\gamma+m\in K[x]$. Then $f^3\left(x+\gamma \right)$ is symmetrically reducible over $K$ if and only if $\gamma=-m$ or there exist $r,s\in K$ such that
\begin{align}
\gamma & = \frac{1}{256r^2}\left( -2 r^5 s^2+9 r^4 s^4-4 r^4 s^2-16 r^3 s^6+32 r^3 s^4+16 r^3 s^2+32 r^3+14 r^2 s^8\right. \label{gamma-in-r-s}\\
&\phantom{=\frac{1}{256r^2}(}-64 r^2 s^6-8 r^2 s^4+96 r^2 s^2+128 r^2-6 r s^{10}+48 r s^8-64 r s^6 \nonumber\\
&\phantom{=\frac{1}{256r^2}(}\left. -160 r s^4+96 r s^2+128 r+s^{12}-12 s^{10}+40 s^8-112 s^4-64 s^2 \right) \nonumber\\
m & = \frac{-r^2-2 r s^2-4 r+s^4-4 s^2-4}{8 r} \label{m-in-r-s}
\end{align}
\end{theorem}
\begin{proof}
First, if $\gamma=-m$ then $f^3(x+\gamma) = \left(x^4 - 2\gamma x^2 + \gamma ^2-\gamma\right)^2$, so $f^3(x+\gamma)$ is symmetrically reducible. Otherwise, suppose that there exist $r,s\in K$ that satisfy \eqref{gamma-in-r-s} and \eqref{m-in-r-s}. Note that $r\neq 0$. Then we compute $f^3\left(x+\gamma\right)$ and substitute the expressions in $r$ and $s$ for $\gamma$ and $m$. This gives
\begin{align*}
f^3\left(x+\gamma\right) &= \frac{h\left(x\right)}{64 r^2}\cdot\frac{h\left(-x\right)}{64r^2},
\end{align*}
where
\begin{align*}
h\left(x\right) &= 64 r^2 x^4 -64 r^2 s x^3 + \left(-16 r^3-64 r^2+16 r s^4-64 r s^2-64 r\right)x^2 \\&\phantom{=} +\left(24 r^3 s+64 r^2 s-8 r s^5+32 r s^3+32 r s\right)x + \left(r^4-12 r^3 s^2+10 r^2 s^4 \right.\\&\phantom{=}\left. -24 r^2 s^2-8 r^2 -4 r s^6+16 r s^4+16 r s^2+s^8-8 s^6+8 s^4+32 s^2+16\right).
\end{align*}
Therefore, $f^3\left(x+\gamma\right)$ is symmetrically reducible.

Conversely, suppose that $f^3\left(x+\gamma\right)$ is symmetrically reducible. Our goal is to show either that $\gamma=-m$ or that \eqref{gamma-in-r-s} and \eqref{m-in-r-s} hold for some $r,s\in K$. We hav
\begin{align*}
f^3\left(x+\gamma\right)
&= x^8 +4 m x^6 +\left(6 m^2+2 m\right)x^4 +\left(4 m^3+4 m^2\right)x^2 + (m^4 + 2 m^3 + m^2 + m + \gamma).
\end{align*}
Since $f^3\left(x+\gamma\right)$ is symmetrically reducible, by definition there exist $h_1,h_2\in K[x]$ such that
\begin{align*}
h_1\left(x+\gamma\right) &= x^4+dx^3+cx^2+bx+a, \\
h_2\left(x+\gamma\right) &= x^4-dx^3+cx^2-bx+a,
\end{align*}
and $f^3\left(x+\gamma\right)=h_1\left(x+\gamma\right)h_2\left(x+\gamma\right)$. Equating coefficients gives
\begin{align}
\gamma +m^4+2 m^3+m^2+m &= a^2 \label{f3eq1}\\
4 m^3+4 m^2 &= 2 a c-b^2  \\
6 m^2+2 m &= 2 a-2 b d+c^2 \label{f3eq3}\\
4 m &= 2 c-d^2 \label{f3eq4}
\end{align}
First, we solve \eqref{f3eq3} for $a$ and substitute the resulting value into the other equations. We then solve \eqref{f3eq4} for $c$, and substitute the resulting value into the other equations. This gives \small
\begin{align}
 \gamma +m^4+2 m^3+m^2+m &= \frac{1}{64} \left(64 b^2 d^2-16 b d^5-128 b d^3 m+128 b d m^2+128 b d m+d^8+16 d^6 m \nonumber
  \right.\\&\phantom{=\frac{1}{64}(6}\left. +48 d^4 m^2 -16 d^4 m-128 d^2 m^3 -128 d^2 m^2+64 m^4+128 m^3+64m^2\right) \label{f3eq5} \\
4\left(m^3+m^2\right) &= \frac{1}{8} \left(-8 b^2+8 b d^3+32 b d m-d^6-12 d^4 m-24 d^2 m^2+8 d^2 m+32 m^3+32 m^2\right) \label{f3eq6}
\end{align}
\normalsize
Next we solve \eqref{f3eq6} for $b$, giving
\begin{align}\label{b-eq}
b &= \frac{1}{4} \left(2 d^3 + 8 d m \pm d \beta \right) \quad \text{where} \quad \beta = \sqrt{2 d^{4} + 8 d^{2} m + 16 m^{2} + 16 m}.
\end{align}
Finally, substituting \eqref{b-eq} into \eqref{f3eq5} and solving for $\gamma$ gives
\begin{equation}\label{gamma}
\gamma=\pm \beta A + \frac{17}{64}d^{8} + \frac{5}{4}md^6 + \frac{11}{4}  m^{2}d^{4} + \frac{7}{4} md^{4} + 2 m^{3}d^{2}  + 2 m^{2}d^{2}  - m,
\end{equation}
where
\begin{equation*}
A = \frac{1}{16} d^2 \left(3 d^4+8 d^2 m+8 m^2+8 m\right).
\end{equation*}
Since $b \in K$, \eqref{b-eq} shows that $d=0$ or $\beta \in K$. If $d=0$, then \eqref{gamma} gives $\gamma = -m$ and we're done.
Otherwise, $\beta\in K$. Then \eqref{b-eq} gives a $K$-rational point on the surface
\begin{align}
S: y^2 &= 2s^4 + 8s^2m + 16m^2 + 16m. \label{mainsurface}
\end{align}
For fixed $s$, this equation is a conic in $y$ and $m$, and we can use rational projection to parametrize its $K$-rational points. The homogeneous form of $S$ can be written
\begin{align*}
\overline{S}: Y^2 &= 16 M^2 + \left(8 s^2 + 16\right) M Z + 2s^4 Z^2.
\end{align*}
Note that the rational point $P = [M:Y:Z] = [1:4:0]\in\mathbb{P}^2(K)$ is on $\overline{S}$; this is the point we will project through. Let $r_0\in K$ be arbitrary. The affine part of the line through $r_0$ and $P$ is given by
$y = 4 m + r_0$. To solve for the other intersection point, we substitute into \eqref{mainsurface} to get
$
\left(4 m + r_0\right)^2 = 16 m^2 + \left(8s^2+16\right)m + 2s^4.
$
This equation is linear in $m$, and we get
\begin{align*}
m &= \frac{2 s^4-r_0^2}{8 \left(r_0-s^2-2\right)}.
\end{align*}
Note that $r_0-s^2-2=0$ corresponds to the line intersecting the point at infinity with multiplicity 2, so with this projection we don't miss any affine rational points. Taking $r=r_0-s^2-2$ now gives
\begin{align*}
m = m\left(r,s\right) &= \frac{2s^4-r^2-2 r s^2-4 r-s^4-4 s^2-4}{8r}
\end{align*}
This is \eqref{m-in-r-s}, so we're halfway done. Since we have $y=4m+r_0=4m+r+s^2+2$,
\begin{align*}
y = y\left(r,s\right) &= \frac{r^2+s^4-4 s^2-4}{2 r}
\end{align*}
Recall that $y=\beta$ and $s=d$. Using these in \eqref{gamma},
\begin{align*}
\gamma=\pm y A + \frac{17}{64}s^{8} + \frac{5}{4}ms^6 + \frac{11}{4}  m^{2}s^{4} + \frac{7}{4} ms^{4} + 2 m^{3}s^{2}  + 2 m^{2}s^{2}  - m.
\end{align*}
If the plus sign satisfies this equation, it simplifies to
\begin{align*}
\gamma = \gamma\left(r,s\right) & = \frac{1}{256r^2}\left(-2 r^5 s^2+9 r^4 s^4-4 r^4 s^2-16 r^3 s^6+32 r^3 s^4+16 r^3 s^2+32 r^3+14 r^2 s^8 \right. \\
&\phantom{=\frac{1}{256r^2}(}-64 r^2 s^6-8 r^2 s^4+96 r^2 s^2+128 r^2-6 r s^{10}+48 r s^8-64 r s^6 \nonumber\\
&\phantom{=\frac{1}{256r^2}(}\left. -160 r s^4+96 r s^2+128 r+s^{12}-12 s^{10}+40 s^8-112 s^4-64 s^2\right). \nonumber
\end{align*}
This is \eqref{gamma-in-r-s}. On the other hand, if we must choose the minus sign, let
\begin{align*}
r' &= \frac{-s^4+4 s^2+4}{r} \\ s' &= s
\end{align*}
Plugging these in, $m\left(r',s'\right)=m\left(r,s\right)$ and $y\left(r',s'\right)=-y\left(r,s\right)$. So in this case, \eqref{gamma-in-r-s} is satisfied by $r'$ and $s'$.
\end{proof}

We make a few observations about Theorem \ref{symmchar} and its proof. First, Theorem \ref{characterization} is now proved, as a consequence of Theorem \ref{symmetric thm} and Theorems \ref{second-iterate-criterion} and \ref{symmchar}. Second, whenever $\gamma=-m$, we have $f(x)=(x-\gamma)^2$, which is obviously reducible. Therefore, Theorem \ref{symmchar} shows that all $f(x)$ with newly reducible third iterate have $m,\gamma$ such that \eqref{gamma-in-r-s} and \eqref{m-in-r-s} are satisfied by some $r,s\in K$.

Third, we emphasize that $K$-rational points on the surface $S$ described in \eqref{mainsurface} play a crucial role in determining the existence of quadratic polynomials over $K$ with newly reducible third iterate. The fibers of the map $S \to \mathbb{A}^1$ given by projection onto the $m$-coordinate are particularly interesting, and worth dwelling on for a moment. In general, these fibers are elliptic curves, unless $m$ divides the discriminant of $2s^4 + 8s^2m + 16m^2 + 16m$ considered as a polynomial in $s$. Observe that
$$
\Disc(2s^4 + 8s^2m + 16m^2 + 16m) = 2^{21}m^3(m+2)^2(m+1).
$$
Hence we obtain a curve of genus zero if and only if $m \in \{0, -1, -2\}$. We remark that for each of these $m$-values, $f(x) = (x-\gamma)^2 + \gamma + m$ is post-critically finite -- that is, the forward orbit of $\gamma$ is finite -- and indeed these are the only such $m$-values for $K = \Q$.

When $m = 0$, we have that $S$ degenerates to $y^2 = 2s^4$, which has no $K$-rational points unless $\sqrt{2} \in K$. Hence if $\sqrt{2} \not\in K$, then $f(x) = (x - \gamma)^2 + \gamma$ has $f^3(x + \gamma)$ symmetrically reducible only if $\gamma = -m$, which implies that $f(x)$ is reducible. If $\sqrt{2} \in K$, then equations \eqref{f3eq1} - \eqref{f3eq4} yield
$
\gamma = d^8 \left( \frac{17}{64} \pm \frac{3 \sqrt{2}}{16} \right),
$
and it follows that
$$\sqrt{2\sqrt{\gamma}} = \frac{d^2}{2}(1 + \sqrt{2}) \in K.$$
By Lemma \ref{second-iterate-criterion} this shows that $f^2(x)$ is reducible over $K$. This is in line with \cite[Theorem 1.2(i) and Remark 4.3]{goksel}, which give that when $K = \Q$ and $m = 0$, either $f^2(x)$ is reducible over $K$, or all iterates of $f$ are irreducible over $K$.

When $m = -2$, $S$ degenerates to $y^2 = 2(s^2-2)^2$, and again one finds that $f^3(x+\gamma)$ is symmetrically reducible if and only if either $f(x)$ or $f^2(x)$ is reducible.

When $m = -1$, matters are different, and the following proposition shows we obtain an infinite family of polynomials with newly reducible third iterate for certain $K$.

\begin{proposition} \label{m-1}
Let $K$ be a field of characteristic not equal to 2, and let $f(x) = (x-\gamma)^2 + \gamma - 1$. Then $f^3(x+\gamma)$ is symmetrically reducible if and only if
\begin{equation} \label{gammaeq}
\gamma = 4\left(\frac{(t^2 - 8t + 8)(t^2 + 8)}{(t^2 - 8)^2}\right)^4 + 1 \quad \text{for some $t \in K$.}
\end{equation}
Suppose moreover that $K$ is a field with a non-trivial discrete valuation $v$ such that $v(5)$ is odd, and that $-1$ is not a square in $K$.
Then $f^3(x)$ is newly reducible over $K$ whenever $t = 25r$ in \eqref{gammaeq}, for any $r \in K$ with $v(r) \geq 0$. In particular, $K \in \mathcal{N}_{2,3}^\infty$.
\end{proposition}

\begin{proof}

When $m = -1$, $S$ degenerates to $y^2 = 2s^2(s^2 - 4)$, whose solutions are given by $s = -2\frac{t^2 + 8}{t^2-8}$ where $t$ is any element of $K$. In the notation of the proof of Theorem \ref{symmchar} we have $s = d$, and we now use equations \eqref{f3eq1} - \eqref{f3eq4} to find \eqref{gammaeq}.

 Observe that for $m = -1$ we have $-m - \gamma = 1-\gamma$, and from \eqref{gammaeq} this is not a square in $K$ because $-1$ is not a square in $K$. We also have $m^2 + m + \gamma = \gamma$, and so if $\gamma$ is not a square in $K$, then Lemma \ref{second-iterate-criterion} gives that $f^2(x)$ is irreducible over $K$. Now from \eqref{gammaeq} we have that the numerator of $\gamma$ as a function of $t$ has constant term $c=5 \cdot 8^8$. Let $v : K^* \to \Z$ be a normalized discrete valuation with $v(5)$ odd. Then $v(c) = v(5) + 8v(8)$, and so $v(c)$ is odd. Taking $t = 25r$ with $v(r) \geq 0$ then shows $v(\gamma)$ is odd, and hence $\gamma$ cannot be a square in $K$.

\end{proof}

\begin{remark}
Proposition \ref{m-1} shows that if $K$ is a number field with ring of integers $\mathcal{O}_K$, $-1$ is not a square in $K$, and the ideal $(5)$ is not the square of an ideal in $\mathcal{O}_K$, then $K \in \mathcal{N}_{2,3}^\infty$. In particular, $\Q \in \mathcal{N}_{2,3}^\infty$.
\end{remark}

The family appearing in Proposition \ref{m-1} is the same family that appears in Lemma 3.9 of \cite{goksel}, though there it is stated only for integer parameters.

In the preceding discussion, we took $m$ to be constant.
It is of interest to give another infinite family with non-constant $m$-values, which we do in the following corollary to Theorem \ref{symmchar}.

\begin{corollary}\label{newfamily}
Let $K$ be a totally real number field and suppose there is a prime ${\mathfrak{p}}$ of $\mathcal{O}_K$ lying over the ideal $(3)$ with residue degree 1.  Let $f(x) = (x - \gamma)^2 + \gamma + m$, where
\begin{align*}
m &=  -1 - 2 k^2 + k^4, \\
\gamma &= 1 + k^2 \left(-2 + k^2\right) \left(-1 - 4 k^2 + 2 k^4\right) \left(1 - 4 k^2 + 2 k^4\right)
\end{align*}
for $k \in \mathcal{O}_K$. If $|k| > \sqrt{2}$ and ${\mathfrak{p}} \nmid (k)$, then $f^3(x)$ has a newly reducible third iterate over $K$. In particular, $K \in \mathcal{N}_{2,3}^\infty$.
\end{corollary}
\begin{proof}
Let $r=2$ and $s=2k$ in Theorem \ref{symmchar}, and observe that it's sufficient to show that neither of $-m-\gamma$ and $m^2 + m + \gamma$ is a square in $K$. One calculates $-m-\gamma = -4 k^6 \left(-2 + k^2\right)^3$.
Because $|k| > \sqrt{2}$, we have $-m-\gamma < 0$, and since $K$ is totally real, we have that $-m-\gamma$ is not a square in $K$.
We also have
\begin{align*}
m^2 + m + \gamma &= \left(-1 + k^2\right)^2 \left(1 + 6 k^2 + 13 k^4 - 16 k^6 + 4 k^8\right)
\end{align*}
Observe that because ${\mathfrak{p}}$ has residue degree 1, we have $\mathcal{O}_K/{\mathfrak{p}} \cong \Z / 3\Z$. Because ${\mathfrak{p}} \nmid (k)$ we conclude $k^2 \equiv 1 \pmod{{\mathfrak{p}}}$, whence
\begin{equation*}
1 + 6 k^2 + 13 k^4 - 16 k^6 + 4 k^8 \equiv 2 \pmod{{\mathfrak{p}}},
\end{equation*}
proving that $m^2 + m + \gamma$ is not a square in $\mathcal{O}_K$, and hence not a square in $K$
\end{proof}
Note that many other infinite families exist; for instance, there is a similar one with $r=-2$.

\section{Rational Quadratics with newly reducible fourth iterate} \label{fourth}

In this section we briefly examine the following question.
\begin{question}
Does there exist $f \in \Q[x]$ of degree $2$ with newly reducible fourth iterate?
\end{question}

Theorem \ref{symmetric thm} tells us that if $f(x)=(x-\gamma)^2+m+\gamma$ has a newly reducible fourth iterate, then $f^4(x+\gamma)$ is symmetrically reducible. That is, $f^4(x+\gamma) = h(x)h(-x),$ where
\begin{align*}
h(x) &= x^8	+a_7 x^7	+a_6 x^6	+a_5 x^5	+a_4 x^4	+a_3 x^3	+a_2 x^2	+a_1 x	+a_0
\end{align*}
Equating coefficients gives the following system of equations:
\begin{align*}
\gamma +m^8+4 m^7+6 m^6+6 m^5+5 m^4+2 m^3+m^2+m &= a_0^2 \\
8 m^7+24 m^6+24 m^5+16 m^4+8 m^3 &= 2 a_0 a_2-a_1^2 \\
28 m^6+60 m^5+36 m^4+16 m^3+4 m^2 &= a_2^2-2 a_1 a_3+2 a_0 a_4 \\
56 m^5+80 m^4+24 m^3+8 m^2 &= -a_3^2+2 a_2 a_4-2 a_1 a_5+2 a_0 a_6 \\
70 m^4+60 m^3+6 m^2+2 m &= a_4^2+2 a_0-2 a_3 a_5+2 a_2 a_6-2 a_1 a_7 \\
56 m^3+24 m^2 &= -a_5^2+2 a_2+2 a_4 a_6-2 a_3 a_7 \\
28 m^2+4 m &= a_6^2+2 a_4-2 a_5 a_7 \\
8 m &= 2 a_6-a_7^2
\end{align*}
This defines a surface, which we denote $S_4$, whose rational points correspond to quadratic $f(x) \in \Q[x]$ such that $f^4(x+\gamma)$ is symmetrically reducible over $\Q$. In particular, this includes all monic quadratic $f(x)$ with newly reducible fourth iterate over $\Q$. At present, obtaining any information about the rational points on $S_4$ seems to be quite difficult.

It is interesting to consider the fibers of the map $\pi : S_4 \mapsto \mathbb{A}^1$ given by projection onto the $m$-coordinate. In general these fibers appear to be very high-genus curves, though it is possible that some particular fibers are sufficiently singular that the genus drops significantly. This was the case for the fibers $m \in \{0, -1, -2\}$ for the similar projection map of the curve $S$ in \eqref{mainsurface} that arose from studying $f$ with symmetrically reducible third iterate. Similarly, the fibers $\pi^{-1}(0)$ and $\pi^{-1}(-2)$ are likely to have smaller genera than other fibers; unfortunately, their rational points cannot lead to $f(x) \in \Q[x]$ with newly reducible fourth iterate (see \cite[Theorem 4.1 and Remark 4.3]{goksel}).

The fiber $\pi^{-1}(-1)$, on the other hand, remains a possible location for rational points on $S_4$. Lemma 3.10 of \cite{goksel} shows that \textit{integer} values of $\gamma$ cannot lead to rational points on this fiber, but the case of non-integer $\gamma \in \Q$ remains open (see Remark 4.3 of \cite{goksel}).

\section{Cubics with newly reducible second iterate} \label{cubics}

In this section we study the question of whether there exist monic irreducible cubic polynomials in $K[x]$ with newly reducible second iterate. Observe that if $K$ does not have characteristic 3, then $f(x) = x^3 + c_2x^2 + c_1x + c_0$ can be conjugated to a cubic with no quadratic term: taking $\gamma = -c_2/3$, $a = c_1 - (c_2^2/3) + (c_2/3)$, and $b = c_0 - (c_2^3/27) -(c_2/3)$ gives
 $f(x) = (x-\gamma)^3 + a(x-\gamma) + b + \gamma$, or equivalently
\begin{equation} \label{cubicform}
f(x + \gamma) = x^3 + ax + b + \gamma \in K[x]
\end{equation}
We emphasize that, unlike in the quadratic case, $\gamma$ is \textit{not} a critical point for $f$ except in the special case that $a=0$. It is this special case that is the source of our main result on newly reducible cubics:

\begin{theorem} \label{cubic2}
Let $K$ be a field of characteristic not equal to 2 or 3, and assume that $2$ is not the cube of an element of $K$. If $f(x) = (x-\gamma)^3 + b + \gamma$, with
\begin{align}
b = 36t^3 & \quad \text{and} \quad \gamma = (-2 \cdot 6^6)t^9 - 36t^3, \quad \text{or}  \label{firstfamily}\\
b = -9t^3& \quad \text{and} \quad  \gamma = 9t^3 + (2\cdot 3^6)t^9 \label{secondfamily}
\end{align}
for some non-zero $t \in K$, then $f$ has a newly reducible second iterate.
\end{theorem}

\begin{proof}
From Lemma \ref{deglem}, we have that if $f^2(x)$ is newly reducible, then $f^2(x) =  p_1(x)p_2(x)$, where $\deg p_1(x) = 3$, $\deg p_2(x) = 6$, and $p_1(x)$ is irreducible.
Write
$$f^2(x + \gamma) = (x^3 + a_2x^2 + a_1x + a_0)(x^6 + b_5x^5 + b_4x^4 + b_3x^3 + b_2x^2 + b_1x + b_0),$$ and equate coefficients to get the system of equations
\begin{align*}
b + b^3 + \gamma & =  a_0b_0 \\
0 & =  a_0b_1 + a_1b_0 \\
0 & =  a_0b_2 + a_1b_1 + a_2b_0 \\
3b^2 & =  a_0b_3 + a_1b_2 + a_2b_1 + b_0 \\
0 & =  a_0 b_4 + a_1 b_3 + a_2 b_2 + b_1 \\
0 & =  a_0b_5 + a_1 b_4 + a_2 b_3 + b_2 \\
3b & =  a_0 + a_1 b_5 + a_2 b_4 + b_3 \\
0 & =  a_1 + a_2 b_5 + b_4 \\
0 & =  a_2 + b_5
\end{align*}
We may solve for $b_5, b_4, b_3, b_2, b_1$, and $b_0$ to obtain
\begin{align}
b + b^3 + \gamma &=  a_0^3 - a_0 a_1^3 - 6 a_0^2 a_1 a_2 + 6 a_0 a_1^2 a_2^2 + 4 a_0^2 a_2^3 -
 5 a_0 a_1 a_2^4 + a_0 a_2^6 \label{cubic_equation_1} \\
 &\phantom{=} - 3 a_0^2 b + 6 a_0 a_1 a_2 b - 3 a_0 a_2^3 b +
 3 a_0 b^2 \nonumber \\
 0 &=  3 a_0^2 a_1 - a_1^4 - 9 a_0 a_1^2 a_2 - 3 a_0^2 a_2^2 + 6 a_1^3 a_2^2 +
 8 a_0 a_1 a_2^3 - 5 a_1^2 a_2^4 \label{cubic_equation_2} \\
 &\phantom{=} - a_0 a_2^5 + a_1 a_2^6 - 6 a_0 a_1 b +
 6 a_1^2 a_2 b + 3 a_0 a_2^2 b \nonumber \\
 &\phantom{=} - 3 a_1 a_2^3 b + 3 a_1 b^2 \nonumber \\
 0 &=  3 a_0 a_1^2 + 3 a_0^2 a_2 - 4 a_1^3 a_2 - 12 a_0 a_1 a_2^2 + 10 a_1^2 a_2^3 +
 5 a_0 a_2^4 - 6 a_1 a_2^5 + a_2^7  \label{cubic_equation_3}\\
 &\phantom{=} - 3 a_1^2 b - 6 a_0 a_2 b + 9 a_1 a_2^2 b -
 3 a_2^4 b + 3 a_2 b^2 \nonumber
\end{align}

Note that \eqref{cubic_equation_2} is quadratic in $a_0$.  Taking the discriminant of this equation in $a_0$ gives
\begin{equation} \label{discriminant1}
9 a_2^4 b^2 + \left(36 a_1^3 a_2 - 42 a_1^2 a_2^3 + 24 a_1 a_2^5 -
 6 a_2^7\right)b + 12 a_1^5 - 3 a_1^4 a_2^2 - 12 a_1^3 a_2^4 + 10 a_1^2 a_2^6 -
 4 a_1 a_2^8 + a_2^{10},
\end{equation}
which in turn is quadratic in $b$.  The discriminant of this equation in $b$ is
\begin{align*}
& -144 \left(-9 a_1^6 a_2^2 + 24 a_1^5 a_2^4 - 25 a_1^4 a_2^6 + 14 a_1^3 a_2^8 -
   5 a_1^2 a_2^{10} + a_1 a_2^{12}\right),
\end{align*}
which has the semi-miraculous factorization
\begin{equation} \label{semimiracle}
12\left(a_1 a_2^2 \left(a_1 - a_2^2\right) \left(9 a_1^4 - 15 a_1^3 a_2^2 + 10 a_1^2 a_2^4 -
   4 a_1 a_2^6 + a_2^8\right)\right).
\end{equation}
Selecting $a_1$ and $a_2$ so that \eqref{semimiracle} vanishes forces \eqref{discriminant1} to be a square, which gives a $K$-rational solution to \eqref{cubic_equation_2}

One such choice is $a_1 = 0$, which reduces \eqref{cubic_equation_2} to $0 = -3a_0^2 a_2^2 - a_0 a_2^5 + 3a_0a_2^2 b$, giving $a_0 = 0$ or $a_0 = -a_2^3/3 + b$. The former is impossible, since $a_0 = 0$ would imply that the degree three factor of our second iterate was reducible. Substituting $a_0 = -a_2^3/3 + b$ into \eqref{cubic_equation_3} produces another semi-miracle, as the $b^2$ terms cancel, giving $b = a_2^3/6$. Then \eqref{cubic_equation_1} gives $\gamma = \frac{-a_2^9}{108} - \frac{a_2^3}{6}$. Observe that $b+\gamma = -a_2^9/108$, and it follows from Proposition \ref{xd-c-fact} that $f(x)$ is irreducible over $K$ if and only if $a_2^9/108$ is not a cube in $K$. But this is equivalent to $2$ not being a cube in $K$, which is true by assumption. Hence $f$ has newly reducible second iterate over $K$. Taking $a_2 = 6t$ clears denominators and gives the family in \eqref{firstfamily}.

Taking $a_1 = a_2^2$ also causes \eqref{semimiracle} to vanish, and reduces \eqref{cubic_equation_2} to
$$0 =   - 2 a_0 a_2^5 + a_2^8 - 3 a_0 a_2^2 b + 3 a_2^5 b + 3 a_2^2 b^2,$$
giving $a_0 = \frac{a_2^8  + 3 a_2^5 b  + 3 a_2^2 b^2}{2a_2^5 + 3 a_2^2 b}$. Substituting this into \eqref{cubic_equation_3} yields a cubic rational function in $b$ whose numerator fortuitously has a factor of $a_2^3 + 3b$, together with an irreducible quadratic in $b$. Taking $b = -a_2^3/3$ and substituting this into \eqref{cubic_equation_1} yields $\gamma = \frac{1}{27}(9a_2^3 + 2a_2^9)$. Then $b + \gamma = \frac{2}{27}a_2^3,$ and we have that $f$ is irreducible over $K$ since $-2$ is not a cube in $K$, as in the previous paragraph. Taking $a_2 = 3t$ gives the family in \eqref{secondfamily}.
\end{proof}

\begin{remark}
We can also force \eqref{semimiracle} to vanish by taking $a_2 = 0$, but this produces a family of polynomials $f$ that are all reducible.
\end{remark}

We obtain the following immediate corollary of Theorem \ref{cubic2}.
\begin{corollary}
If $K$ is an infinite field satisfying the hypotheses of Theorem \ref{cubic2}, then $K \in \mathcal{N}_{3,2}^{\infty}$. In particular, if $K$ is a number field with ring of integers $\mathcal{O}_K$ and $2$ is not the cube of an element of $K$, then there are infinitely many monic $f \in \mathcal{O}_K[x]$ with newly reducible second iterate over $K$.
\end{corollary}

We end this section by briefly addressing Question \ref{cubicquest}, which asks whether there exists $f \in \Q[x]$ with $f$ irreducible but $f^2(x)$ a product of three irreducible cubics. The following proposition addresses this.
\begin{proposition}
Let $f(x) \in \Q[x]$ have the form in \eqref{cubicform}, and assume that $f^2$ is newly reducible with three distinct factors $p_1(x)$, $p_2(x)$, and $p_3(x)$. If $\beta$ is a root of $p_1(x + \gamma)$, and $f(\beta + \gamma) = \alpha$, then
\begin{equation}\label{cubic-root-equation}
\frac{-\beta \pm \sqrt{-3\beta^2 - 4a}}{2}
\end{equation}
are roots of $p_2(x + \gamma)$ and $p_3(x + \gamma)$ whose images under $f(x + \gamma)$ are $\alpha$.
\end{proposition}
\begin{proof}
Observe that each of the factors of $f^2(x)$ must have degree 3. By Lemma \ref{surjective}, the map
\begin{equation*}
    \Phi_i : \{\text{roots of $p_i(x + \gamma)$ in $\overline{\Q}$}\} \to \{\text{roots of $f(x + \gamma)$ in $\overline{\Q}$}\}
\end{equation*}
is surjective and $k$-to-one for each factor $p_i$ of $f^2$.  Since each factor has degree equal to the degree of $f$, $\Phi_i$ is one-to-one for each $i$, and there is one root of each $p_i$ whose image under $\Phi_i$ is $\alpha$.  We will show that the roots of $p_2(x + \gamma)$ and $p_3(x + \gamma)$ satisfying this condition are given by \eqref{cubic-root-equation}.  Let $\beta_2$ be such a root, and observe that $\beta^3 + a(\beta) + b + \gamma = \beta_2^3 + a(\beta_2) + b + \gamma$. This gives a quadratic equation for $\beta_2$ in terms of $\beta$, to which we apply the quadratic formula to obtain $\beta_2  = \frac{-\beta \pm \sqrt{-3\beta^2 - 4a}}{2}$.
\end{proof}
Notice that when $a = 0$, we have that the three preimages of a given root of $f(x)$ differ by a multiple of a cubic root of unity.  As a consequence, the constant terms of our three factors are equal.  We can equate coefficients using this fact to show that there is no monic cubic over $\Q[x]$ with $a = 0$ whose second iterate factors as three irreducible cubics.

\section{Quartics with newly reducible second iterate} \label{quartics}

In this section we will find infinitely many rational quartic polynomials with newly reducible second iterate. We use the following standard fact from field theory to tell when $f(x)$ is irreducible, noting that $K^n$ refers to the set $\{k^n : k \in K\}$.

\begin{proposition}{\cite[Theorem 8.1.6]{karpilovsky}}\label{xd-c-fact}
Let $K$ be a field and $f(x)=x^d-c \in K[x]$ for $d \geq 1$. Then $f(x)$ is irreducible over $K$ if and only if $c\not\in K^p$ for all primes $p\mid d$ and $c\not\in -4K^4$ whenever $4\mid d$.
\end{proposition}

Throughout this section, we only consider polynomials of the form $f(x)=(x-\gamma)^4+m+\gamma$, where $m,\gamma\in\Q$. It turns out that even in this subset of rational quartics, there are infinitely many with newly reducible second iterate.

\begin{theorem}\label{deg-4-theorem}
Let $K$ be a field of characteristic not equal to $2$, and suppose that $f(x)=(x-\gamma)^4+m+\gamma\in K[x]$ is irreducible over $K$. Then $f^2(x+\gamma)$ factors as $p(x^2)p(-x^2)$ for some $p \in K[x]$
if and only if there exist $r,s\in K$ with $r \neq s^2$ such that $m = \frac{2 s^4-r^2}{8 r-8 s^2}$ and
$$
\gamma = \frac{-2 r^5 s^2+19 r^4 s^4-72 r^3 s^6+32 r^3+136 r^2 s^8-32 r^2 s^2-128 r s^{10}-64 r s^4+48 s^{12}+64 s^6}{256\left(r-s^2\right)^2}.
$$
\end{theorem}

\begin{proof}
First, suppose that $m,\gamma,r,s\in K$ are as in the theorem. Writing $r_1=r-s^2$, one computes that
$f^2(x+\gamma) = p(x^2)p(-x^2)$,
where
\begin{align*}
 p(x) &=  x^4
-  s x^3
+ \frac{s^4- r_1^2}{4 r_1}x^2
+ \frac{3 r_1^2 s- s^5}{8 r_1}x
+ \frac{-4 r_1 s^6+10 r_1^2 s^4-12 r_1^3 s^2+r_1^4+s^8}{64 r_1^2}.
\end{align*}

Conversely, suppose that $f(x+\gamma) = p(x^2)p(-x^2)$ for $p(x) = x^4 + dx^3 + cx^2 + bx + a$.
Equating coefficients gives us the following system of equations.
\begin{align}
a^2 &= \gamma +m^4+m \label{deg-4-1}\\
2 a c-b^2 &= 4 m^3 \label{deg-4-2}\\
2 a-2 b d+c^2 &= 6 m^2 \label{deg-4-3}\\
2 c-d^2 &= 4 m\label{deg-4-4}
\end{align}
First we solve \eqref{deg-4-4} for $c$ and substitute this into the other equations. Then we do the same with \eqref{deg-4-3} and $b$, and then with \eqref{deg-4-2} and $a$. This leaves us with equation \eqref{deg-4-1}, which becomes
\begin{align}
\frac{1}{64} \left(3 d^4+8 d^2 m \pm 2 \sqrt{2} d^2 \sqrt{d^4+4 d^2 m+8 m^2}+8 m^2\right)^2 &= \gamma +m^4+m \label{deg-4-large-gamma}
\end{align}
Note that if $d=0$, this equation becomes $m^4=\gamma+m^4+m$, so $\gamma+m=0$. But then $f(x+\gamma)=x^4$ is reducible, which is a contradiction. So $d\neq 0$, and therefore $\sqrt{2(d^4+4 d^2 m+8 m^2)}\in K$. In other words, there is a $K$-rational solution to
\begin{align*}
S:\ y^2 &= 16m^2 + 8s^2m + 2s^4
\end{align*}
where $s=d$. Considering $s$ as a fixed parameter, this is a conic. Now we want to parametrize the $K$-rational points on $S$, so we use the homogeneous form of $S$,
\begin{align*}
\overline{S}:\ Y^2 &= 16M^2 + 8s^2MZ + 2s^4Z^2,
\end{align*}
and project through the point at infinity $[M:Y:Z]=[1:4:0]$. We parametrize our projection line by where it crosses the line $m=0$. Then each line is given by $y=4m+r$, where $r\in\Q$. Then we solve for the intersection point of this line and $S$, giving
\begin{align*}
m = m(r,s) = \frac{2s^4-r^2}{8r-8s^2} \qquad \text{and} \qquad
y = y(r,s) = 4m+r = \frac{2s^4-r^2}{2r-2s^2}-r.
\end{align*}
Then we solve \eqref{deg-4-large-gamma} for $\gamma$. This gives
\begin{align*}
\gamma &= \frac{1}{64} \left(3 s^4 + 8 s^2 m + 8 m^2 \pm 2 s^2 y\right)^2 - m - m^4.
\end{align*}
Note that there is a plus or minus in front of the term containing $y$. The plus sign leads to
\begin{align*}
\gamma &= \frac{-2 r^5 s^2+19 r^4 s^4-72 r^3 s^6+32 r^3+136 r^2 s^8-32 r^2 s^2-128 r s^{10}-64 r s^4+48 s^{12}+64 s^6}{256\left(r-s^2\right)^2},
\end{align*}
as in the statement of the theorem. On the other hand, if we have the minus sign, let
$r^\prime = \frac{s^2 \left(2 s^2-r\right)}{s^2-r}$  and $s^\prime = s$.
Plugging these in gives $m(r^\prime,s^\prime) = m(r,s)$ and $y(r^\prime,s^\prime)=-y(r,s)$, again giving $\gamma$ as in the statement of the theorem.
\end{proof}

\begin{corollary}\label{deg-4-infinitely-many-corollary}
Let $K$ be a field of characteristic not in $\{2,3\}$, and such that $3$ is not the square of an element of $K$, and $-3$ is not the fourth power of an element of $K$. Then
$$
f(x) = \left(x + 192t^8 - 7t^2\right)^4 -192t^8
$$
has a newly reducible second iterate over $K$ for any $t \in K, t \neq 0$.
\end{corollary}
\begin{proof}
Let $t \in K$ with $t \neq 0$, and in the notation of Theorem \ref{deg-4-theorem} take $r=48t^2$ and $s=4t$. Then from Theorem \ref{deg-4-theorem} we have that $f(x) = (x-\gamma)^4 + \gamma + m$ has $f^2(x)$ reducible over $K$, where $\gamma = -192t^8 + 7t^2$ and $m = -7t^2$. We now argue that $f(x)$ is irreducible over $K$. Observe that
$$-m-\gamma = 192t^8 = (2^6 \cdot 3)t^8$$
If $-m-\gamma = k^2$ for $k \in K$, then $3$ is a square in $K$, while if $4(m+\gamma) = k^4$ for $k \in K$, then $-3$ is a fourth power in $K$. Proposition \ref{xd-c-fact} now shows that $f$ is irreducible over $K$.
\end{proof}

\section{Higher-degree polynomials with newly reducible second iterate}\label{largerdegrees}

\begin{theorem} \label{infinited}
Let $d\equiv 2 \pmod{4}$, let $p_1, \ldots, p_r$ be the distinct odd primes dividing $d$, let $K$ be a field, and put $K^n = \{k^n : k \in K\}$. Assume that $-1 \not\in K^2$ and $-4k^4 \not\in K^{p_i}$ for each $i = 1, \ldots, r$. Then
the polynomial $f(x) = \left(x-4k^4\right)^d+4k^4$ has newly reducible second iterate over $K$.
\end{theorem}
\begin{proof}
We have $f^2(x) = \left(x-4k^4\right)^{d^2}+4k^4$.  Because $4 \mid d^2$, Proposition \ref{xd-c-fact} gives that $f^2(x)$ is reducible over $K$.
To show that $f(x)$ is irreducible over $K$, we note that $4 \nmid d$, and so by Proposition \ref{xd-c-fact} it is enough to show that $-4k^4 \not\in K^2$ and $-4k^4 \not\in K^{p_i}$ for each $i = 1, \ldots, r$. The former follows because $-1 \not\in K^2$ and the latter follows by hypothesis.
\end{proof}

\begin{corollary} \label{genbigd}
Let $K$ be a field with a non-trivial discrete valuation. Assume that $K$ has characteristic different from $2$, and that $-1 \not\in K^2$.
Then $K \in \mathcal{N}_{d,2}^\infty$ for all $d \equiv 2 \pmod{4}$.
\end{corollary}

\begin{proof}
Let $v : K^* \to \mathbb{Z}$ be a (normalized) non-trivial discrete valuation on $K$. This map is surjective, and taking $\pi$ with $v(\pi) = 1$, we also have that $v^{-1}(m)$ maps bijectively to $v^{-1}(n)$ via the map $x \mapsto \pi^{m-n}x$.  Observe that $2 \in K^*$ because $K$ has characteristic $\neq 2$, and that if $k \in K^n$ for any $n \geq 2$, then $v(k) \in n\Z$. Fix $d \equiv 2 \pmod{4}$, and let $p_1, \ldots, p_r$ be the odd prime divisors of $d$. By the Chinese Remainder Theorem, the system of congruences
$$
x \equiv  \frac{1-2v(2)}{4}  \pmod{p_i}
$$
has an infinite solution set $S \subset \Z$. If $v(k) \in S$, then $2v(2) + 4v(k) \equiv 1 \pmod{p_i}$ for all $i$, and hence $v(-4k^4) \not\in p_i\Z$, and thus $-4k^4 \not\in K^{p_i}$ for all $i$. Hence by Theorem \ref{infinited} we have that $f(x) = \left(x-4k^4\right)^d+4k^4$ has newly reducible second iterate over $K$.

It remains to argue that $v^{-1}(S)$, and hence $S$, is infinite. We have already established that $v^{-1}(m)$ and $v^{-1}(n)$ are equinumerous for any $m,n \in \Z$. The only way an infinite union of equinumerous sets can be finite is if all are empty, but this contradicts the surjectivity of $v$. Thus $v^{-1}(S)$ is infinite.
\end{proof}

\bibliographystyle{plain}


\section{Introduction}

It is a concensus that pulsar wind nebulae (PWNe) are efficient accelerators for the extremely relativistic leptons (electrons/positrons).
A pulsar located in a PWN releases its rotational energy by driving an ultra-relativistic magnetized wind composed of electron-positron pairs.
The pulsar wind interacts with the ambient medium and slows down at the so-called termination shock where the acceleration of the pulsar wind leptons occurs \citep{Gaensler2006,kargaltsev2015}.
It is generally believed that the accelerated leptons generate the electromagnetic emission from radio to very-high-energy $\gamma$-rays\ via synchrotron and inverse Compton (IC, probably including synchrotron self-Compton process) mechanism, namely the leptonic model \citep[e.g.,][]{Kennel1984a,Venter2007,Zhang2008.PWN,Tanaka2010,Bucciantini2011,Martin2012,Torres2014}.

Although the leptonic model can well explain the observational features, an open question still remains whether, or how large a fraction of, protons (or more generally, ions) may exist in pulsar wind and are further accelerated in PWN.
Theoretically, it has been suggested that protons can also be extracted from the surface of a neutron star and injected into the pulsar wind \citep{Cheng1990, Arons1994}.
Protons carrying a substantial amount of the wind energy were also predicted by the resonant cyclotron absorption model \citep{Hoshino1992,Gallant1994}.
Moreover, ion acceleration in pulsars/PWNe was explored in some recent works \citep[e.g.][]{Fang2013.PSR,Chen2014.PSR,Philippov2014,Kotera2015,Lemoine2015.PWN,Guepin2020}, suggesting that pulsars/PWNs may be the sources of Cosmic Rays.
Especially, \citet{Guepin2020} concluded that the energetic protons accelerated in the pulsar magnetosphere via the magnetic reconnection can take away the spin-down energy with a fraction of 0.2 to 4.0 per cent.

These theoretical works on proton acceleration in pulsars or PWNe indicate that the $\gamma$-ray\ emission of a PWN may contain a contribution from hadrons.
Due to the difficulties of differentiating the leptonic and hadronic origin of the $\gamma$-ray\ emission, the most direct method to test this is the neutrino experiments.
With the current stacked IceCube data toward 35 $\gamma$-ray-bright PWNe, it can only be inferred that PWNe are not the hadronic-dominated sources \citep{Aartsen2020.PWN}.
For some individual sources, however, the hadronic component was indeed invoked to explain the $\gamma$-ray\ emission, e.g. Crab \citep{Atoyan1996}, Vela X \citep{Horns2006,Zhang2009.VelaX}, PWN G54.1+0.3 \citep{Li2010.G54}, and DA 495 \citep{Coerver2019}, although the purely leptonic model can also work, especially for the Crab nebula.

The Crab nebula, powered by pulsar PSR J0534+2200, is the best studied prototype PWN.
Based on the \citeauthor{Kennel1984a} (1984) 1D flow model, \citet{Atoyan1996} calculated the Crab's spectrum extending to the $\gamma$-ray\ band and showed that the pure IC emission underpredicts the observed flux, indicating possible additional contribution from protons.
Nontheless, in a multidimensional simulation for the flow structure in the Crab, \citet{Volpi2008} found that the leptonic model can explain the EGRET data.
With the updated GeV data obtained by {\sl Fermi}\ \citep{PWN.Crab.Fermi.2012}, the broadband spectrum of the Crab was widely modeled using the leptonic model \citep[e.g.][]{Zhang2008.PWN,Tanaka2010,Bucciantini2011,Martin2012}
without any need of a proton contribution.

Recently, $>100$ TeV photons from the Crab were detected by the Tibet AS+MD\ \citep{PWN.Crab.ASg.2019} and HAWC \citep{PWN.Crab.HAWC.2019} experiments.
These new data have helped presenting the most complete spectrum across the electro-magnetic emission window, and may provide information on the proton acceleration.
In this paper, we model the broadest spectral data of the Crab nebula and explore the possibility that the very-high-energy $\gamma$-ray\ emission contains the contribution from the energetic protons.
The spectral evolution model is presented in section~\ref{sec:evol}.
In section~\ref{sec:had}, the hadronic component of $\gamma$-rays\ is constrained via the Markov Chain Monte Carlo (MCMC) method. Some discussions and conclusions are given in section~\ref{sec:dc}.

\section{The spectral evolution model}
\label{sec:evol}

\begin{figure}
\centering
\includegraphics[height=60mm]{f1.eps}
\caption{Fit of the nonthermal spectrum of Crab for the evolutionary model.
The black solid and dashed lines show the total spectrum with and without the $\gamma-\gamma$ absorption, respectively.
The data are taken from \citet{Baldwin1971} and \citet{Macias2010} for radio band, from \citet{Ney1968}, \citet{Grasdalen1979}, \citet{Green2004}, and \citet{Temim2006} for IR, from \citet{Veron1993} for optical, from \citet{Hennessy1992} and \citet{Kuiper2001.crab} for X-rays and soft $\gamma$-rays, and from \citet[][{\sl Fermi}]{PWN.Crab.Fermi.2012}, \citet[][HEGRA]{PWN.Crab.HEGRA.2004}, \citet[][{\rm H.E.S.S.}]{PWN.Crab.HESS.2006}, \citet[][MAGIC]{PWN.Crab.MAGIC.2015}, \citet[][ARGO-YBJ]{PWN.Crab.ARGO.2015}, \citet[][HAWC]{PWN.Crab.HAWC.2019} and \citet[][Tibet AS+MD]{PWN.Crab.ASg.2019} for GeV-TeV $\gamma$-rays.
}
\label{fig:sed_evol}
\end{figure}

The global multiband radiative properties of a PWN are generally described by the ``1-Zone model" \citep[see][for a review]{Bucciantini2014}.
The evolution of the lepton spectrum in the emission zone is given by the continuity equation in the energy space
\begin{equation}
\label{eq:dis_evol}
\frac{\partial N(\gamma,t)}{\partial t} = -\frac{\partial}{\partial\gamma}[\dot{\gamma}(\gamma,t)N(\gamma,t)] -\frac{N(\gamma,t)}{\tau(\gamma,t)} + Q(\gamma,t),
\end{equation}
where $\gamma$ is the Lorentz factor of the leptons, $N(\gamma,t)$ is the lepton distribution function, $\dot{\gamma}(\gamma,t)$ is the summation of the energy losses including the synchrotron radiation, the IC process, and the adiabatic loss \citep[see][]{Li2010.G54},
and $\tau(\gamma,t)$ represents the escape time which can be estimated via Bohm diffusion \citep[e.g.][]{Zhang2008.PWN}.
$Q(\gamma,t)$ represents the injection rate of leptons from shock acceleration process into the PWN emission zone per unit energy at a certain time and is generally assumed to be a broken power-law function
\begin{equation}
\label{eq:Qinj}
Q(\gamma,t)=Q_0(t)\left\{
\begin{array}{ll}
 (\gamma/\gamma_{\rm bre})^{-\alpha_1}\quad {\rm for}\ \gamma<\gamma_{\rm bre} \\
 (\gamma/\gamma_{\rm bre})^{-\alpha_2}\quad {\rm for}\ \gamma>\gamma_{\rm bre},
\end{array}
\right.
\end{equation}
where $\gamma_{\rm bre}$ denotes the break energy, and the parameters $\alpha_1$ and $\alpha_2$ are spectral indices.
The normalization is determined from the total energy in leptons ($W_{\rm e}$) extracted from the spin-down energy of a pulsar ($L_{\rm sd}$) via $W_{\rm e}=\eta_{\rm e} L_{\rm sd}=(1-\eta_{\rm B})L_{\rm sd}$, where $\eta_{\rm e}$ and $\eta_{\rm B}$ are the leptonic and magnetic energy fraction, respectively.
The maximum energy of the injected leptons is obtained by introducing a parameter $\varepsilon$, indicating that the Larmor radius of the lepton must be less than the termination shock radius \citep[see Eq.~(8) in][]{Martin2012}.

Solving equation~\ref{eq:dis_evol}, we can obtain the lepton distribution at a certain age.
Then we can calculate the multiwavelength nonthermal photon spectrum via the synchrotron and IC processes.
We apply this evolutionary model to Crab nebula and fit its observational multiband data including the radio \citep{Baldwin1971, Macias2010}, infra red \citep[IR,][]{Ney1968, Grasdalen1979, Green2004, Temim2006}, optical \citep{Veron1993}, X-ray \citep{Hennessy1992}, soft $\gamma$-ray, \citep{Kuiper2001.crab}, GeV $\gamma$-ray ({\sl Fermi}: \citealt{PWN.Crab.Fermi.2012}), and TeV $\gamma$-ray (HEGRA: \citealt{PWN.Crab.HEGRA.2004}; {\rm H.E.S.S.}: \citealt{PWN.Crab.HESS.2006}; MAGIC: \citealt{PWN.Crab.MAGIC.2015}; ARGO-YBJ: \citealt{PWN.Crab.ARGO.2015}; Tibet AS+MD: \citealt{PWN.Crab.ASg.2019}) data.
For the IC process, the seed photon fields include (i) the 2.7-K cosmic microwave background (CMB) radiation;
(ii) the excess FIR radiation with a temperature of 70 K and an energy density of 0.5 ${\rm eV\ cm^{-3}}$ \citep{Marsden1984}; (iii) the NIR radiation with a temperature of 5000 K and an energy density of 1.0 ${\rm eV\ cm^{-3}}$ \citep{Aharonian1997}; and (iv) the synchrotron radiation photon.

In the calculation, we adopt the distance $d=2.0$~kpc, the radius $R_{\rm PWN}=2.1$~pc, the initial spin-down luminosity $L_0=3.1\times10^{39}\ \ {\rm erg}\ {\rm s}^{-1}$, and the breaking index $n=2.509$ \citep{Trimble1968, Lyne1988,Taylor1993.catalog}.
The eye-fit parameters are $\alpha_1=1.55$, $\alpha_2=2.5$, $\gamma_{\rm bre}=8.5\times10^{5}$, $\eta_{\rm B}=0.012$, and $\varepsilon=0.3$, which are similar to previous results \citep[e.g.][]{Tanaka2010, Martin2012}.
The corresponding spectral energy distribution (SED) is plotted in Figure~\ref{fig:sed_evol}, which includes the $\gamma$-ray absorption based on the cross section of $\gamma$-$\gamma$ interaction \citep{Gould1967.pair} and the interstellar radiation field in Galaxy \citep{Shibata2011}.
As can be seen, the IC process can well explain the $\gamma$-ray data up to $\sim$500 TeV and confirm the leptonic origin for the broadband emission.
According to the fitted results, at the Crab's age of 950 yr, the average magnetic field strength in the emission zone and the maximum energy of injected leptons are about $100\ \mu{\rm G}$ and 4~PeV, respectively.
If protons can also be injected into the acceleration process, their maximum energy may be expected to reach similar order.
Moreover, protons with this energy can be well confined in the nebula at the magnetic field of $100\ \mu{\rm G}$ according to the Larmor formula, and may contribute to the $\ga 100$ TeV photons via p-p interaction.
In the next section, we study the probability that protons contribute some fraction of $\gamma$-ray\ emission from Crab.

\section{Hadronic component?}
\label{sec:had}
Although the broadband emissions from the Crab are explained by the leptonic process, the $\gamma$-ray emission may have hadronic component which can be constrained by the $>$100 TeV photons.
In order to quantitatively investigate this possibility, one should add the hadronic process to the above evolutionary model, and then constrain the model parameters by using a statistical method.
However, for avoiding such time-costing calculation, we can do the estimation in a simplified way.
We directly assume that the lepton distribution at the current age in the emission zone has a smooth broken power-law form with a high-energy cutoff:
\begin{equation}
\label{eq:dis_e}
\frac{dN_{\gamma}(\gamma)}{d\gamma} = A_{\rm e} \cdot \gamma^{-\alpha_1}
              \left[ 1+\left( \frac{\gamma}{\gamma_{\rm bre}} \right)^s  \right]^{\frac{\alpha_1-\alpha_2}{s}} {\rm exp}
              \left[  -\left( \frac{\gamma}{\gamma_{\rm cut}} \right)^{\beta}  \right],
\end{equation}
where $\gamma_{\rm cut}$ is the cut-off energy, $s$ the smooth parameter, and $\beta$ the cutoff shape parameter. The normalization $A_{\rm e}$ can be calculated from the total energy in leptons $W_{\rm e}$.

To check its accurateness, this distribution is firstly employed to fit the electrons in the evolutionary model with the above eye-fit parameters but for three different ages $t=500$, 950, and 5000 yr, which are shown in the black dotted lines in left panel of Figure~\ref{fig:comp}.
As it can be seen, the electrons in the evolutionary model can almost be described by the exponential cutoff smooth broken power law distribution.
\begin{figure}
\centering
\includegraphics[height=60mm]{f2.eps}
\caption{Fit of the electron distribution in evolutionary model at three ages $t=500$, 950, and 5000 yr by using Eq.~\eqref{eq:dis_e} and
comparison of the electron distribution at $t=950$ yr and that in Model A.
The factors enclosed in parenthesis are multiplied in plotting to distinguish the three ages.
}
\label{fig:comp}
\end{figure}

We now apply this simplified leptonic model to refit the broadband data of the Crab and use the MCMC approach to constrain the model parameters (Model A).
There are eight parameters in total: $\alpha_1$, $\alpha_2$, $\gamma_{\rm bre}$, $\gamma_{\rm cut}$, $B_{\rm PWN}$, $W_{\rm e}$, $s$, and $\beta$.
The best-fit parameters are listed in Table~\ref{tab:params} and the corresponding SED is plotted in Figure~\ref{fig:sed_lep}, in which the result of the evolutionary model is also displayed in magenta.
The comparison of the electron distribution in the evolution model and that in Model A is also shown in Figure~\ref{fig:comp}.
We can see that this leptonic model is good enough to explain the broadband data, especial for the $\gtrsim$ 10 TeV data which will play an important role in constraining the hadronic component.
Thus, we will use it instead of the evolutionary model to constrain the contribution of the protons.
In addition, we can directly see that the maximum energy of the leptons at current age needs to be $\sim3$~PeV to explain the data.
This hints that PWNe are one of potential sources of Galactic CRs.
\begin{figure*}
\centering
\includegraphics[height=60mm]{f3a.eps}
\includegraphics[height=57mm]{f3b.eps}
\caption{Left: same as Figure~\ref{fig:sed_evol} but for the simplified leptonic model (Model A). The gray data in IR band are not used in the MCMC fitting. Right: corresponding 1D probability distribution of the parameters. The magenta line represents the fit in Figure~\ref{fig:sed_evol}.
}
\label{fig:sed_lep}
\end{figure*}

Next, we add hadronic photons resulted from the neutral pion-decay process to the simplified leptonic model.
For the proton spectrum, we adopt the exponential cutoff power-law distribution,
\begin{equation}
\label{eq:dis_p}
\frac{dN_{\rm p}(E_{\rm p})}{dE_{\rm p}} = A_{\rm p} \cdot E_{\rm p}^{-\alpha_{\rm p}} {\rm exp}
              \left(  - \frac{E_{\rm p}}{E_{\rm c,p}} \right)
\end{equation}
where $\alpha_{\rm p}$ and $E_{\rm c,p}$ are the the proton index and the cutoff energy, respectively. The normalization $A_{\rm p}$ can be calculated from the total energy in protons $W_{\rm p}$.
For p-p process, we use the analytic photon emissivity developed by Kelner et al. (2006), including the enhancement factor of 1.84 due to contribution from heavy nuclei \citep{Mori2009}.
The density of target gas interacted by relativistic protons is not well constrained.
The mean gas density in the nebula can be estimated from the total nebular mass (gas plus dust) of $7.2 \pm 0.5 M_{\odot}$ \citep{Owen2015} and the nebular radius $R_{\rm PWN}=2.1$ pc, giving a mean density $n\sim5\ {\rm cm}^{-3}$.
Considering the fact that relativistic particles were partially captured in the dense filaments in the Crab nebula, the effective gas density for interactions of relativistic protons may be much higher than the mean density \citep{Atoyan1996}.
Here, we adopt the density of the target gas $n_{\rm t}=10\ {\rm cm}^{-3}$ as a fiducial value.
In addition, the nonthermal bremsstrahlung process can be neglected at this density level \citep[e.g.,][]{Atoyan1996}.

Due to lack of constraint on the hadronic component, we fix the proton index and the cutoff energy.
We set $E_{\rm c,p}=3$~PeV and consider two cases for the proton index:  $\alpha_{\rm p}=2.0$ (Model B, generally assumed particle index) and $\alpha_{\rm p}=1.55$ (Model C, the same as the lepton index below the break).
For this lepto-hadronic model, there are nine parameters in total: $\alpha_1$, $\alpha_2$, $\gamma_{\rm bre}$, $\gamma_{\rm cut}$, $B_{\rm PWN}$, $W_{\rm e}$, $s$, $\beta$, and $W_{\rm p}$.
The best-fit parameters are listed in Table~\ref{tab:params} and the corresponding broadband SED is plotted in Figure~\ref{fig:sed_h1}.
The fitted values of all parameters except $W_{\rm p}$ in Models B and C are almost the same as in Model A.
This means the contribution from hadronic process is negligible, which also can be seen from the result that the hadronic component (the blue line in Figure~\ref{fig:sed_h1}) is far below the observed $\gamma$-ray\ data.
In order to void exceeding the $> 100$~TeV data detected by Tibet AS+MD, it gives a constraint on the energy in protons as $W_{\rm p}\lesssim1.9\times10^{47}$ erg and $W_{\rm p}\lesssim1.0\times10^{47}$ erg at $2\sigma$ level for $\alpha_{\rm p}=2.0$ (Model B) and $\alpha_{\rm p}=1.55$ (Model C), respectively.

\begin{figure*}
\centering
\includegraphics[height=60mm]{f4a.eps}
\includegraphics[height=57mm]{f4b.eps}
\includegraphics[height=60mm]{f4c.eps}
\includegraphics[height=57mm]{f4d.eps}
\caption{Same as Figure~\ref{fig:sed_lep} but for Models B ($\alpha_{\rm p}=2.0$) and C ($\alpha_{\rm p}=1.55$), in which the hadronic components are considered.
}
\label{fig:sed_h1}
\end{figure*}

\begin{table*}
 \centering
 \begin{minipage}{160mm}
  \caption{Best-fit parameters with 1$\sigma$ errors or 95 per cent upper limits. \label{tab:params}
}
  \begin{tabular}{@{}ccccccccccc@{}}
  \hline
   Model    &  $\alpha_1$  & $\alpha_2$ & log$_{10}[\gamma_{\rm bre}$] & log$_{10}[\gamma_{\rm cut}$]
   & log$_{10}[\frac{W_{\rm e}}{\rm erg}$] & log$_{10}[\frac{B_{\rm PWN}}{\rm \mu G}$] & $s$ & $\beta$ & log$_{10}[ \frac{W_{\rm p}}{\rm erg} ]$ & $\chi^2$/dof \\

 \hline
A     & $ 1.54^{+0.01 }_{-0.01 }$ & $ 3.45^{+0.01 }_{-0.01 }$ & $ 6.03^{+0.01 }_{-0.01 }$
      & $ 9.75^{+0.01 }_{-0.01 }$ & $48.85^{+0.00 }_{-0.00 }$ & $ 2.03^{+0.00 }_{-0.00 }$
      & $ 0.88^{+0.02 }_{-0.02 }$ & $ 2.63^{+0.20 }_{-0.20 }$ & --- & 523/270 \\
B     & $ 1.54^{+0.01 }_{-0.01 }$ & $ 3.45^{+0.01 }_{-0.01 }$ & $ 6.03^{+0.01 }_{-0.01 }$
      & $ 9.75^{+0.01 }_{-0.01 }$ & $48.85^{+0.00 }_{-0.00 }$ & $ 2.03^{+0.00 }_{-0.00 }$
      & $ 0.88^{+0.02 }_{-0.02 }$ & $ 2.64^{+0.20 }_{-0.20 }$ & $ \le47.27 $ & 523/269 \\
C     & $ 1.54^{+0.01 }_{-0.01 }$ & $ 3.45^{+0.01 }_{-0.01 }$ & $ 6.03^{+0.01 }_{-0.01 }$
      & $ 9.75^{+0.01 }_{-0.01 }$ & $48.85^{+0.00 }_{-0.00 }$ & $ 2.03^{+0.00 }_{-0.00 }$
      & $ 0.88^{+0.02 }_{-0.02 }$ & $ 2.63^{+0.21 }_{-0.21 }$ & $ \le47.01 $  & 523/269 \\
D     & 1.55$^a$ & $ 3.90^{+0.08 }_{-0.08 }$ & $ 6.25^{+0.07 }_{-0.08 }$ & $\ge8.76$
      & $49.04^{+0.03 }_{-0.04 }$ & 1.8$^a$ & 0.5$^b$ & 3.0$^b$ & $48.42^{+0.15 }_{-0.04 }$ & 79/68\\
E     & 1.55$^a$ & $ 3.83^{+0.06 }_{-0.05 }$ & $ 6.18^{+0.07 }_{-0.06 }$ & $\ge8.67$
      & $49.06^{+0.02 }_{-0.02 }$ & 1.8$^a$ & 0.5$^b$ & 3.0$^b$ & $48.14^{+0.19 }_{-0.06 }$ & 80/68\\
\hline
\end{tabular}
\medskip
$^a$ Eye-fit parameters.\\
$^b$ Fixed in the MCMC fit.
\end{minipage}
\end{table*}

In the calculation, we see that the residual is mainly contributed by the radio to soft $\gamma$-ray\ data.
To bring the $\gamma$-ray\ data into play,
we will only use the GeV-TeV data (i.e., above 1 GeV) to constrain the model parameters.
In doing this, we first fit the data from radio to soft $\gamma$-ray\ band by using a smooth broken power law function, which will describe the seed photons in the SSC process (see the gray line in left panel of Figure~\ref{fig:sed_h2}).
Two cases $\alpha_{\rm p}=2.0$ (Model D) and $\alpha_{\rm p}=1.55$ (Model E) are also considered.

Just using $\gamma$-ray\ data, some of parameters do not converge very well if we do the same MCMC-fitting routine, and hence we make the following attempts.
First, the cutoff shape parameter $\beta$ is the most serious one because it is fully determined by the $\sim$100 MeV {\sl Fermi}\ data.
Considering the results in Models B and C, we round it as $\beta=3$.
In fact, due to the Klein-Nishina effect, this parameter will not obviously affect the spectral shape around 100 TeV.
Secondly, different from the narrow peak of the 1D probability in Models B and C, the smooth parameter $s$ is in a large range 0.3-0.6 and is obviously smaller that in Model B and C.
This means that the break smoothness of the leptonic spectrum is essentially determined by the IR-to-X-ray data, and the IC/$\gamma$-ray\ peak is broader than the synchrotron peak.
It may be the reason that the IC peak is broader than suggested by several previous models which is pointed out by MAGIC group \citep{PWN.Crab.MAGIC.2015,PWN.Crab.MAGIC.2020}.
The broader $\gamma$-ray\ peak not only impacts the convergence of $s$ but also $\alpha_1$, $\alpha_2$, and $\gamma_{\rm bre}$.
In order to constrain the parameters, we fix $\alpha_1=1.55$ for the index below the break and take $s=0.5$\footnote{We have also tried the case of $s=0.6$ which gives similar results.} as a example for the smooth parameter.

Then there are five parameters needed to be determined in Models D and E: $\alpha_2$, $\gamma_{\rm bre}$, $\gamma_{\rm cut}$, $W_{\rm e}$, and $W_{\rm p}$.
The best-fit parameters are listed in Table~\ref{tab:params} and the corresponding SED is plotted in Figure~\ref{fig:sed_h2}.
Different from Models B and C, the hadronic component dominates the $\gamma$-ray\ spectrum above 100~TeV (see the blue lines in left panels of Figure~\ref{fig:sed_h2}).
This implies that the total energy in protons is constrained by the Tibet AS+MD\ data as $W_{\rm p}=2.6\times10^{48}$~erg and $W_{\rm p}=1.4\times10^{48}$~erg for $\alpha_{\rm p}=2.0$ (Model D) and $\alpha_{\rm p}=1.55$ (Model E), respectively.

\begin{figure*}
\centering
\includegraphics[height=60mm]{f5a.eps}
\includegraphics[height=57mm]{f5b.eps}
\includegraphics[height=60mm]{f5c.eps}
\includegraphics[height=57mm]{f5d.eps}
\caption{Same as Figure~\ref{fig:sed_h1} but for Models D ($\alpha_{\rm p}=2.0$) and E ($\alpha_{\rm p}=1.55$) in which only $\gamma$-ray data with energy above 1.0 GeV are used in the MCMC fitting. The gray line represents the fit to data from radio to soft $\gamma$-ray, which describe the seed photons in SSC.
}
\label{fig:sed_h2}
\end{figure*}

\section{Discussion and conclusion}
\label{sec:dc}

As shown above, the broadband spectra of the Crab nebula can be fitted with the leptonic model, in which the $\gamma$-ray\ emission is dominated by the IC process.
In this picture, we also constrain the contribution to $\gamma$-rays\ from the hadronic process.
Fitting the radio to $\gamma$-ray\ data via the MCMC method (namely, Models B and C), we found
the energy in protons confined in the Crab should be less than $1.9\times10^{47}$ erg at $2\sigma$ level.
It is corresponding to the energy conversion efficiency of $\eta_{\rm p}\sim0.5\ (n_{\rm t}/10\ {\rm cm}^{-3})^{-1}$ per cent, considering the total spin-down energy of  $\sim4\times10^{49}$ erg released by Crab during its whole life and the long lifetime of the energetic protons in the Crab's environment.
This fraction is far lower than the required value in the hadronic dominated models for some PWNe in which the spin-down energy are dominantly carried by ions, e.g., Vela X \citep{Horns2006} and G54.1+0.3 \citep{Li2010.G54}.
Of course, the energy carried by protons is dependent to the target gas density $n_{\rm t}$, but it can be inferred that the hadronic emissions play a minor part in explaining the $\gamma$-ray\ emission from the Crab in Models B and C.

On the other hand, in Models D and E in which only $\gamma$-ray\ data were used alone as an extreme case, we find that the $\sim100$~TeV $\gamma$-rays\ are dominated by the hadronic process.
This requires that the energy stored in protons needs to be 1.4 -- $2.6\times10^{48}$~erg, indicating that protons steal a substantial spin-down energy with a fraction of $\eta_{\rm p}\sim7\ (n_{\rm t}/10\ {\rm cm}^{-3})^{-1}$ per cent.
This value is close to the upper limits of $<11\ (n_{\rm t}/10\ {\rm cm}^{-3})^{-1}$ per cent derived by \citet{DiPalma2017} considering the fact that IceCube has not detected any neutrinos from the Crab\footnote{It was scaled to the target density of $10 \ {\rm cm}^{-3}$ based on the material in PWN with the mass of $10\ M_{\odot}$ adopted in the original paper.}.
We also calculate the neutrino spectrum and find that the neutrino fluxes around 30 TeV are easy to exceed the IceCube's differential sensitivity as the increase of $\eta_{\rm p}\cdot n_{\rm t}$.
Our predicted neutrino fluxes at 30 TeV for $\alpha_{\rm p}=2.0$ are displayed in Figure~\ref{fig:eta}, which also shows the dependence of the neutrino flux on the product of $\eta_{\rm p}\cdot n_{\rm t}$.
It should be noted that our results are obtained on the precondition of $E_{\rm c,p}=3$~PeV which approximately equals leptons' maximum energy.
Unlike the leptons, the energy losses are negligible for protons at current age.
So their maximum energy is likely much larger than 3~PeV.
Then, in order to avoid over-estimating the flux at $\sim 100$ TeV, one should decrease the energy in protons.
Thus the conversion fraction $\eta_{\rm p}\sim7$ per cent should be treated as an upper limits.
\begin{figure}
\centering
\includegraphics[height=55mm]{f6.eps}
\caption{The predicted neutrino flux from the Crab at 30 TeV as the function of $\eta_{\rm p}\cdot n_{\rm t}$ for $\alpha_{\rm p}=2.0$ (blue). The black solid line represents the IceCube's differential sensitivity at 30 TeV adopted from \citet{Aartsen2019}. The black dashed lines are corresponding to the constraints for $\eta_{\rm p}=0.5$ and $\eta_{\rm p}=7$ per cent.
}
\label{fig:eta}
\end{figure}

Recently, \citet{Xin2019.G106} performed an analysis of the Fermi data toward the TeV source VER~J2227+608 associated with a SNR-PSR complex.
Considering the fact that the GeV-TeV $\gamma$-ray\ spectra have a hard spectral index ($1.90\pm0.04$) and no high-energy cutoff with energy upto several dozens of TeV, they suggested that the $\gamma$-ray\ emission originate from the energetic protons accelerated by a PWN not a SNR.
Based on the hadronic model, PSR J2229+6114 that is in charge of this PWN needs to convert its spin-down energy to protons with a amount of $6.0\times10^{47}(n_{\rm t}/10\ {\rm cm}^{-3})^{-1}$~erg during its life $\sim 10^{4}$~yr.
Considering the low spin-down luminosity of $2.2\times10^{37}\ {\rm erg\ s^{-1}}$ of PSR J2229+6114 \citep{Halpern2001.PSR}, the corresponding energy conversion efficiency of protons is $\eta_{\rm p}\approx10$.
This value is consistent with our results given by Models D and E, or by Models B and C if the target density in Crab is down to $\sim0.5\ {\rm cm}^{-3}$.

In addition, for Models D and E the population of leptons determined by the $\gamma$-ray\ data also can produce radio to soft-$\gamma$-ray\ emission via synchrotron process.
In order to match the radio data, we obtain $B_{\rm PWN}\approx63\ \mu{\rm G}$ (see red dashed lines in the left panel of Figure~\ref{fig:sed_h2}).
As can be seen, this population underestimates the fluxes in the IR-to-soft-gray bands.
It implies that another population of leptons are needed to fully explain the data, which is consistent with the two-component leptonic models \citep{Zhu2015.Crab,Luo2020.Crab}.

In conclusion, the broadband spectra of Crab can be explained by the leptonic model, giving a constraint on the energy conversion fraction of protons $\eta_{\rm p}<0.5$ per cent.
With this small fraction, however, it can not fully rule out that PWNe are PeVatrons of protons based on the current data.

\section*{Acknowledgements}
X.Z. is indebted to Siming Liu for the helpful discussions.
We thank the support of National Key R\&D Program of China under nos. 2018YFA0404204 and 2017YFA0402600 and NSFC grants under nos. U1931204, 11803011, 11851305, 11773014, 11633007, 11873065, 11673041, and 11533007.

\section*{Data availability}
No new data were generated or analysed in support of this research.