ABSTRACT We examine the scalings of galactic outflows with halo mass across a suite of 20 high-resolution cosmological zoom galaxy simulations covering halo masses in the range . These simulations ...self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables, including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as , with an amplitude and shape that are invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half of the outflow mass across all galaxy masses is later reaccreted. The recycling timescale is typically ∼1 Gyr, virtually independent of halo mass. Recycled material is reaccreted farther out in the disk and with typically ∼2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most cosmic star formation occurs.
We use maximum entropy arguments to derive the phase-space distribution of a virialized dark matter halo. Our distribution function gives an improved representation of the end product of violent ...relaxation. This is achieved by incorporating physically motivated dynamical constraints (specifically on orbital actions) which prevent arbitrary redistribution of energy.
We compare the predictions with three high-resolution dark matter simulations of widely varying mass. The numerical distribution function is accurately predicted by our argument, producing an excellent match for the vast majority of particles.
The remaining particles constitute the central cusp of the halo ( 4 per cent of the dark matter). They can be accounted for within the presented framework once the short dynamical time-scales of the centre are taken into account.
Cold dark matter: Controversies on small scales Weinberg, David H; James S. Bullock; Fabio Governato ...
Proceedings of the National Academy of Sciences - PNAS,
10/2015, Letnik:
112, Številka:
40
Journal Article
Recenzirano
Odprti dostop
The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations ...that probe the innermost regions of dark matter halos and the properties of the Milky Way’s dwarf galaxy satellites. We review the current observational and theoretical status of these “small-scale controversies.” Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: Recent numerical simulations and analytical models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies, where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low-mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.
We present a self-consistent prediction from a large-scale cosmological simulation for the population of "wandering" supermassive black holes (SMBHs) of mass greater than 106 M on long-lived, ...kpc-scale orbits within Milky Way (MW)-mass galaxies. We extract a sample of MW-mass halos from the Romulus25 cosmological simulation, which is uniquely able to capture the orbital evolution of SMBHs during and following galaxy mergers. We predict that such halos, regardless of recent merger history or morphology, host an average of 5.1 3.3 SMBHs, including their central black hole, within 10 kpc from the galactic center and an average of 12.2 8.4 SMBHs total within their virial radius, not counting those in satellite halos. Wandering SMBHs exist within their host galaxies for several Gyr, often accreted by their host halo in the early Universe. We find, with >4 significance, that wandering SMBHs are preferentially found outside of galactic disks.
We present an analysis of star formation and feedback recipes appropriate for galactic smoothed particle hydrodynamics simulations. Using an isolated Milky Way-like galaxy, we constrain these recipes ...based on well-established observational results. Our star formation recipe is based on that of Katz with the additional inclusion of physically motivated supernova feedback recipes. We propose a new feedback recipe in which Type II supernovae are modelled using an analytical treatment of blastwaves. With this feedback mechanism and a tuning of other star formation parameters, the star formation in our isolated Milky Way-like galaxy follows the slope and normalization of the observed Schmidt law. In addition, we reproduce the low-density cut-off and filamentary structure of star formation observed in disc galaxies. Our final recipe will enable better comparison of N-body simulations with observations.
We present a new zoom-in hydrodynamical simulation, ‘ErisBH’, which features the same initial conditions, resolution, and sub-grid physics as the close Milky Way-analogue ‘Eris’ (Guedes et al. 2011), ...but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. This enables a detailed study of black hole evolution and the impact of active galactic nuclei (AGN) feedback in a late-type galaxy. At z = 0, the main galaxy of ErisBH hosts a central black hole of 2.6 × 106 M⊙, which correlates to the bulge mass and the galaxy's central velocity dispersion similarly to what is observed in the Milky Way and in pseudobulges. During its evolution, the black hole grows mostly through mergers with black holes brought in by accreted satellite galaxies and very little by gas accretion (due to the modest amount of gas that reaches the central regions). AGN feedback is weak and it affects only the central
$1\text{--}2 \,\rm {kpc}$
. Yet, it limits the growth of the bulge, which results in a rotation curve that, in the inner ∼ 10 kpc, is flatter than that of Eris. We find that ErisBH is more prone to instabilities than Eris, due to its smaller bulge and larger disc. At z ∼ 0.3, an initially small bar grows to be of a few disc scalelengths in size. The formation of the bar causes a small burst of star formation in the inner few hundred pc, provides new gas to the central black hole and causes the bulge to have a boxy/peanut morphology by z = 0.
We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ...ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) time-scales, temporal behaviour, and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability time-scales, except during the merger proper, lasting ∼0.2–0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus, because of the long time spent in this phase. (iv) The time-scale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data.
Abstract
We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological ...simulations, including explicit baryonic physics treatment. A novel implementation of supermassive black hole (SMBH) formation and evolution is used, as in Tremmel et al., allowing us to explicitly follow the SMBH dynamics at the centre of galaxies. A high SIDM constant cross-section is chosen, σ = 10 cm2gr−1, to amplify differences from CDM models. Milky Way-like galaxies form a shallower DM density profile in SIDM than they do in cold dark matter (CDM), with differences already at 20 kpc scales. This demonstrates that even for the most massive spirals, the effect of SIDM dominates over the adiabatic contraction due to baryons. Strikingly, the dynamics of SMBHs differs in the SIDM and reference CDM case. SMBHs in massive spirals have sunk to the centre of their host galaxy in both the SIDM and CDM run, while in less massive galaxies about 80 per cent of the SMBH population is off-centred in the SIDM case, as opposed to the CDM case in which ∼ 90 per cent of SMBHs have reached their host’s centre. SMBHs are found as far as ∼9 kpc away from the centre of their host SIDM galaxy. This difference is due to the increased dynamical friction time-scale caused by the lower DM density in SIDM galaxies compared to CDM, resulting in core stalling. This pilot work highlights the importance of simulating in a full hydrodynamical context different DM models combined to the SMBH physics to study their influence on galaxy formation.
We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction (‘stochastic’ phase), ...during the ‘merger’ proper, lasting ∼0.2–0.3 Gyr, and in the ‘remnant’ phase. We calculate the bivariate distribution of SFR and BHAR and define the regions in the SFR–BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star-forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bivariate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances.
We make a direct comparison of the derived dark matter (DM) distributions between hydrodynamical simulations of dwarf galaxies assuming a Delta *LCDM cosmology and the observed dwarf galaxies sample ...from the THINGS survey in terms of (1) the rotation curve shape and (2) the logarithmic inner density slope Delta *a of mass density profiles. The simulations, which include the effect of baryonic feedback processes, such as gas cooling, star formation, cosmic UV background heating, and most importantly, physically motivated gas outflows driven by supernovae, form bulgeless galaxies with DM cores. We show that the stellar and baryonic mass is similar to that inferred from photometric and kinematic methods for galaxies of similar circular velocity. Analyzing the simulations in exactly the same way as the observational sample allows us to address directly the so-called cusp/core problem in the Delta *LCDM model. We show that the rotation curves of the simulated dwarf galaxies rise less steeply than cold dark matter rotation curves and are consistent with those of the THINGS dwarf galaxies. The mean value of the logarithmic inner density slopes Delta *a of the simulated galaxies' DM density profiles is ~--0.4 ? 0.1, which shows good agreement with Delta *a = --0.29 ? 0.07 of the THINGS dwarf galaxies. The effect of non-circular motions is not significant enough to affect the results. This confirms that the baryonic feedback processes included in the simulations are efficiently able to make the initial cusps with Delta *a ~--1.0 to --1.5 predicted by DM-only simulations shallower and induce DM halos with a central mass distribution similar to that observed in nearby dwarf galaxies.