Integrable model of a p-wave bosonic superfluid

We present an exactly solvable p-wave pairing model for two bosonic species. The model is solvable in any spatial dimension and shares some commonalities with the p + ip Richardson-Gaudin fermionic model, such as a third-order quantum phase transition. However, contrary to the fermionic case, in the bosonic model the transition separates a gapless fragmented singlet pair condensate from a pair Bose superfluid, and the exact eigenstate at the quantum critical point is a pair condensate analogous to the fermionic Moore-Read state.