Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCosmological SPH simulations: A hybrid multi-phase model for star formation
60
0
0.0
(
0
)
Authors
Volker Springel
Ask ChatGPT about the research
No Arabic abstract
We present a model for star formation and supernova feedback that describes the multi-phase structure of star forming gas on scales that are typically not resolved in cosmological simulations. Our approach includes radiative heating and cooling, the growth of cold clouds embedded in an ambient hot medium, star formation in these clouds, feedback from supernovae in the form of thermal heating and cloud evaporation, galactic winds and outflows, and metal enrichment. Implemented using SPH, our scheme is a significantly modified and extended version of the grid-based method of Yepes et al. (1997), and enables us to achieve high dynamic range in simulations of structure formation. We discuss properties of the feedback model in detail and show that it predicts a self-regulated, quiescent mode of star formation, which, in particular, stabilises the star forming gaseous layers of disk galaxies. The parameterisation of this mode can be reduced to a single free quantity which determines the overall timescale for star formation. We fix this parameter to match the observed rates of star formation in local disk galaxies. When normalised in this manner, cosmological simulations nevertheless overproduce the observed cosmic abundance of stellar material. We are thus motivated to extend our feedback model to include galactic winds associated with star formation. Using small-scale simulations of individual star-forming disk galaxies, we show that these winds produce either galactic fountains or outflows, depending on the depth of the gravitational potential. Moreover, outflows from galaxies in these simulations drive chemical enrichment of the intergalactic medium, in principle accounting for the presence of metals in the Lyman alpha forest. (abridged)
rate research
Read More
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$times10^{11}$ M$_odot$ to 2$times10^{12}$ M$_odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.
We present a multi-epoch analysis of the galaxy populations formed within the cosmological hydrodynamical simulations presented in Vogelsberger et al. (2013). These simulations explore the performance of a recently implemented feedback model which includes primordial and metal line radiative cooling with self-shielding corrections; stellar evolution with associated mass loss and chemical enrichment; feedback by stellar winds; black hole seeding, growth and merging; and AGN quasar- and radio-mode heating with a phenomenological prescription for AGN electro-magnetic feedback. We illustrate the impact of the model parameter choices on the resulting simulated galaxy population properties at high and intermediate redshifts. We demonstrate that our scheme is capable of producing galaxy populations that broadly reproduce the observed galaxy stellar mass function extending from redshift z=0 to z=3. We also characterise the evolving galactic B-band luminosity function, stellar mass to halo mass ratio, star formation main sequence, Tully-Fisher relation, and gas-phase mass-metallicity relation and confront them against recent observational estimates. This detailed comparison allows us to validate elements of our feedback model, while also identifying areas of tension that will be addressed in future work.
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.
We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $left(f_*right)$, the efficiency of thermal energy from stellar feedback $left(epsilonright)$ and the region into which feedback is injected $left(r {rm and} sright)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10},mathrm{M_odot}$ and $10^{12},mathrm{M_odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5,{rm M_odot}$ and $10^7,{rm M_odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10% and 25% in the stellar and baryon mass in simulations evolved from identical initial conditions.
We discuss differences in simulation results that arise between the use of either the thermal energy or the entropy as an independent variable in smoothed particle hydrodynamics (SPH). In this context, we derive a new version of SPH that manifestly conserves both energy and entropy if smoothing lengths are allowed to adapt freely to the local mass resolution. To test various formu- lations of SPH, we consider point-like energy injection and find that powerful explosions are well represented by SPH even when the energy is deposited into a single particle, provided that the entropy equation is integrated. If the thermal energy is instead used as an independent variable, unphysical solutions can be obtained for this problem. We also examine the radiative cooling of gas spheres that collapse and virialize in isolation and of halos that form in cosmological simulations of structure formation. When applied to these problems, the thermal energy version of SPH leads to substantial overcooling in halos that are resolved with up to a few thousand particles, while the entropy formulation is biased only moderately low for these halos. For objects resolved with much larger particle numbers, the two approaches yield consistent results. We trace the origin of the differences to systematic resolution effects in the outer parts of cooling flows. The cumulative effect of this overcooling can be significant. In cosmological simulations of moderate size, we find that the fraction of baryons which cool and condense can be reduced by up to a factor ~2 if the entropy equation is employed rather than the thermal energy equation. We also demonstrate that the entropy method leads to a greatly reduced scatter in the density-temperature relation of the low-density Ly-alpha forest relative to the thermal energy approach, in accord with theoretical expectations.(abridged)
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
