Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH

Federrath, C., Banerjee, R., Clark, P. C., Klessen, R. S., 2010

The Astrophysical Journal, 713, 269  [ ADS link ]  [ PDF ]

Abstract

Star formation is such a complex process that accurate numerical tools are needed to quantitatively examine the mass distribution and accretion of fragments in collapsing, turbulent, magnetized gas clouds. To enable a numerical treatment of this regime, we implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu, and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a density threshold for sink particle creation alone is insufficient for supersonic flows, if the density threshold is below the opacity limit. The density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.

FLASH versus SPH comparison

The movie shows the formation of a star cluster from turbulent initial conditions. We find encouraging agreement between the adaptive mesh refinement code FLASH and a standard smoothed particle hydrodynamics (SPH) code with sink particles.

 

FLASH_vs_SPH_FederrathEtAl2010.mp4, 2.3MB high-res mp4 ]

Sink particle creation checks

Demonstrated in this movie is that physical checks for a bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. If those checks are switched off (right-hand side), many spurious sink particles are created, and the mass accreted by the sink particles is overestimated.

 

NoSinkChecks_FederrathEtAl2010.mp4, 2.3MB high-res mp4 ]

Angular momentum conservation test

This movie shows the formation and fragmentation of a thin gas disk from a rotating one solar mass cloud core. This test is also known as the Boss & Bodenheimer and Bate & Burkert tests. It shows the column density of the disk (face on and edge on) until 68.5% of the initial gas mass was accreted by sink particles and a total of 12 fragments have formed by that time. Below is also a plot of the time evolution of the total angular momentum along the rotation axis. The angular momentum is conserved to within 0.5% for all times.

 

BB_Federrath.mp4, 2.8MB high-res mp4 ]

BB_angmom_z.png

Momentum conservation test

This movie shows the collapse of a (very poorly resolved) gas cloud (0.5 solar masses) on the left, while a sink particle (0.1 solar masses) with initial momentum in positive vertical direction was placed on the right as an initial condition. The gas cloud collapses and a sink particle is formed in its center, while the other particle gets attracted by the gas cloud and goes into highly eccentric orbits. Both sink particles accrete almost all the gas very early on in the dynamical evolution. The two images below show both the momentum of the gas and sink particles in the horizontal and vertical directions, respectively. For momentum conservation to hold both gas and sink particle momenta must be symmetric about the horizontal axis. The discontinuities in the curves are caused by gas accretion events, during which momentum is also conserved.

 

sinks_osci_dens.mp4, 2.2MB high-res mp4 ]

sinks_osci_dens_xmom.png sinks_osci_dens_ymom.png

Acknowledgements

C.F. thanks Daniel Price for useful discussions on the gravitational softening and time integration schemes for sink particles during his visit at Monash University under the GO8/DAAD Joint Research Cooperation Scheme. We thank the anonymous referee for an unbiased and detailed report, which helped to improve the paper. C.F. is grateful for financial support from the International Max Planck Research School for Astronomy and Cosmic Physics (IMPRS-A) and the Heidelberg Graduate School of Fundamental Physics (HGSFP), which is funded by the Excellence Initiative of the Deutsche Forschungsgemeinschaft (DFG) GSC 129/1. R.B. is grateful to the DFG, which funds his Emmy Noether Grant BA 3607/1-1. R.B. is also thankful for subsidies from the FRONTIER initiative of the University of Heidelberg. P.C.C. acknowledges support from the DFG under the Emmy Noether Grant KL 1358/1 and the Priority Program SFB 439 "Galaxies in the Early Universe". R.S.K. thanks for financial resources provided by the DFG under grants No. KL 1358/1, KL 1358/4, and KL 1359/5. R.S.K. furthermore thanks for subsidies from the FRONTIER grant of Heidelberg University sponsored by the German Excellence Initiative and for support from the Landesstiftung Baden-Württemberg via their program International Collaboration II under grant P-LS-SPII/18. All authors acknowledge financial support from the German Bundesministerium für Bildung und Forschung via the ASTRONET project STAR FORMAT (grant 05A09VHA). The FLASH and SPH simulations used computational resources from the HLRB II project grant pr32lo at the Leibniz Rechenzentrum Garching and from a CASPUR project grant. The test simulations were partly run on the GPU cluster KOLOB at the University of Heidelberg. Figure 5 was produced with the visualization tool VISIT. The FLASH code was developed in part by the DOE-supported Alliances Center for Astrophysical Thermonuclear Flashes (ASC) at the University of Chicago.


© C. Federrath 2025