\section{INTRODUCTION}
Impressive progress of the physics of ultracold quantum gases have opened new pathways for the
study of low-dimensional few-body systems (see, for example \cite{Chin2010,Kohler2006})
as well as strongly correlated many-body systems \cite{Bloch2008,Yurovsky2008}).
Specifically, it was shown that the confining geometry of atomic traps can drastically change
the scattering properties of ultracold atoms and induce resonances in the collisions (confinement-induced resonances (CIRs)) \cite{Olshanii1}.
The CIR for bosons has been found to occur in the vanishing energy collisional limit
as a consequence of the coincidence of the binding energy of a diatomic molecular state
with the energy spacing between the levels of the confining (harmonic)
potential \cite{Olshanii2, Melezhik07,Saeidian}. It was shown that this coincidence leads
to a divergence of the effective interatomic coupling constant $g_{1D}$ and to a total atom-atom reflection, which appears as
a broad dip in the transmission coefficient T, thereby creating a gas of impenetrable bosons \cite{Girardeau}.   CIRs have also been extensively
studied e.g. in the context of three-body \cite{Mora1, Mora2}, and four-body \cite{Mora3} scattering in confining traps, fermionic p-wave
scattering \cite{Granger}, distinguishable atom scattering \cite{Kim2, Kim3, Melezhik07, Melezhik09} or multichannel scattering \cite{Saeidian, Melezhik11,Olshanii3} in atomic waiveguides.
Two novel effects were predicted for distinguishable atoms: the so-called dual CIR yielding
a complete suppression of quantum scattering \cite{Kim2}, and the resonant molecule formation in tight waveguides \cite{Melezhik09}.
D-wave resonant scattering of bosons in confining harmonic waveguides has been analyzed very recently in \cite{Giannakeas}.
Remarkable experimental progress has lead to the observation of CIRs,
for identical bosons \cite{Kinoshita, Paredes,Haller2009,Haller2010} and fermions\cite{Guenter,Moritz,Kohl2011}, as well as to
distinguishable atoms \cite{Lamporesi2010}.

However, the recent experimental \cite{Lamporesi2010,Haller2010,Kohl2011,Haller2011}
and theoretical \cite{Peng2010,Zhang2,Melezhik11,Sala2011} investigations of the CIRs clearly
show that, despite the impressive progress, the existing theoretical models of CIRs need to be improved for
a quantitative description of the experiments in this field.
The obvious drawback of the existing theoretical models for CIRs in atomic waveguides is the single-channel character of the
simple interatomic interactions employed. Single-channel potentials with zero-energy
bound states were used for simulating magnetic Feshbach resonances in free space representing a
necessary ingredient for the appearance of the CIR in a confining trap(see, for example \cite{Chin2010,Yurovsky2008}). In the seminal work of
Olshanii \cite{Olshanii1} and in subsequent papers \cite{Olshanii2,LowDSystems}, including the recent works
\cite{Peng2010,Zhang2}, the simple form of a pseudo-potential was used modeling the interatomic interactions. More realistic potentials
were used in our previous works \cite{Kim2,Kim3,Saeidian,Melezhik09,Giannakeas,Melezhik11} as well as
in the works of other authors \cite{Olshanii2,Granger,Sala2011,Grishkevich} which, however, all possess a single-channel character.
The single-channel interatomic interaction approach permits to explore only the main attribute of the
Feshbach resonances in the 3D free space, namely the appearance of a singularity in the s-wave
scattering length $a_s \rightarrow\pm\infty$ when the molecular bound state with energy $E_B$ crosses the atom-atom
scattering threshold at energy $E=0$ in the entrance channel. However, other important
parameters of the Feshbach resonance, such as the rotational and spin structure of
the molecular bound state in the closed channel as well as the width
$\Delta$ of the resonance characterizing the coupling $\Gamma$ of the molecular state with the
entrance channel, were ignored. The main goal of the present work is to extend these theoretical approaches developed
earlier for the CIRs and transverse excitations/deexcitation processes for collisions
in harmonic waveguides \cite{Saeidian} to the case of a tensorial interaction
as well as taking into account the width of the magnetic Feshbach resonances responsible for the CIR. The parameters
obtained from the experimental analysis of the magnetic Feshbach resonances in free space,
namely the resonant energies $E_{c,i}$ (or the corresponding values of the field strengths $B_{c,i}$ of the external magnetic field),
the widths of the resonances $\Delta_i (\Gamma_i)$, spin characteristics and the background scattering length $a_{bg}$,
are used as input parameters of our model.

It is well-known that the most adequate computational schemes
including all the above mentioned parameters of the Feshbach
resonances are multi-channel scattering
\cite{Chin2010,Chin2004,Lange} and multi-channel quantum-defect
\cite{Mies2000} approaches providing a quantitative description
for a broad regime of experimental parameters for the Feshbach
resonances in ultracold atomic gases. In the present work the
two-channel potential scheme \cite{Lange}, permitting the
efficient modeling of the interatomic interaction near the known
magnetic Feshbach resonances of the Cs gas in free space
\cite{Chin2004}, is included in our multichannel approach
\cite{Saeidian} instead of the simplified interparticle
single-channel potential we have used so far for analyzing CIRs
and transverse excitations in atomic waveguides. With this
approach we model s-, d- and g-wave Feshbach resonances in
ultracold Cs (s-,d- and g- indicate here the "exit" channels, i.e.
they correspond to rotational quantum numbers of the molecular
states in the closed channels), which were observed in free space
scattering experiments \cite{Chin2004}, and quantitatively analyze
the shifts and widths of the resonances in the harmonic waveguides
at experimental conditions closed to the ones encountered in the
works \cite{Haller2009,Haller2010}. We note that selected aspects
of the modeling of magnetic Feshbach resonances in  optical traps and lattices
were also considered in \cite{Dickerscheid,Diener,Yurovsky2006,Schneider,Peng}.

\section{FESHBACH RESONANCE MODEL}

We consider the collision of two identical bosonic atoms in a harmonic waveguide.  This two-body problem permits the separation of the center-of-mass
and relative motion yielding the following Hamiltonian for the relative atomic motion
\begin{equation}
\hat{H}(r,\theta)=[-\frac{\hbar ^2}{2\mu}\nabla^2+ \frac{1}{2}\mu\omega_{\perp}^2\rho^2]\hat{I} + \hat{V}(r)
\end{equation}
with $\rho = r\sin\theta$ and the trap potential $1/2\mu\omega_{\perp}^2\rho^2$.
$\hat{V}(r)$ is the 4-channel interatomic potential and $\hat{I}$ is the unit matrix,
$r$ is the relative radial coordinate and $\mu = m/2$ is the reduced mass of the atoms.

Following the scheme suggested in \cite{Lange} for describing the three magnetic Feshbach resonances in an ultracold Cs gas let us suppose
that initially  the scattering atoms are prepared in one spin configuration $|e\rangle$ (called the
``entrance channel") and the ``closed channels" $|c,i\rangle$ (i=1, 2, 3) support s-, d- and g-wave molecular bound states.
at -11.1G, 47.78G and 53.449G.
The quantum state of an atomic pair with energy $E$ is described as
$$|\psi\rangle=\sum_{i=1}^3\psi_{c,i}({\bf r})|c,i\rangle+\psi_{e}({\bf r})|e\rangle$$
satisfying the Schr\"{o}dinger's equation with the Hamiltonian (1).  A four-channel square-well potential
\begin{equation}
\hat{V}=\left (\begin{array}{cccc} -V_{c,3}& 0 & 0 & \hbar\Omega_3\\ 0 & -V_{c,2} & 0 &  \hbar\Omega_2\\
  0&0&-V_{c,1}  & \hbar\Omega_1\\ \hbar\Omega_3 &  \hbar\Omega_2 & \hbar\Omega_1& -V_e
\end{array}\right )  \begin{array}{cc} (if & r<\overline{a})\end{array}\\
\end{equation}
\begin{equation}\nonumber
=\left (\begin{array}{cccc} \infty& 0 & 0 &0\\ 0 & \infty & 0 & 0\\
  0&0& \infty&0\\0 & 0&0& 0
\end{array}\right )  \begin{array}{cc} (if & r>\overline{a})\end{array}
\end{equation}
is employed to describe the colliding atoms in the ``entrance channel"
$|e\rangle$ and the weakly-bound molecules in the ``closed channels" $|c,i\rangle$ near a Feshbach
resonance. For $r< \overline{a}=4\pi\Gamma(1/4)^{-2}R_{vdW}$, we assume the attractive potential can support multiple molecular states - that is,
$ V_{e}$,$ V_{c_{i}}\gg E_{vdW}=\hbar^2/2mR_{vdW}^2$ - and $\hbar\Omega_i$ induce Feshbach couplings
between the channels \cite{Lange}. The size $\overline{a}$ of the potential action is chosen to account for the interatomic interaction
determined by the van-der-Waals (vdW) tail with the length $R_{vdW}=1/2(2mC_6/\hbar^2)^{1/4}$ and the energy
scale $E_{vdW}=\hbar^2/2mR_{vdW}^2$ \cite{Gribakin},
where $C_6$ is the corresponding van der Waals coefficient and $\Gamma(x)$ is the gamma function. For $r> \overline{a}$, entrance- and closed-channel
thresholds are set to be $ E=0$ and $E=\infty$, respectively (See Fig.1 in Ref.\cite{Lange}).

Such a choice of the interatomic interaction permits a simple
parametrization of the atom-atom scattering in universal terms
of the energy of the i-th bare bound state $E_i$,  the Feshbach
coupling strength of the bound molecular state with the entrance
channel $\Gamma_i$ and the background scattering length $a_{bg}$,
which is convenient for an analysis of experimental data near
magnetic Feshbach resonances \cite{Chin2010,Lange}. When the
mixing between the closed channels and the entrance channel is
weak and the background scattering length $ |a_{bg}|$ considerably
exceeds the range of the interatomic interaction $\overline{a}$, the
s-wave scattering length $a$ and the binding energy $E_b$ in free space are given by \cite{Lange}:
\begin{equation}
\frac{1}{a-\overline{a}}=\frac{1}{a_{bg}-\overline{a}}+\frac{1}{\overline{a}}\sum_{i=1}^{3}\frac{\Gamma_i/2}{E_i},
\end{equation}
and
\begin{equation}\label{binding energy}
E_b=\frac{\hbar^2k_m^2}{2\mu}\,, \,\,\, k_m=\frac{1}{a_{bg}-\overline{a}}+\frac{1}{\overline{a}}\sum_{i=1}^{3}\frac{\Gamma_i/2}{E_b+E_i}
\end{equation}
respectively.
Then, assuming that the bare bound  states can be linearly tuned magnetically by a linear Zeeman shift -namely, $ E_{i}=\delta\mu_{i}(B-B_{c,i})$, where $\delta\mu_{i}$
is the relative magnetic moment between the entrance and i-th channels and $ B_{c,i}$ is the crossing field value of the $i$th bare bound
state, the scattering length $a(B)$ can be represented as
\begin{equation}\label{Scatt length eq}
\frac{a}{a_{bg}}=\prod^{3}_{i=1}\frac{B-B^*_i}{B-B_{0,i}}\,\,.
\end{equation}
Here $B_{0,i}$ is the $i$th lowest pole of $a$ and $B^*_{i}$ the $i$th lowest zero.  The width of the $i$th Feshbach resonance can be defined as
$\Delta_{i}=B_{i}^*-B_{0,i}$.  The binding energy of Feshbach molecules can be measured in the
experiment by e.g. radio frequency and microwave
spectroscopy \cite{Regal, Bartenstein, Claussen, Thompson, Papp}.  Using the fitting parameters $a_{bg}$, $\delta\mu_{i}$, $\Gamma_{i}$ and $B_{c,i}$,
one can fit the experimental data \cite{Regal, Bartenstein, Claussen, Thompson, Papp} with  (\ref{binding energy}), from which one can calculate
$a(B)$ (3), $B^*_i$ and $B_{0,i}$ (5) \cite{Lange}.

\begin{figure}[!t]
\includegraphics[width=1.1\columnwidth]{1.eps}
\caption{(Color online) The s-wave scattering length $a$ of  $|F=3, m_F=3\rangle$ Cesium atoms as a function of the magnetic field  $B$.  The solid curve shows the analytical result (5) and the
dots show the numerical result (see the text).}
\label{fig1}
\end{figure}

\begin{figure}[hbt]
\includegraphics[width=1.1\columnwidth]{2.eps}
\caption{(Color online) The transmission coefficient $T$ for the harmonic trap with $\omega_{\perp}=14.9$ kHz as a function of the
magnetic field $B$ for the case of the 2-channel (solid curve) and
4-channel (dots) potential (2).} \label{fig2}
\end{figure}

The result of the procedure is given in the  Table I of ref.
\cite{Lange} for ultracold $Cs$ Feshbach molecules where  $
a_{bg}=1875 a_{0}$ and $ \overline{a} = 95.7 a_{0}$ ($a_0$ is the
Bohr radius) are also given. With these values for $B_{0,i}$ and $B_{i}^*$ we
show the scattering length $a(B)$ as a function of the
magnetic field $B$ according to Eq.(\ref{Scatt length eq}) in Fig.1.

The parameters $ a_{bg}$, $\delta\mu_{i}$, $\Gamma_{i}$, $B_{c,i}$, $B^*_i$ and $B_{0,i}$ from
\cite{Lange} together with $a(B)$ defined by Eq.(5) are used for fitting
the diagonal terms $V_{c,i}$ and $V_e$ in the tensor potential (2). The nondiagonal terms $\hbar\Omega_i$
are defined by the formulas \cite{Lange}
$$\Gamma_{i}/2=2\theta^{2}V_{c,i}\,,\,\,\,
\tan 2\theta_i=\frac{2\hbar\Omega_i}{V_e-V_{c,i}}\,\,.$$ The
scattering length $a(B)$ is then calculated for different $B$ and varying parameters of the potential $\hat{V}$ by solving the
Schr\"odinger equation
\begin{equation}
\left ([-\frac{\hbar^2}{2\mu}\nabla^2+ \frac{1}{2}\mu\omega_{\perp}^2\rho^2]\hat{I}+\hat{B} +\hat{V}(r)\right )|\psi\rangle=E|\psi\rangle
\end{equation}
in free space ($\omega_{\perp}=0$) with the scattering boundary conditions
$$\psi_e({\bf r})\rightarrow \exp\{ikz\} + f(k,\theta)/r\exp\{ikr\}\,,\,\,\, \psi_{c,i}({\bf r})\rightarrow 0$$
at $kr\rightarrow\infty$ for the fixed $E\rightarrow 0$  $ (k=\sqrt{2\mu E}/\hbar\rightarrow 0)$ \cite{MelHu}. The diagonal
matrix $\hat{B}$ in (6) is defined as $B_{ii}=\delta\mu_i(B-B_i)$ $(i=1,2,3)$ and $B_{ee}=0$.
After separation of the angular part in (6) we come to the system of four coupled radial equations
\begin{equation}
[-\frac{\hbar^2}{2\mu}\frac{d^2}{dr^2} -\frac{\hbar^2l_{\alpha}(l_{\alpha}+1)}{2\mu r^2} +B_{\alpha\alpha}]\phi_{\alpha}(r) + \sum_{\beta}V_{\alpha\beta}(r)\phi_{\beta}(r)= E\phi_{\alpha}(r)
\end{equation}
for the radial part $\phi_{\alpha}(r)$ of the desired wave
function
$|\psi\rangle=\sum_{\alpha}\psi_{\alpha}({\bf{r}})|\alpha\rangle =
\sum_{\alpha}\phi_{\alpha}(r)Y_{l_{\alpha}0}(\hat{r})|\alpha\rangle$,
where $\alpha=\{e,i=1,2,3\}$, $(l_e=0, l_1=0, l_2=2, l_3=4)$ and the matrix
elements $V_{\alpha\beta}(r)$ are defined by Eq.(2). The
centrifugal barrier  $-\frac{\hbar^2l_{\alpha}(l_{\alpha}+1)}{2\mu
r^2}$ in (7) models at $r\rightarrow 0$ the correct asymptotic
behavior $\phi_{i}(r) \sim r^{l_i+1}$ of the molecular bound
states $|c,i\rangle$ in the closed channels which couple to the
entrance s-wave ($l_e=0$) channel $|e\rangle$ by the nondiagonal terms
$V_{\alpha\beta}(r) (\alpha\not=\beta)$.

By varying the
$V_{c,i}$, $V_e$ and $\Omega_i$ we obtain an excellent agreement of the
calculated s-wave scattering length $a(B)$ with the analytical results from
\cite{Lange} for Cesium atoms in the hyperfine state $|F=3, m_F=3\rangle$ for the considered magnetic
field regime $-40G<B<60G$ (see dots in Fig.1).  In this regime,
we observe three resonance terms, which correspond to the coupling to the s-, d-, and g-wave molecular states.

Next we analyze the scattering properties of the s-, d- and g-wave
magnetic Feshbach resonances in harmonic waveguides by integrating the Schr\"odinger equation (6)
for $\omega_{\perp}\neq 0$ with the scattering boundary conditions
\begin{equation}
\psi_e({\bf r}) = \left (\cos (k_0 z) + f_e \exp\{ik_0 \mid z\mid\}\right
)\Phi_{0}(\rho)\,\,,\,\,\, \psi_{c,i}({\bf r})\rightarrow 0
\end{equation}
at $\mid z\mid = \mid r\cos\theta\mid \rightarrow \infty$ adopted
for a confining trap \cite{Saeidian}. $f_e(E)$ is the
scattering amplitude, corresponding to the symmetry with respect to the exchange $z
\rightarrow -z$ (we consider collisions of identical bosonic Cs
atoms), $\Phi_0(\rho)$ is the wave function of the ground-state of
the  two-dimensional harmonic oscillator  and $k_0=\sqrt{2\mu
(E-\hbar\omega_{\perp})}/\hbar=\sqrt{2\mu E_{\parallel}}/\hbar$.
In the presence of of the harmonic trap ($\omega_{\perp}\not=0$) the problem (6),(8) becomes non-separable
in the plane $\{\rho,z\}$, i.e. the azimuthal angular part is
separated in the wave function $|\psi\rangle$ and Eq.(6) is
reduced to the coupled system of four 2D Schr\"{o}dinger-type
equations. To integrate  this coupled
channel 2D scattering problem in the plane $\{r,\theta\}$ we have extended the computational
scheme \cite{Saeidian}.

The computations have been performed in a range of variation of
$\omega_{\perp}$ close to the experimental values of the trap
frequencies $\sim 2\pi \times 14.5$kHz \cite{Haller2010}. We have
integrated Eq.(6) for varying $B$ and fixed longitudinal colliding
energy $E_{||}=E-\hbar\omega_{\perp}$. In the main part of
computations the energy $E_{||}$ was chosen very low
$E_{||}=1.0\times 10^{-16}(\frac{\hbar^2}{\mu\overline{a}^2})$
$\rightarrow 0$ $(k_0=1.0\times 10^{-8}(\frac{1}{a_0})\rightarrow
0)$ to have a possibility for direct comparison with existing
pseudopotential estimates obtained in zero-energy limit. The
integration was performed in the units of the problem leading to
the scale transformation: $r \rightarrow \frac{r}{\overline{a}}$,
$E \rightarrow \frac{E}{E_0}$, $V \rightarrow \frac{V}{E_0}$, and
$\omega_{\perp} \rightarrow \frac{\omega_{\perp}}{\omega_0}$ with
$E_0=\frac{\hbar^2}{\mu\overline{a}^2}$, and
$\omega_0=\frac{E_0}{\hbar}$.

\section{RESULTS AND DISCUSSION}

\subsection{Transmission coefficient}
In Fig.2 we present the transmission coefficient $T(B)=|1+
f_e(B)|^2$ as a function of B calculated for the harmonic trap
$\omega_{\perp}=14.9$kHz with the two-channel (i=1) and
four-channel (i=1,2,3) character of the tensorial
interaction $\hat{V}(r)$ (2).
\begin{figure}[t]
\includegraphics[width=1.\columnwidth]{3.eps}
 \caption{(Color online) The coupling constant $g_{1D}$ as a function
of the scattering length $a$ in free space in the region
near the s-wave Feshbach resonance, calculated for the harmonic trap with $\omega_{\perp}=14.9$kHz. Dots indicate the numerical
data and the solid curve the analytical results with
the formulas derived in the pseudopotential approach
\cite{Olshanii1}.} \label{fig3}
\end{figure}
The two-channel potential $\hat{V}(i=1)$ supports only one broad
s-wave Feshbach resonance
while the four-channel potential
$\hat{V}(i=1,2,3)$ supports all three s-, d- and g-wave magnetic
Feshbach resonances in free space ($\omega_{\perp}=0$) at
the fields $B_{0,1}=-11.1$G, $B_{0,2}=47.78$G and $B_{0,3}=53.449$G,
respectively \cite{Lange}. All these resonances are also developed
in the calculated curve $T(B)$ of the transmission coefficient in
the trap ($\omega_{\perp}\neq 0$).

\begin{figure}[hbt]
\includegraphics[width=1.\columnwidth]{4.eps}
\caption{(Color online) The transmission coefficient $T(B)$ as a function of the
strength of the magnetic field for several trap frequencies
$\omega_{\perp}$. Here the s-wave scattering length $a(B)$ from
Fig.1 is also given (solid line).} \label{fig4}
\end{figure}
\begin{figure}[hbt]
\includegraphics[width=1.\columnwidth]{5.eps}
\caption{(Color online) The dependence of the s-wave scattering length $a(B)$ on
the width $a_{\perp}$ of the waveguide at the points $B_{min}$ of the minimum of the
transmission coefficient $T(a(B),a_{\perp})$ (see Fig.4).
The dots, pluses and cycles correspond to the calculated points
near the s-,d- and g-wave magnetic Feshbach resonances, respectively. The solid
curve corresponds to the formula $a=a_{\perp}/C$
\cite{Olshanii1}.} \label{fig5}
\end{figure}

\begin{figure}[hbt]
\includegraphics[width=1.1\columnwidth]{6a.eps}
\includegraphics[width=1.\columnwidth]{6b.eps}
\caption{(Color online) (a) The dependence of the width $\Gamma_i^*$ on the
longitudinal momentum  $k_0$ near the d-wave Feshbach resonance.
The solid line has been obtained via Eq.(10), solid circles
indicate the widths $\Gamma_i^*(k_0)$ extracted from the
numerically calculated $T(B,k_0)$.(b) The transmission coefficient
T as a function of the magnetic field B calculated for a few
$k_0$. For both subfigures $\omega_{\perp}$=14.9kHz.}
\label{fig55a}
\end{figure}

\begin{figure}[hbt]
\includegraphics[width=1.\columnwidth]{7.eps}
\caption{(Color online) The dependence of the width $\Gamma_i^*$ on the trap
frequency $\omega_{\perp}$ near the d-wave Feshbach resonance at
$k_0$=2.0899$\times 10^{-8}(\frac{1}{k_0})$. The solid curve has been obtained via Eq.(10), solid circles indicate the widths
$\Gamma_i^*(\omega_{\perp})$ extracted from the numerically calculated
$T(B,\omega_{\perp})$.   The analytical (dashed line) and numerical (squares) results for the right half width $1/2\Gamma^{(i)}_{CIR}$
versus $\omega_{\perp}$ are also given.} \label{fig55b}
\end{figure}

First, we have analyzed the region of $B\simeq-11.1$ G near the
s-wave resonance. In Fig.3 we show the calculated 1D coupling
constant $g_{1D}=\lim_{k_0\rightarrow
0}Re\{f_e(k_0)\}/Im\{f_e(k_0)\}k_0/\mu$ \cite{Olshanii2} as a
function of the scattering length $a(B)$ in this region.  By
tuning $B$ in the interval $-60G \leq B\leq 40G$ , the s-wave
scattering length $a(B)$ changes from $-\infty$ to $+\infty$. At
$B=-49.78$G corresponding to the point $a(B)=a_{\perp}/C $ (where
$C=1.4603..$ and $a_{\perp}=\sqrt{\hbar/\mu\omega_{\perp}}$) of
the appearance of the CIR \cite{Olshanii1}, the coupling constant
$g_{1D}(a/a_{\perp})$ diverges and the behavior is in very good
agreement with our previous computation of $g_{1D}(a/a_{\perp})$
performed with the single channel screened Coulomb potential
\cite{Saeidian} and the formulas $g_{1D}=2\hbar a/(\mu
a_{\perp}^2)/(1-C a/a_{\perp})$ derived in a pseudopotential
approach \cite{Olshanii1}.

Let us next analyze the region of the d-wave magnetic Feshbach
resonance near the point $47.78$G (see Fig.1) - the region of
major experimental interest due to atomic loss and the
formation of Cs molecules in the confined ultracold gas of Cs
atoms \cite{Haller2009,Haller2010}. The results of our
computations are illustrated by the curves $T(B)$ calculated for
different $\omega_{\perp}$ (Fig.4) corresponding to the transverse
frequencies of the optical trap being used in the
experiment \cite{Haller2010}. It is clearly shown that the
position $B_{min}$ of the transmission coefficient $T(B)$ minimum
is dependent on the trap width
$a_{\perp}=\sqrt{\hbar/(\mu\omega_{\perp})}$ and the corresponding
scattering length $a(B_{min})$ at the point $B_{min}$ is
accurately described by the formulas $a(B_{min}) = a_{\perp}/C$
obtained by Olshanii \cite{Olshanii1} for the position of the CIR
in a harmonic waveguide. The latter fact is illustrated in Fig.5. Here
the results of the calculation of the dependence of $a(B_{min})$
on  $a_{\perp}$ for s-wave and g-wave Feshbach resonances are also
given. This analysis clearly demonstrates that the law
$a=a_{\perp}/C$ for the position of the CIR in a harmonic
waveguide is fulfilled with high accuracy for the Feshbach
resonances of different tensorial structure, although the law was
initially obtained at zero-energy limit in the framework of an s-wave single-channel
pseudopotential approach to the interatomic interaction
\cite{Olshanii1}. We also see that the positions $B_{max}$ of the
maximums of the coefficients $T(B)$ calculated for different
$\omega_{\perp}$ are independent
of $\omega_{\perp}$ and coincide at the point $B_{2}^*=47.944$G of the zero of the scattering length in
free space ($a(B_{2}^*)=0$). This fact is in agreement with
the analytic result $T=\mid 1+f_e\mid^2\rightarrow 1$ obtained in
the pseudopotential approach at $a\rightarrow 0$. The same
behavior of the $T(B)$ coefficients has been found in the
vicinity of 18.1G and 53.46G (points $B_1^*$ and $B_3^*$ of s-wave
scattering length zeros near the s- and g-wave Feshbach
resonances).

\subsection{Width of the resonant enhancement of the transmission}

While the coefficients $T(B)$ in the vicinity of the points $B_i^*$ (positions of the zeros of the scattering length $a(B)$)
show a resonant behavior one can define the width $\Gamma_{i}^*$ of this resonance as the width at half $T$ maximum.
By using the formulas for the even scattering amplitude in the trap \cite{Schmiedmayer,Olshanii3}
$$
f_e=-\frac{1}{1+i\cot\delta_{1D}} = -\frac{1}{1+i k_0 a_{1D}}
$$
valid near the CIR in the zero-energy limit and the definitions $a_{1D}=\frac{a_{\perp}}{2}(C-\frac{a_{\perp}}{a})$ and
$a=a_{bg}\prod^{3}_{i=1}\frac{(B-B_i^*)}{(B-B_{0,i})}$,
we obtain
\begin{equation}
\Gamma_{i}^* = \frac{4 a_{bg}\gamma_i a_{\perp}^2k_0\Delta_i}{4
a_{bg}^2 \gamma_i^2-a_{\perp}^4k_0^2(1- C\gamma_i
\frac{a_{bg}}{a_{\perp}})^2} \,\,,
\end{equation}
where $\gamma_i=\prod_{j\neq
i}^3\frac{(B_i^*-B_j^*)}{(B_i^*-B_{0,j})}$. The above formula is
valid for the condition $\Gamma_i^*\ll B_i^*$ of narrow resonances,
which is fulfilled with high accuracy for the d- and g-wave resonances
and less accurately for s-wave resonance (see Fig.2). In the low
energy limit $k_0\rightarrow 0$ the expression is reduced to
\begin{equation}
 \Gamma_{i}^* = \Delta_i \frac{a_{\perp}^2 k_0}{a_{bg}\gamma_i} =
 \Delta_i\frac{\sqrt{2E_{\parallel}}}{\sqrt{\mu}a_{bg}\omega_{\perp}\gamma_i} \,\,\,,
\end{equation}
where $\Delta_i = B_i^*-B_{0,i}$ is the width of the Feshbach
resonance in free space and the dimensionless  $\gamma_i$- factor
$\approx\frac{1}{2}$ for the case of d- and g-wave resonances and
$\approx 1$ for s-wave resonance. Figs.6 and 7 demonstrate the
very good agreement of these formulas with the numerical
computations. Actually, in Fig.6(a) one can see in a broad range
of $k_0$ variation the very good agreement of the width $\Gamma_i^*(k_0)$,
extracted from numerically calculated $T(B,k_0)$, with the linear
functional dependence on the longitudinal momentum $k_0$, following from
Eq.(10). Also, the $\Gamma_i^*(\omega_{\perp})$ has inverse dependence on
$\omega_{\perp}$. It is confirmed by results of numerical calculations of
$\Gamma_i^*(\omega_{\perp})$ given in Fig.7. Thus, Eq.(10) can be
used for extracting important information about the analyzed
system. Indeed, by measuring the width $ \Gamma_{i}^*$ one can
extract from Eq.(10) the longitudinal momentum $k_0=\sqrt{2\mu
E_{\parallel}}/\hbar$ (longitudinal colliding energy
$E_{\parallel}$) and estimate the ``longitudinal" temperature of
the atomic cloud in the trap. This expression also shows that one
can control the width $\Gamma_{i}^*$ of the resonance by varying
the trap frequency $\omega_{\perp}$. Increasing of
$\omega_{\perp}$ leads to a narrowing of the resonance (see Figs.4
and 7), the effect of which can be used experimentally. Fig.6(b)
demonstrates the stability of the position of the CIR (the minimum
of the transmission coefficient $T$) with respect to $k_0$ variation and the linear growing of the
transmission coefficient $T$ with $k_0$ increasing. It is shown
that already at very low $k_0\sim 10^{-5}(\frac{1}{a_0})$ the
$T$-coefficient becomes large enough and consequently can be
experimentally ``visible" near CIR.
\begin{figure}[hbt]
\includegraphics[width=1.\columnwidth]{8a.eps}
\includegraphics[width=1.\columnwidth]{8b.eps}
\caption{(Color online) The relative populations $P_i/P_e$ (13) of the molecular
states $i=s,d$ and $g$  calculated as a function of $B$ near the
d-wave Feshbach resonance in Cs for the pair collisions in free
space (a) and in the harmonic waveguide with
$\omega_{\perp}=59.6$kHz (b). The s-wave scattering length $a(B)$
in free space from Fig.1 is also given.} \label{fig6}
\end{figure}
\begin{figure}[t]
\includegraphics[width=1.\columnwidth]{9.eps}
\caption{(Color online) Schematic picture (presented in \cite{Olshanii2,Melezhik11} for a single-channel potential model of
interatomic interaction) of
the bound state $E_b$ of two atoms as a function of $a_{\perp}/a$
in free space (solid line) and in the harmonic waveguide (dashed
line). The first excited resonance state $E_r$ is also presented.}
\label{fig7}
\end{figure}

To conclude this subsection we also give the expression for the position of the CIR with respect to the resonance at $a=0$ (to the point
$B_i^*$)
\begin{equation}
\Delta_{CIR}^{(i)}= \frac{\Delta_i}{1-C\gamma_i\frac{a_{bg}}{a_{\perp}}} \,\,.
\end{equation}
Since the resonance maximum at $B_i^*$ corresponds to the unit value
$T(B_i^*)\rightarrow 1$ and the minimum of the CIR approaches to
zero $T(B_i^*-\Delta_{CIR}^{(i)})\rightarrow 0$, we can define the ``right"
half width $\frac{1}{2}\Gamma_{CIR}^{(i)}$ of the CIR as the distance of the zero of the
coefficient $T$ to the right situated point where $T=\frac{1}{2}$. Using the formulas for
$\Delta_{CIR}^{(i)}$ and Eqs.(10) for $\Gamma_{i}^*$ we obtain
\begin{equation}
\frac{1}{2}\Gamma_{CIR}^{(i)}=\Delta_{CIR}^{(i)}-\frac{1}{2}\Gamma_i^*
\ _{\overrightarrow{(k_0\rightarrow 0)}}
\Delta_i[\frac{1}{1-C\gamma_i\frac{a_{bg}}{a_{\perp}}}-\frac{a_{\perp}^2
k_0}{2a_{bg}\gamma_i}]\,\,.
\end{equation}
The above expression describes the narrowing of the CIR width
$\Gamma_{CIR}^{(i)}$ with decreasing trap frequency
$\omega_{\perp}$, opposite to $\Gamma_i^*$ (see Fig.4).
This expression might also be of interest to the experimental analysis of the CIRs.

\subsection{Population of molecular states in free space and in the waveguide}

We have also calculated the relative populations $P_i/P_e$ of the
molecular states $|c,i\rangle$ in the process of pair atomic
collisions for harmonic traps of different frequencies
$\omega_{\perp}$ as well as in free space ($\omega_{\perp}=0$).
The populations $P_i$ are defined as
\begin{equation}
P_i=2\pi\int_0^{\infty}\int_0^{\pi}|\psi_{c,i}(r,\theta)|^2
 r^2 d r \sin\theta d \theta\,\,,
\end{equation}
where integration over the infinite region $0 \leq r <\infty$
gives a convergent results due to the decaying tails of the
molecular bound-state wave-functions $\psi_{c,i}({\bf
r})\rightarrow 0$ at distances of the order $\sim\overline{a}$ in
the closed channels $|c,i\rangle$. The population $P_e$ is defined
in a region near the origin $r\rightarrow 0$ of the molecular
dimension $\sim\overline{a}$  in the entrance channel by using
Eq.(13) where the upper limit of the integration over $r$ is
$\overline{a}$. The relative populations $P_i/P_e$ of the
molecular states for ultracold atomic collisions in free space and
in the harmonic waveguide with $\omega_{\perp}=59.6$kHz are shown
in Fig.8 in the region of B near the d-wave Feshbach resonance. In
free space the pair collision leads to a relatively low population
of the d-wave molecular state $P_2/P_e\sim 10^{-4}-10^{-5}$. The
populations of non-resonant s- and g-wave states are $P_1 \ll P_2$
and $P_3 \ll P_2$ essentially suppressed here. The dependence of
the d-wave molecular state population on $B$ repeats the
dependence on $B$ of the population of the region
$\sim\overline{a}$ in the entrance channel $P_2(B)/P_e(B)\sim
const$. However, in the confined geometry of the waveguide
($\omega_{\perp}\neq 0$) the pair collision leads to a resonant
enhancement of the relative population $P_2(B)/P_e(B)$ of the
d-wave molecular state near the point $B_2^*=47.944G$ where the
free space scattering length is  zero $a(B_{2}^*)=0$. In the
waveguide the point of appearance of the bound state is shifted
from the position defined by $1/|a|\rightarrow 0$ for free space
scattering (see point A in the illustrative scheme of the bound
and resonant states of an atomic pair given in Fig.9) to the point
$1/|a|\rightarrow \infty$ (the point B in Fig.9). This is why we
observe in Fig.8(b) the strong resonant enhancement of the
population $P_2(B_{2}^*)$ at the magnetic field $B_{2}^*=47.944$G.
The populations $P_1/P_e$ and $P_3/P_e$ of other molecular states
also show some enhancement with respect to the entrance channel
due to the coupling of the states with the ``resonance" channel
$|c,2\rangle$, which, however, are a few orders of magnitude less
than the enhancement of the population $P_2/P_e$. We do not
observe a resonant behavior of $P_2(B)$ at the point $B=47.57$G of
the CIR (see Fig.4 and the point C in Fig.9). We suspect that this
is the case due to the rather weak coupling in our model potential
$\hat{V}$(2) between the entrance channel $|e\rangle$ and the
resonant state in the closed channel $|c,2\rangle$ for stimulating
considerable transition to the molecular state $|c,2\rangle$ in
the closed excited channel of the waveguide.

In Figs.10  and 11 we present the results of our computation of
the molecular state populations near the s- and g-wave Feshbach
resonances -11.1G and 53.449G, respectively. We observe the
qualitatively analogous effect of resonant enhancement of the
relative populations $P_1(B)/P_e(B)$ and $P_3(B)/P_e(B)$ at the
points $B_{1}^*=18.1$G and $B_{3}^*=53.457$G in the harmonic
waveguide, which correspond to the positions of the zero of the
s-wave scattering length in the vicinity of s- and g-wave Feshbach
resonances, respectively. Here we notice $P_1(B)/P_e(B)\simeq const$, i.e.
the same dependence on $B$ of the s-wave molecular state population as
the population of the region $\leq\overline{a}$ in the entrance
channel $P_1(B)/P_e(B)\simeq const$ (see Fig.10(a)). In Fig.11(a) we
observe the resonant enhancement of the relative population
$P_3(B)/P_e(B)$ of the g-wave molecular state near the Feshbach
resonance in free space at $B_{0,3}=53.449G$ where $1/|a|\rightarrow
0$. A similar resonance enhancement near the
point $1/|a|\rightarrow 0$ of the appearance of the near threshold
resonance or weakly-bound state in free space (see point A in
Fig.9) is observable via the values $P_1(B)$,$P_2(B)$ and $P_e(B)$ but,
because $P_1(B)/P_e(B)\simeq const$ and $P_2(B)/P_e(B)\sim const$,
we do not observe this effect in the relative populations
$P_1(B)/P_e(B)$ and $P_2(B)/P_e(B)$ in Figs.10(a) and 8(a).

\begin{figure}[hbt]
\includegraphics[width=1.\columnwidth]{10a.eps}
\includegraphics[width=1.\columnwidth]{10b.eps}
\caption{(Color online) The relative populations $P_i/P_e$ (13) of the molecular
states $i=s,d$ and $g$  calculated as functions of $B$ near the
s-wave Feshbach resonance for Cs for the pair collisions in
free space (a) and in the harmonic waveguide with
$\omega_{\perp}=59.6$kHz (b). The s-wave scattering length $a(B)$
from Fig.1 is also provided.} \label{fig8}
\end{figure}

\begin{figure}[hbt]
\includegraphics[width=1.\columnwidth]{11a.eps}
\includegraphics[width=1.\columnwidth]{11b.eps}
\caption{(Color online) The relative populations $P_i/P_e$ (13) of the molecular
states $i=s,d$ and $g$  calculated as functions of $B$ near the
g-wave Feshbach resonance in Cs for the pair collisions in
free space (a) and in the harmonic waveguide with
$\omega_{\perp}=59.6$kHz (b). The s-wave scattering length $a(B)$
from Fig.1 is also presented.} \label{fig9}
\end{figure}

\section{CONCLUSION}
We have developed a theoretical model for a
quantitative analysis of the Feshbach resonance shift and width
induced by an atomic waveguide. It is based on our multichannel
approach for confinement-induced resonances and atomic
transitions in the waveguides in the multi-mode regime
\cite{Saeidian}. In this scheme the single-channel
(scalar) interatomic interaction is replaced by a four-channel (tensorial)
potential modeling resonances of different structure according
to the two-channel parametrization of A.D.Lange et. al.
\cite{Lange}. The experimentally known parameters of
the Feshbach resonance in the absence of the waveguide are used
as an input in our approach.
We have calculated the shifts and widths of s-, d- and g-wave
magnetic Feshbach resonances of Cs atoms emerging in harmonic
waveguides as CIRs and resonant enhancement of the transmission at zeros of the free space scattering length.

In particular we find that the relationship $a=a_{\perp}/C$ for the position of
the CIR in a harmonic waveguide is fulfilled with high accuracy
for the Feshbach resonances of different tensorial structure which
holds in spite of the fact that this property was originally obtained in the framework of a s-wave
single-channel pseudopotential approach \cite{Olshanii1}. Note,
that this property was experimentally confirmed for d-wave Feshbach
resonances in a gas of Cs atoms \cite{Haller2010}. The maximum
of the transmission, corresponding to the zero of the scattering
amplitude, is shown to be independent of the trap field strength
and then again corresponds to the zero $B_{i}^*$ of the s-wave scattering length $a(B)$
in free space. In a nutshell, the Feshbach resonance in
free space develops in the harmonic waveguide into a minimum of $T$ (position
of CIR), defined by the formulas $a=a_{\perp}/C$, and a maximum,
coinciding with the position of zero of the s-wave scattering
length a. The 'distance' between these extrema is equal to $\Delta_{CIR}^{(i)} =\Delta_i/(1-C\gamma_i
a_{bg}/a_{\perp})$.

We have derived expressions for the widths $ \Gamma_{i}^*,\frac{1}{2}\Gamma_{CIR}^{(i)}$ of the resonant enhancement
of $T$ at $B_i^*$ and the ``right" side half-width of the minimum of the
$T$-coefficient i.e. at the position of CIR and confirmed its
validity by numerical results for $k_0\rightarrow 0$. By
measuring the width $ \Gamma_{i}^*$ one can, in principle, extract from these expressions
the longitudinal collision energy and estimate the ``longitudinal" temperature of
the atomic cloud in the trap. In other words the
width of the atomic loss resonance observed in the experiment
\cite{Haller2010} at the point of CIR might contain important
information about the longitudinal atomic momentum $k_0$ and the
temperature of the gas. It also shows that one can control the
width $\Gamma_{i}^*$ of the resonance at $a(B_i^*)=0$ by varying the trap
frequency $\omega_{\perp}$. An increase of $\omega_{\perp}$ leads
to a narrowing of the resonance, an effect which could potentially
be used experimentally.

Finally, the molecule formation rates in a waveguide show an enhancement
for the case of a corresponding zero of the s-wave scattering
length $a(B_{i}^*)=0$. We have shown that the positions of these
resonances are stable w.r.t. the variation of the confining frequency
$\omega_{\perp}$ of the waveguide.

Our model adds to the possible studies of
scattering processes of ultracold atomic gases in waveguides
beyond the framework of s-wave resonant scattering. Our model
might be extended to the cases of fermions or distinguishable atom
scattering, including transverse excitation/diexcitation processes \cite{Saeidian}.
It permits also for the investigation of other trap geometries
\cite{MelHu} and more realistic interatomic interactions.

\section{ACKNOWLEDGEMENTS}
We thank E. Haller and H.-C. N\"agerl for fruitful discussions. Authors
V.S.M. and P.S. acknowledge financial support by the Deutsche
Forschungsgemeinschaft and the Heisenberg-Landau Program. V.S.M.
thanks the Zentrum f\"ur Optische Quantentechnologien of the
University of Hamburg and Sh.S. thanks the Bogoliubov Laboratory of
Theoretical Physics of JINR at Dubna for their warm hospitality.


\section{Introduction}

The study of the fractional quantum Hall effect (FQHE) of electrons in semi-conductor materials \cite{prangebook} and, potentially, of bosons in rotating Bose gasses \cite{cooper2008} or in artificial gauge fields \cite{cooper2011} has been a fertile ground for the conception of topological phases of matter, which exhibit exotic and exciting theoretical properties. \cite{nayak2008}  In this paper we shall discuss the construction of a FQHE wave function exhibiting a ``multicomponent'' degree of freedom such as a spin, \cite{du1995} valley, \cite{bishop2007,dean2011} isospin, \cite{eng2007} layer, \cite{eisensteinbook} or subband. \cite{hormozi2012,palmer2006} The wave function describes the low-energy behaviour of a system possessing a vanishing energy gap to quasihole excitations in the thermodynamic limit, and is thus not itself a candidate to describe a topological phase. Nevertheless the study of this wave function is an interesting exploration that enhances the understanding of the construction of ``multicomponent'' FQHE wave functions in general.

An ideal theoretical description of a FQHE state comprises at least the following three ingredients: A local Hamiltonian which describes the ground state and excitation spectrum of the Hilbert space; relatively simple analytic wave functions which are the highest density zero energy states corresponding to the Hilbert space of that Hamiltonian; and a (rational) two-dimensional (2D) conformal field theory (CFT) which generates these wave functions. If a plasma analogy is available, then a gapped state is associated with the analogous plasma being in a screening phase. \cite{laughlin1983a} In general, the quantum Hall Hilbert space is built from a basis of monomials in the complex particle coordinates $z_i$.  In certain instances---the Laughlin series, the Read--Rezayi series, the Gaffnian and Haffnian wave functions, the Halperin wave functions and the non-Abelian spin-singlet (NASS) states \cite{laughlin1983a,read1999, simon2007b,green2002,halperin1983,ardonne1999,ardonne2001,ardonne2002} ---it has been possible to i) determine simple analytical expressions for the polynomial wave functions, and ii) construct a  local Hamiltonian whose zero energy eigenstates are in one-to-one correspondence with those polynomial wave functions. A key feature in each of these special cases is that the wave functions are uniquely defined by their vanishing properties.

Changing wave functions describing a FQHE state in a seemingly innocuous way, may in fact have severe consequences. The wave function of the Moore--Read state for spinless bosons vanishes quadratically when three of the constituent bosons are coincident. The Gaffnian wave function is obtained by changing this behaviour, such that the wave function vanishes as a third power, when three constituent particles are coincident. \cite{simon2007b} This simple change, however, results in a nonunitary, compressible wave function, which does not describe a topological phase. Nevertheless nonunitary wave functions such as the Gaffnian are still of interest as they are thought to correspond to critical points between other unitary, incompressible topological phases. This scenario is well understood in the case of the Haldane--Rezayi\cite{haldane1988} wave function, which describes the phase transition between a $d$-wave spin-singlet phase and a strongly paired state. \cite{read2000,green2002} Similar scenarios have been suggested for the Gaffnian.\cite{simon2007b} It is noteworthy that the Gaffnian wave function has large overlap with an incompressible composite fermion state thus suggesting that the Gaffnian is a critical point between the composite fermion phase and some other phase. \cite{jain1989,simon2007b,regnault2008,regnault2009}

In this paper, we present evidence for a quantum Hall wave function, which inherits properties from both the NASS state---a unitary, spin-singlet quantum Hall state---and the nonunitary Gaffnian. Such a wave function has also been considered in Ref.~\onlinecite{estienne2011}. Most notably, we argue that the wave function can be written as a special polynomial and that it also corresponds to the highest density zero energy state of a certain local Hamiltonian. We christen it the ``spin-singlet Gaffnian'' (\wf{}). (Throughout this work, we shall describe how the wave function corresponds to a ``spin'' degree of freedom, however, note that all of what follows applies equally well to any other type of multicomponent degree of freedom such as valley or layer index.)

\subsection*{Statement of results}

Our first key result is a model local Hamiltonian describing the \wf{} wave function. In Sec.~\ref{secPseudopotential} we shall explain how it can be written in terms of generalized Haldane pseudopotentials. \cite{haldane1983,simon2007a,davenport2012a} In Secs.~\ref{secSqueezing} and \ref{secPauli}, we give an overview of how both the squeezing algorithm \cite{bernevig2008b,ardonne2011} and the generalized Pauli principle for spinful states \cite{estienne2011}  apply to the \wf{} wave function.  In Sec.~\ref{secNumerics} we shall present the results of numerical exact diagonalization of the \wf{} Hamiltonian. We have determined that the counting of zero energy eigenstates in the quasihole spectrum and also the entanglement spectrum \cite{li2008,sterdyniak2011} matches the counting predicted by the squeezing algorithm. The counting of zero energy eigenstates is also found to be identical to the result generated by the spinful version of the generalized Pauli principle.

Our second key result is an analytic form for the ground-state wave function. The \wf{} wave functions are constructed from conformal blocks in a CFT associated with the semidirect product of nonunitary minimal models expressible as $\mathcal{M}(3,5) \ltimes \mathcal{M}(5,7)$ (see also Ref. \onlinecite{estienne2011}). We shall discuss the CFT in Sec.~\ref{secCft}. In a Bosonic incarnation, the \wf{} ground state occurs at filling factor $\nu=4/5$ (and on the sphere, it has a shift $\delta=3$). The proposed analytic construction for the \wf{} is consistent with the vanishing properties required by the CFT and the ground state of the Hilbert space generated by squeezing.  We have determined that our proposed ground-state wave function is the unique, highest density zero energy ground state of our model Hamiltonian at $\nu=4/5$ and $\delta=3$. We shall discuss an explicit form of the ground-state wave function in Sec.~\ref{secwave function}.

\section{Model Hamiltonian and the Characterization of its Zero Modes}
\label{secMethods}

In this section we shall describe in detail three independent methods by which we are able to study the \wf{} wave function. The cornerstone of our argument lies with the proposal of a model local Hamiltonian. First, we shall explain how this local Hamiltonian is constructed and then we shall present two alternative methods to count the number of zero energy modes of the proposed Hamiltonian. The first method employs the spin-dressed squeezing algorithm (which can also be used to obtain the wave functions) and the second method employs the generalized Pauli principle. In Sec.~\ref{secNumerics} we shall present numerical evidence demonstrating that the zero modes generated by diagonalizing the local Hamiltonian presented here are in precise agreement with the zero modes generated by the squeezing algorithm and the counting obtained from the generalized Pauli principle.

\subsection{Pseudopotential construction of Hamiltonian}
\label{secPseudopotential}

In prior investigations, most notably for the Laughlin and Moore--Read wave functions, it has been determined that trial quantum Hall wave functions correspond to unique, zero energy ground states of certain model Hamiltonians. \cite{haldane1983,greiter1992} Often, these model Hamiltonians are most simply expressible in terms of Haldane pseudopotentials \cite{haldane1983} and their generalizations. \cite{simon2007a,davenport2012a} Given a system with a certain $M$-body interaction potential $V(z_1,\ldots, z_M)$, the action of a pseudopotential is to project out a particular component of that interaction.  Those components are labelled by a convenient set of quantum numbers, which describe all possible few-particle interactions. The vector space of few-particle interactions is spanned by specifying both the relative angular momentum $L$ and, if we consider a spin degree of freedom as well, the spin quantum number $S$ also. (Note that in general it is necessary to further distinguish between distinct interaction components with the same $L,S$. In other words, there exists in general a subvector space of dimension $d_{L,S}$ for each $L$ and $S$ sector. \cite{davenport2012a,simon2010} In this paper, however, we shall only be concerned with subspaces of dimension $d_{L,S}$=1 or 0, and so, for clarity, we omit any additional notation.) Pseudopotentials specify projection operators in the Hamiltonian and are thus labelled by two distinct sets of such quantum numbers. The pseudopotential $V^M_{L,S ; L' S' }$ is expressed in terms of the vector space $\left|  L, S \right \rangle$ and the $M$-body interaction potential $V(z_1,\ldots, z_M)$ as
\begin{equation}
\label{eqMulltiParticlePseudopotential}
V^M_{L,S ; L' S' } = \left\langle L, S \right |      V(z_1,\ldots, z_M)        \left|  L', S' \right \rangle.
\end{equation}
The general form of the Hamiltonian is then given by
\begin{equation}
\label{eqGeneralHamiltonian}
H=\sum_{L, S, L' S'}    \left|  L, S \right \rangle  V^M_{L,S ; L'S' } \left\langle L', S' \right |.
\end{equation}
Such a basis of pseudopotentials is particularly convenient when the interaction potential is rotationally and/or spin rotationally invariant, since in that case the pseudopotentials are diagonal in $L$, $S$, or both.

At our convenience, we can pick certain special many-body interactions (such as $\delta$-function-type interactions) for which only a small set of pseudopotentials remain nonzero. The impact of specifying a positive value of a given pseudopotential in a model Hamiltonian is to assign energy to the corresponding component of the interaction, therefore, if such a component is present in a given trial wave function, then that wave function will not be a zero energy state of our Hamiltonian. Conversely, if a component is not present in a given trial wave function, then we can include the corresponding pseudopotential in the Hamiltonian without introducing an extra energy cost. In this way we can tailor the Hamiltonian to correspond to the desired properties of a given trial ground-state wave function (these properties might come, for instance, from a CFT description of the state; see Sec. \ref{secCft}).

To give an example, briefly, it is known that the Moore--Read wave function for bosons at $\nu=1$ is the unique, highest density zero energy ground state of a spin-polarized three-body contact interaction. \cite{greiter1992} The space of pseudopotentials for spin-polarized three-body interactions is spanned by a relative angular momentum $L$. $S$ takes only its maximal value for a three-body interaction, $S=3/2$. \cite{simon2007a} (Note a slight abuse of notation here: We are denoting bosons as if they are spin-1/2 objects, e.g., three bosons can have $S=3/2$ and $S=1/2$ interaction channels. In actual fact we are describing systems such as cold atomic gasses for which bosons can be engineered to have internal two-state degrees of freedom, and we have simply mapped that two-state degree of freedom onto a spin degree of freedom.) In the language of pseudopotentials,  the Moore--Read state corresponds to $ V^3_{0,3/2; 0,3/2}$ being positive and all other pseudopotentials being zero and, in addition, the corresponding Hilbert space is restricted to only spin-polarized sectors. Along similar lines, the spin-polarized Gaffnian wave function for bosons corresponds to  $ V^3_{0 ,3/2; 0,3/2}$ and  $ V^3_{2,3/2 ;2,3/2}$ being positive and all other pseudopotentials being zero, and again  the Hilbert space is restricted to spin-polarized sectors (note that for $L=1$, $S=3/2$, and $M=3$, we have a 0-dimensional vector space of wave functions, so no corresponding $L=1$, $S=3/2$, and $M=3$ pseudopotential can occur). In these examples we observe that the pseudopotentials are in fact diagonal in the $L,S$ sectors. This is a direct consequence of the interaction being, respectively, rotationally and spin-rotationally invariant. This feature will remain present in our application of pseudopotentials to the \wf{} wave function, and so at this point we shall drop the repeated indices and denote diagonal pseudopotentials by $V^M_{L,S} \equiv V^M_{L,S ; L S } $.

The Bosonic NASS state, like the Moore--Read state, corresponds to a spin-polarized three-body contact interaction, however the Hilbert space now includes additional spin sectors (not just the spin-polarized sector).  In the language of pseudopotentials,  the Bosonic NASS state is the highest density zero energy ground state of a Hamiltonian with the positive, diagonal pseudopotential $ V^3_{0,3/2}$. Motivated by the generalization of the Moore--Read Hamiltonian to the Gaffnian Hamiltonian, our proposal for the \wf{} Hamiltonian is to keep $ V^3_{0,3/2}$, $ V^3_{2,3/2}$, and $ V^3_{1,1/2}$ positive. The proposed Hamiltonian for the \wf{} wave function is thus expressed as
\begin{align}
\label{eqHamiltonian}
H_{\mbox{\wf{}}} &= \left|  0, 3/2 \right \rangle  V^3_{0,3/2}  \left\langle 0, 3/2 \right | \nonumber \\ &+ \left|  2, 3/2 \right \rangle  V^3_{2,3/2}  \left\langle 2, 3/2 \right | \nonumber \\&+ \left|  1, 1/2 \right \rangle  V^3_{1,1/2}  \left\langle 1,  1/2 \right |.
\end{align}

From this pseudopotential Hamiltonian, we can already infer some properties of its zero energy eigenstates. Because the three-body interaction in the $S=3/2$ channel is identical to the Hamiltonian generating the (polarized) Gaffnian, the zero energy eigenstates will vanish as at least a third power when three particles with the same spin are brought to the same point. In addition, when three particles have overall spin $S=1/2$, the wave functions vanish at least quadratically. These are indeed the vanishing properties consistent with the CFT description of the \wf{}, as we shall describe in Sec.~\ref{secCft}. It is worth stating here that the \wf{} wave function likely also corresponds to the ground state of other, more complicated, local Hamiltonians involving $M$-body terms with $M>3$, following the line of reasoning discussed in Ref.~\onlinecite{simon2007a}.

\subsection{Spin-dressed squeezing algorithm}
\label{secSqueezing}

There is another way to characterize the spin-singlet Gaffnian, apart from by means of the Hamiltonian we introduced above. This method is the so-called ``squeezing method''. The idea underlying the method is that for many model states, the wave functions have a large number of zero coefficients if expressed in terms of the space of all possible monomials. These zeros are closely related to the vanishing properties of the wave functions. The method was pioneered in Ref.~ \onlinecite{bernevig2008b}. The spinful case, which we will employ here, was described in great detail in Ref.~\onlinecite{ardonne2011}.

We consider quantum Hall states on the sphere, \cite{haldane1983} in the presence of $N_\Phi$ flux quanta. In the lowest Landau level, this construction gives rise to $N_\Phi+1$ orbitals, whose occupation numbers will be denoted by $(n_0,n_1,\ldots,n_{N_\Phi})$. These orbitals have angular momentum $(N_\Phi/2, N_\Phi/2-1,\ldots,-N_\Phi/2)$. We will use the orbital occupation numbers to label the states in the Hilbert space. Because total angular momentum is a good quantum number, we can split the Hilbert space into sectors with different values of total $L_z$, and we will always consider the number of particles $N$ to be fixed. For spinless fermions, the Pauli principle specifies that the occupation of any orbital $n_i$ can be 0 or 1, while for bosons, there is no constraint.

To obtain all states in a particular $L_z$ sector, we divide the $N$ particles over the orbitals, such that one obtains the correct value of total $L_z$, and such that the particles are ``desqueezed'' as much as possible. This means that for an even number of bosons in the $L_z=0$ sector, one considers the occupation $(N/2,0,\ldots,0,N/2)$. To obtain the other states in this sector, one generates all possible pairwise rearrangements of particles occupying orbitals, keeping the total $L_z$ fixed. In the general case, this entails $(\ldots,n_{i},n_{i+1},\ldots,n_{j-1},n_{j},\ldots)$ transforming into   $(\ldots,n_{i}-1,n_{i+1}+1,\ldots,n_{j-1}+1,n_{j}-1,\ldots)$, where all the other occupation numbers remain unchanged. In this way, one obtains the {\em full} Hilbert space in each total $L_z$ sector.

As we alluded to above, the wave functions of many model states have a large number of zero components, if expressed in the Hilbert space described above. For instance, if one considers the $\nu=1/2$ Bosonic Laughlin state, it suffices to construct a {\em reduced} Hilbert space, by starting to squeeze from the following so-called {\em root configuration}, $(1,0,1,0,\ldots,0,1,0,1)$, instead of the completely desqueezed configuration $(N/2,0,\ldots,0,N/2)$. All states in the {\em full} Hilbert space which do not appear in the {\em reduced} Hilbert space have zero coefficient in the Laughlin state. Moreover, it turns out that there is only one $L=0$ state one can construct in the {\em reduced} Hilbert space. Thus, to obtain the Laughlin state for a certain number of particles, one constructs the reduced Hilbert space from the root configuration $(1,0,1,0,\ldots,0,1,0,1)$, and demands that $L^{+}$ on a general state in this reduced Hilbert space gives zero. This procedure gives a set of equations for the coefficients, whose solution gives the Laughlin state, expressed in the monomial basis on the sphere. In a similar way, many other model states can be obtained. For instance, to obtain the level-$k$ Read-Rezayi states \cite{read1999}, one only has to change the root configuration to $(k,0,k,\ldots,k,0,k)$. Interestingly, the root configurations correspond to the ``thin-torus'' limit of the states. \cite{bergholtz2005,bergholtz2006,seidel2006}

To describe spinful wave functions, the squeezing method was augmented in Ref.~\onlinecite{ardonne2011}. The idea is to start with a {\rm root configuration}, for which one at first ignores the spin degrees of freedom. To construct the reduced Hilbert space, one creates a set of orbital occupations, by squeezing in all possible ways. Finally, one assigns spin to all the particles, in all possible ways. We will be interested in spin-singlet states, for which total $S$ is a good quantum number. So to construct the ground state, we work in the total $S_z=0$ sector. After distributing the spin in all possible ways for each configuration, we have constructed the reduced Hilbert space.  The spinful orbital occupations are now denoted by $(n_{0,\uparrow},n_{0,\downarrow},n_{1,\uparrow},n_{1,\downarrow},\ldots,n_{N_{\Phi},\uparrow},n_{N_{\Phi},\downarrow})$ .

We still have to specify how to obtain the correct wave function in the reduced Hilbert space. In general, we will be interested in ground-state wave functions which have $L=S=0$, so to obtain those, we work in the $L_z=S_z=0$ subsectors, and demand that the action of $L^{+}$ and $S^{+}$ on the state gives zero. It turns out that this is in general not enough to completely specify the wave functions. In addition, one has to demand that certain states in the reduced Hilbert space have zero coefficient. The SSG is an example where this happens, as we shall presently describe.

Let us now discuss how the squeezing algorithm applies to the construction of the \wf{} Hilbert space. We have already seen, in the construction of the Hamiltonian in Eq.~(\ref{eqHamiltonian}), that the zero energy eigenstates of the \wf{} wave function vanish as at least a third power when three particles with the same spin are brought to the same point and, in addition, when three particles have overall spin $S=1/2$, the wave functions vanish at least quadratically. These observations motivate the use of the root configurations based on the pattern $(2,0,2,0,0,2,0,2,0,0,\ldots,0,0,2,0,2) $, because like the Moore--Read state, the resulting wave function will vanish at least quadratically when three particles are at the same location. These root configurations have the property that two neighbouring orbitals can maximally be occupied by two particles, while five consecutive orbitals can be occupied by at most four particles. We note in advance that these root configurations satisfy the spinful generalized Pauli principle, which we discuss in the next subsection.

To construct the ground-state wave function, we start with the root configuration described above, and then construct the reduced Hilbert space. To obtain the ground state, we then demand that both $L^{+}$ and $S^{+}$ act to give zero on all states in the reduced Hilbert. This construction, however, only guarantees that the obtained wave function vanishes as a second power when three particles of the same spin are coincident; but we really wanted the wave function to vanish as a third power in that case. We therefore enforce an additional constraint: That all basis states in the reduced Hilbert space have zero coefficient, on the condition that either one has for spin-up particles $n_{0,\uparrow} = 2$ and $n_{2,\uparrow} > 0$, or similarly, for spin-down particles, one has $n_{0,\downarrow} = 2$ and $n_{2,\downarrow} > 0$. With that additional constraint, the wave function now vanishes as a third power when three particles are coincident with total $S=3/2$.

\subsection{Spinful generalized Pauli principle}
\label{secPauli}

A further method to characterize the zero energy  space of the SSG Hamiltonian is to use the generalized Pauli principle for quantum Hall wave functions. \cite{bernevig2008a} The generalized Pauli principle  is most readily expressed in terms of partitions (or occupied orbitals).  A partition ${\lambda}$ is defined to be an ordered set of $N$ integers, \{$\lambda_1,\ldots,\lambda_N$\}, where $\lambda_i\geq \lambda_{i+1}$.

To describe the Hilbert space of spinless particles in flux $N_\Phi$, we restrict the integers to the set $\lambda_i \in \{0,1,\ldots,N_{\Phi} \}$. We then have the following relation between the orbital occupation numbers $n_{0},n_{1},\ldots,n_{N_\Phi}$, which were introduced in the previous subsection, and the
$\lambda_i$ forming the partition $\lambda$: Namely, $n_{N_\Phi}$ is the number of $i$ such that $\lambda_{i} = N_\phi$ etc.  Thus in general,
$n_{j}$ is the number of $i$ such that $\lambda_i=j$.

To characterize the quasiholes of the $(k,r)$ clustered states\cite{bernevig2008c} we introduce the notion of $(k,r)$-admissible partitions, which are partitions obeying the following condition for all $i$:
\begin{equation}
 \lambda_i - \lambda_{i+k} \geq r.  \label{PauliSpinless}
\end{equation}

For given fixed values of $N$ and $N_\Phi$, one identifies the number of admissible partitions with the number of quasihole states in the corresponding spectrum. For instance, the $r=2$ series corresponds to the level-$k$ Read--Rezayi states, the Laughlin state $(k=1)$, and the Moore--Read state $(k=2)$ being the simplest cases. The Gaffnian wave function is associated with the $(k=2,r=3)$ generalized Pauli principle.

In order to generalize this method to the spinful case, we follow the argument presented in Ref.~\onlinecite{estienne2011}. First, the partition is replaced by a spinful partition that mixes momentum and spin. A spinful partition $(\lambda,\sigma)$ is specified by $N$ integers, \{$\lambda_1,\ldots,\lambda_N$\}, and now in addition, a set of spin indices  $\{\sigma_1, \ldots \sigma_N\}$, where $\sigma_i \in \{-1, 1\} $. If either the condition
$\lambda_i>\lambda_{i+1}$ or $\lambda_i=\lambda_{i+1}$ and $\sigma_i \geq \sigma_{i+1}$ holds, we say that $(\lambda,\sigma)$ constitutes a spinful partition.

The spinful generalization of the $(k,r)$ admissible partitions is given by the following conditions, which have to hold for all $i$:
\begin{eqnarray}
&& \lambda_i - \lambda_{i+k} \geq r \nonumber,\\
&\mathrm{~or~}& \lambda_i - \lambda_{i+k} = r - 1 \mathrm{~and~} \sigma_i < \sigma_{i+k}\label{Paulispinful}.
\end{eqnarray}

In Ref.~\onlinecite{estienne2011} it is shown that the densest Bosonic state that leads to a spinful $(k,r)$ admissible possible corresponds to a filling factor $\nu = 2k/(2r-1)$ and a shift $\delta = r$. Note that on the sphere, the relation between the filling factor $\nu$ and shift $\delta$ is
\[
\nu=\frac{N}{N_\Phi+\delta}.
\]

It has been proposed \cite{estienne2011} that nonsymmetric Jack polynomials can represent spin-singlet states that are the spin generalization of the spinless clustered states. In that case, the root partition is chosen to be a spinful $(k,r)$ admissible partition. The number of quasihole states can be obtained in a similar way as the spinless case, counting the number of admissible partitions. Indeed, these properties have been checked for the Halperin spin-singlet states $(k=1,r)$ and the NASS state $(k=2,r=2)$. In our application of the spinful generalized Pauli principle to the \wf{}, we employ the $(k=2,r=3)$ Pauli principle (as for the spin-polarized Gaffnian), but now using the spinful partitions.

\section{Numerical results}
\label{secNumerics}

In this section we shall present numerical evidence demonstrating that the zero energy states of the local Hamiltonian presented in Sec.~\ref{secMethods}, Eq. (\ref{eqHamiltonian}), are correctly reproduced by both the squeezing algorithm and by the generalized Pauli principle which were also described in Sec.~\ref{secMethods}. We shall present evidence derived from both the quasihole spectrum and particle entanglement spectrum.

\subsection{Quasihole spectrum}
\label{secSpectrum}

The quasihole spectrum is determined by numerical exact diagonalization of the Hamiltonian for finite-sized systems in the sphere geometry. This is done for a variety of system sizes, $N$, and for a variety of fluxes $N_\Phi$. Eigenstates are labelled by the quantum numbers $L_z$ and $S_z$ and fall into $(L,S)$ multiplets. (Note that these quantum numbers are completely separate from the $L$ and $S$ describing the pseudopotentials unless one has only $M$ particles with an $M$-body interaction.) Given constraints on the dimension of the spinful Hilbert space, we were able to study systems of up to $N=12$ and $N_\Phi=12$. An example of such a spectrum is plotted in Fig.~\ref{figSpectrum}.

For comparison, we generated the same Hilbert space via the spin-dressed squeezing and generalized Pauli algorithms. For the spin-dressed squeezing algorithm we started from a root partition $(2,0,2,0,0)$, e.g.,  $(2,0,2,0,0,2,0,2,0,0,2,0,2)$ for $N=12$, and we applied the procedure described in Sec.~\ref{secSqueezing} to generate the Hilbert space. The ground-state wave functions, corresponding to the densest root configurations with $L_z=0$, occur at filling $\nu=4/5$ and shift $\delta=3$. For the generalized Pauli algorithm, we specified a spinful admissible partition $(k=2,r=3)$.

Our key observations are as follows: First, in the spectrum of the Hamiltonian there is a unique zero energy state occurring only in the $L=0$, $S=0$ sector for $N=4$ at $N_\Phi=2$, for $N=8$ at $N_\Phi=7$ (see Fig.~\ref{figSpectrum}) and for $N=12$ at $N_\Phi=12$, which all correspond to $\nu=4/5$ and $\delta=3$, consistent with both the squeezing and generalized Pauli approaches; second, for the quasihole spectrum generated for $N$ between 3 and 12 and for $N_\Phi$ up to 20, we have checked that the counting of zero energy states in each $L_z$ and $S_z$ sector precisely matches the counting predicted by both the spin-dressed squeezing algorithm and the generalized Pauli principle.

We have further determined for $N=4$ and $N=8$ that the ground-state monomial expansion of the Hilbert space generated by exact diagonalization of Eq.~(\ref{eqHamiltonian}) precisely matches the ground-state monomial expansion generated by the squeezing algorithm as described in Sec.~\ref{secSqueezing}. In addition we have found that the ground-state monomial expansion also precisely matches that of the proposed analytic form of the ground-state \wf{} wave function (we shall discuss the analytical form of the ground-state wave function in Sec.~\ref{secwave function}).

\begin{figure}[t,b]
\includegraphics[width=0.5\textwidth]{spinGaffSpectrum_n_8_flux_7_s_0.eps}
\caption{(Colour online) The $S=0$ sector of the quasihole spectrum of the Hamiltonian $H_{\mbox{\wf{}}} $, defined in Eq.~(\ref{eqHamiltonian}), for $N=8$ and $N_\Phi=7$, obtained by numerical exact diagonalization in the sphere geometry. The inset zooms in on the bottom left corner of the spectrum and shows the unique zero energy state located in the $L=0$, $S=0$ sector.}
\label{figSpectrum}
\end{figure}

\subsection{Particle entanglement spectrum}
\label{secPES}

The so-called particle entanglement spectrum (PES) is determined from the reduced density matrix of a subsystem that results from the partition of the whole system into two (or more) parts $A$ and $B$. \cite{li2008, sterdyniak2011} During this partitioning we also keep the overall geometry unchanged. The reduced density matrix $\rho_A$ is given in terms of the full density matrix $\rho=\sum \left | \psi \right \rangle \left \langle \psi \right |$ by tracing out the $N_B$ particles in the $B$ partition: $\rho_A = \mbox{Tr}_B (\rho) $. The PES arises from diagonalizing $\rho_A$ and then classifying the resulting eigenstates according to the symmetries of the problem, in our case by $L_A$ and $S_A$.

For model states for topological phases such as the Laughlin and Moore--Read states, it has been observed that there is a characteristic PES: The number of nonzero eigenvalues for $\rho_A$ is identical to the number of quasihole states for a system with identical geometry but only $N_A$ particles.\cite{ li2008, sterdyniak2011} The number is usually exponentially lower than the dimension of $\rho_A$, or equivalently, there is an infinite so-called ``entanglement gap'' to the remaining eigenvalues. \cite{li2008,zozuly2009,thomale2010,hermanns2011}

For the \wf{} wave function, our key observations are as follows: We find evidence of an infinite entanglement gap for $N=4$ at $N_\Phi=2$, for $N=8$ at $N_\Phi=7$ and for $N=12$ at $N_\Phi=12$ (see Fig.~\ref{figPesSpectrum}), which we associate with our model ground-state wave function for the \wf{}; and we find that the counting predicted for the quasihole excitations up to $N=12$, obtained by considering the PES for all possible partitions of the system into $A$ and $B$ subsystems, precisely matches the counting predicted by the generalized Pauli principle in each $S$ and $L$ sector.

With the evidence presented in this section in mind, we conclude that the proposed \wf{} Hamiltonian in Eq.~(\ref{eqHamiltonian}) generates the correct zero energy Hilbert space for the \wf{} wave function, at least for system sizes up to $N=12$. Based on this evidence, we also expect our conclusion to hold for any other system size. Thus we have a Hamiltonian which describes the ground state and excitation spectrum of the Hilbert space.

An important question, left unanswered by our numerical calculations here, is what is the size of the gap for quasiparticle excitations in the thermodynamic limit? For finite-sized systems, looking at Fig.~\ref{figSpectrum}, the gap is clearly finite (as it is in the Gaffnian case). There is good reason, however, to argue that the gap will in fact vanish in the thermodynamic limit, and hence the \wf{} will be gapless and the wave function therefore compressible (see Ref.~\onlinecite{read2009}). That argument is based on the fact that the CFT corresponding to the \wf{} is nonunitary, as we shall presently discuss.

\begin{figure}[t,b]
\includegraphics[width=0.5\textwidth]{spinGaffPes_n_12_flux_12_s_0_na_6.eps}
\caption{(Colour online) The particle entanglement spectrum (PES) for $N=12$, $N_\Phi=12$ and $N_A=6$, obtained by numerical exact diagonalization in the sphere geometry. The negative log of the eigenvalues, $\xi$, of the reduced density matrix $\rho_A$ (defined in the text) are plotted.}
\label{figPesSpectrum}
\end{figure}

\section{Conformal Field Theory}
\label{secCft}

Let us now discuss the guiding principles of the derivation of the CFT for the \wf{} wave function (see also related considerations in Ref. \onlinecite{estienne2011}). We shall expand further on the discussion in the Appendix. For an introduction on CFT, we refer to the book, Ref. ~\onlinecite{byb}, and the seminal paper by Belavin, Polyakov, and Zamolodchikov (BPZ), Ref~\onlinecite{belavin1984}.

\subsection{Coset constructions for minimal models}

The minimal models introduced by Belavin, Polyakov, and Zamolodchikov can be written in terms of coset models of the $su(2)_k$ WZW model. \cite{goddard1985} In the unitary case, one has
\begin{equation}
\mathcal{M} (k+1,k+2) = \frac{su(2)_1\times su(2)_{k-1}}{su(2)_k} \ .
\end{equation}
For $k=2$, this gives the Ising CFT, $\mathcal{M} (3,4)$. For integer $k$, the coset theories are unitary, a property which is inherited from the $su(2)_k$ WZW model, which is unitary for $k$ integer.

General minimal models are labeled $\mathcal{M}(p',p)$ for arbitrary (non-negative) co-prime integers $(p',p)$. The minimal model is unitary only if $| p - p'| = 1$, otherwise it is nonunitary. Nevertheless, there exists a coset description of nonunitary minimal models in terms of fractional level WZW models. \cite{mathieu1990,mathieu1992} In particular, the parameter $k$ is given by $k=\frac{3p'-2p}{p-p'}$. The nonunitary minimal model $\mathcal{M}(3,5)$, featuring in the CFT description of the Gaffnian wave function, corresponds to $k=-1/2$, while the Yang-Lee model $\mathcal{M}(2,5)$ has $k=-4/3$.

The central charge of the minimal models is given by $c(p',p) = 1-\frac{6(p'-p)^2}{pp'}$. The primary fields of the minimal models are labeled by integers $(r,s)$, which take the values $1\leq r < p'$ and $1\leq s < p$. The labels $(r,s)$ and $(p'-r,p-s)$ correspond to the same primary field. Finally, the conformal dimensions of the fields are given by
\begin{equation}
h(r,s) = \frac{(rp-sp')^2-(p'-p)^2}{4pp'} \ .
\end{equation}

\subsection{Gepner parafermions in terms of minimal models}

The non-Abelian part of the CFT describing the NASS state, was originally written in terms of Gepner parafermions, \cite{gepner1987} which can be expressed in terms of the coset,
\[
{\rm GPf}[su(3)_k]= \frac{su(3)_k}{u(1)_{2k} u(1)_{6k}}.
\]
For our present purposes, we shall focus on $k=2$, in which case the above coset is equivalently written as
\begin{equation}
{\rm GPf}[su(3)_2] = \frac{su(2)_1\times su(2)_1 \times su(2)_1}{su(2)_3}
\end{equation}
[for arbitrary $su(3)_k$ parafermions, $su(2)$ has to be replaced with $su(k)$ in the above]. If one multiplies the numerator and denominator of the coset above by $su(2)_2$, one can factorize the coset into the (semidirect) product of the minimal models
$\mathcal{M}(3,4)$ and $\mathcal{M}(4,5)$,
\begin{align}
{\rm GPf}[su(3)_2] &\approx \frac{su(2)_1\times su(2)_1 \times su(2)_1\times su(2)_2}
{su(2)_2 \times su(2)_3}  \nonumber \\ &\approx \mathcal{M}(3,4) \ltimes \mathcal{M}(4,5).
\end{align}

Let us be more precise about this correspondence. The {\em direct} product
$\mathcal{M}(3,4) \times \mathcal{M}(4,5)$ does not correspond to the $su(3)_2$ Gepner parafermions. Instead, one has to consider
a so-called different modular invariant, \cite{cappelli1987a,cappelli1987b} constructed from the fields present in the direct product.
To do that, we will follow the logic presented in Refs.~\onlinecite{bais2009,fuchs1996,frohlich2004}. This amounts to identifying a boson in the CFT, which is then said to be ``condensed'' (or added to the chiral algebra). To identify a suitable boson, we give the fields of the models
$\mathcal{M}(3,4)$ and $\mathcal{M}(4,5)$ in Table~\ref{Tab:KacM4534}. (Fusion rules for the minimal models are given in Refs. \onlinecite{byb,belavin1984}.)

\begin{table}[ht]
  \begin{tabular}{r|ccc}
  \hline \hline
  \multicolumn{4}{c}{$\mathcal{M}(3,4)$} \\
  \hline
  $h(r,s)$ & $s=1$ & $2$ & $3$\\
  \hline
  $r=1$ & 0 & 1/16 & 1/2 \\
  $2$     & 1/2 & 1/16 & 0 \\
  \hline \hline
  \end{tabular}
  \hskip 8mm
  \begin{tabular}{r|cccc}
  \hline \hline
  \multicolumn{5}{c}{$\mathcal{M}(4,5)$} \\
  \hline
  $h(r,s)$ & $s=1$ & $2$ & $3$ & $4$\\
  \hline
  $r=1$ & $0$ & $1/10$ & $3/5$ & $3/2$\\
  $2$ & $7/16$ & $3/80$ & $3/80$ & $7/16$\\
  $3$ & $3/2$ & $3/5$ & $1/10$ & $0$\\
   \hline \hline
  \end{tabular}
   \caption{Kac table of conformal weights for the unitary minimal models $\mathcal{M}(3,4)$ and
   $\mathcal{M}(4,5)$.}
  \label{Tab:KacM4534}
\end{table}

We will label the fields in the product theory in terms of the conformal dimensions (or conformal weights) of the contributing fields of the original CFTs. We find that there is indeed a Bosonic field (i.e., a field with integer scaling dimension) in the product, namely $(1/2,3/2)$. The condensation picture amounts to the following procedure: \cite{bais2009} Particles which can be obtained from one another by fusion of the boson, are ``identified'' (in the same way as the boson itself is identified with the vacuum, or trivial particle); in addition, particles which are not mutually local with the boson, are ``confined''; a particle which is mutual local with the boson does not generate a phase factor when transported around the boson; finally, it can happen that particles ``split''. We will see an example of splitting below. In this section, we will be rather brief; more details on how one constructs the correct theory can be found in the Appendix.

To get started, we list the particles in the product theory which are mutually local with this boson, namely
\begin{align*}
e &= (0,0) = (1/2,3/2),\\
(1/10)_a &= (1/16,3/80),\\
(1/10)_b &= (0,1/10) = (1/2,3/5),\\
(1/2)_a &= (1/16,7/16),\\
(1/2)_b &= (0,3/2) = (1/2,0),\\
3/5 &= (0,3/5) = (1/2,1/10).
\end{align*}
The equalities within one line signify that the particles corresponding to the two fields are identified. As an example, one obtains $(1/2,3/5)$ from
$(0,1/10)$ by fusion with the boson, namely $(0,1/10) \times (1/2,3/2) = (1/2,3/5)$.
We also note that the particles $(\frac{1}{2})_a$ and $(\frac{1}{10})_a$ have to split into two particles, thereby giving the full $Z_3$ symmetry present in the $su(3)_2$ parafermion theory. That this splitting is necessary can be seen by taking, for instance, the fusion product $(\frac{1}{2})_a\times (\frac{1}{2})_a = (0,0) + (1/2,0) + (0,3/2) + (1/2,3/2)$. Because both $(0,0)$ and $(1/2,3/2)$ correspond to the identity, splitting of the field $(\frac{1}{2})_a$ is necessary. A similar argument applies for $(\frac{1}{10})_a$, and one can convince oneself that in the end, one indeed obtains the correct fusion rules of the $su(3)_2$ parafermion theory. We conclude this section by mentioning that the example of the equivalence
${\rm GPf} [su(3)_2] \approx \mathcal{M}(3,4) \ltimes \mathcal{M}(4,5)$ featured prominently in a paper relating the Moore--Read and NASS states.\cite{grosfeld2009}

\subsection{CFT for the spin-singlet Gaffnian}

The spin-singlet Gaffnian that we aim to construct should have the same vanishing properties as the Gaffnian when considered as a function of only one spin species of particles. Motivated by the product of minimal models describing the $su(3)_2$ Gepner parafermions, we take as a starting point the coset
\begin{equation*}
\frac{su(2)_1 \times su(2)_1 \times su(2)_{k} }{su(2)_{k+2}} \ ,
\end{equation*}
with $k$ some fraction. Multiplying numerator and denominator by $su(2)_{k+1}$, we can write
\begin{equation*}
\frac{su(2)_1 \times su(2)_{k} }{su(2)_{k+1}}
\times
\frac{su(2)_1 \times su(2)_{k+1} }{su(2)_{k+2}}  .
\end{equation*}
The CFT describing the Gaffnian is the nonunitary minimal model $\mathcal{M}(3,5)$. To retain $\mathcal{M}(3,5)$ in the above construction,  we will choose $k=-1/2$. With this choice, our coset takes the form $\mathcal{M}(3,5)\times \mathcal{M}(5,7)$. In analogy with the Gepner parafermion case above, we look for a field with integer scaling dimension in the product theory, and condense it. The field content of the models is given in Table~\ref{Tab:KacM5735}.  (Fusion rules for the minimal models are given in Refs. \onlinecite{byb,belavin1984}.)

\begin{table}[ht]
  \begin{tabular}{r|cccc}
  \hline \hline
  \multicolumn{5}{c}{$\mathcal{M}(3,5)$} \\
  \hline
  $h(r,s)$ & $s=1$ & $2$ & $3$ & $4$\\
  \hline
  $r=1$ & 0 & $-$1/20 & 1/5 & 3/4 \\
  $2$     & 3/4 & 1/5 & $-$1/20 & 0 \\
  \hline \hline
  \end{tabular}
  \vskip 8mm
  \begin{tabular}{r|cccccc}
  \hline \hline
  \multicolumn{7}{c}{$\mathcal{M}(5,7)$} \\
  \hline
  $h(r,s)$ & $s=1$ & $2$ & $3$ & $4$ & $5$ & $6$\\
  \hline
  $r=1$ & 0 & 1/28 & 3/7 & 33/28 & 16/7 & 15/4 \\
  $2$ & 11/20 & 3/35 & $-$3/140 & 8/35 & 117/140 & 9/5 \\
  $3$ & 9/5 & 117/140 & 8/35 & $-$3/140 & 3/35 & 11/20 \\
  $4$ & 15/4 & 16/7 & 33/28 & 3/7 & 1/28 & 0 \\
   \hline \hline
  \end{tabular}
   \caption{Kac table of conformal weights for the unitary minimal models $\mathcal{M}(3,5)$ and
   $\mathcal{M}(5,7)$.}
  \label{Tab:KacM5735}
\end{table}

It turns out that the only boson that we can condense (or add to the chiral algebra), is the field $(9/5,1/5)$, which curiously has nontrivial fusion rules. The only fields that survive this condensation are found to correspond to product of fields of the first column of the Kac table for $\mathcal{M}(3,5)$ and the fields of the first row of the Kac table for $\mathcal{M}(5,7)$. In addition, none of these product fields are split. It follows that the resulting theory can be thought of as the product of a nonunitary semion theory (two fields with $Z_2$ fusion rules and dimensions $0$ and $3/4$) and a theory of six fields, satisfying $su(2)_5$ fusion rules. This second theory is also nonunitary, and the scaling dimensions are again simply read off from the Kac table, $(0,1/28,3/7,33/28,16/7,15/4)$. We note that this theory is modular (inherited via the coset construction). Here, we were rather brief in our description of the construction of the theory $\mathcal{M}(3,5) \ltimes \mathcal{M}(5,7)$. In Sec. 2 of the Appendix, we will give the details of the construction.

Before we start with the construction of the electron and quasihole operators, in Table~\ref{Tab:spingaffcft} we give the Kac table of the fields present in the CFT we constructed,
$\mathcal{M}(3,5) \ltimes \mathcal{M}(5,7)$, which will constitute the non-Abelian part of the CFT describing the \wf{}. The full CFT also contains the $u(1)$ vertex operators (see, for instance, Ref.~\onlinecite{moore1991}). It is important to note that one has to be careful in determining the scaling dimensions of the fields in the product theory. The dimensions obtained by simply adding the scaling dimensions of the constituent fields can in fact correspond to the scaling dimensions of descendant fields, differing from the scaling dimensions of the primaries by integers. In particular, the field $(15/4,0)$ is identified with $(9/5,1/5) \times (15/4,0) = (11/20,1/5)$, which has scaling dimension $3/4$. Similarly, $(9/5,1/5) \times (15/4,3/4) = (11/20,-1/20)$, which has scaling dimension $1/2$.

The fields occupying the corners of the Kac table, namely $(0,0), (0,3/4), (15/4,0), (15/4,3/4)$, are special, because they are simple currents. A simple current is a field, which when fused with any other field, always gives a single field as the result. Therefore, the particles corresponding to simple currents are Abelian. To construct the ``electron'' operator, one is therefore only allowed to use one of these four fields. Before we get started, we note that the Abelian sector of the theory at hand has less symmetry than the one describing the NASS state, which has three primary fields all with the same conformal weight.

\begin{table}
\begin{tabular}{r|cccccc}
\hline
\hline
\multicolumn{7}{c}{$\mathcal{M}(3,5) \ltimes \mathcal{M}(5,7)$} \\
\hline
$h(i,j)$ & $j=0$ & $1$ & $2$ & $3$ & $4$ & $5$\\
\hline
$i=0$ & 0 & 1/28 & 3/7 & 5/28 & 2/7 & 3/4 \\
$1$ & 3/4 & 11/14 & 5/28 & -1/14 & 1/28 & 1/2\\
\hline
\hline
\end{tabular}
\caption{Kac table of conformal weights
of the CFT describing the non-Abelian part of the spin Gaffnian.}
\label{Tab:spingaffcft}
\end{table}

We will now work with the assumption that we do have a $Z_2$ symmetry between the fields $(15/4,0)$ and $(0,3/4)$, [i.e., the fields with labels  $(i,j)=(0,5),(1,0)$] which we need, if the correlator is to describe a spin-singlet state (in the Appendix, we will give evidence supporting this statement). In addition, we want the wave function to look like a Gaffnian when viewed as a wave function for either spin-up or spin-down particles alone. We will introduce the following notation for the fields:
\begin{align}
\mathbf 1 &= (0,0), & \psi_\uparrow &= (0,3/4), \nonumber \\ \psi_\downarrow &= (15/4,0), & \psi_{\uparrow\downarrow} &= (15/4,3/4).
\end{align}
The fusion rules of these fields read $\psi_\uparrow \times \psi_\uparrow = \psi_\downarrow \times \psi_\downarrow = \psi_{\uparrow\downarrow} \times \psi_{\uparrow\downarrow} = \mathbf 1$ and $\psi_\uparrow \times \psi_\downarrow = \psi_{\uparrow\downarrow}$.

We will continue with adding the appropriate vertex operators, giving charge and spin to the constituent particles, which is the usual procedure (see, for instance, Ref. ~\onlinecite{moore1991}). The guiding principle will be to construct electron operators, which give rise to the vanishing properties we want. The ansatz for the operators is
\begin{align}
V_{\uparrow} (z_\uparrow) &= \psi_\uparrow (z_\uparrow) e^{i \alpha \phi_c + i \beta \phi_s} (z_\uparrow) \nonumber, \\
V_{\downarrow} (z_\downarrow) &= \psi_\downarrow (z_\downarrow) e^{i \alpha \phi_c - i \beta \phi_s} (z_\downarrow),
\end{align}
where $\alpha$ and $\beta$ are constants to be determined. The fields $\phi_c$ and $\phi_s$ are $u(1)$ compactified bosons. As we already pointed out, when we bring several up-particles together, we demand that the wave function behaves in the same way as the Gaffnian wave function. The electron operator in the CFT description of the Gaffnian reads $V (z) = \psi(z) e^{i \sqrt{\frac{3}{2}} \phi}(z)$, with $\psi(z)$ a field with conformal dimension $3/4$. Equivalences with this polarized case gives us the constraint $\alpha^2 + \beta^2 = 3/2$.

Moreover, bringing a spin-up and a spin-down particle together should not lead to a divergence. Using the fusion rule $V_\uparrow \times V_\downarrow = \psi_{\uparrow\downarrow} e^{2 i \alpha \phi_c}$ gives rise to a factor $(z_\uparrow-z_\downarrow)^{-1 + \alpha^2-\beta^2}$ in the operator product expansion. To avoid a pole, one obtains the constraint $\alpha^2-\beta^2 = p$, with $p$ an integer greater then or equal to 1. Picking the minimal choice $p=1$, one finds that $\alpha = \sqrt{5/4}$ and $\beta = 1/2$. This value for $\beta$ is in fact equal to the value this parameter takes in the case of the NASS state.

With the electron operators
$V_{\uparrow} = \psi_\uparrow e^{i \sqrt{5/4}\phi_c + 1/2 i \phi_s}$
and
$V_{\downarrow} = \psi_\downarrow e^{i \sqrt{5/4}\phi_c - 1/2 i \phi_s}$
in place, we conclude that the corresponding wave function does not vanish when any two particles come together. It vanishes as three powers when three up (or three down) particles come together. In addition, when two up and one down particle come together, the wave function generically vanishes quadratically. We then assume that one can use logic similar to the NASS case, and thus, exploiting the SU(2) symmetry, \cite{ardonne2001} one finds that the wave function vanishes as three powers when three particles in the $S=3/2$ channel are coincident, and quadratically when three particles in the $S=1/2$ channel are coincident.

The filling fraction of the wave function is given by $\nu = 4/5$. The scaling dimension $h_e = 3/4 + 3/4 = 3/2$ gives rise to a shift $\delta = 3$, the same as for the Gaffnian. We will defer the construction of the quasihole operators to the Appendix, but mention here that they are constrained by the fact that the wave functions for the electrons should not have poles, even in the presence of quasiholes. For $N$ particles the ground-state \wf{} wave function can be written as a conformal block containing equal numbers of spin-up and spin-down electron fields:
\[
\Psi_{\mbox{\wf{}}} =\left \langle \psi_{\uparrow} (z^{\uparrow}_1) \ldots  \psi_{\uparrow} (z^{\uparrow}_{\frac{N}{2}})  \psi_{\downarrow} (z^{\downarrow}_{\frac{N}{2}+1}) \ldots  \psi_{\uparrow} (z^{\downarrow}_{N}) \right \rangle.
\]

Finally, we would like to mention that if we would have assumed that the field $\psi_{\uparrow,\downarrow}$ would have had conformal dimension $3/2$ instead of $1/2$ (which is not consistent, because of the field identification above), one would have found $\alpha = \beta = \sqrt{3/4}$. The resulting wave function would then factorize as Gaffnian(up) $\times$ Gaffnian(down). The CFT description for such a wave function should be $\mathcal{M}(3,5)\times \mathcal{M}(3,5)$, which is at odds with the CFT description we have used here.

\section{Ground-state wave function}
\label{secwave function}

In this section we shall present an explicit construction for the ground-state \wf{} wave function for the Bosonic case. Our construction is conjectured with the view to satisfying the vanishing properties as well as the constraints of filling factor $\nu=4/5$ and shift $\delta=3$, arising from the CFT considerations put forward in Sec.~\ref{secCft}. In addition, the wave function must describe a spin-singlet state. That condition specifies the requirement to satisfy the Fock cyclic symmetry conditions (see Ref. \onlinecite{ham}).

Before we describe how to construct the \wf{}, it will be useful to motivate our methodology by briefly reviewing the construction of the spin-polarized Gaffnian wave function. \cite{simon2007b} For the Gaffnian, the CFT requirements imposed on the vanishing properties are that for any three particles coincident the wave function must vanish as three powers and that the wave function must not vanish for any two particles coincident [the CFT describing the Gaffnian is the nonunitary minimal model $\mathcal{M}(3,5)$]. Further, it was determined that in its Bosonic form, the Gaffnian occurs at filling $2/3$ and $\delta=3$.

To give the expression of the Gaffnian wave function that can be generalized to the \wf{}, we divide the particles into two groups $A$ and $B$ of equal size. The Gaffnian wave function can then be written as
\begin{widetext}
\[
\Psi_{\mbox{Gaffnian}} = \hat S \left [
\prod_{i<j \in A}(z^{A}_i -z^{A}_j)^2
\prod_{i<j \in B}(z^{B}_i -z^{B}_j)^2
\prod_{i \in A, j \in B} (z^{A}_i -z^{B}_j) \;
\mathbf{Per} \left [ \frac{1}{z^{A}_i-z^{B}_j} \right ]
\right ],
\]
\end{widetext}
where $\hat S$ represents a symmetrization operation over all $N$ particle coordinates, and
$\mathbf{Per} [ M_{ij}] $ denotes the ``permanent'' of a matrix $M$ whose elements are in this case given by $M_{ij} = (z^{A}_i-z^{B}_j)^{-1}$. On its own, this permanent factor contains a certain pattern of poles.
When placed within the full construction, these poles conspire to ensure that the overall wave function does not vanish as two particles become coincident. Thus the polynomial vanishes only when three particles are coincident. The vanishing power can be tuned by adjusting the exponents of each Jastrow-type factor in the construction. For example,  $\prod_{i<j \in A}(z^{A}_i -z^{A}_j)^2$ could be adjusted to, say, $\prod_{i<j \in A}(z^{A}_i -z^{A}_j)^4$ (and similar for the particles in group $B$) to alter the vanishing power from 3 to 5 in this example. An entirely nontrivial step is to determine whether or not these vanishing properties are retained once the overall symmetrization operation has been completed. For the Gaffnian, it was found that the vanishing properties are retained.

Now let us turn to the construction of the ground-state wave function for the \wf{}. We have determined that the following construction, given in Eq.~(\ref{eqwave function}), gives a wave function at filling factor $\nu=4/5$ and $\delta=3$ and it satisfies the CFT vanishing constraints for $N=4,8,12$. We have reason to believe that it will work for all other values of $N$ (i.e., 16, 20, etc.), although an explicit check is not possible. Thus, we conclude that it must be proportional to the conformal block giving the \wf{} ground state. In the construction, we have used a clustering principle: The particle coordinates are first divided into equal sets of spin up and spin down, and then each of these sets is further subdivided into two equal subsets labelled by either $A$ or $B$, giving four sets in total ($A\uparrow$, $A\downarrow$, $B\uparrow$, $B\downarrow$). With the label $A$ or $B$ on its own, we refer to the two sub-sets with either spin direction,
\begin{widetext}
\begin{equation}
\label{eqwave function}
\Psi_{\mbox{\wf{}}} = \hat Y_{S=0}  \left [   \prod_{a=A,B} \left\{ \prod_{i<j \in a\uparrow}(z^{a\uparrow}_i -z^{a\uparrow}_j)^2   \prod_{i<j \in a\downarrow}(z^{a\downarrow}_i -z^{a\downarrow}_j)^2  \prod_{i \in a\uparrow, j \in a\downarrow}(z^{a\uparrow}_i -z^{a\downarrow}_j)  \right \}   \prod_{i \in A, j \in B} (z^{A}_i -z^{B}_j) \; \mathbf{Per} \left [ \frac{1}{z^{A}_i-z^{B}_j} \right ] \right ],
\end{equation}
\end{widetext}
where, as above, $\mathbf{Per} [ M_{ij}] $ denotes the permanent of a matrix $M$ whose elements are in this case given by $M_{ij} = (z^{A}_i-z^{B}_j)^{-1} $ and where $\hat Y_{S=0}$ is the Young operator for a spin-singlet representation of the symmetric group---this operation is required in order to guarantee that the wave function satisfies the correct Fock cyclic symmetry conditions for a spin-singlet state. \cite{ham}
\[
\hat Y_{S=0} = \hat S_{z^{\uparrow}_1 \ldots z^{\uparrow}_{N/2} } \hat S_{z^{\downarrow}_1 \ldots z^{\downarrow}_{N/2} } \hat A_{z^{\uparrow}_1 z^{\downarrow}_1}\ldots  \hat A_{z^{\uparrow}_{N/2} z^{\downarrow}_{N/2}},
\]
that is, the Young operator for a spin-singlet representation corresponds to the operation of antisymmetrizing over ordered pairs of spin-up and spin-down coordinates, followed by symmetrizating over all spin-down and then all spin-up coordinates. Crucially, it is important to note that the wave function does not completely vanish when $\hat Y_{S=0}$ is applied!

In Sec.~\ref{secNumerics} we provided strong evidence showing that the spectrum of the Hamiltonian for the \wf{} contains a unique zero energy state with $L=0$ and $S=0$, for precisely the values of particle number and flux corresponding to filling $\nu=4/5$ and $\delta=3$. Therefore we conclude that this trial ground-state wave function is unique, and further, it is precisely the highest density ground state of the \wf{} Hamiltonian proposed in Eq.~(\ref{eqHamiltonian}).

\section{Discussion}
\label{discussion}

To summarize our findings, we have presented evidence to demonstrate that the proposed \wf{} wave function satisfies many of the ``ingredients of an ideal theory of a FQHE state'': A local Hamiltonian, a relatively simple analytic expression for the wave function, at least for the ground state, and a correspondence of that wave function to a 2D rational CFT. We have verified the quasihole spectrum of the \wf{} Hamiltonian, checking against both the spin dressed squeezing algorithm and the spinful version of the generalized Pauli principle. Indeed, this study is an interesting test case for the application of such techniques.

Concerning quasiholes, although we have made progress on the CFT description of the quasihole operators (see the Appendix), we have not yet determined simple, analytic expressions for the corresponding quasihole wave functions. For the spin-polarized Gaffnian, the quasihole wave functions have been constructed; \cite{simon2007b} the main stumbling block here is the additional complexity due to the spin degree of freedom.

In Sec.~\ref{secMethods}  we described how to generate the \wf{} wave functions by means of a local Hamiltonian written in terms of spin-dependent pseudopotentials. An interesting question left unaddressed is whether other states could be constructed with faster vanishing properties than the \wf{}. For instance, in the spin-polarized case, adding the next highest $L$ pseudopotential to the Gaffnian Hamiltonian is known to produce the Haffnian Hamiltonian. \cite{green2002} Might we be able to generate a ``spin-Haffnian''  state with a Hamiltonian containing positive
$ V^3_{0,3/2}$, $ V^3_{2,3/2,}$, $ V^3_{1,1/2}$ and now, in addition, $V^3_{3,3/2}$ (and possibly $ V^3_{2,1/2}$)? Wave functions of this type were also considered in Ref.~\onlinecite{estienne2011}.
Presently we lack a corresponding CFT description with which to conduct the same checks as for the \wf{}. It is also unclear at the present time whether that Hilbert space could identically be constructed by a squeezing algorithm or generalized Pauli principle approach.

Concerning CFT coset constructions, it is noteworthy that we were able to derive a self-consistent CFT from the product of two nonunitary minimal models by constructing a nondiagonal modular invariant. Specifically, the nonunitary models in the product can be thought of as cosets of fractional level affine Lie algebras. The resulting CFT corresponds to the nonunitary coset $su(2)_1\times su(2)_1 \times su(2)_{k-2}/su(2)_k$, with $k$ fractional. In Sec.~2 of the Appendix, we deal with a more general case.

In general, for \emph{any} unitary coset construction with simple current extensions there is a proven procedure to generate the CFT. \cite{fuchs1996,frohlich2004} For nonunitary cosets, no such general procedure exists, and each CFT must be constructed on a case-by-case basis. \cite{byb} In the CFT construction presented in this work, we have provided a further example of a case where the construction of non-diagonal modular invariants for nonunitary CFTs is possible.

Another powerful tool for the analysis of the quantum Hall wave functions has been the concept of Jack polynomials. \cite{bernevig2008b,estienne2011} The Jack polynomials provide a convenient basis in which to describe polynomial wave functions with precisely defined vanishing properties.  Based on this new insight, it has been determined that states such as that of Moore and Read---where there exists a corresponding local Hamiltonian, a relatively simple analytic form for the wave function, and a CFT description of the state---are in fact simply special cases within a much broader classification of FQHE states in terms of these Jack polynomials. It has been further realized that such states are rather atypical: Out of all the Jacks, there exist only a handful of special cases for which all three of these ideal ingredients of a theory of a FQHE state are present. In this paper we have presented a further example, albeit that the CFT corresponding to the SSG is nonunitary implying a gapless, compressible wave function.

\acknowledgements

We thank J.~K.~Slingerland, B.~Schellekens, B.~Estienne and B.~A.~Bernevig for helpful discussions. We also thank B.~Estienne for useful comments on the manuscript. S.C.D. and S.H.S. were supported by EPSRC Grants No. EP/I032487/1 and No. EP/I031014/1. N.R. was supported by  NSF CAREER DMR-095242, ONR-N00014-11-1-0635,  ARMY-245-6778, MURI-130-6082, ANR-12-BS04-0002-02, Packard Foundation, and a Keck grant. We thank NORDITA, the Institute Henri Poincar\'e and the Aspen Centre for Physics for hospitality during our collaboration.