Convergence analysis of a stochastic interacting particle-field algorithm for 3D parabolic-parabolic Keller-Segel systems

Chemotaxis models describe the movement of organisms in response to chemical gradients. In this paper, we present a stochastic interacting particle-field algorithm with a random batch approximation (SIPF-$r$) for the three-dimensional (3D) parabolic-parabolic Keller-Segel (KS) system, also referred to as the fully parabolic KS system. The SIPF-$r$ method approximates the KS system by coupling particle-based representations of the density with a smooth field variable computed using spectral methods. By incorporating the random batch method (RBM), we bypass the mean-field limit and significantly reduce computational complexity. Under mild assumptions on the regularity of the original KS system and the boundedness of numerical approximations, we prove that the empirical measure of the SIPF-$r$ particle system converges, with high probability, to the exact measure of the limiting McKean-Vlasov process in the $1$-Wasserstein distance. Finally, we present numerical experiments to validate the theoretical convergence rates, and demonstrate the performance and robustness of the SIPF-$r$ method as a diagnostic tool for intense focusing and potential finite-time singularity in 3D, subject to critical initial mass thresholds in the system.