Snake instability of dark solitons in fermionic superfluids

We present numerical calculations of the snake instability in a Fermi superfluid within the Bogoliubov-de Gennes theory of the Bose-Einstein condensate (BEC) to BCS crossover using the random-phase approximation complemented by time-dependent simulations. We examine the snaking behavior across the crossover and quantify the time scale and length scale of the instability. While the dynamics shows extensive snaking before eventually producing vortices and sound on the BEC side of the crossover, the snaking dynamics is preempted by decay into sound due to pair breaking in the deep BCS regime. At the unitarity limit, hydrodynamic arguments allow us to link the rate of snaking to the experimentally observable ratio of inertial to physical mass of the soliton. In this limit we witness an unresolved discrepancy between our numerical estimates for the critical wave number of suppression of the snake instability and recent experimental observations with an ultracold Fermi gas.