Just as a supplement to ACuriousMind's answer, it is worth noting that buried in the bottom of their paper they actually show what the "spin 1/2" eigenstates are in terms of the regular basis:
$|j=1/2\rangle=\frac{1}{\sqrt{2}}(|1, -1 \rangle + |0,1\rangle$)
$|j=-1/2\rangle=\frac{1}{\sqrt{2}}(|-1, 1 \rangle + |0,-1\rangle$)
where $|l, \sigma\rangle$ is the angular momentum in the normal $|l,s\rangle$ basis. Written out explicitly, it is clear:
- That these are eigenstates of $L+S/2$, and
- That this trick could only work for an integer or half-integer $\gamma$, and
- That there's nothing too special going on here.
Still, it can still be interesting to frame an old system in a new way. I knew about the possibility of anyons in low dimensions, but I still would not have guessed that the very natural and common reduction in symmetry caused by picking a propagation axis might be sufficient for this effect. However, that might is an important qualifier: since the authors don't actually demonstrate fractional statistics or a procedure to measure them, this remains to be seen.
Edit: Emilio asks for a concrete demonstration of point two:
We want
$(L+\gamma S)(\alpha|l_1,1\rangle+\beta|l_2,-1\rangle)=j(\alpha|l_1,1\rangle+\beta|l_2,-1\rangle)$ .
This is the most general angular momentum superposition possible with a fixed total $j$, since there are only two spin possibilities. Furthermore, we know from beginning QM classes that an eigenstate of $j$ will have all these possible elements with some Clebsch-Gordon coefficients.
Applying the operators:
$((l_1+\gamma)\alpha|l_1,1\rangle+(l_2-\gamma)\beta|l_2,-1\rangle)=j(\alpha|l_1,1\rangle+\beta|l_2,-1\rangle)$
giving the conditions
$l_1+\gamma=j$
$l_2-\gamma=j$,
or $l_2-l_1=2\gamma$, $l_1\neq l_2$ .
Since $l$ are integers, this implies that $\gamma$ must be a half-integer.