Closures of moment expansion of anisotropic active Brownian particles

We study analytically the dynamics of anisotropic active Brownian particles (ABPs), and more precisely their intermediate scattering function (ISF). To this end, we develop a systematic closure scheme for the moment expansion of their Fokker-Planck equation. Starting from the coupled evolution of translational and orientational degrees of freedom, we derive equations for the density, polarization, and nematic tensor fields, which naturally generate an infinite hierarchy of higher-order moments. To obtain explicit solutions, we investigate truncation strategies and analyze closures at different orders. While the closure at lowest order yields Gaussian dynamics with an effective translational diffusion, closures at higher orders incorporate orientational correlations and reproduce non-Gaussian features in the ISF. By confronting these approximations with exact solutions based on spheroidal wave functions and with Brownian dynamics simulations, we identify their range of validity in terms of Péclet number, wavenumber, and observation timescales. An advantage of this method is its ability to yield approximate yet explicit expressions not only for the ISF but also for polarization and nematic fields, which are often neglected but relevant in scattering experiments and theoretical modeling. Beyond providing a practical guide to select the appropriate closure according to the spatiotemporal regime, our framework highlights the efficiency of moment-based approaches compared to exact yet implicit formulations. This strategy can be systematically extended to more complex situations, such as propulsion switching, confinement, or external fields, where functional bases for exact solutions are generally unavailable.