\section{Introduction}
Geometry and quantum mechanics are inextricably linked.
Whenever a quantum system evolves in Hilbert space, its wavefunction acquires a phase that solely depends on the path the quantum system has taken.
This concept of geometric phases was introduced by Sir Michael Berry \cite{Berry} who pointed out that there are instances in which these emerging phase factors cannot be removed by some gauge transformation.
A famous example is the Aharonov--Bohm effect \cite{Aharonov} in which the wavefunction of a charged particle travelling in a loop around a solenoidal magnetic field acquires a phase proportional to the flux through the surface enclosed by the loop, i.e. a line integral over the (Abelian) vector potential.
Here the emerging (Abelian) phase is merely a complex number.

However, any quantum system with degenerate energy levels possesses a far
richer structure as matrix-valued geometric phases \cite{Wilczek}, known as
non-Abelian holonomies. These can occur as soon as the emerging phase depends
on the order of consecutive paths. In contrast to the vector potential in the
Aharonov--Bohm scenario, a non-Abelian gauge field is responsible for the
appearance of these matrix-valued phases.
Such phases are crucial for topological quantum computation
\cite{Zanardi,Pachos}, non-Abelian anyon statistics \cite{Nayak}, and the
quantum simulation of Yang--Mills theories.
Non-Abelian synthetic gauge fields are usually realized in systems where
coupled energy levels naturally lead to the required degeneracy such as in cold
atomic samples \cite{Dalibard} and artifical atoms in superconducting circuits
\cite{Abdumalikov}.
In the case of electromagnetic fields, due to their intrinsic Abelian nature
the required degeneracy is more intricate to design.
One succesful scheme utilized the spin-orbit coupling of polarized light
in asymmetric microcavities \cite{Ma}.

In our work, we focus on a different degree of freedom and synthesize a
non-Abelian geometric phase by implementing adiabatic population transfer
\cite{Unanyan} of light. We employ an integrated photonic structure possessing
dark states, in which a non-Abelian geometric phase associated with a U(2) group
transformation is realized. The quantum metric spanned by the subspace of the
dark states induces a geodesic that defines the optimal adiabatic evolution of
these non-Abelian phases.

\section{Theory}
Gauge fields naturally arise in the context of field theories when demanding the invariance of a field under some transformation.
Invariance under multiplication by a scalar phase factor leads to the concept of Abelian gauge fields, such as the four-vector potential of electromagnetism, with commuting components.
In contrast, matrix-valued phases entail non-Abelian gauge fields where the commutator of the individual components involves the structure constants of the underlying Lie algebra.
A famous example are Yang--Mills theories of particle physics \cite{Cheng}.

A wavepacket evolving in the presence of a gauge field acquires a geometric phase.
For Abelian gauge fields, this is the famous Berry--Pancharatnam phase \cite{Berry,Pancharatnam}.
The non-Abelian generalization is known as the Wilson loop
\begin{equation}
\label{eq:wilson}
W_{\mathcal{C}}=\mathrm{Tr} \left[ \mathcal{P} \mathrm{exp} \left( -\oint_{\mathcal{C}} \mathbf{A}_{\nu} \mathrm{d}x^{\nu} \right) \right],
\end{equation}
which is the trace of the path-ordered ($\mathcal{P}$) exponential of the non-Abelian gauge field $\mathbf{A}_{\nu}$ \cite{Wilson}.
Our system is a non-trivial collection of interacting Abelian subsystems
and, thus, non-Abelian, being
characterized by a Wilson loop $W_{\mathcal{C}} < 2$~\cite{Goldman,GoldmanSpielman}.

In order to implement this concept, one seemingly requires a non-Abelian gauge field.
However, this is not necessary, as a non-Abelian structure naturally appears
whenever the evolution of a quantum system is confined to a degenerate subspace
of some Hamiltonian \cite{Wilczek}.
As a consequence, generating non-Abelian geometric phases is not connected to
the presence of a non-Abelian gauge fields, but to the existence of a
degenerate subspace.

\begin{figure}[h!]
 \center
 \includegraphics[width=8cm]{Fig1_Projection.pdf}
 \caption{\textbf{a} Scheme of the four-waveguide system.
 \textbf{b} Visualization of changing dark subspace projector along curve $\mathcal{C}$.
 The change of the dark subspace projectors of the system is given by the quantum metric \textbf{$g$}.}
\end{figure}

In our work, we consider the system sketched in Fig.~1 that consists of four potential wells $C, U, R, L$ that are coupled  with time-dependent hopping constants $Q, S, P$.
The Schr\"odinger equation for the field amplitudes $a_{\textrm{\footnotesize{L}}}, a_{\textrm{\footnotesize{U}}}, a_{\textrm{\footnotesize{R}}}, a_{\textrm{\footnotesize{C}}}$ in the individual wells reads, therefore,
\begin{equation}
\label{eq:Hamiltonian}
\mathrm{i} \, \partial_t \left( {\begin{array}{c} a_{\textrm{\footnotesize{L}}}\\a_{\textrm{\footnotesize{R}}}\\a_{\textrm{\footnotesize{U}}}\\a_{\textrm{\footnotesize{C}}}\\ \end{array}} \right)
=
\left( {\begin{array}{cccc}0 & 0 & 0 & S\\0 & 0 & 0 & P\\0 & 0 & 0 & Q\\S & P & Q & 0\\ \end{array}} \right)
\cdot
\left( {\begin{array}{c} a_{\textrm{\footnotesize{L}}}\\a_{\textrm{\footnotesize{R}}}\\a_{\textrm{\footnotesize{U}}}\\a_{\textrm{\footnotesize{C}}}\\ \end{array}} \right)
\end{equation}

This Hamiltonian supports two dark states with zero eigenvalue:
\begin{align}
\label{eq:darkstate1}
|\mathcal{D}_1\rangle &= \sin \theta |w_{\textrm{\footnotesize{L}}}\rangle - \cos \theta |w_{\textrm{\footnotesize{R}}}\rangle, \\
|\mathcal{D}_2\rangle &= \cos \theta \sin \phi |w_{\textrm{\footnotesize{L}}}\rangle + \sin \theta \sin \phi |w_{\textrm{\footnotesize{R}}}\rangle - \cos \phi |w_{\textrm{\footnotesize{U}}}\rangle,
\end{align}
where $|w_{\textrm{\footnotesize{U,R,L}}}\rangle$ are the eigenmodes of the potential wells  $U, R, L$ respectively, with the angle parameterization $\theta=\arctan(P/S)$ and $\phi= \arctan(Q/\sqrt{S^2 + P^2})$.
Notably, they do not involve the eigenstate $|w_{\textrm{\footnotesize{C}}}\rangle$ to which all other states are coupled.
These dark states span a (dark) subspace in which the adiabatic evolution of a wavefunction along a closed path can be described by a non-Abelian geometric phase (\ref{eq:wilson}) with the gauge field
\begin{equation}
\label{eq:gaugeField}
(\mathbf{A}_{\nu})_{ki} = \langle \mathcal{D}_k | \frac{\partial}{\partial x^{\nu}} | \mathcal{D}_i \rangle,
\end{equation}
written in the coordinates $\{x^{\nu}\}=\{S,P,Q\}$ (for details, see Supplementary Materials).
In the context of adiabatic evolution, the Hamiltonian (\ref{eq:Hamiltonian}) is a generalization \cite{Unanyan} of the STIRAP (\textbf{STI}mulated \textbf{R}aman \textbf{A}diabatic \textbf{P}assage) protocol \cite{STIRAP}.
What is required is to ensure adiabaticity of the evolution.

Interestingly, adiabatic transport is equivalent to parallel transport in a curved (metric) space via vanishing covariant derivative \cite{Liu}, i.e. along a geodesic defined in our parameter manifold.
The quantum metric $g_{\mu \nu} = \mathrm{Tr} \left( \partial_{\mu} P \partial_{\nu} P \right)$ is constructed from infinitesimal changes of the dark subspace projector $P(x^\nu)=\sum_i | \mathcal{D}_i (x^\nu) \rangle \langle \mathcal{D}_i (x^\nu) |$ (see Fig.~1~B).
This is the real part of a quantity known as the quantum geometric tensor
\cite{Provost,Rezakhani,Nakahara,Liu}, whose imaginary part is the field
strength tensor of the (non-Abelian) gauge field, $\mathbf{F}_{\mu\nu} =
\partial_{\mu} \mathbf{A}_{\nu} - \partial_{\nu} \mathbf{A}_{\mu} - [
\mathbf{A}_{\mu} , \mathbf{A}_{\nu} ]$.

The coordinates $x^\nu(z)$ in parameter space themselves are a function of the propagation distance $z$.
The quantum metric $g_{\mu \nu}$ defines a path length (action) along the curve $\mathcal{C}$ in the parameter manifold from the input facet $z_{\textrm{\footnotesize{i}}}$ to the output facet $z_{\textrm{\footnotesize{f}}}$,
\begin{equation}
\label{eq:Ldef}
L = \int_{z_{\textrm{\footnotesize{i}}}}^{z_{\textrm{\footnotesize{f}}}} \sqrt{g_{\mu \nu} \partial_z x^{\mu} \partial_z x^{\nu}} \mathrm{d} z.
\end{equation}
The principle of least action defines a geodesic that describes the evolution with the least diabatic error through parameter space \cite{Rezakhani}.
As a consequence, the notion of adiabaticity is intimately connected to the concept of the quantum metric.
This defines the optimal strategy for determining the time dependence of the parameters for adiabatic evolution in parameter space.
Starting from the desire to realize non-Abelian geometric phases, one first has to find a Hamiltonian with a degenerate subspace \cite{Wilczek} on which a metric can be defined.
The geodesic induced by this metric then specifies the variation of the parameters of the Hamiltonian such that the evolution through parameter space occurs with the least diabatic error.
A closed path along the geodesic in parameter space then necessarily results in a non-Abelian geometric phase.
In our experimental implementation, we minimize $L$ under the constraint of a given pulse shape, which provides the curves with the least diabatic error for a given total length $z_{\textrm{\footnotesize{f}}} - z_{\textrm{\footnotesize{i}}}$ (see Supplementary Materials).

\section{Experiment}
In order to implement our findings, we employ a photonic platform manifested in the form of integrated coupled waveguides.
Using the analogy between the quantum evolution of a wavefunction and the propagation of an optical wavepacket in the paraxial approximation \cite{Longhi}, the quantum wells in our structure can be replaced by optical waveguides.
The temporal evolution of the light amplitudes in those waveguides is governed by Eq.~(\ref{eq:Hamiltonian}) with the sole difference that the time evolution is replaced by the evolution along the waveguides described by the spatial coordinate $z$ (see Fig.~2~A).
Our design protocol yields a spatial evolution of the intersite hoppings $Q$, $S$, $P$ with an example depicted in Fig.~2~B.
The hoppings $S$ and $P$ resemble the Stokes and pump pulses of Gaussian shape in the counterintuitive sequence known from STIRAP \cite{STIRAP}, with $Q$ as an additional constant coupling.
The evolution of the parameters in parameter space is chosen to form a closed-loop trajectory as shown in Fig.~2~C.
Therefore, this evolution results in a geometric phase.

\begin{figure}
 \center
 \includegraphics[width=8cm]{Fig2_couplings.pdf}
 \caption{\textbf{a} 3D-rendering of the waveguides for one realization. \textbf{b} Coupling variation along propagating distance $z$. \textbf{c} Curve in the parameter manifold $\{S,P,Q\}$ ($Q=\textrm{const}$).}
\end{figure}

In the following, we will describe the measurement protocol for retrieving the Wilson loop. From the choice of the temporal evolution of the parameters $S$ and $P$, we have at the input facet $z_{\textrm{\footnotesize{i}}}$ and the output facet $z_{\textrm{\footnotesize{f}}}$ the relations $P(z_{\textrm{\footnotesize{i}}})/S(z_{\textrm{\footnotesize{i}}})\simeq 0\simeq S(z_{\textrm{\footnotesize{f}}})/P(z_{\textrm{\footnotesize{f}}})$.
Hence, the dark states at both facets simply become $|\mathcal{D}_1(z_{\textrm{\footnotesize{i}}})\rangle = -|w_{\textrm{\footnotesize{R}}}\rangle$, $|\mathcal{D}_2(z_{\textrm{\footnotesize{i}}})\rangle = |w_{\textrm{\footnotesize{L}}}\rangle$ and $|\mathcal{D}_1 (z_{\textrm{\footnotesize{f}}}) \rangle = |w_{\textrm{\footnotesize{L}}} \rangle$, $|\mathcal{D}_2 (z_{\textrm{\footnotesize{f}}}) \rangle = |w_{\textrm{\footnotesize{R}}} \rangle$.
As a consequence, launching light into the waveguides $L$ and $R$ excites only the dark states of the system.
Also, measuring the light intensity emanating from the waveguides $L$ and $R$ at the output facet provides the information about the population transfer between the dark states.
An initial superposition of the dark states evolves according to a unitary evolution $\mathbf{\mathcal{U}}$.
As we show in the Supplementary Materials, the elements of this unitary matrix can be expressed in terms of the amplitudes of the dark states at the input and output facet.
Moreover, the value of the Wilson loop is given by $W_{\mathcal{C}}=\mathrm{Tr}\left[\mathbf{\mathcal{U}}\right]$.
Measuring the field intensities yields the absolute values of the matrix elements $|\mathcal{U}_{ki}|$, and hence the absolute value $|W_{\mathcal{C}}|$ as shown in the Supplementary Materials.

For the fabrication of our samples, we use the femtosecond laser writing technique \cite{Szameit}.
Details of the fabrication are given in the Methods section.
We realize several structures with varying temporal profiles of the coupling parameters $S$ and $P$, resulting in different values of the Wilson loop.
An example of the evolution along the waveguides recorded by fluorescence microscopy (see Methods) is shown in Fig.~3~A.
Launching light into waveguide $L$ excites only the dark state $|\mathcal{D}_2 (z_{\textrm{\footnotesize{i}}}) \rangle$.
During the evolution, the light is coupled to the waveguides $R$ and $U$ without ever populating $C$ (see Fig.~3~B).
Hence, the evolution indeed remains in the dark subspace for all times as required for adiabaticity.

\begin{figure}
 \center
 \includegraphics[width=8cm]{Fig3_Evolution.pdf}
 \caption{Intensity distribution in the four waveguides along the propagation distance.
 \textbf{a} Experimentally measured fluorescence signal which is proportional to the intensity in the waveguides.
 The waveguide L was excited, thus the second dark state.
 The central waveguide C is located between the dashed lines, highlighting the almost vanishing intensity.
 \textbf{b} Theoretically predicted intensity from the coupled-mode theory using Eq.~(\ref{eq:Hamiltonian}).
 Red dots at the end facet are the experimentally measured intensities (compare Fig.~4).}
\end{figure}

For retrieving the elements of $\mathbf{\mathcal{U}}$, we measure the intensities at the output facet.
A representative example is shown in Fig.~4~A. In our experiments, we realized Wilson loops by implementing three different sets of parameters (details of which are given in the Supplementary Materials).
The results are summarized in Tab.~\ref{tab:results}, where the theoretical predictions and the experimental results are shown to agree well.
In all three cases, the (absolute) value of the  Wilson loop is well above $0$, thus proving the non-Abelian character of the underlying contour.

\begin{table}
\caption{\label{tab:results} \textbf{Experimental results and  corresponding theoretical Wilson loops for three different coupling variations}.
For detailed description of the used parameters and pulse shapes see the Supplementary Materials.}
\centering
\begin{tabular}{cc}
\hline\hline $|W_{\mathcal{C}}^{\textrm{\footnotesize{theo}}}|$ & $|W_{\mathcal{C}}^{\textrm{\footnotesize{exp}}}|$ \\
 \hline
 0.88 & 0.87\\
 0.97 & 1.00\\
 1.07 & 1.13\\\hline\hline
\end{tabular}
\end{table}

In order to prove that waveguide $C$ is part of the full eigenspace, we specifically launched light into $C$ and excited states in the bright subspace that extend over all waveguides (see Fig.~4~B).
From our measurements, we find that at the output facet, all waveguides are bright such that one can conclude that, indeed, waveguide $C$ has to reside within the bright subspace.

\begin{figure}
 \center
 \includegraphics[width=8cm]{Fig4_Endfacets.pdf}
 \caption{Measured intensities at the end facet. \textbf{a} Intensity distribution for dark-state excitation.
 This is for the same set of parameters as seen in Fig.~3.
 The result is clearly restricted to the dark subspace.
 \textbf{b} Intensity distribution for bright-state distribution.
 The waveguide C and U, and with that the bright states, are dominantly excited.}
\end{figure}

\section{Conclusion}
We employed evanescently coupled photonic waveguides to simulate the action of a non-Abelian gauge field on the dark subspace of the associated Hamiltonian.
The non-Abelian nature of the process was verified by measurement of the gauge invariant Wilson loop.
As the present implementation of the Wilson loop requires an adiabatic evolution within the dark subspace, the quantum metric is the appropriate tool to quantify the diabatic error.

Our results lay the foundations for the simulation of non-Abelian gauge fields using Abelian systems such as sound waves, matter waves or, in our case, light.
In particular, within this construction principle, the implementation of non-Abelian gauge fields that transform under SU($N$) could be realized with $N+2$ coupled sites containing an $N$-dimensional dark subspace.
Moreover, the use of geodesics of the quantum metric to quantify adiabaticity sheds new light on the optimization of all STIRAP-type processes.

The implementation of non-Abelian Abelian gauge fields prompts various important questions.
One of them concerns the simulation of lattice gauge field theories such as Yang--Mills theories where Wilson loops are the observable quantities.
In another context, using nonclassical light, our proposed setting is conducive to realize holonomic quantum operations as quantum logical gates can be defined as the action of non-Abelian geometric phases on the space of degenerate states, i.e. the dark subspace \cite{Pachos}.
Finally, the definition of a quantum metric induced by parametric changes of the waveguide couplings allows to study the evolution of a quantum system on curved manifolds.

\section{Acknowledgments}
The authors acknowledge funding from the Deutsche Forschungsgemeinschaft (grants SCHE 612/6-1, SZ 276/20-1, SZ 276/15-1, BL 574/13-1, SZ 276/9-2) and the Krupp von Bohlen and Halbach foundation.

\newpage


\section{I. introduction}

Spintronics relies on the transport of spins in condensed matter~\cite{Zutic,Maekawa,oxfordspin}. Spin transport has been investigated in a variety of materials, including metals, semiconductors, and insulators. In metals and semiconductors, spins are transported by the diffusion of conduction electrons~\cite{oxfordspin}. In contrast, in magnetically-ordered materials, spins can be transported even in the absence of conduction electrons; spins are carried by the elementary excitations of magnetic moments, magnons~\cite{kajiwara}. The magnonic spin current in insulators is of particular recent interest because this sets a new direction for experimental and theoretical studies of the physics of spin transport~\cite{cornelissen2015long,wesenberg2017long}.

Antiferromagnetic insulators is a new class of materials for spin transport~\cite{jungwirth2016antiferromagnetic,JUNGFLEISCH2018865,baltz2018antiferromagnetic}. This class of materials potentially entails a number of advantages as compared to ferromagnets: antiferromagnets are robust against external magnetic fields, produce no stray fields, and display ultrafast dynamics. Since the first observation of the transmission of spins through an antiferromagnetic insulator NiO~\cite{wang2014antiferromagnonic,hahn2014conduction,wang2015spin}, intense experimental and theoretical efforts have been invested in unraveling the physics of the spin transport in antiferromagnetic insulators~\cite{wang2014antiferromagnonic,hahn2014conduction,wang2015spin,doi:10.1063/1.4918990,PhysRevB.92.020409,PhysRevLett.115.266601,lin2016enhancement,takei2014superfluid,PhysRevB.93.054412,bender2017enhanced,qaiumzadeh2017spin,lebrun2018tunable,yuan2018experimental}. In antiferromagnetic insulators, the spin-decay length is known to be typically limited to only a few nanometers~\cite{baltz2018antiferromagnetic}, although theories predict long-distance spin transport in antiferromagnets~\cite{khymyn2016transformation}. This is in stark contrast to the situation for ferromagnetic insulators, where long-distance spin propagation has been observed~\cite{kajiwara,cornelissen2015long}.

In this Letter, we reveal the intrinsic character of magnonic spin transport in an antiferromagnetic insulator. We found that, in the conventional spin-injector/antiferromagnetic-insulator/spin-detector structure, the spin-transmission signal is strongly suppressed by two-magnon scattering. By eliminating the two-magnon contribution in the spin-transmission signal, we show that the spin decay length of a prototypical antiferromagnetic insulator NiO changes by two-orders of magnitude through the paramagnetic to antiferromagnetic transition. This result shows that the intrinsic spin decay length of the antiferromagnetic NiO is an order of magnitude longer than the previously believed, providing an important information for the fundamental understanding of antiferromagnetic spintronics.

\begin{figure}[tb]
	\center\includegraphics[scale=1]{Fig1.eps}
	\caption{(a) A schematic illustration of the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayer. ${\bf H}$ denotes the external magnetic field. ${\bf M}$ and $\theta_M$ represent the equilibrium direction of the magnetization when ${\bf H}$ is applied at an angle of $\theta_H$ from the film normal. (b) Magnetic field $H$ dependence of the microwave absorption signal $dI/dH$ and voltage signal $V$ for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers with $d_{\mathrm{NiO}}=0$~nm (black) and 4.1~nm (red). $I$ is the microwave absorption intensity and $H_\mathrm{res}$ is the FMR field.}
	\label{fig1}
\end{figure}

\begin{figure*}[tb]
	\center\includegraphics[scale=1]{Fig2.eps}
	\caption{(a) Out-of-plane magnetic field angle $\theta_H$ dependence of the ISHE voltage $V_\mathrm{ISHE}$ for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers with $d_\mathrm{NiO}=0$ nm, 4.1 nm, 7.2 nm, and 10.5 nm. The solid curve is the theoretical prediction based on the model of the spin pumping and ISHE, $j_\mathrm{s}(\theta_M)\sin\theta_M$. (b) $\theta_H$ dependence of the FMR field $H_{\mathrm{res}}$ for $d_{\mathrm{NiO}}=$~0 and 4.1~nm. The solid curves are the fitting results. The inset sows the $d_{\mathrm{NiO}}$ dependence of the effective demagnetization field $M_\text{eff}$. (c) Out-of-plane angle of the magnetization-precession axis, $\theta_M$, dependence of  $V_\mathrm{ISHE}$ for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers. The solid curve is a function proportional to $\sin\theta_M$. The inset shows $\theta_H$ dependence of $\theta_M$ for $d_\mathrm{NiO}=0$. }
	\label{fig2}
\end{figure*}

To quantify the intrinsic spin decay length of NiO, we prepared Ni$_{81}$Fe$_{19}$(8)/NiO($d_{\mathrm{NiO}}$)/Pt(5) trilayers on thermally oxidized Si substrates by RF magnetron sputtering at room temperature [see Fig.~\ref{fig1}(a)]. The numbers in brackets represent the thickness of each layer in nm unit, where $d_{\mathrm{NiO}}=0$ to 10.5~nm. The Ni$_{81}$Fe$_{19}$ layer, capped by 4-nm-thick SiO$_2$, is a $1\times 1.5$ mm$^2$ rectangular shape.
For the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers, we measured the spin pumping by varying a magnetic field ${\bf H}$ applied at an angle of $\theta_H$ from the film normal at room temperature [see Fig.~\ref{fig1}(a)]. The spin pumping from the Ni$_{81}$Fe$_{19}$ layer injects a spin current into the NiO layer~\cite{tserkovnyak2002spin}. The spin current reaching the Pt layer is converted into an electric voltage $V_\text{ISHE}$ through the inverse spin Hall effect (ISHE) in the Pt layer~\cite{Saitoh}, and thus the spin-current decay in the NiO layer can be characterized by measuring the $d_\mathrm{NiO}$ dependence of  $V_\text{ISHE}$. In Fig.~\ref{fig1}(b), we show the $H$ dependence of the microwave absorption intensity $I(H)$ and voltage $V(H)$ signals for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers with $d_{\mathrm{NiO}}=0$ and 4.1~nm at  $\theta_{{H}}=90^\circ$. For the measurement, the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayer was placed at the center of a TE$_{011}$ cavity with the frequency of $f=9.43$~GHz and power of $P=200$ mW, and we measured dc electric voltage $V$ between electrodes attached to the edges of the film [see Fig.~\ref{fig1}(a)].
Figure~\ref{fig1}(b) shows that the ISHE voltage $V_\text{ISHE}$ is generated around the FMR field $H=H_\mathrm{res}$. This result also shows that $V_\text{ISHE}\equiv V(H_\mathrm{res})$ is strongly suppressed by inserting the NiO layer, as expected for the spin-current decay in the antiferromagnet.

Our finding is that magnetic-field angle $\theta_H$ dependence of $V_\text{ISHE}$ strongly depends on the NiO thickness $d_\mathrm{NiO}$. In Fig.~\ref{fig2}(a), we show the $\theta_H$ dependence of $V_\text{ISHE}$ for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers with various $d_\mathrm{NiO}$. This result shows that the $\theta_H$ dependence of $V_\text{ISHE}$ for the trilayers with different $d_\mathrm{NiO}$ is the same only around $\theta_H=0$. Here, the variation of $V_\mathrm{ISHE}$ for the film with $d_\mathrm{NiO}=0$ nm is consistent with the standard model of the spin pumping and ISHE~\cite{andoJAPfull}. In this model, when the magnetic damping constant $\alpha$ is independent of $\theta_H$, the spin current generated by the spin pumping is expressed as~\cite{andoJAPfull}
 \begin{equation}
{j_\mathrm{s}}(\theta_M)=\frac{gh^2\hbar\gamma^2\omega }{4\pi\alpha^2A(\theta_M)\left[(4\pi M_\mathrm{s})^2\gamma^2\sin^4\theta_M+4\omega^2\right]},\label{spinc}
\end{equation}
where $g_\text{eff}^{\uparrow\downarrow}$ is the effective spin-mixing conductance, $h$ is the microwave magnetic field, $\gamma$ is the gyromagnetic ratio, $M_\mathrm{s}$ is the saturation magnetization, and $\omega=2\pi f$. $\theta_M$ is the out-of-plane angle of the magnetization-precession axis [see Fig.~\ref{fig1}(a)]. $A(\theta_M)={2\omega}\left[{4\pi M_\mathrm{s}\gamma \sin^2\theta_M+\sqrt{(4\pi M_\mathrm{s}\gamma)^2 \sin^4\theta_M+4\omega^2}}\right]^{-1}$ is the precession ellipticity factor. When the magnetization-precession axis is oblique to the film plane, the ISHE voltage $V_\mathrm{ISHE}$ is proportional to ${j_\mathrm{s}}(\theta_M)\sin\theta_M $ because of ${\bf j}_\mathrm{s}^\mathrm{Pt}\parallel {\bf j}_\mathrm{c}^\mathrm{Pt}\times {\bm \sigma}$~\cite{andoJAPfull}, where $j_\mathrm{s}^\mathrm{Pt}$ is the spin current density injected into the Pt layer and $j_\mathrm{c}^\mathrm{Pt}$ is the charge current density generated by the ISHE. $\bm{\sigma}$ is the spin-polarization direction of the spin current, which is parallel to the magnetization-precession axis. As shown in Fig.~\ref{fig2}(a), this model well reproduces the experimental data only for $d_\mathrm{NiO}=0$ nm [see the solid curve]. For the calculation, we determined $\theta_M$ and $M_\mathrm{s}$ from measured $\theta_H$ dependence of $H_\mathrm{res}$, shown in Figs.~\ref{fig2}(b), by solving $\omega={\gamma}\sqrt{H_XH_Y}$, where $\omega= 2\pi f$, $H_X=H_\mathrm{res}\cos(\theta_H-\theta_M)-M_{\mathrm{eff}}\cos^2\theta_M$, and $H_Y=H_\mathrm{res}\cos(\theta_H-\theta_M)-M_{\mathrm{eff}}\cos(2\theta_M)$~\cite{mizukami2002effect,lindner2009two,landeros2008two,arias1999extrinsic,arias2000extrinsic} [see the inset to Figs.~\ref{fig2}(b) and \ref{fig2}(c)]. $M_\mathrm{eff}\simeq M_\mathrm{s}$ is the effective demagnetization field.

To clarify the origin of the anomaly in the $\theta_H$ dependence of $V_\mathrm{ISHE}$ for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers with $d_\mathrm{NiO}\neq 0$ nm, we plot $\theta_M$ dependence of $V_\mathrm{ISHE}$ in Fig.~\ref{fig2}(c). Since $j_\mathrm{c}^\mathrm{Pt}(\theta_M)$ does not change drastically with $\theta_M$, $V_\text{ISHE}$ is approximately proportional to $\sin \theta_M$. In fact, the $\theta_M$ dependence of $V_\text{ISHE}$ is consistent with this scenario for the Ni$_{81}$Fe$_{19}$/Pt bilayer ($d_\mathrm{NiO}=0$ nm). However, for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers, the measured $V_\mathrm{ISHE}$ values are proportional to $\sin\theta_M$ only at $|\theta_M|<45^\circ$ as shown in Fig.~\ref{fig2}(c); $V_\mathrm{ISHE}$ deviates from $\sin\theta_M$ at $|\theta_M|>45^\circ$ with increasing the thickness of the NiO layer.

\begin{figure}[tb]
	\center\includegraphics[scale=1]{Fig3.eps}
	\caption{(a) $\theta_H$ dependence of the peak-to-peak FMR linewidth $\Delta H$ for $d_{\text{NiO}}=0$ nm and 4.1 nm. The open circles are the experimental data and the solid curves are the fitting results~\cite{supplementary}. (b) The relation between the amplitude of the two-magnon scattering $C_\mathrm{TMS}$ and $\Delta \bar{V}_\mathrm{ISHE} \equiv V_\mathrm{ISHE}^\mathrm{calc}(\theta_H=90^\circ)/V_\mathrm{ISHE}^\mathrm{calc}(\theta_H=10^\circ)-V_\mathrm{ISHE}^\mathrm{exp}(\theta_H=90^\circ)/V_\mathrm{ISHE}^\mathrm{exp}(\theta_H=10^\circ)$, where $V_\mathrm{ISHE}^\mathrm{calc(exp)}(\theta_H)$ is the calculated(measured) ISHE voltage at $\theta_H$. (c) $d_{\text{NiO}}$ dependence of the magnetic damping constant $\alpha$. (d) $d_{\text{NiO}}$ dependence of $V_{\text{ISHE}}$ at ${\theta_M = 90^\circ}$ (open circles) and $\theta_M = 40^\circ$ (solid circles). The solid line in black is the fitting result using an exponential function, $\exp(-d_\mathrm{NiO}/\lambda_\mathrm{NiO})$, for $d_\mathrm{NiO}<3$ nm. The solid lines in blue and red are the fitting result for the data at $\theta_M=90^\circ$ and $\theta_M=40^\circ$ for $d_\mathrm{NiO}>3$ nm, respectively. }
	\label{fig3}
\end{figure}

The drastic change in $V_\mathrm{ISHE}$ at $|\theta_M|>45^\circ$ indicates that the nontrivial variation of $V_\mathrm{ISHE}$ is caused by two-magnon scattering in the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers. The two-magnon scattering can be induced only when $|\theta_M|>45^\circ$ because the degenerated states with $k=0$ mode disappear at $|\theta_M|<45^\circ$~\cite{kurebayashi2013uniaxial,arias1999extrinsic,arias2000extrinsic}. Here, as shown in Fig.~\ref{fig1}(b), the peak-to-peak FMR linewidth $\Delta H$ is clearly enhanced by inserting the NiO layer, despite the negligible change in the effective demagnetization field $M_\mathrm{eff}$ [see the inset to Fig.~\ref{fig2}(b)]. To quantitatively study the damping enhancement induced by the NiO insertion, we plot $\theta_H$ dependence of $\Delta H$ in Fig.~\ref{fig3}(a). Figure~\ref{fig3}(a) shows $\Delta H(\theta_H=\theta_M=0^\circ)\simeq \Delta H(\theta_H=\theta_M=90^\circ)$ for $d_\mathrm{NiO}=0$ nm, while $\Delta H(\theta_H=\theta_M=0^\circ)<\Delta H(\theta_H=\theta_M=90^\circ)$ for $d_\mathrm{NiO}=4.1$ nm. This result indicates that $\Delta H(\theta_H=90^\circ)$ for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayer is influenced by the two-magnon scattering.

The two-magnon scattering is known to be activated by the random fluctuation of uniaxial anisotropy, surface/interface roughness, and defects~\cite{azevedo2000extrinsic,arias1999extrinsic,arias2000extrinsic,landeros2008two,lindner2009two}. We note that in the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers, the NiO layer is polycrystalline, as evidenced by the X-ray diffractometry~\cite{supplementary}. This suggests that the two-magnon scattering can be induced by the random fluctuation of uniaxial anisotropy due to randomly oriented exchange bias fields~\cite{sakimura2018}.
In fact, the measured $\theta_H$ dependence of $\Delta H$ is well reproduced by a calculation which takes into account the additional damping due to the two-magnon scattering as shown in Fig.~\ref{fig3}(a)~\cite{sakimura2018,lindner2009two} [for details, see \cite{supplementary}]. In the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayers, the random fluctuation of uniaxial anisotropy due to the randomly oriented exchange bias increases with $d_\mathrm{NiO}$~\cite{sakimura2018}; although the surface roughness of the NiO layer is almost unchanged with $d_\mathrm{NiO}$~\cite{supplementary}, the amplitude of the two-magnon scattering $C_\mathrm{TMS}$ increases with $d_\mathrm{NiO}$, which is reminiscent of the increased suppression of $V_\mathrm{ISHE}$ with $d_\mathrm{NiO}$ shown in Fig.~\ref{fig2}(c).
Here, we characterize the suppression of $V_\mathrm{ISHE}$ induced by the NiO insertion as the difference between the measured $V_\mathrm{ISHE}$ and $V_\mathrm{ISHE}$ calculated using the conventional spin-pumping model, $\Delta \bar{V}_\mathrm{ISHE} \equiv V_\mathrm{ISHE}^\mathrm{calc}(\theta_H=90^\circ)/V_\mathrm{ISHE}^\mathrm{calc}(\theta_H=10^\circ)-V_\mathrm{ISHE}^\mathrm{exp}(\theta_H=90^\circ)/V_\mathrm{ISHE}^\mathrm{exp}(\theta_H=10^\circ)$, where $V_\mathrm{ISHE}^\mathrm{calc}(\theta_H)$ and $V_\mathrm{ISHE}^\mathrm{exp}(\theta_H)$ are the calculated and measured ISHE voltage at $\theta_H$, respectively [see Fig.~\ref{fig2}(a)].
To clarify the relation between $C_\mathrm{TMS}$ and the voltage suppression, we plot $\Delta \bar{V}_\mathrm{ISHE} $ with respect to $C_\mathrm{TMS}$, extracted by the calculation shown in Fig.~\ref{fig3}(a). As shown in Fig.~\ref{fig3}(b), $\Delta \bar{V}_\mathrm{ISHE} $ increases with $C_\mathrm{TMS}$, supporting that the suppressed $V_\mathrm{ISHE}$ signals at $|\theta_M|>45^\circ$ is caused by the two-magnon scattering.

From the calculation of the $\theta_H$ dependence of $\Delta H$, we also extracted the damping constant $\alpha= \mu _ { 0 } [\Delta H-(\Delta H_\mathrm{inhomo}+\Delta H_\mathrm{TMS})] ( { \sqrt { 3 }  }/ { 2})  ({  \gamma \Xi} /{ \omega })$, where $\Delta H_\mathrm{inhomo}$ and $\Delta H_\mathrm{TMS}$ are the linewidth due to inhomogeneity and two-magnon scattering, respectively. $\Xi$ is the dragging function~\cite{supplementary}. Figure~\ref{fig3}(c) shows that $\alpha$ decreases at $d_{\text{NiO}}=2.0$~nm, while $\alpha$ increases above $d_{\text{NiO}}=4.1$~nm, consistent with previous reports~\cite{wang2014antiferromagnonic, wang2015spin}; $\alpha$ decreases due to the decoupling of the Ni$_{81}$Fe$_{19}$ and Pt layers by the insulating and non-N\'eel-ordered NiO layer because the N\'eel temperature of 2-nm-thick NiO is below the room temperature~\cite{gruyters2002structural,baruth2008enhanced,wang2014spin}. Above $d_{\text{NiO}}=4.1$~nm, $\alpha$ increases because of the enhanced antiferromagnetic correlation due to the thickness growth~\cite{wang2014antiferromagnonic,lang2007dependence}.

Commonly, the spin decay length $\lambda_\mathrm{NiO}$ of NiO is obtained from the thickness $d_\mathrm{NiO}$ dependence of $V_\mathrm{ISHE}$ at $\theta_H=\theta_M=90^\circ$~\cite{wang2014antiferromagnonic, wang2015spin,hahn2014conduction}. Following this procedure, we plot the $d_\mathrm{NiO}$ dependence of $V_\mathrm{ISHE}$ at $\theta_M=90^\circ$ in Fig.~\ref{fig3}(d). This result shows that the spin decay length is increased from  $\lambda_\mathrm{NiO}=1.8$ nm for $d_\mathrm{NiO}<3$ nm to $\lambda_\mathrm{NiO}=8.8$ nm for $d_\mathrm{NiO}>3$ nm. The increase of $\lambda_\mathrm{NiO}$ can be attributed to the paramagnetic to antiferromagnetic transition; for $d_\mathrm{NiO}<3$~nm, the N\'eel temperature is lower than the room temperature, while the NiO layer with $d_{\text{NiO}}>3$~nm is antiferromagnetic at room temperature ~\cite{wang2014antiferromagnonic,baruth2008enhanced,gruyters2002structural}.
$\lambda_\mathrm{NiO}=8.8$ nm in the antiferromagnetic state is consistent with previous reports~\cite{wang2014antiferromagnonic,wang2015spin}. However, we note that, as is clear from Fig.~\ref{fig2}(a), the $V_\mathrm{ISHE}$ signals at $\theta_H=90^\circ$ are strongly suppressed by the two-magnon scattering. This results in under estimation of the spin decay length in the antiferromagnetic state because the voltage suppression increases with $d_\mathrm{NiO}$.

The intrinsic spin decay length, where the two-magnon contribution is excluded, can be determined only from the $d_\mathrm{NiO}$ dependence of $V_\mathrm{ISHE}$ at $|\theta_M|<45^\circ$, where the voltage suppression due to the two-magnon scattering is absent. As shown in Fig.~\ref{fig3}(d), the $d_\mathrm{NiO}$ dependence of $V_\mathrm{ISHE}$ at $\theta_M=40^\circ$ is clearly different from that at
$\theta_M=90^\circ$. From the data at $\theta_M=40^\circ$, for the antiferromagnetic NiO, we obtain $ \lambda_\mathrm{NiO}=109$ nm, which is almost ten times longer than previously reported values~\cite{wang2014antiferromagnonic,wang2015spin}. We also note that the characteristic length of spin decay in NiO changes by two-orders of magnitude through the paramagnetic to antiferromagnetic transition, illustrating the crucial role of the antiferromagnetic order for efficient spin transport in antiferromagnetic insulators.

In summary, we investigated magnonic spin transport in an antiferromagnetic insulator NiO. We found that in the in-plane magnetic field geometry, the spin transport signal is strongly suppressed by the two-magnon scattering. By changing the magnetic-field angle, the two-magnon scattering contribution can be eliminated, which enables to determine the intrinsic spin decay length of the antiferromagnetic insulator.
Although the spin transport signal for the Ni$_{81}$Fe$_{19}$/NiO/Pt trilayer with much thicker $d_\mathrm{NiO}$ is difficult to measure because the surface roughness of the NiO layer increases with $d_\mathrm{NiO}$, our result shows that the intrinsic spin decay length of the prototypical antiferromagnetic insulator NiO is longer than 100 nm, which is an order of magnitude longer than previously believed. The result shows that the spin decay length changes by two-orders of magnitude through the paramagnetic to antiferromagnetic transition. Our results therefore demonstrate the crucial role of the antiferromagnetic order for efficient spin transport in antiferromagnetic insulators, as well as the two-magnon scattering in quantifying the spin transport in antiferromagnets.

\begin{acknowledgments}
This work was supported by JSPS KAKENHI Grant Numbers 26220604, 26103004, the Asahi Glass Foundation, and JGC-S Scholarship Foundation. H.S. is supported by JSPS Grant-inAid for Research Fellowship for Young Scientists (DC1) No. JP17J03624.
\end{acknowledgments}


\section*{Introduction}

Let $\cal{O}_K$ be an excellent Henselian DVR with perfect residue field $k$ and fraction field $K$. Let $K^\mathrm{s}$ be a separable closure of $K$ and let $G_K=\mathrm{Gal}(K^{\rm s}/K)$ be the corresponding absolute Galois group. Then, the residue field $\bar k$ of $K^{\rm s}$ is an algebraic closure of $k$ and we let $G_k$ denote the corresponding absolute Galois group. Let $X$ be a K3 surface over $K$.

\begin{definitionu} A model for $X$ over $\mathcal{O}_K$ is a flat and proper morphism $\mathcal{X}\rightarrow \spec{\mathcal{O}_K}$ of algebraic spaces together with an isomorphism $X\overset{\cong}{\longrightarrow} \mathcal{X}_K$.
\end{definitionu}

We say that $X$ has good reduction if it admits a smooth model over $\mathcal{O}_K$. It follows from the smooth and proper base change theorem that if $X$ has good reduction, then the $G_K$-action on $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified for any prime $\ell\neq \mathrm{char}(k)$. In order to study the extent to which the converse of this is true, we consider the following (strong) notion of semistable reduction.

\begin{AssumptionStar}
  A K3 surface $X$ over $K$ satisfies ($\star$) if
  there exists a finite field extension $L/K$ and a model $\mathcal{X}$ for $X_L$ over $\mathcal{O}_K$ such that:
  \begin{enumerate}
  \item the special fibre of $\mathcal{X}$ is a strict normal crossing divisor;
  \item the relative canonical divisor $\omega_{\mathcal{X}/\mathcal{O}_L}$ is trivial.
  \end{enumerate}
\end{AssumptionStar}

In honour of Kulikov's results \cite{Kul77}, such a model $\mathcal{X}$ is called a \emph{Kulikov} model.
In equicharacteristic zero, Assumption $(\star)$ is always
satisfied \cite{Kul77, PP81}. Also, if the residue characteristic $p$ is at least $5$, Assumption $(\star)$ is satisfied whenever $X$ has potential scheme theoretic semistable reduction \cite{LM18, Mau14}. For example, this is satisfied if $X$ admits an ample invertible sheaf $\cal{L}$, whose self-intersection number
satisfies $\cal{L}^2+4<p$. We refer to \cite[\S3]{LM18} for further results and discussion.

\subsection*{Good reduction in characteristic zero}
If $\mathrm{char}(K)=0$, then converse results, that is, N\'eron-Ogg-Shafarevich criteria for good reductions
of K3 surfaces have been established in \cite{CLL17, LM18, Mat15}. More precisely, the main result (Theorem 1.3) of \cite{LM18} reads as follows.

\begin{theoremn}[Liedtke--Matsumoto] \label{theo: LM18 main} Assume that $\mathrm{char}(K)=0$, let $X$ be a K3 surface over $K$, satisfying Assumption $(\star)$, and let $\ell$ be a prime $\neq \mathrm{char}(k)$. Consider the following three statements:
\begin{enumerate}
\item the $G_K$-action on $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified; \label{cond: unram}
\item $X$ has good reduction over a finite and unramified extension $L/K$;
\item $X$ admits a model $\mathcal{X}$ that is projective over $\mathcal{O}_K$ and whose special fibre is a K3 surface with at worst RDP singularities.
\end{enumerate}
Then, we have $(1)\Leftrightarrow(2)\Rightarrow (3)$.
\end{theoremn}

\begin{remarku} There exist examples of K3 surfaces over ${\mathbb Q}_p$ with $p\geq5$ that only admit good reduction over a non-trivial unramified extension, see \cite[Theorem 1.6]{LM18}. We will show how to modify these to produce similar examples over Laurent series fields ${\mathbb F}_p\lser{t}$  in \S\ref{sec: counter} below.
\end{remarku}

This begs the question of whether there exists a criterion that sees good reduction `on the nose'. This was the question investigated in \cite{CLL17} and to explain the answer, we note that any K3 surface $X/K$ satisfying $(1)$, $(2)$ and $(3)$ in Theorem \ref{theo: LM18 main} has a `canonical reduction' $Y$, which is a K3 surface over $k$. Concretely, this is given as the minimal resolution of singularities of the special fibre of any RDP model $\mathcal{X}$ as in $(3)$.
We note that $Y$ is unique up to canonical isomorphism and does not depend on the choice of model. The main result (Theorem 1.6) of \cite{CLL17} is then the following.

\begin{theoremn}[Chiarellotto--Lazda--Liedtke] \label{theo: CLL17 main}
Assume that $\mathrm{char}(K)=0$, let $X$ be a K3 surface over $K$, satisfying Assumption $(\star)$ and let $\ell$ be a prime $\neq \mathrm{char}(k)$. Then, then following are equivalent:
\begin{enumerate}
\item $X$ has good reduction over $K$;
\item \label{cond: NOS} the $G_K$-representation $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified, and there is an isomorphism
\[ H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell) \overset{\cong}{\longrightarrow} H^2_\mathrm{\acute{e}t}(Y_{\bar k},{\mathbb Q}_\ell)\]
of $G_k$-representations.
\end{enumerate}
\end{theoremn}

\begin{remarku} When $\mathrm{char}(k)=p>0$, there are also $p$-adic versions of both results in terms of crystalline cohomology, see \cite[Theorem 1.1 and Theorem 1.6]{CLL17}.\end{remarku}

\subsection*{Good reduction in positive equicharacteristic}
The classical $\ell$-adic good reduction criterion of Serre--Tate for abelian varieties \cite{ST68} makes no restriction on the generic characteristic of the base. Neither does its $p$-adic analogue \cite[Exp. IX, Thm. 5.13]{SGA7i}, phrased in terms of the $p$-divisible group associated to an abelian variety. The purpose of this note is therefore to explain how to extend the above results on good reduction of K3 surfaces to the case when $\mathrm{char}(K)=p>0$, both in terms of $\ell$-adic and $p$-adic cohomology. As we already observed in \cite{CLL17}, essentially everything in the $\ell$-adic case from \cite{LM18,CLL17} goes through if Assumption ($\star$) is strengthened to require the given finite extension $L/K$ to be \emph{separable} - this is Assumption ($\star^{\rm s}$) in \S\ref{sec: recap}. The key point, then, is to deal with the problem of descending good reduction over a finite and purely inseparable extension $L/K$. Our main result in this direction is the following.

\begin{theoremu}[Corollary \ref{cor: gr desc}] Let $L/K$ be a finite and purely inseparable extension and $X/K$ a K3 surface. Then, $X$ has good reduction over $K$ if and only if it has good reduction over $L$.
\end{theoremu}

In fact, the main property of K3 surfaces needed is that they do not admit non-zero global vector fields.
Thus, we actually prove a more general statement about purely inseparable descent of models (Theorem \ref{thm: gr desc}) that may be of independent interest.
This result already gives the desired extensions of the above results to $\ell$-adic cohomology.

\begin{corollaryu} Theorems \ref{theo: LM18 main} and \ref{theo: CLL17 main} hold without any restriction on $\mathrm{char}(K)$.
\end{corollaryu}

In order to provide $p$-adic versions of these results in equicharacteristic $p>0$, we need to use the Robba ring valued version of rigid cohomology, as was developped in reasonable generality in \cite{LP16}.\footnote{That this cohomology theory can be used to rephrase the `classical' $p$-adic criterion of Grothendieck for abelian varieties is explained in \cite[\S5.3]{LP16}} Let $W=W(k)$ denote the ring of Witt vectors of $k$, write $\kappa=W[1/p]$ for its fraction field, and $\mathcal{R}$ for the Robba ring over $\kappa$. Then for any variety $X/K$, there are rigid cohomology groups $H^i_\mathrm{rig}(X/\mathcal{R})$,
which are $(\varphi,\nabla)$-modules over $\mathcal{R}$, and behave completely analogously to the $\ell$-adic \'etale cohomology groups for $\ell\neq p$. In fact, since we will only be working with smooth and proper varieties it is possible to phrase everything purely in terms of crystalline cohomology, see \S\ref{sec: robba} for a detailed discussion. In the $p$-adic case in equicharacteristic, the analogue of Theorem \ref{theo: LM18 main}
is the following.

\begin{theoremn} \label{theo: LM18 main equi} Assume that $\mathrm{char}(K)=p>0$, and let $X$ be a K3 surface over $K$, satisfying Assumption $(\star)$. Consider the following statements:
\begin{enumerate}
\item for some prime (resp. all primes) $\ell\neq p$, the $G_K$-action on $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified;
\item the $(\varphi,\nabla)$-module $H^2_\mathrm{rig}(X/\mathcal{R})$ admits a basis of horizontal sections;
\item $X$ has good reduction over a finite, unramified extension $L/K$;
\item $X$ admits a model $\mathcal{X}$ that is projective over $\mathcal{O}_K$, and whose special fibre is a K3 surface with at worst RDP singularities.
\end{enumerate}
Then, we have $(1)\Leftrightarrow(2)\Leftrightarrow (3) \Rightarrow (4)$.
\end{theoremn}

Similarly, the analogue of Theorem \ref{theo: CLL17 main} is the following.

\begin{theoremn} \label{theo: CLL17 main equi}
Assume that $\mathrm{char}(K)=p>0$, and let $X$ be a K3 surface over $K$, satisfying Assumption $(\star)$. Then then following are equivalent:
\begin{enumerate}
\item $X$ has good reduction over $K$;
\item for some prime (resp. all primes) $\ell \neq p$, the $G_K$-representation $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified, and there is an isomorphism
\[ H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell) \overset{\cong}{\longrightarrow} H^2_\mathrm{\acute{e}t}(Y_{\bar k},{\mathbb Q}_\ell)\]
of $G_k$-representations.
\item the $(\varphi,\nabla)$-module $H^2_\mathrm{rig}(X/\mathcal{R})$ admits a basis of horizontal sections, and there is an isomorphism
\[ H^2_\mathrm{rig}(X/\mathcal{R})^{\nabla=0} \overset{\cong}{\longrightarrow} H^2_\mathrm{rig}(Y/\kappa) \]
of $F$-isocrystals over $\kappa$.
\end{enumerate}
\end{theoremn}

Finally, there also exist examples similar to the ones of \cite[Theorem 1.6]{LM18} over Laurent series fields
that show that in general one \emph{cannot} expect to choose the finite and unramified extension $L/K$ from (3)
of Theorem \ref{theo: LM18 main equi} to be trivial.

\begin{theoremu}[\ref{theo: counter example}]
For every prime $p\geq5$, there exists a smooth K3 surface over ${\mathbb F}_p\lser{t}$
  that admits good reduction over ${\mathbb F}_{p^2}\lser{t}$ but not over ${\mathbb F}_p\lser{t}$.
\end{theoremu}

\subsection*{Organisation}
Let us now give a brief summary of the contents of this note. In \S\ref{sec: robba} we will give a brief overview of $p$-adic cohomology over positive characteristic local fields. In \S\ref{sec: recap} we will roughly outline the arguments from \cite{LM18} and \cite{CLL17}, and explain why they carry over to the equicharacteristic $p>0$ case under a suitable strengthening of Assumption $(\star)$ that we call Assumption ($\star^{\rm s}$). In \S\ref{sec: descent} we will discuss descent under purely inseparable fpqc covers, and prove our main result that smooth models of K3 surfaces descend under purely inseparable field extensions. In \S\ref{sec: together} we will bring everything together to prove our main results, Theorems \ref{theo: LM18 main equi} and \ref{theo: CLL17 main equi}. Finally, in \S\ref{sec: counter} we will construct counter-examples asserted by Theorem \ref{theo: counter example}.

\subsection*{Acknowledgements} We would like to thank Yuya Matsumoto for explaining to us why arithmetic threefold flops still exist in equicharacteristic $2$. The second named author is support by the Netherlands Organisation for Scientific Research (NWO). The third named author is supported by the ERC Consolidator Grant 681838 ``K3CRYSTAL''.

\section{Review of \texorpdfstring{$p$}{p}-adic cohomology in equicharacteristic \texorpdfstring{$p>0$}{pg0}}\label{sec: robba}

In this section we will briefly review the facts needed on $p$-adic cohomology when $\mathrm{char}(K)=p>0$.
Thus, $\cal{O}_K$ is an excellent and Henselian DVR of characteristic $p>0$ with fraction field $K$ and perfect residue field $k$. We let $W=W(k)$ denote the ring of Witt vectors of $k$, $\kappa$ its fraction field, and $\sigma$ the $p$-power Frobenius on $W$ and $\kappa$.

\begin{remark} The choice of $\kappa$ for the fraction field of $W(k)$ is distinctly ``non-standard''. However, the more usual choices of $K_0$ or $K$ are unfortunately excluded by the use of $K$ for our given equicharacteristic local field.
\end{remark}

In this situation, a choice of uniformiser $\varpi\in\cal{O}_K$ induces an isomorphism $\widehat{K}\cong k\lser{\varpi}$ between the completion of $K$ and a Laurent series field over $k$. We will denote by $\mathcal{R}$ the Robba ring over $\kappa$, that is, the ring consisting of those Laurent series $\sum_i a_it^i$ with $a_i\in\kappa$
such that:
\begin{itemize} \item for all $\rho<1$, $\norm{a_i}\rho^i\rightarrow 0$ as $i\rightarrow \infty$;
\item for some $\eta<1$, $\norm{a_i}\eta^i\rightarrow 0$ as $i\rightarrow -\infty$.
\end{itemize}
In other words, $\mathcal{R}$ is the ring of functions convergent on some semi-open annulus $\eta\leq \norm{t}<1$. The ring of integral elements $\mathcal{R}^{\mathrm{int}}$ (that is, those with $a_i\in W$) is therefore a lift of $\widehat{K}$ to characteristic $0$, in the sense that mapping $t\mapsto \varpi$ induces an isomorphism $\mathcal{R}^{\mathrm{int}}/(p)\cong \widehat{K}$. We will denote by $\sigma$ a Frobenius on $\mathcal{R}$, that is, a continuous $\sigma$-linear endomorphism preserving $\mathcal{R}^{\mathrm{int}}$ and lifting the absolute $p$-power Frobenius on $\widehat{K}$.
We will moreover assume that $\sigma(t)=ut^p$ for some $u\in (W\pow{t}\otimes_W \kappa)^\times$. The reader is welcome to assume that $\sigma(\sum_i a_it^i)=\sum_i \sigma(a_i)t^{ip}$. Let $\partial_t:\mathcal{R}\rightarrow \mathcal{R}$ denote the derivation given by differentiation with respect to $t$.

\begin{definition} \label{pnm1} A $(\varphi,\nabla)$-module over $\mathcal{R}$ is a finite free $\mathcal{R}$-module $M$ together with:
\begin{itemize} \item a connection, that is, a $\kappa$-linear map $\nabla:M\rightarrow M$ such that
\[ \nabla(rm)=\partial_t(r)m+r\nabla(m)\;\;\;\;\text{for all}\;\; r\in \mathcal{R} \;\; \text{and} \;\;m\in M;\]
\item a horizontal Frobenius $\varphi:\sigma^*M:=M\otimes_{\mathcal{R},\sigma} \mathcal{R} \overset{\sim}{\rightarrow} M$.
\end{itemize}
\end{definition}

Then, $(\varphi,\nabla)$-modules over $\mathcal{R}$ should be considered as $p$-adic analogues of Galois representations. For example, they satisfy a local monodromy theorem, see \cite{Ked04a}. More specifically, the connection $\nabla$ should be viewed as an analogue of the action of the inertia subgroup $I_K$ and the Frobenius $\varphi$ the action of some Frobenius lift in $G_K$. The analogue for $(\varphi,\nabla)$-modules of a Galois representation being unramified is therefore the connection acting trivially, or, in other words, the $(\varphi,\nabla)$-module admitting a basis of horizontal sections.

The $p$-adic completion $\widehat{\mathcal{R}}^\mathrm{int}$ of the integral Robba ring is a Cohen ring for $\widehat{K}$ and hence, any smooth and proper $\widehat{K}$-variety $Y$ has crystalline cohomology groups
\[ H^i_\mathrm{cris}(Y/\widehat{\mathcal{R}}^\mathrm{int}), \]
which are $(\varphi,\nabla)$-modules over $\widehat{\mathcal{R}}^\mathrm{int}$, see for example \cite[\S.3.1]{GM87}. It follows from \cite[Theorem 7.0.1]{Ked00} (where the notations $\Gamma= \widehat{\mathcal{R}}^\mathrm{int}$ and $\Gamma^\dagger=\mathcal{R}^\mathrm{int}$ are used) that the crystalline cohomology groups of any smooth and proper variety $Y/\widehat{K}$ descend uniquely to $(\varphi,\nabla)$-modules
\[ H^i_\mathrm{cris}(Y/\mathcal{R}^\mathrm{int}) \]
over $\mathcal{R}^\mathrm{int}$. Therefore, we may define for any smooth and proper variety $X/K$
\[ H^i_\mathrm{rig}(X/\mathcal{R}) := H^i_\mathrm{cris}(X_{\widehat{K}}/\mathcal{R}^\mathrm{int})\otimes_{\mathcal{R}^\mathrm{int}} \mathcal{R} \]
as $(\varphi,\nabla)$-modules over $\mathcal{R}$. We have the following analogue of the smooth and proper base change theorem.

\begin{proposition} \label{prop: spbc} Let $\mathcal{X}\rightarrow \spec{\mathcal{O}_K}$ be a smooth and proper morphism of algebraic spaces, whose generic fibre is a scheme. Then, for any $i\geq 0$ there exists a canonical isomorphism
\[ H^i_\mathrm{rig}(\mathcal{X}_K/\mathcal{R})^{\nabla=0} \overset{\cong}{\longrightarrow} H^i_\mathrm{rig}(\mathcal{X}_k/\kappa), \]
of $F$-isocrystals over $K$.
\end{proposition}

\begin{proof}
Since $\mathcal{X}_k$ is smooth over $k$, the right hand side is isomorphic to the rational log-crystalline cohomology of $\mathcal{X}_k$, equipped with the log structure
\[ {\mathbb N} \rightarrow \mathcal{O}_{\mathcal{X}_k} ,\;\;\;\;1\mapsto 0 \]
so we may apply \cite[Proposition 2.3]{CL16}.
\end{proof}

We will also need cycle class maps in $p$-adic cohomology. For any smooth and proper variety $Y/\widehat{K}$, homomorphisms
\[ [-]:\mathrm{CH}^d(Y) \rightarrow H^{2d}_\mathrm{cris}(Y/\widehat{\mathcal{R}}^\mathrm{int})_{\mathbb Q} \]
are constructed in \cite[Theorem 4.3.1]{GM87}, and since their image lands in the subspace of horizontal sections on which Frobenius acts as $p^d$, it follows from Kedlaya's full faithfullness theorem \cite[Theorem 5.1]{Ked04c} that for any smooth and proper variety $Y/\widehat{K}$, the codimension-$d$ crystalline cycle class map for $Y$ actually takes values in $H^{2d}_\mathrm{cris}(Y/\mathcal{R}^\mathrm{int})_{{\mathbb Q}}^{\nabla=0,\varphi=p^d}$. This allows us to define, for any smooth and proper variety $X/K$, cycle class maps
\[ [-]:\mathrm{CH}^d(X) \rightarrow H^{2d}_\mathrm{rig}(X/\mathcal{R}) \]
landing in the subspace of horizontal sections on which Frobenius acts as $p^d$. On the other hand, for any smooth and proper variety $Y/k$, we also have crystalline cycle class map
\[ [-]:\mathrm{CH}^d(Y) \rightarrow H^{2d}_\mathrm{rig}(Y/\kappa). \]
The following is the $p$-adic analogue of \cite[Lemma 5.6]{LM18}, see also \cite[\S1]{CCM13}.

\begin{proposition} \label{prop: cyccomp} Let $\mathcal{X}\rightarrow \spec{\mathcal{O}_K}$ be a smooth and proper morphism of algebraic spaces, whose special and generic fibres are schemes. Then, for any closed subspace $\mathcal{Z}\hookrightarrow \mathcal{X}$ of constant codimension $d$ and flat over $\mathcal{O}_K$, the isomorphism
\[  H^{2d}_\mathrm{rig}(\mathcal{X}_K/\mathcal{R})^{\nabla=0} \overset{\cong}{\longrightarrow} H^{2d}_\mathrm{rig}(\mathcal{X}_k/\kappa) \]
sends $[\mathcal{Z}_K]$ to $[\mathcal{Z}_k]$.
\end{proposition}

\begin{proof} We may assume that $K=\widehat{K}$ is complete. By combining Chow's lemma with de\thinspace Jong's theorem on alterations \cite{dJ96} we may choose a finite extension $L/K$, a projective and strictly semistable scheme $\mathcal{Y}\rightarrow \spec{\mathcal{O}_L}$ and an alteration $\pi: \mathcal{Y}\rightarrow \mathcal{X}$ of algebraic spaces over $\mathcal{O}_K$.

Let $k_L$ denote the residue field of $L$, $W_L=W(k_L)$ and fix a lift of $\mathcal{O}_K\rightarrow \mathcal{O}_L$ of the form $W\pow{t}\rightarrow W_L\pow{t_L}$. Let $\mathcal{R}_L$ denote a copy of the Robba ring over $\kappa_L=W_L[1/p]$ with parameter $t_L$, and $\mathcal{R}\rightarrow \mathcal{R}_L$ the induced finite flat extension. This sends $\mathcal{R}^\mathrm{int}$ into $\mathcal{R}_L^\mathrm{int}$. Let $\mathcal{Y}^\times$ denote the scheme $\cal{Y}$ equipped with the fs log structure induced by the special fibre, and $\mathcal{Y}_{k_L}^\times$ the special fibre of this log scheme. Let $W_{k_L}^\times$ denote the ring $W_{k_L}$ equipped with the log structure \begin{align*}
{\mathbb N} &\rightarrow W_{k_L} \\
1 &\mapsto 0
\end{align*}
and $W_L\pow{t_L}^\times$ for $W_L\pow{t_L}$ equipped with the log structure
\begin{align*}
{\mathbb N} &\rightarrow W_{k_L}\pow{t_L} \\
1 &\mapsto t_L.
\end{align*}
By \cite[Theorem 5.46]{LP16} there is, for all $i\geq0$, an isomorphism
\[ H^i_\mathrm{rig}(\mathcal{Y}_L/\mathcal{R}_L)^{\nabla=0} \overset{\cong}{\longrightarrow} H^i_{\log\text{-}\mathrm{cris}}(\mathcal{Y}_{k_L}^\times/W_L^\times)\otimes {\mathbb Q}, \]
fitting into a commutative diagram
\[ \xymatrix{ H^i_\mathrm{rig}(\mathcal{Y}_L/\mathcal{R}_L)^{\nabla=0} \ar[r]^-{\cong} &  H^i_{\log\text{-}\mathrm{cris}}(\mathcal{Y}_{k_L}^\times/W^\times_L)\otimes {\mathbb Q} \\
H^i_\mathrm{rig}(\mathcal{X}_K/\mathcal{R})^{\nabla=0} \ar[u] \ar[r]^-{\cong} & H^i_{\mathrm{rig}}(\mathcal{X}_k/\kappa). \ar[u]  } \]
Both vertical maps are injective by Poincar\'e duality (in the log crystalline case, see \cite{Tsu99a}). In particular, it suffices to verify that $[\mathcal{Z}_K]$ and $[\mathcal{Z}_k]$ map to the same element in $H^{2d}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}_{k_L}^\times/W_L)\otimes {\mathbb Q}$.

For any algebraic space over $\mathcal{O}_K$, we will let $K(-)$ denote the $K$-theory of vector bundles. We therefore have total Chern class maps
\begin{align*} c: K_0(\mathcal{X}_K) &\rightarrow \bigoplus_{n\geq0} H^{2n}_\mathrm{cris}(\mathcal{X}_K/\widehat{\mathcal{R}}^\mathrm{int}) \\
c: K_0(\mathcal{Y}_K) &\rightarrow \bigoplus_{n\geq0} H^{2n}_\mathrm{cris}(\mathcal{Y}_K/\widehat{\mathcal{R}}_L^\mathrm{int}) \\
c: K_0(\mathcal{X}_k) &\rightarrow \bigoplus_{n\geq0} H^{2n}_\mathrm{cris}(\mathcal{X}_k/W)
\end{align*}
taking values in crystalline cohomology, as constructed in \cite{BI70}. Using exactly the same method as in \cite{BI70} we can also construct total Chern class maps
\begin{align*} c: K_0(\mathcal{Y}) &\rightarrow \bigoplus_{n\geq0} H^{2n}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}^\times/W_L\pow{t_L}^\times) \\
c: K_0(\mathcal{Y}_{k_L}) &\rightarrow \bigoplus_{n\geq0} H^{2n}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}^\times_{k_L}/W_L^\times)
\end{align*}
in log-crystalline cohomology. Indeed, if we have a $p$-adic log PD base $S$, (e.g. $W^\times$ or $W\pow{t}^\times$) and an fs log scheme $X$ equipped with a log smooth and proper morphism $X\rightarrow S$, then for any rank $r$ vector bundle $\mathcal{E}$ on $X$, we can endow its projectivisation $g:\P(\mathcal{E})\rightarrow X$ with the pullback log-structure from $X$. The first Chern class
\[ \xi:= c_1(\mathcal{O}_{\P(\mathcal{E})}(1))
\in H^2_{\log\text{-}\mathrm{cris}}(\P(\mathcal{E})/W)\]
can be defined directly, and in order to follow the contruction in \cite{BI70} we need to know that the cohomology $H^*_{\log\text{-}\mathrm{cris}}(\P(\mathcal{E})/W)$ decomposes as a direct sum
\[ H^*_{\log\text{-}\mathrm{cris}}(\P(\mathcal{E})/S) \cong \bigoplus_{i=0}^{r-1}H^*_{\log\text{-}\mathrm{cris}}(X/S)\cdot \xi^i.	 \]
We can turn this into a local statement, namely that the higher direct image $\mathbf{R}g_*\mathcal{O}^{\log\text{-}\mathrm{cris}}_{\P(\mathcal{E})/S}$ decomposes as
\[ \mathbf{R}g_*\mathcal{O}^{\log\text{-}\mathrm{cris}}_{\P(\mathcal{E})/S} \cong \bigoplus_{i=0}^{r-1}\mathcal{O}^{\log\text{-}\mathrm{cris}}_{X/S} \cdot \xi^i,  \]
where we have abused notation and also written $\xi\in H^0_{\log\text{-}\mathrm{cris}}(X/S,\mathbf{R}^2g_* \mathcal{O}^{\log\text{-}\mathrm{cris}}_{\P(\mathcal{E})/S} )$ for the image of the first Chern class of $\mathcal{O}_{\P(\mathcal{E})}(1)$ under the edge morphism of the Leray spectral sequence. Since this statement is now local, we can assume that in fact $\mathcal{E}=\mathcal{O}_{X}^{\oplus r}$, thus $\P(\mathcal{E})=\P^{r-1}_X$. In this case the claim follows from the smooth and proper base change theorem in log-crystalline cohomology \cite[Theorem 6.10]{Kat89}, together with the standard computation of the crystalline cohomology of projective space.

Now, since $\mathcal{Z}$ is flat, the fibres $\mathcal{Z}_K$ and $\mathcal{Z}_k$ are also of codimension $d$ in $\mathcal{X}_K$ and $\mathcal{X}_k$ respectively. Since $\mathcal{X}_K$ and $\mathcal{X}_k$ are both regular and  projective schemes, the $K$-theory of vector bundles coincides with that of coherent sheaves, so we have well-defined classes $[\mathcal{O}_{\mathcal{Z}_K}]\in K_0(\mathcal{X}_K)$ and $[\mathcal{O}_{\mathcal{Z}_k}]\in K_0(\mathcal{X}_k)$. Essentially by definition \cite[Theorem 4.3.1]{GM87}, we have
\[ [\mathcal{Z}_K]=\frac{(-1)^{d-1}}{(d-1)!}c_d([\mathcal{O}_{\mathcal{Z}_K}])  \mbox{ \quad and \quad } [\mathcal{Z}_k] =\frac{(-1)^{d-1}}{(d-1)!}c_d([\mathcal{O}_{\mathcal{Z}_k}]). \]
Since $\mathcal{O}_{\mathcal{Z}_K}$ and $\mathcal{O}_{Z_k}$ are strictly perfect complexes on $\mathcal{X}_K$ and $\mathcal{X}_k$ respectively, we get well-defined classes $[\mathbf{L}\pi^*\mathcal{O}_{\mathcal{Z}_K}]\in K_0(\mathcal{Y}_{L})$ and $[\mathbf{L}\pi^*\mathcal{O}_{\mathcal{Z}_k}]\in K_0(\mathcal{Y}_{k_L})$. By functoriality of Chern classes, it suffices to show that the isomorphism
\[ H^{2d}_\mathrm{rig}(\mathcal{Y}_L/\mathcal{R}_L)^{\nabla=0} \overset{\cong}{\longrightarrow} H^{2d}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}_{k_L}^\times/W_L^\times)\otimes {\mathbb Q} \]
sends $c_d([\mathbf{L}\pi^*\mathcal{O}_{\mathcal{Z}_K}])$ to $c_d([\mathbf{L}\pi^*\mathcal{O}_{\mathcal{Z}_k}])$. Since $\mathcal{X}$ is regular we know that $\mathcal{O}_\mathcal{Z}$ is a perfect complex of $\mathcal{O}_\mathcal{X}$-modules, and since $\mathcal{Y}$ is a projective scheme, $\mathbf{L}\pi^*\mathcal{O}_\mathcal{Z}$ is actually a \emph{strictly} perfect complex on $\mathcal{Y}$. It therefore has a well defined class in $K_0(\mathcal{Y})$, which restricts to the class of $\mathbf{L}\pi^*\mathcal{O}_{\mathcal{Z}_K}$ (resp. $\mathbf{L}\pi^*\mathcal{O}_{\mathcal{Z}_k}$) on the generic (resp. special) fibre. Given the construction of the isomorphism
\[ H^{2d}_\mathrm{rig}(\mathcal{Y}_L/\mathcal{R}_L)^{\nabla=0} \overset{\cong}{\longrightarrow} H^{2d}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}_{k_L}^\times/W_L^\times)\otimes {\mathbb Q} \]
it now simply suffices to note that the diagram
\[ \xymatrix{ K_0(\mathcal{Y}_K) \ar[d]^c & K_0(\mathcal{Y})\ar[d]^c \ar[l]  \ar[r] & K_0(\mathcal{Y}_k)\ar[d]^c \\ \bigoplus_{n\geq0} H^{2n}_{\mathrm{cris}}(\mathcal{Y}_L/\widehat{\mathcal{R}}_L^\mathrm{int}) &
\bigoplus_{n\geq0} H^{2n}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}^\times/ W_L\pow{t_L}^\times) \ar[r]\ar[l] & \bigoplus_{n\geq0} H^{2n}_{\log\text{-}\mathrm{cris}}(\mathcal{Y}_{k_L}^\times/ W_L^\times)  } \]
commutes.
\end{proof}

\begin{remark} We need to argue on the alteration $\mathcal{Y}$ rather than on the model $\mathcal{X}$ itself because of the potential difference between the $K$-group of vector bundles (equivalently: strictly perfect complexes) and that of perfect complexes. They coincide on regular schemes (like $\mathcal{Y}$), but it is not known whether the same is true for algebraic spaces (like $\mathcal{X}$). If it were, the above proof could be significantly simplified.
\end{remark}

\section{Descending good reduction under separable extensions} \label{sec: recap}

In this section, we will explain how and why the main results of \cite{LM18} and \cite{CLL17} carry over in equicharacteristic $p$ under a strengthening of Assumption ($\star$), which we will call Assumption ($\star^{\rm s}$). Under this stronger assumption, we also establish a version of these results in terms of the Robba ring valued $p$-adic cohomology discussed in \S\ref{sec: robba}.

\begin{AssumptionStarSep}
A K3 surface $X/K$ satisfies ($\star^{\rm s}$) if it satisfies Assumption ($\star$) and we can moreover take the finite extension $L$ to be separable over $K$.
\end{AssumptionStarSep}

Then we can use the arguments of \cite{LM18} word-for-word to prove the following.

\begin{theorem} \label{theo: LM18 main sep} Let $X$ be a K3 surface over $K$, satisfying Assumption $(\star^{\rm s})$, and let $\ell$ be a prime $\neq p$. Consider the following three statements:
\begin{enumerate}
\item \label{num: urQl} the $G_K$-action on $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified;
\item \label{num: grux} $X$ has good reduction over a finite, unramified extension $L/K$;
\item \label{num: RDPmod} $X$ admits a model $\mathcal{X}$ that is projective over $\mathcal{O}_K$, and whose special fibre is a K3 surface with at worts RDP singularities.
\end{enumerate}
Then, we have $(\ref{num: urQl})\Leftrightarrow(\ref{num: grux})\Rightarrow (\ref{num: RDPmod})$.
\end{theorem}

\begin{proof}
 The point is that the characteristic $0$ hypothesis is only used in \cite{LM18} in order to ensure that $X$ admits a Kulikov model over a Galois extension $L/K$,
which enables arguments via Galois descent to be used.
Thus, if we admit this from the start, the proofs in \cite{LM18} go through unchanged, with one important exception.

This exception occurs in the proof of \cite[Proposition 4.2]{LM18}, showing the existence of flops on relative surfaces over $\mathcal{O}_K$ with numerically trivial relative canonical divisor. One step of the construction uses the fact that a deformation
\[ \spf{\frac{\mathcal{O}_K\pow{x,y,z}}{z^2-H_1(x,y)z-H_0(x,y)}}\]
of a rational double point singularity over $k$ admits a non-trivial involution
\[ z\mapsto H_1(x,y)-z.\]
This will clearly remain true provided $\mathrm{char}(K)\neq 2$, and to show that it also remains true when $\mathrm{char}(K)=2$, we need to explain why we will always have $H_1(x,y)\neq 0$.

This was pointed out to us by Yuya Matsumoto. Indeed, suppose that we have $H_1(x,y) = 0$. Then, since the equation
\[ \frac{\partial H_0}{\partial x} = \frac{\partial H_0}{\partial y} =
0\]
has the solution $(x,y) = (0,0)$ on the special fibre, after possibly replacing $K$ by a finite extension, it also has a solution $(x,y) = (\xi, \eta)$ on the generic fibre. This would mean that the generic fibre has a singular point $(x,y,z) = (\xi, \eta, \sqrt{H_0(\xi,\eta)})$. Therefore smoothness of the generic fibre means that we must have $H_1(x,y)\neq 0$ after all.
\end{proof}

In order to have a version of this result also in terms of the Robba ring valued $p$-adic cohomology discussed in \S\ref{sec: robba}, the first step is to have an analogue of \cite[Theorem 1.1]{Mat15} in this cohomology theory. This is provided by \cite[Theorem 6.4]{CL16}. The next key point will be to show that smooth models descend under finite, \emph{totally ramified} extensions, under suitable cohomological assumptions. More precisely, we need to establish the analogue of \cite[Proposition 5.8]{LM18} in Robba-ring value rigid cohomology.

Thus, let $X$ be a K3 surface over $K$. Suppose that there exists a totally ramified Galois extension $L/K$ with Galois group $G$ such that $X_L$
admits a smooth model $\mathcal{X}\rightarrow \spec{\mathcal{O}_L}$. Let $\mathcal{R}$ be the Robba ring over $\kappa$ and let $\mathcal{R}_L$ denote the unique unramified extension of $\mathcal{R}$ corresponding to $L/K$ as in \cite{Mat95}, see also the proof of Proposition \ref{prop: cyccomp}. Then, there is a natural $G$-action on the $(\varphi,\nabla)$-module $H^2_\mathrm{rig}(X_L/\mathcal{R}_L)$ and hence, on its subspace $H^2_\mathrm{rig}(X_L/\mathcal{R}_L)^{\nabla=0}$
of horizontal sections.

\begin{proposition} \label{prop: extension}
Under the previous assumptions, assume that the $G$-action on $H^2_\mathrm{rig}(X_L/\mathcal{R}_L)^{\nabla=0}$ is trivial.
 Then, the $G$-action on $X_L$ extends to $\mathcal{X}$ and the induced $G$-action on the special fibre $\mathcal{X}_{k}$ is trivial.
\end{proposition}

\begin{proof}
The key ingredient here is Proposition \ref{prop: cyccomp}. Indeed, given this we can copy the proof of \cite[Proposition 5.5]{LM18} word for word to show that the $G$-action on $X_L$ extends to $\mathcal{X}$. To show that the induced action on the special fibre is trivial, we use \cite[Corollary 2.5]{Ogu79} and \cite[Theorem 1.4]{Keu16}, which give injectivity of the natural map
\[\mathrm{Aut}(Y) \rightarrow \mathrm{GL}(H^2_\mathrm{rig}(Y/\kappa))  \]
for K3 surfaces over $k$.
\end{proof}

After this preparation, we arrive at the following $p$-adic version of Theorem \ref{theo: LM18 main sep}.

\begin{theorem} \label{theo: LM18 main padic} Assume that $\mathrm{char}(K)=p>0$, let $X$ be a K3 surface over $K$, satisfying  Assumption $(\star^{\rm s})$, and let $\ell$ be a prime $\neq p$. Then $(\ref{num: urQl})$ and $(\ref{num: grux})$ in Theorem \ref{theo: LM18 main sep} are both equivalent to the following:
\begin{enumerate}\item[$(1)'$] the $(\varphi,\nabla)$-module $H^2_\mathrm{rig}(X/\mathcal{R})$ admits a basis of horizontal sections.
\end{enumerate}
\end{theorem}

\begin{proof}
The implication $(\ref{num: grux})\Rightarrow (1)'$ follows from Proposition \ref{prop: spbc}, together with the fact that a $(\varphi,\nabla)$-module $M$ over $\mathcal{R}$ admits a basis of horizontal sections if and only if it does so after replacing $\kappa$ by a finite extension.

For the implication $(1)'\Rightarrow (\ref{num: grux})$, the main extra input we need is the correct analogue of \cite[Proposition 5.8]{LM18}, which is provided by
Proposition \ref{prop: extension}. The completion of the proof of Theorem \ref{theo: LM18 main padic} is now geometric and identical to the $\ell$-adic case.
\end{proof}

We can also extend the results of \cite{CLL17} to equicharacteristic with essentially no difficulties.

\begin{theorem}
Assume that $\mathrm{char}(K)=p>0$ and let $X$ be a K3 surface over $K$, satisfying  Assumption $(\star^{\rm s})$. Then, then following are equivalent:
\begin{enumerate}
\item $X$ has good reduction over $K$;
\item for some prime (resp. all primes) $\ell \neq p$, the $G_K$-representation $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified, and there is an isomorphism
\[ H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell) \overset{\cong}{\longrightarrow} H^2_\mathrm{\acute{e}t}(Y_{\bar k},{\mathbb Q}_\ell)\]
of $G_k$-representations.
\item the $(\varphi,\nabla)$-module $H^2_\mathrm{rig}(X/\mathcal{R})$ admits a basis of horizontal sections, and there is an isomorphism
\[ H^2_\mathrm{rig}(X/\mathcal{R})^{\nabla=0} \overset{\cong}{\longrightarrow} H^2_\mathrm{rig}(Y/\kappa) \]
of $F$-isocrystals over $\kappa$.
\end{enumerate}
\end{theorem}

\begin{proof}
The existence of flops that is required in \cite[Proposition 7.5]{CLL17}, follows from the same argument we used during the proof of Theorem \ref{theo: LM18 main sep}. Apart from this,
everything in \cite[\S7]{CLL17} works exactly the same in positive equicharacteristic.
The same is true for \cite[\S8]{CLL17}, replacing $\mathbb{D}_\mathrm{cris}(H^2_\mathrm{\acute{e}t}(X_{\overline{K}},{\mathbb Q}_p))$ everywhere by $H^2_\mathrm{rig}(X/\mathcal{R})^{\nabla=0}$, and using Proposition \ref{prop: cyccomp} in order to prove the correct analogue of \cite[Lemma 8.1]{CLL17}. Everything in \cite[\S9]{CLL17} goes through unchanged.
\end{proof}

\section{Descending good reduction under purely inseparable extensions} \label{sec: descent}

In this section, we shall show that smooth models of K3 surfaces descend under purely inseparable extensions of $K$,
although we will start off in rather more generality than that.
Let $S$ be a scheme and let $p:S'\rightarrow S$ be a purely inseparable (that is, radicial, or universally injective) fpqc covering.

\begin{lemma} The map $p:S'\rightarrow S$ is separated and the diagonal embedding $\Delta: S'\rightarrow S'\times_S S'$ is a nilpotent closed immersion.
\end{lemma}

\begin{proof}
Since $p$ is surjective and universally injective, it is universally bijective.
In particular, we see that either projection $p_i:S'\times_S S'\rightarrow S'$ is bijective.
Since $\Delta$ is a section of $p_i$, we deduce that $\Delta$ and $p_i$ are in fact both homeomorphisms.
As $\Delta$ is a locally closed embedding with closed image, it must therefore be a closed immersion, and $p$ is therefore separated.
Since $\Delta$ is a closed immersion inducing a homeomorphism on the underlying topological spaces, it is defined by a nilpotent sheaf of ideals.
\end{proof}

We denote by $\mathcal{I}$ the ideal of the diagonal $S'\hookrightarrow S'\times_S S'$, and write
\[ S'^{(n)}:= \underline{\mathrm{Spec}}_{S'\times_S S'}(\mathcal{O}_{S'\times_S S'}/\mathcal{I}^{n+1})\]
for the $n$th infinitesimal neighbourhood of $S'$ in $S'\times_S S'$.
Let $f:Y\rightarrow S'$ be an algebraic space with relative tangent sheaf $\mathcal{T}_{Y/S'}:=\mathcal{H}\mathrm{om}_{\mathcal{O}_Y}(\Omega^1_{Y/S'},\mathcal{O}_Y)\in \mathrm{QCoh}(Y)$.
For any morphism of schemes $U\rightarrow S$ we write $(-)_U$ for the fibre product $(-)\times_S U$.

\begin{proposition} \label{prop: extending descent}Assume that:
\begin{enumerate}
\item the ideal $\mathcal{I}$ is of finite type;
\item each $S'^{(n)}$ is flat over $S'$ (under either projection);
\item $Y$ is smooth and separated over $S'$;
\item $f_*\mathcal{T}_{Y/S'}=0$, and $\mathbf{R}^1f_*\mathcal{T}_{Y/S'}$ is torsion free.
\end{enumerate}
Then, for any dense open subscheme $U\subset S$, the function
\[ \left\{\text{descent data on }Y\text{ along }S'\rightarrow S \right\} \rightarrow \left\{\text{descent data on }Y_U\text{ along }S_U'\rightarrow U \right\}\]
obtained by restriction is a bijection.
\end{proposition}

\begin{remark} Note that we do not assume $Y$ proper over $S'$.
However, the condition that $f_*\mathcal{T}_{Y/S'}=0$ will generally force this in practice.
\end{remark}

\begin{proof} Write $p_i^*$ for pullback along either projection $p_i:S'\times_S S'\rightarrow S'$.
The given function is injective since $Y$ is separated over $S'$, thus an isomorphism $p_1^*Y_U \rightarrow p_2^*Y_U$ can extend
to at most one isomorphism $p_1^*Y\overset{\sim}{\rightarrow} p_2^*Y$.
To see that it is surjective, it suffices to show that we can always extend such an isomorphism
(the fact that such an extension will define a descent datum again
follows from separatedness of $Y\rightarrow S'$ and the corresponding fact for $Y_U$). Since $S'\rightarrow S'\times_S S'$ is a homeomorphism, $S'\times_S S'$ is covered by open sets of the form $V\times_S V$ with $V\subset S'$ open and affine. Therefore appealing once more to the fact that $U\subset S$ is dense and $Y\rightarrow S'$ is separated, we can see that the problem of extending $\alpha_U$ is local on both $S$ and $S'$, which we may therefore assume to be affine, say $S=\spec{R}$ and $S'=\spec{R'}$. (We may have lost surjectivity of $S'\rightarrow S$ but this doesn't matter). We may also assume that $U=\spec{R[g^{-1}]}$ is a basic open affine in $S$.

We will prove that any such $\alpha_U:p_1^*Y_U \rightarrow p_2^*Y_U$ extends by using deformation theory. Let $R'':=R'\otimes_R R'$ and let $I=\ker \left( R''\rightarrow R'\right)$.
Thus, $I$ is a nilpotent ideal. By assumption, it is finitely generated and hence, $I^N=0$ for $N$ large enough.
If we write $R'^{(n)}:= R''/I^{n+1}$, then by assumption each $R'^{(n)}$ is flat over $R$, hence so is the kernel $I^n/I^{n+1}$
of the surjection $R'^{(n)}\rightarrow R'^{(n-1)}$.
Let $p_i^{(n)*}Y$ denote the base change to $R'^{(n)}$ along either `projection' map $p_i^{(n)}:R'\rightarrow R^{'(n)}$.
We will show by induction on $n\geq 1$ that $\alpha_U^{(n)}: p_1^{(n)*}Y_U \overset{\sim}{\rightarrow} p_2^{(n)*}Y_U$ extends to an isomorphism
$p_1^{(n)*}Y \overset{\sim}{\rightarrow} p_2^{(n)*}Y$.
Since $I^N=0$ for large enough $N$, this suffices to prove the result.

The case $n=0$ is trivial and we may therefore suppose that $n\geq 1$ and that we are given $p_1^{(n-1)*}Y \overset{\sim}{\rightarrow} p_2^{(n-1)*}Y$ extending $\alpha_U^{(n-1)}$.
We can thus view both $p_1^{(n)*}Y$ and $p_2^{(n)*}Y$ as deformations of $p_1^{(n-1)*}Y$ along $R'^{(n)}\twoheadrightarrow R'^{(n-1)}$.
By flatness of $I^n/I^{n+1}$ over $R'$, the isomorphism classes of such deformations are controlled by
\[ H^1( Y, \mathcal{T}_{Y/R'})\otimes_{R'} I^n/I^{n+1} ,\]
and we are given that these two classes map to the same element in
\[ H^1( Y, \mathcal{T}_{Y/R'})\otimes_{R'} I^n/I^{n+1} \otimes_R R[g^{-1}]. \]
But by the torsion-free hypothesis on $\mathbf{R}^1f_*\mathcal{T}_{Y/S'}$, we know that
\[ H^1( Y, \mathcal{T}_{Y/R'})\otimes_{R'} I^n/I^{n+1} \hookrightarrow H^1( Y, \mathcal{T}_{Y/R'})\otimes_{R'} I^n/I^{n+1} \otimes_R R[g^{-1}] \]
is injective, from which we deduce that $p_1^{(n)*}Y$ and $p_2^{(n)*}Y$ are isomorphic as deformations of $p_1^{(n-1)*}Y$.
Moreover, the automorphism group of a given deformation of $p_1^{*(n-1)}Y_U$ is
\[ H^0( Y, \mathcal{T}_{Y/R'})\otimes_{R'} I^n/I^{n+1} \otimes_R R[g^{-1}] =0, \]
since $f_*\mathcal{T}_{Y/S'}=0$, from which we deduce that \emph{any} isomorphism $p_1^{(n)*}Y \overset{\sim}{\rightarrow} p_2^{*(n)}Y$ as deformations of $p_1^{(n-1)*}Y$
restricts to the given isomorphism $\alpha_U^{(n)}: p_1^{(n)*}Y_U \overset{\sim}{\rightarrow} p_2^{(n)*}Y_U$. \end{proof}

Concerning effectivity of descent, we have the following, presumably well-known, result.

\begin{proposition} \label{prop: effective}
 For any algebraic space $Y\rightarrow S'$, all descent data on $Y$ along $p:S'\rightarrow S$ are effective.
\end{proposition}

\begin{proof}
This is entirely straightforward for schemes and not much harder for algebraic spaces.
Let $Y\cong [R \rightrightarrows U]$ be a presentation of $Y$ via an \'etale equivalence relation of affine schemes over $S'$.
Let $\alpha:p_1^*Y \overset{\sim}{\rightarrow} p_2^*Y$ be an isomorphism defining a descent datum on $Y$.
Since the diagonal $S'\rightarrow S'\times_S S'$ is a nilpotent closed immersion, it is a universal homeomorphism, and we can deduce from topological invariance
of the small \'etale site (in particular, from the equivalence $\mathrm{\acute{E}t}(Y)\cong \mathrm{\acute{E}t}(p_i^*Y)$ for each $i$)
that we can uniquely extend $\alpha$ to an isomorphism
\[ p_1^*R\rightrightarrows p_1^*U \cong p_2^* R\rightrightarrows p_2^*U \]
of \'etale equivalence relations. This isomorphism is compatible on $S'\times_S S'\times_S S'$ and hence (since $U$ and $R$ are affine), descends to an \'etale equivalence
relation $R_S \rightrightarrows U_S $ of affine schemes over $S$.
Thus, taking $X\cong [R_S\rightrightarrows U_S]$ to be the corresponding quotient we get the required descent of $Y$ to $S$.
\end{proof}

We will apply these results as follows.
Let $\mathcal{O}_K$ be our fixed excellent and Henselian DVR with fraction field $K$.
Let $X$ be a smooth, proper, and geometrically connected scheme over $K$ and $L/K$ a finite and purely inseparable extension.
The integral closure $\mathcal{O}_L$ of $\mathcal{O}_K$ inside $L$ is then also an excellent and Henselian DVR.
Suppose that we have a smooth model $\mathcal{Y}$ for $X_L$ over $\mathcal{O}_L$.

\begin{theorem}\label{thm: gr desc}
Assume that $H^0(X,\mathcal{T}_{X/K})=0$ and that $H^1(\mathcal{Y},\mathcal{T}_{\mathcal{Y}/\mathcal{O}_L})$ is torsion free.
Then, $\mathcal{Y}$ descends uniquely to a smooth model for $X$ over $\mathcal{O}_K$.
\end{theorem}

\begin{proof}
Since $X_L$ comes with a canonical descent datum, we apply Proposition \ref{prop: extending descent}
with $S=\spec{\mathcal{O}_K}$, $S'=\spec{\mathcal{O}_L}$ and $U=\spec{K}$ to extend this descent datum to $\mathcal{Y}$,
and then apply Propostion \ref{prop: effective} to deduce that this descent datum is effective.
\end{proof}

Since K3 surfaces do not admit non-zero global vector fields \cite{RS76}, the previous result applies to them:

\begin{corollary} \label{cor: gr desc}
 Suppose that $X$ is a K3 surface over $K$.
 Then, any smooth model for $X_L$ descends uniquely to a smooth model for $X$.
 In particular, $X$ has good reduction over $K$ if and only if $X_L$ has good reduction over $L$.
\end{corollary}

\section{Main results} \label{sec: together}

Using the results of \S\ref{sec: descent}, it is now straightforward to deduce our main results,
via the following proposition.

\begin{proposition}
 Suppose that $\mathrm{char}(K)=p>0$ and let $X$ be a K3 surface over $K$.
 Assume that either of the following conditions holds:
\begin{enumerate}
\item for some prime $\ell\neq p$ the $G_K$-representation $H^2_\mathrm{\acute{e}t}(X_{K^{\rm s}},{\mathbb Q}_\ell)$ is unramified;
\item the $(\varphi,\nabla)$-module $H^2_\mathrm{rig}(X/\mathcal{R})$ admits a basis of horizontal sections.
\end{enumerate}
Then, $X$ satisfies Assumption $(\star)$ if and only if it satisfies Assumption $(\star^{\rm s})$.
\end{proposition}

\begin{proof}
We clearly have that Assumption $(\star^{\rm s})$ implies Assumption $(\star)$.

For the converse direction, let $L/K$ be a finite extension over which $X$ admits a Kulikov model.
Note that either condition of the proposition also holds for the base change $X_L$ and hence, by \cite[Theorem 6.4]{CL16} we know that $X_L$ has good reduction.
If we let $K\subset K_1\subset  L$ denote the maximal separable subextension, it therefore follows from Corollary \ref{cor: gr desc}
that $X_{K_1}$ has good reduction. In particular, $X$ admits a Kulikov model over $K_1$ and thus, satisfies Assumption $(\star)$.
\end{proof}

Applying the results of \S\ref{sec: recap}, Theorems \ref{theo: LM18 main equi} and \ref{theo: CLL17 main equi} now follow immediately.

\section{A counter-example}\label{sec: counter}

In this section, we will explain how to modify the example given in \cite[\S7]{LM18} to produce K3 surfaces in equicharacteristic $p>0$ that only admit good reduction after a \emph{non-trivial} finite and unramified extension. Let $p\geq 5$, choose $c\in {\mathbb F}_p^*\setminus ({\mathbb F}_p^*)^2$, $a\in {\mathbb F}_p^*\setminus \left\{\frac{16}{27}\right\}$ and let $F\in {\mathbb F}_p[t][x,y,z,w]$ be the polynomial
\[ F:= w(x^3+y^3+z^3+(t+a(1-t))w^3)+(tz^2+xy+tyz)^2-(c-t)t^2y^2z^2. \]
Define $\mathcal{X}:= V(F) \subset \P^3_{{\mathbb F}_p\pow{t}}$ and let $X=\mathcal{X}_{{\mathbb F}_p\lser{t}}$ be its generic fibre.

\begin{theorem} \label{theo: counter example}
$X$ is a smooth K3 surface over ${\mathbb F}_p\lser{t}$ that admits good reduction over ${\mathbb F}_{p^2}\lser{t}$ but not over ${\mathbb F}_p\lser{t}$.
\end{theorem}

\begin{proof}
Exactly as in \cite[Theorem 7.2]{LM18}, we can verify smoothness of $X$ by considering the subscheme cut out by $F$ in $\P^3_{{\mathbb F}_p[t]}$ and reducing modulo $1-t$. The special fibre of $\mathcal{X}$ is isomorphic to the special fibre of the example $\mathcal{X}(p)$ constructed in \cite{LM18}, it is therefore a singular K3 surface with 6 RDP singularities, two of which are defined over ${\mathbb F}_p$ and the other four over ${\mathbb F}_p[\zeta_3]$. Since $\sqrt{c-t}\in {\mathbb F}_{p^2}\pow{t}$, after extending scalars to ${\mathbb F}_{p^2}\pow{t}$ we can resolve the singularities of $\mathcal{X}$ by blowing up either ideal
\[ \mathcal{I}_{\pm}=(w,tz^2+xy+tyz\pm tyz\sqrt{c-t}).\]
However, these two different ideals give rise to distinct elements of the strictly local Picard group $\mathrm{Pic}(\mathcal{O}^{\rm sh}_{\mathcal{X},\bar{x}})$, which are interchanged by the action of  $\mathrm{Gal}({\mathbb F}_p\lser{t}^{\rm{s}}/{\mathbb F}_p\lser{t})$. Thus, exactly as in \cite[Theorem 7.2]{LM18}, $X$ cannot have a smooth model over ${\mathbb F}_p\pow{t}$.
\end{proof}

\bibliographystyle{../../Templates/bibsty}