\section{Introduction}
Helical structures lie at the heart of several feats of innovation.
The Archimedean screw and spiral staircases are inventions
dating back to ancient history. Simple bolts and spring coils stand as hallmarks
of practicality in modern everyday life. Helical growth is a clear demonstration
that such structures also exist naturally without human intervention.
On the other hand, helices are vital for life itself as we know it,
the DNA polymer being a famous example of a double-helix structure~\cite{Watson1953}.

The study of helical spin textures has recently drawn attention in the field of
gaseous spinor Bose--Einstein condensates (BECs). In a quantum degenerate gas of ferromagnetic spin-1 ${}^{87}{\rm Rb}$,
a helical magnetization texture was observed to decay into small spin domains~\cite{Vengalattore2008}.
This effect was argued to result from the weak interatomic magnetic dipole forces in the system.
Closely related to the experiment, the dynamical instability of an $XY$ spiral state has been investigated theoretically~\cite{Cherng2008}.

It has been predicted that a sufficiently strong magnetic dipolar forces can spontaneously give rise to
intriguing spin textures in ferromagnetic condensates: the long-range dipolar potential stabilizes
spin-vortex states in various geometries, as was demonstrated using a \mbox{spin-1} model~\cite{Yi2006,Kawaguchi2006}
and a classical spin approach~\cite{Takahashi2007,Huhtamaki2010}.
Even weak dipolar interactions can lead to dynamic formation
of a helical spin texture in a ferromagnetic spin-1 BEC~\cite{Kawaguchi2007}.
In the absence of external magnetic fields, the spin helix can appear as the ground-state texture in a suitable geometry and with strong enough dipolar
interactions~\cite{Huhtamaki2010}.

In the present article, we investigate helical spin textures
in ferromagnetic dipolar condensates using a spin-1 model. We find stationary helical
states under the assumption that the system is infinitely long in the direction of the helix axis and find solutions hosting different types
of topological line defects, i.e., quantized vortices.
Due to the symmetry of the order-parameter field, these line defects
encircle the condensate in a helical pattern. The resulting vortex structure resembles a pair of vortices
excited in Kelvin modes in a stationary configuration.
Direct experimental evidence of Kelvin waves has been observed in a spin-polarized BEC of ${}^{87}{\rm Rb}$
by studying two transverse quandrupole modes of the atomic cloud~\cite{Bretin2003}. Kelvin-wave excitations of quantized vortices
have also been studied theoretically in a similar setup~\cite{Mizushima2003,Fetter2004,Simula2008}. Moreover, helical vortices in a two-component
condensate have been investigated~\cite{Cho2005}.

Although we focus on a spin-1 BEC, the study can also shed light on phenomena in strongly dipolar systems with more complex order parameters.
The energy-minimizing helical spin textures that we obtain are in qualitative agreement with those found earlier using the classical
spin approach~\cite{Huhtamaki2010}. Since
the classical spin model is supposed to be accurate for ferromagnetic systems in the limit of large magnetic moments, it is plausible to expect that similar
structures exist independent of the value of the atomic spin number $F$.
The most timely example of a condensate with larger $F$ and significant dipolar interactions is the spin-3 gas of ${}^{52} {\rm Cr}$~\cite{Griesmaier2005} which
has been produced by purely optical means~\cite{Beaufils2008}. The chromium atoms have magnetic moments of $6\, \mu_\mathrm{B}$, whereas
the maximal atomic magnetic moment for an alkali-metal condensate is $1\, \mu_\mathrm{B}$. Recently, there has been progress in cooling and trapping
vapors of $\mathrm{Tm}$ $(4\, \mu_\mathrm{B})$~\cite{Sukachev2010}, $\mathrm{Er}$ $(7\, \mu_\mathrm{B})$~\cite{Berglund2008},
and $\mathrm{Dy}$ $(10\, \mu_\mathrm{B})$~\cite{Lu2010}, the last having the largest atomic magnetic moment of all known elements.

\section{Theory}
In this work, we study ferromagnetic spin-1 BECs, such as
${}^{87}{\rm Rb}$, using a zero-temperature mean-field model. In
addition to the density--density and spin--spin interatomic forces,
we also include long-range dipolar interactions in the model.

Stationary states of the system are solutions to the
time-independent Gross--Pitaevskii (GP) equation
\begin{eqnarray}
\label{GP} \hat{h} \Psi({\bf r}) &+& g n({\bf r}) \Psi({\bf r}) + g_{\mathrm s} \sum_\alpha
M_\alpha({\bf r}) \hat{S}_\alpha \Psi({\bf r}) \\
&+&g_{\mathrm d} \sum_{\alpha,\beta} \int d^3r' D_{\alpha \beta}({\bf r}-{\bf r}')
M_\alpha({\bf r}') \hat{S}_\beta \Psi({\bf r}) = 0, \nonumber
\end{eqnarray}
where $\Psi=\left(\psi_1, \psi_0, \psi_{-1} \right)^T$ is a three-component spinor order parameter,
$\hat{h}=-\hbar^2 \nabla^2/2m+V_{\rm trap}({\bf r})-\mu$ is the
single-particle Hamiltonian, and $\hat{S}_\alpha$ denotes the
$\alpha$th component of the dimensionless $F=1$ spin operator
whose spin-space expectation value gives the $\alpha$th component of magnetization,
$M_\alpha({\bf r})=\Psi^\dagger({\bf r}) \hat{S}_\alpha \Psi({\bf r})$. The density of particles is given by
$n({\bf r})=\sum \psi_k^*({\bf r}) \psi_k({\bf r})$, where the components
$\left(\psi_1, \psi_0, \psi_{-1} \right)$ are the projections of
$\Psi$ onto the eigenbasis of $\hat{S}_z$. The total number of
particles per unit length, $\int d^2 r\, n({\bf r})=N$, is controlled through the chemical
potential $\mu$ acting as a Lagrange multiplier.

The coupling constants $g$, $g_{\mathrm s}$ and $g_{\mathrm d}$
measure the strengths of the local density--density, local
spin--spin, and non-local magnetic dipole--dipole interactions,
respectively. The first two are related to the scattering lengths $a_0$
and $a_2$ into spin channels with total spin $0$ and $2 \hbar$
through $g=4\pi \hbar^2 (a_0+2a_2)/3m$ and $g_{\mathrm s}=4\pi
\hbar^2 (a_2-a_0)/3m$. Throughout the work, we use
$g_{\mathrm s}=-0.01\,g$, which is roughly the coupling
constant for ${}^{87}{\rm Rb}$
\cite{Klausen2001,Kempen2002,Widera2006} and a value previously used,
e.g., in \cite{Yi2006,Simula2010}. The dipolar coupling constant is
given by $g_{\mathrm d}=\mu_0 \mu_{\mathrm B}^2 g_{\mathrm
F}^2/4\pi$ with $\mu_0$, $\mu_{\mathrm B}$, and $g_{\mathrm F}$
being the permeability of vacuum, the Bohr magneton, and the Land\'e
factor, respectively. Rather than fixing $g_{\mathrm d}$ to some
particular value, e.g., $g_{\mathrm d} \sim 10^{-3} g$ as for
${}^{87}{\rm Rb}$, we present results for various interaction
strengths in order to emphasize the role of dipolar effects. Moreover,
by using greater values of $g_\mathrm{d}/g$,
the study should also provide useful information for systems subject to
strong dipolar forces, such as gases of ${}^{52}{\rm Cr}$
\cite{Lahaye2009}. In experiments, the ratio $g_\mathrm{d}/g$ may be controlled
with an optical Feshbach resonance~\cite{Fedichev1996}, which has been
demonstrated for ${}^{87}{\rm Rb}$~\cite{Theis2004}.

The long-range dipolar interactions are characterized by the
functions $D_{\alpha \beta}({\bf R})=\left(\delta_{\alpha \beta}
R^2-3 R_\alpha R_\beta \right)/R^5$, where $\{ R_\alpha\}$ denote
the components of the argument ${\bf R}$ and $R=\sqrt{\sum_\alpha R_\alpha^2}$.
The elements of the traceless symmetric tensor $D_{\alpha \beta}$
have a $d$-wave symmetric form and are thus simply expressed in
cylindrical coordinates. The interaction integral in Eq.~(\ref{GP})
can be viewed as an effective potential for the order parameter
$\Psi$ at ${\bf r}$ arising from the magnetization throughout the
system.

In the present work, we concentrate on elongated helical solutions
to Eq.~(\ref{GP}), similar to the helical textures studied recently
in cigar-shaped systems using a classical spin
approximation~\cite{Huhtamaki2010}. For simplicity, we assume an
infinitely long system in the axial direction, confined radially by
a harmonic potential $V_{\rm trap}=\frac{1}{2}m\omega_r^2 r^2$, where
$r=\sqrt{x^2+y^2}$ and $\omega_r$ is the radial trapping frequency.
The axial symmetry of the trap allows a well-defined wave vector
$\kappa$ for the helical texture. The aim is to fix $\kappa$ and
calculate the energy-minimizing texture for a given radial plane,
say, for $z=0$, subject to the condition that the planar texture is
mapped along the axial direction according to the helical structure.
Therefore, we write the order parameter in the general form
\begin{equation}
\label{ansatz} \Psi(r,\varphi,z)=e^{i q_z z}e^{i \kappa z
\hat{S}_z}\Psi(r,\varphi+\kappa z,0)
\end{equation}
expressed in cylindrical coordinates $(r,\varphi,z)$. The second
exponential factor describes a spin-rotation of angle $\kappa z$
about the $z$-axis in the clockwise direction.
Due to the non-linear terms in Eq.~(\ref{GP}),
the order parameter feels an effective potential with a period of
$2\pi/\kappa$ in the axial direction, and thus the Ansatz is written
in the form of a Bloch-wave including the first exponential factor.

By substituting Eq.~(\ref{ansatz}) into Eq.~(\ref{GP}) and
setting $z=0$, we obtain a GP equation reduced to
polar coordinates $(r,\varphi)$,
\begin{eqnarray}
\label{reduced_GP} \hat{h}' \Psi'(r,\varphi) &+& g n'(r,\varphi) \Psi'(r,\varphi) + g_{\mathrm s} \sum_\alpha
M_\alpha'(r,\varphi) \hat{S}_\alpha \Psi'(r,\varphi) \nonumber \\
&+&g_{\mathrm d} \sum_\beta I_\beta(r,\varphi,\kappa) \hat{S}_\beta \Psi'(r,\varphi) = 0,
\end{eqnarray}
where the primes denote that the quantities are evaluated at $z=0$. The kinetic
energy in the transformed single-particle operator $\hat{h}'$ is
given by
\begin{equation}
\label{kinetic_energy} -\frac{\hbar^2}{2m}\nabla^2 +
\frac{\hbar^2}{2m}\left[q_z + \kappa \left( \hat{L}_z + \hat{S}_z
\right) \right]^2,
\end{equation}
where $\nabla^2=\partial_r^2+\partial_r/r+\partial_\varphi^2/r^2$
and $\hat{L}_z=-i\partial_\varphi$. The functions
$I_\beta(r,\varphi,\kappa)$ arising from the dipolar interactions
are evaluated in Appendix A.

Only the kinetic-energy term in the energy functional depends
explicitly on the axial wave vector $q_z$, and hence minimization of
the total energy with respect to $q_z$ yields
\begin{equation}
q_z^{\rm min} = -\frac{\kappa}{N} \int d^2 r\, \Psi^\dagger \left( \hat{L}_z +
\hat{S}_z \right) \Psi.
\end{equation}
Substitution back into Eq.~(\ref{kinetic_energy}) reveals that for
a fixed planar texture, the kinetic energy is a quadratic function
of the helix wave vector $\kappa$. Moreover, the kinetic energy
becomes independent of $\kappa$ exactly when ${\rm Var} \big(
\hat{L}_z + \hat{S}_z \big)=0$, implying that $\Psi$ must be an
eigenstate of $\hat{L}_z + \hat{S}_z$. Denoting the (integer)
eigenvalue by $\lambda$, we find that in such a case the order
parameter is of the form
$\psi_k(r,\varphi,z)=e^{i(\lambda-k)\varphi}\psi_k(r,0,z)$, which is an
integer-spin vortex state. This describes the general form of states
for which the magnetization is invariant with respect to rotations
about the $z$-axis, implying that the dipolar energy and hence also
the total energy are independent of $\kappa$.

Previous studies have shown that for very elongated condensates with
finite dipolar interaction strengths, $g_{\mathrm d} > 0$, the spin-polarized
textures with the magnetization lying along the long axis are the
energetically favored ones~\cite{Yi2006,Huhtamaki2010}. Such
states are of the general form discussed above and hence do not
depend on $\kappa$. In this study, however, we are interested in
helical textures which are found by requiring the symmetry
conditions~\cite{Huhtamaki2010}
\begin{eqnarray}
M'_{x,y}(r,\varphi+\pi) = M'_{x,y}(r,\varphi), \nonumber \\
\label{Symmetry_Conditions} M'_z(r,\varphi+\pi) = -M'_z(r,\varphi).
\end{eqnarray}
Equations~(\ref{Symmetry_Conditions}) are satisfied, e.g., when
$\psi_{k}(r,\varphi+\pi,z)=\psi^*_{-k}(r,\varphi,z)$. The textures
presented in the next section are solutions to Eq.~(\ref{GP})
subject to these symmetry conditions with respect to inversion
about the $z$-axis. A typical spin helix is shown schematically in Fig.~\ref{FIG1}.

\begin{figure}[!t]
\includegraphics[width=180pt]{FIG1.eps}
\caption{\label{FIG1} (Color online) A schematic illustration of a
single period of a typical spin helix state. Here, $p=2\pi/\kappa$ denotes
the pitch of the helix whose wave vector is $\kappa$. The solid curves depict the helical trajectories
of two points on the $z=0$ plane, and the arrows point along the local direction of magnetization. If the
state hosts vortices, their cores will also follow similar trajectories.}
\end{figure}

\section{Results}

We have solved the reduced GP equation, Eq.~(\ref{reduced_GP}), for
various values of the helical wave vector $\kappa$ and the dipolar
coupling constant $g_{\mathrm d}$. In the numerical simulations, we
choose the harmonic oscillator length $a_r=\sqrt{\hbar/m\omega_r}$
for the unit of length and $\hbar \omega_r$ for the unit of energy.
For definiteness, we fix the value of the
dimensionless density--density coupling constant to $\tilde{g}=g N /\hbar
\omega_r a_r^2=10^3$, which corresponds roughly to $1.5 \times 10^4$
${}^{87}{\rm Rb}$ ${\rm atoms}/\mu {\rm m}$ in a trap with $a_r=1\,
\mu {\rm m}$.

One particular class of stationary solutions is illustrated in
Fig.~\ref{FIG2} with $g_{\mathrm d}=0.10\, g$. On the left-hand side,
Figs.~\ref{FIG2}(a)--\ref{FIG2}(c), the state is shown in the
axially homogeneous case, $\kappa=0$, and on the right-hand side,
Figs.~\ref{FIG2}(d)--\ref{FIG2}(f), with finite helical pitch,
$\kappa=0.25\, {\rm rad}/a_r$. The left column in each
subfigure shows the amplitudes of the order-parameter components
$|\psi_1|$, $|\psi_0|$, and $|\psi_{-1}|$, from up to down. The complex
phases of the corresponding components are given in the right
column in each subfigure. For both states, the magnetization is
pointing predominantly in the positive $z$ direction in the upper
half $(y>0)$ and in the negative $z$ direction in the lower half
$(y<0)$ of the $xy$ plane.

The complex-phase plots in Figs.~\ref{FIG2}(a)--\ref{FIG2}(c) reveal
that the system is hosting two integer-spin vortices with phase
windings $(0,1,2)$ (upper half) and $(-2,-1,0)$ (lower half) in the
components ($\psi_1$, $\psi_0$, $\psi_{-1}$), respectively. The
ferromagnetic cores of the vortices are deformed into elliptic
shapes, which is also indicated by the separation of the phase
singularities in the $\psi_{\pm 1}$ components, cf.
Ref.~\cite{Simula2010}. The two vortices carry both spin- and mass
currents: the spin currents flow in the same direction whereas the
mass currents flow in the opposite directions, canceling the total
mass current in the state.

For finite helical wave vector, $\kappa>0$, the vortices move closer
to the surface of the cloud, as illustrated in
Figs.~\ref{FIG2}(d)--\ref{FIG2}(f). Not only do the cores of the
vortex pair separate farther apart for larger $\kappa$, but also the
separation between the singularities in the $\psi_{\pm 1}$ components
increases. However, one should note that for larger $\kappa$, the
vortex lines are tilted steeper with respect to the axial direction.
Due to the structure of the topological defects hosted by the order
parameter, we refer to this state as the Mermin--Ho vortex (MHV) helix.
This type of solution is found to exist in the
whole stability range of the system, $0 \le g_{\mathrm d}/g \lesssim 0.24$.

\begin{figure}[!t]
\includegraphics[width=240pt]{FIG2.eps}
\caption{\label{FIG2} (Color online) Illustration of the MHV
helix state for $g_{\mathrm d}=0.10\, g$ with $\kappa=0$ in (a)--(c)
and $\kappa=0.25\, {\rm rad}/a_r$ in (d)--(f). For both states, the
amplitudes of the order-parameter components ($\psi_1$,
$\psi_0$, $\psi_{-1}$) are given in the left column and
the complex phases in the right column. This class of
states hosts a pair of MHVs with aligned spin
currents and opposite mass flows. For each state separately, the
color-map range in the left column is $\big[0, {\rm max}
\left\{n({\bf r})\right\} \big]$ and in right column $\big[0,
2\pi\big]$. The field of view in each panel is $[14\, a_r \times
14\, a_r]$.}
\end{figure}

The MHV helix is the ground-state texture within the symmetry requirements of
Eqs.~(\ref{Symmetry_Conditions}) for all parameter values considered
in this work. However, there exists an interesting class of excited
states which is illustrated in Fig.~\ref{FIG3} for $g_{\mathrm d}=0.10\, g$
with $\kappa=0$ and $\kappa=1.0\, {\rm rad}/a_r$ on the left- and
right-hand side, respectively. The complex-phase plots of
Figs.~\ref{FIG3}(a) and \ref{FIG3}(c) reveal that the state hosts
two spin vortices with the same phase windings $(-1,0,1)$ in the components
$(\psi_{1}, \psi_{0}, \psi_{-1})$, respectively. These
defects are spin vortices with polar (nonmagnetized) core regions.
Such vortices carry a spin current, whereas the mass current about the
vortex line vanishes. The polar vortex cores within the
ferromagnetic cloud are energetically analogous to air bubbles in
water, hence increasing the total energy of the system through
spin--spin interactions. The third singularity visible in
Figs.~\ref{FIG3}(a) and \ref{FIG3}(c) shows that the state hosts
also a pair of fractional half-quantum vortices close to the
surface of the cloud. However, these defects have very little effect
on the texture because the related phase gradients lie in regions of
nearly vanishing component amplitude $|\psi_{\pm 1}|$ for all values of
$\kappa$.

\begin{figure}[!t]
\includegraphics[width=240pt]{FIG3.eps}
\caption{\label{FIG3} (Color online) Illustration of the PCV
helix state for $g_{\mathrm d}=0.10\, g$ with $\kappa=0$ in (a)--(c)
and $\kappa=1.0\, {\rm rad}/a_r$ in (d)--(f). For both states, the
amplitudes of the order-parameter components ($\psi_1$,
$\psi_0$, $\psi_{-1}$) are given in the left column and
the complex phases in the right column. This class of
states hosts a pair of PCVs with aligned spin
currents and vanishing mass flows. For each state separately, the
color-map range in the left column is $\big[0, {\rm max}
\left\{n({\bf r})\right\} \big]$ and in right column $\big[0,
2\pi\big]$. The field of view in each panel is $[14\, a_r \times
14\, a_r]$.}
\end{figure}

Solutions with similar vortex structures also exist for finite
helical wave vectors, as depicted in
Figs.~\ref{FIG3}(d)--\ref{FIG3}(f) for $\kappa=1.0\, {\rm rad}/a_r$.
Thus, we refer to this state as the polar-core vortex (PCV) helix. This
class of stationary states is found within the dipolar stability
range of the system, except for weak dipolar interaction strengths,
$g_{\mathrm d}/g \lesssim 10^{-3}$, likely due to numerical difficulties.

If one considers the vortex pairs in Figs.~\ref{FIG2} and \ref{FIG3}
as single entities, the total phase windings in both states are
$(-2,0,2)$ in the order-parameter components $(\psi_1, \psi_0,
\psi_{-1})$, respectively. By tracing the order parameter about the
$z$-axis along a path close to the surface of the condensate, the expectation
value of the spin rotates a total angle of $4\pi$ about the axial
direction. This suggests that both states represent doubly quantized
spin vortices that have split into different kinds of singly quantized
spin defects. The splitting of doubly quantized mass vortices in scalar
condensates has been studied previously both experimentally and
numerically~\cite{Shin2004,Mottonen2003,Huhtamaki2006,Mateo2006}.

The spin textures of the states depicted in Figs.~\ref{FIG2} and
\ref{FIG3} are shown in Fig.~\ref{FIG4}. Figures \ref{FIG4}(a) and
\ref{FIG4}(b) illustrate the planar magnetization ${\bf M}'(r,\varphi)$ for
the MHV helix and
Figs.~\ref{FIG4}(c) and \ref{FIG4}(d) for the PCV
helix. The projection of magnetization on the $xy$ plane, ${\bf
M}'_{xy}(r,\varphi)=M'_x(r,\varphi) {\bf e}_x + M'_y(r,\varphi) {\bf e}_y$, is depicted
by the cones, with the length of each cone being proportional to the
local magnitude $\sqrt{M'^2_x(r,\varphi)+M'^2_y(r,\varphi)}$. The axial
magnetization, $M'_z(r,\varphi)$, is shown by color. At the center of the
system in the axially homogeneous textures ($\kappa=0$) in
Figs.~\ref{FIG4}(a) and \ref{FIG4}(c), the magnetization is pointing
perpendicular to $\nabla M'_z(r,\varphi)$ in the MHV and parallel to $\nabla M'_z(r,\varphi)$ in the
PCV state. For the MHV helix, this relative orientation is
maintained for all values of the helical wave vector $\kappa$.
However, for the PCV helix, the relative orientation
twists continuously with increasing $\kappa$ and finally locks into
the symmetric configuration shown in Fig.~\ref{FIG4}(d) for $\kappa
\sim 1\, {\rm rad}/a_r$.

\begin{figure}[!t]
\includegraphics[width=245pt]{FIG4.eps}
\caption{\label{FIG4} (Color online) Illustration of the spin
textures for the states shown in Figs.~\ref{FIG2} and \ref{FIG3}.
Panels (a) and (b) correspond to the MHV helix state with
$\kappa=0$ and $\kappa=0.25\, {\rm rad}/a_r$, respectively. Panels (c) and
(d) refer to the PCV helix state with $\kappa=0$ and
$\kappa=1.0\, {\rm rad}/a_r$, respectively. The cones show the
projection of magnetization onto the $xy$ plane and the color
indicates the local value of axial magnetization $M'_z(x,y)$.
The field of view in each panel is $[13\, a_r \times 13\, a_r]$.}
\end{figure}

In Fig.~\ref{FIG5}, the total energy per particle, $E_{\rm tot}$, is shown as a function of the
helical wave vector $\kappa$ for the
MHV and the PCV helices. As previously, the coupling constants have the
values $\tilde{g}=1000$ and $g_{\mathrm d}/g=0.10$ for both states. The MHV helix
is the ground-state texture, within the constraints of
Eq.~(\ref{Symmetry_Conditions}), for all values of $\kappa$. Whereas
$E_{\rm tot}$ for PCV helix is minimized for the
axially homogeneous state, $\kappa=0$, the MHV helix
is stabilized for a finite wave vector $\kappa_{\rm min}$. This
minimum persists in the total energy for all positive interaction
strengths $g$ and $g_{\mathrm d}$ considered. However, the relative
magnitude of the dip, $\left[E_{\rm tot}(0)-E_{\rm
tot}(\kappa_{\rm min})\right]/E_{\rm tot}(0)$, decreases for weaker dipolar
interaction strengths. The curve $E_{\rm tot}(\kappa)$ is
roughly a shifted parabola, closely resembling the result
that long-period helical structures in MnSi and FeGe can become
stable due to ferromagnetic Dzyaloshinskii
instability~\cite{Bak1980}.

The energy-minimizing helical wave vector $\kappa_{\rm min}$ is
shown in Fig.~\ref{FIG5}(b) as a function of the dipolar interaction
strength $g_{\mathrm d}/g$ for two values of the dimensionless
coupling constant, $\tilde{g}=1000$ (solid curve) and $\tilde{g}=200$ (dashed curve).
In the absence of dipolar interactions, the total energy is
minimized by the axially homogeneous texture for which the kinetic
energy is minimal. The value of $\kappa_{\rm min}$ increases rapidly
as a function of $g_{\mathrm d}/g$ because the dipolar interactions,
which overwhelm the kinetic energy already for $g_{\mathrm d}/g
\approx 0.01$, favor a finite helical pitch. In general,
$\kappa_{\rm min}$ increases for stronger dipolar couplings. However,
the pitch of the energy-minimizing helix, $p_{\rm
min}=2\pi/\kappa_{\rm min}$, is not determined directly by the
dipole--dipole coherence length $\xi_{\rm dd}=\hbar/\sqrt{2m
g_{\mathrm d} n(0)}$, where $n(0)$ is the particle density at the
trap center. In fact, the coherence length $\xi_{\rm dd}$ shrinks
for larger particle numbers, whereas $p_{\rm min}$ typically
increases. This is likely due to the overall expansion of the condensate.

\begin{figure}[!t]
\includegraphics[width=250pt]{FIG5.eps}
\caption{\label{FIG5} (Color online) (a) Total energy per particle,
$E_{\rm tot}$, as a function of the helical wave vector $\kappa$ for
the MHV helix (solid curve) and the PCV helix (dashed
curve) with $\tilde{g}=1000$ and $g_{\mathrm d}/g=0.10$. The total
energy of the MHV helix exhibits a local minimum with
a finite wave vector $\kappa_{\rm min}$. (b) The energy-minimizing
wave vector $\kappa_{\rm min}$ for the MHV helix as a
function of the dipolar interaction strength $g_{\mathrm d}/g$ for
two values of the density--density coupling constant, $\tilde{g}=1000$ (solid
curve) and $\tilde{g}=200$ (dashed curve).}
\end{figure}

As a measure of how large the axial magnetization $M_z({\bf r})$ is on
average, we define the integrated axial magnetization per unit length as
$\mathcal{M}_z=\int d^2r |M'_z(r,\varphi)|$. Similarly, the integrated
transversal magnetization is given by $\mathcal{M}_r=\int d^2r
\sqrt{M'^2_x(r,\varphi)+ M'^2_y(r,\varphi)}$. These quantities are plotted in
Figs.~\ref{FIG6}(a)--\ref{FIG6}(d) as functions of the helical
wave vector $\kappa$ for $\tilde{g}=1000$ and four different values of the
dipolar coupling constant $g_{\mathrm d}/g$. The solid curve
corresponds to $\mathcal{M}_r$ and the dashed curve to $\mathcal{M}_z$.
As shown in Fig.~\ref{FIG6}(a), the axial magnetization
vanishes when $g_{\mathrm d}/g=0$. For $\kappa \gtrsim
0.7\, {\rm rad}/a_r$, the system becomes nonmagnetized because the
additional kinetic energy due to finite $\kappa$ for magnetized states exceeds
the energy gain from the spin--spin interaction, cf.
Ref.~\cite{Cherng2008}. Already for weak dipolar coupling,
$g_{\mathrm d}/g=10^{-3}$, which is roughly the value for
${}^{87}{\rm Rb}$, the axial magnetization becomes significant, as
shown in Fig.~\ref{FIG6}(b). Also, for finite dipolar coupling, the
value of $\kappa$ at which the system enters a polar state is
increased. The axial component of magnetization becomes larger for increasing $g_\mathrm{d}$, as illustrated in
Figs.~\ref{FIG6}(c) and \ref{FIG6}(d).

\begin{figure}[!t]
\includegraphics[width=250pt]{FIG6.eps}
\caption{\label{FIG6} (Color online) Integrated transversal $\mathcal{M}_r$ (solid curve) and axial $\mathcal{M}_z$ (dashed curve)
magnetization per unit length as a function of the helical wave vector
$\kappa$ for the MHV helix. The strength of the
dipolar coupling $g_{\mathrm d}/g$ is indicated in each
panel. The overall magnetization is reduced for larger values of
$\kappa$, and the axial magnetization tends to increase with $g_\mathrm{d}/g$.}
\end{figure}

\section{Discussion}
In summary, we have studied helical spin textures in a spin-1
BEC subject to long-range dipolar
interactions. The axially elongated system was assumed to be confined radially by a
cylindrically symmetric harmonic potential. By using a helical
Ansatz, we reduced the zero-temperature GP equation
to a two-dimensional problem, in which the helical wave vector $\kappa$
appeared as a parameter. This allowed us to
investigate states with variable values of $\kappa$. We found
two classes of helical solutions which we refer to as the Mermin--Ho
and the polar-core vortex helices, according to the structure of their
topological defects. Whereas the total energy
of the polar-core vortex helix is minimized for the axially
homogeneous case, $\kappa=0$, the Mermin--Ho vortex helix is
stabilized for a finite pitch.

The helical spin textures studied in this work are naturally most
transparent in elongated systems, such as cigar-shaped
BECs. One difficulty in observing them
as stable configurations in a condensate subject to weak dipolar interactions
is that the spins tend to align predominantly
parallel to the weak axis of the trap due to the head-to-tail
attraction of the dipoles~\cite{Yi2006,Takahashi2007,Huhtamaki2010}.
However, we have performed additional three-dimensional simulations indicating that this problem can be overcome
by placing the cigar-shaped system in a one-dimensional optical
lattice potential. A strong enough lattice deforms the elongated condensate into a series of
oblate clouds. Within each cloud, the preferred direction of
magnetization lies perpendicular to the weak axis of the trap, and
the relative orientation of magnetization between neighboring
systems is determined by the long-range part of the interaction
potential, i.e., by dipolar forces.

In the three-dimensional simulations, we calculated the total energy of an axially spin-polarized state
and a helical spin texture subject to the symmetry constraint in
Eq.~(\ref{Symmetry_Conditions}) as a function of the strength of the
optical lattice potential. We chose the aspect ratio of the
confining harmonic potential such that $\omega_z/\omega_r = 0.10$,
where $\omega_z$ and $\omega_r$ are the axial and radial trapping
frequencies, respectively. The values of the dimensionless coupling
constants were fixed to $\tilde{g}=5 \times 10^4$, $g_{\mathrm
s}/g=-10^{-2}$, and $g_{\mathrm d}/g=10^{-3}$, corresponding to $8 \times 10^5$ ${}^{87}{\rm Rb}$ atoms in a trap with $\omega_r=2\pi \times 100\, {\rm Hz}$. The optical
lattice potential was of the form $V_{\rm opt}=E_{\mathrm s} \sin^2
\left( k_{\mathrm s} z \right)$, where the constant $E_{\mathrm s}$
is the strength of the lattice potential, and the wave vector
$k_{\mathrm s}$ determining the distance between neighboring clouds
was fixed to $k_{\mathrm s}=\pi/10\, a_r$. The energy difference
between the helical and polarized states was found to decrease monotonously as a
function of $E_{\mathrm s}$, reaching degeneracy at
$E_{\mathrm s} \approx 17\, \hbar\omega_r$.

\begin{acknowledgments}
We acknowledge financial support from Japan Society for the Promotion of Science
(JSPS), the Emil Aaltonen Foundation, the V\"ais\"al\"a Foundation,
and Finnish Academy of Science and Letters.
V.~Pietil\"a, T.~P.~Simula, T.~Mizushima, and K.~Machida are appreciated for
useful comments and discussion.
\end{acknowledgments}


\section{Introduction}
When two projections $p$ and $q$ in $B(H)$, whose difference is compact, are given, an integer $[p:q]$ is defined as the Fredholm index of $v^*w$ where $v,w$ are isometries on $H$ with $vv^*=p$ and $ww^*=q$. This number is called the essential codimension because it gives the codimension of $p$ in $q$ if $p\leq q$ \cite{BDF}. A modern interpretation of this essential codimension is provided using the Kasparov group $\KK(\mathbb{C},\mathbb{C})$. Indeed, a $*$-homomorphism from $\mathbb{C}$ to $B(H)$ is determined by the image of $1$ which is a projection. Thus we can associate to the essential codimension  a Cuntz pair.
An important result of the essential codimension is the following: $[p:q]=0$ if and only if there is a unitary $u$ of the form identity + compact such that $upu^*=q$. Motivated by this result, Dadarlat and Eilers defined a new equivalence relation on $\KK$-group \cite{DE}. When  $\pi,\sigma: A \to \mathcal{L}(E)$ are two representations, with $E$ is a Hilbert $B$-module, we say $\pi$ and $\sigma$ are properly asymptotically unitarily equivalent and write $\pi \approxeq \sigma$ if there is a continuous path of unitaries $u:[0,\infty) \to \mathcal{U}(\mathcal{K}(E)+\mathbb{C}1_E)$, $u=(u_t)_{t\in[0,\infty)}$, such that
\begin{itemize}
\item $lim_{t\to \infty}\| \sigma(a)-u_t \pi(a)u_t^*\|=0$ for all $a \in A$,
\item $ \sigma(a)-u_t \pi(a)u_t^* \in \mathcal{K}(E)$ for all $t\in [0,\infty )$, and $a\in A$.
\end{itemize}

Note that the word `proper' reflects the fact that implementing unitaries are of the form `identity+compact'. The main result of them is \cite[Theorem 3.8]{DE} which asserts that if $\phi, \psi: A \to M(B\otimes K(H))$ is a Cuntz pair of representations, then the class $[\phi,\psi]$ vanishes in $\KK(A,B)$ if and only if there is another representation $\gamma:A \to M(B\otimes K(H))$ such that $\phi\oplus\gamma \approxeq \psi\oplus \gamma$. When $B=\mathbb{C}$, which corresponds to $K$-homology, the result is improved as a non-stable version. In fact, if $(\phi,\psi)$ is a Cuntz pair of faithful, non-degenerate representations from $A$ to $B(H)$ such that both images do not contain any non-trivial compact operator, then the cycle $[\phi,\psi]=0 $ in $\KK(A,\mathbb{C})$ if and only if $\phi \approxeq \psi$ \cite[Theorem 3.12]{DE}. This fits nicely with the above aspect of the essential codimension. An abstract version of this is proved given a Cuntz pair of absorbing representations (See Theorem \ref{T:DE}). Thus the proper asymptotic unitary equivalence must be the right notion and tool for further developments of the non-stable K-theory. Our intrinsic interest lies in when this non-stable version of proper asymptotic unitary equivalence happens as shown in $K$-homology case. We show a similar result for $K$-theory. In fact, we prove that if $(\phi,\psi) $ is a Cuntz pair of faithful representations from $\mathbb{C} \to M(B\otimes K)$ whose images are not in $B\otimes K$, then $[\phi,\psi]=0$ in $K(B)$ if and only if $\phi \approxeq \psi$ provided that $B$ is non-unital, separable, purely infinite simple $C\sp*$-algebra such that $M(B)$ has real rank zero (See Theorem \ref{T:BDF}).

Besides our intrinsic interest, Theorem \ref{T:BDF} was motivated by the projection lifting problem from the corona algebra to the multiplier algebra of a $C\sp*$-algebra of the form $C(X)\otimes B$. To lift a projection from a quotient algebra to a projection  has been a fundamental question related to K-theory (See \cite{Da}). We show that  a projection in the corona algebra is `locally' liftable to a projection in the multilpler algebra  but not `globally' in general. In other words, it can be represented by finitely many  projection valued functions so that  its discontinuities are described in terms of  Cuntz pairs. They give rise to $K$-theoretical obstructions. We show that these discontinuities  can be resolved if corresponding $K$-theoretical terms are vanishing. In this process, the crucial point of proper asymptotic unitary equivalence is exploited as a key step (See Theorem \ref{T:lifting}).

\section{Proper asymptotic unitary equivalence}
Let $E$ be a (right) Hilbert $B$-module. We denote $\mathcal{L}(E,F)$ by the $C\sp{*}$-algebra of adjointable, bounded  operators from $E$ to $F$. The ideal of `compact' operators from $E$ to $F$ is denoted by $\mathcal{K}(E,F)$. When $E=F$, we write $\mathcal{L}(E)$ and $\mathcal{K}(E)$ instead of  $\mathcal{L}(E,E)$ and $\mathcal{K}(E,E)$. Throughout the paper, $A$ is a separable $C\sp{*}$-algebra, and all Hilbert modules are assumed to be countably generated over a separable $C\sp*$-algebra. We use the term representation for a $*$-homomorphism from $A$ to $\mathcal{L}(E)$. We let $H_B$ be the standard Hilbert module over $B$ which is $H\otimes B $ where $H$ is a separable infinite dimensional Hilbert space. We denote $M(B)$ by the multiplier algebra of $B$. It is well-known that $\mathcal{L}(H_B)=M(B\otimes K)$ and $\mathcal{K}(H_B)=B\otimes K$ where $K$ is the $C\sp{*}$-algebra of the compact operators on $H$ \cite{Kas80}.
\begin{defn}\cite[Definition 2.1]{DE}
Let $\pi,\sigma$ be two representations from $A$ to $E$ and $F$ respectively. We say $\pi$ and $\sigma$ are approximately unitarily equivalent and write $\pi \sim \sigma$, if there exists a sequence of unitaries $u_n \in \mathcal{L}(E,F)$ such that for any $a\in A$
\begin{itemize}
\item[(i)] $lim_{n\to \infty}\| \sigma(a)-u_n \pi(a)u_n^*\|=0$,\\
\item[(ii)] $ \sigma(a)-u_n \pi(a)u_n^* \in \mathcal{K}(F)$ for all $n$.
\end{itemize}
\end{defn}
\begin{defn}\cite[Definition 2.5]{DE}
A representation $\pi:A \to \mathcal{L}(E)$ is called absorbing if $\pi \oplus \sigma \sim \pi$ for any representation $\sigma :A \to \mathcal{L}(F)$.
\end{defn}
  We say that $\pi$ and $\sigma$ are asymptotically unitarily equivalent, and write $\pi \underset{\text{asym}}{\sim}\sigma$ if there is a unitary valued norm continuous map $u:\left[0,\infty \right) \to \mathcal{L}(E,F)$ such that $t \to \sigma(a)- u_t\pi(a)u_t^* $ lies in $C_0([0,\infty))\otimes \mathcal{K}(E)$ for any $a \in A$, or if
\begin{itemize}
\item[(i)]$lim_{t\to \infty}\| \sigma(a)-u_t \pi(a)u_t^*\|=0$,\\
\item[(ii)] $ \sigma(a)-u_t \pi(a)u_t^* \in \mathcal{K}(F)$ for all $t\in [0,\infty)$.
\end{itemize}
If $\pi:A \to \mathcal{L}(E)$ is a representation, we define $\pi^{(\infty)}:A \to \mathcal{L}(E^{(\infty)})$ by as $\pi^{(\infty)}=\pi\oplus \pi\oplus \cdots$ where $E^{(\infty)}=E\oplus E\oplus \cdots$.
\begin{lem}\label{L:absorbing}
 Let $\psi$ be an absorbing representation, and $\phi$ be a representation, of a separable $C\sp{*}$-algebra $A$ on the standard Hilbert $B$-module $H_B$. Then there exists a sequence of isometries $\{v_n\} \subset \mathcal{L}(H_B^{(\infty)},H_B)$ such that for each $a \in A$
\begin{align*}
   &v_n \phi^{(\infty)}(a) -\psi(a) v_n \in \mathcal{K}(H_B^{(\infty)}, H_B), \\
   &\| v_n \phi^{(\infty)}(a) -\psi(a) v_n  \| \to 0 \quad \text{as $n \to \infty$},\\
   &v_j^* v_i=0 \quad \text{for $i \neq j$}.
\end{align*}
\end{lem}
\begin{proof}

Let $S_i, \, i=1,2,3 \cdots $, be a sequence of isometries of $\mc{L}(H_B)$ such that $S_i^*S_j=0, \,i\neq j$, and $\sum_iS_iS_i^*=1$ in the strict topology. Let $\phi_{\infty}(a)=\sum_i S_i \phi(a) S_i^*$. Since $\psi$ is absorbing, there is a unitary $U \in \mathcal{L}(H_B, H_B)$ such that
\begin{equation}\label{E:(1)}
   U^* \psi(a)U- \phi_{\infty}(a)  \in \mc{K}(H_B) \quad a\in A.
\end{equation}
Define $T:H_B^{(\infty)} \to H_B$ by $T=(S_1,S_2,\cdots )$. Then  $$\phi_{\infty}(a)=T \phi^{(\infty)}(a)T^*. $$
Thus equation (\ref{E:(1)}) is rewritten as
\begin{equation}\label{E:(2)}
  T^*U^* \psi(a)UT- \phi^{(\infty)}(a)  \in \mc{K}(H_B^{(\infty)}) \quad a\in A.
 \end{equation}
 If we identify $\phi^{(\infty)}$ as $(\phi^{(\infty)})^{(\infty)}$, there is a partition $N_i$, $i=1,2,3,\cdots$, of $\mathbb{N}$ so that we generate a sequence of isometries $v_i \in \mc{L}(H_B^{(\infty)},H_B)$ from $UT=(US_1,US_2,\cdots,)$. More concretely, if we let $\nu_i: N_i \to \mathbb{N}$ be bijections, we can define
 $v_i=(US_{\nu_i^{-1}(1)},US_{\nu_i^{-1}(2)},\cdots)$. It is easily checked that $v_iv_j^*=0$ for $i\neq j$.
 Equation (\ref{E:(2)}) implies that
 \begin{align*}
 & v_i^*\psi(a) v_i -\phi^{(\infty)}(a) \in \mc{K}(H_B^{(\infty)}), \\
 &\left\| v_i^*\psi(a) v_i -\phi^{(\infty)}(a) \right\| \to 0 \quad \text{as}\, i \to \infty.
 \end{align*}
 Finally, our claim follows from
 \[
 \begin{split}
 &(v_n \phi^{(\infty)}(a) -\psi(a) v_n)^*(v_n \phi^{(\infty)}(a) -\psi(a) v_n)=\phi^{(\infty)}(a^*)(\phi^{(\infty)}(a)-v_n^*\psi(a)v_n)\\
 &+(\phi^{(\infty)}(a^*)-v_n^*\psi(a) v_n)\phi^{(\infty)}(a)-(\phi^{(\infty)}(a^*a)-v_n^*\psi(a^*a)v_n).
 \end{split}
 \]
\end{proof}

\begin{lem}\cite[Lemma 2.6]{DE}\label{L:asym}
Let $\pi:A \to \mathcal{L}(E)$ and $\sigma:A \to \mathcal{L}(F)$ be two representations. Suppose that there is a sequence of
isometries $v_i:F^{(\infty)} \to E$ such that for  $a\in A$
$$ v_i\sigma ^{(\infty)}(a)-\pi(a)v_i \in \mathcal{K}(F^{(\infty)}, E), \quad \lim_{i \to \infty}\|v_i\sigma ^{(\infty)}(a)-\pi(a)v_i\|\to 0, $$
and $v_j^*v_i=0$ for $i\neq j$. Then $\pi \oplus\sigma \underset{\text{asym}}{\sim}\pi$.
\end{lem}
 We say $\phi:A \to B(H)$ is admissible if $\phi$ is faithful, non-degenerate, and $\phi(A)\cap K = \{0 \}$. The main result in \cite {Voi}
states that any pair of admissible representations $\phi$ and $\psi$ satisfies that $\phi \sim \psi$. Dadarlat and Eilers proved a much stronger version which states that any pair of admissible representations $\phi$ and $\psi$ satisfies  $\phi \underset{\text{asym}}{\sim} \psi$ \cite[Theorem 3.11]{DE}. Since the admissible representation is absorbing, the following result is the appropriate generalization of Voiculescu's result.
\begin{thm}\label{T:Voiculescu}
If two representations $\psi$, $\phi$ of a separable $C\sp{*}$-algebra $A$ on the standard Hilbert $B$-module $H_B$ are absorbing,
then we have $\phi \underset{\text{asym}}{\sim} \psi$.
\end{thm}
\begin{proof}
By Lemma \ref{L:absorbing} and Lemma \ref{L:asym}, we have $\psi\oplus\phi  \underset{\text{asym}}{\sim} \psi$, and the proof is complete by symmetry.
\end{proof}

\begin{defn}
Let $\phi$ be a representation from $A$ to $M(B\otimes K)$. Then we define a $C\sp{*}$-algebra by
$$D_{\phi}(A,B)=\{x\in M(B\otimes K)\mid x\phi(a)-\phi(a)x \in B \otimes K, \quad a \in A \}.$$
\end{defn}

\begin{lem}
If $M(B\otimes K)$ has real rank zero, then $D_{\phi}(\mathbb{C},B)$ has real rank zero for any representation $\phi:\mathbb{C} \to M(B\otimes K)$.
\end{lem}
\begin{proof}
The proof of the lemma is essentially based on the argument due to Brown and Pedersen \cite{BP}.

Note that any representation $\phi: \mathbb{C} \to M(B\otimes K)$ is determined by $\phi(1)$, which is a projection in $M(B\otimes K)$. Say $\phi(1)=p$. Then we see that $D_{\phi}(\mathbb{C},B)=\{x\in M(B\otimes K)\mid x p- p x \in B \otimes K \}$.

To show $D_{\phi}(\mathbb{C},B)$ has real rank zero, it is enough to show any self-adjoint element in $D_{\phi}(\mathbb{C},B)$ is approximated by a self-adjoint, invertible element. Let $x$ be a self-adjoint element. Using the obvious matrix notation
\[x=\left(
                        \begin{matrix}
             a & c \\
             c^* & b
                        \end{matrix}
                      \right),\]
$x p- p x \in B\otimes K$ implies that $c$ is `compact', i.e., it is in $B\otimes K$.
Since $M(B\otimes K)$ has real rank zero, $pM(B\otimes K)p$ and $(1-p)M(B\otimes K)(1-p)$ have real rank zero. Given $\epsilon>0$ we can find $b_0$ invertible in $(1-p)M(B\otimes K)(1-p)$ with $b_0=b_0^*$ and $\|b-b_0 \|<\epsilon$. Then considering $a-cb_0^{-1}c^*$, we can find $a_0$ in $pM(B\otimes K)p$ with $a_0=a_0^*$ and $\| a-a_0\|<\epsilon$, such that $a_0-cb_0^{-1}c^*$ is invertible in $pM(B\otimes K)p$. Then
$\left(\begin{matrix}
             p & cb_0^{-1} \\
             0 & 1-p
                        \end{matrix}
                      \right) $, $\left(
                        \begin{matrix}
             p & 0 \\
             b_0^{-1}c^* & 1-p
                        \end{matrix}
                      \right)$ are in $D_{\phi}(\mathbb{C},B)$ since $ cb_0^{-1}$ is `compact'. Thus
\[x_0=\left(
                        \begin{matrix}
             a_0 & c \\
             c^* & b_0
                        \end{matrix}
                      \right)= \left(
                        \begin{matrix}
             p & cb_0^{-1} \\
             0 & 1-p
                        \end{matrix}
                      \right)\left(
                        \begin{matrix}
             a_0-cb_0^{-1}c^* & 0  \\
             0 & b_0
                        \end{matrix}
                      \right)\left(
                        \begin{matrix}
             p & 0 \\
             b_0^{-1}c^* & 1-p
                        \end{matrix}
                      \right)\]
is invertible in $D_{\phi}(\mathbb{C},B)$. Evidently $\|x-x_0 \| <\epsilon$, so we are done.
\end{proof}

Let us recall the definition of Kasparov group $\KK(A,B)$. We refer the reader to \cite{Kas81} for the general introduction of the subject. A $\KK$-cycle is a triple $(\phi_0,\phi_1,u)$, where $\phi_i:A \to \mathcal{L}(E_i)$ are representations and $u \in \mathcal{L}(E_0,E_1)$ satisfies that
\begin{itemize}
\item[(i)] $u\phi_0(a)-\phi_1(a)u \in \mathcal{K}(E_0,E_1)$,
\item[(ii)]$\phi_0(a)(u^*u-1)\in \mathcal{K}(E_0)$, $\phi_1(a)(uu^*-1)\in \mathcal{K}(E_1)$.
\end{itemize}
The set of all $KK$-cycles will be denoted by $\mathbb{E}(A,B)$.  A cycle is degenerate if
\[ u\phi_0(a)-\phi_1(a)u=0, \quad\phi_0(a)(u^*u-1)=0,\quad \phi_1(a)(uu^*-1)=0.\]
An operator homotopy through $KK$-cycles is a homotopy $(\phi_0,\phi_1,u_t)$, where the map $t \to u_t$ is norm continuous.
The equivalence relation $\underset{\text{oh}}{\sim}$ is generated by operator homotopy and addition of degenerate cycles up to unitary equivalence.
Then $\KK(A,B)$ is defined as the quotient of $\mathbb{E}(A,B)$ by $\underset{\text{oh}}{\sim}$. When we consider  non-trivially graded $C\sp{*}$-algebras, we define a triple $(E,\phi,F)$, where $\phi:A \to \mathcal{L}(E)$ is a graded representation, and $F\in \mathcal{L}(E)$ is of odd degree such that $F\phi(a)-\phi(a)F$, $(F^2-1)\phi(a)$, and $(F-F^*)\phi(a)$ are all in $\mathcal{K}(E)$ and call it a Kasparov $(A,B)$-module. Other definitions like degenerate cycle and operator homotopy are defined in  similar ways.
Let $v$ be a unitary in $M_n(D_{\phi}(A,B))$. Define $\phi^n:A \to \mathcal{L}_B(B^n)$ by $\phi^n(a)(b_1,b_2,\cdots,b_n)=(\phi(a)b_1,\phi(a)b_2,\cdots,\phi(a)b_n)$. Let $B^n\oplus B^n$ be graded by $(x,y)\mapsto (x,-y)$. Then
\[\left( B^n\oplus B^n,\left( \begin{matrix}
                           \phi^n & 0\\
                           0  & \phi^n \end{matrix}\right), \left( \begin{matrix}
                           0 & v\\
                           v^*  & 0 \end{matrix}\right)      \right) \]
 is a Kasparov $(A,B)$-module. The class of this module depends only on the class of $v$ in $K_1(D_{\phi}(A,B))$ so that the construction gives rise to a group homomorphism $\Omega:K_1(D_{\phi}(A,B)) \to \KK(A,B)$.
\begin{lem}\label{L:Thomsen}
Let $\phi$ an absorbing representation from $A$ to $\mathcal{L}(H_B)=M(B)$ where $B$ is a stable $C\sp*$-algebra. Then
$\Omega: K_1(D_{\phi}(A, B)) \to \KK(A,B)$ is an isomorphism.
\end{lem}
\begin{proof}
See \cite[Theorem 3.2]{Th}. In fact, Thomsen proved $K_1(\frac{D_{\phi}(A,B)}{(D_{\phi}(A,A,B))})$ is isomorphic to $\KK(A,B)$ via a map $\Theta$ where $D_{\phi}(A,A,B)=\{x\in D_{\phi}(A,B)\mid x\phi(A) \subset B\}$ is the ideal of $D_{\phi}(A,B)$. However, the same proof shows $\Omega$ is an isomorphism. Alternatively we can show that $K_i(D_{\phi}(A,A,B))=0$ for $i=0,1$ by the argument of \cite[Lemma 1.6]{H} with the fact that $K_*(M(B))=0$.
Thus, using the six term exact sequence,  $K_*(D_{\phi}(A,B))$ is isomorphic to $K_*\left(\frac{D_{\phi}(A,B)}{(D_{\phi}(A,A,B))}\right)$. This implies the map $\Omega$ which is the composition with $\Theta$ and $q_1$ is an isomorphism. Here $q_1$ is the induced map between K-groups from the quotient map from $D_{\phi}(A, B)$ onto $\frac{D_{\phi}(A,B)}{(D_{\phi}(A,A,B))}$.
\end{proof}
\begin{defn}\cite[Definition 3.2]{DE}
If $\pi,\sigma:A \to \mathcal{L}(E)$ are representations, we say that $\pi$ and $\sigma$ are properly asymptotically unitarily equivalent  and write $\pi \approxeq \sigma$ if  there is a continuous path of unitaries $u: [0,\infty ) \to \mathcal{U}(\mathcal{K}(E)+\mathbb{C}I_E)$, $u=(u_t)_{t\in [0,\infty )}$ such that for all $a\in A$
\begin{itemize}
\item[(i)] $lim_{t\to \infty}\| \sigma(a)-u_t \pi(a)u_t^*\|=0$,\\
\item[(ii)] $ \sigma(a)-u_t \pi(a)u_t^* \in \mathcal{K}(E)$ for all $t\in [0,\infty ).$
\end{itemize}
\end{defn}
In the above, we introduced the Fredholm picture of $\KK$-group. There is an alternative way to describe the element of $KK$-group. The Cuntz picture is described by a pair of representations $\phi,\psi:A \to \mathcal{L}(H_B)=M(B\otimes K)$ such that $\phi(a)-\psi(a)\in \mathcal{K}(H_B)=B\otimes K$. Such a pair is called a Cuntz pair. They form a set denoted by $\mathbb{E}_h(A,B)$. A homotopy of Cuntz pairs consists of a Cuntz pair $(\Phi,\Psi):A \to M(C([0,1])\otimes (B\otimes K))$. The quotient of $\mathbb{E}_h(A,B)$ by homotopy equivalence  is a group $\KK_h(A,B)$ which is isomorphic to $\KK(A,B)$ via the mapping sending $[\phi,\psi]$ to $[\phi,\psi,1]$.

 Dadarlat and Eilers proved that $[\phi,\psi]=0$ in $\KK_h(A,B)$ if and only if there is a representation $\gamma:A \to M(B\otimes K)=\mathcal{L}(H_B)$ such that $\phi\oplus\gamma \approxeq \psi\oplus\gamma$ \cite[Proposition 3.6]{DE}. The point is that the equivalence is implemented by unitaries of the form  compact $+$ identity. Sometimes, we can have a non-stable equivalence keeping this useful point.

\begin{defn}
 Let  $A$ be a $C\sp*$-algebra. Denote by ${\widetilde A}$ its unitization. We say that $A$ has $K_1$-injectivity if the map from $\mathcal{U}(\widetilde{A})/\mathcal{U}_0(\widetilde{A})$ to $K_1(A)$ is injective where $\mathcal{U}(\widetilde{A})$ is the unitary group and $\mathcal{U}_0(\widetilde{A})$ is the connected component of the identity. We note that H. Lin proved in \cite[Lemma 2.2]{Lin96} that real rank zero implies $K_1$-injectivity.
\end{defn}

\begin{thm}\label{T:DE}
Let $A$ be a separable $C\sp*$-algebra and let $\psi, \phi: A \to H_B$ be a Cuntz pair of absorbing representations. Suppose that the composition of $\phi$ with the natural quotient map $\pi:M(B\otimes K) \to M(B\otimes K)/ B\otimes K$, which will be denoted by $\dot{\phi}$, is faithful. Further, we suppose that $D_{\phi}(A,B)$ satisfies $K_1$-injectivity. If $[\phi,\psi]=0$ in $KK(A,B)$, then $\phi \approxeq \psi$.
\end{thm}
\begin{proof}
The proof of this theorem is almost identical to the one given in \cite[Theorem 3.12]{DE}. We just give the proof to illustrate how our assumptions play the roles.

By Theorem \ref{T:Voiculescu}, we get a continuous family of unitaries $(u_t)_{t\in\left[0,\infty\right)}$ in $M(B\otimes K)$ such that
\begin{equation}\label{E:asymt}
u_t\phi(a)u_t^*- \psi(a) \in C_0(\left[0,\infty \right))\otimes (B\otimes K).
\end{equation}
Note that (\ref{E:asymt}) implies $[\phi,\psi]=[\phi,u_1\phi u_1^*]$ (See \cite[Lemma 3.1]{DE}).
We assume that $[\phi,\psi]=0$ and we conclude that $[\phi,u_1 \phi u_1^*]=0$. Since $(\phi,\phi,u_1^*)$ is unitarily equivalent to $(\phi,u_1\phi u_1^*,1)$, \[[(\phi,\phi,u_1)]=[(\phi,\phi,u_1^*)]=0.\]
Since the isomorphism $\Omega:K_1(D_{\phi}(A,B)) \to \KK(A,B)$  sends $[u_1]$ to $[\phi,\phi,u_1]$  by Lemma \ref{L:Thomsen}, $K_1$-injectivity implies that $u_1$ is homotopic to $1$ in $D_{\phi}(A,B)$. Thus we may assume that $u_0=1$ in (\ref{E:asymt}).

Let $E_{\phi}$ be a $C\sp*$-algebra $\phi(A)+B\otimes K$.
We define $(\alpha_t)_{t\in \left[ 0, \infty \right)}$ in $\Aut_{0}(E_{\phi})$ by $\Ad (u_t)$. Note that $\alpha_0=\id $ and $(\alpha_t)$ is a uniform continuous family of automorphisms. Thus we apply Proposition 2.15 in \cite{DE} and get a continuous family $(v_t)_{\left[ 0, \infty\right)}$ of unitaries in $E_{\phi}$ such that
\begin{equation}\label{E:Aut}
 \lim_{t \to \infty}\| \alpha_{t}(x) -\Ad v_t(x)\|=0
\end{equation}
for any $x \in E_{\phi}$.

Combining (\ref{E:Aut}) with (\ref{E:asymt}), we obtain $(v_t)_{\left[ 0, \infty\right)}$ of unitaries in $E_{\phi}$ such that
\[  \lim_{t \to \infty}\| v_t\phi(a) v_t^*-\psi(a)\|=0 \]
for any $a \in A$. Since $\dot{\phi}$ is faithful, we can replace  $(v_t)_{\left[ 0, \infty\right)}$ by a family of unitaries in
$B\otimes K +\mathbb{C}1$ by the argument shown in the STEP 1 of the proof of Proposition 3.6 in \cite{DE}.
\end{proof}
Recall the definition of the esssential codimension of Brown, Douglas, and Fillmore defined by  two projections $p,q$ in $B(H)$  whose difference is compact as we have defined  in Introduction. Using KK-theory, or K-theory, we generalize this notion as follows, keeping the same notation.
\begin{defn}\label{D:BDF}
Given two projections $p,q \in M(B\otimes K)$ such that $p-q \in B\otimes K$, we consider representations $\phi,\psi $ from $\mathbb{C}$ to $M(B\otimes K)$ such that $\phi(1)=p,\psi(1)=q$. Then $(\phi,\psi)$ is a Cuntz pair so that
we define $[p:q]$ as the class $[\phi,\psi]\in \KK(\mathbb{C},B) \simeq K(B)$.
\end{defn}

\begin{lem}\cite{Lin1}\label{L:Lin}
Let $B$ be a non-unital ($\sigma$-unital) purely infinite simple $C\sp{*}$-algebra. Let $\phi, \psi$ be two monomorphisms from $C(X)$ to $M(B\otimes K)$ where $X$ is a compact metrizable space. If $\dot{\phi},\dot{\psi}$ are still injective, then they are approximately unitarily equivalent.
\end{lem}
The following theorem is a sort of generalization of BDF's result about the essential codimension.
\begin{thm}\label{T:BDF}
Let $B$ be a non-unital ($\sigma$-unital) purely infinite simple $C\sp{*}$-algebra such that $M(B\otimes K)$ has real rank zero.
Suppose two projections $p$ and $q$ in $M(B\otimes K)=\mathcal{L}(H_B)$ such that $p-q \in B\otimes K$ and neither of them are in $B\otimes K$. If $[p,q] \in K_{0}(B)$ vanishes, then there is a unitary $u$ in $ \id + B\otimes K$ such that $upu^*=q $.
\end{thm}
\begin{proof}
Step1: Let $\phi,\psi:\mathbb{C} \to M(B\otimes K)$ be representations from $p$ and $q$ respectively.
Evidently $\phi$ is injective. Moreover, it does not contain any ``compacts'' since $p$ does not belong to $B\otimes K$. Thus $\dot{\phi}$ is faithful.
Recall $\psi_{\infty}$ is defined by $\psi_{\infty}(a)=\sum S_i \psi(a) S_i^*$ where $\{S_i\}$ is a sequence of isometries in $M(B\otimes K)$ such that $S_iS_j^*=0$ for $i\neq j$. Suppose that $\psi_{\infty}(\lambda)=0$ for $\lambda \in \mathbb{C}$. Then $S_i^*\psi_{\infty}(\lambda) S_i=\psi(\lambda)=0$ or  $\lambda q=0$. Thus $\lambda=0$. Similarly, $\dot{\psi}_{\infty}$ is injective.  Then they are approximately unitarily equivalent by  applying Lemma \ref{L:Lin} to $X=\{x_0\}$. Thus we have a unitary $U$ in
$\mathcal{L}(H_B)$ such that
\begin{equation} \label{E:unitaryequivalecne}
U^* \phi(a)U - \psi_{\infty}(a)
\end{equation}
 for $a \in \mathbb{C}$.

Note that to get a sequnce of isometries $\{v_i\} \in \mathcal{L}(H_B^{(\infty)},H_B)$ satisfying the conditions of  Lemma \ref{L:absorbing}, what we needed was the equation (\ref{E:unitaryequivalecne}).  Following the same argument in the proof of  Theorem \ref{T:Voiculescu}, we get $\phi \underset{\text{asym}}{\sim} \psi$. In other words, we have a continuous family of unitaries $(u_t)_{t\in\left[0,\infty\right)}$ in $M(B\otimes K)$ such that
\begin{equation*}\label{E:asym}
u_t\phi(a)u_t^*- \psi(a) \in C_0(\left[0,\infty \right))\otimes (B\otimes K) \quad\text{for any $a$ in A}.
\end{equation*}
Since $D_{\phi}(\mathbb{C},B)$ has real rank zero, it satisfies $K_1$-injectivity. Thus it follows that $\phi \approxeq \psi$ as in the proof of Theorem \ref{T:DE}.

Step2: For large enough $t$, we can take $u_t=u$ of the form `identity + compact' such that $\|upu^*-q\| <1$. For the moment we write $upu^*$ as $p$. Thus $\|p-q\| <1$.
Note that $p-q \in B\otimes K$. Then  $z=pq+(1-p)(1-q)\in 1+ B\otimes K $ is invertible and
 $pz=zq$. If we consider the polar decomposition of $z$ as $z=v|z|$.
 It is easy to check that  $v \in 1+B\otimes K$ and $vpv^{\ast}=q$. Now $w=vu$ is also a unitary of the form `identity + compact' such that
 \[ wpw^*=q.\]

\end{proof}

\section{Application: projection lifting}
In this section, we show an application of proper asymptotic unitary equivalence of two projections. In this application, with an additional real rank zero property, the unitary of the form `identity + compact' plays a crucial role as we shall see.

Let $B$ be a stable $C\sp*$-algebra such that the multiplier algebra $M(B)$ has real rank zero.
Let $X$ be $[0,1], [0,\infty)$, $(-\infty,\infty)$ or
$\mathbb{T}=[0,1]/\{0,1\}$. When $X$ is compact, let $I=C(X)\otimes
   B$ which is the $C\sp{*}$-algebra of (norm continuous) functions from $X$ to $B$. When $X$ is not compact,  let $I=C_0(X)\otimes B$ which is the $C\sp*$-algebra of continuous functions from $X$ to $B$ vanishing at infinity. Then $M(I)$ is given by $C_b(X, M(B)_s)$, which is the set of
   bounded functions from $X$ to $B(H)$, where $M(B)$ is given the strict
   topology. Let $\mathcal{C}(I)=M(I)/I$ be the corona algebra of $I$ and also
   let $\pi:M(I) \to \mathcal{C}(I)$ be the natural quotient map. Then an element
   $\mathbf{f}$ of the corona algebra can be represented as follows:  Consider a
 finite partition of $X$, or $X \smallsetminus \{0,1\}$ when $X=\mathbb{T}$ given by partition points $x_1 < x_2 < \cdots
 < x_n $ all of which are in the interior of $X$ and divide $X$ into
 $n+1$ (closed) subintervals $X_0,X_1,\cdots,X_{n}$. We can take $f_i \in
 C_b(X_i, M(B)_s)$ such that $f_i(x_i) -f_{i-1}(x_i)\in B$
 for $i=1,2,\cdots,n$ and $f_0(x_0)-f_n(x_0) \in B$  where $x_0=0=1$ if $X$ is $\mathbb{T}$.

 \begin{lem}
 The coset in $\mathcal{C}(I)$ represented by
 $(f_0,\cdots,f_n)$  consists of functions $f$ in $M(I)$ such that $f- f_i \in
 C(X_i)\otimes B$ for every $i$ and $f-f_i $ vanishes (in norm) at any
 infinite end point of $X_i$.
 \end{lem}
\begin{proof}
 If $X$ is compact, then we set $x_0=0$, $x_{n+1}=1$. Otherwise, we set $x_0=x_1-1$ when $X$ contains $-\infty$, and $x_{n+1}=x_n +1$ when $X$ contaions $+\infty$.
Then we define a function in $C(X)\otimes B$ by
  \begin{equation*}
m_i(x)=
\begin{cases}
\frac{x-x_{i-1}}{x_i-x_{i-1}}(f_i(x_i) -f_{i-1}(x_i)), &\text{if $x_{i-1}\leq x \leq x_i$}\\
\frac{x-x_{i+1}}{x_i-x_{i+1}}(f_i(x_i) -f_{i-1}(x_i)), &\text{if $x_{i}\leq x \leq x_{i+1}$}\\
0,                               &\text{otherwise}
\end{cases}
\end{equation*}
for each $i=1,\cdots n$. In addition, we set $m_0=m_{n+1}=0$.
Then we define a function $\widetilde{f}$ from $f_i$'s by
\[\widetilde{f}(x)=f_i(x)-m_i(x)/2+m_{i+1}(x)/2\] on each $X_i$.
It follows that $f_i(x_i)-m_i(x_i)/2+m_{i+1}(x_i)/2=f_{i-1}(x_i)-m_{i-1}(x_i)/2+m_{i}(x_i)/2$. Thus $\widetilde{f}$ is well defined.
 The conditions $f-f_i \in C(X_i)\otimes B$ for each $i$ imply that $f-\widetilde{f}$ is  norm continuous function from $X$ to $B$ since
$f|_{X_i}(x_i) -\widetilde{f}|_{X_i}(x_i)= f|_{X_{i-1}}(x_i)-\widetilde{f}|_{X_{i-1}}(x_i)$.

\end{proof}
 Similarly $(f_0,\cdots,f_n)$ and $(g_0,\cdots,g_n)$
 define the same element of $\mathcal{C}(I)$ if and only if $f_i - g_i \in
 C(X_i)\otimes B$ for $i=0,\cdots,n$ if $X$ is compact.  $(f_0,\cdots,f_n)$ and $(g_0,\cdots,g_n)$
 define the same element of $\mathcal{C}(I)$ if and only if  $f_i - g_i \in
 C(X_i)\otimes B$ for $i=0,\cdots,n-1$, $f_n -g_n \in
 C_0([x_n,\infty))\otimes B$  if $X$ is $[0.\infty)$.  $(f_0,\cdots,f_n)$ and $(g_0,\cdots,g_n)$
 define the same element of $\mathcal{C}(I)$ if and only if
  $f_i - g_i \in
 C(X_i)\otimes B$ for $i=1,\cdots,n-1$, $f_n -g_n \in
 C_0([x_n,\infty))\otimes B$, $f_0-g_0 \in
 C_0((-\infty,x_1])\otimes B$ if $X=(-\infty,\infty)$.\\

The following theorem says that any projection in the corona algebra of $C(X)\otimes B$ for some $C\sp{*}$-algebras $B$ is described by a ``locally trivial fiber bundle'' with the fibre $H_{B}$ in the sense of Dixmier and Duady \cite{DixDua}.
\begin{thm}\label{T:locallift}
Let $I$ be $C(X)\otimes B$ or $C_0(X)\otimes B$ where $B$ is a stable $C\sp{*}$-algebra such that $M(B)$  has real rank zero.
Then a projection $\mathbf{f}$ in $M(I)/I$ can be represented by $(f_0,f_1,\cdots, f_n)$ as above where $f_i$ is a projection valued function in
$C(X_i)\otimes M(B)_s$ for each $i$.
\end{thm}
\begin{proof}
Let $f$ be the element of $M(I)$ such that $\pi(f)=\mathbf{f}$.
 Without loss of generality, we can assume $f$ is self-adjoint and $0 \leq f \leq 1$.
\begin{itemize}
 \item [(i)]Suppose $X$ does not contain any infinite point.
 Choose a point $t_0 \in X $. Then there is a self-adjoint element $T \in M(B)$ such that $T-f(t_0) \in B$ and
 the spectrum of $T$ has a gap around $1/2$ by \cite[Theorem 3.14]{BP}. So we consider $f(t)+T-f(t_0)$ which is still self-adjoint whose image is $\mathbf{f}$.
 Thus we may assume $f(t_0)$ is a self-adjoint element whose spectrum has a gap around $1/2$.

 Since $r(f(t)):t \to f(t)-f(t)^{2}$ is norm continuous where $r(x)=x-x^2$,
 if we pick a point $z$ in $\left( 0,\frac{1}{4}\right)$ such that $z \notin \sigma(f(t_0)-f(t_0)^2)$,
 then  $\sigma(f(s))$ omits $r^{-1}(J)$  for $s$ sufficiently
 close to $t$ where $J$ is an interval containing $z$. In other
 words, there is  $\delta>0$ and $b> a >0 $ such that if $|t_0-s|<\delta $, then $ \sigma(f(s)) \subset
 [0, a)\cup (b,1]$. \\
  If we let $f_{t_0}(s)=\chi_{(b,1]}(f(s))$ for $s$ in $(t_0 - \delta,
  t_0+\delta)$ where $\chi_{(b,1]}$ is the characteristic function on $(b,1]$, then it is a continuous projection valued function such that
  $f_{t_0} - f \in C(t_0-\delta,t_0+\delta)\otimes B$.

  By repeating the above procedure, since $X$ is compact, we can find $n+1$ points $t_0,\cdots, t_n$, $n+1$ functions $f_{t_0}, \cdots,  f_{t_n}$,
  and an open covering $\{O_i\}$ such that $t_i \in O_i$, $O_i \cap O_{i-1} \ne
  \varnothing$, and  $f_{t_i}$ is projection valued function on $O_i$. Now
   let $f_i=f_{t_i}$ as above. Take the point $x_i \in O_{i-1} \cap O_i$ for
  $i=1, \cdots, n$. Then $f_{i}(x_i) - f_{i-1}(x_i)= f_{i}(x_i) - f(x_i) +
  f(x_i)- f_{i-1}(x_i) \in B$ and $f_0(x_0)-f_n(x_0) \in B$ if applicable. Let $X_{i}=[x_{i},x_{i+1}]$ for $i=1,
  \cdots, n-1$, $X_0=[0,x_1]$, and $X_n=[x_n,1]$. Since each $f_i$ is also
  defined on $X_{i}$, $(f_0, \cdots, f_n)$ is what we
  want.
  \item[(ii)] let $X$ be $[0, \infty)$.  Since $f^2(t) - f(t) \to 0 $ as $t$
  goes to $\infty$, for given $ \delta$ in $(0, 1/2) $, there is $M >0$ such that
  whenever $t \geq M$ then $\|f^2(t)-f(t) \| < \delta -\delta^2 $. It follows
  that $\sigma(f(t)) \subset [0, \delta) \cup (1-\delta,1]$ for $t \geq M$. Then
  again $\chi_{(1-\delta,1]}(f(t))$ is a continuous projection valued function for $t \geq M$
  such that $f(t) - \chi_{(1-\delta,1]}(f(t))$ vanishes in norm as $t$ goes to $\infty$.  By applying
  the argument in $(i)$ to $[0,M]$, we get a closed sub-intervals $X_i$ for $i=0,
  \cdots, n-1$ of $[0,M]$ and $f_{i} \in C_{b}(X_{i}, B(H))$. Now if we let
  $X_n=[M, \infty)$ and $f_n(t) = \chi_{(1-\delta,1]}(f(t)) $, we are done.
  \item[(iii)]The case $X=(-\infty, \infty)$ is similar to (ii).
  \end{itemize}
\end{proof}

When  a projection $\mathbf{f}\in C(I)$ is represented by $(f_0,f_1,\cdots,f_n)$  by Theorem \ref{T:locallift}, we note that  $f_i(x)$ is a projection in $M(B\otimes K)$ for each $x \in X_i$ and $f_i(x_i)-f_{i-1}(x_i)\in B$. Applying Definition \ref{D:BDF} we have $K$-theoretical terms $k_i=[f_i(x_i):f_{i-1}(x_i)]\in \KK(\mathbb{C},B)$ for $i=1,2,\cdots,n$. The following theorem shows that if all $k_i$'s are vanishing, then  a projection $\mathbf{f}$ in $C(I)$ lifts to a projection in $M(I)$.
\begin{thm}\label{T:lifting}
Let $I$ be $C(X)\otimes B$ where $B$ is a $\sigma$-unital, non-unital, purely infinite simple $C\sp{*}$-algebra such that $M(B)$  has real rank zero or $K_1(B)=0$ (See \cite{Zh}).
Let a projection $\mathbf{f}$ in $M(I)/I$ be represented by $(f_1,f_2,\cdots, f_n)$, where $f_i$ is a projection valued function in
$C(X_i)\otimes M(B)_s$ for each $i$, as in Theorem \ref{T:locallift}.
If $k_i=[f_i(x_i):f_{i-1}(x_i)]=0$ for all $i$, then the projection $\mathbf{f}$ in $M(I)/I$ lifts.
\end{thm}
\begin{proof}
Note that, by Zhang's dichotomy, $B$ is stable \cite[Theorem 1.2]{Zh}.
By induction, assume that $f_{j}(x_j)=f_{j-1}(x_j)$ for $j=1,2, \cdots, i-1$.

Let $f_i(x_i)=p_i, f_{i-1}(x_i)=p_{i-1}$. Since $[p_i:p_{i-1}]=0$, we have a unitary $u$ of the form ` identity + compact ' such that
$\|p_{i}-u^*p_{i-1}u\|<1/2$ by Theorem \ref{T:BDF}. Since $B$ has real rank zero,  given $0< \epsilon <1/4 $ there is a unitary $v \in \mathbb{C}1+B$ with finite spectrum such that $\|u-v \|<\epsilon$ \cite{Lin93}, \cite{Lin96}. Then \[ \|p_{i}-vp_{i-1}v^* \| \leq \|p_{i}-up_{i-1}u^* \|+\|up_{i-1}u^*-vp_{i-1}v^*\|<1.\]
Note that  $p_{i}-vp_{i-1}v^* \in B$. Thus we have $wp_{i}w^*=vp_{i-1}v^*$ for some unitary $w \in \id+B$. (Recall that Step 2 of the proof of Theorem \ref{T:BDF}.) Let $g_i=wf_iw^*$, then $f_i-g_i\in C(X_i)\otimes B $ since $w$ is of the form ` identity + compact'.

On the other hand, we can write $v$ as $e^{ih}$ where $h$ is a self-adjoint element in $B$ since $v$ has the finite spectrum. A homotopy of unitaries $t \to e^{ith}$, which are of the form `` identity + compact'', connects $1$ to $v$. Now we define $g_{i-1}(t)$ as
 \[\exp \left(i \frac{t-x_{i-1}}{x_i-x_{i-1}}h\right)f_{i-1}(t)\overline{\exp\left(i\frac{t-x_{i-1}}{x_i-x_{i-1}}h\right)}\] for $t \in [x_{i-1},x_i]$. Then
we see that $g_{i-1}(x_i)=g_{i}(x_i)$, $g_{i-1}-f_{i-1}\in C(X_{i-1})\otimes K$, and $g_{i-1}(x_{i-1})=f_{i-1}(x_{i-1})$. Moreover, if we let $g_{i+1}=wf_{i+1}w^*$, then $f_{i+1}-g_{i+1} \in C(X_{i+1})\otimes B $, and
\[
\begin{split}
[g_{i+1}(x_{i+1}):g_{i}(x_{i+1})]&=[wf_{i+1}(x_{i+1})w^*:wf_i(x_i+1)w^*]\\
&=[f_{i+1}(x_{i+1}):f_i(x_{i+1})]=0.
\end{split}
\]
Then $(f_0,f_1,\cdots, f_n)$ and  $(f_0,f_1,\cdots,g_{i-1},g_i,g_{i+1},f_{i+2},\cdots, f_n)$  define the same element  $\mathbf{f}$ while the $k_i$'s are unchanged and $i$-th discontinuity is resolved. So we take the latter as $(f_0,\cdots,f_n)$ such that $f_j(x_j)=f_{j-1}(x_j)$ for $j=1,\dots,i$.
We can repeat the same procedure until we have $f_i(x_i)=f_{i-1}(x_i)$ for all $i$. It follows that $(f_0,\cdots,f_n)$ is a projection in $M(C(X)\otimes B)$ which lifts $\mathbf{f}$.
\end{proof}
\begin{rem}
When $I=C_0(X)\otimes B$ where $X$ is $ [0,\infty)$ or $(-\infty,\infty)$, the similar result holds replacing $C(X_i)\otimes B$ with $C_0(-\infty,x_1]\otimes B$ or $C_0[x_n,\infty) \otimes B$  for $i=0$ or $i=n$ respectively.
\end{rem}
\section{Acknowledgements}
 Although this work was not carried out at Purdue, a significant influence on the auther was made by Larry Brown and Marius Dadarlat  who have aquainted him with geometric ideas in operator algebras.
He also would like to thank Huaxin Lin for answering the question related to Lemma 2.14.


\section{introduction}

The Standard Model explains the mass spectrum of fermions via their interaction with the Higgs field. The heavier the fermion is the stronger is its interaction with the Higgs field. Since the Higgs is a scalar this interaction has an attractive nature, which makes it tempting to look for new phenomena related to the heaviest known fermion, the top quark, as well as consider what may happen if heavier fermions of the next (fourth) generation exist in nature. It was suggested previously by Froggatt et al \cite{Froggatt:etal} that a system of 12 top quarks, 6 top and 6 anti-top occupying the lowest $S_{1/2}$ orbital can be bound by the forces originating from the Higgs exchange. This interesting suggestion was based on a simplified description, which neglected the Higgs mass. In addition, it was supported by estimations based on the Hydrogen-atom model. However, our more accurate calculations based on the nonrelativistic mean field  approach \cite{Kuchiev:2008fd}, which were supported by \cite{Richard:2008uq}, indicated that the binding is nonzero, albeit very small, only in the mentioned simplified approach, which neglects the Higgs mass. An account of the Higgs mass, which is limited by the current experimental bounds,  makes this 12 top-antitop system unbound. Recently Froggatt and Nielsen \cite{Froggatt:2008hc} discussed the issue anew.  This latter work claims again that the binding state of 12 top quarks exists and, moreover, the binding is very strong. This claim was supported by model-type approaches, which present a modification of the ones used previously in \cite{Froggatt:etal}.

Having in mind to clarify the issue this papers examines the fermion bound state for 12 fermions occupying $S_{1/2}$ shell considering the mass of the heavy fermion $m$ as a parameter. The  aim is  to find a critical value for this mass $m_{\mathrm cr}$ below which, when $m\le m_{\mathrm cr}$, the mentioned bound state cannot exist. The mean field approximation for the relativistic system that contains fermions interaction with the Higgs field is employed. The critical mass found depends strongly on the Higgs mass $m_{\mathrm H}$.
Taking the later in the interval $100\le m_{\mathrm H}\le 200$ Gev$/c^2$ we find the following restriction on the critical fermion mass $320 < m_{\mathrm cr}< 410$ Gev$/c^2$.

Consider $N$ heavy fermions that interact with the Higgs field $\Phi$. Take  the conventional unitary gauge where $\Phi$ is represented by the real field $ \xi$
\begin{equation}
\Phi\,=\,\frac{v}{\sqrt{2}}\left( \begin{array}{c} 0 \\ \xi \end{array}\right)~.
\label{Phi-xi}
\end{equation}
Here $v=246$ Gev is the VEV, which is achieved when $\xi=1$. The part of the Standard Model Lagrangian, which describes the interacting Higgs $\xi$ and fermion fields $\psi$ reads ($\hbar=c=1$)
\begin{align}
\!{\cal L}=
\frac{v^2}{2}\Big(\partial^\mu {\xi} \partial_\mu {\xi}\!-\!\frac{m_{\mathrm H}^2}{4}(\xi^2\!-\!1)^2\Big)\!+
\bar \psi( i\gamma^\mu \partial_\mu \!-m \xi ) \psi\,.
\label{L}
\end{align}
The important feature of the problem is that the fermion mass $m$ is presumed large, which allows us to neglect in the first approximation interactions with all other fields, gluons, $Z^{\,0}$ and $W$ mesons etc, which are omitted in (\ref{L}).
Moreover, the large fermion mass and the large number of fermions considered, $N=12$, allow one to rely in the lowest approximation on the mean field approach. Using it we replace the fields in Eq.(\ref{L}) by stationary wave functions calling them $\xi$ for the Higgs field and $\psi$ for the single-particle wave functions of fermions. Searching for the spherically symmetrical solution, take $\xi=\xi(r)$ and assume that all fermions occupy the same shell with total and angular momenta $j$ and $l$ (ultimately we will concentrate on the $S_{1/2}$ state). In this state the fermions are described by the same large $F(r)$ and small $G(r)$ components of the Dirac spinor. Using Eq.(\ref{L}) one writes the Hamiltonian $H$ of the system
\begin{align}
H\,=\,&\int_0^\infty\,
\Big[ \,\,\frac{v^2}{2}\,\Big(\,\xi'^2+\frac{1}{4}\,m_{\mathrm H}^2\,(\xi^2-1)^2
\,\Big)\,(4 \pi\,r^2)
\label{HFG}
\\
& + N\,\Big(\,2\big(F'+\frac{\varkappa}{r} F\,\big)\,G+m \xi\, (\,F^2-\,G^2\,)\,\Big)\,\Big]\,dr~.
\nonumber
\end{align}
Here $\varkappa=\pm\,(j+1/2)$ for $l=j\,\pm 1/2$, and normalization $\int_0^\infty(F^2+G^2)\,dr=1$ is taken.
From Eq.(\ref{HFG})  one derives the
mean-field equations for $\xi,F,G$
\begin{align}
\xi''+\frac{2}{r}\,\xi'&+\frac{m_{\mathrm H}^2}{2}\,\,\xi\,(1-\xi^2)\,=\,\frac{(N-1)m}{4\pi \,v^2}\,\frac{F^2-G^2}{r^2}~,
\label{xi}
\\
&(\varepsilon-m\,\xi)\,F~\,=\,-G'\,+\,(\varkappa/r)\,G~,
\label{F}
\\
&(\varepsilon+m\,\xi)\,G~\,=\,\,~~F'\,+\,(\varkappa/r)\,F~,
\label{G}
\end{align}
where the eigenvalue $\varepsilon$ is positive, and the factor $N-1$ in Eq.(\ref{xi}) takes into account the fact that each fermion interacts only with $N-1$ other fermions. These equations are supplemented by the boundary conditions, which for the localized fermion wave function have a conventional form
\begin{equation}
F(0)=F(\infty)=G(0)=G(\infty)=0~.
\label{FFGG}
\end{equation}
For the Higgs field the boundary conditions read
\begin{equation}
\xi'(0)\,=\,0~,\quad \quad \xi(\infty)\,=\,1~.
\label{xi0}
\end{equation}
The second one guarantees that outside the fermion bag the Higgs field takes its VEV,
while the first one makes the Higgs field regular at the origin by suppressing a singularity
at $r=0$ in the term $(2/r) \xi'$ in Eq.(\ref{xi}).
\begin{figure}[h]
\centering
\includegraphics[height=5.5 cm,keepaspectratio = true,
]{fig1}
\caption{ \label{ras}
Binding energy per fermion derived from the total energy $H$ in Eq.(\ref{HFG}) versus the mass $m$ of heavy fermions when 12 fermions
(6 fermions and 6 antifermions) occupy the same $S_{1/2}$ shell. The Higgs mass is $m_\text{H}=100$ Gev$/c^2$. The dots present calculations based on Eqs.(\ref{HFG})-(\ref{G}). An interpolating solid line shows that the bound state exists only for $m\ge m_{\mathrm cr}=324.2$ Gev$/c^2$.}
\end{figure}

The behaviour of the system described by Eqs. (\ref{HFG}) - (\ref{G}) was studied in Ref. \cite{Kuchiev:2008gt} for large fermion masses $m>1$ Tev$/c^2$, when the ultrarelativistic description is valid. It was found that the system reacts to the presence of the large fermion mass in such a way as to eradicate its influence on its physical parameters. In particular, the mass of the bag proves be independent on $m$. This phenomenon was described both by numerical methods and by the semi-analytical solution that is valid in the ultrarelativistic region.
Ref. \cite{Kuchiev:2008fd} studied the case of top quarks, when the nonrelativistic description is valid. In this situation the Higgs field proves be close to its VEV, $\xi\approx 1$, and consequently the linearized in $\xi$ version of Eq.(\ref{xi}) is applicable, which allows one to resolve Eq.(\ref{xi}) reducing the system (\ref{xi})-(\ref{G}) to a single equation on the nonrelativistic wave function of the fermion.

Here we consider the case when $m$ can be larger than the top quark mass, but remains below the ultrarelativistic limit,  $170<m<1000$ Gev$/c^2$. In this region the nonrelativistic approximation can become unreliable, while the mentioned analytical approach for the ultrarelativistic case is also not applicable. We are therefore to rely on  numerical methods solving Eqs.(\ref{xi})-(\ref{G}) directly, without additional simplifications. Fig. \ref{ras} shows the results of these calculations presenting the binding energy $E_{\text B}$ of the bag of $N=12$ fermions (6 fermions and 6 antifermions), which occupy the same $S_{1/2}$ shell. The binding energy per fermion shown in Fig. \ref{ras} equals $E_{\text B}/N =m-H/N$ where $H$ is the total energy (\ref{HFG}) calculated with the wave functions found from solution of Eqs. (\ref{xi})-(\ref{G}).
Fig. \ref{ras} presents this energy versus the mass of the fermion $m$. The positive binding indicates that the bound state is stable. From data presented in Fig. \ref{ras} we derive that the bag can exists provided $m\ge m_{\mathrm cr}\simeq 324.2$ Gev$/c^2$.
\begin{figure}[t]
\centering
\includegraphics[height=5.7 cm,keepaspectratio = true,
]{fig2}
\caption{
 \label{dva} Fields inside the fermion bag.
The Higgs and fermion masses are $m_\text{H}=100$, $m=328$ Gev$/c^2$.
Solid and dashed lines - large $F(r)$ and small $G(r)$ components of the Dirac spinor (scaled to make $\max(F)=1$), dot-dashed line - Higgs field $\xi(r)$, all solutions of Eqs.(\ref{xi})-\ref{G}), dotted line - nonperturbed Higgs $\xi_0=1$. The distance from the center of the bag $r$ is scaled by the VEV of the Higgs field $v=246$ Gev.}\end{figure}
\noindent

The wave functions $F(r)$, $G(r)$ and $\xi(r)$, which describe
the fermions and Higgs field inside the bag, are shown in Fig. \ref{dva}, which presents them for $m_\text{H}=100$, and the fermions mass $m=328$ Gev$/c^2$. This behaviour of the wave functions is typical when fermion mass just above the critical value.
In this region the Higgs field is only slightly modified inside the bag, while the large component of the fermion radial wave function $f(r)=F(r)/r$ is prominent at the center of the bag and smoothly decreases in the outside region. It is instructive to compare this behaviour with the case when the fermion is very heavy, $m>1$ Tev, which was considered in \cite{Kuchiev:2008gt}. In this ultrarelativistic case the Higgs field develops a node on the surface of the bag becoming negative inside, while the fermion wave function is localized on the surface of the bag, in the vicinity of the node of the Higgs field. Thus, for very heavy fermions the fermion wave function is suppressed at the center of the bag, the bag can be called `hollow'. In contrast, in the considered presently case of large, though moderately large fermion masses, the fermion wave function $f(r)$ is large at the origin, the bag is `solid'. Observe that Fig. \ref{dva} shows that the bag has no sharp, well defined surface, which contrasts the nuclear structure. This property complies with the fact that the range of forces between fermions $r\sim 1/m_{\mathrm H}$ induced by the Higgs exchange inside the bag is comparable with the size of the bag, while nuclear forces have a short-range nature in comparison to the size of large nuclei.

The latter comment reminds us that the  fermion binding energy may depend on the Higgs mass. To study this dependence the calculations were repeated for different values of $m_{\mathrm H}$ in the most interesting region $100\le m_{\mathrm H}\le 200$ Gev$/c^2$. For each $m_{\mathrm H}$ from this region a set of the fermion masses $m$ was taken and the critical fermion mass $m_{\mathrm cr}$ was found using a procedure illustrated by Fig. \ref{ras}. The found from these calculations critical fermion mass is plotted in Fig. \ref{tri} versus the Higgs mass. We see that the critical mass quickly rises with the Higgs mass, which complies with the fact that the heavy Higgs reduces the range in which an effective attraction between fermions exists. From data plotted i Fig. \ref{tri} we derive that if the Higgs mass is restricted by the interval $100\le m_{\mathrm H} \le 200$ Gev$/c^2$ then the critical fermion mass lies in the interval  $320 <m_\mathrm{cr}<410$ Gev$/c^2$.

\begin{figure}[t]
\centering
\includegraphics[height=5.5 cm,keepaspectratio = true,
]{fig3}
\caption{
 \label{tri} The critical fermion mass $m_{\mathrm cr}$ versus the Higgs mass $m_{\mathrm H}$. For heavy fermions with $m>m_{\mathrm cr}$ there can exist a bag of 12 fermions
occupying $S_{1/2}$ shell, which are bound via the Higgs field.
The dots present data from numerical solution of Eqs.(\ref{xi})-(\ref{G}), the solid line is an interpolation.}
\end{figure}
\noindent

Assessing the accuracy of these results one should keep in mind errors related to several phenomena neglected in the approximation considered. An important class among them represent radiative corrections, which importance was emphasized previously in
\cite{Dimopoulos:1990at,Bagger:1991pg,Farhi:1998vx,Crichigno:2009kk}.
It was alleged that these corrections are able to make the bag shallow
\cite{Dimopoulos:1990at,Crichigno:2009kk},  deflate it into inexistence \cite{Farhi:1998vx}, or ultimately even make the vacuum unstable \cite{Bagger:1991pg}. It seems that the issue is important and need further consideration; we are planning to address it in some detail in the near future. Presently, one can state that the reported results \cite{Dimopoulos:1990at,Bagger:1991pg,Farhi:1998vx,Crichigno:2009kk} indicate on a repulsive nature of radiative corrections. Therefore in any case they cannot reduce the estimates for the critical mass of heavy fermions.

There exist a number of other corrections. One class of them represent the many-body corrections, which include deviation of the many-body wave function from the mean field approximation, as well as the retardation and recoil corrections. One should anticipate though that the considered system of large number of tightly bound fermions, 12 fermions located in the region of the size of $r\sim \hbar c/v$, $v=246$ Gev$/c^2$, makes the mean field approach sensible, many body corrections should not be very pronounced. There exit also corrections that arise due to interactions with gluons, $Z^0$ and $W$ bosons etc, but for the considered prominent fermion masses $>300$ Gev$/c^2$ these interactions are certainly much weaker than the one induced by the Higgs exchange. Speaking of which one could remember that there exists also an interaction between fermions, which is induced by the virtual fermion-antifermion annihilation into the Higgs. However, for the considered heavy fermions this channel is strongly suppressed due to the fact that the virtual Higgs, which arises from an annihilation, propagates with large virtual frequency $\omega\simeq 2m$, which greatly reduces contribution of this process(similar situation reveals itself in the positronium, where the annihilation channel provides only a small correction.)

This discussion shows that there is no reason to unticipate that any of the mentioned phenomena would be able to change the main qualitative conclusion of this work.  The critical mass for fermions below which they are certainly not able to produce a bag (12 fermions in $S_{1/2}$ shell bound by the Higgs exchange) is large, $m>320$ Gev$/c^2$. This bound greatly exceeds the mass of the top quark.  Consequently 12 top quarks are not able to form the bound state. This important qualitative conclusion complies with previous calculations \cite{Kuchiev:2008fd,Richard:2008uq} and contradicts results of Froggatt and Nielsen \cite{Froggatt:2008hc}. In order to find a bridge  between conclusions of \cite{Kuchiev:2008fd,Richard:2008uq} and \cite{Froggatt:2008hc}
Froggatt and Nielsen put forward an idea that there may exist a particular phase transition, a phenomenon which manifests itself as an existence of two different states, which have the same quantum numbers. These authors alleged that one of these states, which is bound for top quarks, is studied in \cite{Froggatt:2008hc}, while another one, which remains unbound, presents itself in calculations of \cite{Kuchiev:2008fd,Richard:2008uq} and, developing this argument, in the present work. This idea is supported in \cite{Froggatt:2008hc} by arguments derived from the proposed in this work `toy model', which makes a number of simplifications.

In order to check the idea that there may exist some other bound state, in which
the Higgs field is also spherically symmetrical and the fermion wave function has the same $S_{1/2}$ symmetry, the following numerical approach was used. A set of initial approximations for the Higgs field was taken. This set included, in particular, the wave functions for the Higgs field, which described it as being mostly flat inside the bag, in accord with presumptions of the `toy model' of \cite{Froggatt:2008hc}. Each initial approximation was used as input data for a numerical procedure, which generated the self-consistent numerical solution of the mean field equations. It was found that for any sufficiently  smooth initial approximation the numerical procedure converges, and converges to one and only self-consistent solution, regardless of the initial approximation. Thus, there was found no indications of a possible existence of a different bound state with the given symmetry. Summarizing, we have a situation when direct numerical calculations of this work, as well as of \cite{Kuchiev:2008fd,Richard:2008uq} show that top quarks are too light to be able to produce the bound state in question.

In conclusion, we investigated a possibility that 12 heavy fermions, which occupy the same $S_{1/2}$ shell can be bound by an attraction, which originates from the Higgs field exchange. The fermions in question need to be heavy to make this attraction sufficiently strong. The critical mass for fermions, below which this bound state cannot exist, is found using the relativistic mean field approximation. A prominent dependence of this critical mass on the mass of the Higgs field is found. Taking the later in the interval $100\le m_{\mathrm H}\le 200 $ Gev$/c^2$ the critical fermion mass is found in the interval $320\le m_{\mathrm H}\le 410 $ Gev$/c^2$. The fact that even the lowest estimate from this set is almost twice as heavy as the mass  of top quarks makes it unrealistic to expect that top quarks are able to produce the mentioned bound state. The only hope to observe experimentally the multi-fermion bound states, which exists due to the Higgs exchange, could rest on a hope that there exist heavy fermions of the next, fourth generation.

This work is supported by the Australian Research Council