Kate Pilkington

Hierarchical formation of bulgeless galaxies: why outflows have low angular momentum

Using high resolution, fully cosmological smoothed particle hydrodynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how high redshift gas outflows can modify the baryon angular momentum distribution and allow pure disc galaxies to form. We outline how galactic outflows preferentially remove low angular momentum material due a combination of (a) star formation peaking at high redshift in shallow dark matter potentials, an epoch when accreted gas has relatively low angular momentum, (b) the existence of an extended reservoir of high angular momentum gas in the outer disc to provide material for prolonged SF at later times and (c) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of ‘classical’, merger driven bulges in low-mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component similar to those detected in nearby dwarfs.

Metallicity gradients in disks - Do galaxies form inside-out?

Aims. We examine radial and vertical metallicity gradients using a suite of disk galaxy hydrodynamical simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient slope and reconcile the differences existing between extant models and observations within the canonical “inside-out” disk growth paradigm. Methods. A suite of 25 cosmological disks is used to examine the evolution of metallicity gradients; this consists of 19 galaxies selected from the RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code ramses ,i ncluding eight drawn from the “field” and six from “loose group” environments. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code gasoline. Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. Results. We first show that generically flatter gradients are observed at redshift zero when comparing older stars with those forming today, consistent with expectations of kinematically hot simulations, but counter to that observed in the Milky Way. The vertical abundance gradients at ∼1−3 disk scalelengths are comparable to those observed in the thick disk of the Milky Way, but significantly shallower than those seen in the thin disk. Most importantly, we find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. Conclusions. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical “inside-out” growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z ∼ 1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

Constraining sub-grid physics with high-redshift spatially-resolved metallicity distributions

Astronomy and Astrophysics 544 (2013): A47 reproduced with permission from Astronomy & Astrophysics

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

By means of high-resolution cosmological hydrodynamical simulations of Milky Way (MW) 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 subcomponents, 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 MW. 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 MW analogue, 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.

The cold gas content of bulgeless dwarf galaxies

We present an analysis of the neutral hydrogen (H i) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al., the high-resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the H i distribution and kinematics of this simulated bulgeless disc with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g. velocity ellipsoid and turbulence) and the power spectrum of structure within the cold interstellar medium (ISM) from the simulations. The highest resolution dwarf, when using a high-density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with the observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the Small Magellanic Cloud. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within parsec-scale molecular clouds.

Decoding the message from meteoritic stardust silicon carbide grains

Micron-sized stardust grains that originated in ancient stars are recovered from meteorites and analyzed using high-resolution mass spectrometry. The most widely studied type of stardust is silicon carbide (SiC). Thousands of these grains have been analyzed with high precision for their Si isotopic composition. Here we show that the distribution of the Si isotopic composition of the vast majority of stardust SiC grains carries the imprints of a spread in the age-metallicity distribution of their parent stars and of a power-law increase of the relative formation efficiency of SiC dust with the metallicity. This result offers a solution for the long-standing problem of silicon in stardust SiC grains, confirms the necessity of coupling chemistry and dynamics in simulations of the chemical evolution of our Galaxy, and constrains the modeling of dust condensation in stellar winds as a function of the metallicity.

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 Chemical Evolution of Globular Clusters - II. Metals and Fluorine

In the first paper in this series, we proposed a new framework in which to model the chemical evolution of globular clusters. This model, is predicated upon the assumption that clusters form within an interstellar medium enriched locally by the ejecta of a single Type Ia supernova and varying numbers of asymptotic giant branch stars, superimposed on an ambient medium pre-enriched by low-metallicity Type II supernovae. Paper I was concerned with the application of this model to the observed abundances of several reactive elements and so-called non-metals for three classical intermediate-metallicity clusters, with the hallmark of the work being the successful recovery of many of their well-known elemental and isotopic abundance anomalies. Here, we expand upon our initial analysis by (a) applying the model to a much broader range of metallicities (from the factor of three explored in Paper I, to now, a factor of �50; i.e., essentially, the full range of Galactic globular cluster abundances, and (b) incorporating a broader suite of chemical species, including a number of ironpeak isotopes, heavier �-elements, and fluorine. While allowing for an appropriate fine tuning of the model input parameters, most empirical globular cluster abundance trends are reproduced, our model would suggest the need for a higher production of calcium, silicon, and copper in low-metallicity (or so-called “prompt”) Type Ia supernovae than predicted in current stellar models in order to reproduce the observed trends in NGC 6752, and a factor of two reduction in carbon production from asymptotic giant branch stars to explain the observed trends between carbon and nitrogen. Observations of heavy-element isotopes produced primarily by Type Ia supernovae, including those of titanium, iron, and nickel, could support/refute unequivocally our proposed framework, although currently the feasibility of the proposed observations is well beyond current instrumental capabilities. Hydrodynamical simulations would be necessary to study its viability from a dynamical point of view.

The Lowest Metallicity Stars in the LMC: Clues from MaGICC Simulations

Using a cosmological hydrodynamical simulation of a galaxy of similar mass to the Large Magellanic Cloud (LMC), we examine the predicted characteristics of its lowest metallicity populations. In particular, we emphasise the spatial distributions of first (Pop III) and second (polluted by only immediate Pop III ancestors) generation stars. We find that primordial composition stars form not only in the central galaxys progenitor, but also in locally collapsed subhaloes during the early phases of galaxy formation. The lowest metallicity stars in these subhaloes end up in a relatively extended distribution around the host, with these accreted stars possessing present-day galactocentric distances as great as similar to 40 kpc. By contrast, the earliest stars formed within the central galaxy remain in the inner region, where the vast majority of star formation occurs, for the entirety of the simulation. Consequently, the fraction of stars that are from the earliest generation increases strongly with radius.

The Chemical and Dynamical Evolution of Isolated Dwarf Galaxies

Using a suite of simulations [2] which successfully produce bulgeless (dwarf) disk galaxies, we provide an analysis of their associated cold interstellar media (ISM) and stellar chemical abundance patterns. A preliminary comparison with observations is undertaken, in order to assess whether the properties of the cold gas and chemistry of the stellar components are recovered successfully. To this end, we have extracted the radial and vertical gas density profiles, neutral hydrogen velocity dispersion, and the power spectrum of structure within the ISM. We complement this analysis of the cold gas with a brief examination of the simulations’ metallicity distribution functions and the distribution of α-elements-to-iron.