Multiparticle collision framework for active polar fluids

Sufficiently dense intrinsically out-of-equilibrium suspensions, such as those observed in biological systems, can be modeled as active fluids characterized by their orientational symmetry. While mesoscale numerical approaches to active nematic fluids have been developed, polar fluids are simulated as either ensembles of microscopic self-propelled particles or continuous hydrodynamic-scale equations of motion. To better simulate active polar fluids in complex geometries or as a solvent for suspensions, mesoscale numerical approaches are needed. In this work, the coarse-graining multiparticle collision dynamics (MPCD) framework is applied to three active particle models to produce mesoscale simulations of polar active fluids. The first active-polar MPCD (AP-MPCD) is a variant of the Vicsek model, while the second and third variants allow the speed of the particles to relax towards a self-propulsion speed subject to Andersen and Langevin thermostats, respectively. Each of these AP-MPCD variants exhibits a flocking transition at a critical activity and banding in the vicinity of the transition point. We leverage the mesoscale nature of AP-MPCD to explore flocking in the presence of external fields, which destroys banding, and anisotropic obstacles, which act as a ratchet that biases the flocking direction. These results demonstrate the capacity of AP-MPCD to capture the known phenomenology of polar active suspensions and its versatility to study active polar fluids in complex scenarios.