Search

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We study the time evolution of molecular clouds across three Milky Way-like isolated disc galaxy simulations at a temporal resolution of 1 Myr and at a range of spatial resolutions spanning two orders of magnitude in spatial scale from ∼10 pc up to ∼1 kpc. The cloud evolution networks generated at the highest spatial resolution contain a cumulative total of ∼80 000 separate molecular clouds in different galactic-dynamical environments. We find that clouds undergo mergers at a rate proportional to the crossing time between their centroids, but that their physical properties are largely insensitive to these interactions. Below the gas-disc scale height, the cloud lifetime τ life obeys a scaling relation of the form τ life ∝ℓ -0.3 with the cloud size ℓ, consistent with over-densities that collapse, form stars, and are dispersed by stellar feedback. Above the disc scale height, these self-gravitating regions are no longer resolved, so the scaling relation flattens to a constant value of ∼13 Myr, consistent with the turbulent crossing time of the gas disc, as observed in nearby disc galaxies., (© 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society.)