Theoretical Challenges in Galaxy Formation
TL;DR: In this paper, the authors focus on one major challenge for galaxy formation theory: to understand the underlying physical processes that regulate the structure of the interstellar medium, star formation and the driving of galactic outflows.
Abstract: Numerical simulations have become a major tool for understanding galaxy formation and evolution. Over the decades the field has made significant progress. It is now possible to simulate the formation of individual galaxies and galaxy populations from well defined initial conditions with realistic abundances and global properties. An essential component of the calculation is to correctly estimate the inflow to and outflow from forming galaxies since observations indicating low formation efficiency and strong circum-glactic presence of gas are persuasive. Energetic 'feedback' from massive stars and accreting super-massive black holes - generally unresolved in cosmological simulations - plays a major role for driving galactic outflows, which have been shown to regulate many aspects of galaxy evolution. A surprisingly large variety of plausible sub-resolution models succeeds in this exercise. They capture the essential characteristics of the problem, i.e. outflows regulating galactic gas flows, but their predictive power is limited. In this review we focus on one major challenge for galaxy formation theory: to understand the underlying physical processes that regulate the structure of the interstellar medium, star formation and the driving of galactic outflows. This requires accurate physical models and numerical simulations, which can precisely describe the multi-phase structure of the interstellar medium on the currently unresolved few hundred parsecs scales of large scale cosmological simulations. Such models ultimately require the full accounting for the dominant cooling and heating processes, the radiation and winds from massive stars and accreting black holes, an accurate treatment of supernova explosions as well as the non-thermal components of the interstellar medium like magnetic fields and cosmic rays.
Citations
More filters
TL;DR: In this article, an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO is introduced.
Abstract: We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multi-mode feedback. In this paper we give a comprehensive description of the physical and numerical advances which form the core of the IllustrisTNG (The Next Generation) framework. We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content, and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black hole driven kinetic feedback at low accretion rates, magnetohydrodynamics, and improvements to the numerical scheme. Using a suite of (25 Mpc $h^{-1}$)$^3$ cosmological boxes we assess the outcome of the new model at our fiducial resolution. The presence of a self-consistently amplified magnetic field is shown to have an important impact on the stellar content of $10^{12} M_{\rm sun}$ haloes and above. Finally, we demonstrate that the new galactic winds promise to solve key problems identified in Illustris in matching observational constraints and affecting the stellar content and sizes of the low mass end of the galaxy population.
1,105 citations
TL;DR: Simba as discussed by the authors is the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics, which includes updates to Mufaa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angles-Alcazar et al.
Abstract: We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angles-Alcazar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies.
518 citations
Cites background from "Theoretical Challenges in Galaxy Fo..."
...While this framework has been broadly successful at reproducing many key observed characteristics of the galaxy population at a range of cosmic epochs, the physical understanding of most of these small-scale feedback processes remains coarse and heuristic (see Naab & Ostriker 2017, for a review)....
[...]
TL;DR: In this paper, the authors provide an overview of the methods used to constrain the chemical enrichment in galaxies and their environment, and discuss the observed scaling relations between metallicity and galaxy properties, the observed relative chemical abundances, how the chemical elements are distributed within galaxies, and how these properties evolve across the cosmic epochs.
Abstract: The evolution of the content of heavy elements in galaxies, the relative chemical abundances, their spatial distribution, and how these scale with various galactic properties, provide unique information on the galactic evolutionary processes across the cosmic epochs. In recent years major progress has been made in constraining the chemical evolution of galaxies and inferring key information relevant to our understanding of the main mechanisms involved in galaxy evolution. In this review we provide an overview of these various areas. After an overview of the methods used to constrain the chemical enrichment in galaxies and their environment, we discuss the observed scaling relations between metallicity and galaxy properties, the observed relative chemical abundances, how the chemical elements are distributed within galaxies, and how these properties evolve across the cosmic epochs. We discuss how the various observational findings compare with the predictions from theoretical models and numerical cosmological simulations. Finally, we briefly discuss the open problems and the prospects for major progress in this field in the nearby future.
257 citations
TL;DR: In this article, the authors investigate outflows utilizing the solar neighborhood model of the TIGRESS simulation suite, which is self-consistently treated and well resolved in the multiphase, turbulent, magnetized interstellar medium.
Abstract: Gas blown away from galactic disks by supernova (SN) feedback plays a key role in galaxy evolution. We investigate outflows utilizing the solar neighborhood model of our high-resolution, local galactic disk simulation suite, TIGRESS. In our numerical implementation, star formation and SN feedback are self-consistently treated and well resolved in the multiphase, turbulent, magnetized interstellar medium. Bursts of star formation produce spatially and temporally correlated SNe that drive strong outflows, consisting of hot (T>5x10^5K) winds and warm (5050K 1kpc from the midplane has mass and energy fluxes nearly constant with d. The hot flow escapes our local Cartesian box barely affected by gravity and is expected to accelerate up to the terminal velocity of v_wind~350-500km/s. The mean mass and energy loading factors of the hot wind are 0.1 and 0.02, respectively. For warm gas, the mean outward mass flux through d=1kpc is comparable to the mean star formation rate, but only a small fraction of this gas is at velocity >50km/s. Thus, the warm outflows eventually fall back as inflows. The warm fountain flows are created by expanding hot superbubbles at d< 1kpc; at larger d neither ram pressure acceleration nor cooling transfers significant momentum or energy flux from the hot wind to the warm outflow. The velocity distribution at launching near d~1kpc better represents warm outflows than a single mass loading factor, potentially enabling development of subgrid models for warm galactic winds in arbitrary large-scale galactic potentials.
161 citations
Cites background from "Theoretical Challenges in Galaxy Fo..."
...…2005; Heckman & Thompson 2017, for reviews) and believed to be essential to distribution of the gas and metals in galaxies and the circumgalactic/intergalactic medium (CGM/IGM) and hence to regulating cosmic star formation history (see Somerville & Davé 2015; Naab & Ostriker 2017, for reviews)....
[...]
TL;DR: The findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven life cycles that vary with the galactic environment and collectively define how galaxies form stars.
Abstract: The physics of star formation and the deposition of mass, momentum and energy into the interstellar medium by massive stars (‘feedback’) are the main uncertainties in modern cosmological simulations of galaxy formation and evolution1,2. These processes determine the properties of galaxies3,4 but are poorly understood on the scale of individual giant molecular clouds (less than 100 parsecs)5,6, which are resolved in modern galaxy formation simulations7,8. The key question is why the timescale for depleting molecular gas through star formation in galaxies (about 2 billion years)9,10 exceeds the cloud dynamical timescale by two orders of magnitude11. Either most of a cloud’s mass is converted into stars over many dynamical times12 or only a small fraction turns into stars before the cloud is dispersed on a dynamical timescale13,14. Here we report high-angular-resolution observations of the nearby flocculent spiral galaxy NGC 300. We find that the molecular gas and high-mass star formation on the scale of giant molecular clouds are spatially decorrelated, in contrast to their tight correlation on galactic scales5. We demonstrate that this decorrelation implies rapid evolutionary cycling between clouds, star formation and feedback. We apply a statistical method15,16 to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, which drives cloud dispersal on short timescales (around 1.5 million years). The rapid feedback arises from radiation and stellar winds, before supernova explosions can occur. This feedback limits cloud lifetimes to about one dynamical timescale (about 10 million years), with integrated star formation efficiencies of only 2 to 3 per cent. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven life cycles that vary with the galactic environment and collectively define how galaxies form stars. Observations that molecular gas in NGC 300 is spatially uncorrelated with high-mass stars are attributed to rapid evolution, with molecular clouds quickly destroyed by stellar feedback, and low star-formation efficiency.
151 citations
References
More filters
TL;DR: In this paper, the authors present the first cosmological results based on Planck measurements of the cosmic microwave background (CMB) temperature and lensing-potential power spectra, which are extremely well described by the standard spatially-flat six-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations.
Abstract: This paper presents the first cosmological results based on Planck measurements of the cosmic microwave background (CMB) temperature and lensing-potential power spectra. We find that the Planck spectra at high multipoles (l ≳ 40) are extremely well described by the standard spatially-flat six-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations. Within the context of this cosmology, the Planck data determine the cosmological parameters to high precision: the angular size of the sound horizon at recombination, the physical densities of baryons and cold dark matter, and the scalar spectral index are estimated to be θ∗ = (1.04147 ± 0.00062) × 10-2, Ωbh2 = 0.02205 ± 0.00028, Ωch2 = 0.1199 ± 0.0027, and ns = 0.9603 ± 0.0073, respectively(note that in this abstract we quote 68% errors on measured parameters and 95% upper limits on other parameters). For this cosmology, we find a low value of the Hubble constant, H0 = (67.3 ± 1.2) km s-1 Mpc-1, and a high value of the matter density parameter, Ωm = 0.315 ± 0.017. These values are in tension with recent direct measurements of H0 and the magnitude-redshift relation for Type Ia supernovae, but are in excellent agreement with geometrical constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent level precision using Planck CMB data alone. We use high-resolution CMB data together with Planck to provide greater control on extragalactic foreground components in an investigation of extensions to the six-parameter ΛCDM model. We present selected results from a large grid of cosmological models, using a range of additional astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured over the standard six-parameter ΛCDM cosmology. The deviation of the scalar spectral index from unity isinsensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find an upper limit of r0.002< 0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles beyond the three families of neutrinos in the standard model. Using BAO and CMB data, we find Neff = 3.30 ± 0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the sum of neutrino masses. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of Neff = 3.046. We find no evidence for dynamical dark energy; using BAO and CMB data, the dark energy equation of state parameter is constrained to be w = -1.13-0.10+0.13. We also use the Planck data to set limits on a possible variation of the fine-structure constant, dark matter annihilation and primordial magnetic fields. Despite the success of the six-parameter ΛCDM model in describing the Planck data at high multipoles, we note that this cosmology does not provide a good fit to the temperature power spectrum at low multipoles. The unusual shape of the spectrum in the multipole range 20 ≲ l ≲ 40 was seen previously in the WMAP data and is a real feature of the primordial CMB anisotropies. The poor fit to the spectrum at low multipoles is not of decisive significance, but is an “anomaly” in an otherwise self-consistent analysis of the Planck temperature data.
7,060 citations
6,206 citations
"Theoretical Challenges in Galaxy Fo..." refers methods in this paper
...The Lagrangian particle based SPH method (Gingold & Monaghan (1977), Lucy (1977), see also Springel (2010b) and Somerville & Davé (2015) for recent reviews) is relatively simple to implement and due to its adaptive spatial resolution and good conservation properties has been very popular for…...
[...]
TL;DR: GADGET-2 as mentioned in this paper is a massively parallel tree-SPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics.
Abstract: We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our implementation of SPH manifestly conserves energy and entropy in regions free of dissipation, while allowing for fully adaptive smoothing lengths. Gravitational forces are computed with a hierarchical multipole expansion, which can optionally be applied in the form of a TreePM algorithm, where only short-range forces are computed with the ‘tree’ method while long-range forces are determined with Fourier techniques. Time integration is based on a quasi-symplectic scheme where long-range and short-range forces can be integrated with different time-steps. Individual and adaptive short-range time-steps may also be employed. The domain decomposition used in the parallelization algorithm is based on a space-filling curve, resulting in high flexibility and tree force errors that do not depend on the way the domains are cut. The code is efficient in terms of memory consumption and required communication bandwidth. It has been used to compute the first cosmological N-body simulation with more than 10 10 dark matter particles, reaching a homogeneous spatial dynamic range of 10 5 per dimension in a three-dimensional box. It has also been used to carry out very large cosmological SPH simulations that account for radiative cooling and star formation, reaching total particle numbers of more than 250 million. We present the algorithms used by the code and discuss their accuracy and performance using a number of test problems. GADGET-2 is publicly released to the research community. Ke yw ords: methods: numerical ‐ galaxies: interactions ‐ dark matter.
6,196 citations
"Theoretical Challenges in Galaxy Fo..." refers methods in this paper
...The newly developed moving mesh code AREPO (Springel 2010a) similarly suffers from numerical diffusion but combines advantages of the Lagrangian and Eulerian approaches and performs much better than traditional SPH codes like GADGET on mixing problems with a high convergence rate (Sijacki et al. 2012, Springel 2010b)....
[...]
...The recent GIZMO implementation is based on the GADGET framework and shows some significant improvements on idealised test problems, in particular for low Mach number gas (Hopkins 2015)....
[...]
...…recent SPH work on cosmological galaxy formation is based on derivatives of either the GASOLINE (Wadsley, Stadel & Quinn 2004) code or the GADGET (Springel 2005) code and include updated implementations to treat the mixing problem better (see e.g. Hopkins et al. 2014; Hu et al. 2014; Read &…...
[...]
...Bottom middle: SPH (GADGET ) simulation of Aumer et al. (2013)....
[...]
...Most of the recent SPH work on cosmological galaxy formation is based on derivatives of either the GASOLINE (Wadsley, Stadel & Quinn 2004) code or the GADGET (Springel 2005) code and include updated implementations to treat the mixing problem better (see e.g. Hopkins et al. 2014; Hu et al. 2014; Read & Hayfield 2012; Schaller et al. 2015; Schaye et al. 2015; Wadsley, Veeravalli & Couchman 2008 and references therein)....
[...]
TL;DR: A finite-size particle scheme for the numerical solution of two-and three-dimensional gas dynamical problems of astronomical interest is described and tested in this article, which is then applied to the fission problem for optically thick protostars.
Abstract: A finite-size particle scheme for the numerical solution of two- and three-dimensional gas dynamical problems of astronomical interest is described and tested. The scheme is then applied to the fission problem for optically thick protostars. Results are given, showing the evolution of one such protostar from an initial state as a single, rotating star to a final state as a triple system whose components contain 60% of the original mass. The decisiveness of this numerical test of the fission hypothesis and its relevance to observed binaries are briefly discussed.
5,508 citations
"Theoretical Challenges in Galaxy Fo..." refers methods in this paper
...The Lagrangian particle based SPH method (Gingold & Monaghan (1977), Lucy (1977), see also Springel (2010b) and Somerville & Davé (2015) for recent reviews) is relatively simple to implement and due to its adaptive spatial resolution and good conservation properties has been very popular for…...
[...]
TL;DR: In this paper, the Schmidt law was used to model the global star formation law over the full range of gas densities and star formation rates observed in galaxies, and the results showed that the SFR scales with the ratio of the gas density to the average orbital timescale.
Abstract: Measurements of Hα, H I, and CO distributions in 61 normal spiral galaxies are combined with published far-infrared and CO observations of 36 infrared-selected starburst galaxies, in order to study the form of the global star formation law over the full range of gas densities and star formation rates (SFRs) observed in galaxies. The disk-averaged SFRs and gas densities for the combined sample are well represented by a Schmidt law with index N = 1.4 ± 0.15. The Schmidt law provides a surprisingly tight parametrization of the global star formation law, extending over several orders of magnitude in SFR and gas density. An alternative formulation of the star formation law, in which the SFR is presumed to scale with the ratio of the gas density to the average orbital timescale, also fits the data very well. Both descriptions provide potentially useful "recipes" for modeling the SFR in numerical simulations of galaxy formation and evolution.
5,299 citations