\section{Introduction}
\label{sctn:intro}

The goal of this article is to study some phenomena in the
$C_2$-equivariant stable homotopy groups.
Let $\ensuremath{\mathbb{R}}^{n,k}$ be the $n$-dimensional real representation of $C_2$
in which the nonidentity element of
$C_2$ acts as $-1$ on the last $k$ coordinates (and trivially on the first $n-k$),
and let $S^{n,k}$ be its one-point compactification.
Then $\pi^{C_2}_{n,k}$ is the set of
$2$-completed $C_2$-equivariant
stable homotopy classes of maps $S^{n,k} \longrightarrow S^{0,0}$.
In this article, we are primarily concerned with the groups
$\pi^{C_2}_{k,k}$.

The classical Hopf map $\eta_\mathrm{cl}:S^3\longrightarrow S^2$ can be modeled as the defining quotient map $\ensuremath{\mathbb{C}}^2-\{0\} \longrightarrow \ensuremath{\mathbb{C}}\ensuremath{\mathbb{P}}^1$ for complex projective space. When we remember the action of $C_2$ via complex conjugation, this represents a $C_2$-equivariant stable map $\eta$ in $\pi^{C_2}_{1,1}$.
Classically, $\eta_\mathrm{cl}$ is nilpotent in the stable homotopy ring, as is every element in positive stems \cite{N}.
However, the equivariant Hopf map $\eta$ is not nilpotent
because $\eta$ induces the non-nilpotent element $-2$ on
geometric fixed points.
(The distinction between $2$ and $-2$ depends on choices of
orientations and is inconsequential to the argument.)
We will concern ourselves with phenomena associated to the non-zero
elements $\eta^k$ in $\pi^{C_2}_{k,k}$.

Because the fixed points of the representation sphere $S^{k,k}$ consist
of two points, the geometric fixed point homomorphism takes the form
$\phi:\pi_{k,k}^{C_2} \longrightarrow \pi_0\cong \ensuremath{\mathbb{Z}}$.
Bredon \cite{Br} and Landweber \cite{L} proved that the image of $\phi$
is not in general generated by $\phi(\eta^k)$. For instance, $\phi(\eta^5) = (-2)^5=-32$,
but $\phi(\pi^{C_2}_{5,5})=16\ensuremath{\mathbb{Z}}$.
In fact, the higher powers of $\eta$ are increasingly divisible by 2
in the $C_2$-equivariant stable homotopy groups (\autoref{etadiv2}).

Let $\rho: S^{0,0} \longrightarrow S^{1,1}$ be the inclusion of fixed points.
Our main result describes to what extent the powers of $\eta$ are divisible by $\rho$, from which we will deduce several other results.

\begin{thm}\label{mainthm}
Let $k = 4j + \varepsilon \geq 1$, where $0 \leq \epsilon \leq 3$.
If $\varepsilon = 0$, then
the $C_2$-equivariant stable homotopy class $\eta^k$ is divisible by $\rho^{k-1}$ and no higher power of $\rho$.
Otherwise, the $C_2$-equivariant stable homotopy class $\eta^{k}$ is divisible by $\rho^{4j}$ and no higher power of $\rho$.
\end{thm}

Our primary tool for studying $C_2$-equivariant
stable homotopy groups is the equivariant Adams spectral sequence
\cite{G} \cite{HK}.
The {\bf Bredon-Landweber region} refers to
the subgroups of $\pi^{C_2}_{k,k}$ that are detected in
Adams filtration greater than $\frac{1}{2}k - 1$.
This region is displayed in the top part of \autoref{AdamsEinfChart}.
We will show that the Bredon-Landweber region
is additively generated by
the elements $\rho^i \eta^k$, together
with elements $\alpha$ such that
$\rho^i \alpha = \eta^k$ for some $i$.

We recover the following theorem of Landweber \cite[Theorem~2.2]{L}, which was originally conjectured by Bredon \cite{Br}.

\begin{cor}
\label{cor:fixed-pt-image}
Let $k = 8j + \varepsilon \geq 1$, with $0 \leq \varepsilon \leq 7$.
The image of the geometric fixed points homomorphism $\pi^{C_2}_{k,k}\xrtarr{\phi} \pi_0$ is generated by
\[\begin{cases}
2^{4j + 1} & \text{if} \quad \varepsilon = 0. \\
2^{4j + \varepsilon} & \text{if} \quad 1 \leq \varepsilon \leq 4. \\
2^{4j + 4} & \text{if} \quad 5 \leq \varepsilon \leq 7. \\
\end{cases}\]
\end{cor}

Similarly, we have:

\begin{cor}\label{etadiv2}
Let $k = 8j + \varepsilon \geq 5$, with $0 \leq \varepsilon \leq 7$.
The $C_2$-equivariant stable homotopy class $\eta^k$ is divisible by
\[ \begin{cases}
2^{4j-1} & \text{if} \quad \varepsilon = 0, \\
2^{4j} & \text{if} \quad 1 \leq \varepsilon \leq 4, \\
2^{4j+\varepsilon-4} & \text{if} \quad 5 \leq \varepsilon \leq 7, \\
\end{cases} \]
and no higher power of $2$.
\end{cor}

The proofs of \autoref{mainthm}, \autoref{cor:fixed-pt-image},
and \autoref{etadiv2} appear in \autoref{sctn:Adams}.
First, we must carry out some $C_2$-equivariant Adams spectral
sequence calculations.

Our calculations can also be used to compute the classical
Mahowald invariants of $2^k$ for all $k \geq 0$
(see \autoref{thm:Mahowald-inv}).
We directly apply the Bruner-Greenlees formulation \cite{BG}
of the Mahowald invariant that uses $C_2$-equivariant homotopy groups.
These invariants were previously established by
Mahowald and Ravenel \cite{MR} using entirely different methods.

The charts were created using Hood Chatham's {\tt spectralsequences} package.

\section{Notation}
\label{sec:background}

We continue with notation from \cite{DI} and \cite{GHIR} as follows.
\begin{enumerate}
\item
$\ensuremath{\mathbb{M}_2}^\ensuremath{\mathbb{R}}=\ensuremath{\mathbb{F}}_2[\tau,\rho]$
is the motivic cohomology of $\ensuremath{\mathbb{R}}$ with $\ensuremath{\mathbb{F}}_2$ coefficients, where $\tau$ and $\rho$ have bidegrees $(0,1)$ and $(1,1)$, respectively.
\item
$\ensuremath{\mathbb{M}_2}^{C_2}$ is the bigraded $C_2$-equivariant Bredon cohomology of a point with coefficients in the constant Mackey functor $\underline{\ensuremath{\mathbb{F}}}_2$.
\item
$\ensuremath{\mathcal{A}}^{C_2}$ is the
$C_2$-equivariant mod 2 Steenrod algebra, using coefficients $\underline{\ensuremath{\mathbb{F}}}_2$, and
$\ensuremath{\mathcal{A}}^{C_2}(1)$ is the $\ensuremath{\mathbb{M}_2}^{C_2}$-subalgebra of $\ensuremath{\mathcal{A}}^{C_2}$
generated by $\mathrm{Sq}^1$ and $\mathrm{Sq}^2$.
\item
$\Ext_\mathrm{cl}$, $\Ext_\ensuremath{\mathbb{C}}$, $\Ext_\ensuremath{\mathbb{R}}$, and $\Ext_{C_2}$ are the
cohomologies of the classical, $\ensuremath{\mathbb{C}}$-motivic, $\ensuremath{\mathbb{R}}$-motivic, and
$C_2$-equivariant mod 2 Steenrod algebras respectively.
These objects are the $E_2$-pages of Adams spectral sequences.
\item
$\pi_{*,*}^{C_2}$ and $\pi_*$ are the stable homotopy rings of the
2-completed $C_2$-equivariant sphere spectrum and the
2-completed classical sphere spectrum respectively.
\end{enumerate}

We will use some specific familiar elements of the Adams $E_2$-page.
These elements lie near the ``Adams edge" at the top of the Adams
chart along a line of slope $1/2$.
Our notation for these elements is standard.  They include
$P^k h_1$, $P^k h_1^2$, $P^k h_1^3$, $P^k c_0$, $P^k h_1 c_0$,
$P^k h_2$, and $P^k h_0 h_2$.
In addition, we will consider the elements
$P^k h_0 h_3$, $P^k h_0^2 h_3$, and $P^k h_0^3 h_3$.
These slightly non-standard (but technically correct)
names conveniently refer to a well-understood,
regular family of elements in the Adams $E_2$-page.
They are the top three elements in a tower of $h_0$-multiplications
in stems congruent to $7$ modulo $8$.
For more details,
we refer the reader to any Adams chart,
such as \cite{Icharts} or \cite[Figure A3.1]{R}.

We follow \cite{Istems} in grading $\Ext$ groups
according to  $(s,f,w)$, where:
\begin{enumerate}
\item
$f$ is the Adams filtration, i.e., the homological degree.
\item
$s+f$ is the internal degree, i.e., corresponds to the
first coordinate in the bidegree of the Steenrod algebra.
\item
$s$ is the stem, i.e., the internal degree minus
the Adams filtration.
\item
$w$ is the weight.
\end{enumerate}

Following this grading convention,
the elements $\tau$ and $\rho$, as elements of $\Ext_\ensuremath{\mathbb{R}}$, have degrees
$(0,0,-1)$ and $(-1,0,-1)$ respectively.
We will also often refer to the \textbf{coweight},
which is defined to be $c = s-w$.
Since both $\eta$ and $\rho$ have coweight $0$,
the Bredon-Landweber region consists entirely of elements of
coweight $0$. The coweight is also called the Milnor-Witt degree in the motivic context (\cite{DI},\cite{GI}).

\section{The $\rho$-Bockstein spectral sequence}
\label{sctn:rho-Bock}

As an $\ensuremath{\mathbb{M}_2}^\ensuremath{\mathbb{R}}$-module,
the equivariant cofficient ring $\ensuremath{\mathbb{M}_2}^{C_2}$ splits as
$\ensuremath{\mathbb{M}_2}^{C_2} \cong \ensuremath{\mathbb{M}_2}^\ensuremath{\mathbb{R}} \oplus NC$, where $NC$ is the ``negative cone''.
The negative cone has $\ensuremath{\mathbb{F}}_2$-basis $\{ \frac{\gamma}{\rho^j \tau^{k+1}}\}$, where $j,k\geq 0$ and $\frac{\gamma}{\tau}$ lives in degree $(0,0,2)$.
See \cite[Section 2.1]{GHIR} for more details.
This splitting of $\ensuremath{\mathbb{M}_2}^{C_2}$ induces a splitting
\[
\Ext_{C_2} \cong \Ext_\ensuremath{\mathbb{R}} \oplus \Ext_{NC},
\]
where $\Ext_{NC}=\Ext_\ensuremath{\mathbb{R}}(NC,\ensuremath{\mathbb{M}_2}^\ensuremath{\mathbb{R}})$.
The splitting of $\ensuremath{\mathbb{M}_2}^{C_2}$
also yields a splitting for the Bockstein spectral sequence, and we
follow \cite[Proposition 3.1]{GHIR} in writing $E_1^+$ for the summand of the Bockstein $E_1$-term which converges to $\Ext_\ensuremath{\mathbb{R}}$ and $E_1^-$ for the summand  which converges to $\Ext_{NC}$.

The $\ensuremath{\mathbb{R}}$-motivic $\rho$-Bockstein spectral sequence (\cite{H,DI}) takes the form
\[ E_1^+ = \Ext_\ensuremath{\mathbb{C}}[\rho] \Rightarrow \Ext_\ensuremath{\mathbb{R}}.\]
The groups $\Ext_\ensuremath{\mathbb{R}}$ are computed for low coweights in \cite{DI}. In coweight $0$, $E_1^+$ is $\ensuremath{\mathbb{F}}_2[h_0,h_1,\rho]/(h_0h_1)$, and the only relevant Bockstein differential is $d_1(\tau) = \rho h_0$, giving

\begin{lem}[\cite{DI}]
\label{lem:ExtR-mw0}
 $\Ext_\ensuremath{\mathbb{R}}$ in coweight 0 is $\ensuremath{\mathbb{F}}_2[h_0,h_1,\rho]/(h_0h_1,\rho h_0)$.
\end{lem}

The calculation of $\Ext_{NC}$ in coweight $0$ is much
more complicated.
We have a short exact sequence \cite[Prop 3.1]{GHIR}
\begin{equation}
\bigoplus_{s\geq 0} \frac{\ensuremath{\mathbb{F}}_2[\tau]}{\tau^\infty}
\left\{\frac\gamma{\rho^s}\right\} \otimes_{\ensuremath{\mathbb{F}}_2[\tau]} \Ext_{\ensuremath{\mathbb{C}}}
\rightarrow E_1^- \rightarrow
\bigoplus_{s\geq 0} \Tor_{\ensuremath{\mathbb{F}}_2[\tau]}
\left( \frac{\ensuremath{\mathbb{F}}_2[\tau]}{\tau^\infty} \left\{\frac\gamma{\rho^s}\right\},
\Ext_{\ensuremath{\mathbb{C}}}\right),
\label{E1NCses}
\end{equation}
which we abbreviate as
\[ \gamma E_1^- \longrightarrow E_1^- \longrightarrow Q E_1^-.\]
As explained in \cite[Remark~3.5]{GHIR},
for each class $x$ in $\Ext_\ensuremath{\mathbb{C}}$ such that $\tau x = 0$,
we get a class $Q x$ in $Q E_1^-$, and this element is infinitely
divisible by $\rho$.
The $\tau$-torsion elements of $\Ext_\ensuremath{\mathbb{C}}$ in coweight 0 are $h_1^k$ for $k\geq 4$, so these give rise to infinitely $\rho$-divisible classes $Qh_1^k$.  This describes $QE_1^-$ in coweight $0$.

We now describe $\gamma E_1^-$ in coweight $0$.
First note that $\frac{\gamma}{\tau^i}$ has coweight
$-i-1$.
Now let $x$ be a class in $\Ext_\ensuremath{\mathbb{C}}$ that is $\tau$-free
and not divisible by $\tau$, and let $c \geq 0$ be the coweight
of $x$.
If $c \geq 2$,
then $\frac{\gamma}{\tau^{c-1}} x$ is an element of $\gamma E_1^-$
in coweight $0$
that is infinitely divisible by $\rho$.
When $c \leq 1$, there is no corresponding
element of $\gamma E_1^-$ in coweight $0$.

This description of $E_1^-$ is incomplete in the sense that it depends
on the $\tau$-free part of $\Ext_\ensuremath{\mathbb{C}}$, which is only known in a range.
The $\tau$-free part of $\Ext_\ensuremath{\mathbb{C}}$ corresponds precisely to
classical $\Ext_{\mathrm{cl}}$ \cite[Proposition 2.10]{Istems}.
In a range, information about
$\Ext_{\mathrm{cl}}$ can be obtained from an $\Ext$ chart, such as
\cite{Icharts} or \cite[Figure A.3.1]{R}.

In order to rule out certain Bockstein differentials later, we need some
structural results for the Bockstein spectral sequence.

\begin{prop}
\label{prop:rho-cofree-diff}
Let $x$ be an element of $E_r^-$ such that
$d_r(x)$ is non-zero.  Then $x$ and $d_r(x)$ are both
infinitely divisible by $\rho$ in $E_r^-$.
\end{prop}

\begin{proof}
Let $E_r^-[k]$ be the part of $E_r^-$ in Bockstein filtration
$k$.  Note that $E_1^-[k]$ is zero if $k > 0$, and that
$\rho: E_1^-[k] \longrightarrow E_1^-[k+1]$ is an isomorphism if $k < 0$.
The $d_r$ differential takes the form
$E_r^-[k] \longrightarrow E_r^-[k+r]$.

By induction,
diagram chases show that
$\rho: E_r^-[k] \longrightarrow E_r^-[k+1]$ is injective if
$0 > k > -r$, and it is an isomorphism if $-r \geq k$.
In particular, if $-r \geq k$, then every element
of $E_r^-[k]$ is infinitely divisible by $\rho$.

Now let $x$ be an element of $E_r^-[k]$ such that $d_r(x)$ is
non-zero.  This implies that $E_r^-[k+r]$ is non-zero, so
$-r \geq k$, and $x$ is infinitely divisible by $\rho$.
Finally, the multiplicative structure implies that $d_r(x)$
must also be infinitely divisible by $\rho$.
\end{proof}

\begin{rem}
\label{rem:rho-cofree-diff}
\autoref{prop:rho-cofree-diff} is dual to
\cite[Lemma 3.4]{DI}, which shows that if $d_r(x)$ is
non-zero in $E_r^+$, then $\rho^k x$ and $\rho^k d_r(x)$
are non-zero in $E_r^+$ for all $k \geq 0$.
In fact, the proof of \autoref{prop:rho-cofree-diff} dualizes
line by line.
\end{rem}

\begin{prop}\label{NegConeVanish}
The $C_2$-equivariant Bockstein $E_1$-page is zero
in degrees $(s,f,w)$ such that
the coweight $s - w$ is negative,
the stem $s$ is positive, and
$f > \frac{1}{2} s + \frac{3}{2}$.
\end{prop}

\begin{pf}
The summand $E_1^+$ vanishes when $s - w < 0$, i.e., in negative
coweights.
Similarly, $QE_1^-$ vanishes in negative coweights.

It only remains to consider $\gamma E_1^-$.
Consider a non-zero element of $\gamma E_1^-$ in degree $(s,f,w)$
with $s > 0$.
This element has the form
$\frac{\gamma}{\rho^j \tau^k} x$, where $x$ is $\tau$-free in $\Ext_\ensuremath{\mathbb{C}}$.
Moreover, the degree of $x$ is $(s - j, f, w - j - k - 1)$.

Using a vanishing result for $\Ext_\ensuremath{\mathbb{C}}$
\cite[Theorem~1.1]{GIVanish},
we know that
\[
f \leq \frac{1}{2} (s-k) + \frac{3}{2}.
\]
Since $k$ is non-negative, it follows that
$f \leq \frac{1}{2} s + \frac{3}{2}$.
\end{pf}

\begin{lem}\label{h1PerNC1}
In coweight $1$, the
localization $E_1[h_1^{-1}]$ of the Bockstein $E_1$-page
vanishes.
\end{lem}

\begin{pf}
We know that $\Ext_\ensuremath{\mathbb{C}}[h_1^{-1}]$ vanishes in
coweight $1$ \cite[Theorem 1.1]{GI}.
Therefore $E_1^+[h_1^{-1}]$ and $QE_1^-[h_1^{-1}]$ both vanish.
Finally, $\gamma E_1^- [h_1^{-1}]$ also vanishes
because there are no $\tau$-free classes in
$\Ext_\ensuremath{\mathbb{C}}[h_1^{-1}]$.
\end{pf}

\section{Some Bockstein differentials}

The goal of this section is to compute some explicit Bockstein
differentials that we will need for our analysis of the
Bredon-Landweber region.

\begin{lem}
\label{lem:gamma-diff}
For $k \geq 0$,
\[
d_1\left( \frac{\gamma}{\rho \tau^{2k+1}} \right) =
\frac{\gamma}{\tau^{2k+2}} h_0.
\]
\[
d_2\left( \frac{\gamma}{\rho^2 \tau^{4k+2}} \right) =
\frac{\gamma}{\tau^{4k+3}} h_1.
\]
\[
d_3\left( \frac{\gamma}{\rho^3 \tau^{4k+4}} \right) = 0.
\]
\end{lem}

\begin{proof}
These formulas follow from the Leibniz rule and the
$\ensuremath{\mathbb{R}}$-motivic Bockstein differentials
$d_1(\tau^{2k+1}) = \rho \tau^{2k} h_0$,
$d_2(\tau^{4k+2}) = \rho^2 \tau^{4k+1} h_1$,
and $d_3(\tau^{4k+4}) = 0$
\cite[Proposition 3.2]{DI}.

More specifically, start with the relation
$\tau^{2k+1} \cdot \frac{\gamma}{\rho \tau^{2k+1}} = 0$.
Apply the $d_1$ differential to obtain
\[
0 =
\rho \tau^{2k} h_0 \cdot \frac{\gamma}{\rho \tau^{2k+1}} +
\tau^{2k+1} \cdot d_1 \left( \frac{\gamma}{\rho \tau^{2k+1}} \right) =
\frac{\gamma}{\tau} h_0 +
\tau^{2k+1} \cdot d_1 \left( \frac{\gamma}{\rho \tau^{2k+1}} \right).
\]
Therefore,
$d_1 \left( \frac{\gamma}{\rho \tau^{2k+1}} \right)$ must equal
$\frac{\gamma}{\tau^{2k+2}} h_0$.

The second and third formulas follow from a similar argument, starting with
the relations $\tau^{4k+2} \cdot \frac{\gamma}{\rho^2 \tau^{4k+2}} = 0$
and $\tau^{4k+4} \cdot \frac{\gamma}{\rho^2 \tau^{4k+4}} = 0$.
\end{proof}

\begin{lem}
\label{lem:v1-perm}
For $k \geq 0$, the elements
$\tau P^k h_1$, $P^k h_2$, and $\tau P^k c_0$, are
permanent cycles in the $\ensuremath{\mathbb{R}}$-motivic $\rho$-Bockstein
spectral sequence.
\end{lem}

\begin{proof}
We can express the classes $\tau P^k h_1$ recursively
as matric Massey products \cite{Quig}:
\[
\tau P^k h_1 = \left\langle [ h_3 \  c_0 ], \begin{bmatrix} h_0^4 \\ \rho^3 h_1^2\end{bmatrix}, \tau P^{k-1} h_1\right\rangle.
\]
The May Convergence Theorem \cite[Theorem~4.1]{May}
\cite[Theorem~2.2.1]{Istems},
applied to the $\rho$-Bockstein spectral sequence,
shows that $\tau P^k h_1$ is a permanent cycle.

Similarly, we have recursive matric Massey products
\[
P^k h_2 = \left\langle [ h_3 \  c_0 ], \begin{bmatrix} h_0^4 \\ \rho^3 h_1^2\end{bmatrix}, P^{k-1} h_2 \right\rangle
\]\[
\tau P^k c_0 = \left\langle [ h_3 \  c_0 ], \begin{bmatrix} h_0^4 \\ \rho^3 h_1^2\end{bmatrix}, \tau P^{k-1} c_0\right\rangle.
\]
\end{proof}

\begin{lem}
\label{lem:d3}
For $k \geq 0$,
\[
d_3(\tau^3 P^k h_0^3 h_3) = \rho^3 \tau P^{k+1} h_1.
\]
\[
d_3(\tau^3 P^k h_1 c_0) = \rho^3 P^{k+1} h_2.
\]
\end{lem}

\begin{proof}
By \cite[Theorem~4.1]{DI}, the only classes in $\Ext_\ensuremath{\mathbb{R}}[\rho^{-1}]$ that survive the $\rho$-inverted Bockstein spectral sequence are those satifsying $s+f-2w=0$.
Since $\tau P^k h_1$  does not satisfy this equation,
either it supports a Bockstein differential, or
$\rho^r \tau P^k h_1$  is hit by a Bockstein differential for some $r$.
But \autoref{lem:v1-perm} shows that $\tau P^k h_1$ does not
support a differential.
Therefore, $\rho^r \tau P^k h_1$ must be hit by some differential.
By inspection, there is only one possibility.
This establishes the first formula.

The same argument applies to establish the second formula.
\end{proof}

\begin{lem}\label{d4km1Qh14}
For $k\geq 1$,
\[d_{4k-1}\left( \frac{Q}{\rho^{4k-1}} h_1^{4k}\right) = \frac{\gamma}{ \tau^{4k-1}} P^{k-1}h_0^3 h_3.\]
\[d_{4k}\left( \frac{Q}{\rho^{4k}} h_1^{4k+1}\right) = \frac{\gamma}{ \tau^{4k}} P^{k}h_1.\]
\end{lem}

\begin{pf}
The class $ \frac{Q}{\rho^{4k}} h_1^{4k}$ restricts to a class of the same name in $\Ext$ over $\ensuremath{\mathcal{A}}(1)^{C_2}$. There,
we have
\[
d_{4k}\left( \frac{Q}{\rho^{4k}} h_1^{4k}\right) =
\frac{\gamma}{\tau^{4k}} b^k
\]
by \cite[Proposition~7.9]{GHIR}.
The second formula follows immediately from this, since
$P^k h_1$ restricts to $h_1 b^k$.

On the other hand,
the Bockstein class $\frac{\gamma}{\tau^{4k}} b^k$ is not in the image of the restriction from the $\rho$-Bockstein spectral sequence over $\ensuremath{\mathcal{A}}^{C_2}$,
so $\frac{Q}{\rho^{4k}} h_1^{4k}$ must support a shorter Bockstein differential over $\ensuremath{\mathcal{A}}^{C_2}$. The claimed differential is the only possibility.
\end{pf}

\autoref{tbl:BockDiff} summarizes the key differential calculations
that we will need later.

\begin{table}[ht]
\captionof{table}{Some Bockstein differentials}
\label{tbl:BockDiff}
\begin{center}
\begin{tabular}{LLLLLl}
\hline
c &  (s,f, w) & \text{element} & r & d_r & proof \\ \hline
-2k-2 & (1, 0, 2k+3) & \frac{\gamma}{\rho \tau^{2k+1}} & 1 &
	\frac{\gamma}{\tau^{2k+2}} h_0 & \autoref{lem:gamma-diff} \\
-4k-3 & (2, 0, 4k+5) & \frac{\gamma}{\rho^2 \tau^{4k+2}} & 2 &
	\frac{\gamma}{\tau^{4k+3}} h_1 & \autoref{lem:gamma-diff} \\
4k+6 & (8k+7, 4k+4, 4k+1) & \tau^3 P^k h_0^3 h_3 & 3
	& \rho^3 \tau P^{k+1} h_1 & \autoref{lem:d3} \\
4k+6 & (8k+9, 4k+4, 4k+3) & \tau^3 P^k h_1 c_0 & 3
	& \rho^3 P^{k+1} h_2 & \autoref{lem:d3} \\
0 & (8k, 4k-1, 8k) & \frac{Q}{\rho^{4k-1}} h_1^{4k} & 4k-1
  & \frac{\gamma}{\tau^{4k-1}}P^{k-1}h_0^3 h_3 & \autoref{d4km1Qh14} \\
0 & (8k+2, 4k, 8k+2) & \frac{Q}{\rho^{4k+1}} h_1^{4k+1} & 4k
  & \frac{\gamma}{\tau^{4k}}P^kh_1 & \autoref{d4km1Qh14} \\
\hline
\end{tabular}
\end{center}
\end{table}

\section{$\Ext_{C_2}$ in coweight $0$}

\autoref{prop:ExtC2-BL} explicitly describes a large part
of $\Ext_{C_2}$ in coweight $0$.  This result is more
easily understood in the $\Ext$ chart in \autoref{ExtChart},
where we are considering only elements above the shaded region.

\begin{prop}
\label{prop:ExtC2-BL}
In degrees $(s,f,w)$ satisfying $s- w = 0$ and
$f > \frac{1}{2} s - 1$,
$\Ext_{C_2}$ consists of the following classes:
\begin{enumerate}
\item
$h_0^k$ for $k \geq 0$.
\item
$\rho^j h_1^k$ for $j \geq 0$ and $k \geq 0$.
\item
$\frac{Q}{\rho^j} h_1^{4k + \varepsilon}$,
with $k \geq 1$, $0 \leq \varepsilon \leq 3$ and:
\begin{enumerate}
\item
$j \leq 4k - 2$ when $\varepsilon = 0$.
\item
$j \leq 4k - 1$ when $1 \leq \varepsilon \leq 3$.
\end{enumerate}
\end{enumerate}
\end{prop}

\begin{proof}
\autoref{lem:ExtR-mw0} explains how the classes
$h_0^k$ and $\rho^j h_1^k$ arise in $\Ext_\ensuremath{\mathbb{R}}$.
It remains to study $\Ext_{NC}$.

The desired elements of the form $\frac{Q}{\rho^j}{h_1^{4k+\varepsilon}}$
arise from the differentials in \autoref{d4km1Qh14}.
\autoref{prop:rho-cofree-diff} implies that
these elements cannot be involved in any further differentials.

There are several additional Adams periodic families of elements in the
Bockstein $E_1$-page that lie above the line $f = \frac{1}{2}s-1$
in coweight $0$.  All of these elements are the
targets of Bockstein differentials,
as shown in \autoref{tbl:BockDiff-MW1},
so they do not appear in $\Ext_{C_2}$.
Each differential in \autoref{tbl:BockDiff-MW1} follows from the
Leibniz rule and the differentials in \autoref{tbl:BockDiff}.

However, the last three calculations are not entirely obvious.
In these cases, write
\[
\frac{\gamma}{\rho^2 \tau^{4k+1}} P^k c_0 =
\frac{\gamma}{\rho^2 \tau^{4k+2}} \cdot \tau P^k c_0
\]
\[
\frac{\gamma}{\rho^3 \tau^{4k+1}} P^k h_0^3 h_3 =
\frac{\gamma}{\rho^3 \tau^{4k+4}} \cdot \tau^3 P^k h_0^3 h_3
\]
\[
\frac{\gamma}{\rho^3 \tau^{4k+1}} P^k h_1 c_0 =
\frac{\gamma}{\rho^3 \tau^{4k+4}} \cdot \tau^3 P^k h_1 c_0,
\]
and then apply the Leibniz rule to these products.

\begin{table}[ht]
\captionof{table}{Some Bockstein differentials in coweight $1$}
\label{tbl:BockDiff-MW1}
\begin{center}
\begin{tabular}{LLLLL}
\hline
c &  (s,f, w) & \text{element} & r & d_r \\
\hline
1 & (8k+3, 4k+1, 8k+2) & \frac{\gamma}{\rho^2\tau^{4k-2}}P^k h_1 & 2 &
	\frac{\gamma}{\tau^{4k-1}} P^k h_1^2 \\
1 & (8k+4, 4k+1, 8k+3) & \frac{\gamma}{\rho\tau^{4k-1}}P^k h_2 & 1 &
	\frac{\gamma}{\tau^{4k}} P^k h_0 h_2 \\
1 & (8k+4, 4k+2, 8k+3) & \frac{\gamma}{\rho\tau^{4k-1}}P^k h_0h_2 & 1 &
	\frac{\gamma}{\tau^{4k-1}} P^k h_1^3 \\
1 & (8k+8, 4k+2, 8k+7) & \frac{\gamma}{\rho\tau^{4k+1}} P^k h_0 h_3 & 1 &
	\frac{\gamma}{\tau^{4k+2}} P^k h_0^2h_3 \\
1 & (8k+8, 4k+3, 8k+7) & \frac{\gamma}{\rho\tau^{4k+1}} P^k h_0^2 h_3
	& 1 & \frac{\gamma}{\tau^{4k+2}} P^k h_0^3h_3 \\
1 & (8k+10, 4k+3, 8k+9) & \frac{\gamma}{\rho^2\tau^{4k+1}} P^k c_0 & 2 &
	\frac{\gamma}{\tau^{4k+2}} P^k h_1 c_0 \\
1 & (8k+10, 4k+4, 8k+9) & \frac{\gamma}{\rho^3\tau^{4k+1}} P^k h_0^3h_3
	& 3 & \frac{\gamma}{\tau^{4k+3}} P^{k+1} h_1 \\
1 & (8k+12, 4k+4, 8k+11) & \frac{\gamma}{\rho^3\tau^{4k+1}}P^k h_1c_0
	& 3 & \frac{\gamma}{\tau^{4k+4}} P^{k+1} h_2 \\
\hline
\end{tabular}
\end{center}
\end{table}
\end{proof}

\autoref{ExtChart} also shows some
classes that are not part of the Bredon-Landweber region.
These classes arise in the shaded part of the chart.
The structure there is quite complicated, and it will be analyzed
in a range in future work.

\section{The Adams spectral sequence}
\label{sctn:Adams}

We  show in \autoref{NoAdamsDiffs}
that the entire Bredon-Landweber region described in
\autoref{prop:ExtC2-BL} survives
the $C_2$-equivariant Adams spectral sequence.

\begin{prop}\label{NoAdamsDiffs} No element
listed in \autoref{prop:ExtC2-BL}
is either the target or the source of an Adams differential.
\end{prop}

\begin{pf}
Except for the elements $h_0^k$,
all of the classes in the Bredon-Landweber region are $h_1$-periodic. Therefore, any class supporting an Adams differential into the Bredon-Landweber region would be $h_1$-periodic and of coweight 1.
\autoref{h1PerNC1} shows that there are no such classes.

On the other hand,
the elements $h_0^k$ are $h_0$-periodic.
By inspection in low dimensions, there are no
$h_0$-periodic elements that could support differentials whose
values are $h_0^k$.

Adams differentials on the classes in the Bredon-Landweber region
lie in the vanishing region of \autoref{NegConeVanish}.
Therefore, these classes must be permanent cycles.
\end{pf}

\autoref{AdamsEinfChart} shows the Bredon-Landweber region
in the Adams $E_\infty$-page.
Similarly to the Adams $E_2$-page in \autoref{ExtChart}, there are additional classes
in the shaded part of the chart that we will consider in future work.

We now analyze hidden $\rho$ extensions in the Bredon-Landweber region.
The key tool is \autoref{prop:rho-div}.

\begin{prop}
\label{prop:rho-div}
The kernel of the underlying homomorphism
$U: \pi^{C_2}_{n,k} \longrightarrow \pi_n$ is
the image of $\rho:\pi^{C_2}_{n-1,k-1} \longrightarrow \pi^{C_2}_{n,k}$.
\end{prop}

\begin{proof}
This follows immediately from the
cofiber sequence
\[
(C_2)_+ \longrightarrow S^{0,0} \xrtarr{\text{ }\rho\text{ }} S^{1,1},
\]
using the free-forgetful adjunction between
equivariant homotopy classes $(C_2)_+ \longrightarrow X$ and
classical homotopy classes from $S^0$ to the
underlying spectrum of $X$.
\end{proof}

\begin{lem}
\label{lem:rho-hidden}
For $k \geq 4$,
there is a hidden $\rho$ extension from
$Q h_1^k$ to $h_1^k$.
\end{lem}

\begin{proof}
The classical Hopf map $\eta_\mathrm{cl}$ in $\pi_1$
is the underlying map of the equivariant Hopf map
$\eta$ in $\pi^{C_2}_{1,1}$.
Let $k \geq 4$.
Since $(\eta_\mathrm{cl})^k=0$ in $\pi_k$,
\autoref{prop:rho-div} implies that $\eta^k$ must be
a multiple of $\rho$.
The only possibility is that there is a hidden
$\rho$ extension from $Q h_1^k$ to $h_1^k$.
\end{proof}

The hidden extensions of \autoref{lem:rho-hidden} appear in
\autoref{AdamsEinfChart} as dashed lines of negative slope.

In principle, it would be possible for there to be additional
hidden $\rho$ extensions whose sources lie in the
shaded part of \autoref{AdamsEinfChart} and whose targets
lie in the Bredon-Landweber region.
\autoref{lem:no-rho-extn} eliminates this possibility.

\begin{lem}
\label{lem:no-rho-extn}
There are no hidden $\rho$-extensions
whose targets lie in the Bredon-Landweber region.
\end{lem}

\begin{pf}
We use the unit map $S^{0,0} \longrightarrow ko_{C_2}$ for the $C_2$-equivariant connective real $K$-theory spectrum, as studied in \cite{GHIR}.
The entire Bredon-Landweber region is detected
by the Adams $E_\infty$-page for $ko_{C_2}$.
Therefore, a hidden $\rho$ extension into the
Bredon-Landweber region would be detected by a
$\rho$ extension in the homotopy of $ko_{C_2}$.

There are in fact some elements in the homotopy of
$ko_{C_2}$ that support $\rho$ extensions into the image of
the Bredon-Landweber pattern.
These elements are detected by
$\frac{Q}{\rho^{4k-4}} h_1^{4k-1}$ and
$\frac{Q}{\rho^{4k-1}} h_1^{4k}$ in the
Adams $E_\infty$-page for $ko_{C_2}$.
We must show that they do not lie in the image
of the unit map.

These elements support infinite towers of $h_0$-multiplications
in the Adams $E_\infty$-page for $ko_{C_2}$.
Therefore, they cannot lie in the image of the unit map, since
the Adams $E_\infty$-page in \autoref{AdamsEinfChart}
does not include elements that support infinite towers
of $h_0$-multiplications in the relevant degrees.
\end{pf}

We have now collected enough results to prove \autoref{mainthm},
\autoref{cor:fixed-pt-image}, and \autoref{etadiv2}.

\begin{proof}[Proof of \autoref{mainthm}]
The $C_2$-equivariant element $\eta^k$ is detected by $h_1^k$ in the
Adams $E_\infty$-page.
\autoref{lem:rho-hidden} and \autoref{prop:ExtC2-BL} give a lower bound on the
power of $\rho$ that divides $\eta^k$, and
\autoref{lem:no-rho-extn} gives an upper bound on this
power of $\rho$.
\end{proof}

\begin{proof}[Proof of \autoref{cor:fixed-pt-image}]
The geometric fixed points homomorphism
$\phi$ takes the values
$\phi (\rho) = 1$ and $\phi(\eta) = -2$.
(The minus sign in $\phi(\eta)$ depends on choices of orientations
and is inconsequential to the proof.)
Consequently, $\phi(\alpha) = 0$ if $\alpha$ is not $\rho$-periodic,
i.e., if $\alpha$ is detected below the Bredon-Landweber region.

The corollary now follows from \autoref{mainthm}.
\end{proof}

\begin{proof}[Proof of \autoref{etadiv2}]
Recall \cite[Section 8]{DI} that $h_0$ detects $2+\rho\eta$.
Therefore, for homotopy classes detected by $h_0$-torsion classes,
multiplication by $2$ coincides with multiplication by $-\rho\eta$.
Thus $\eta^k$ is divisible by $2^m$ if and only if
$\eta^{k-m}$ is divisible by $\rho^m$.
This latter condition can be determined by \autoref{mainthm}.
\end{proof}

\section{The Mahowald invariant of $2^k$}

The goal of this section is compute the Mahowald invariant of $2^k$
for all $k \geq 0$.
We begin by determining the values of the underlying homomorphism
$U: \pi^{C_2}_{s,s} \longrightarrow \pi_s$
on classes in the Bredon-Landweber region.
\autoref{prop:rho-div} implies that $U(\alpha)$ is necessarily zero
when $\alpha$ is divisible by $\rho$.  On the other hand,
$U(\alpha)$ is always non-zero when $\alpha$ is not divisible by $\rho$.

\begin{thm}
\label{thm:U}
The underlying homomorphism $U: \pi^{C_2}_{s,s} \longrightarrow \pi_s$
takes values as described in \autoref{tbl:u0}.
\end{thm}

\begin{table}[ht]
\captionof{table}{Some values of the underlying homomorphism}
\label{tbl:u0}
\begin{center}
\begin{tabular}{LLL}
\hline
s & \alpha \text{ detected by} & U(\alpha) \text{ detected by} \\
0 & 1 & 1 \\
1 & h_1 & h_1 \\
2 & h_1^2 & h_1^2 \\
3 & h_1^3 & h_1^3 \\
8k-1 & \frac{Q}{\rho^{4k-2}} h_1^{4k} & P^{k-1} h_0^3 h_3 \\
8k+1 & \frac{Q}{\rho^{4k-1}} h_1^{4k+1} & P^k h_1 \\
8k+2 & \frac{Q}{\rho^{4k-1}} h_1^{4k+2} & P^k h_1^2 \\
8k+3 & \frac{Q}{\rho^{4k-1}} h_1^{4k+3} & P^k h_1^3 \\
\hline
\end{tabular}
\end{center}
\end{table}

\begin{proof}
The underlying homomorphism $U$ induces a map of Adams spectral sequences $\Ext_{C_2} \longrightarrow \Ext_{\mathrm{cl}}$.
This map of spectral sequences detects
the first four values in \autoref{tbl:u0}.

The summand $\Ext_{NC}$ lies in the kernel of this map
of spectral sequences.
Therefore,
if $\alpha$ in $\pi^{C_2}_{n,k}$ is detected by an element
of $\Ext_{NC}$, then $U(\alpha)$ must be detected in strictly
higher Adams filtration.
For each of the last four entries in \autoref{tbl:u0},
there is only one possible element in higher Adams filtration.
\end{proof}

The underlying homomorphism $U$ plays a central role in the Mahowald root invariant \cite{MR}.
The Bruner-Greenlees  formulation of the
Mahowald invariant of a homotopy class is given as follows \cite{BG}.
 First, recall that the geometric fixed points homomorphism
$\phi:\pi^{C_2}_{n,k} \longrightarrow \pi_{n-k}$ gives rise to an isomorphism
\cite[Proposition 7.0]{AI}
\[
\pi^{C_2}_{*,*} [\rho^{-1}] \cong
\pi_* [\rho^{\pm 1}].
\]
Then  the Mahowald invariant is defined via the diagram
\[
\begin{tikzcd}
\pi_* \arrow[r]
\arrow[dr, dashed, bend right=5] \arrow[ddr,  bend right=10, "M(-)"'] & \pi^{C_2}_{*,*}[\rho^{-1}] \\
 & \pi^{C_2}_{*,*} \arrow[u] \arrow[d, "U"] \\
  & \pi_*
\end{tikzcd}\]
Here the dashed arrow picks out an element
from the largest possible stem.
More precisely, given a classical stable homotopy class
$\alpha$ in $\pi_k$,
one first chooses an equivariant stable homotopy class $\beta$
in $\pi^{C_2}_{n, n-k}$
such that $\phi(\beta) = \alpha$ and such that $n$ is as large
as possible.  In particular, $\beta$ is not divisible by $\rho$,
for otherwise $n$ would not be maximal.
Then $M(\alpha)$ contains the element $U(\beta)$.
Beware that there can be more than one choice for $\beta$,
so the Mahowald invariant has indeterminacy in general.

Our $C_2$-equivariant calculations allow us to easily
recover the Mahowald invariants  of $2^k$ (\cite[Theorem~2.17]{MR}).

\begin{thm}
\label{thm:Mahowald-inv}
Let $k = 4j + \varepsilon \geq 4$ with $0 \leq \varepsilon \leq 3$.
The Mahowald invariant of $2^k$ contains
an element that is detected by
\begin{center}
\begin{tabular}{ll}
$P^{j-1} h_0^3 h_3$ & if $\varepsilon = 0$. \\
$P^j h_1$ & if $\varepsilon = 1$. \\
$P^j h_1^2$ & if $\varepsilon = 2$. \\
$P^j h_1^3$ & if $\varepsilon = 3$.
\end{tabular}
\end{center}
\end{thm}

\begin{pf}
\autoref{mainthm} determines the value of $\beta$ in the
Bruner-Greenlees formulation of the Mahowald invariant.
Then \autoref{thm:U} gives the value of $U(\beta)$.
\end{pf}

\begin{rem}
\label{rem:MI-indet}
The indeterminacy in $M(2^k)$ is determined by the values
of the underlying map on classes that are detected in the
shaded region of \autoref{AdamsEinfChart}.
It does not seem possible to predict the indeterminacy
of $M(2^k)$ in general.  However, inspection in low dimensions
shows that the indeterminacy of
$M(2^5)$ is generated by
elements detected by $h_1^2 h_3$ and $h_1 c_0$,
while $M(2^k)$ has no indeterminacy for all other values of
$k \leq 8$.
This indeterminacy calculation depends on a detailed analysis of the shaded
region of \autoref{AdamsEinfChart} and will be justified in
future work.
\end{rem}

\section{Charts}

Here is a key for reading the charts of \autoref{ExtChart}, \autoref{AdamsEinfChart}, and \autoref{AdamsEinfkoChart}:
\begin{enumerate}
\item
Blue dots indicate copies of $\ensuremath{\mathbb{F}}_2$ from $\Ext_\ensuremath{\mathbb{R}}$ (or $\Ext$
over the $\ensuremath{\mathbb{R}}$-motivic version of $\ensuremath{\mathcal{A}}(1)$
in the case of \autoref{AdamsEinfkoChart}).
\item
Gray dots indicate copies of $\ensuremath{\mathbb{F}}_2$  from $\Ext_{NC}$ (or
from the negative cone part of $\Ext$ over $\ensuremath{\mathcal{A}}^{C_2}(1)$
in the case of \autoref{AdamsEinfkoChart}).
\item
Horizontal lines indicate multiplications by $\rho$.
\item
Dashed lines of negative slope indicate $\rho$ extensions that
are hidden in the Adams spectral sequence.
\item
Vertical lines indicate multiplications by $h_0$.
\item
Vertical arrows indicate infinite sequences of multiplications
by $h_0$.
\item
Lines of slope $1$ indicate multiplications by $h_1$.
\end{enumerate}

\begin{figure}[h]
\caption{}
\label{ExtChart}
\printpage[name=MW0]
\end{figure}

\begin{figure}
\caption{}
\label{AdamsEinfChart}
\printpage[name=MW{0}Einf]
\end{figure}
\medskip

\begin{figure}
\caption{}
\label{AdamsEinfkoChart}
\printpage[name=MW{0}Einfko]
\end{figure}

\vfill


\section{\label{sec:intro} Introduction:}
Magnetic confinement approaches to fusion rely on the ordered topology of nested magnetic surfaces to prevent the plasma from escaping. Devices such as tokamaks and stellerators are designed towards this; however, in reality, magnetic fields are not perfect and are subject to error fields. These change the magnetic topology and result in magnetic islands appearing at rational surfaces. \par
Magnetic islands are characterised by flat density and temperature profiles due to enhanced transport through them. This reduction of pressure gradients suppresses bootstrap current within the island, generally making the island grow until saturation\cite{Chang,LaHaye97,ZohmEPS97,Gates97}. As a result, the confinement ability of the system decreases, paired with the occurrence of so called neoclassical tearing modes (NTMs)\cite{Buttery_2000}. \par
NTMs were recognised as a source of major disruptions in experiments such as JET\cite{devries,devries14}, and thus their stabilisation  is central. Amongst proposed stabilisation approaches, driving current\cite{reiman,Yoshioka} at the islands with RF waves has stimulated a long list of added efforts,\cite{Sauter2004,Kamendje2005,LaHaye2006,Lahaye08,Henderson08,Lazzari,Volpe,Sauter2010,Bertelli2011,Hennen2012,Smolyakov_2013,Ayten,Borgogno,Volpe2015,Fevrier2016,Wang2015,JC_Li_2017,
Grasso_JPP2016,Grasso18,Poli15}
including many experimental demonstrations\cite{Bernabei,Warrick2000,Gantenbein,Zohm2001,Isayama2000,LaHaye2002,Petty04}. By driving current at the centre of the island using ECCD\cite{fisch80,karney81} or LHCD\cite{fisch78,karney79} one may balance the lack of bootstrap current, and prevent the island growth.\cite{Lazzari} This technique is, however, limited to its application to islands of smaller size due to available power constraints. This makes driving current precisely at their centre difficult, as the deposition width is comparable to the island size. \par
It has been recently suggested\cite{Reiman18} that some of these stringent requirements may be relaxed due to the so called "RF current condensation effect". This effect takes into consideration the non-linear feedback of temperature variations\cite{Westerhof07} resultant from RF wave heating onto the deposition itself. It was found that condensation could improve mitigation as well as reduce radial sensitivity.
 \par
Formally, in Ref.~[\onlinecite{Reiman18}] this non-linear feedback was modeled using a simplified diffusion energy balance equation that included resonant power deposition for a prescribed profile. Yet, the lack of a dissipation mechanism (eg. radiation) or the unlimited absorbable power from RF waves, gave, as they observed, an nonphysical temperature blow-up beyond a bifurcation point. Here that model is extended to include the damping of the RF wave.  \par
In what follows, this extended model is first introduced in detail. Next, the equation is analytically and numerically solved and a hysteresis effect related to the island heating is described. Having introduced this phenomena, the effects of the non-linearities on different RF deposition schemes are explored. In order to evaluate these effects fairly, comparisons are made to the analytic linear solution, which is taken as representative of current approaches.

\section{Fundamental theoretical model}

The model describes both the temperature variations of a magnetic island and the RF wave deposition. The latter may be described as a wave that is being damped along a ray trajectory, so as to provide the plasma with energy. The temperature of the plasma, driven by the RF waves, is described as part of an energy balance model with thermal diffusion. To construct such a model, various approximations are introduced. \par
 Take, as a starting point, the transport equation\cite{Braginski}, representing the second moment of the Boltzmann equation, to describe temperature, $T$, dynamics. Under the assumption of no significant flows, and considering fast equilibration between electrons and ions, a single equation may be written combining the two-fluid Braginskii equations
\begin{equation}
\frac{3}{2}nk_B\partial_t T-\nabla\cdot(\mathbb{\kappa}\nabla T)=P \label{eqn::difEqn}
\end{equation}
where $n$ is the plasma density, $\kappa$ is the generalised heat conductivity tensor, and $P$ a volumetric power deposition which will be later related to the RF power input. \par
It has been stated as an assumption for Eq.~(\ref{eqn::difEqn}) that electrons and ions are effectively equilibrated. Thus, the equation should only apply to those time scales larger than the typical equilibration time $\tau_\mathrm{eq}$; this constitutes the first temporal constraint: $t\gg\tau_\mathrm{eq}$, where $t$ denotes the time scales that the model is suited to describe.
\par
A second point is related to the use of temperature as a measure of plasma energy. The concept of temperature customarily applies only to thermalised systems, in which the populations in $v$-space are Maxwellian distributed. However, the continuous injection of RF waves distorts the distribution function locally so as to deposit energy resonantly onto a small fraction of faster electrons (eg. for LH waves $v\approx 4.5v_{Te}$ and for EC $v\approx 3v_{Te}$)\cite{fisch87,karney81}. Hence, the plasma is made up of a Maxwellian bulk with a well defined $T$, but also a resonant minority population. If the bulk $T$ is to represent the total internal energy of the plasma, then the energy drawn locally in $v$-space needs to be redistributed quickly. Here a second time ordering is introduced: collisional thermalisation ($\tau_\mathrm{se}$) and isotropisation ($\tau_\perp$) processes must be faster than the time scales of interest ($t\gg\tau_\mathrm{se},\tau_\perp$). In that case non-Maxwellian features that affect a minute fraction of the total population may be ignored to leading order. Kinetic details will still prove important for $P$. \par
As it stands, Eq.~(\ref{eqn::difEqn}) remains a three dimensional problem, but it may be cast into an approximate reduced 1D problem assuming the following. Consider as it happens (see Table \ref{tab::timeScales})\cite{Spitz53,Westerhof07}, transport over a given magnetic flux surface to be much faster than perpendicular to it; that is, $f_\kappa\equiv\kappa_\parallel/\kappa_\perp\gg1$. In that case, magnetic surfaces will be approximately isothermal, simplifying the derivatives in the diffusive term (i.e. the second term in Eq.~(\ref{eqn::difEqn})). For analytic simplicity, the geometrical factor associated with the particular shape of island flux surfaces will be ignored, using instead a single slab coordinate $x$ (one may think of making a cut to an elongated, narrow island). This reduced form is,
\begin{equation}
\frac{3}{2}nk_B\partial_t T-\partial_x(\mathbb{\kappa_\perp}\partial_x T)=P \label{eqn::dif2}
\end{equation}
As a result of this slab adaptation, areal weighting is made equal for all points, though in reality this should be larger for the edges. A more complete treatment considering the flux coordinate is left for future work, though it was shown itn [\onlinecite{Reiman18}]  that the slab geometry shared the qualitative physics with the more realistic geometry.\par
	 Nonetheless, the model retains an important feature of the magnetic island geometry: the closed nature of magnetic surfaces about the island centre. Consequently, in the reduced 1D model,  temperature solutions are required to be even about the centre $x=0$. \par
Proceed now to linearise Eq.~(\ref{eqn::dif2}). Assume that the temperature changes ($\widetilde{T}$) in the island due to directing RF heating ($P$) to it are small; ie. $\epsilon=\widetilde{T}/T_0\ll 1$, where $T_0$ is the equilibrium temperature. Seeking precision in this definition, $T_0$ is defined to be the island temperature when the RF power is not aimed directly at the island, but is rather part of the total power budget that heats the centre of the tokamak (see Fig.~\ref{fig::figBry}). Because $T_0$ is constant over the island (as is the density), any island inhomogeneity that develops will be at least $O(\epsilon)$; and with this ordering we drop derivatives with respect to $\kappa_\perp$ in Eq.~(\ref{eqn::dif2}), \par
\begin{equation}
\frac{3}{2}nk_B\partial_t \widetilde{T}-\kappa_\perp\partial^2_x \widetilde{T}=P \label{eqn::difEqn1D}
\end{equation}
The consideration of stiff temperature profiles that could modify $\kappa_\perp$ non-smoothly is left for future work.
\par
	To drop the time dependence of Eq.~(\ref{eqn::difEqn1D}), the time scales of concern should exceed those of energy diffusion ($\kappa_\perp\partial^2_x \widetilde{T}$ term) and the driving times ($P$ term). Consider the former; clearly for the dominance of the diffusion term $ t\gg\tau_D=W_i^2/\chi_\perp$. Here $W_i$ represents the island width and $\chi_\perp\approx\kappa_\perp/nk_B$ is the heat diffusion coefficient. \par
For the latter, $nk_B\partial_t T\sim P$ suggests that the time derivative may be dropped provided the power density of the wave has had enough time to deposit all the needed thermal energy; i.e., $ t\gg \tau_E= nk_B\Delta T/P$, where $\Delta T\sim\epsilon T_0
$ is the characteristic variation of the island temperature.  Then, with $t$ in this regime, and being consistent with all of previous requirements, \par
\begin{equation}
-\kappa_\perp\partial^2_x \widetilde{T}=P \label{eqn::difEqn1DSimp}
\end{equation}
Prior to detailing the form of $P$, the problem should be closed by both defining the spatial domain and setting appropriate boundary conditions. Naturally, one defines the last closed surface of the island including the X-points as boundaries of the domain of $x$, ie. $|x|\leq W_i/2$. The island width, $W_i$, will be kept constant. That means that $t$ must be shorter than the typical island growth $\tau_i=\left(\partial \ln W_i/\partial t\right)^{-1}\approx\tau_\eta$, where $\tau_\eta$ is the global resistive time scale.\cite{Ruther73} Table \ref{tab::timeScales} shows this last requirement is consistent with previous time orderings. Width changes may then be treated adiabatically; ie. the steady state Eq.~(\ref{eqn::difEqn1DSimp}) may be taken to be satisfied at all times as $W_i$ is changed artificially. \par

\begin{table}
\begin{tabular}{cc}
    Time scale & s \\\hline \hline
    $\tau_\mathrm{se}$ & $1\times 10^{-3}$ \\
    $\tau_\perp$ & $1\times 10^{-3}$ \\
    $\tau_D$ & $3\times 10^{-3}$ \\
    $\tau_E$ & $8\times 10^{-3}$\\
    $\tau_\mathrm{eq}$ & $4\times 10^{-1}$\\
    $\tau_i,\tau_\eta$ & $2\times10^2$\\
 \\
\end{tabular}
\begin{tabular}{cc}
    Dimensionless & scales  \\\hline \hline\\
    $f_\kappa$ & $10^{13}$ \\
    $\epsilon$ & $10^{-1}$ \\
\\
\\
\\
\\
\end{tabular}
\caption{Summary of the various scales relevant to obtain the form of Eqs. (\ref{eqn::approx1}) and (\ref{eqn::approx2}) in typical tokamak parameters of $T=10~$keV, $a=1~$m, $R=5~$m, $W_i=0.05~$m, $n=10^{14}~$cm$^{-3}$, $\chi_{e\perp}=0.5~$m$^2$s$^{-1}$, $Z=1$, $q=2$.}
\label{tab::timeScales}
\end{table}
\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.4\textwidth]{bryTemp.png}
\caption{Schematic sequence of temperature profiles of the full plasma for: a) No island present b) An island present at $R=R_i$ c) RF energy deposition displaced to within the island. The reddish band represents the region of RF deposition, while the gray one corresponds to other heating sources.}\label{fig::figBry}
\end{figure}

The boundary condition on temperature encodes the influence of the island on the remaining of the plasma (and vice versa). To specify it, a simplified treatment of the energy dynamics of the rest of the tokamak is done, using a steady state diffusion model like Eq.~(\ref{eqn::difEqn}). It is convenient to apply Gau{\ss}' theorem to magnetic flux surfaces, so that $\int_\psi PdV=-\partial_\psi T\int_{\partial \psi}\kappa_t\hat{n}\cdot\nabla\psi~dS$, where $\psi$ is the magnetic flux coordinate. This shows that the slope of the temperature profile at a particular flux surface is determined by the power deposited inside it. For simplicity, let us associate a spatial 1D coordinate $R$, in the absence of islands monotonic with $\psi$, and take the heat conductivity, $\kappa_t$, to be constant; then, the slope of temperature at some $R_0$ is determined by the power deposited at $R<R_0$. With this in mind, let there be some heating in the tokamak centre: $P_c$ (fusion power, Ohmic heating, etc.) and $P_\mathrm{RF}$ (RF heating). The temperature profile is then determined by the heat flux and the fixed plasma edge temperature (see Fig.~\ref{fig::figBry}a). \par

Now, let there be an island of size $W_i$ at a distance $R_i$ from the core over which the temperature profile is flat (see Fig.~\ref{fig::figBry}b). Because heat sources have not changed, the temperature slope remains unchanged elsewhere. Let then $P_\mathrm{RF}$ be redirected to the island (ie. the case of interest). For those magnetic flux surfaces at $R>R_i+W_i/2$, the enclosed total power does not change, and thus the slope of $T$ should neither (see Fig. \ref{fig::figBry}c). Given that the tokamak plasma edge temperature is fixed, the temperature of at  the edge of the island, $T_0$, remains unchanged.
The boundary condition for our island temperature may then be taken to be $\widetilde{T}(x=\pm W_i/2)=0$.\par

It is now the turn of specifying $P$ in Eq.~(\ref{eqn::difEqn1DSimp}) to represent the energy deposition from RF waves. Adopting a geometrical optics (GO) description of the wave envelope\cite{Tracy}, the evolution of the energy of the wave may be written as
\begin{equation}
d_t\Bar{V}=[-(\nabla\cdot v_g)+\frac{\omega_t}{\omega}+2\gamma]\Bar{V}\approx 2\gamma \Bar{V}, \label{eqn::waveEnrgy}
\end{equation}
where $\Bar{V}$ represents the wave energy density, $v_g$ is the group velocity of the wave, $\omega_t$ represents the time derivative of frequency due to a time dependent medium, $d_t$ represents the total time derivative following a wave along a ray and $\gamma$ represents collisionless damping rate. Assuming the medium to be stationary in the wave damping time scale ($\omega_t/\omega\ll \gamma$) and the spatial inhomogeneity to be much smaller than the variation resultant from the damping ($\nabla\cdot v_g\ll\gamma$), the last approximated equality follows. This condition is not difficult to satisfy considering only small variations are created within the island. \par
Expressing Eq.~(\ref{eqn::waveEnrgy}) in terms of $x$, the distance along the ray $d_tx=v_g$,
\begin{equation}
d_x\Bar{V}=2\frac{\gamma}{v_g}\Bar{V}, \label{eqn::waveEnerSimp}
\end{equation}
which has the form of damped propagation.
The factor $\gamma$ may be obtained under the assumption of a Maxwellian magnetised background\cite{Stix,karney81}, one may show that for EC and LH, $\gamma\propto \exp(-\chi^2)$ where $\chi=(\omega-n\Omega_e)/kv_{Te}=v_\parallel/v_{Te}$ where $v_\parallel$ is the phase velocity of the wave, $v_{Te}$ is the electron thermal speed, $\Omega_e$ is the electron cyclotron frequency, and $n=0$ corresponds to LH waves and $n=-1$ to EC waves. The power deposition both for electron cyclotron and lower hybrid waves occurs on the tail of the Maxwellian velocity distribution, with damping exponentially small in the lowest resonant velocity.  \par
It is this exponential factor which makes deposition highly sensitive to variations in temperature. Indeed, considering the phase velocity of the wave to remain constant over the extent of the island,
 $$\gamma\propto e^{-(v_\parallel/v_{Te_0})^2}\exp\left(\frac{v_\parallel^2\widetilde{T}}{v_{Te,0}^2T_0}\right)\rightarrow \frac{2\gamma}{v_g}=-\alpha e^u.$$
 It is convenient here to define the dimensionless variable $u\equiv v_\parallel^2\widetilde{T}/v_{Te,0}^2T_0=w^2\widetilde{T}/T_0$, where $w=v_\parallel/v_{T_e,0}$.
The location of the wave damping within the whole plasma is generally dependent on $T$, $B$, $\omega$ and $k$. This location can be determined using ray tracing\cite{Prater08,bonoli1986simulation}. Here the picture is simplified by artificially restricting the damping to a particular defined region within the island , while keeping $v_\parallel$ constant, and that way allowing for the $\widetilde{T}$ expansion. A fully self consistent, full GO analysis is left for future work.
\par
Let us express Eq.~(\ref{eqn::waveEnerSimp}) as
\begin{equation}
d_x\Bar{V}=-\alpha(x)e^u\Bar{V} \label{eqn::damp}
\end{equation}
where the damping strength $\alpha(x)\equiv2\alpha_0f(x)/W_i$. The factor $\alpha_0$ represents the strength of the damping, but it is also defined in a dimensionless way to include the island width $W_i$.  For example, for EC waves\cite{karney81} $\alpha_0\approx W_i\sqrt{\pi}\omega_{pe}^2 \exp(-w^2)/2ckv_{Te}$. It is helpful to introduce a more physically motivated interpretation of $\alpha_0$.  If a linear limit is taken of Eq.~(\ref{eqn::damp}), the power deposition profile takes the form $|V'|\propto \exp\left(-2\alpha_0x/W_i\right)$; i.e., $\alpha_0$ is the ratio of the island half-width to the characteristic deposition width. Note that the deposition has an exponential shape, and not its usual Gaussian form generally considered for electron-cyclotron waves\cite{Prater08}; however, both schemes are peaked and of finite width, and ultimately quite similar.  \par
Now, going back to the original question: how is $P$ related to this wave energy $\Bar{V}$? From the damping of the wave along a ray, it is easily seen that the volumetric power deposition at a given point is given by $P=-d_t\Bar{V}\approx -v_gd_x\Bar{V}\equiv -v_g\Bar{V}'(x)$. But, because in this particular geometry the points $\pm x$ are linked together (recall this is true due to them belonging to the same flux surface) these points share the total deposition at $x$ and $-x$. All things considered,
\begin{subequations}
\begin{align}
-\kappa_\perp\partial^2_x \widetilde{T}(x)=-v_g\frac{\Bar{V}'(x)+\Bar{V}'(-x)}{2} \label{eqn::approx1}\\
\Bar{V}'(x)=-\alpha(x)e^u\Bar{V}(x)\label{eqn::approx2}
\end{align}
\end{subequations}

These equations may be non-dimensionalised, reducing them to
\begin{subequations}
\begin{align}
V'(\widetilde{x})&=-f(\widetilde{x})e^uV(\widetilde{x})  \label{eqn::WaveEnergy} \\
u''&=\frac{V'(\widetilde{x})+V'(-\widetilde{x})}{2} \label{eqn::diffus}
\end{align}
\end{subequations}
where the new $V(\widetilde{x})=\bar{V}(x)W_iv_gv_p^2/2\alpha_0\kappa_\perp T_0v_T^2=\bar{V}(x)W_i\Upsilon^2/\alpha_0$ and $\widetilde{x}=2\alpha_0x/W_i$. Note that the edges are now at $x=\pm\alpha_0$. \par
In order to complete the setting of the problem, an initial value must be taken for Eq.~(\ref{eqn::WaveEnergy}). Let $V_X\equiv V(-\alpha_0)\equiv V_0/\alpha_0$, where $V_0$ is a constant representing some wave energy density input. For interpreting solutions, it is important to bear in mind that $V_X$ is independent of island width, but will however scale as $1/\alpha_0$ with the deposition strength. \par
For clarity in the following Sections III and IV, one may refer to Appendix B as a quick reference for the variables employed.

\section{Hysteresis phenomena}

To investigate the effect of the non-linear wave deposition, consider the tractable basic problem of wave damping occurring everywhere within the island. This case, represented by $f(\widetilde{x})=1$, allows for an analytic solution of Eqs.~(\ref{eqn::WaveEnergy}) and (\ref{eqn::diffus}) (see Appendix A for a detailed derivation). Implementing the appropriate boundary and initial conditions,
\begin{equation}
u(\widetilde{x})=2\log\gamma -\log\left[\sqrt{(\lambda+1)^2-\gamma^2}\cosh\gamma\widetilde{x}+(\lambda+1)\right] \label{eqn::tempProf}
\end{equation}
where the paramters $\lambda$ and $\gamma$ are determined by,
\begin{subequations}
\begin{align}
\left[\gamma^2-(1+\lambda)\right]^2=\cosh ^2(\alpha_0\gamma) \left[(\lambda+1)^2-\gamma^2\right] \label{eqn::param1} \\
\gamma^2=(2\lambda+1)+(V_X-\lambda)^2 \label{eqn::param2}
\end{align}
\end{subequations}
 \par
The integration constant is $\lambda=(V_X+V_f)/2$, where $V_f$ is the energy density when exiting the island. Eqs. (\ref{eqn::param1}-\ref{eqn::param2}) solve $\lambda$ implicitly, which ultimately determines the temperature of the island as a function of $\alpha_0$ and $V_X$ from Eq.~(\ref{eqn::tempProf}). \par

The dependence of perturbed central island temperature on these two parameters will be represented as contour curves (see for example Fig.~\ref{fig::fig2}). Two main representations are of particular physical interest. First, contours of constant $\alpha_0$ (ie. fixed deposition strength and island width) in the $u(0)-V_0$ plane. These contours show the effects of the wave power on temperature (see Fig.~\ref{fig::fig2}a). The second interesting picture is related to how the heating of the island evolves as its width or the wave profile width changes. This is captured by curves of constant $V_X$ or $V_0$ at fixed deposition strength in the $u(0)-\alpha_0$  plane (see Fig.~\ref{fig::fig2}b for an example).  \par

\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.5\textwidth]{hysteresisExamples.jpg}
\caption{a) Island temperature at fixed deposition strength for different island widths, as a function of absorbed RF power, showing the appearance of the bifurcation point (shown bigger points). Dashed lines correspond to the linear limit of the solution; the dotted line to the asymptotic form of the solution for almost complete power deposition.
b) Central island temperature with varying island width for constant incident wave energy densities $V_X$. The broken lines represent the solution to the linear problem, while the dotted ones the asymptotic limit solution for complete energy deposition.}
\label{fig::fig2}
\end{figure}

There are a number of general features in the solutions to Eqs.~(\ref{eqn::tempProf})-(\ref{eqn::param2}) worth highlighting. The first of those is the existence of bifurcation points. As previously observed\cite{Reiman18}, for sufficiently broad depositions (small $\alpha_0$), saddle-node bifurcations appear, at which the two lower temperature solutions disappear. Such points are marked in Fig.~\ref{fig::fig2}a. The appearance of these points may be linked to the action of a self-focusing mechanism affecting RF waves. Schematically, below the bifurcation, significant energy leakage takes place. As the bifurcation is approached, island temperature perturbations become larger, while the power deposited by the wave increases accordingly. This positive feedback eventually extracts all RF energy effectively, reaching a higher temperature steady state and thus jumping into an upper branch of the solution. \par

This hot stable solution may be seen in Fig.~\ref{fig::fig2}, along with the asymptotic form of the solution as $V_f\rightarrow 0$ (dotted line). The proximity of the two solutions demonstrates that the upper branch indeed corresponds to nearly complete deposition of the wave energy in the island. Such a solution branch is also, immediately after the bifurcation point, significantly larger than the linear prediction (see broken lines). \par

A consequence of the solution structure obtained is the hysteresis behaviour of island heating. To illustrate such a process,
 take as a starting point the system to be in equilibrium at the lower temperature branch in Fig.~\ref{fig::fig2}a, and increase the absorbed energy of the incoming wave ($V_0$) gradually. As a result, the temperature of the island will grow until the bifurcation point is reached. Once at this point, and driven by the self-focusing feedback, the temperature of the island will rapidly increase  towards the upper branch, which is the only stable solution at high absorbed power.\par
 The hot island exhibits, at this point, a large temperature difference between the centre and the separatrix. A priori, this would help to absorb RF power closer to the O-point, and thus also drive thte central current\cite{fisch87}. Here we ignore current associated with the DC electric field, which could be due to Spitzer conductivity or due to the hot electron conductivity\cite{fisch85a,karney_fisch_jobes}. These currents are less important than the directly driven RF current.\cite{Reiman18} With such presumed centred current drive, the island would tend to shrink and stabilise, as governed by the Rutherford equation\cite{Ruther73}. This size reduction corresponds to a leftwards displacement towards smaller values of $\alpha_0$ in Fig.~\ref{fig::fig2}b. In such a case, and if the energy available to the island is maintained, the reduction does not imply a return back to the original low temperature, but instead remains in the more effective current driving upper branch for some time. Similarly, once in the upper branch, driving power requirements are relaxed, and lower $V_0$ would still keep the plasma hot. This constitutes the hysteresis effect. A more careful discussion on the usefulness, accessibility and consistency of this sketched simplified picture for particular deposition schemes is the concern of following sections. \par

 \par

\section{Effects of non-linear feedback}

\subsection{Typical parameters}
	Before proceeding further, a brief estimate and collection of typical values for both $V_0$ and $\alpha_0$ is presented. We emphasise that our slab model provides a physical qualitatively correct picture of the problem, but only a rough guide into the quantitative behaviour of the more realistic geometry, as previous calculations suggest\cite{Reiman18}. Other simplified features, such as the exponential form of the linear deposition profile, are also different when compared to actual experiments\cite{Prater08}, but do however share the fundamental characteristics. Thus, the linear case will be taken as reference in guiding conclusions in the following sections, as well as orientative comparison standard to existing experimental parameters.    \par
	Focus first on the values for the parameter $\alpha_0$. Two different routes are taken at this point. One possible method uses, given the definition of $\alpha_0$ as the size of the linear power deposition width, typical deposition widths in tokamak experiments could be used to obtain $\alpha_0\sim 0.5-3$.\cite{Volpe15} It has been recently reported that current drive profiles are in experiment subject to broadening\cite{Brook17} by factors of 2--3 due to effects unaccounted for in ray tracing routines, such as edge density fluctuations. This effective broadening could make typical $\alpha_0$ values even lower, down to $\sim 0.2$.   \par
	Alternatively, one could use the form for $\alpha_0$ given before and obtained in the context  of GO. In the case of ECCD, for instance, using typical approximated hydrogen tokamak values (see caption of Table \ref{tab::timeScales}), with $w^2\sim10$,\cite{karney81} wavenumber\cite{Prater04} $k\sim2\pi/(5~\mathrm{mm})$, density $n\sim10^{20}$~m$^{-3}$ and $W_i\sim10~$cm; $\alpha_0\approx W_i\sqrt{\pi}\omega_{pe}^2 \exp(-w^2)/2ckv_{Te}\sim10^{-1}$.
\par
	Now consider the wave power density $V_0$. First, we may compare existing literature\cite{Westerhof07} where island temperature variations are computed to our linear toy model. Comparing values of $u(0)$,orientative typical power parameters on the order of $V_0\sim5$ for 20~MW RF power are found. In a more first principle approach, we might use the rescaled definition of $V_0$ introduced before. Using an RF power on the order of $P\sim10$~MW, with a beam of cross section $A\sim1~$m$^2$, with $\chi_\perp\sim1$m$^2$/s, temperature $T_0\sim10$~keV, density $n\sim10^{20}$~m$^{-3}$ and $w^2\sim10$, then $V_0\approx PW_iw^2/2A\chi_\perp nk_BT_0\sim10^1$. \par
	Summarising:
\begin{equation}
	\alpha_0\sim0.1-3 \hspace{1cm}  V_0\sim 10^1
\end{equation}

\subsection{Central deposition}
	The spatial distribution of the RF deposition strongly affects the final temperature of the island, as well as the island mitigation efficiency. In this section, the best case scenario is first analysed; i.e., deposition starting from the island centre. To formally emulate this ideal case, $f(\widetilde{x})=H(\widetilde{x})$, where $H$ is the Heaviside step function. Given this newly introduced asymmetry, well defined parity is lost from the equations and the solution to the equation is only found numerically. \par
	Consider first the occurrence of bifurcation points, and in particular, how they depend on $\alpha_0$. To illustrate these points, Figures \ref{fig::fig3} and \ref{fig::fig4} are presented. \par
	In the broad linear deposition limit, with $3\alpha_0\ll1$, there always exists a bifurcation point (see Fig.~\ref{fig::fig3}). The turning point, however, occurs at increasingly larger wave energy densities $V_0$. One may understand this result by referring to the analytic asymptotic form of the solution at low $V_0$. In that limit, the system takes the form of the linear problem, for which $u(0)\sim V_0\left(\alpha_0-1+e^{-\alpha_0}\right)/2\alpha_0$. This shows that the temperature of the island becomes decreasingly responsive as $\alpha_0\rightarrow 0$ (see the decreasing initial slopes of curves in Fig.~\ref{fig::fig3} inset), which is ultimately related to there being a significant wave energy leakage ($V_f=V_Xe^{-\alpha_0}$). \par
	Only for those cases for which the initial energy leakage $V_f$ is significant will a bifurcation occur. The bifurcation is a result of the system being able to access all that previously lost energy when the damping $e^u$ factor becomes significant. This jump will be associated with a narrowing of the deposition and a current that is more efficiently utilised in stabilisation. \par

\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.5\textwidth]{centrDepoBifPoint.jpg}`
\caption{Wave power density value for the bifurcation at a given deposition strength $\alpha_0$, for deposition starting at the island centre. The larger scatter point represents the limiting value of $\alpha_0$ over which no bifurcation exists. The inset shows curves of constant deposition strength in the $u(0)-V_0$ plane. The broken line represents the asymptotic form of the non-linear solution.}
\label{fig::fig3}
\end{figure}

\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.5\textwidth]{centrDeposDiffNLLinComplete.jpg}
\caption{a) Wave power density value for 50\% difference between the linear and non-linear solutions for centred wave energy deposition. b) Upper bound to ratio of non-linear to linear island temperature. c) Location of middle point of the wave deposition $x_\mathrm{cent}|V(x_\mathrm{cent}\alpha_0)=(V_X+V_f)/2$ for the points in the curve in a) (continuous line), and middle point for the linear equivalent problem (dotted line).  Solutions for the region $\alpha_0$ to the left of the broken line represent solutions with existing bifurcation points.}
\label{fig::fig4}
\end{figure}

	As the initial deposition is reduced by increasing $\alpha_0$, $V_f$ in the lower branch decreases, and it eventually becomes too small to sustain a bifurcation. Fig.~\ref{fig::fig3}) shows the boundary value $\alpha_0=0.32$ beyond which no bifurcation occurs. \par
	Where no bifurcation occurs, the temperature of the island only undergoes a smooth transition in temperature between the linear solution and the high temperature asymptote (see Fig.~\ref{fig::fig3}). That limiting form of the nonlinear solution as $V_f\rightarrow0$ is $u(0)\sim V_0/2$ (see dashed line in Fig.~\ref{fig::fig3} inset and Appendix C), which is also the limit as $\alpha_0\rightarrow\infty$ of the linear deposition. That is, the non-linear response serves as a short cut via self-focusing to the linear ideal infinitely narrow deposition.  Thus, for the case of centralized deposition, the nonlinear mechanism always leads to an enhanced temperature increase.
 \par
	However, the wave power required to obtain a substantial improvement exhibits strong dependence on $\alpha_0$ as shown in Fig.~\ref{fig::fig4}a. The large energy leakage and the small difference between linear and non-linear solutions at low and high $\alpha_0$ respectively leave a most easily accessible (lower $V_0$) central region at values $\alpha_0\sim1-1.5$. To emphasise the second of these limitations Fig.~\ref{fig::fig4}b shows the ratio of the analytic asymptotic forms of the non-linear and linear solutions. Evidently, the differences become marginal (ie. the ratio tends to one) for stronger depositions, which explains why Fig.~\ref{fig::fig4}a diverges at $\alpha_0\sim2.8$. \par
	Finally, we examine the extent to which self focusing narrows the RF deposition profile. From Fig.~\ref{fig::fig4}c the self-focusing effect is apparent, and will undoubtedly improve stabilisation by bringing current drive closer to the O point. This reduction by a factor of $\sim2-3$ in width opens the possibility of previously disregarded regimes of island stabilisation. Note also that for $\alpha_0\approx2.8$, the non-linear deposition profile becomes similar to that of a delta function (within discretisation). \par
	This analysis suggests that the region of interest and current experimental relevance may in some subset of cases (for the broadest depositions) show some hysteresis behaviour, but most will just show significant temperature variation. In addition, as a result of deposition narrowing due to the non-linearity, islands could be stabilised when traditionally predicted not to. This opens the door to experimental verification of the non-linear effect, as well as extension of mitigation schemes.  \par

\subsection{Edge deposition}
	The scenario adopted for analysis before was that of central deposition. This is the ideal case, and so it presumes that one is experimentally capable of aiming perfectly at the centre of the island without depositing any energy before that. But, what would happen if the deposition departs from this idealised case? The worst case scenario is now presented. To that end, we recover the analytic solution from Eqs. (\ref{eqn::tempProf})-(\ref{eqn::param2}) which represents wave deposition from the very edge of the magnetic island. \par
With this analytic result at hand, let us explore first the limit of complete wave deposition: $u(0)\sim \log\left(V_0/2\alpha_0+1\right)$. The linear solution gives $u(0)=V_0e^{-\alpha_0}(\cosh\alpha_0-1)/\alpha_0$. It is remarkable that the non-linear model gives a logarithmic growth of the island temperature as the power input is increased, while the linear case grows linearly. It necessarily follows that some non-linear inhibition mechanism must take part. Indeed, one may relate this to the deposition profile becoming localised ever closer to the edge of the island. \par
 	It is the same self-focusing that narrowed the deposition closer to the centre when central deposition was considered, which displaces deposition towards the island edge (see Fig.~\ref{fig::fig6}). Physically, the RF wave becomes so strongly damped that it runs out of energy very close to the edge. There, the temperature slope is large, and thus as heat gets to the centre of the island, it is quickly lost across the edges, $u(0)$ becoming limited. For typical values $\alpha_0\sim1$ and $V_0\sim5-10$, the peak of deposition $x_\mathrm{peak}\sim-0.7$; ie. power is deposited somewhere between the X- and O-points. \par

\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.5\textwidth]{normalDeposChngStrngth.jpg}
\caption{Power at which the non-linear and linear solutions give the same central temperature with deposition starting from the edge, as a function of the deposition strength. The region to the left of the curve represent the case for larger non-linear solution. The inset shows examples of constant deposition curves as a function of power. The broken lines represent the linear solution, while the dotted curves show the logarithmic asymptotic behaviour of the non-linear one. }
\label{fig::fig5}
\end{figure}

\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.5\textwidth]{normalDeposWidthAnalysisPk.jpg}
\caption{Location of the power deposition peak as a function of power for various deposition strengths (complementary to Fig.~\ref{fig::fig5}). A value of $x_\mathrm{peak}=0$ corresponds to a centred deposition, while a value of $-1$ represents the edge.}
\label{fig::fig6}
\end{figure}
	As a result of this detrimental displacement, only over a limited regime will the non-linear solution be hotter than the linear one. As shown in Fig.~\ref{fig::fig5}, this interval is larger for the broader depositions, but tends to disappear as $\alpha_0\rightarrow\infty$. When the deposition is broader, the non-linear focusing increases the amount of deposited power significantly, a benefitial addition that outweights the inhibiting displacement for a more extended range of powers.  \par
	In addition to temperature, this deposition shift will also bring the driven current closer to the X-point. This displacement can be catastrophic when trying to mitigate the growth of magnetic islands. The proximity of the deposition to the edge may be seen in Fig.~\ref{fig::fig6}. The plot shows that indeed, for large $V_0$ values, the non-linear self-focusing brings the deposition ever closer to the edge. This is the result of the large damping $e^u$ which drains the incoming wave faster than the linearly increasing $V_0$. Nevertheless, there is an initial region in which the focusing is starting to affect the system, and clearly becomes beneficial in terms of drawing power towards the centre. \par

	To further explore the sensitivity of stability on $x_\mathrm{peak}$, the procedure in [\onlinecite{Lazzari}] is followed. The relative Fourier weighting to $\Delta '$ in the Rutherford equation due to driving current at a particular flux coordinate, $\psi$ ($\Omega$ in the reference), may be estimated taking the 1D spatial variable $x$ in our model to match the spatial $x$ in [\onlinecite{Lazzari}], looking at $\xi=0$ (see Fig.~\ref{fig::fig7}). The calculation shows that within approximately 90\% of the island the drive is stabilising. Therefore, looking back at Fig.~\ref{fig::fig6}. The current drive will still be central enough to be stabilising for a significant fraction of the cases, even when the non-linear solution is colder than the linear one. No definitive conclusion may however be drawn on the precise fraction of the island that is truly stabilising, as the treatment of the island geometry in [\onlinecite{Lazzari}] is different from that of our model. A fully consistent treatment is left for future work.  \par
	We have shown that an initially broad RF profile may then be used to stabilise islands, so long as the input power remains below some upper bound. This result is promising and an idea to further explore. \par

\begin{figure}
\hspace*{-0.3cm}
\includegraphics[width=0.5\textwidth]{CDcontribDeltaMinus.jpg}
\caption{Scaled contribution to $\Delta'$ in the Rutherford equation due to current drive at different positions within the island. Current drive is stabilising for $\Delta'<0$, ie. when driving is roughly within 90\% of the island extent. The inset shows the deposition profiles for the linear and non-linear solutions corresponding to the point shown in Fig.~\ref{fig::fig6}, represented in log scale (the linear deposition is almost uniform, but largest at the X point). }
\label{fig::fig7}
\end{figure}

\par
More generally, if deposition was to start midway between the X- and O-points, then deposition would be driven and narrowed towards that location instead. The trends will then fit between the extreme centre and edge cases shown explicitly, with less constraining requirements as the centre is reached. In this context, RF overshooting scenarios will never suffer from the inhibition that takes place when undershooting with increasing $V_0$.  \par
Yet another possibility to circumvent inhibition might be to look for means to amplify the wave power within the island, such as through an $\alpha$-channeling effect\cite{fisch1992interaction}. This could give a power not peaked at the island periphery. Unfortunately, this volumetric amplification does not occur for electron cyclotron waves, but might be exploited when using lower hybrid waves\cite{ochs_2015a}. Despite lacking this possible enhancement, EC waves do however have the benefit of having a $B$-dependent resonance that allows for deposition starting at a particular point in space.

\section{Conclusion}
The possibility of hysteresis involving the heating of magnetic islands with RF waves is shown for a symmetrised, 1D slab model. The wave power deposition self focuses mediated by island temperature, leading to higher temperatures. Past a bifurcation point, the island may remain in this high temperature solution even as it shrinks or power is reduced. Exploiting the hysteresis effect could thus provide an easier and improved way to eliminate magnetic islands. \par
In typical parameter regimes of current experiments, though, bifurcation is likely to occur only for the broadest profiles. Yet, centred deposition scenarios show that important temperature increments (on the order of $\sim20\%$) occur. These differences are most significant in broader deposition schemes, where the self-focusing mechanism of the non-linear model is most different from the linear model. Alongside these thermal variations, the non-linear narrowing of the RF deposition profile by a factor of up to $\sim2-3$ for typical values will improve the utility of current drive for purposes of stabilisation.\par
Broad schemes that deviate from the centre still lead to not dissimilar self-narrowing and stabilisation under certain circumstances. These circumstances involve the form of EC power density profiles, and must be considered in designing deposition scenarios. In particular, for deposition profiles that peak before the O-point, there is a threshold power density above which a self-inhibition mechanism is encountered; beyond the O-point, this threshold does not exist. This opens the door to exploring previously disregarded broad RF deposition stabilisation schemes.  \par

\begin{acknowledgments}
Thanks to Suying Jin for fruitful discussions and help. \par
This work was supported by US DOE DE-AC02-09CH11466 and DE-SC0016072.

\end{acknowledgments}