The Role of Cold Flows in the Assembly of Galaxy Disks

We use high resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy as it enters the virial radius. Low mass galaxies are instead dominated by accretion of gas that stays well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation in the disk at all times to the present. Even for galaxies at higher masses, cold flows dominate the growth of the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although star formation in the disk is mitigated by supernovae feedback, the short cooling times allow for the growth of stellar disks at higher redshifts than predicted by the standard model.Comment: accepted to Ap

MaGICC baryon cycle: The enrichment history of simulated disc galaxies

Using cosmological galaxy formation simulations from the MaGICC (Making Galaxies in a Cosmological Context) project, spanning stellar mass from ∼107 to 3 × 1010 M⊙, we trace the baryonic cycle of infalling gas from the virial radius through to its eventual participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and circumgalactic medium. We derive the distributions of time between gas crossing the virial radius and being accreted to the star-forming region (which allows for mixing within the corona), as well as the time between gas being accreted to the star-forming region and then ultimately forming stars (which allows for mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star-forming region. Gas entering high-mass galaxies is pre-enriched in low-mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z ∼ 5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star-forming region, with gas that is feeding the star-forming region surpassing 0.1 Z⊙ by z = 0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulated galaxy (Mvir ∼ 2 × 1010 M⊙, with M⋆ ∼ 107 M⊙), by contrast, accretes primordial gas through the virial radius and on to the disc, throughout its history. Much like the case for classical analytical solutions to the so-called ‘G-dwarf problem’, overproduction of low-metallicity stars is ameliorated by the interplay between the time of accretion on to the disc and the subsequent involvement in star formation – i.e. due to the inefficiency of star formation. Finally, gas outflow/metal removal rates from star-forming regions as a function of galactic mass are presented

NIHAO XI: Formation of Ultra-Diffuse Galaxies by outflows

We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses $10^{10-11}\rm M_{\odot}$, have stellar masses of $10^{7-8.5}\rm M_{\odot}$, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colors, an average S\'ersic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible HI gas mass of $10^{7-9}\rm M_{\odot}$, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes that form UDGs: feedback driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more HI gas than less extended dwarfs of similar $\rm M_{\star}$.Comment: matches accepted version, MNRAS Letter 2016-10-1

The role of feedback in shaping the structure of the interstellar medium

We present an analysis of the role of feedback in shaping the neutral hydrogen (H I) content of simulated disc galaxies. For our analysis, we have used two realizations of two separate Milky Way-like (similar to L star) discs - one employing a conservative feedback scheme (McMaster Unbiased Galaxy Survey), the other significantly more energetic [Making Galaxies In a Cosmological Context (MaGICC)]. To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred H I distribution; construct power spectra associated with the underlying structure of the simulated cold interstellar medium, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from The H I Nearby Galaxy Survey (THINGS). Single power-law fits (P proportional to k(gamma)) to the power spectra of the stronger feedback (MaGICC) runs (over spatial scales corresponding to similar to 0.5 to similar to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma similar to -2.5). The weaker feedback (MUGS) runs exhibit shallower power-law slopes (gamma similar to -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e. gamma similar to -1.4 for scales in excess of similar to 2 kpc) and a steeper slope on scales below similar to 1 kpc (gamma similar to -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial H I surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs

The stellar metallicity distribution of disc galaxies and bulges in cosmological simulations

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.Comment: 20 pages, 14 figures, MNRAS, accepte

ECOLOGICAL-ECONOMIC MODELS OF SUSTAINABLE HARVEST FOR AN ENDANGERED BUT EXOTIC MEGAHERBIVORE IN NORTHERN AUSTRALIA

ABSTRACT. How can one manage wildlife under a suite of competing values? In isolation, the ecological economics of native wildlife harvest, threatened species conservation and control of exotic species are all well established sub-disciplines of wildlife management. However, the wild banteng (Bos javanicus) population of northern Australia represents an interesting combination of these aspirations. A native bovid of Southeast Asia now 'endangered' in its native range, banteng were introduced into northern Australia in 1849. Today, a population of 8,000 10,000 resides on one small, isolated peninsula in western Arnhem Land, Northern Territory and is harvested by both recreational (trophy) and aboriginal subsistence hunters. Indigenous, industry and conservationist stakeholders differ in their requirements for population management. Here we analyze the ecological and economic costs/benefits of a series of potential harvest management options for Australia's banteng population, with the aim being either to: (1) maximize sustainable yield (MSY); (2) maximize harvest of trophy males; (3) maximize indigenous off-take; (4) suppress density or completely eradicate the population; (5) minimize risk of extinction whilst limiting range expansion; (6) scenarios incorporating two or more of options 1 5. The modeling framework employed stochastic, density-regulated matrix population models with life-history parameters derived from (i) allometric relationships (for estimating r max , generation length, fecundity and densities for a banteng-sized mammal) and (ii) measured For each management option, we present a simple economic analysis that incorporates estimated costs of management implementation and associated profits projected. Results demonstrate that revenue of >AU$200,000 is possible from meat production and safari hunting without compromising long-term population stability or the conservation status of this endangered bovid