  Recent experiments (G. Ariel, et al., Nature Comm. 6, 8396 (2015)) revealed
an intriguing behavior of swarming bacteria: they fundamentally change their
collective motion from simple diffusion into a superdiffusive L\'{e}vy walk
dynamics. We introduce a nonlinear non-Markovian persistent random walk model
that explains the emergence of superdiffusive L\'{e}vy walks. We show that the
alignment interaction between individuals can lead to the superdiffusive growth
of the mean squared displacement and the power law distribution of run length
with infinite variance. The main result is that the superdiffusive behavior
emerges as a nonlinear collective phenomenon, rather than due to the standard
assumption of the power law distribution of run distances from the inception.
At the same time, we find that the repulsion/collision effects lead to the
density dependent exponential tempering of power law distributions. This
qualitatively explains experimentally observed transition from superdiffusion
to the diffusion of mussels as their density increases (M. de Jager et al.,
Proc. R. Soc. B 281, 20132605 (2014)).
