Star Formation Mitosis: The Effects of Particle Splitting on the Star Formation Outcome

Particle splitting can increase the dynamic range in Lagrangian astrophysical simulations by several orders of magnitude. However, such methods can introduce numerical noise that may affect the rate of convergence, particularly for collapse problems where fragmentation could be artificially suppressed or enhanced. Thus it is necessary to determine how these errors propagate to key quantities of interest, and to correct or mitigate them where possible. Using GIZMO's MFM-MHD method, we perform a set of about 220 simulations of a uniform-density turbulent gas sphere undergoing collapse, systematically varying the cloud mass and radius, feedback physics, initial resolution, the pace and cadence at which we split particles and the sink (accreting) or single-age stellar population (SSP; non-accreting) particle formation prescription. Crucially, we compare with control runs that start with uniformly-high resolution without splitting. While we find no evidence for artificial fragmentation caused by splitting particles, particle splitting runs have a smoother density field than their non-splitting counterparts. Due to this, particle splitting runs exhibit minor differences in the density and turbulence statistics. However, they form about the same total stellar mass as the control non-splitting runs, albeit with a slightly different mass spectrum. As a result, particle splitting runs have a different evolutionary history, compared to runs without particle splitting, especially when stellar feedback is present. Despite the changes seen in runs where we split particles, the numerical solution lies within the scatter inherently present in these chaotic systems, which we model by modifying the initial turbulent velocity field. Our results here argue that particle splitting generally ``works" for star formation applications, despite some systematic differences in the specifics.