\section{Introduction}
We study the dynamics of magnetic \emph{domain walls} (DWs) in nanowires for a constant applied field and current aligned with the material anisotropy. The dynamics of a magnetization field $\boldsymbol{m}=(m_1,m_2,m_3)^{\mkern-1.5mu\mathsf{T}}\in\mathbb{S}^2$ in case of an applied spin-polarized current is governed by the \emph{Landau-Lifschitz-Gilbert-Slonczewski} equation
\begin{equation*}\label{eq:LLGS}
\partial_t \boldsymbol{m} - \alpha\boldsymbol{m}\times\partial_t\boldsymbol{m}=-\boldsymbol{m}\times\boldsymbol{h}_\textnormal{eff}+\boldsymbol{m}\times\left(\boldsymbol{m}\times \boldsymbol{J}\right). \tag{LLGS}
\end{equation*}
Here, $\alpha>0$ is the \emph{Gilbert} damping parameter and $\boldsymbol{J}$ the applied current. The magnetization $\boldsymbol{m}$ as well as the effective field $\boldsymbol{h}_\textnormal{eff}$ are normalized by the spontaneous magnetization~\cite{berger1996emission, gilbert2004phenomenological, hubert2008magnetic, landautheory, slonczewski1996current, slonczewski2002currents}.

The effective field is defined as the negative variational derivative of the total magnetic free energy $\boldsymbol{E}(\boldsymbol{m})$~\cite{abert2019micromagnetics, hubert2008magnetic}. We consider uniaxial anisotropy, with material parameter $\mu<0$, applied field $h\in\mathbb{R}$, both aligned along the $\boldsymbol{e}_3$-axis, and the simplified form
\[\boldsymbol{h}_\textnormal{eff}\coloneqq -\frac{\delta \boldsymbol{E}}{\delta \boldsymbol{m}}=\partial_x^2 \boldsymbol{m}+h\boldsymbol{e}_3-\mu m_3\boldsymbol{e}_3.\]
In the considered axially symmetric configuration, the applied current is given by
\begin{equation}\label{eq:Current}
\boldsymbol{J}\coloneqq \frac{\beta}{1+c_{\rm cp} m_3}\boldsymbol{e}_3,
\end{equation}
where $\beta\ge0$ and $c_{\rm cp}\in(-1,1)$ describe the current density and the ratio of spin-polarization, respectively~\cite{bertotti2008spin, bertotti2005magnetization, mayergoyz2009nonlinear}. We provide more details of the modeling assumptions in \S~\ref{sec:Preliminaries}.

In contrast to the well-known \emph{Landau-Lifschitz-Gilbert}~\eqref{eq:LLG} equation, where $\boldsymbol{J}\equiv 0$ and which can be written as the gradient of the total magnetic free energy, the generalized~\eqref{eq:LLGS} equation is no longer variational. This makes parts of the analysis more challenging, but all results presented in this paper include the~\eqref{eq:LLG} equation by setting $c_{\rm cp}=0$.

The basic states of~\eqref{eq:LLGS} are the constant up- and down-magnetization states $\boldsymbol{m}=\pm\boldsymbol{e}_3$, which form spatially homogeneous equilibria of~\eqref{eq:LLGS}. We are interested in coherent structures connecting these states, i.e., magnetic \emph{domain walls} (DWs), and their dynamics. This fundamentally depends on the stability properties of $\pm\boldsymbol{e}_3$, and we refer to the situation in which both asymptotic states $\pm \boldsymbol{e}_3$ are dynamically stable as \emph{bistable}, while \emph{monostable} refers to the case in which only of $\pm\boldsymbol{e}_3$ is stable.

In the considered axially symmetric setting, coherent structures occur as relative equilibria with respect to translation as well as rotation and are solutions of the form
\begin{equation}\label{eq:CoherentStructureAnsatz}
\boldsymbol{m}(\xi,t)=\boldsymbol{m}_\ast(\xi)\mathrm{e}^{\mathbf{i} \varphi(\xi,t)},\quad \xi=x-st \quad\textnormal{and}\quad \varphi(\xi,t)=\phi(\xi)+\Omega t,
\end{equation}
with \emph{speed} $s\in\mathbb{R}$ and \emph{frequency} $\Omega\in\mathbb{R}$. Here the complex exponential acts on $\boldsymbol{m}_\ast\in\mathbb{S}^2$ by rotation around the $\boldsymbol{e}_3$-axis. DWs are solutions of the form~\eqref{eq:CoherentStructureAnsatz} with asymptotics $\boldsymbol{m}_\ast(\xi)\to\pm\boldsymbol{e}_3$ for $\xi\to\pm\infty$ or $\xi\to\mp\infty$.

The only known explicit DW solutions, denoted $\boldsymbol{m}_0$, occur at $\beta\cdot c_{\rm cp}=0$, where \eqref{eq:LLGS} reduces to \eqref{eq:LLG}. This family of solutions was discovered in~\cite{goussev2010domain} for applied field with arbitrary strength, cf.\ Figure~\ref{fig:DWsSphere} (a), and will be discussed in \S~\ref{sec:ExplicitFamily}.

\begin{figure}
	\centering
	\begin{subfigure}[b]{0.3\textwidth}
		\includegraphics[trim=0 0 0 150, clip, width=\textwidth]{figures/atansphere.eps}
		\caption*{(a)}
	\end{subfigure}
	\hfill
	\begin{subfigure}[b]{0.3\textwidth}
		\includegraphics[trim=0 0 0 150, clip, width=\textwidth]{figures/freezedsmallfield.eps}
		\caption*{(b)}
	\end{subfigure}
	\hfill
	\begin{subfigure}[b]{0.3\textwidth}
		\includegraphics[trim=0 0 0 150, clip, width=\textwidth]{figures/freezedlargefield.eps}
		\caption*{(c)}
	\end{subfigure}
	\caption{Snapshots of DWs projected onto $\mathbb{S}^2$. (a) explicit homogeneous DW profile given by~\eqref{eq:AnalyticHomogeneousSolution}, independent of $h$. (b)-(c) profiles of inhomogeneous DWs obtained by direct long-time simulations of the full PDE~\eqref{eq:LLGS} for $\alpha=1,\,\beta=0.75,\, \mu=-1$, and $c_{\rm cp}=0.5$ with $h=1$ (bistable regime) in (b) and $h=10$ (monostable regime) in (c). Numerical asymptotic speeds and frequencies are $s=0.081$ and $\Omega=0.923$ in (b) as well as $s=3.868$ and $\Omega=8.965$ in (c), where in (c) the linear predictions in the monostable regime are $s^\textnormal{lin}=3.873$ and $\Omega^\textnormal{lin}=9$.}
	\label{fig:DWsSphere}
\end{figure}

A first classification of more general DWs is based on the \emph{local wavenumber} $q\coloneqq\phi'(\xi)$, where a DW with $q\equiv 0$ is called \emph{homogeneous} and \emph{inhomogeneous} otherwise~\cite{siemer2020}; in particular $\boldsymbol{m}_0$ is homogeneous. It is known that DW solutions to~\eqref{eq:LLGS} for $c_{\rm cp}\neq 0$ cannot be homogeneous~\cite{Melcher2017}. An additional classification of DWs is based on the convergence behavior of the local wavenumber $q$ as $\xi\to\pm\infty$: we refer to DWs as \emph{flat}, if $|q(\xi)|$ has a limit in $\overline{\mathbb{R}}$ for $|\xi|\to\infty$, and \emph{non-flat} otherwise; homogeneous DWs are flat by definition.

Utilizing perturbations methods, it has been proven in~\cite{siemer2020} that non-flat DWs also exists in case $c_{\rm cp}=0$, and more importantly that flat as well as non-flat DWs exists for $|c_{\rm cp}|\ll 1$ sufficiently small, where the applied field strength can be arbitrary. Examples of inhomogeneous DWs for different applied fields and $c_{\rm cp}$ away from zero are presented in Figures~\ref{fig:DWsSphere} (b) and (c), which were obtained by direct long-time simulations of~\eqref{eq:LLGS} starting from a Heaviside step function initial condition in (b) and a compact perturbation of the unstable state $-\boldsymbol{e}_3$ in (c).

In this paper, we study stability properties and the selection of DWs analytically and numerically. The selection mainly refers to the monostable case for which -- in contrast to the bistable case -- a family of DWs coexists for given parameters. Our main interest lies in understanding selection mechanisms and determining conditions on spreading speeds and frequencies of DWs for an arbitrary strength of the applied field, both in the bistable as well as monostable parameter regime.

Within the bistable regime, DWs can be orbitally asymptotically stable, i.e., perturbations relax back to a translate of the DW. This occurs for stable point spectrum (up to translation as well as rotation modes), which has been proven in~\cite{goussev2020dynamics} for field strength $h$ below a threshold. Here, we slightly improve this threshold analytically and show that stability is robust for $0<|c_{\rm cp}|\ll 1$. Moreover, we discuss the selection of spreading speed and its sign by heuristically extending the micromagnetic free energy~\eqref{eq:FreeEnergy}. Together with numerical simulations of the full PDE and numerical continuation in the $h$-$c_{\rm cp}$ parameter plane for standing DWs, i.e. $s=0$, we corroborate that: given any $\alpha>0, \beta\ge 0,\mu<0$, and $c_{\rm cp}\in(-1,1)$, the selected DW is approximately, though not exactly, stationary for the applied field
\[h=\frac{\beta}{2\alpha c_{\rm cp}}\big( \ln(1+c_{\rm cp})-\ln(1-c_{\rm cp})\big).\]

In the monostable case, one expects that DWs propagate towards the unstable state~\cite{ducrot2019spreading, holzer2014anomalous, van2003front} and thus invade it, as studied in case $c_{\rm cp}=0$ for $\boldsymbol{m}_0$ in~\cite{goussev2020dynamics}.
In order to determine the asymptotic speed and frequency of this invasion process, we follow in detail the approach of linear spreading speeds and pointwise growth, which is based on refined linear stability properties of the unstable state~\cite{holzer2014criteria}. Briefly, one considers the stability problem in spaces with weighted norms, which moves the essential spectrum, thus changing stability properties. A barrier for this motion is the so-called \emph{absolute spectrum}~\cite{sandstede2000absolute} which is defined by the spatial roots of the dispersion relation, that relates the exponential rate in spatial direction with that in time. The absolute spectrum contains in particular those spectral values for which (certain) two spatial roots coincide that are at the same time singularities of the so-called \emph{pointwise Green's function}~\cite{beck2017stability, holzer2014criteria, howard1999pointwise, zumbrun1998pointwise}. This function is the resolvent kernel associated to the linearization around the unstable state and pointwise growth occurs if and only if a singularity has positive real part.

We will present an explicit formula for the absolute spectrum of certain (complex) linear matrix operators, which include the operators that arise from linearizing~\eqref{eq:LLGS} around the invaded state. This also includes the operators arising in the well-known complex Ginzburg-Landau equation (CGL) around the zero state, which is not surprising since~\eqref{eq:LLGS} can be written as a generalized (CGL) via stereographic projection. We will further utilize information about the explicit location of the absolute spectrum in order to locate singularities of the pointwise Green's function and to describe the transition from pointwise decay to pointwise growth of perturbations localized around the unstable state. This will enable us to determine the linear spreading speed and the associated linear spreading frequency, which are predictors for the dynamically selected DWs in~\eqref{eq:LLGS}. While the linear spreading speed is determined by the transition from pointwise decay to pointwise growth, the linear spreading frequency is not associated to a selection in the sense of stability. We substantiate our analytic results via direct numerical simulations, which corroborate that the explicit family is selected also inside the monostable region, but only until the linear spreading speed exceeds that of the explicit family, cf.\ Figure~\ref{fig:FreezingZeroCCPLight}. From the perspective of linear versus nonlinear speed selection, these DWs are so-called \emph{pushed} fronts in the monostable regime~\cite{paquette1994structural, stokes1976two, van2003front}, since they propagate faster than the linear prediction. Thus, their selection is necessarily driven by nonlinear effects. When the linear spreading speed exceeds the speed of the explicit family, the selected DWs lie outside the explicit family, propagate according to the linear prediction, and hence are so-called \emph{pulled} fronts~\cite{paquette1994structural, stokes1976two, van2003front}. For large applied fields, the explicit family has speeds that are again lower than that the linear spreading speed. However, the explicit family is not selected in this regime due to unstable point spectrum, which we verify numerically.

Our results match the prediction from~\cite{goussev2020dynamics} when $c_{\rm cp}=0$, i.e., for~\eqref{eq:LLG}, which is based on the so-called saddle-point approximation~\cite{ebert2000front, miller2006applied, van2003front}. Our more general approach is based on weighted $L^2$-spaces, the explicit expression of the absolute spectrum associated to the unstable states, and its relation to the pointwise Green's function. Moreover, utilizing the weighted spectrum paves the way to tackle problems beyond those considered in this paper, see e.g.~\cite{faye2020remnant}.

\begin{figure}
	\centering
	\begin{subfigure}[b]{0.49\textwidth}
		\includegraphics[trim=30 0 170 20, clip, width=\textwidth]{figures/zeroccpspeedlight6.eps}
		\caption*{(a)}
	\end{subfigure}
	\hfill
	\begin{subfigure}[b]{0.49\textwidth}
		\includegraphics[trim=30 0 170 20, clip, width=\textwidth]{figures/zeroccpfreqlight6.eps}
		\caption*{(b)}
	\end{subfigure}
	\caption{Numerically observed spreading speed (crosses) in (a) and spreading frequency (circles) in (b) in the absence of spin torque compared with the explicit speed and frequency of $\boldsymbol{m}_0$ (dashed lines) and the linear predictions (solid lines). Parameters are $\alpha=1, \beta=0.75, \mu=-1, c_{\rm cp}=0$ and an applied field $h$ between $\beta/\alpha$ and $\beta/\alpha-3\mu$. Vertical dotted line at $h=1.75$ separates bistable and monostable regions. Pushed fronts are numerically observed for $1.75<h<1.92$ (red shaded region), whereas pulled fronts are observed for $1.92<h$ (blue shaded region).}
	\label{fig:FreezingZeroCCPLight}
\end{figure}

The remainder of the paper is organized as follows. We motivate the problem, describe first properties, introduce the linearized eigenvalue problem, and discuss the stability of the asymptotic states $\pm\boldsymbol{e}_3$ in \S~\ref{sec:Preliminaries}. Next, propagation dynamics of DWs in the \emph{bistable} regime are discussed in \S~\ref{sec:Bistable}, where we also extend the stability threshold for the explicit family analytically and present the heuristic relation between $h$ and $c_{\rm cp}$ for standing DWs. The \emph{monostable} regime is discussed in \S~\ref{sec:Monostable}, where we present the explicit formula of the absolute spectrum for the more general class of linear operators and relate this to the singularities of the pointwise Green's function. From this, we compute the \emph{linear spreading speed} and \emph{linear spreading frequency} and show that non-flat DWs are favored for large applied fields. In \S~\ref{sec:Numerics}, we numerically verify stable point spectrum of $\boldsymbol{m}_0$ into the monostable region. Additionally, we provide some details behind the numerical continuation and simulation results which are presented throughout the paper. Finally, we conclude with a discussion and outlook in \S~\ref{sec:DiscussionAndOutlook}.

\section{Preliminaries on homogeneous domain walls and asymptotic states}\label{sec:Preliminaries}

We start with a brief outline of the model background. The classical equation which describes the magnetization dynamics for a unit vector field $\boldsymbol{m}\in\mathbb{S}^2$ was proposed by L.\ Landau and E.\ Lifschitz~\cite{landautheory} and later extended by T.\ Gilbert~\cite{gilbert2004phenomenological}. For normalized time, it reads
\begin{equation*}\label{eq:LLG}
\partial_t \boldsymbol{m} - \alpha \boldsymbol{m} \times \partial_t \boldsymbol{m}=-\boldsymbol{m}\times \boldsymbol{h}_\textnormal{eff} , \tag{LLG}
\end{equation*}
where $\alpha>0$ is the Gilbert damping parameter. It describes a damped precession of the magnetization $\boldsymbol{m}=(m_1,m_2,m_3)$ around the effective field $\boldsymbol{h}_\textnormal{eff}$, both normalized by the spontaneous magnetization $M_S$~\cite{hubert2008magnetic, lakshmanan2011fascinating}. Based on normalization, time is measured in units of $(\gamma M_S)^{-1}$ ($\gamma$ is the absolute value of the gyromagnetic ration) and space is measured in units of $2A/\mu_0 M_S^2$ ($A$ is the exchange stiffness and $\mu_0$ the permeability)~\cite{hubert2008magnetic}. The effective field itself is defined as the negative variational derivative of the micromagnetic free energy $\boldsymbol{E}(\boldsymbol{m})$ and thus collects the exchange, the anisotropy, the Zeeman or external field, the demagnetization or stray field, and also the magnetoelastic energy~\cite{abert2019micromagnetics}. The leading order contributions of uniaxial anisotropy with parameter $\mu$ and easy axis $\boldsymbol{e}_3$, as well as local approximation of the demagnetization field, without effects caused by mechanical stress or deformation of the ferromagnetic object, lead to free energies of the form
\begin{equation}\label{eq:FreeEnergy}
\boldsymbol{E}(\boldsymbol{m})=\frac{1}{2}\int (\nabla \boldsymbol{m})^2 + \mu m_3^2\, \textnormal{d}x - \int \boldsymbol{h} \cdot \boldsymbol{m} \, \textnormal{d}x,
\end{equation}
where $\boldsymbol{h}\in\mathbb{R}^3$ is the applied field. We choose $\boldsymbol{h}$ to be aligned along the $\boldsymbol{e}_3$ axis, i.e., $\boldsymbol{h}=h\boldsymbol{e}_3$, with a time-independent and uniform field parameter $h\in\mathbb{R}$. Because the anisotropy energy is rotationally symmetric around $\boldsymbol{e}_3$ in our case, \emph{easy-axis} and \emph{easy-plane} anisotropy correspond to $\mu<0$ and $\mu>0$, respectively, according to the energetically preferred directions. Since we are interested in DWs in a nanowire geometry, we make the standing assumption $\mu<0$ throughout this paper.

In order to take effects caused by spin polarized currents into account, J.\ Slonczewski derived a model in~\cite{slonczewski1996current} in which the term $-\boldsymbol{m}\times \boldsymbol{J}$ was added to the effective field $\boldsymbol{h}_\textnormal{eff}$. We consider the axially symmetric case, i.e., rotation symmetric around the $\boldsymbol{e}_3$-axis, which yields a current of the form~\eqref{eq:Current} and where the parameters $\beta\ge 0$ and $c_{\rm cp}\in(-1,1)$ describe the strength of spin-transfer and ratio of polarization in the considered axially symmetric situation~\cite{berger1996emission, bertotti2008spin, mayergoyz2009nonlinear, slonczewski1996current, slonczewski2002currents}.

In one space dimension, where $\boldsymbol{m}=\boldsymbol{m}(x,t), x\in\mathbb{R}$, we consider the~\eqref{eq:LLGS} equation, which in detail reads
\[\partial_t \boldsymbol{m} -\alpha \boldsymbol{m}\times \partial_t \boldsymbol{m}= -\boldsymbol{m}\times \left( \partial_x^2\boldsymbol{m} +(h-\mu m_3)\boldsymbol{e}_3-\frac{\beta}{1+c_{\rm cp} m_3}\boldsymbol{m} \times \boldsymbol{e}_3\right).\]

The~\eqref{eq:LLG} equation, with its Lyapunov (variational) structure, arises as a special case of~\eqref{eq:LLGS} if $\beta=0$ or $c_{\rm cp}=0$, where in the latter case this follows from shifting the applied field $h\to h-\beta/\alpha$~\cite{Melcher2017}. Therefore, while our focus is on the non-variational spintronic case $\beta\cdot c_{\rm cp}\neq 0$, all results presented in this paper that contain the case $\beta\cdot c_{\rm cp}=0$ carry over to~\eqref{eq:LLG}.

Additionally, \eqref{eq:LLGS} possesses a skew variational structure (see for example~\cite{Melcher2017}) as it can be written in the form
\[\partial_t\boldsymbol{m}-\alpha\boldsymbol{m}\times \partial_t\boldsymbol{m}=\boldsymbol{m}\times \frac{\delta \boldsymbol{E}}{\delta \boldsymbol{m}} + \boldsymbol{m} \times \left(\boldsymbol{m} \times \frac{\delta \boldsymbol{\Psi}}{\delta \boldsymbol{m}} \right),\]
with potential for $\boldsymbol{J}$ given by
\begin{equation}\label{eq:Potential}
\boldsymbol{\Psi}(\boldsymbol{m})\coloneqq \frac{\beta}{c_{\rm cp}}\ln\left(1+c_{\rm cp} (\boldsymbol{m}\cdot\boldsymbol{e}_3)\right).
\end{equation}
It is also well-known that~\eqref{eq:LLGS} admits an equivalent form as an explicit evolution equation of quasilinear parabolic type,  which reads
\begin{equation}\label{eq:QuasilinearLLGS}
\partial_t \boldsymbol{m}=D(\boldsymbol{m})\partial_x^2\boldsymbol{m}+ \frac{\alpha}{1+\alpha^2} | \partial_x \boldsymbol{m}|^2 \boldsymbol{m}+f(\boldsymbol{m})
\end{equation}
with $| \partial_x \boldsymbol{m}|^2=(\partial_x m_1)^2 + (\partial_x m_2)^2 + (\partial_x m_3)^2$,
\begin{align*}
D(\boldsymbol{m})&\coloneqq\frac{1}{1+\alpha^2}\begin{pmatrix} \alpha & m_3 & -m_2 \\ -m_3 & \alpha & m_1 \\ m_2 & -m_1 & \alpha \end{pmatrix},\\
f(\boldsymbol{m})&\coloneqq \frac{\boldsymbol{\beta}(m_3)}{1+\alpha^2}\begin{pmatrix} m_1m_3-\alpha m_2 \\ m_2m_3+\alpha m_1 \\ m_3^2-1 \end{pmatrix} -\frac{h-\mu m_3}{1+\alpha^2} \begin{pmatrix}\ \alpha m_1m_3+m_2 \\ \alpha m_2m_3-m_1 \\\alpha \left( m_3^2-1 \right) \end{pmatrix},
\end{align*}
and where we introduced
\begin{align}\label{e:defbeta}
\boldsymbol{\beta}(m_3)=\boldsymbol{\beta}\coloneqq\frac{\beta}{1+c_{\rm cp} m_3} \quad \textnormal{as well as} \quad \boldsymbol{\beta}^\pm\coloneqq \boldsymbol{\beta}(\pm 1)=\frac{\beta}{1\pmc_{\rm cp}}.
\end{align}
For solutions of the form~\eqref{eq:CoherentStructureAnsatz}, which satisfy $\boldsymbol{m}(\pm\infty) =\pm \boldsymbol{e}_3$, or vice versa, we now consider the third component only. Since rotation acts on the first and second component only, we may write $m_3(x,t)=m_3(\xi)$ in a co-moving frame with $\xi=x-st$.

Following~\cite{komineas2019traveling}, for such solutions consider the integral
\[M\coloneqq \int_{-\infty}^\infty m_3(x-st) \,\mathrm{d} x \]
(in the Cauchy principle sense), which relates to the speed $s$ via
\[\frac{d}{dt} M = -s\int_{-\infty}^\infty m_3'\,\mathrm{d} x. \]
Since for a DW $\partial_t m_3=-s m'_3$ with $'=\textnormal{d}/\textnormal{d}\xi$, such a solution to~\eqref{eq:QuasilinearLLGS} fulfills
\[(1+\alpha^2)\partial_t m_3=\alpha m''_3 + m''_1 m_2 - m_1 m''_2 + \alpha|\boldsymbol{m}'|^2 m_3 - (\alpha(h-\mu m_3)-\boldsymbol{\beta}(m_3))(m_3^2-1),\]
and we obtain, via integration by parts,
\begin{align*}
\Delta_3\cdot  s &= -\int_{-\infty}^\infty \partial_t m_3 \,\mathrm{d} x\\
&= \frac{\alpha}{1+\alpha^2}\int_{-\infty}^\infty |\boldsymbol{m}'|^2 m_3 \,\mathrm{d} x + \frac{1}{1+\alpha^2}\int_{-\infty}^\infty \boldsymbol{\beta}(m_3)(m_3^2-1) \,\mathrm{d} x \\
&\quad -\frac{\alpha h}{1+\alpha^2}\int_{-\infty}^\infty m_3^2-1 \,\mathrm{d} x- \frac{\alpha \mu}{1+\alpha^2}\int_{-\infty}^\infty m_3(m_3^2-1) \,\mathrm{d} x.
\end{align*}
This simplifies further for symmetric DWs. In particular for even $m_1, m_2$ and odd $m_3$ the first and last integral vanish. Moreover, for odd $m_3$ we obtain to leading order
\begin{equation}\label{eq:PropagationSignIntegral}
\Delta_3 \cdot s = \frac{\beta - \alpha h}{1+\alpha^2}\int_{-\infty}^\infty m_3^2-1 \,\mathrm{d} x + \mathcal{O}(c_{\rm cp}^2).
\end{equation}
In particular, this defines a curve in the $h$-$c_{\rm cp}$-plane with quadratic tangency at $c_{\rm cp}=0$, cf.\ \S~\ref{sec:StandingDomainWalls}. Moreover, for the sign of propagation we obtain the leading order relation
\begin{equation}\label{eq:velsym}
\mathrm{sgn}(s) = \mathrm{sgn}\big(\Delta_3\big)\cdot  \mathrm{sgn}(\beta - \alpha h).
\end{equation}

\subsection{Coherent structure form and explicit DW}\label{sec:ExplicitFamily}
As we are interested in the dynamics of and near DWs, it is natural to apply the coherent structure ansatz~\eqref{eq:CoherentStructureAnsatz}. In a co-moving frame $\xi=x-st$ with wavespeed $s\in\mathbb{R}$ as well as co-rotating frame with (superimposed) frequency $\Omega\in\mathbb{R}$ around $\boldsymbol{e}_3$ system \eqref{eq:QuasilinearLLGS} turns into
\begin{equation}\label{eq:QuasilinearSpeedFrequencyLLGS}
\partial_t \boldsymbol{m}=D(\boldsymbol{m})\boldsymbol{m}''+(1+\alpha^2)s\boldsymbol{m}'+ R\boldsymbol{m}+\frac{\alpha}{1+\alpha^2} | \boldsymbol{m}'|^2 \boldsymbol{m}+f(\boldsymbol{m}),
\end{equation}
where
\[R\coloneqq (1+\alpha^2)\Omega\begin{pmatrix}
0&1&0\\
-1&0&0\\
0&0&0
\end{pmatrix}.\]
Since $\boldsymbol{m}(\xi,t)$ lies on the unit sphere, it is convenient to rewrite~\eqref{eq:LLGS} in spherical coordinates via
\begin{equation}\label{eq:SphericalCoordinates}
\boldsymbol{m}=\sin(\theta)\left(\cos(\varphi)\boldsymbol{e}_1 + \sin(\varphi)\boldsymbol{e}_2 \right) + \cos(\theta) \boldsymbol{e}_3.
\end{equation}
Applying the coherent structure ansatz~\eqref{eq:CoherentStructureAnsatz} within these coordinates, \eqref{eq:LLGS} transforms into \begin{equation}\label{eq:SphericalCoherentStructurePDE}
\begin{split}
\begin{pmatrix}
\alpha & -1\\
1 & \alpha
\end{pmatrix}
\begin{pmatrix}\dot{\phi}\sin(\theta)\\-\dot{\theta}
\end{pmatrix}&=
\begin{pmatrix}
2\phi ' \theta '\cos(\theta) + s\theta'\\
-\theta '' -\alpha s \theta'
\end{pmatrix}\\
&\quad +\sin(\theta)\begin{pmatrix}
\phi '' +\alpha s \phi'-\alpha\Omega+ \beta/(1+c_{\rm cp}\cos(\theta))\\
(\phi ')^2 \cos(\theta) + s\phi'- \Omega + h-\mu\cos(\theta)
\end{pmatrix}
\end{split}
\end{equation}
with derivatives $\dot{\,}=\,\mathrm{d}/\!\,\mathrm{d} t$ and again $'=\textnormal{d}/\textnormal{d}\xi$. The existence of domain walls asymptotically connecting $\pm\boldsymbol{e}_3$ in~\eqref{eq:LLGS} now translates into the existence of heteroclinic orbits between $\theta=0$ and $\theta=\pi$ in the steady state equation to~\eqref{eq:SphericalCoherentStructurePDE}, which is given by the system of ODEs
\begin{equation}\label{eq:SphericalCoherentStructureODE}
\begin{split}
\theta''&=-\alpha s \theta' + \sin(\theta)\left(s\phi'+(\phi')^2\cos(\theta)+\alpha(h-\mu\cos(\theta))-\Omega\right)\\
\phi''\sin(\theta)&=-s\theta'-2\phi'\theta'\cos(\theta)-\sin(\theta)\left(\alpha s \phi'+\beta/(1+c_{\rm cp} \cos(\theta))-\alpha \Omega\right)
\end{split}
\end{equation}
In case $c_{\rm cp}=0$, the aforementioned explicit family of solutions $\boldsymbol{m}_0$ from~\cite{goussev2010domain, Melcher2017} in the original coordinates reads (recall $\mu<0$)
\begin{equation}\label{eq:AnalyticHomogeneousSolution}
    \boldsymbol{m}_0=\begin{pmatrix}
    \sech\left(\sigma\sqrt{-\mu} \xi\right)\\0\\ -\tanh\left(\sigma\sqrt{-\mu} \xi\right)
    \end{pmatrix}.
\end{equation}
In spherical coordinates, we denote this solution by $(\theta_0, \phi_0)$, where $\phi_0=\textnormal{constant}$ and $\theta_0=2\arctan(\exp(\sigma \sqrt{-\mu} \xi))$. The sign of propagation is, consistent with~\eqref{eq:velsym}, given by $\sigma=\mathrm{sgn}(s\cdot(\alpha h-\beta))$ in terms of speed
\begin{equation}\label{eq:HomogeneousSpeed}
s^2=-\frac{(\alpha h-\beta)^2}{\mu(1+\alpha^2)^2},
\end{equation}
and the frequency is
\begin{equation}\label{eq:HomogeneousFrequency}
\Omega=\frac{h+\alpha\beta}{1+\alpha^2}.
\end{equation}
Note that the shape as well as width of $\boldsymbol{m}_0$ depend on the anisotropy only and that these DWs are homogeneous since $\phi_0$ is constant. We further note that two explicit solutions of the form~\eqref{eq:AnalyticHomogeneousSolution} exist for any value of the applied field: for $s\neq0$ one left- and one right-moving, and for $s=0$ two stationary ones.

\subsection{Eigenvalue Problem} In order to study spectral stability properties of solutions to~\eqref{eq:LLGS}, we now introduce the associated eigenvalue problem. The linearization around a (steady) solution $\overline{m}=(\overline{m}_1(\xi), \overline{m}_2(\xi), \overline{m}_3(\xi))^{\mkern-1.5mu\mathsf{T}}$ of~\eqref{eq:LLGS} can be obtained by substituting $\boldsymbol{m}(\xi,t)=\overline{m}(\xi)+\varepsilon \boldsymbol{n}(\xi)\mathrm{e}^{\lambda t}+\mathcal{O}(\varepsilon^2)$ for $|\varepsilon| \ll 1$ into~\eqref{eq:QuasilinearSpeedFrequencyLLGS}, which leads (at order $\varepsilon$) to the matrix eigenvalue problem
\begin{equation}\label{eq:LinearEigenvalueProblem}
    \mathcal{L}(\overline{m})\boldsymbol{n}=\lambda \boldsymbol{n},
\end{equation}
where $\boldsymbol{n}$ is pointwise tangent to the sphere at $\overline{m}$ and
\[(1+\alpha^2)\mathcal{L}(\overline{m})=\]
\begin{align*}
    &\quad\begin{pmatrix}
    \alpha\partial_\xi^2 & \overline{m}_3\partial_\xi^2-\overline{m}_3'' & -\overline{m}_2\partial_\xi^2+\overline{m}_2''\\ -\overline{m}_3\partial_\xi^2+\overline{m}_3'' & \alpha\partial_\xi^2 & \overline{m}_1\partial_\xi^2-\overline{m}_1'' \\ \overline{m}_2\partial_\xi^2-\overline{m}_2'' & -\overline{m}_1\partial_\xi^2+\overline{m}_1'' & \alpha\partial_\xi^2
    \end{pmatrix}+ (1+\alpha^2)\begin{pmatrix}
    s\partial_\xi &\Omega&0\\
    -\Omega&s\partial_\xi&0\\
    0&0&s\partial_\xi
    \end{pmatrix}\\
    &+\left(\alpha|\overline{\boldsymbol{m}}'|^2+2\alpha\overline{\boldsymbol{m}}\left( \overline{\boldsymbol{m}}'\right)^{\mkern-1.5mu\mathsf{T}}\partial_\xi\right)\textnormal{Id}+\begin{pmatrix}
    \overline{\boldsymbol{\beta}}\,\overline{m}_3 & -\alpha\overline{\boldsymbol{\beta}}\, & \overline{\boldsymbol{\beta}}\,\overline{m}_1-\widetilde{\boldsymbol{\beta}}(\overline{m}_1\overline{m}_3-\alpha\overline{m}_2)
    \\
    \alpha\overline{\boldsymbol{\beta}}\,& \overline{\boldsymbol{\beta}}\,\overline{m}_3 & \overline{\boldsymbol{\beta}}\,\overline{m}_2-\widetilde{\boldsymbol{\beta}}(\overline{m}_2\overline{m}_3+\alpha\overline{m}_1)\\
    0&0 & 2\overline{\boldsymbol{\beta}}\,\overline{m}_3 + \widetilde{\boldsymbol{\beta}}(1-\overline{m}_3^2)
    \end{pmatrix}\\
    &+(h-\mu\overline{m}_3)\begin{pmatrix}
    -\alpha\overline{m}_3 & -1 & -\alpha\overline{m}_1\\
    1 & -\alpha\overline{m}_3 & -\alpha\overline{m}_2\\
    0 & 0 & -2\alpha\overline{m}_3
    \end{pmatrix}
    +\begin{pmatrix}
    0 & 0 & \mu(\alpha\overline{m}_1\overline{m}_3+\overline{m}_2)\\
    0 & 0 & \mu(\alpha\overline{m}_2\overline{m}_3-\overline{m}_1)\\
    0 & 0 & \alpha\mu(\overline{m}_3^2-1)
    \end{pmatrix}
\end{align*}
with $\overline{m}'=\partial_\xi \overline{m}$ and where we set $\overline{\boldsymbol{\beta}}\,\coloneqq \boldsymbol{\beta}(\overline{m}_3)$ as well as $\widetilde{\boldsymbol{\beta}}\coloneqq\betac_{\rm cp}/(1+c_{\rm cp} \overline{m}_3)^2$.

The spectrum of $\mathcal{L}(\overline{m})$ for DWs, e.g.\ in $L^2$ (tangent to the sphere at $\overline{m}$), decomposes in terms of Fredholm properties into \emph{essential} and \emph{point} spectrum, denoted by $\Sigma_{\textnormal{ess}}(\mathcal{L})$ and $\Sigma_{\textnormal{pt}}(\mathcal{L})$, respectively. Moreover, the boundary of $\Sigma_{\textnormal{ess}}(\mathcal{L})$ is contained in the (purely essential) spectra of $\pm\boldsymbol{e}_3$, cf.\ e.g.~\cite{Melcher2017}.

\subsection{Stability of Steady States}\label{sec:StabilitySteadyStates}
To start our analysis, we first review the stability of the steady states $\pm\boldsymbol{e}_3$ as homogeneous solutions to~\eqref{eq:LLGS} and its dependence on our main parameters $h$ and $c_{\rm cp}$. It is well-known that the constant up-magnetization state $+\boldsymbol{e}_3$ loses stability for an applied field below a threshold, whereas the constant down-magnetization state $-\boldsymbol{e}_3$ loses it above a certain threshold, both via a Hopf instability of the (purely) essential spectrum with frequency $\boldsymbol{\beta}^\pm/\alpha$, cf.\ e.g.~\cite{Melcher2017}. Moreover, due to the presence of a spin torque term, there exists an additional parameter regime, mainly characterized by the anisotropy, in which both states are (linearly) unstable (cf.\ Figure~\ref{fig:StabilityRegionsUp-AndDown-Magnetization}), cf.\ e.g.~\cite{Melcher2017, siemer2020}. More precisely, allowing only for perturbations tangential to the sphere at the poles, i.e., $\boldsymbol{n}=(n_1,n_2,0)\in \textnormal{T}_{\boldsymbol{m}} \mathbb{S}^2$, we obtain from~\eqref{eq:LinearEigenvalueProblem} with $s=\Omega=0$ the eigenvalue problems $\widetilde{\mathcal{L}}^+ \boldsymbol{n} = \lambda \boldsymbol{n}$ associated to $+\boldsymbol{e}_3$ as well as $\widetilde{\mathcal{L}}^- \boldsymbol{n} = \lambda \boldsymbol{n}$ associated to $-\boldsymbol{e}_3$, where
\[\widetilde{\mathcal{L}}^+\coloneqq\begin{pmatrix}
\frac{\alpha}{1+\alpha^2}\partial_\xi^2+\frac{\boldsymbol{\beta}^+-\alpha(h-\mu)}{1+\alpha^2} & \frac{1}{1+\alpha^2}\partial_\xi^2 - \frac{\alpha\boldsymbol{\beta}^++h-\mu}{1+\alpha^2} \\ -\frac{1}{1+\alpha^2}\partial_\xi^2+ \frac{\alpha\boldsymbol{\beta}^++h-\mu}{1+\alpha^2} & \frac{\alpha}{1+\alpha^2}\partial_\xi^2 +  \frac{\boldsymbol{\beta}^+-\alpha(h-\mu)}{1+\alpha^2}
\end{pmatrix}\]
and
\[\widetilde{\mathcal{L}}^-\coloneqq\begin{pmatrix}
\frac{\alpha}{1+\alpha^2}\partial_\xi^2+\frac{\alpha(h+\mu)-\boldsymbol{\beta}^-}{1+\alpha^2} & -\frac{1}{1+\alpha^2}\partial_\xi^2 - \frac{\alpha\boldsymbol{\beta}^-+h+\mu}{1+\alpha^2} \\ \frac{1}{1+\alpha^2}\partial_\xi^2+ \frac{\alpha\boldsymbol{\beta}^-+h+\mu}{1+\alpha^2} & \frac{\alpha}{1+\alpha^2}\partial_\xi^2 +  \frac{\alpha(h+\mu)-\boldsymbol{\beta}^-}{1+\alpha^2}
\end{pmatrix}.\]
The spectra of the (closed) linear matrix operators $\widetilde{\mathcal{L}}^\pm : H^2(\mathbb{R})\times H^2(\mathbb{R})\subset L^2(\mathbb{R})\times L^2(\mathbb{R}) \rightarrow L^2(\mathbb{R})\times L^2(\mathbb{R})$ determine the stability of the states $\pm\boldsymbol{e}_3$ for these function spaces. The associated dispersion relations read
\begin{equation}\label{eq:DispersionRelationHomogeneousStates}
\begin{split}
\widetilde{d}^\pm(\lambda, \mathbf{i} k)&=\left(-\frac{\alpha}{1+\alpha^2}k^2+\frac{\boldsymbol{\beta}^\pm\tau-\alpha(h-\mu\tau)\tau}{1+\alpha^2}-\lambda\right)^2\\
&\quad+\left(\frac{\tau}{1+\alpha^2}k^2-\frac{\alpha\boldsymbol{\beta}^\pm +h-\mu\tau}{1+\alpha^2} \right)^2,
\end{split}
\end{equation}
where $\tau=1$ for $+\boldsymbol{e}_3$ and $\tau=-1$ for $-\boldsymbol{e}_3$. Since each of $\widetilde{\mathcal{L}}^\pm$ has constant coefficients, $\Sigma_{\textnormal{pt}}(\widetilde{\mathcal{L}}^\pm)=\emptyset$, and the essential spectra parameterized by $k$ read
\begin{equation}\label{eq:EssentialSpectrumUp}
\Sigma_{\textnormal{ess}}(\widetilde{\mathcal{L}}^+)=\left\lbrace \frac{\boldsymbol{\beta}^+ -\alpha(k^2+ h-\mu)}{1+\alpha^2}\pm\mathbf{i}\frac{\alpha\boldsymbol{\beta}^+ +k^2+h-\mu}{1+\alpha^2}, k\in\mathbb{R} \right\rbrace
\end{equation}
as well as
\begin{equation}\label{eq:EssentialSpectrumDown}
\Sigma_{\textnormal{ess}}(\widetilde{\mathcal{L}}^-)=\left\lbrace \frac{\alpha( h+\mu-k^2)-\boldsymbol{\beta}^-}{1+\alpha^2}\pm\mathbf{i}\frac{\alpha\boldsymbol{\beta}^- +h+\mu-k^2}{1+\alpha^2}, k\in\mathbb{R} \right\rbrace.
\end{equation}
In particular, $k=0$ gives the most unstable points of $\Sigma_{\textnormal{ess}}(\widetilde{\mathcal{L}}^\pm)$, so that the homogeneous state $\boldsymbol{m}=+\boldsymbol{e}_3$ is (strictly) stable in case $h>\boldsymbol{\beta}^+/\alpha+\mu$ and $\boldsymbol{m}=-\boldsymbol{e}_3$ in case $h<\boldsymbol{\beta}^-/\alpha-\mu$. Note that the stability regions in the $h$-$c_{\rm cp}$-plane crucially depend on the sign of the anisotropy parameter. Although we are only interested in $\mu<0$, for clarity we include $\mu>0$, and hence define the boundaries $\Gamma^+\coloneqq(\boldsymbol{\beta}^+/\alpha)/(h-\mu)-1$ and $\Gamma^-\coloneqq 1-(\boldsymbol{\beta}^-/\alpha)/(h+\mu)$. For each fixed $\mu$, this results in the stability regions in the $h$-$c_{\rm cp}$-plane labelled in terms of the stability of $\pm\boldsymbol{e}_3$ as the aforementioned \emph{bistable} and \emph{monostable}, as well as \emph{biunstable}, where both are unstable. The monostable region may be further divided into two regions: $+\boldsymbol{e}_3$ stable and $-\boldsymbol{e}_3$ unstable, or vice versa. All regions are separated by the curves $\Gamma^\pm$ (cf.~Figure~\ref{fig:StabilityRegionsUp-AndDown-Magnetization} and see~\cite{siemer2020} for more details).

\begin{figure}[h]
    \centering
    \begin{subfigure}[b]{0.4\textwidth}
    \includegraphics[trim= 80 50 50 80, clip, width=\textwidth]{figures/steady_states_negative_mu_2.pdf}
    \caption*{\small $\mu<0$}
    \end{subfigure}
    \hspace*{15mm}
    \begin{subfigure}[b]{0.4\textwidth}
    \includegraphics[trim= 80 50 50 80, clip, width=\textwidth]{figures/steady_states_positive_mu_2.pdf}
    \caption*{\small $\mu>0$}
    \end{subfigure}
    \caption{Stability diagram of homogeneous states $\pm\boldsymbol{e}_3$ in $h$ and $c_{\rm cp}$ for $\alpha=1, \beta=3/4$, and $\mu=-1$ (left) as well as $\mu=1$ (right). State $+\boldsymbol{e}_3$ unstable to the left and stable to the right of $\Gamma^+$, $-\boldsymbol{e}_3$ stable to the left and unstable to the right of $\Gamma^-$.}
    \label{fig:StabilityRegionsUp-AndDown-Magnetization}
\end{figure}

\begin{Remark}
The terms \emph{bistable} and \emph{monostable} have been coined in the study of scalar parabolic PDEs, for which the theory of fronts, their stability, invasion and spreading has been originally developed. In the bistable case two homogeneous states are separated by an unstable state in the corresponding ODE for spatially constant states, and prototypical examples are the Nagumo equation in neuroscience~\cite{nagumo1962active} and the Allen-Cahn equation in materials science~\cite{allen1979microscopic}. The generic and extensively studied monostable paradigm for scalar parabolic PDEs is the Fisher-Kolmogorov-Petrovskii-Piscounov (FKPP) equation~\cite{aronson1975nonlinear, faye2019asymptotic, fisher1937wave, kolmogorov1937bulletin, tikhomirov1991study}.
\end{Remark}

In order to capture the monostable parameter regime in case $\mu<0$, we define:
\begin{equation}\label{eq:MonostableSetPlus}
\mathcal{M}^+\coloneqq\lbrace h\in\mathbb{R},\, c_{\rm cp}\in(-1,1): \boldsymbol{\beta}^-/\alpha-\mu<h \textnormal{ and } \boldsymbol{\beta}^+/\alpha+\mu<h \rbrace
\end{equation}
\begin{equation}\label{eq:MonostableSetMinus}
\mathcal{M}^-\coloneqq\lbrace h\in\mathbb{R},\, c_{\rm cp}\in(-1,1): h<\boldsymbol{\beta}^-/\alpha-\mu \textnormal{ and } h<\boldsymbol{\beta}^+/\alpha+\mu \rbrace
\end{equation}
Parameters from $\mathcal{M}^+$ lie to the `right' of $\Gamma^\pm$ in the $(h,\,c_{\rm cp})$-plane, whereas parameters which lie to the `left' of $\Gamma^\pm$ correspond to $\mathcal{M}^-$.
\begin{Remark}\label{rem:AlphaDependence}
The bistable parameter region increases for $\alpha \to 0$, whereas the monostable region with stable $\boldsymbol{e}_3$ moves to the left and that with $-\boldsymbol{e}_3$ stable to the right, respectively. On the other hand, the bistable region shrinks for $\alpha \to +\infty$, whereas the monostable regions as well as the biunstable region increase. This is in accordance with the underlying model and the description of damped precession around the applied field aligned with the $\boldsymbol{e}_3$-axis.
\end{Remark}

\section{Domain wall motion: bistable case}\label{sec:Bistable}
In case both asymptotic states are stable, the existence of a DW \emph{and} the selection of the speed $s$ and frequency $\Omega$ is equivalent to the construction of a heteroclinic orbit in the underlying coherent structure ODE~\eqref{eq:SphericalCoherentStructureODE}. It was shown in~\cite{siemer2020} that such a heteroclinic connection is generically `codimension-2' in the sense that both $s$ and $\Omega$ need to be selected in order to close the connection. It was also proven there that these DWs can be continued to inhomogeneous DWs in the parameter $c_{\rm cp}$ and it was further shown numerically that this holds for the entire interval $c_{\rm cp}\in(-1,1)$, see also Figure~\ref{fig:DWsSphere} (b)-(c). As mentioned above, this means that -- in contrast to the monostable case -- in general the selection of the speed and frequency are determined by the full nonlinear problem for the `interior' profile of a DW and not by the asymptotic states alone.

To start, we recall that~\eqref{eq:LLGS} reduces to~\eqref{eq:LLG} in case $\beta=0$ which is equivalent to $c_{\rm cp}=0$ with the additional shift in the constant applied field $h\to h-\beta/\alpha$~\cite{Melcher2017, siemer2020}. Let us now denote the shifted applied field by $\tilde{h}$ and set $c_{\rm cp}=0$. Concerning the explicit family of homogeneous DWs~\eqref{eq:AnalyticHomogeneousSolution}, denoted by $\boldsymbol{m}_0$, an upper bound on the applied field for linear stability was obtained in~\cite{gou2011stability}, where it was already mentioned that this is not optimal. To simplify notations, we rewrite the free energy~\eqref{eq:FreeEnergy} in the spherical coordinates~\eqref{eq:SphericalCoordinates} and without applied field as
\begin{equation}\label{eq:FreeEnergySpherical}
\boldsymbol{E}(\theta,\phi)=\frac{1}{2}\int \left(\left(\partial_x\theta\right)^2+\left(\partial_x\phi\right)^2\sin^2(\theta)+\mu\cos^2(\theta)\right)\, \textnormal{d}x.
\end{equation}

The idea of the proof presented in~\cite{gou2011stability} is not study the $\mathcal{O}(\varepsilon)$ correction, i.e., the spectrum of the linear operator arising from linearizing around $\boldsymbol{m}_0$, but rather show that to leading order the difference in the energies of the perturbed profile $\boldsymbol{m}_\varepsilon=\boldsymbol{m}_0 + \varepsilon \boldsymbol{m}_1 + \mathcal{O}(\varepsilon^2)$ and the explicit DW $\boldsymbol{m}_0$, denoted by $\Delta \boldsymbol{E}_\varepsilon$, decays to zero asymptotically in time.

One difference in the setups is that~\cite{gou2011stability} considers a different direction of the easy-axis anisotropy as well as applied field, which (for a nanowire geometry) requires an anisotropy parameter larger than zero that is normalized to one. However, our choice of applied field and anisotropy axis, and the requirement $\mu<0$ for a nanowire geometry, does not affect the computations compared with~\cite{gou2011stability}, since the temporal evolution of the energy $\!\,\mathrm{d}/\!\,\mathrm{d} t \, \boldsymbol{E}(\boldsymbol{m})$ and thus also $\!\,\mathrm{d}/\!\,\mathrm{d} t \, \Delta \boldsymbol{E}_\varepsilon$ does not change. We obtain the following result.

\begin{lemma}\label{lem:LinearStabilityArctan}
For any $\alpha>0,\,\beta\ge 0$, and $\mu<0$, the homogeneous DW given by~\eqref{eq:AnalyticHomogeneousSolution} for $c_{\rm cp}=0$ is linearly stable for $|h|<\beta/\alpha-9\mu/(9+2\sqrt{3})$. Moreover, this also holds for its unique continuation along $c_{\rm cp}$ for any sufficiently small $|c_{\rm cp}|$.
\end{lemma}

\begin{proof}We linearize around the explicit DW given by~\eqref{eq:AnalyticHomogeneousSolution} for $\sigma=-1$ only, since the result also holds for $\sigma=+1$ due to symmetry. We will obtain the bound on the applied field by improving the estimate in~\cite[Equation (38)]{gou2011stability}, therefore we use the notation as in~\cite{gou2011stability}, and thus set $\mu=-1$ for the moment.
Again, we denote by $\theta_0$ the function of the polar angle of $\boldsymbol{m}_0$ in spherical coordinates, see~\eqref{eq:AnalyticHomogeneousSolution}. Our simple observation is that $|\sin^2(\theta)\cos(\theta)|=|\cos(\theta)-\cos^3(\theta)|\le 2/(3\sqrt{3})$
for any $\theta$. This means the estimates in~\cite[Equation (38)]{gou2011stability} can be improved to
\[|F|\le \frac{2}{\varepsilon^2}\left(1+\frac{2}{3\sqrt{3}}\right)\Delta \boldsymbol{E}_\varepsilon,\]
where
\[F\coloneqq\int \cos(\theta_0) f_0 + \cos(\theta_0)\sin^2(\theta_0)\theta_1^2 \, \textnormal{d}x,\]
$f_0\coloneqq (\theta_1')^2 + \cos(2\theta_0)\theta_1^2 + \sin^2(\theta_0)(\phi_1')^2$, and where $\phi_1$ as well as $\theta_1$ are the angles corresponding to $\boldsymbol{m}_1$. Therefore, we obtain analogous to~\cite[Equation (39)]{gou2011stability} the slightly improved estimate
\begin{equation}\label{eq:EnergyEstimate}
\frac{\,\mathrm{d}}{\,\mathrm{d} t}\Delta \boldsymbol{E}_\varepsilon \le -2\alpha\left(1-\left(1+\frac{2}{3\sqrt{3}}\right)|\tilde{h}| \right)\Delta \boldsymbol{E}_\varepsilon .
\end{equation}
Although the estimate on $\sin^2(\theta)\cos(\theta)$ is indeed sharp, the one in~\eqref{eq:EnergyEstimate} is certainly not, as this includes several (non-sharp) estimates (see~\cite{gou2011stability} for detailed computations).

For general $\mu<0$, we obtain from~\eqref{eq:EnergyEstimate} that linear stability holds now for $|\tilde{h}|< -9\mu /(9+2\sqrt{3})$ and by applying the backward shift in the applied field, it holds for $|h|<\beta/\alpha-9\mu /(9+2\sqrt{3})$.

Therefore, in this regime, the linearization in $\boldsymbol{m}_0$ admits a spectral gap to the simple zero eigenvalues from translation as well as rotation. It directly follows from compact perturbation theory, see e.g.~\cite{kapitula2013spectral}, that the continuation of $\boldsymbol{m}_0$ to DWs for $0<|c_{\rm cp}|\ll 1$ established in~\cite{siemer2020} is also spectrally stable in this sense, as long as $|c_{\rm cp}|$ is sufficiently small.
\end{proof}
In direct comparison, the stability result presented in~\cite{gou2011stability} holds for an applied field less than $0.5$ with $\mu=-1$, while the improved bound is $9/(9+2\sqrt{3})\approx 0.722$. We emphasize that~\cite{gou2011stability} have studied nonlinear stability of $\boldsymbol{m}_0$ numerically via time-evolution of $\Delta \boldsymbol{E}_\varepsilon$ and have obtained that $\boldsymbol{m}_0$ is stable for an applied field $|h|<\beta/\alpha-\mu$ in our notation. This threshold is exactly the transition point to the monostable parameter regime $\mathcal{M}^\pm$ discussed in \S~\ref{sec:StabilitySteadyStates}, where one of the two stable asymptotic states $\pm\boldsymbol{e}_3$ becomes unstable.

\begin{Remark}\label{rem:HomogeneousDW}
In~\cite{gou2011stability}, also stability results for the case of a small transverse applied field and also for a small hard-axis anisotropy have been obtained, based on their stability bound. Therefore, one could potentially improve these results utilizing Lemma~\ref{lem:LinearStabilityArctan}.
\end{Remark}

\subsection{Standing Domain Walls}\label{sec:StandingDomainWalls}
We turn to the asymptotic propagation behavior in the bistable parameter regime. In particular, we introduce a heuristic approach for standing DWs to describe the higher order terms in~\eqref{eq:PropagationSignIntegral}. This approach utilizes~\eqref{eq:Potential} and although the result is not exact, it provides a surprisingly good approximation even for large $|c_{\rm cp}|<1$.

Bistability for a free energy corresponds to (at least) two local minima and one expects front propagation towards the state with lower energy. In analogy, we consider the free energy~\eqref{eq:FreeEnergy} in case $c_{\rm cp}=0$ and apply again the shift $h\mapsto h-\beta/\alpha$. Upon neglecting the Dirichlet energy term, the shifted free energy depends on $m_3$ only and possesses the potential
\[V(m_3)\coloneqq-\left(h-\beta/\alpha\right)m_3+\frac{1}{2}\mu m_3^2.\]
From this, we compute that $V(+1)>V(-1)$ in case $h<\beta/\alpha$ and $V(+1)<V(-1)$ if $\beta/\alpha<h$. Therefore, any DW connecting $+\boldsymbol{e}_3$ $[-\boldsymbol{e}_3]$ on the left with $-\boldsymbol{e}_3$ $[+\boldsymbol{e}_3]$ on the right will invade the state $-\boldsymbol{e}_3$ $[+\boldsymbol{e}_3]$ for $\beta/\alpha<h<\beta/\alpha-\mu$ $[\textnormal{for }\beta/\alpha+\mu<h<\beta/\alpha]$ and will invade the state $+\boldsymbol{e}_3$ $[-\boldsymbol{e}_3]$ for $\beta/\alpha+\mu<h<\beta/\alpha$ $[\textnormal{for }\beta/\alpha<h<\beta/\alpha-\mu]$. Note that these directions of propagation are in complete agreement with~\eqref{eq:velsym} and the definition of $\sigma$ as well as the sign of $s$ in~\eqref{eq:HomogeneousSpeed} of the (explicit) homogeneous DW family~\eqref{eq:AnalyticHomogeneousSolution}.

In order to extend the expression of $V$ to cover also the case $c_{\rm cp}\neq0$ and thus approximate higher order terms in~\eqref{eq:PropagationSignIntegral}, we include~\eqref{eq:Potential}, set $\beta=\boldsymbol{\beta}$ from \eqref{e:defbeta}, and obtain
\[\widetilde{V}(m_3)\coloneqq -hm_3 + \frac{1}{2}\mu m_3^2 + \frac{\beta}{\alpha c_{\rm cp}}\ln\left(1+c_{\rm cp} m_3 \right),\]
which coincides with $V$ in the limit $c_{\rm cp}\to0$. Concerning the energy difference and direction of motion, we define the quantity
\[\mathcal{H}\coloneqq\frac{\beta}{2\alphac_{\rm cp}}\left(\ln(1+c_{\rm cp})-\ln(1-c_{\rm cp})\right)\]
and observe that
$\widetilde{V}(+1)>\widetilde{V}(-1)$ holds for $h<\mathcal{H}$ and $\widetilde{V}(+1)<\widetilde{V}(-1)$ if $h>\mathcal{H}$. From this, we presume that a DW connecting $+\boldsymbol{e}_3$ on the left and $-\boldsymbol{e}_3$ on the right will invade $+\boldsymbol{e}_3$ if $\boldsymbol{\beta}^+/\alpha+\mu<h<\mathcal{H}$ and $-\boldsymbol{e}_3$ if $\mathcal{H}<h<\boldsymbol{\beta}^-/\alpha-\mu$. On the other hand, DWs connecting $-\boldsymbol{e}_3$ on the left and $+\boldsymbol{e}_3$ on the right should propagate towards $+\boldsymbol{e}_3$ for $\mathcal{H}<h<\boldsymbol{\beta}^-/\alpha-\mu$ and towards $-\boldsymbol{e}_3$ if $\boldsymbol{\beta}^+/\alpha+\mu<h<\mathcal{H}$. This is indeed in good agreement with numerical simulations for $h$ and $c_{\rm cp}$ above $\Gamma^-$ as well as above $\Gamma^+$, as illustrated in Figure~\ref{fig:SpeedSignBistableCase} (a). Here we have simulated the PDE dynamics with \textsc{Pde2Path}\xspace from step function-like initial conditions connecting the state $+\boldsymbol{e}_3$ on the left and $-\boldsymbol{e}_3$ on the right, and have tracked the asymptotic propagation speed over time by the freezing method, which is described in more detail in \S~\ref{sec:NumericalSimulationMethod}. We corroborate these predictions for $s=0$ by numerical continuation in \textsc{AUTO-07P} with $s=0$ fixed in order to obtain the loci of (numerically) exact stationarity. We illustrate the results in Figure~\ref{fig:SpeedSignBistableCase} (b) for the range $-0.598<c_{\rm cp}\le 0.9$, where the lower bound represents the point of transition into the $\mathcal{M}^-$ regime. We further remark that $\mathcal{H}$ does not give precise predictions: it is symmetric in the parameter $c_{\rm cp}$, whereas system~\eqref{eq:SphericalCoherentStructureODE} is not, and indeed the numerical results in Figure~\ref{fig:SpeedSignBistableCase} (b) obtained from continuation (red dashed) are slightly asymmetric. Nevertheless, we conclude that $\mathcal{H}-h$ approximates the sign of the asymptotic spreading speed surprisingly well for a broad range of $c_{\rm cp}\neq 0$.

\begin{figure}
    \centering
    \begin{subfigure}[b]{0.56\textwidth}
    	\includegraphics[trim=50 120 70 130,clip,width=\textwidth]{figures/speed_freezed2.pdf}
    	\caption*{(a)}
    \end{subfigure}
    \hspace*{10mm}
    \begin{subfigure}[b]{0.3\textwidth}
    	\includegraphics[trim=10 0 20 35,clip,width=\textwidth]{figures/zerospeedcontinuation2.eps}
    	\caption*{(b)}
    \end{subfigure}
    \caption{Comparison of $h=\mathcal{H}$
     with numerical computations in part of the bistable region (above $\Gamma^+$ and $\Gamma^-$) in left panel of Figure~\ref{fig:StabilityRegionsUp-AndDown-Magnetization}; the other parameters are taken as there. (a) $h\in[\beta/\alpha+\mu, \beta/\alpha-\mu]$ and $c_{\rm cp}\in[-0.9,0.9]$. Green curve marks $h=\mathcal{H}$; left of it $h>\mathcal{H}$,
    right of it $h<\mathcal{H}$.
    Black shaded region marks numerically observed approach of DWs with negative speed. In gray region the speed is within $\pm$\num{e-6}. In white bistable region DWs attain positive speed. Initial data is step function-type connecting $+\boldsymbol{e}_3$ on the left and $-\boldsymbol{e}_3$ on the right. (b) Magnification around green curve in (a) and red dashed line results from numerical continuation for $-0.598<c_{\rm cp}\le 0.9$ with $s=0$ fixed.}
    \label{fig:SpeedSignBistableCase}
\end{figure}

\section{Domain wall motion: monostable case}\label{sec:Monostable}
The study of interfaces invading an unstable state, also in the context of ferromagnetic dynamics, has a long history, see e.g.~\cite{aronson1975nonlinear, ebert2000front, elmer1993front, elmer1994dual, fisher1937wave, garnier2012inside, goussev2020dynamics, holzer2014anomalous, kolmogorov1937bulletin, van2003front}. The most accurate method for predicting the asymptotic behavior in systems of parabolic equations is based on the so-called pointwise Green's function associated with the unstable state~\cite{beck2017stability, holzer2014anomalous, holzer2014criteria}, rather than the saddle-point analysis that is sufficient for the classical scalar case~\cite{ebert2000front, van2003front}. It turns out that in the present case this can be studied via the so-called absolute spectrum as defined below.

We will discuss a more general problem and present an explicit expression for the absolute spectrum, denoted by $\Sigma_{\textnormal{abs}}(\mathscr{L})$, of certain complex linear $2\times 2$ matrix operators $\mathscr{L}$. This class also includes the operator arising from linearizing the complex Ginzburg-Landau equation around the origin, as discussed below. Based on the explicit formula for $\Sigma_{\textnormal{abs}}(\mathscr{L})$, we will show that its most unstable elements are \emph{simple pinched double roots} and correspond to singularities of the pointwise Green's function. From this, we will infer the so-called \emph{linear spreading speed} and also define the associated \emph{linear spreading frequency}, where both provide predictions for the asymptotic behavior of DWs in the nonlinear system~\eqref{eq:LLGS}.

Although the linear spreading speed can be obtained from values of $s$ for which the critical pinched double root is placed on the imaginary axis, the spreading frequency is not selected in the sense of (linear) stability properties. It is selected \emph{passively} from the oscillatory dynamics, since -- due to rotation symmetry -- changes in frequency modify the imaginary part only.

To simplify notations, we will focus on the parameter regime $\mathcal{M}^+$ introduced in~\eqref{eq:MonostableSetPlus}; all results presented in this section also transfer to the regime $\mathcal{M}^-$.

\subsection{Pointwise growth, absolute spectra, and double roots}\label{sec:Pointwise growth, absolute spectra, and double roots}

The pointwise Green's function in our situation can be defined as follows. Again, we allow for perturbations tangential to the sphere at the pole $-\boldsymbol{e}_3$ only and obtain from linearization of~\eqref{eq:QuasilinearSpeedFrequencyLLGS}
\[\partial_t \boldsymbol{n}=\mathcal{L}^- \boldsymbol{n} + \mathcal{N}(\boldsymbol{n}),\]
where the linear matrix operator $\mathcal{L}^-$ is given via~\eqref{eq:LinearEigenvalueProblem} for $\overline{m}=-\boldsymbol{e}_3$ and $\mathcal{N}(\boldsymbol{n})$ describes the nonlinear terms. Considering the linear part only and after a Laplace transform in time, the equation reduces to a system of second order ODEs which depends on the spectral parameter $\lambda$. The associated pointwise (matrix) Green's function $\boldsymbol{G}_\lambda(\xi-y)$ is defined as the solution to
\[(\mathcal{L}^--\lambda)\boldsymbol{G}_\lambda(\xi-y)=\delta(y)\]
in the distributional sense. Applying the inverse Laplace transform with a contour $\mathcal{C}$ to the right of the essential spectrum of $\mathcal{L}^-$, which we determine later, the solution to the initial problem can now be written as
\begin{equation}\label{eq:LinearSolution}
\boldsymbol{n}=-\frac{1}{2\pi\mathbf{i}}\int_\mathcal{C}\mathrm{e}^{\lambda t}\int_\mathbb{R} \boldsymbol{G}_\lambda(\xi-y)\boldsymbol{n}(y,0)\,\textnormal{d}y\, \textnormal{d}\lambda,
\end{equation}
where we supposed that $\boldsymbol{n}(y,0)$ has compact support. We notice that the contour $\mathcal{C}$ can be deformed beyond the essential spectrum up to singularities of $\boldsymbol{G}_\lambda$~\cite{gardner1998gap, kapitula2013spectral}. In this way, if $\boldsymbol{G}_\lambda$ can be analytically continued into the left half-plane, decay bounds on the solution can be inferred.

Analytical properties of $\boldsymbol{G}_\lambda$ depend on the frame of reference and thus on $s$ and $\Omega$ in our situation. It turns out that for $s=0$, $\boldsymbol{G}_\lambda$ possesses singularities at the complex conjugated pair of points of the essential spectrum with maximal real part. For more general operators it is well known that upon changing $s$, the singularities can be moved into the stable half-plane, which describes precisely the transition from pointwise growth to pointwise decay. As mentioned before, the frequency $\Omega$ has no influence on the real part of the singularities. Viewing $s$ now as a parameter, the transition from pointwise growth to pointwise decay determines the change from so-called \emph{convective} to \emph{absolute} instability of $-\boldsymbol{e}_3$. Here, we call an instability convective, if perturbations grow in norm, but decay pointwise, while for an absolute instability a generic perturbation grows exponentially at each point where it was applied~\cite{kapitula2013spectral, sandstede2002stability, sandstede2000absolute}.

As mentioned, relevant singularities of $\boldsymbol{G}_\lambda$ in our situation can be identified via the absolute spectrum associated to $\mathcal{L}^-$. Here, we consider the more general (complex) linear operators
\begin{equation}\label{eq:GeneralMatrixOperator}
\mathscr{L}\coloneqq \begin{pmatrix} a_2\partial_\xi^2+a_1\partial_\xi + a_0 & -b_2\partial_\xi^2 -b_1\partial_\xi- b_0\\ b_2\partial_\xi^2 +b_1\partial_\xi  +b_0 & a_2\partial_\xi^2+a_1\partial_\xi + a_0\end{pmatrix},
\end{equation}
with $a_j, b_j\in \mathbb{C}, j=0,1,2$, $a_2, b_2\neq 0$, and whose dispersion relation reads
\begin{equation}\label{eq:GeneralDispersionRelation}
\mathscr{D}(\lambda,\nu)\coloneqq \left(a_2\nu^2+a_1\nu+a_0-\lambda\right)^2+\left( b_2\nu^2+b_1\nu+b_0\right)^2,
\end{equation}
For what follows it is convenient to denote
\[a_2^r\coloneqq\mathrm{Re}(a_2),\, b_2^r\coloneqq\mathrm{Re}(b_2),\, a_2^i\coloneqq\mathrm{Im}(a_2),\, b_2^i\coloneqq\mathrm{Im}(b_2).\]
We shall assume the following, which relates to uniform ellipticity of $\mathscr{L}$.
\begin{Hypothesis}\label{hyp:Coefficients} The coefficients $a_2, b_2\in\mathbb{C}$ satisfy $a^r_2>|b_2^i|^2$.
\end{Hypothesis}

In order to define the absolute spectrum, cf.~\cite{sandstede2000absolute}, we solve $\mathscr{D}(\lambda,\nu)=0$ for given $\lambda$, which leads to four so-called spatial eigenvalues $\nu_j(\lambda),\, j=1,\dots,4$, including multiplicity, that we order by descending real parts as
\begin{equation}\label{eq:Ordering}
\mathrm{Re}(\nu_1(\lambda))\ge \mathrm{Re}(\nu_2(\lambda))\ge\mathrm{Re}(\nu_3(\lambda))\ge \mathrm{Re}(\nu_4(\lambda)).
\end{equation}
Next, we determine the Morse index~\cite{sandstede2000absolute} associated to the spatial eigenvalues, i.e., the number of unstable ones for sufficiently large $\mathrm{Re}(\lambda)$.
\begin{lemma}\label{lem:MorseIndexGeneral}
For $a_2, b_2$ satisfying Hypothesis~\ref{hyp:Coefficients}, the Morse index associated to $\mathscr{L}$ is $i_\infty=2$.
\end{lemma}
\begin{proof}
We consider $\mathrm{Re}(\lambda)\gg 1$ by taking $0<\varepsilon\ll1$ in the parabolic scaling $\lambda=\tilde{\lambda}/\varepsilon^2$ as well as $\nu=\tilde{\nu}/\varepsilon$ for arbitrary fixed $\tilde{\lambda}\in \mathbb{C}$ with $\mathrm{Re}(\tilde{\lambda})>0$. Multiplication of $\mathscr{D}(\tilde{\lambda},\tilde{\nu})=0$ by $\varepsilon^4$ leads to
\[(a^2_2+b^2_2)\tilde{\nu}^4-2a_2\tilde{\lambda}\tilde{\nu}^2+\tilde{\lambda}^2+\mathcal{O}(\varepsilon)=0,\]
whose solutions for $\varepsilon=0$ with four choices of signs are
\[\tilde{\nu}_\pm^\pm(\tilde{\lambda})=\pm\sqrt{\frac{\tilde{\lambda}}{a_2\pm\mathbf{i} b_2}}.\]
Under the assumptions on $a_2$ and $b_2$, we show that $\tilde{\nu}_\pm^\pm(\tilde{\lambda})\notin\mathbf{i}\mathbb{R}$, which implies the claim, i.e., $\mathrm{Re}(\nu_{1})=\mathrm{Re}(\nu_{2})>0> \mathrm{Re}(\nu_{3})=\mathrm{Re}(\nu_{4})$ in terms of the notation \eqref{eq:Ordering}.
Since the sign of $b_2$ will be irrelevant, it suffices to consider the case $+ \mathbf{i} b_2$ and assume that $\tilde{\lambda}/(a_2+ \mathbf{i} b_2)\in \mathbb{R}^-$ under Hypothesis~\ref{hyp:Coefficients}. From this, $\mathrm{Im}(\tilde{\lambda}/(a_2+ \mathbf{i} b_2))=0$ is equivalent to $(a_2^r-b_2^i)\tilde{\lambda}^i=(a_2^i+b_2^r)\tilde{\lambda}^r$, where $\tilde{\lambda}^r\coloneqq \mathrm{Re}(\tilde{\lambda})$ and $\tilde{\lambda}^i\coloneqq \mathrm{Im}(\tilde{\lambda})$. In case $a_2^r=b_2^i$, it follows that $a_2^i=-b_2^r$, which together leads to $a_2+\mathbf{i} b_2=0$, contradicting Hypothesis~\ref{hyp:Coefficients}. Thus $a_2^r\neq b_2^i$ and we obtain
\[\tilde{\lambda}^i=\tilde{\lambda}^r \frac{a_2^i+b_2^r}{a_2^r-b_2^i}.\]
Now $\mathrm{Re}(\tilde{\lambda}/(a_2\pm \mathbf{i} b_2))<0$ with $\mathrm{Re}(\tilde{\lambda})>0$ implies
\[a_2^r-b_2^i + \frac{(a_2^i+b_2^r)^2}{a_2^r-b_2^i}<0.\]
This inequality only holds for $a_2^r<b_2^i$ and hence contradicts Hypothesis~\ref{hyp:Coefficients}.
\end{proof}

It follows that under Hypothesis \ref{hyp:Coefficients} the absolute spectrum of $\mathscr{L}$ with~\eqref{eq:Ordering} is given by the set
\[\Sigma_{\textnormal{abs}}(\mathscr{L})=\lbrace \lambda \in \mathbb{C} : \mathscr{D}(\lambda,\nu)=0 \textnormal{ and } \mathrm{Re}(\nu_2(\lambda))=\mathrm{Re}(\nu_3(\lambda))\rbrace.\]

We now turn to the aforementioned points in the absolute spectrum that can give rise to singularities of the pointwise Green's function~\cite{holzer2014criteria, sandstede2000absolute}.
\begin{definition}\label{def:DoubleRoots}
A pair $(\lambda^\textnormal{dr},\nu^\textnormal{dr})\in\mathbb{C}^2$ is called \emph{double root} if
\[\mathscr{D}(\lambda^\textnormal{dr},\nu^\textnormal{dr})=0\quad \textnormal{and} \quad \partial_\nu \mathscr{D}(\lambda^\textnormal{dr},\nu^\textnormal{dr})=0.\]
A double root is called \emph{simple} if $\partial_\lambda \mathscr{D}(\lambda^\textnormal{dr},\nu^\textnormal{dr})\neq 0$ and $\partial_\nu^2 \mathscr{D}(\lambda^\textnormal{dr},\nu^\textnormal{dr})\neq 0$.
\end{definition}

To obtain a direct relation between double roots and singularities of $\boldsymbol{G}_\lambda$, one has to further require a \emph{pinching condition}, where non-pinched double roots correspond to removable singularities and an analytic extension of $G_\lambda$ beyond these roots is available~\cite{holzer2014anomalous, holzer2014criteria, sandstede2000absolute}.
\begin{definition}\label{def:PinchingCondition}
A double root $(\lambda^\textnormal{dr},\nu^\textnormal{dr})$ is called \emph{pinched} if there exists a continuous curve $\Lambda(\tau)\in\mathbb{C}, \tau\ge 0$, with $\Lambda(0)=\lambda^\textnormal{dr}$, $\mathrm{Re}(\Lambda(\tau))\to +\infty$ for $\tau\to +\infty$, and continuous curves of roots $\nu_\pm(\Lambda(\tau))$ to $\mathscr{D}(\Lambda(\tau),\nu)=0$ with $\nu_\pm(\lambda^\textnormal{dr})=\nu^\textnormal{dr}$ and $\mathrm{Re}(\nu_\pm(\Lambda(\tau)))\to\pm\infty$ for $\tau\to +\infty$.
\end{definition}

Decisive for the identification of pinched double roots in our setting is the following observation. In preparation, denote by $\mathcal{G}_\mathrm{abs}\subset \mathbb{C}$ the connected component of $\mathbb{C} \setminus\Sigma_{\textnormal{abs}}$ that contains an unbounded interval of $\mathbb{R}^+$. Since the Morse index $i_\infty$ is well-defined, there exists a $\rho\in \mathbb{R}$ such that $\Sigma_{\textnormal{abs}}\cap \lbrace \lambda\in\mathbb{C} : \mathrm{Re}(\lambda)\ge \rho \rbrace=\emptyset$, so that $\mathcal{G}_\mathrm{abs}$ is well-defined.

\begin{lemma}\label{lem:PinchingAndMorseIndex}
A double root $(\lambda^\textnormal{dr},\nu^\textnormal{dr})$ for $\mathscr{D}$ with $\lambda^\textnormal{dr}\in \partial \mathcal{G}_\mathrm{abs}$ is pinched. In particular, this is the case if $\mathrm{Re}(\lambda^\textnormal{dr}) = \max\mathrm{Re} (\Sigma_{\textnormal{abs}})$.
\end{lemma}
\begin{proof}
Since $\lambda^\mathrm{dr}\in\partial \mathcal{G}_\mathrm{abs}$, one finds a continuous curve $\Lambda(\tau)\in\mathbb{C}, \tau\ge 0$, with $\Lambda(0)=\lambda^\mathrm{dr}$, into the region $\lbrace \lambda\in\mathbb{C}:\mathrm{Re}(\lambda)>\max \mathrm{Re}(\Sigma_{\textnormal{ess}})\rbrace\subset \mathbb{C}$, such that $\Sigma_{\textnormal{abs}}\cap\Lambda=\lbrace \lambda^\mathrm{dr}\rbrace$. We may assume there is no crossing with $\Sigma_{\textnormal{abs}}$ along $\Lambda$ for $\tau>0$ and thus the roots $\nu(\Lambda(\tau))$ of $\mathcal{D}$ are ordered as in~\eqref{eq:Ordering}. The fact that the real parts grow unboundedly with opposite signs follows from the scaling result in Lemma~\ref{lem:MorseIndexGeneral} (and similarly for general parabolic well-posed problems).
\end{proof}
\begin{Remark}
If all curves that connect $\lambda^\mathrm{dr}$ with a point in $\lbrace \lambda\in\mathbb{C}:\mathrm{Re}(\lambda)>\max \mathrm{Re}(\Sigma_{\textnormal{ess}})\rbrace\subset \mathbb{C}$ intersect $\Sigma_{\textnormal{abs}}$, then Lemma~\ref{lem:PinchingAndMorseIndex} is not applicable and it is not known under what (additional) conditions the double root would be pinched.
\end{Remark}

Due to the index conditions, explicit representations of $\Sigma_{\textnormal{abs}}$ for higher order or non-scalar cases, e.g.\ the Swift-Hohenberg equation, are rarely possible. An exception is the complex Ginzburg-Landau (CGL) equation and due to the similar structure in the dispersion relation, it is possible to obtain explicit expressions for the class of operators $\mathscr{L}$. Our approach is based on the fact that $\Sigma_{\textnormal{abs}}$ is the union of curves that must lie fully in the region bounded by the essential spectrum that includes an unbounded part of $\mathbb{R}^-$~\cite{kapitula2013spectral, rademacher2006geometric}. Moreover, each point of the absolute spectrum arises as an intersection point of (at least) two curves of the essential spectrum in exponentially weighted $L^2$-spaces~\cite{fiedler2003spatio, rademacher2006geometric}. Here the essential spectrum of $\mathscr{L}$ in such a space with exponential weight $\eta\in\mathbb{R}$ consists of solutions to the dispersion relation~\eqref{eq:GeneralDispersionRelation} with $\nu=\mathbf{i} k+\eta$, i.e., the set
\[\Sigma_{\textnormal{ess}}^\eta(\mathscr{L})=\lbrace \lambda\in\mathbb{C}:\mathscr{D}(\lambda,\mathbf{i} k+\eta)=0, k\in\mathbb{R} \rbrace;\]
note $\Sigma_{\textnormal{ess}}(\mathscr{L})=\Sigma_{\textnormal{ess}}^0(\mathscr{L})$. The key to locate pinched double roots in our case is the following.
\begin{theorem}\label{theo:AbsolutSpectrumGeneralLinearMatrixOperator}
For $a_2, b_2$ satisfying Hypothesis~\ref{hyp:Coefficients}, the absolute spectrum of $\mathscr{L}$ is given by the following two (straight) half-lines parameterized by $r\in[0,\infty)$,
\[-(a_2+\mathbf{i} b_2) r-\frac{(a_1+\mathbf{i} b_1)^2}{4(a_2+\mathbf{i} b_2)}+a_0+\mathbf{i} b_0\quad \textnormal{and}\quad -(a_2-\mathbf{i} b_2) r-\frac{(a_1-\mathbf{i} b_1)^2}{4(a_2-\mathbf{i} b_2)}+a_0-\mathbf{i} b_0.\]
The rightmost points lie at $r=0$ and are pinched double roots of~\eqref{eq:GeneralDispersionRelation}.
\end{theorem}
\begin{proof} To shorten notations, we set $c_j=a_j+\mathbf{i} b_j, \tilde{c}_j=a_j-\mathbf{i} b_j$ and
\[c_j^r=\mathrm{Re}(c_j), \,c_j^i=\mathrm{Im}(c_j), \,\tilde{c}_j^r=\mathrm{Re}(\tilde{c}_j), \,\tilde{c}_j^i=\mathrm{Im}(\tilde{c}_j)\]
for $j=0,1,2$. Starting from $\mathscr{D}(\lambda,\mathbf{i} k +\eta)=0$ with wavenumber $k\in\mathbb{R}$ and weight $\eta\in\mathbb{R}$, the weighted essential spectrum of $\mathscr{L}$ is given by the two (parabolic) curves
\[\lambda_1(k;\eta)=-c_2k^2 + \mathbf{i}(c_1+2c_2\eta)k+c_2\eta^2+c_1\eta+c_0\]
and
\[\lambda_2(k;\eta)=-\tilde{c}_2k^2 + \mathbf{i}(\tilde{c}_1+2\tilde{c}_2\eta)k+\tilde{c}_2\eta^2+\tilde{c}_1\eta+\tilde{c}_0\]
for $k\in\mathbb{R}$. The maximal real part of $\lambda_1$ satisfies $\mathrm{Re}\left(\partial_k \lambda_1(k;\eta)\right)=0$ which is solved by $k^\ast\coloneqq-\frac{c_1^i+2c_2^i\eta}{2c_2^r}$.
Additionally,
\[\mathrm{Re}\left(\partial_\eta \lambda_1(k^\ast; \eta)\right)=2\frac{(c_2^r)^2 + (c_2^i)^2}{c_2^r}\eta+c_1^r+\frac{c_2^i c_1^i}{c_2^r},\]
which vanishes at
\[\eta^\ast=-\frac{c_2^r c_1^r + c_2^i c_1^i}{2\left((c_2^r)^2+(c_2^i)^2\right)}.\]
For this value of $\eta$, we obtain
\begin{align*}\lambda_1(k;\eta^\ast)&=-c_2k^2+\mathbf{i}\frac{c_2^i c_1^r - c_2^rc_1^i}{c_2^i+\mathbf{i} c_2^r}k\\ &\quad + \frac{\left(c_2^rc_1^2 +  c_2^i c_1^i - 2\mathbf{i} c_2^ic_1^r + 2\mathbf{i} c_2^rc_1^i \right) \left( c_2^r c_1^r+ c_2^i c_1^i \right)}{4c_2 \left(c_2^i + \mathbf{i} c_2^r \right)^2} + c_0\\
&=-c_2\overline{k}^2 -\frac{c_1^2}{4c_2} + c_0,
\end{align*}
where we introduced the shift
\[\overline{k}\coloneqq k+\frac{c_2^r c_1^i-c_2^i c_1^r}{2(c_2^r)^2+2(c_2^i)^2}\in\mathbb{R}.\]
This leads to a straight half-line with slope $-c_2$ and maximal real part at
\[\frac{c_2^r(c_1^i)^2-c_2^r(c_1^r)^2-2c_2^ic_1^rc_1^i}{4((c_2^r)^2+(c_2^i)^2)}+c_0^r \quad \textnormal{for} \quad k=\frac{c_2^ic_1^r-c_2^rc_1^i}{2(c_2^r)^2 + 2(c_2^i)^2}.\]
Analogously, for $\lambda_2$ we obtain
\[\lambda_2(k;\eta^\ast)=-\tilde{c}_2 \check{k}^2-\frac{\tilde{c}_1^2}{4\tilde{c}_2}+\tilde{c}_0, \quad \check{k}\coloneqq k+\frac{\tilde{c}_2^r \tilde{c}_1^i-\tilde{c}_2^i \tilde{c}_1^r}{2(\tilde{c}_2^r)^2+2(\tilde{c}_2^i)^2}\in\mathbb{R}.\]
To locate double roots, we determine the roots of the resultant function of $\mathscr{D}(\lambda,\nu)$ and $\partial_\nu \mathscr{D}(\lambda,\nu)$ in $\lambda$, which indeed results in the two endpoints. Since $\Sigma_{\textnormal{abs}}(\mathscr{L})$ has to lie in the region bounded by $\Sigma_{\textnormal{ess}}^\eta(\mathscr{L})$ (see~\cite{rademacher2006geometric} for full statement) which forms (doubly covered) straight lines precisely at $\eta=\eta^*$, it follows that $\Sigma_{\textnormal{abs}}(\mathscr{L})$ is given by $\lambda_{1,2}(k;\eta^\ast)$. One readily checks that the double roots $\lambda^\pm_\textnormal{dr}$ are the rightmost points of $\lambda_{1,2}(k;\eta^\ast)$ and thus also of $\Sigma_{\textnormal{abs}}(\mathscr{L})$. Lemma~\ref{lem:PinchingAndMorseIndex} implies that these are pinched.
\end{proof}

Theorem~\ref{theo:AbsolutSpectrumGeneralLinearMatrixOperator} allows to explicitly represent the absolute spectrum of $\mathcal{L}^-$ and its most unstable modes. See Figure~\ref{fig:SpecsOfDownMagnetization} for illustrations. The dispersion relation of $\mathcal{L}^-$ reads
\begin{equation}\label{eq:DispersionRelationDownState}
\begin{split}
d^-(\lambda, \nu)&\coloneqq\left(\frac{\alpha}{1+\alpha^2}\nu^2+s\nu+\frac{\alpha(h+\mu)-\boldsymbol{\beta}^-}{1+\alpha^2}-\lambda\right)^2\\
&\quad+\left(\frac{1}{1+\alpha^2}\nu^2-\Omega+\frac{h+\mu+\alpha\boldsymbol{\beta}^-}{1+\alpha^2} \right)^2
\end{split}
\end{equation}
which gives the form~\eqref{eq:GeneralDispersionRelation} with
\begin{equation}\label{eq:GeneralCoefficientsA}
a_2=\frac{\alpha}{1+\alpha^2},\quad a_1=s,\quad a_0=\frac{\alpha(h+\mu)-\boldsymbol{\beta}^-}{1+\alpha^2},
\end{equation}
as well as
\begin{equation}\label{eq:GeneralCoefficientsB}
b_2=\frac{1}{1+\alpha^2},\quad b_1=0, \quad b_0=\frac{h+\mu+\alpha\boldsymbol{\beta}^-}{1+\alpha^2}-\Omega.
\end{equation}
In particular, Hypothesis~\ref{hyp:Coefficients} is always satisfied for~\eqref{eq:GeneralCoefficientsA} and~\eqref{eq:GeneralCoefficientsB}, and thus we obtain the following statement.

\begin{corollary}\label{cor:AbsoluteSpectrum}
The absolute spectrum $\Sigma_{\textnormal{abs}}(\mathcal{L}^-)$ of the linearization around the unstable state $-\boldsymbol{e}_3$ in the $\mathcal{M}^+$ regime is given by
\[\Sigma_{\textnormal{abs}}(\mathcal{L}^-)=\Sigma_{\textnormal{abs}}^+(\mathcal{L}^-) \cup \Sigma_{\textnormal{abs}}^-(\mathcal{L}^-),\]
where $\Sigma_{\textnormal{abs}}^\pm(\mathcal{L}^-)$ are the straight half-lines
\[\Big\lbrace -\frac{\alpha\pm\mathbf{i}}{1+\alpha^2}k-\frac{\alpha\mp\mathbf{i}}{4}s^2+\frac{\alpha(h+\mu)-\boldsymbol{\beta}^-}{1+\alpha^2}\pm\mathbf{i}\left( \frac{h+\mu+\alpha\boldsymbol{\beta}^-}{1+\alpha^2}-\Omega \right), k\in[0,\infty) \Big\rbrace\]
in the complex plane. Moreover, the most unstable points are pinched double roots that are simple if $s\neq 0$, and are given by
\begin{equation}\label{eq:DoubleRootLambda}
\lambda^\textnormal{dr}_\pm=-\alpha s^2/4\pm \mathbf{i} s^2/4+\frac{\alpha(h+\mu)-\boldsymbol{\beta}^-}{1+\alpha^2}\pm\mathbf{i}\left(\frac{h+\mu+\alpha\boldsymbol{\beta}^-}{1+\alpha^2}-\Omega\right),
\end{equation}
with $\nu^\textnormal{dr}_+=\nu(\lambda_\textnormal{dr}^+)=-\alpha s/2+\mathbf{i} s/2$ and $\nu^\textnormal{dr}_-=\nu(\lambda_\textnormal{dr}^-)=-\alpha s/2-\mathbf{i} s/2$.
\end{corollary}
\begin{proof}
The first part is a direct consequence of Theorem~\ref{theo:AbsolutSpectrumGeneralLinearMatrixOperator} via~\eqref{eq:GeneralCoefficientsA} and~\eqref{eq:GeneralCoefficientsB}, from which it also follows that the rightmost points are pinched double roots given by $(\lambda^\textnormal{dr}_\pm,\nu^\textnormal{dr}_\pm)$. Solving for $d_\lambda(\lambda^\textnormal{dr}_\pm,\nu^\textnormal{dr}_\pm)=0$ and $d_{\nu\nu}(\lambda^\textnormal{dr}_\pm,\nu^\textnormal{dr}_\pm)=0$ requires $s=0$, which completes the proof.
\end{proof}
\begin{Remark}\label{rem:OptimalWeight}
The absolute spectrum $\Sigma_{\textnormal{abs}}(\mathcal{L}^-)$ coincides with the weighted essential spectrum $\Sigma_{\textnormal{ess}}^{\eta^*}(\mathcal{L}^-)$ for $\eta^*\coloneqq-(\alpha s)/2$.
\end{Remark}
\begin{Remark}\label{rem:ComplexGinzburgLandau}
We highlight the structural relation to the extensively studied complex Ginzburg-Landau equation (CGL), whose spectra are well-known, see e.g.~\cite{aranson2002world} for details. The CGL in a co-moving frame as well as co-rotating frame with speed $s$ and frequency $\Omega$ reads
\begin{equation*}\label{eq:ComplexGinzburgLandauEquation}
\partial_t A=(1+\mathbf{i}\gamma)\partial_\xi^2 A +s\partial_\xi A + A - (1+\mathbf{i}\kappa)|A|^2 A-\mathbf{i} \Omega A, \tag{CGL}
\end{equation*}
where $\xi=x-st$, $A=A(\xi,t)\in\mathbb{C}$ and the real parameters $\gamma, \kappa$ characterize the linear and nonlinear dispersion. In terms of real and imaginary parts, $A=u+\mathbf{i} v$, the linearization around the zero state reads
\begin{align*}
\partial_t \begin{pmatrix}
u\\ v
\end{pmatrix}=
\begin{pmatrix}
\partial_\xi^2 + s\partial_\xi + 1 & - (\gamma\partial_\xi^2 -\Omega)\\
\gamma \partial_\xi^2 -\Omega & \partial_\xi^2 +s\partial_\xi + 1
\end{pmatrix}\begin{pmatrix}
u\\ v
\end{pmatrix},
\end{align*}
which is of the form~\eqref{eq:GeneralMatrixOperator} with $a_2=a_0=1$, $a_1=s$, $b_2=\mathbf{i}\gamma $, $b_1=0$, and $b_0=-\mathbf{i}\Omega$. Since $a_2$ and $b_2$ satisfy the conditions of Hypothesis~\ref{hyp:Coefficients}, we directly obtain from Theorem~\ref{theo:AbsolutSpectrumGeneralLinearMatrixOperator} that the absolute spectrum of (CGL) linearized around the zero state $A=0$ is given by
\begin{align*}
\Sigma_{\textnormal{abs}}&=\bigg\lbrace -(1+\mathbf{i} \gamma)k - \frac{s^2}{4(1+\mathbf{i} \gamma)}+1-\mathbf{i} \Omega, \, k\in[0,\infty) \bigg\rbrace
\\ &\,\cup \bigg\lbrace -(1-\mathbf{i} \gamma)k - \frac{s^2}{4(1-\mathbf{i} \gamma)}+1+\mathbf{i} \Omega, \, k\in[0,\infty) \bigg\rbrace
\end{align*}
with maximal real part located at $-s^2/(4(1+\gamma^2))+1$. It is well-known (see e.g.~\cite{goh2020spectral, van2003front}), that $\Sigma_{\textnormal{abs}}$ is marginally stable for $|s|=2\sqrt{1+\gamma^2}$.
\end{Remark}

\subsection{Linear spreading speed and frequency}\label{sec:LinearSpreadingSpeed}

It has been shown in~\cite{holzer2014criteria} that for any simple pinched double root $(\lambda^\textnormal{dr}, \nu^\textnormal{dr})$ it holds that $\lambda^\textnormal{dr}$ is a singularity of $\boldsymbol{G}_\lambda$. Indeed there are examples that simplicity is in general necessary~\cite{holzer2014criteria}. Thus, the most definite indicator for the transition from convective to absolute instability is the location of the most unstable simple pinched double root~\cite{holzer2014criteria}.
\begin{definition}\label{def:LinearSpreadingSpeed}
The (positive) \emph{linear spreading speed} is
\[s^\textnormal{lin}\coloneqq\sup \,\lbrace s\,:\, \mathscr{D} \textnormal{ possesses a simple pinched double root with } \mathrm{Re}(\lambda^\textnormal{dr})>0 \rbrace.\]
\end{definition}
The set in this definition is non-empty in our case of unstable essential spectrum and the spreading speed is bounded, cf.~\cite{holzer2014criteria}. One analogously defines a negative spreading speed via the infimum and for reflection symmetric problems, as ours, these have the same absolute value. Based on Corollary~\ref{cor:AbsoluteSpectrum}, we directly obtain the following result.
\begin{lemma}\label{lem:LinearSpreadingSpeed}
For any $\alpha>0, \beta>0, \mu<0, c_{\rm cp}\in(-1,1)$, the positive linear spreading speed associated to $\mathcal{L}^-$ is
\begin{equation}\label{eq:SelectedSpeedMonostable}
    s^{\textnormal{lin}}=2\sqrt{\frac{h+\mu-\boldsymbol{\beta}^-/\alpha}{1+\alpha^2}}
\end{equation}
in the $\mathcal{M}^+$ regime, where in particular $\boldsymbol{\beta}^-/\alpha-\mu<h$ holds.
\end{lemma}
\begin{proof}
Due to Corollary~\ref{cor:AbsoluteSpectrum}, the two points in $\Sigma_{\textnormal{abs}}(\mathcal{L}^-)$ with maximal real part are simple pinched double roots for $s\neq 0$; $\lambda^\textnormal{dr}_\pm$ with $\lambda^\textnormal{dr}_+=\overline{\lambda^\textnormal{dr}_-}$. Moreover, we readily see that $\mathrm{Re}(\lambda^\pm_\textnormal{dr})=0$ has positive unique solution $s=s^{\textnormal{lin}}$. Based on Definition~\ref{def:LinearSpreadingSpeed}, we obtain the statement.
\end{proof}

\begin{Remark}\label{rem:SymmetryLinearSpreadingSpeed}
Recall that fronts, in particular DWs, that travel with the linear spreading speed are called pulled. By reflection symmetry, the negative spreading speed is $-s^\textnormal{lin}$. A pulled DW with stable state $+\boldsymbol{e}_3$ to the left and unstable state $-\boldsymbol{e}_3$ to the right will propagate with $s^\textnormal{lin}$, whereas it will propagate with $-s^\textnormal{lin}$ if the stable state is to the right and the unstable state to the left, both in the $\mathcal{M}^+$ parameter regime (see also \S~\ref{sec:StandingDomainWalls} for discussion on the spreading direction). For parameters in $\mathcal{M}^-$, one has to interchange the directions of propagation.
\end{Remark}

Setting $s=s^\textnormal{lin}$, all singularities of $\boldsymbol{G}_\lambda$ lie in the closed left half-plane and one can compute~\eqref{eq:LinearSolution}. Thus, we directly obtain the following definition for the solution to be also stationary regarding rotation.
\begin{definition}\label{def:LinearSpreadingFrequency}
The \emph{linear spreading frequency} for $\mathcal{L}$ is $\Omega^\textnormal{lin}\coloneqq \mathrm{Im}\left(\lambda^\textnormal{dr}\big\vert_{\Omega=0}\right)$.
\end{definition}

This reflects the fact that in a frame moving with speed $s^\textnormal{lin}$, critical perturbations in the weighted space with weight $\eta^*$ as specified in Theorem~\ref{theo:AbsolutSpectrumGeneralLinearMatrixOperator} oscillate for $\Omega=0$ with frequency $\mathrm{Im}(\lambda^\textnormal{dr})$. Hence, in a frame moving and rotating with the linear spreading speed as well as linear spreading frequency, critical modes of~\eqref{eq:LinearSolution} are stationary with respect to the two symmetries; translation and rotation, see Figure~\ref{fig:SpecsOfDownMagnetization}.

Again based on Corollary~\ref{cor:AbsoluteSpectrum}, we obtain the  following statement.
\begin{lemma}\label{lem:LinearSpreadingFrequency}
For any $\alpha>0, \beta>0, \mu<0, c_{\rm cp}\in(-1,1)$, the linear spreading frequency associated to $\mathcal{L}^-$ in the $\mathcal{M}^+$ regime is
\begin{equation}\label{eq:SelectedFrequencyMonostable}
\Omega^\textnormal{lin}=\frac{2h+2\mu+\alpha \boldsymbol{\beta}^- - \boldsymbol{\beta}^-/\alpha}{1+\alpha^2}
\end{equation}
\end{lemma}
\begin{proof}
By Lemma~\ref{lem:LinearSpreadingSpeed}, the linear spreading speed is given by~\eqref{eq:SelectedSpeedMonostable} and from Corollary~\ref{cor:AbsoluteSpectrum}, we compute that $\mathrm{Im}(\lambda_\pm^\textnormal{dr}\big\vert_{\Omega=0})=0$ for $s=s^\textnormal{lin}$ has unique solution $\Omega^\textnormal{lin}$. The statement follows now from Definition~\ref{def:LinearSpreadingFrequency}.
\end{proof}

\begin{Remark}
Although the sign of the linear spreading speed depends on whether the invaded state is on the left or on the right (see Remark~\ref{rem:SymmetryLinearSpreadingSpeed}), the linear spreading frequency is positive within the parameter space to $\mathcal{M}^+$, where in particular $\boldsymbol{\beta}^-/\alpha-\mu<h$ holds. Additionally, for $\alpha=1$, $\Omega^\textnormal{lin}$ is independent of the parameter $c_{\rm cp}$ (recall $\beta>0$ in case of spin-torque), whereas $s^\textnormal{lin}_\pm$ depends on $c_{\rm cp}$. Moreover, we remark that in case $s=s^\textnormal{lin}_\pm$ and $\Omega=\Omega^\textnormal{lin}$ one obtains from~\eqref{eq:DoubleRootLambda} that $\lambda^\textnormal{dr}_\pm=0$ and that solving $d^-(0,\nu)=0$ leads to the four spatial solutions $\nu_j=-\alpha s^\textnormal{lin}_\pm/2 \pm \mathbf{i} s^\textnormal{lin}_\pm/2$. This is a `doubled' double root in the sense of the two symmetries: translation and rotation.
\end{Remark}

\begin{figure}
    \centering
    \begin{subfigure}[b]{0.24\textwidth}
    \includegraphics[width=\textwidth]{figures/specs_large_s_large_Omega.eps}
    \caption*{(a)}
    \end{subfigure}
    \hspace*{20mm}
    \begin{subfigure}[b]{0.32\textwidth}
    \includegraphics[trim=0 70 0 70, clip, width=\textwidth]{figures/absolut_specs.pdf}
    \caption*{(b)}
    \end{subfigure}
    \caption{Spectra of $\mathcal{L}^-$ in the $\mathcal{M}^+$ parameter regime with $s>s^\textnormal{lin}_+$ and $\Omega>\Omega^\textnormal{lin}$ in (a), where blue solid line is essential and blue dashed the weighted essential spectrum with weight $\eta^*<\eta<0$. Red lines illustrate the absolute spectrum. (b) absolute spectrum for $s=s^\textnormal{lin}_+$ as defined in~\eqref{eq:SelectedSpeedMonostable} and different frequencies: red for $\Omega>\Omega^\textnormal{lin}$ and blue for $\Omega=\Omega^\textnormal{lin}$ as defined in~\eqref{eq:SelectedFrequencyMonostable}.}
    \label{fig:SpecsOfDownMagnetization}
\end{figure}

As usual in spatial dynamics, multiple spatial eigenvalues at the origin indicate the bifurcation of coherent structures. Indeed, one readily verifies that $\Omega^\textnormal{lin}=(s^\textnormal{lin})^2/2+\boldsymbol{\beta}^-/\alpha$ is the degenerate \emph{center case} in \cite{siemer2020} that implies special properties for the `spatial' coherent structure ODE~\eqref{eq:SphericalCoherentStructureODE}. In case $\Omega=\Omega^\textnormal{lin}$, it was proven (see~\cite[Theorem 2]{siemer2020} for full statement) that if $|c_{\rm cp}|\ll1$ then a family of inhomogeneous DWs (flat as well as non-flat ones) co-exist.

We again note that the linear spreading speed and linear spreading frequency of DWs in the monostable parameter regime were determined in~\cite{goussev2020dynamics} for the case $c_{\rm cp}=0$ utilizing a saddle-point approximation. While this is generally valid for scalar reaction-diffusion equations and the~\eqref{eq:ComplexGinzburgLandauEquation} equation, the saddle-point approach may fail for systems of PDEs, as already conjectured in~\cite{ebert2000front}. Roughly speaking, each component leads to a candidate for the spreading speed and possibly additional ones in case of (linear) coupling, and one has to pick the right one~\cite{ebert2000front, van2003front}. However, the situation of~\eqref{eq:LLGS}, i.e., the computation of the linear spreading speed and frequency associated to $\mathcal{L}^-$, is analogous to the~\eqref{eq:ComplexGinzburgLandauEquation} case, and both approached yields the same prediction. Heuristically, this is not surprising since~\eqref{eq:LLGS} can be written as a generalized~\eqref{eq:ComplexGinzburgLandauEquation} by pointwise stereographic projection, cf. e.g.~\cite{Melcher2017}. Here we have generalized this to operators $\mathcal{L}$ that correspond to a general scalar complex equation, which means the dispersion relation~\eqref{eq:GeneralDispersionRelation} has the form $d(\lambda,\nu)=d_1^2(\lambda,\nu) + d_2^2(\nu)$ with polynomials $d_1$ and $d_2$. As shown above, it follows that there is only one candidate for the linear spreading speed (in absolute value) as well as frequency. Beside these observations, we again remark that the study of the pointwise Green's function, and in particular the locations of its singularities, can be generalized from scalar equations to systems of PDEs in a straightforward way. Moreover, the description of the absolute spectrum via weighted spectral curves can possibly be generalized to even more complicated equations, see e.g.~\cite{faye2020remnant}.

\subsection{Relation to homogeneous domain walls}\label{sec:RelationHomogeneousDW}
We return to the full nonlinear~\eqref{eq:LLGS} and the question whether an initial interface will propagate asymptotically with $s^\textnormal{lin}$ and $\Omega^\textnormal{lin}$ or not. For this, we first focus on the situation of zero spin polarization, $c_{\rm cp}=0$, and thus assume $\beta/\alpha<h$. We further select the positive speed in~\eqref{eq:HomogeneousSpeed}, denoted by $s_+^\textnormal{hom}$, by taking the unstable state $-\boldsymbol{e}_3$ on the right side. In this parameter regime, $s_+^\textnormal{hom}$ equals $s^\textnormal{lin}$ for
\begin{equation}\label{eq:AbsoluteSpectrumAppliedField}
h=h_\pm^s\coloneqq\frac{\beta}{\alpha}-2\mu-\frac{2\mu}{\alpha^2}\pm\frac{2\mu}{\alpha^2}\sqrt{1+\alpha^2},
\end{equation}
where $h^s_+<h^s_-$, since we consider the easy-axis case $\mu<0$. Concerning the boundary between the bistable and monostable regimes given by $\beta/\alpha-\mu$, i.e., the left boundary of $\mathcal{M}^+$ in case $c_{\rm cp}=0$, one readily verifies that $\beta/\alpha-\mu<h^s_+$. This means that the linear spreading speed is smaller than the speed of $\boldsymbol{m}_0$ for $h\in(\beta/\alpha-\mu,h^s_+)$, so that $\boldsymbol{m}_0$ may be selected rather than a DW with linear spreading speed~\eqref{eq:SelectedSpeedMonostable}. Indeed, we verify numerically in \S~\ref{sec:PointSpectrumNumerics} that $\boldsymbol{m}_0$ has stable point spectrum up to $h^s_+$, which means that in this regime $\boldsymbol{m}_0$ is convectively unstable only, since the weighted essential spectrum is stable by the above analysis. Moreover, also in direct numerical simulations, we observe the selection of $\boldsymbol{m}_0$ from a localized perturbation around $-\boldsymbol{e}_3$, so that $\boldsymbol{m}_0$ would be a pushed front in the applied field regime $h\in(\beta/\alpha-\mu,h^s_+)$, see \S~\ref{sec:NumericalSimulationMethod} for details.

Analogously, we compute that~\eqref{eq:HomogeneousFrequency} is equal to $\Omega^\textnormal{lin}$ in case
\[h=h^\Omega\coloneqq \frac{\beta}{\alpha}-2\mu.\]
Since $\alpha>0$, it follows that $h_\pm^s\neq h^\Omega$ for all $\beta\ge 0$ and $\mu<0$, and thus $h^s_+<h^\Omega<h^s_-$. We summarize these findings as follows.
\begin{corollary}\label{cor:DifferenceToHomogeneousDW}
In the $\mathcal{M}^+$ regime in case $c_{\rm cp}=0$, a DW that propagates with the linear spreading speed and frequency given by~\eqref{eq:SelectedSpeedMonostable} and~\eqref{eq:SelectedFrequencyMonostable}, respectively, differs from $\boldsymbol{m}_0$ given by~\eqref{eq:AnalyticHomogeneousSolution}. In particular, such a DW is non-flat, due to~\cite[Theorem 2]{siemer2020}. By symmetry, this also holds in the $\mathcal{M}^-$ regime given by~\eqref{eq:MonostableSetMinus}.
\end{corollary}

\section{Numerical Studies}\label{sec:Numerics}
In this section, we study on the one hand the location of the point spectrum of the linear operator arising from linearization around the explicitly known family of homogeneous DWs~\eqref{eq:AnalyticHomogeneousSolution} and its dependence on the applied field. On the other hand, we describe the method of numerical time integration of~\eqref{eq:LLGS} used in this paper, and also the method of freezing from which the long-time behavior within the bistable as well as monostable parameter regime were obtained.

\subsection{Point Spectrum}\label{sec:PointSpectrumNumerics}

We are interested in the location of the point spectrum of $\boldsymbol{m}_0$ and thus set $c_{\rm cp}=0$. Due to symmetry, we only consider the case $\beta/\alpha<h$ and hence $\sigma=1$ in~\eqref{eq:AnalyticHomogeneousSolution}. The associated speed and frequency is given by~\eqref{eq:HomogeneousSpeed} as well as~\eqref{eq:HomogeneousFrequency}.

Stability of $\boldsymbol{m}_0$ was numerically verified in~\cite{gou2011stability} for $h<\beta/\alpha-\mu$, where the aforementioned energy difference $\Delta E_\varepsilon$ was used. For applied fields beyond $\beta/\alpha-\mu$, however, $\boldsymbol{m}_0$ is unstable in $L^2$ -- also in case of stable point spectrum -- due to the presence of unstable essential spectrum. As discussed in Sections~\ref{sec:Pointwise growth, absolute spectra, and double roots} and~\ref{sec:RelationHomogeneousDW}, in the convectively unstable regime $|h|\notin(h^s_+,h_-^s)$, the essential spectrum is stable in the exponentially weighted spaces $L^2_\eta$ for suitable $\eta<0$. In particular, the absolute spectrum $\Sigma_{\textnormal{abs}}$ is stable, and in the following we distinguish the parameter regimes $\beta/\alpha-\mu<h<h^s_+$ and $h^s_-<h$.

Numerical computations of the point spectrum were performed using the \textsc{Matlab}\xspace toolbox \textsc{Stablab}\xspace . This toolbox uses the so-called \emph{Evans function} to numerically locate eigenvalues, and we refer to~\cite{barker2018evans, barker2009stablab} for further details. For implementation in \textsc{Stablab}\xspace , we rewrite~\eqref{eq:LinearEigenvalueProblem}, including exponential weights, as a first order system of the form
\begin{equation}\label{eq:WeightedFirstOrderSystem}
W'=A_\eta(\xi;\lambda)W,
\end{equation}
where $A_\eta(\xi;\lambda)$ is defined in the Appendix. The fixed parameters are $\alpha=1$, $\beta=3/4$, and $\mu=-1$, which yields $\beta/\alpha-\mu=1.75$ and $h^s_+= 1.92$. The domain is set to $\xi \in [-100,100]$ and since we know the locations of $\Sigma_{\textnormal{ess}}$, $\Sigma_{\textnormal{ess}}^\eta$, and $\Sigma_{\textnormal{abs}}$ explicitly (see \S~\ref{sec:StabilitySteadyStates} and~\ref{sec:Pointwise growth, absolute spectra, and double roots}), we can choose an integration contour that excludes zero eigenvalues which arise from the translation and rotation invariances.

Within the first parameter region, we take the samples $h\in\lbrace 1.75, 1.8, 1.85, 1.9\rbrace$ with $\eta\in\lbrace 0, -0.1, -0.15, -0.25 \rbrace$. Due to the relatively small range, $1.75\le h\le 1.9$, the individual Evans function images of the chosen contour do not vary significantly. Thus, we only illustrate one exemplary Evans function output for $h=1.9$ with $\eta=-0.25$ in Figure~\ref{fig:EvansFunctionOutput} (b). Here we utilize the normalization option within \textsc{Stablab}\xspace for illustration purposes. In all these cases, the winding number is zero, demonstrating stability of the point spectrum of $\boldsymbol{m}_0$ for $h<h^s_+$ in $L^2_\eta$.

For the second region, we choose $h\in\lbrace 8, 15, 20\rbrace$ and $\eta\in\lbrace -1.5, -1.2, -1.1 \rbrace$. The images of the same contour as before (cf.\ Figure~\ref{fig:EvansFunctionOutput} (a)) under the resulting Evans functions are illustrated for $h=8$ as well as $h=20$ in Figure~\ref{fig:EvansFunctionOutputLargeH}. In all cases, the winding number is non-zero, indicating unstable point spectrum. For an applied field in this region, i.e., beyond $h^s_-$, we also observe numerically in \S~\ref{sec:NumericalSimulationMethod} that not $\boldsymbol{m}_0$, but a critical pulled front is dynamically selected (cf.\ Figure~\ref{fig:SpeedAndFrequencyComparison}). This selection behavior persists for $c_{\rm cp}\neq 0$ (cf.\ Figure~\ref{fig:SelectedSpeedAndFrequency}), where pushed DWs are asymptotically selected also for small applied fields, and where pulled DWs are selected for $h$ beyond a certain threshold. It is clear that this threshold lies above $\boldsymbol{\beta}^-/\alpha-\mu$, but more specific analytic results do not seem to be known.

\begin{figure}
    \centering
    \begin{subfigure}[b]{0.37\textwidth}
    \includegraphics[trim=50 0 90 30, clip, width=\textwidth]{figures/contour2.eps}
    \caption*{(a)}
    \end{subfigure}
    \hspace*{10mm}
    \begin{subfigure}[b]{0.37\textwidth}
    \includegraphics[trim=50 30 120 30, clip, width=\textwidth]{figures/evans2.eps}
    \caption*{(b)}
    \end{subfigure}
    \caption{Shown are the Evans-function computations for~\eqref{eq:WeightedFirstOrderSystem} with $\alpha=1, \beta=0.75, \mu=-1, h=1.9, \eta=-0.29, s=0.575,$ and $\Omega=1.325$. Semi-circular contour of radius $100$ in (a) using $1500$ mesh points, excluding the half-circle of radius $0.1$ around the origin shown in the inset. The normalized image of the contour under the Evans function is shown in (b), where the red cross indicates the origin and the inset zooms in around the origin. The winding number is zero.}
    \label{fig:EvansFunctionOutput}
\end{figure}

\begin{figure}
    \centering
    \begin{subfigure}[b]{\textwidth}
    \centering
    \includegraphics[trim=20 0 100 30, clip, width=0.3\textwidth]{figures/evans8.eps}
    \hfill
    \includegraphics[trim=20 0 100 30, clip, width=0.3\textwidth]{figures/evans8zoom.eps}
	\hfill
    \includegraphics[trim=20 0 100 30, clip, width=0.3\textwidth]{figures/evans8zoomzoom.eps}
    \caption*{$h=8$}
    \end{subfigure}
    \vskip\baselineskip
    \begin{subfigure}[b]{\textwidth}
    \centering
    \includegraphics[trim=20 0 100 30, clip, width=0.3\textwidth]{figures/evans20.eps}
    \hfill
    \includegraphics[trim=20 0 100 30, clip, width=0.3\textwidth]{figures/evans20zoom.eps}
	\hfill
    \includegraphics[trim=20 0 100 30, clip, width=0.3\textwidth]{figures/evans20zoomzoom.eps}
    \caption*{$h=20$}
    \end{subfigure}
    \caption{The left panels show the images of the contour Figure~\ref{fig:EvansFunctionOutput} (a) under the Evans function for $h=8$ (upper) and $h=20$ (lower). The middle and right panels show enlargements of the region around the origin, which is indicated by a red cross. Other parameters are as in Figure~\ref{fig:EvansFunctionOutput}. The winding number is six for $h=8$, and fourteen for $h=20$.}
    \label{fig:EvansFunctionOutputLargeH}
\end{figure}

\begin{Remark}
Since we are interested in the location of the point spectrum of solutions on the unit sphere, one could also start from linearizing~\eqref{eq:SphericalCoherentStructurePDE} around $\boldsymbol{m}_0$ in spherical coordinates which reads $(\theta_0,\phi_0)=(2\arctan(\exp(\sqrt{-\mu}\xi)),0)$. However, due to the singular behavior of the term $\sin^{-1}(\theta_0)$ for $|\xi|\to\infty$, convergence is not guaranteed and \textsc{Stablab}\xspace cannot be applied in a straightforward way. Therefore, we chose to study the location of the point spectrum in the original coordinates.
\end{Remark}

\subsection{Numerical simulation method and freezing}\label{sec:NumericalSimulationMethod}~

Before discussing our results from numerical simulations, we comment on the methods used. The results of numerical time integration of~\eqref{eq:LLGS} rely on a semi-implicit method with a finite element scheme in space from the \textsc{Matlab}\xspace package \textsc{Pde2Path}~\cite{dohnal2014pde2path}. In order to capture the spreading speeds and frequencies throughout the simulations, we extend the implementation of the `freezing' method from~\cite{BeynThuemmler2004} as discussed in~\cite{rademacher2017symmetries}, such that both symmetries are tracked. Briefly, here we first project onto the two generators and secondly solve the two component equation system for the current translation as well as current rotation, both within each time step. In order to ensure $\boldsymbol{m} \in \mathbb{S}^2$, in each time step we scale the resulting predictor pointwise onto the unit sphere. In this way, we perform simulations with a typical domain size $[-50,50]$ with grid spacing $\Delta_x=\num{e-2}$, homogeneous Neumann boundary conditions. The final time within the bistable parameter regime is $t=50$, and $t=100$ within the monostable parameter regime, both with step size $\Delta_t=\num{e-4}$. The final time is increased for $h$ close to $h^s_+$. The initial conditions within the bistable parameter regime are of step function-type, connecting the state $+\boldsymbol{e}_3$ on the left and $-\boldsymbol{e}_3$ on the right, whereas the initial conditions in the monostable parameter regime are perturbations locally around the unstable state $-\boldsymbol{e}_3$. Parameters are again $\alpha=1, \beta=0.75$, and $\mu=-1$. For an applied field less than $h^\Omega=2.75$ and close to $h^s_+=1.92$, we observe modulated fronts which are characterized by oscillating frozen speed and frequency (cf.\ Figure~\ref{fig:ModulatedFront}). Although these oscillations mostly decay for a final time extended to $t=500$, the decay significantly slowed down for an applied field quite close to $h^s_+$. Therefore, we do not plot a few points in the figures for which there seems to be no convergence at $t=500$.

Our results from these numerical simulations confirm the analytical predictions and also the persistence for $|c_{\rm cp}|\ll 1$, even for $c_{\rm cp}$ far away from zero. The results for $c_{\rm cp}=0$ are presented in Figure~\ref{fig:SpeedAndFrequencyComparison} and for $c_{\rm cp}=\pm 0.5$ in Figure~\ref{fig:SelectedSpeedAndFrequency}.

The frozen speeds and frequencies in case $c_{\rm cp}=0$, but also for $c_{\rm cp}\neq 0$, approach the linear predictions~\eqref{eq:SelectedSpeedMonostable} and~\eqref{eq:SelectedFrequencyMonostable} from below for $h^\Omega< h$, where no modulated fronts are observed. We suspect that this is a consequence of the logarithmic correction, which is known to be associated with the spreading speed selection in scalar reaction-diffusion equations (see for example~\cite{ebert2000front, hamel2013}).
In our numerical results also the spreading frequency appears to be approached monotonically from below, which we suspect is due to the functional relation of spreading speed and frequency. For $h^s_+<h<h^\Omega$, this behavior does not persists in a strict sense, due to the aforementioned oscillations in $s$ and $\Omega$ (see Figure~\ref{fig:ModulatedFront} for an example).

\begin{figure}
\centering
\begin{tabular}{cc}
  \begin{minipage}[c]{.33\textwidth}
  \centering
  \includegraphics[trim=0 0 0 150, clip, width=\textwidth]{figures/modulatedfrontsphere.eps}
\end{minipage}
\begin{minipage}[c]{.52\textwidth}
  \begin{subfigure}{\textwidth}
  \centering
  \includegraphics[trim=0 0 10 0, clip,width=0.6\textwidth]{figures/modulatedfrontspeed.eps}
  \end{subfigure}

  \vspace*{0mm}
  \begin{subfigure}{\textwidth}
  \centering
  \includegraphics[width=0.6\textwidth]{figures/modulatedfrontfrequency.eps}
  \vspace*{-5mm}
  \end{subfigure}
\end{minipage}

\end{tabular}
\vspace*{5mm}
\caption{Left panel: DW at final time $t=500$ projected onto $\mathbb{S}^2$. Right panels: temporal evolution of the speed (upper) and frequency (lower) for $t\in[100,500]$. Parameters are $h=2.4$ and $c_{\rm cp}=0$, other parameters are as in Figure~\ref{fig:StabilityRegionsUp-AndDown-Magnetization}. Linear predictions are $s^\mathrm{lin}=1.14$ as well as $\Omega^\mathrm{lin}=1.4$.}
\label{fig:ModulatedFront}
\end{figure}

For $h<\beta/\alpha-\mu$, we find that $\boldsymbol{m}_0$ is selected, which is expected due to stable spectrum (see \S~\ref{sec:StabilitySteadyStates} and \S~\ref{sec:PointSpectrumNumerics}). For $\beta/\alpha-\mu<h<h^s_+$ in the monostable regime, a nonlinear selection mechanism appears to dominate the linear spreading speed, which is plausible due to stable point spectrum numerically observed in \S~\ref{sec:PointSpectrumNumerics}. This is a well-known phenomenon in scalar reaction diffusion equations, cf.\ e.g.~\cite{bramson1983convergence, ebert2000front, garnier2012inside, kirchgassner1992nonlinear, van2003front}, and so these DWs appear to be pushed fronts, as mentioned before. For $h^s_+<h$, fronts with speed and frequency given by the linear predictions, i.e., $s^\mathrm{lin}$ and $\Omega^\mathrm{lin}$, are numerically selected over time. These dynamically selected DWs must differ from $\boldsymbol{m}_0$ due to Corollary~\ref{cor:DifferenceToHomogeneousDW}. This behavior persists also for $c_{\rm cp}$ far away from zero (cf.\ Figure~\ref{fig:SelectedSpeedAndFrequency}).

\begin{figure}[h]
    \centering
    \includegraphics[trim=50 0 0 0, clip, width=0.8\textwidth]{figures/freezingzeroccpfull.eps}
    \caption{Speed and frequency selection for $h\in(\beta/\alpha, \beta/\alpha-10\mu)$ in case $c_{\rm cp}=0$. Other parameters are as in Figure~\ref{fig:StabilityRegionsUp-AndDown-Magnetization}. Speed and frequency of $\boldsymbol{m}_0$, cf.~\eqref{eq:HomogeneousSpeed} (dashed blue) as well as~\eqref{eq:HomogeneousFrequency} (dashed red), and linear spreading speed~\eqref{eq:SelectedSpeedMonostable} (solid blue) and spreading frequency~\eqref{eq:SelectedFrequencyMonostable} (solid red). Numerically observed (asymptotic) speeds and frequencies are indicated by blue crosses and red circles, respectively. Vertical lines mark $h=\beta/\alpha-\mu=1.75$ (solid), $h^s_+=1.92$ (dashed), $h^\Omega=2.75$ (dotted dashed), and $h^s_-=7.58$ (dotted).}
    \label{fig:SpeedAndFrequencyComparison}
\end{figure}

\begin{figure}[h]
\centering
   	\begin{subfigure}[b]{0.49\textwidth}
		\includegraphics[width=\textwidth]{figures/negativeonehalfccp2.eps}
	\caption*{\small $c_{\rm cp}=-0.5$}
	\end{subfigure}
	\hfill
	\begin{subfigure}[b]{0.49\textwidth}
		\includegraphics[width=\textwidth]{figures/positiveonehalfccp2.eps}
	\caption*{\small $c_{\rm cp}=0.5$}
	\end{subfigure}
	\caption{Linear spreading speed~\eqref{eq:SelectedSpeedMonostable} (solid blue) and spreading frequency~\eqref{eq:SelectedFrequencyMonostable} (solid red) for $h\in(\beta/\alpha, \beta/\alpha-5\mu)$ in case $c_{\rm cp}=-0.5$ in (a) and $c_{\rm cp}=0.5$ in (b). Other parameters are as in Figure~\ref{fig:StabilityRegionsUp-AndDown-Magnetization}. Numerically observed (asymptotic) speed and frequency is indicated by blue crosses and red circles, respectively. Black vertical line separates the bistable region to the left from the monostable region to the right, cf.\ Figure~\ref{fig:FreezingZeroCCPLight}.}
	\label{fig:SelectedSpeedAndFrequency}
\end{figure}

\section{Discussion and Outlook}\label{sec:DiscussionAndOutlook}
In this paper, we have investigated the dynamics of DWs that connect the up- and down-magnetization states in axially symmetric models of spintronic nanowires. Taking the applied field as the main parameter, we have initially separated the parameter range into the bistable and the monostable region.

For the bistable regime, and without spin-torque effect, we have analytically extended the threshold for stability of an explicitly known family of homogeneous DWs and also into the region where spin-torque effects are present. We have also verified the stability of this family beyond this threshold numerically up to (and beyond) the transition into the monostable regime. Moreover, we have presented a heuristic formula which relates the applied field and polarization ratio that allows to identify parameters for (nearly) standing DWs. We have numerically corroborated that this prediction approximates the direction of propagation well for a broad range of $c_{\rm cp}\neq 0$. We believe that this can be a basis to study weak interactions of multiple stacked DWs in a nanowire that individually would be stationary, but that interact through their tails and thus are expected to form metastable states.

Concerning the monostable regime, we have first presented an explicit formula for the absolute spectrum for a class of complex matrix operators for two-component systems. From this, we were able to relate the most unstable points of the absolute spectrum to the location of the essential singularities of the associated pointwise Green's functions. This allowed us to explicitly determine the linear spreading speed and the associated frequency. We have numerically verified that, for increasing $h$, the explicit DWs are merely convectively unstable until the point of absolute instability of an asymptotic state. Specifically, we find that until this point the explicit family possesses stable point spectrum and also stable essential spectrum in a weighted $L^2$ space. At the same time, this provides numerical evidence that the explicit family yields pushed fronts in this parameter regime, whereas for larger values of $h$ pulled fronts of unknown analytic form are selected. The latter ones are certainly of inhomogeneous DW type and we expect these to be also non-flat, which is the case for zero spin-torque. Moreover, while the explicit family of DWs propagates faster than the linear prediction for applied fields beyond a larger threshold, we find that in this regime the explicit family possesses unstable point spectrum and thus cannot be dynamically selected.

In order to illustrate and corroborate these results, we have presented numerical long-time simulations for applied fields in both regimes; bistable and monostable, using the method of freezing that filters motion in the continuous symmetry groups by selecting the appropriate speed and rotation frequency. We have found that these principles of selecting DWs persist in case of spin-torque effects, i.e., $0<|c_{\rm cp}|<1$.

A natural step towards a deeper understanding of domain wall motion in nanowires is on the one hand to study the selection mechanism of fronts that propagate with the linear predictions at least for the case $c_{\rm cp}=0$ via the energy functional. On the other hand, an Evans function approach could be utilized to analytically verify the stability of the explicit family of homogeneous DWs. A further interesting problem to study, for which one could start from the case $c_{\rm cp}=0$, are the properties and selection principles of modulated fronts which numerically seem to appear close to $h^\Omega$.

Since the main results are also valid for $c_{\rm cp}=0$, they transfer over to the well-known~\eqref{eq:LLG} equation.

\section*{Acknowledgments}
The authors acknowledge support by the Deutsche \linebreak Forschungsgemeinschaft (German Research Foundation), Projektnummer \linebreak 281474342/GRK2224/1. The authors gratefully thank Matt Holzer for stimulating discussions related to \S~\ref{sec:Pointwise growth, absolute spectra, and double roots}, as well as Blake Barker for support on \textsc{Stablab}\xspace .