We propose a model of superdiffusive Lévy walk as an emergent nonlinear phenomenon in systems of interacting individuals. The aim is to provide a qualitative explanation of recent experiments [G. Ariel et al., Nat. Commun. 6, 8396 (2015)] revealing an intriguing behavior: swarming bacteria fundamentally change their collective motion from simple diffusion into a superdiffusive Lévy walk dynamics. We introduce microscopic mean-field kinetic equations in which we combine two key ingredients: (1) alignment interactions between individuals and (2) non-Markovian effects. Our interacting run-and-tumble model leads 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 cooperative effect without using the standard assumption of the power-law distribution of run distances from the inception. At the same time, we find that the collision and repulsion interactions lead to the density-dependent exponential tempering of power-law distributions. This qualitatively explains the experimentally observed transition from superdiffusion to the diffusion of mussels as their density increases [M. [15][16][17][18]. Recently an intriguing behavior of swarming bacteria was found: they fundamentally change their collective motion from simple diffusion into a superdiffusive Lévy walk dynamics [19]. The extraordinary feature of this movement is that the emergence of a superdiffusive motility is a result of the interactions between bacteria rather than the standard mechanism of controlling the individual frequency of tumbling. However, it is still an open question how Lévy walks emerge in systems of interacting self-propelled particles. The current theory of Lévy walk [8] assumes noninteracting particles and power-law distribution of traveled distances from the inception.The collective behavior of large groups of interacting individuals such as bird flocks, fish schools, and the collective migration of cells or bacteria is another rapidly growing area of active matter research [20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37][38]. There exist two main types of models used for a collective behavior: (1) Lagrangian models describing the movements of self-propelled particles in terms of nonlinear equations for the positions and velocities of all particles [20][21][22][23], and (2) kinetic models involving partial differential equations for the population densities [24][25][26][27][28][29][30]. In this Rapid Communication we are only concerned with the nonlinear kinetic and macroscopic equations. They have been used to describe interactions between individuals and investigate the formation of a large variety of spatiotemporal self-organized aggregations. Most of these models converge to the density-dependent diffusive transport and do not lead to non-Markovian Lévy walks.In this Rapid Communication we propose a kinetic nonlinear Lévy walk model for interacting individuals in which...