Variational wavefunction for multi-species spinful fermionic superfluids and superconductors

We present numerical solutions for two electron species in 2 dimensions. Based on these solutions, we show that, for equivalent spin up and down fermions, the Fermi occupation factor (per spin) equals 1/2 deep in the Fermi sea. This constitutes a unique experimental prediction of the theory, both for the normal and superfluid states. Interestingly, this result, obtained in the thermodynamic limit, is consistent with Fermi occupation factor (in-)equalities for finite systems of electrons, derived (in a different context) by Borland and Dennis (1972) and by Altunbulak and Klyachko (2008).