Mitigating density fluctuations in particle-based active nematic simulations

Understanding active matter has led to new perspectives on biophysics and non-equilibrium dynamics. However, the development of numerical tools for simulating active fluids capable of incorporating non-trivial boundaries or inclusions has lagged behind. Active particle-based methods, which typically excel at this, suffer from large density fluctuations that affect the dynamics of inclusions. To this end, we advance the Active-Nematic Multi-Particle Collision Dynamics algorithm, a particle-based method for simulating active nematics, by addressing the large density fluctuations that arise from activity. This paper introduces three activity formulations that mitigate the coupling between activity and local density. Local density fluctuations are decreased to a level comparable to the passive limit while retaining active nematic phenomenology and increasing the active turbulence regime four-fold in two dimensions. These developments extend the technique into a flexible tool for modeling active systems, including solutes and inclusions, with broad applications for the study of biophysical systems. Particle-based models of soft active matter are known to exhibit large density fluctuations, which may be problematic for considering homogeneous systems. Through considering the coupling of activity to local density, we reduce density fluctuations while retaining active turbulent dynamics.