ABSTRACT
We present a model for the formation of the first galaxies before and after the reionization of hydrogen in the early universe. In this model, galaxy formation can only take place in dark ...matter haloes whose mass exceeds a redshift-dependent critical value, which, before reionization, is equal (in the simplest case) to the mass at which atomic hydrogen cooling becomes effective and, after reionization, is equal to the mass above which gas cannot remain in hydrostatic equilibrium. We define the Halo Occupation Fraction (HOF) as the fraction of haloes that host a luminous galaxy as a function of halo mass. The HOF is established by the interplay between the evolution of the critical mass and the assembly history of haloes and depends on three factors: the minimum halo mass for galaxy formation before reionization, the redshift of reionization, and the intensity of the (evolving) external photoheating rate. Our fiducial model predicts a cutoff in the galaxy mass function at a present-day halo mass, $M_{200} \sim 3\times 10^{8} \, \mathrm{M}_{\odot }$; 100 per cent occupation at $M_{200} \gt 5\times 10^9 \, \mathrm{M}_{\odot }$; and a population of starless gaseous haloes of present-day mass in the range 106 ≲ M200/M⊙ ≲ 5 × 109, in which the gas is in thermal equilibrium with the ultraviolet background radiation and in hydrostatic equilibrium in the gravitational potential of the halo. The transition between HOF = 0 and HOF = 1 reflects the stochastic nature of halo mass growth. We explore how these characteristic masses vary with model assumptions and parameter values. The results of our model are in excellent agreement with cosmological hydrodynamic simulations of galaxy formation.
ABSTRACT
We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the eagle model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density ...threshold (ρth) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. LT runs show little effect on the inner DM profile, even in systems with extended and bursty SF, two ingredients often cited as critical for core formation. For higher thresholds, gas is able to dominate the gravitational potential before being ejected by feedback. This can lead to a substantial reduction in the inner DM content, but only if the gas is gravitationally important over an extended period of time, allowing the halo to contract before gas removal. Rapid assembly and removal of gas in short SF bursts is less effective at altering the inner DM content. Subsequent gas accretion may draw DM back in and reform a cusp, unless SF is bursty enough to prevent it, preserving the core. Thus, for the eagle SF + feedback model, there is no simple relation between core formation and SF history, contrary to recent claims. The dependence of the inner DM content of dwarfs on ρth hinders robust predictions and the interpretation of observations. A simulation of a $(12 \rm \ Mpc)^3$ volume with high ρth results in dwarfs with sizeable cores over a limited halo mass range, but with insufficient variety in mass profiles to explain the observed diversity of dwarf galaxy rotation curves.
ABSTRACT
We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical ...simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r ≲ 20 kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2, and 4 at radial distances of 20, 8, and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well.1 The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.97_{-0.19}^{+0.24}\times 10^{12}\,\mathrm{M}_\odot$, and concentration before baryon contraction of $9.4_{-2.6}^{+1.9}$, which lie close to the median halo mass–concentration relation predicted in ΛCDM. The inferred total mass, $M_{200}^{\rm total}=1.08_{-0.14}^{+0.20} \times 10^{12}\,\mathrm{M}_\odot$, is in good agreement with recent measurements. The model gives an MW stellar mass of $5.04_{-0.52}^{+0.43}\times 10^{10}\,\mathrm{M}_\odot$ and infers that the DM density at the Solar position is $\rho _{\odot }^{\rm DM}=8.8_{-0.5}^{+0.5}\times 10^{-3}\,\mathrm{M}_\odot \,\mathrm{pc}^{-3}\equiv 0.33_{-0.02}^{+0.02}\,\rm {GeV}\,\rm {cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo.
ABSTRACT
We examine the gas content of field dwarf galaxies in a high-resolution cosmological simulation. In agreement with previous work, we find that galaxies inhabiting dark matter haloes with ...mass below a critical value, M200 ≲ Mcrit ≈ 5 × 109 M⊙, are quiescent at the present day. The gas content of these galaxies is thus insensitive to feedback from evolving stars. Almost half of these quiescent systems today have gas masses much smaller than that expected for their mass. We find that gas-deficient galaxies originate from (1) past interactions with massive hosts, in which a dwarf loses gas and dark matter via tidal and ram-pressure forces; and (2) from hydrodynamic interactions with the gaseous filaments and sheets of the cosmic web, in which a dwarf loses gas via ram pressure. We refer to these systems as ‘flybys’ and ‘COSWEBs’. Flybys locate in high-density regions, tracing the location of the most massive galaxies in the simulation. In contrast, COSWEBs are dispersed throughout the volume and trace the cosmic web. For sub-critical systems, M200 < Mcrit, the fraction of COSWEB galaxies can be as high as $35 {{\ \rm per\ cent}}$, and much higher for flybys, which make up 100 per cent of the galaxies with $M_{200}\lt 3\times 10^8 \ \rm M_{\odot }$. The deficit of gas caused by these mechanisms may preclude the detection of a large fraction of field dwarfs in future H i surveys. For galaxies inhabiting haloes with mass M200 > Mcrit, we find that cosmic web stripping, on average, shuts down star formation in more than $70{{\ \rm per\ cent}}$ of the affected systems.
The core–cusp problem: a matter of perspective Genina, Anna; Benítez-Llambay, Alejandro; Frenk, Carlos S ...
Monthly Notices of the Royal Astronomical Society,
02/2018, Letnik:
474, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Abstract
The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their ...matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Peñarrubia have used this approach to argue that data for these galaxies are consistent with constant-density ‘cores’ in their inner regions and rule out ‘cuspy’ Navarro–Frenk–White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Λ cold dark matter cosmological hydrodynamic simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Peñarrubia to a sample of 50 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed.
ABSTRACT
An ionization front (I-front) that propagates through an inhomogeneous medium is slowed down by self-shielding and recombinations. We perform cosmological radiation hydrodynamics simulations ...of the I-front propagation during the epoch of cosmic reionization. The simulations resolve gas in mini-haloes (halo mass 104 ≲ MhM⊙ ≲ 108) that could dominate recombinations, in a computational volume that is large enough to sample the abundance of such haloes. The numerical resolution is sufficient (gas-particle mass ∼20 M⊙ and spatial resolution <0.1 ckpc) to allow accurate modelling of the hydrodynamic response of gas to photoheating. We quantify the photoevaporation time of mini-haloes as a function of Mh and its dependence on the photoionization rate, Γ−12, and the redshift of reionization, zi. The recombination rate can be enhanced over that of a uniform medium by a factor ∼10−20 early on. The peak value increases with Γ−12 and decreases with zi, due to the enhanced contribution from mini-haloes. The clumping factor, cr, decreases to a factor of a few at ∼100 Myr after the passage of the I-front when the mini-haloes have been photoevaporated; this asymptotic value depends only weakly on Γ−12. Recombinations increase the required number of photons per baryon to reionize the Universe by 20 per cent–100 per cent, with the higher value occurring when Γ−12 is high and zi is low. We complement the numerical simulations with simple analytical models for the evaporation rate and the inverse Strömgren layer. The study also demonstrates the proficiency and potential of sph-m1rt to address astrophysical problems in high-resolution cosmological simulations.
ABSTRACT
We use a compilation of disc galaxy rotation curves to assess the role of the luminous component (‘baryons’) in the rotation curve diversity problem. As in earlier work, we find that ...rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly rising rotation curves consistent with cuspy cold dark matter haloes; slowly rising rotation curves (characteristic of galaxies with inner mass deficits or ‘cores’) occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from ‘cuspy’ systems where baryons are unimportant in the inner mass budget to ‘cored’ galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the ‘mass discrepancy–acceleration relation’ (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modelling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter haloes remains an open issue.
Abstract
We study the vertical structure of polytropic centrifugally supported gaseous discs embedded in cold dark matter (CDM) haloes. At fixed radius, R, the shape of the vertical density profile ...depends weakly on whether the disc is self-gravitating (SG) or non-self-gravitating (NSG). The disc ‘characteristic’ thickness, zH, set by the midplane sound speed and circular velocity, zNSG = (cs/Vc)R, in the NSG case, and by the sound speed and surface density, $z_{\rm SG} = c_{\rm s}^2/G\Sigma$, in SG discs, is smaller than zSG and zNSG. SG discs are typically Toomre unstable, NSG discs are stable. Exponential discs in CDM haloes with roughly flat circular velocity curves ‘flare’ outwards. Flares in mono abundance or coeval populations in galaxies like the Milky Way are thus not necessarily due to radial migration. For the polytropic equation of state of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations, discs that match observational constraints are NSG for Md < 3 × 109 M⊙ and SG at higher masses, if fully gaseous. We test these analytic results using a set of idealized smoothed particle hydrodynamic simulations and find excellent agreement. Our results clarify the role of the gravitational softening on the thickness of simulated discs, and on the onset of radial instabilities. EAGLE low-mass discs are NSG so the softening plays no role in their vertical structure. High-mass discs are expected to be SG and unstable, and may be artificially thickened and stabilized unless gravity is well resolved. Simulations with spatial resolution high enough to not compromise the vertical structure of a disc also resolve the onset of their instabilities, but the converse is not true.
ABSTRACT
We use the APOSTLE cosmological simulations to examine the role of the cosmic UV background in regulating star formation (SF) in low-mass Lambda cold dark matter (ΛCDM) haloes. In agreement ...with earlier work, we find that after reionization SF proceeds mainly in haloes whose mass exceeds a redshift-dependent ‘critical’ mass, Mcrit, set by the structure of the haloes and by the thermal pressure of UV-heated gas. Mcrit increases from $\sim 10^{8}\, \mathrm{M}_\odot$ at z ∼ 10 to $M_{\rm crit}\sim 10^{9.7}\, \mathrm{M}_\odot$ at z = 0, roughly following the average mass growth of haloes in that mass range. This implies that haloes well above or below critical at present have remained so since early times. Haloes of luminous dwarfs today were already above-critical and star forming at high redshift, explaining naturally the ubiquitous presence of ancient stellar populations in dwarfs, regardless of luminosity. The SF history of systems close to the critical boundary is more complex. SF may cease or reignite in dwarfs whose host halo falls below or climbs above the critical boundary, suggesting an attractive explanation for the episodic nature of SF in some dwarfs. Also, some subcritical haloes today may have been above critical in the past; these systems should at present make up a sizable population of faint field dwarfs lacking ongoing star formation. Although few such galaxies are currently known, the discovery of this population would provide strong support for our results. Our work indicates that, rather than stellar feedback, it is the ionizing UV background and mass accretion history what regulates SF in the faintest dwarfs.