Theoretical Challenges in Galaxy Formation
TL;DR: In this article, a variety of plausible subresolution models were proposed to estimate the inflow to and outflow from forming galaxies because observations indicating low formation efficiency and strong circumgalactic presence of gas are persuasive.
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 because observations indicating low formation efficiency and strong circumgalactic presence of gas are persuasive. Energetic “feedback” from massive stars and accreting supermassive black holes—generally unresolved in cosmological simulations—plays a major role in driving galactic outflows, which have been shown to regulate many aspects of galaxy evolution. A surprisingly large variety of plausible subresolution models succeeds in this exercise. They capture the essential characteristics of the problem, i.e., outflows regulating galactic gas flows, but their predicti...
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: Giacconi Fellowship from the Space Telescope Science Institute; NASA through a Hubble Fellowship grant from NASA's HST-HF2-51353.001-A; NASANational Aeronautics & Space Administration (NAS5-26555); NSFNational Science Foundation (NSF) [1066293]; National Science Foundation(NSF)'s National Research Foundation (NRF) [PHY11-25915]; Munich Institute for Astro-and Particle Physics (MIAPP) of the DFG cluster of excellence 'Origin and Structure of the Universe'
Abstract: Giacconi Fellowship from the Space Telescope Science Institute; NASA through a Hubble Fellowship grant from the Space Telescope Science Institute [HST-HF2-51353.001-A]; NASANational Aeronautics & Space Administration (NASA) [NAS5-26555]; NSFNational Science Foundation (NSF) [1066293]; National Science Foundation (NSF)National Science Foundation (NSF) [PHY11-25915]; Munich Institute for Astro-and Particle Physics (MIAPP) of the DFG cluster of excellence 'Origin and Structure of the Universe; Office of Science of the U.S. Department of EnergyUnited States Department of Energy (DOE) [DE-AC02-05CH11231]; DOEUnited States Department of Energy (DOE) [DE-AC02-76SF00515]; NASA High-EndComputing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center; PRACE [012060963]; Alfred P. Sloan FoundationAlfred P. Sloan Foundation; U.S. Department of Energy Office of ScienceUnited States Department of Energy (DOE); University of Arizona; Brazilian Participation Group; Brookhaven National LaboratoryUnited States Department of Energy (DOE); Carnegie Mellon University; French Participation Group; German Participation Group; Harvard University; Instituto de Astrofisica de Canarias; Michigan State/Notre Dame/JINA Participation Group; Johns Hopkins UniversityJohns Hopkins University; Lawrence Berkeley National LaboratoryUnited States Department of Energy (DOE); Max Planck Institute for Astrophysics; Max Planck Institute for Extraterrestrial Physics; Pennsylvania State University; Princeton UniversityPrinceton University; Spanish Participation Group; Yale University; University of FloridaUniversity of Florida; NewMexico State University; New York University; Ohio State UniversityOhio State University; University of Portsmouth; University of Tokyo; University of Utah; Vanderbilt University; University of Virginia; University of WashingtonUniversity of Washington
612 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
01 Jan 2020
TL;DR: Cosmological simulations of galaxy formation have been instrumental in advancing our understanding of structure and galaxy formation in the Universe as discussed by the authors, and have also proven useful to study alternative cosmological models and their impact on the galaxy population.
Abstract: Over recent decades, cosmological simulations of galaxy formation have been instrumental in advancing our understanding of structure and galaxy formation in the Universe. These simulations follow the nonlinear evolution of galaxies, modelling a variety of physical processes over an enormous range of time and length scales. A better understanding of the relevant physical processes, improved numerical methods and increased computing power have led to simulations that can reproduce a large number of the observed galaxy properties. Modern simulations model dark matter, dark energy and ordinary matter in an expanding space-time starting from well-defined initial conditions. The modelling of ordinary matter is most challenging due to the large array of physical processes affecting this component. Cosmological simulations have also proven useful to study alternative cosmological models and their impact on the galaxy population. This Technical Review presents a concise overview of the methodology of cosmological simulations of galaxy formation and their different applications. Cosmological computer simulations of galaxy formation emerged as the primary tool to study structure formation in the Universe. This Technical Review describes the main techniques and ingredients of such simulations and their application to develop and constrain galaxy formation theories.
279 citations
TL;DR: The EMERGE model as mentioned in this paper describes the evolution of individual galaxies in large volumes from z similar to 10 to the present day, and assigns a star formation rate to each dark matter halo based on its growth rate, and the instantaneous baryon conversion efficiency, which determines how efficiently this material is converted to stars.
Abstract: We present EMERGE, an Empirical ModEl for the foRmation of GalaxiEs, describing the evolution of individual galaxies in large volumes from z similar to 10 to the present day. We assign a star formation rate to each dark matter halo based on its growth rate, which specifies how much baryonic material becomes available, and the instantaneous baryon conversion efficiency, which determines how efficiently this material is converted to stars, thereby capturing the baryonic physics. Satellites are quenched following the delayed-then-rapid model, and they are tidally disrupted once their subhalo has lost a significant fraction of its mass. The model is constrained with observed data extending out to high redshift. The empirical relations are very flexible, and the model complexity is increased only if required by the data, assessed by several model selection statistics. We find that for the same final halo mass galaxies can have very different star formation histories. Galaxies that are quenched at z = 0 typically have a higher peak star formation rate compared to their star-forming counterparts. EMERGE predicts stellar-to-halo mass ratios for individual galaxies and introduces scatter self-consistently. We find that at fixed halo mass, passive galaxies have a higher stellar mass on average. The intracluster mass in massive haloes can be up to eight times larger than the mass of the central galaxy. Clustering for star-forming and quenched galaxies is in good agreement with observational constraints, indicating a realistic assignment of galaxies to haloes.
274 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
TL;DR: In this article, a simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, BH 2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) was proposed to fit the three-year WMAP temperature and polarization data.
Abstract: A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.
6,295 citations
6,206 citations
TL;DR: In this paper, the authors present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra, which are extremely well described by the standard spatially-flat six-parameter LCDM cosmology.
Abstract: We present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra. The Planck spectra at high multipoles are extremely well described by the standard spatially-flat six-parameter LCDM cosmology. In this model Planck data determine the cosmological parameters to high precision. We find a low value of the Hubble constant, H0=67.3+/-1.2 km/s/Mpc and a high value of the matter density parameter, Omega_m=0.315+/-0.017 (+/-1 sigma errors) in excellent agreement with 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 present results from an analysis of extensions to the standard cosmology, using astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured significantly over standard LCDM. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find a 95% upper limit of r<0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles. Using BAO and CMB data, we find N_eff=3.30+/-0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the summed neutrino mass. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of N_eff=3.046. We find no evidence for dynamical dark energy. Despite the success of the standard LCDM model, this cosmology does not provide a good fit to the CMB power spectrum at low multipoles, as noted previously by the WMAP team. While not of decisive significance, this is an anomaly in an otherwise self-consistent analysis of the Planck temperature data.
6,201 citations
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