scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Cosmological simulations of the formation of the stellar haloes around disc galaxies

TL;DR: In this article, the formation of stellar spheroids of the Milky Way mass disc galaxies is studied using the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations.
Abstract: We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ≈400 simulated ∼L∗ disc galaxies. The simulated galaxies are surrounded by low-mass, low surface brightness stellar haloes that extend out to ∼100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in the slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalization of the spheroid density profile within r 30 kpc (e.g. when fitted by a spherical power-law profile, the indices range from −2.6 to −3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically, older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the long-standing failure of pure accretion-based models to reproduce the observed properties of the inner spheroid.
Citations
More filters
Journal ArticleDOI
TL;DR: In this article, a multi-epoch abundance matching (MEAM) model was proposed to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z � 4 to the present.
Abstract: We present a new statistical method to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z � 4 to the present. This multi-epoch abundance matching (MEAM) model self-consistently takes into account that satellite galaxies first become satellites at times earlier than they are observed. We employ a redshift-dependent parameterization of the stellar-to-halo mass relation to populate haloes and subhaloes in the Millennium simulations with galaxies, requiring that the observed stellar mass functions at several redshifts be reproduced simultaneously. We show that physically meaningful growth of massive galaxies is consistent with these data only if observational mass errors are taken into account. Using merger trees extracted from the dark matter simulations in combination with MEAM, we predict the average assembly histories of galaxies, separating into star formation within the galaxies (in-situ) and accretion of stars (ex-situ). Our main results are: The peak star formation efficiency decreases with redshift from 23 per cent at z = 0 to 9 per cent at z = 4 while the corresponding halo mass increases from 10 11.8 M⊙ to 10 12.5 M⊙. The star formation rate of central galaxies peaks at a redshift which depends on halo mass; for massive haloes this peak is at early cosmic times while for low-mass galaxies the peak has not been reached yet. In haloes similar to that of the Milky-Way about half of the central stellar mass is assembled after z = 0.7. In low-mass haloes, the accretion of satellites contributes little to the assembly of their central galaxies, while in massive haloes more than half of the central stellar mass is formed ex-situ with significant accretion of satellites at z < 2. We find that our method implies a cosmic star formation history and an evolution of specific star formation rates which are consistent with those inferred directly. We present convenient fitting functions for stellar masses, star formation rates, and accretion rates as functions of halo mass and redshift.

1,327 citations


Cites background from "Cosmological simulations of the for..."

  • ...In lower-mass disk-like galaxies the accreted stars contribute to the build-up of the outer stellar halo and can also directly contribute on the order of 10 percent to the stel- lar mass of the disk (e.g. Abadi et al. 2006; Font et al. 2011; Tissera et al. 2012)....

    [...]

Journal ArticleDOI
TL;DR: In this paper, the authors review the key integrated, structural and kinematic parameters of the Galaxy, and point to uncertainties as well as directions for future progress, and show that the Galaxy is a luminous (L⋆) barred spiral with a central box/peanut bulge, a dominant disk, and a diffuse stellar halo.
Abstract: Our Galaxy, the Milky Way, is a benchmark for understanding disk galaxies. It is the only galaxy whose formation history can be studied using the full distribution of stars from faint dwarfs to supergiants. The oldest components provide us with unique insight into how galaxies form and evolve over billions of years. The Galaxy is a luminous (L⋆) barred spiral with a central box/peanut bulge, a dominant disk, and a diffuse stellar halo. Based on global properties, it falls in the sparsely populated “green valley” region of the galaxy color-magnitude diagram. Here we review the key integrated, structural and kinematic parameters of the Galaxy, and point to uncertainties as well as directions for future progress. Galactic studies will continue to play a fundamental role far into the future because there are measurements that can only be made in the near field and much of contemporary astrophysics depends on such observations.

1,084 citations

Journal ArticleDOI
TL;DR: In this article, the authors compare the results of various cosmological gas-dynamical codes used to simulate the formation of a galaxy in the Λ cold dark matter structure formation paradigm.
Abstract: We compare the results of various cosmological gas-dynamical codes used to simulate the formation of a galaxy in the Λ cold dark matter structure formation paradigm. The various runs (13 in total) differ in their numerical hydrodynamical treatment [smoothed particle hydrodynamics (SPH), moving mesh and adaptive mesh refinement] but share the same initial conditions and adopt in each case their latest published model of gas cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z= 0, due mainly to the different implementations of star formation and feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disc is discernible in most simulations, although its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observed spirals. In general, models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disc component. Although numerical convergence is not particularly good for any of the codes, our conclusions hold at two different numerical resolutions. Some differences can also be traced to the different numerical techniques; for example, more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high angular momentum stellar discs remains a challenge.

502 citations


Cites methods from "Cosmological simulations of the for..."

  • ...A full description of the model can be found in Crain et al. (2009) and Schaye et al. (2010) and further applications of this model can be found in Crain et al. (2010), Cui et al. (2011), Deason et al. (2011), Font et al. (2011) and McCarthy et al. (2011)....

    [...]

Journal ArticleDOI
TL;DR: The Subaru Prime Focus Spectrograph (PFS) as discussed by the authors is a massively multiplexed fiber-fed optical and near-infrared three-arm spectrograph with a large field of view.
Abstract: The Subaru Prime Focus Spectrograph (PFS) is a massively multiplexed fiber-fed optical and near-infrared three-arm spectrograph (N_fiber = 2400, 380 ≤ λ ≤ 1260 nm, 1 $_{.}^{\circ}$3 diameter field of view). Here, we summarize the science cases in terms of provisional plans for a 300-night Subaru survey. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 h^−3 Gpc^3 in the redshift range 0.8 < z < 2.4. In each of six redshift bins, the cosmological distances will be measured to 3% precision via the baryonic acoustic oscillation scale, and redshift-space distortion measures will constrain structure growth to 6% precision. In the near-field cosmology program, radial velocities and chemical abundances of stars in the Milky Way and M 31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo, and tidal streams as faint as V ∼ 22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple α-element abundances and more precise velocities for Galactic stars. For the galaxy evolution program, our simulations suggest the wide wavelength range of PFS will be powerful in probing the galaxy population and its clustering over a wide redshift range. We plan to conduct a color-selected survey of 1 < z < 2 galaxies and AGN over 16 deg^2 to J ≃ 23.4, yielding a fair sample of galaxies with stellar masses above ∼10^10 M_⊙ at z ≃ 2. A two-tiered survey of higher redshift Lyman break galaxies and Lyman alpha emitters will quantify the properties of early systems close to the reionization epoch.

495 citations

Journal ArticleDOI
TL;DR: The Subaru Prime Focus Spectrograph (PFS) as mentioned in this paper is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph, offering unique opportunities in survey astronomy.
Abstract: The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380

478 citations


Cites background from "Cosmological simulations of the for..."

  • ...All of these processes are predicted to be controlled by hierarchical clustering of CDM via self-gravity on galactic and sub-galactic scales (e.g., Bullock & Johnston 2005; Cooper et al. 2010; Font et al. 2011)....

    [...]

References
More filters
Journal ArticleDOI
TL;DR: In this article, the spectral evolution of stellar populations at ages between 100,000 yr and 20 Gyr at a resolution of 3 A across the whole wavelength range from 3200 to 9500 A for a wide range of metallicities.
Abstract: We present a new model for computing the spectral evolution of stellar populations at ages between 100,000 yr and 20 Gyr at a resolution of 3 A across the whole wavelength range from 3200 to 9500 A for a wide range of metallicities. These predictions are based on a newly available library of observed stellar spectra. We also compute the spectral evolution across a larger wavelength range, from 91 A to 160 micron, at lower resolution. The model incorporates recent progress in stellar evolution theory and an observationally motivated prescription for thermally-pulsing stars on the asymptotic giant branch. The latter is supported by observations of surface brightness fluctuations in nearby stellar populations. We show that this model reproduces well the observed optical and near-infrared colour-magnitude diagrams of Galactic star clusters of various ages and metallicities. Stochastic fluctuations in the numbers of stars in different evolutionary phases can account for the full range of observed integrated colours of star clusters in the Magellanic Clouds. The model reproduces in detail typical galaxy spectra from the Early Data Release (EDR) of the Sloan Digital Sky Survey (SDSS). We exemplify how this type of spectral fit can constrain physical parameters such as the star formation history, metallicity and dust content of galaxies. Our model is the first to enable accurate studies of absorption-line strengths in galaxies containing stars over the full range of ages. Using the highest-quality spectra of the SDSS EDR, we show that this model can reproduce simultaneously the observed strengths of those Lick indices that do not depend strongly on element abundance ratios [abridged].

10,384 citations

Journal ArticleDOI
TL;DR: A review of the present-day mass function and initial mass function in various components of the Galaxy (disk, spheroid, young, and globular clusters) and in conditions characteristic of early star formation is presented in this paper.
Abstract: We review recent determinations of the present-day mass function (PDMF) and initial mass function (IMF) in various components of the Galaxy—disk, spheroid, young, and globular clusters—and in conditions characteristic of early star formation. As a general feature, the IMF is found to depend weakly on the environment and to be well described by a power-law form forM , and a lognormal form below, except possibly for m!1 early star formation conditions. The disk IMF for single objects has a characteristic mass around M , m!0.08 c and a variance in logarithmic mass , whereas the IMF for multiple systems hasM , and . j!0.7 m!0.2 j!0.6 c The extension of the single MF into the brown dwarf regime is in good agreement with present estimates of L- and T-dwarf densities and yields a disk brown dwarf number density comparable to the stellar one, n!n! BD " pc !3 .T he IMF of young clusters is found to be consistent with the disk fi eld IMF, providing the same correction 0.1 for unresolved binaries, confirming the fact that young star clusters and disk field stars represent the same stellar population. Dynamical effects, yielding depletion of the lowest mass objects, are found to become consequential for ages!130 Myr. The spheroid IMF relies on much less robust grounds. The large metallicity spread in the local subdwarf photometric sample, in particular, remains puzzling. Recent observations suggest that there is a continuous kinematic shear between the thick-disk population, present in local samples, and the genuine spheroid one. This enables us to derive only an upper limit for the spheroid mass density and IMF. Within all the uncertainties, the latter is found to be similar to the one derived for globular clusters and is well represented also by a lognormal form with a characteristic mass slightly larger than for the disk, M , ,e xcluding as ignif icant population of m!0.2-0.3 c brown dwarfs in globular clusters and in the spheroid. The IMF characteristic of early star formation at large redshift remains undetermined, but different observational constraints suggest that it does not extend below!1M , .T hese results suggest a characteristic mass for star formation that decreases with time, from conditions prevailing at large redshift to conditions characteristic of the spheroid (or thick disk) to present-day conditions.Theseconclusions,however, remain speculative, given the large uncertainties in the spheroid and early star IMF determinations. These IMFs allow a reasonably robust determination of the Galactic present-day and initial stellar and brown dwarf contents. They also have important galactic implications beyond the Milky Way in yielding more accurate mass-to-light ratio determinations. The mass-to-light ratios obtained with the disk and the spheroid IMF yield values 1.8-1.4 times smaller than for a Salpeter IMF, respectively, in agreement with various recent dynamical determinations. This general IMF determination is examined in the context of star formation theory. None of the theories based on a Jeans-type mechanism, where fragmentation is due only to gravity, can fulfill all the observational constraints on star formation and predict a large number of substellar objects. On the other hand, recent numerical simulations of compressible turbulence, in particular in super-Alfvenic conditions, seem to reproduce both qualitatively and quantitatively the stellar and substellar IMF and thus provide an appealing theoretical foundation. In this picture, star formation is induced by the dissipation of large-scale turbulence to smaller scales through radiative MHD shocks, producing filamentary structures. These shocks produce local nonequilibrium structures with large density contrasts, which collapse eventually in gravitationally bound objects under the combined influence of turbulence and gravity. The concept of a single Jeans mass is replaced by a distribution of local Jeans masses, representative of the lognormal probability density function of the turbulent gas. Objects below the mean thermal Jeans mass still have a possibility to collapse, although with a decreasing probability.

8,218 citations


"Cosmological simulations of the for..." refers methods in this paper

  • ...This corresponds to using approximately 80% of the total energy available from supernovae for a Chabrier (2003) IMF, which is assumed in the simulation....

    [...]

Journal ArticleDOI
TL;DR: In this article, high-resolution N-body simulations show that the density profiles of dark matter halos formed in the standard CDM cosmogony can be fit accurately by scaling a simple universal profile.
Abstract: High resolution N-body simulations show that the density profiles of dark matter halos formed in the standard CDM cosmogony can be fit accurately by scaling a simple “universal” profile. Regardless of their mass, halos are nearly isothermal over a large range in radius, but significantly shallower than r -2 near the center and steeper than r -2 in the outer regions. The characteristic overdensity of a halo correlates strongly with halo mass in a manner consistent with the mass dependence of the epoch of halo formation. Matching the shape of the rotation curves of disk galaxies with this halo structure requires (i) disk mass-to-light ratios to increase systematically with luminosity, (ii) halo circular velocities to be systematically lower than the disk rotation speed, and (iii) that the masses of halos surrounding bright galaxies depend only weakly on galaxy luminosity. This offers an attractive explanation for the puzzling lack of correlation between luminosity and dynamics in observed samples of binary galaxies and of satellite companions of bright spiral galaxies, suggesting that the structure of dark matter halos surrounding bright spirals is similar to that of cold dark matter halos.

7,622 citations

Journal ArticleDOI
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


"Cosmological simulations of the for..." refers methods in this paper

  • ...The simulations were evolved to z = 0 using the TreePM-SPH code gadget-3 (last described in Springel 2005)....

    [...]

Journal ArticleDOI
TL;DR: In this article, the Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data were used to constrain the physics of cosmic inflation via Gaussianity, adiabaticity, the power spectrum of primordial fluctuations, gravitational waves, and spatial curvature.
Abstract: The Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data provide stringent limits on deviations from the minimal, six-parameter Λ cold dark matter model. We report these limits and use them to constrain the physics of cosmic inflation via Gaussianity, adiabaticity, the power spectrum of primordial fluctuations, gravitational waves, and spatial curvature. We also constrain models of dark energy via its equation of state, parity-violating interaction, and neutrino properties, such as mass and the number of species. We detect no convincing deviations from the minimal model. The six parameters and the corresponding 68% uncertainties, derived from the WMAP data combined with the distance measurements from the Type Ia supernovae (SN) and the Baryon Acoustic Oscillations (BAO) in the distribution of galaxies, are: Ω b h 2 = 0.02267+0.00058 –0.00059, Ω c h 2 = 0.1131 ± 0.0034, ΩΛ = 0.726 ± 0.015, ns = 0.960 ± 0.013, τ = 0.084 ± 0.016, and at k = 0.002 Mpc-1. From these, we derive σ8 = 0.812 ± 0.026, H 0 = 70.5 ± 1.3 km s-1 Mpc–1, Ω b = 0.0456 ± 0.0015, Ω c = 0.228 ± 0.013, Ω m h 2 = 0.1358+0.0037 –0.0036, z reion = 10.9 ± 1.4, and t 0 = 13.72 ± 0.12 Gyr. With the WMAP data combined with BAO and SN, we find the limit on the tensor-to-scalar ratio of r 1 is disfavored even when gravitational waves are included, which constrains the models of inflation that can produce significant gravitational waves, such as chaotic or power-law inflation models, or a blue spectrum, such as hybrid inflation models. We obtain tight, simultaneous limits on the (constant) equation of state of dark energy and the spatial curvature of the universe: –0.14 < 1 + w < 0.12(95%CL) and –0.0179 < Ω k < 0.0081(95%CL). We provide a set of WMAP distance priors, to test a variety of dark energy models with spatial curvature. We test a time-dependent w with a present value constrained as –0.33 < 1 + w 0 < 0.21 (95% CL). Temperature and dark matter fluctuations are found to obey the adiabatic relation to within 8.9% and 2.1% for the axion-type and curvaton-type dark matter, respectively. The power spectra of TB and EB correlations constrain a parity-violating interaction, which rotates the polarization angle and converts E to B. The polarization angle could not be rotated more than –59 < Δα < 24 (95% CL) between the decoupling and the present epoch. We find the limit on the total mass of massive neutrinos of ∑m ν < 0.67 eV(95%CL), which is free from the uncertainty in the normalization of the large-scale structure data. The number of relativistic degrees of freedom (dof), expressed in units of the effective number of neutrino species, is constrained as N eff = 4.4 ± 1.5 (68%), consistent with the standard value of 3.04. Finally, quantitative limits on physically-motivated primordial non-Gaussianity parameters are –9 < f local NL < 111 (95% CL) and –151 < f equil NL < 253 (95% CL) for the local and equilateral models, respectively.

5,904 citations


"Cosmological simulations of the for..." refers background in this paper

  • ...The value of σ8 is approximately 2-sigma higher than has been inferred from the most recent CMB data (Komatsu et al. 2009), which will affect the abundance of Milky Way-mass systems somewhat, but should not significantly affect their individual properties....

    [...]