Showing 1–50 of 99 results for author: Brook, C

A high fidelity Milky Way simulation with Kraken, Gaia-Enceladus, and Sequoia analogues: clues to their accretion histories

Authors: Guacimara García-Bethencourt, Chris B. Brook, Robert J. J. Grand, Daisuke Kawata

Abstract: Within a simulated Milky Way-like galaxy, we identify and analyse analogues of the Gaia-Enceladus (GE), Kraken and Sequoia mergers that each matches remarkably well observational results, including in velocity and chemical abundance space, and their distributions in the $j_{z}$-Energy plane. The Kraken analogue is the earliest merger and has the highest total mass ratio. Consistent with previous s… ▽ More Within a simulated Milky Way-like galaxy, we identify and analyse analogues of the Gaia-Enceladus (GE), Kraken and Sequoia mergers that each matches remarkably well observational results, including in velocity and chemical abundance space, and their distributions in the $j_{z}$-Energy plane. The Kraken analogue is the earliest merger and has the highest total mass ratio. Consistent with previous studies, it is chemically indistinguishable from old in-situ stars at the time of its accretion. The GE and Sequoia analogue events accrete at similar times in our simulation, both along filaments but from opposite sides of the main galaxy. The mean stellar ages of the GE and Sequoia analogues are both similar and, from our simulation results, we see that they can be separate entities and still naturally reproduce the observed properties of their stellar remnants at the present day, including the significant retrograde velocities of the Sequoia analogue remnant stars and the difference in the tracks of the two galaxies through chemical abundance space. Our results provide supporting information about the properties of these three merger events, and show for the first time that they can all be reproduced with a fully cosmological simulation, providing a possible self consistent evolutionary pathway for the Milky Way's formation. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 8 pages, 6 figures, published in MNRAS (accepted version)

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 526, Issue 1, November 2023, Pages 1190-1197

Constraining gas metal mixing strength in simulations using observations of the Milky Way's disc

Authors: Jorge Sarrato-Alós, Christopher Brook, Arianna Di Cintio

Abstract: This work explores the mixing rate of metals in the interstellar medium (ISM), comparing observational constraints from our solar neighbourhood to high resolution cosmological hydrodynamical simulations of Milky Way (MW)-like galaxies. The mixing rate, described by the coefficient C, is varied in simulations between 0 and 0.05, with resultant simulated galaxies compared to observations of metallic… ▽ More This work explores the mixing rate of metals in the interstellar medium (ISM), comparing observational constraints from our solar neighbourhood to high resolution cosmological hydrodynamical simulations of Milky Way (MW)-like galaxies. The mixing rate, described by the coefficient C, is varied in simulations between 0 and 0.05, with resultant simulated galaxies compared to observations of metallicity dispersion in young star clusters, HII regions and neutral gas in the disc of the MW. A value of C between 0.003125 and 0.0125 is found to self-consistently match a range of observables, with a best estimate of C=0.0064$\pm$0.0004. We demonstrate that the relationship between metal dispersion in young stars, HII regions and neutral gas, versus the coefficient C, can be described by a power law. These constrained mixing rates infer a comparatively well mixed ISM in the solar neighbourhood, at odds with some recent observations that have reported a highly inhomogeneous ISM. The degree of mixing suggested by this work is lower than what often employed in many hydrodynamical simulations. Our results have implications for studying the metallicity distribution of stars as well as of gas in the interstellar and circumgalactic media. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 7 pages, 3 figures, accepted for publication in MNRAS

A probabilistic deep learning model to distinguish cusps and cores in dwarf galaxies

Authors: J. Expósito-Márquez, C. B. Brook, M. Huertas-Company, A. Di Cintio, A. V. Macciò, R. J. J. Grand, G. Battaglia

Abstract: Numerical simulations within a cold dark matter (DM) cosmology form halos whose density profiles have a steep inner slope (`cusp'), yet observations of galaxies often point towards a flat central `core'. We develop a convolutional mixture density neural network model to derive a probability density function (PDF) of the inner density slopes of DM halos. We train the network on simulated dwarf gala… ▽ More Numerical simulations within a cold dark matter (DM) cosmology form halos whose density profiles have a steep inner slope (`cusp'), yet observations of galaxies often point towards a flat central `core'. We develop a convolutional mixture density neural network model to derive a probability density function (PDF) of the inner density slopes of DM halos. We train the network on simulated dwarf galaxies from the NIHAO and AURIGA projects, which include both DM cusps and cores: line-of-sight velocities and 2D spatial distributions of their stars are used as inputs to obtain a PDF representing the probability of predicting a specific inner slope. The model recovers accurately the expected DM profiles: $\sim$82$\%$ of the galaxies have a derived inner slope within $\pm$0.1 of their true value, while $\sim$98$\%$ within $\pm$0.3. We apply our model to four Local Group dwarf spheroidal galaxies and find results consistent with those obtained with the Jeans modelling based code GravSphere: the Fornax dSph has a strong indication of possessing a central DM core, Carina and Sextans have cusps (although the latter with large uncertainties), while Sculptor shows a double peaked PDF indicating that a cusp is preferred, but a core can not be ruled out. Our results show that simulation-based inference with neural networks provide a innovative and complementary method for the determination of the inner matter density profiles in galaxies, which in turn can help constrain the properties of the elusive DM. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 14 pages, 11 figures, submitted to MNRAS. comments welcome

NIHAO-LG: The uniqueness of Local Group dwarf galaxies

Authors: Nikhil Arora, Andrea V. Macció, Stéphane Courteau, Tobias Buck, Noam I. Libeskind, Jenny G. Sorce, Chris B. Brook, Yehuda Hoffman, Gustavo Yepes, Eduardo Carlesi, Connor Stone

Abstract: Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulat… ▽ More Recent observational and theoretical studies of the Local Group (LG) dwarf galaxies have highlighted their unique star-formation history, stellar metallicity, gas content, and kinematics. We investigate the commonality of these features by comparing constrained LG and field central dwarf halo simulations in the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project. Our simulations, performed with NIHAO-like hydrodynamics which track the evolution of the Milky Way (MW) and M31 along with ~100 dwarfs in the LG, reveal the total gas mass and stellar properties (velocity dispersion, evolution history, etc.) of present-day LG dwarfs to be similar to field systems. However, relative to field galaxies, LG dwarfs have more cold gas in their central parts and more metal-rich gas in the halo stemming from interactions with other dwarfs living in a high-density environment like the LG. Interestingly, the direct impact of massive MW/M31 analogues on the metallicity evolution of LG dwarfs is minimal; LG dwarfs accrete high-metallicity gas mostly from other dwarfs at late times. We have also tested for the impact of metal diffusion on the chemical evolution of LG dwarfs, and found that it does not affect the stellar or gaseous content of LG dwarfs. Our simulations suggest that the stellar components of LG dwarfs offer a unique and unbiased local laboratory for galaxy-formation tests and comparisons, especially against the overall dwarf population in the Universe. △ Less

Submitted 14 April, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: 16 pages, 12 figures, Accepted for publication in MNRAS

A shallow dark matter halo in Ultra Diffuse Galaxy AGC 242019: are UDGs structurally similar to low surface brightness galaxies?

Authors: Chris B. Brook, Arianna Di Cintio, Andrea V. Maccio, Marvin Blank

Abstract: A central question regarding Ultra Diffuse Galaxies (UDGs) is whether they are a separate category to Low Surface Brightness (LSB) galaxies, or just their natural continuation towards low stellar masses. In this letter, we show that the rotation curve of the gas rich UDG AGC 242019 is well fit by a dark matter halo with inner slope that asymptotes to -0.54, and that such fit provides a concentrati… ▽ More A central question regarding Ultra Diffuse Galaxies (UDGs) is whether they are a separate category to Low Surface Brightness (LSB) galaxies, or just their natural continuation towards low stellar masses. In this letter, we show that the rotation curve of the gas rich UDG AGC 242019 is well fit by a dark matter halo with inner slope that asymptotes to -0.54, and that such fit provides a concentration parameter that matches theoretical expectations. This finding, together with previously works in which shallow inner profiles are derived for UDGs, shows that the structural properties of these galaxies are like other observed LSBs. UDGs show slowly rising rotation curves and this favours formation scenarios in which internal processes, such as SNae driven gas outflows, are acting to modify UDGs profiles. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJL

NIHAO XXIV: Rotation or pressure supported systems? Simulated Ultra Diffuse Galaxies show a broad distribution in their stellar kinematics

Authors: Salvador Cardona-Barrero, Arianna Di Cintio, Christopher B. A. Brook, Tomas Ruiz-Lara, Michael A. Beasley, Jesus Falcón-Barroso, Andrea V. Macciò

Abstract: In recent years a new window on galaxy evolution opened, thanks to the increasing discovery of galaxies with a low surface brightness, such as Ultra Diffuse Galaxies (UDGs). The formation mechanism of these systems is still a much debated question, and so are their kinematical properties. In this work, we address this topic by analyzing the stellar kinematics of isolated UDGs formed in the hydrody… ▽ More In recent years a new window on galaxy evolution opened, thanks to the increasing discovery of galaxies with a low surface brightness, such as Ultra Diffuse Galaxies (UDGs). The formation mechanism of these systems is still a much debated question, and so are their kinematical properties. In this work, we address this topic by analyzing the stellar kinematics of isolated UDGs formed in the hydrodynamical simulation suite NIHAO. We construct projected line-of-sight velocity and velocity dispersion maps to compute the projected specific angular momentum, $λ_{\rm R}$, to characterize the kinematical support of the stars in these galaxies. We found that UDGs cover a broad distribution, ranging from dispersion to rotation supported galaxies, with similar abundances in both regimes. The degree of rotation support of simulated UDGs correlates with several properties such as galaxy morphology, higher HI fractions and larger effective radii with respect to the dispersion supported group, while the dark matter halo spin and mass accretion history are similar amongst the two populations. We demonstrate that the alignment of the infalling baryons into the protogalaxy at early $z$ is the principal driver of the $z$=0 stellar kinematic state: pressure supported isolated UDGs form via mis-aligned gas accretion while rotation supported ones build-up their baryons in an ordered manner. Accounting for random inclination effects, we predict that a comprehensive survey will find nearly half of field UDGs to have rotationally supported stellar disks. △ Less

Submitted 16 July, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

The mass of our Galaxy from satellite proper motions in the Gaia era

Authors: T. K. Fritz, A. Di Cintio, G. Battaglia, C. Brook, S. Taibi

Abstract: We use Gaia DR2 systemic proper motions of 45 satellite galaxies to constrain the mass of the Milky Way using the scale free mass estimator of Watkins et al. (2010). We first determine the anisotropy parameter $β$, and the tracer satellites' radial density index $γ$ to be $β$=$-0.67^{+0.45}_{-0.62}$ and $γ=2.11\pm0.23$. When we exclude possible former satellites of the Large Magellanic Cloud, the… ▽ More We use Gaia DR2 systemic proper motions of 45 satellite galaxies to constrain the mass of the Milky Way using the scale free mass estimator of Watkins et al. (2010). We first determine the anisotropy parameter $β$, and the tracer satellites' radial density index $γ$ to be $β$=$-0.67^{+0.45}_{-0.62}$ and $γ=2.11\pm0.23$. When we exclude possible former satellites of the Large Magellanic Cloud, the anisotropy changes to $β$=$-0.21^{+0.37}_{-0.51}$. We find that the index of the Milky Way's gravitational potential $α$, which is dependent on the mass itself, is the parameter with the largest impact on the mass determination. Via comparison with cosmological simulations of Milky Way-like galaxies, we carried out a detailed analysis of the estimation of the observational uncertainties and their impact on the mass estimator. We found that the mass estimator is biased when applied naively to the satellites of simulated Milky Way halos. Correcting for this bias, we obtain for our Galaxy a mass of $0.58^{+0.15}_{-0.14}\times10^{12}$M$_\odot$ within 64 kpc, as computed from the inner half of our observational sample, and $1.43^{+0.35}_{-0.32}\times10^{12}$M$_\odot$ within 273 kpc, from the full sample; this latter value extrapolates to a virial mass of $M_\mathrm{vir\,Δ=97}$=$1.51^{+0.45}_{-0.40} \times 10^{12}M_{\odot}$ corresponding to a virial radius of R$_\mathrm{vir}$=$308\pm29$ kpc. This value of the Milky Way mass lies in-between other mass estimates reported in the literature, from various different methods. △ Less

Submitted 10 April, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 17 pages, 10 figures, 2 tables. Accepted for publication in MNRAS

Explaining the chemical trajectories of accreted and in-situ halo stars of the Milky Way

Authors: Chris B. Brook, Daisuke Kawata, Brad K. Gibson, Carme Gallart, Andrés Vicente

Abstract: The Milky Way underwent its last significant merger ten billion years ago, when the Gaia-Enceladus-Sausage (GES) was accreted. Accreted GES stars and progenitor stars born prior to the merger make up the bulk of the inner halo. Even though these two main populations of halo stars have similar $durations$ of star formation prior to their merger, they differ in [$α$/Fe]-[Fe/H] space, with the GES po… ▽ More The Milky Way underwent its last significant merger ten billion years ago, when the Gaia-Enceladus-Sausage (GES) was accreted. Accreted GES stars and progenitor stars born prior to the merger make up the bulk of the inner halo. Even though these two main populations of halo stars have similar $durations$ of star formation prior to their merger, they differ in [$α$/Fe]-[Fe/H] space, with the GES population bending to lower [$α$/Fe] at a relatively low value of [Fe/H]. We use cosmological simulations of a 'Milky Way' to argue that the different tracks of the halo stars through the [$α$/Fe]-[Fe/H] plane are due to a difference in their star formation history and efficiency, with the lower mass GES having its low and constant star formation regulated by feedback whilst the higher mass main progenitor has a higher star formation rate prior to the merger. The lower star formation efficiency of GES leads to lower gas pollution levels, pushing [$α$/Fe]-[Fe/H] tracks to the left. In addition, the increasing star formation rate maintains a higher relative contribution of Type~II SNe to Type~Ia SNe for the main progenitor population that formed during the same time period, thus maintaining a relatively high [$α$/Fe]. Thus the different positions of the downturns in the [$α$/Fe]-[Fe/H] plane for the GES stars are not reflective of different star formation durations, but instead reflect different star formation efficiencies. We argue that cosmological simulations match a wide range of independent observations, breaking degeneracies that exist in simpler models. △ Less

Submitted 16 April, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: Accepted version MNRAS

NIHAO XXI: The emergence of Low Surface Brightness galaxies

Authors: Arianna Di Cintio, Chris B. Brook, Andrea V. Macciò, Aaron A. Dutton, Salvador Cardona-Barrero

Abstract: The existence of galaxies with a surface brightness $μ$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rmΛCDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within h… ▽ More The existence of galaxies with a surface brightness $μ$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rmΛCDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within hydrodynamical cosmological simulations from the NIHAO suite. The simulated and observed LSBs share similar properties, having large HI reservoir, extended star formation histories and effective radii, low Sérsic index and slowly rising rotation curves. The formation mechanism of these objects is explored: simulated LSBs form as a result of co-planar co-rotating mergers and aligned accretion of gas at early times, while perpendicular mergers and mis-aligned gas accretion result in higher $μ$ galaxies by $z$=0. The larger the merger, the stronger the correlation between merger orbital configuration and final $μ$. While the halo spin parameter is consistently high in simulated LSB galaxies, the impact of halo concentration, feedback-driven gas outflows and merger time only plays a minor-to-no role in determining $μ$. Interestingly, the formation scenario of such `classical' LSBs differs from the one of less massive, M$_{\star}$$\sim$10$^{7-9}$M$_{\odot}$, Ultra-Diffuse Galaxies, the latter resulting from the effects of SNae driven gas outflows: a M$_{\star}$ of $\sim$10$^9$M$_{\odot}$ thus represents the transition regime between a feedback-dominated to an angular momentum-dominated formation scenario in the LSB realm. Observational predictions are offered regarding spatially resolved star formation rates through LSB discs: these, together with upcoming surveys, can be used to verify the proposed emergence scenario of LSB galaxies. △ Less

Submitted 9 April, 2019; v1 submitted 24 January, 2019; originally announced January 2019.

Comments: MNRAS Accepted 2019 April 02; in original form 2019 January 24. 15 pages,17 figures. main results in Fig 12 and 14. Added appendix A

The birth of the Milky Way as uncovered by accurate stellar ages with Gaia

Authors: Carme Gallart, Edouard J. Bernard, Chris B. Brook, Tomas Ruiz-Lara, Santi Cassisi, Vanessa Hill, Matteo Monelli

Abstract: Knowledge of ages for stars formed over a galaxy's lifetime is fundamental to understand its formation and evolution. However, stellar ages are difficult to obtain since they cannot be measured from observations, being comparison with stellar models (Soderblom 2010) required. Alternatively, age distributions can be derived applying the robust technique of colour-magnitude diagram fitting (Gallart… ▽ More Knowledge of ages for stars formed over a galaxy's lifetime is fundamental to understand its formation and evolution. However, stellar ages are difficult to obtain since they cannot be measured from observations, being comparison with stellar models (Soderblom 2010) required. Alternatively, age distributions can be derived applying the robust technique of colour-magnitude diagram fitting (Gallart et al. 2005), till now mainly employed to study nearby galaxies. The new distances to individual Milky Way stars from the Gaia mission (Brown et al. 2018) have allowed us to use this technique to derive ages from a thick disk colour-magnitude diagram, and from the enigmatic, two-sequenced colour-magnitude diagram of the kinematically hot local halo (Babusiaux et al. 2018), which blue-sequence has been linked to a major accretion event (Haywood et al. 2018, Helmi et al. 2018). Because accurate ages were lacking, the time of the merger and its role on our Galaxy's early evolution remained unclear. We show that the stars in both halo sequences share identical age distributions, and are older than the bulk of thick disc stars. The sharp halo age cut 10 Gyr ago can be identified with the accretion of Gaia-Enceladus. Along with state-of-the-art cosmological simulations of galaxy formation (Brook et al. 2012), these robust ages allow us to order the early sequence of events that shaped our Galaxy, identifying the red-sequence as the first stars formed within the Milky Way progenitor which, because of their kinematics, can be described as its long sought in-situ halo. △ Less

Submitted 29 July, 2019; v1 submitted 9 January, 2019; originally announced January 2019.

Comments: 18 pages, 4 figures, accepted by Nature Astronomy, July 2019; online version available at https://rdcu.be/bLifw

CLUES about M33: the reversed radial stellar age gradient in the outskirts of Triangulum galaxy

Authors: Robert Mostoghiu, Arianna Di Cintio, Alexander Knebe, Noam I. Libeskind, Ivan Minchev, Chris Brook

Abstract: HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we inves… ▽ More HST/ACS observations along the major axis of M33 show that the mean age of its stars decreases with increasing distance from the galaxy center. Such a behavior is consistent with an inside-out growth of the disc. However, in the outermost observed field, at r$\simeq$11.6 kpc, a reversal of this gradient is detected, with old stars found in high percentages beyond this radius. In this work we investigate the origin of such a reversal in stellar age gradient, by using a simulated M33 analogue from the Constrained Local UniversE Simulations (CLUES). The simulated M33 is similar to the observed one in terms of mass, rotation velocity, surface brightness and, similar to what has been reported in observations, shows a stellar age turnaround at large radii. We demonstrate that this reversal is mostly a result of stellar accretion from old satellite galaxies and, to a lesser extent, of stellar migration of in-situ stars. The old accreted stars, with formation times tf < 4 Gyrs, are kinematically hot and can be differentiated from the in-situ stars by their high velocity dispersion and the fact that they do not have rotationally-supported orbits. In the future, obtaining kinematic information of the stars in the outskirt of M33 will help to verify this scenario. △ Less

Submitted 9 August, 2018; originally announced August 2018.

Comments: accepted for publication in MNRAS; 13 pages, 13 figures

Journal ref: MNRAS 480, Issue 4 (2018) 4455-4467

Stellar Feedback and the Energy Budget of Late-Type Galaxies: Missing Baryons and Core Creation

Authors: Harley Katz, Harry Desmond, Federico Lelli, Stacy McGaugh, Arianna Di Cintio, Chris Brook, James Schombert

Abstract: In a $Λ$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems.… ▽ More In a $Λ$CDM cosmology, galaxy formation is a globally inefficient process: it is often the case that far fewer baryons are observed in galaxy disks than expected from the cosmic baryon fraction. The location of these "missing baryons" is unclear. By fitting halo profiles to the rotation curves of galaxies in the SPARC data set, we measure the "missing baryon" mass for individual late-type systems. Assuming that haloes initially accrete the cosmological baryon fraction, we show that the maximum energy available from supernovae is typically not enough to completely eject these "missing baryons" from a halo, but it is often sufficient to heat them to the virial temperature. The energy available from supernovae has the same scaling with galaxy mass as the energy needed to heat or eject the "missing baryons", indicating that the coupling efficiency of the feedback to the ISM may be constant with galaxy virial mass. We further find that the energy available from supernova feedback is always enough to convert a primordial cusp into a core and has magnitude consistent with what is required to heat the "missing baryons" to the virial temperature. Taking a census of the baryon content of galaxies with ${\rm 10^9<M_{vir}/M_{\odot}<10^{12}}$ reveals that $\sim86\%$ of baryons are likely to be in a hot phase surrounding the galaxies and possibly observable in the X-ray, $\sim7\%$ are in the form of cold gas, and $\sim7\%$ are in stars. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: 16 pages, 8 figures, MNRAS accepted

On the early evolution of Local Group dwarf galaxy types: star formation and supernova feedback

Authors: Jose R. Bermejo-Climent, Giuseppina Battaglia, Carme Gallart, Arianna Di Cintio, Chris B. Brook, Luis Cicuendez, Matteo Monelli, Ryan Leaman, Lucio Mayer, Jorge Penarrubia, Justin I. Read

Abstract: According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are p… ▽ More According to star formation histories (SFHs), Local Group dwarf galaxies can be broadly classified in two types: those forming most of their stars before $z=2$ (${\it fast}$) and those with more extended SFHs (${\it slow}$). The most precise SFHs are usually derived from deep but not very spatially extended photometric data; this might alter the ratio of old to young stars when age gradients are present. Here we correct for this effect and derive the mass formed in stars by $z=2$ for a sample of 16 Local Group dwarf galaxies. We explore early differences between ${\it fast}$ and ${\it slow}$ dwarfs, and evaluate the impact of internal feedback by supernovae (SN) on the baryonic and dark matter (DM) component of the dwarfs. ${\it Fast}$ dwarfs assembled more stellar mass at early times and have larger amounts of DM within the half-light radius than ${\it slow}$ dwarfs. By imposing that ${\it slow}$ dwarfs cannot have lost their gas by $z=2$, we constrain the maximum coupling efficiency of SN feedback to the gas and to the DM to be $\sim$10%. We find that internal feedback alone appears insufficient to quench the SFH of ${\it fast}$ dwarfs by gas deprivation, in particular for the fainter systems. Nonetheless, SN feedback can core the DM halo density profiles relatively easily, producing cores of the sizes of the half-light radius in ${\it fast}$ dwarfs by $z=2$ with very low efficiencies. Amongst the "classical" Milky Way satellites, we predict that the smallest cores should be found in Draco and Ursa Minor, while Sculptor and Fornax should host the largest ones. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 15 pages, 3 figures, 2 tables. Accepted for publication in MNRAS

Gaia DR2 Proper Motions of Dwarf Galaxies within 420 kpc: Orbits, Milky Way Mass, Tidal Influences, Planar Alignments, and Group Infall

Authors: T. K. Fritz, G. Battaglia, M. S. Pawlowski, N. Kallivayalil, R. van der Marel, T. S. Sohn, C. Brook, G. Besla

Abstract: A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectr… ▽ More A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectroscopically observed in the literature. We derive systemic PMs for 39 galaxies and galaxy candidates out to 420 kpc, and generally find good consistency for the subset with measurements available from other studies. We derive the implied Galactocentric velocities, and calculate orbits in canonical MW halo potentials of "low" ($0.8 \times 10^{12} M_\odot$) and "high" mass ($1.6 \times 10^{12} M_\odot$). Comparison of the distributions of orbital apocenters and 3D velocities to the halo virial radius and escape velocity, respectively, suggests that the satellite kinematics are best explained in the high-mass halo. Tuc III, Crater II, and additional candidates have orbital pericenters small enough to imply significant tidal influences. Relevant to the missing satellite problem, the fact that fewer galaxies are observed to be near apocenter than near pericenter implies that there must be a population of distant dwarf galaxies yet to be discovered. Of the 39 dwarfs: 12 have orbital poles that do not align with the MW plane of satellites (given reasonable assumptions about its intrinsic thickness); 10 have insufficient PM accuracy to establish whether they align; and 17 satellites align, of which 11 are co-orbiting and (somewhat surprisingly, in view of prior knowledge) 6 are counter-orbiting. Group infall might have contributed to this, but no definitive association is found for the members of the Crater-Leo group. △ Less

Submitted 29 August, 2018; v1 submitted 2 May, 2018; originally announced May 2018.

Comments: 18 pages, 7 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 619, A103 (2018)

Spectroscopic characterisation of the stellar content of ultra diffuse galaxies

Authors: T. Ruiz-Lara, M. A. Beasley, J. Falcón-Barroso, J. Román, F. Pinna, C. Brook, A. Di Cintio, I. Martín-Navarro, I. Trujillo, A. Vazdekis

Abstract: Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star for… ▽ More Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star formation histories (SFH) and mean population parameters, such as ages and metallicities, for a sample of five UDG candidates in the Coma cluster. From the radial velocities, we confirm the Coma membership of these galaxies. We find that their rotation properties, if detected at all, are compatible with dwarf-like galaxies. The SFHs of the UDG are dominated by old (~ 7 Gyr), metal-poor ([M/H] ~ -1.1) and alpha-enhanced ([Mg/Fe] ~ 0.4) populations followed by a smooth or episodic decline which halted ~ 2 Gyr ago, possibly a sign of cluster-induced quenching. We find no obvious correlation between individual SFH shapes and any UDG morphological properties. The recovered stellar properties for UDGs are similar to those found for DDO44, a local UDG analogue resolved into stars. We conclude that the UDGs in our sample are extended dwarfs whose properties are likely the outcome of both internal processes, such as bursty SFHs and/or high-spin haloes, as well as environmental effects within the Coma cluster. △ Less

Submitted 6 July, 2018; v1 submitted 16 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS

The Radial Distribution of Mono-Metallicity Populations in the Galactic Disk as Evidence for Two-Phase Disk Formation

Authors: R. Domínguez-Tenreiro, A. Obreja, C. B. Brook, F. J. Martínez-Serrano, A. Serna

Abstract: Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [$α$/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin - thick disk dichotomy. The analysis of these observations lead to the non-$[α$/Fe] enhanced populations splitting into MMPs with different surface densities according to thei… ▽ More Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [$α$/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin - thick disk dichotomy. The analysis of these observations lead to the non-$[α$/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, $[α$/Fe] enhanced (i.e., old) populations show an homogeneous behaviour. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the $\rm ΛCDM$ model, we have found that the two phases of halo mass assembly (an early, fast phase, followed by a slow one, with low mass assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing has only a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients, and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars to suffer from strong radial mixing. By linking the two-component disk concept with the two-phase halo mass assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event marking periods characterized by different physical conditions under which thick and thin disk stars were born. △ Less

Submitted 15 June, 2017; originally announced June 2017.

Comments: accepted for publication in ApJ

NIHAO XIV: Reproducing the observed diversity of dwarf galaxy rotation curve shapes in LCDM

Authors: Isabel M. Santos-Santos, Arianna Di Cintio, Chris B. Brook, Andrea Macciò, Aaron Dutton, Rosa Domínguez-Tenreiro

Abstract: The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO gala… ▽ More The significant diversity of rotation curve (RC) shapes in dwarf galaxies has recently emerged as a challenge to LCDM: in dark matter (DM) only simulations, DM halos have a universal cuspy density profile that results in self-similar RC shapes. We compare RC shapes of simulated galaxies from the NIHAO project with observed galaxies from the homogeneous SPARC dataset. The DM halos of the NIHAO galaxies can expand to form cores, with the degree of expansion depending on their stellar-to-halo mass ratio. By means of the $V_{\rm 2kpc}-V_{\rm Rlast}$ relation (where $V_{\rm Rlast}$ is the outermost measured rotation velocity), we show that both the average trend and the scatter in RC shapes of NIHAO galaxies are in reasonable agreement with SPARC: this represents a significant improvement compared to simulations that do not result in DM core formation, suggesting that halo expansion is a key process in matching the diversity of dwarf galaxy RCs. Note that NIHAO galaxies can reproduce even the extremely slowly rising RCs of IC 2574 and UGC 5750. Revealingly, the range where observed galaxies show the highest diversity corresponds to the range where core formation is most efficient in NIHAO simulations, 50$<V_{\rm Rlast}$\km s$^{-1}<$100. A few observed galaxies in this range cannot be matched by any NIHAO RC nor by simulations that predict a universal halo profile. Interestingly, the majority of these are starbursts or emission-line galaxies, with steep RCs and small effective radii. Such galaxies represent an interesting observational target providing new clues to the process/viability of cusp-core transformation, the relationship between starburst and inner potential well, and the nature of DM. △ Less

Submitted 16 October, 2017; v1 submitted 13 June, 2017; originally announced June 2017.

Comments: 13 pages, 6 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society 9 October 2017

Metallicity Gradient of the Thick Disc Progenitor at High Redshift

Authors: D. Kawata, C. Allende Prieto, C. B. Brook, L. Casagrande, I. Ciucă, B. K. Gibson, R. J. J. Grand, M. R. Hayden, J. A. S. Hunt

Abstract: We have developed a novel Markov Chain Mote Carlo (MCMC) chemical "painting" technique to explore possible radial and vertical metallicity gradients for the thick disc progenitor. In our analysis we match an N-body simulation to the data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We assume that the thick disc has a constant scale-height and has completed its f… ▽ More We have developed a novel Markov Chain Mote Carlo (MCMC) chemical "painting" technique to explore possible radial and vertical metallicity gradients for the thick disc progenitor. In our analysis we match an N-body simulation to the data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We assume that the thick disc has a constant scale-height and has completed its formation at an early epoch, after which time radial mixing of its stars has taken place. Under these assumptions, we find that the initial radial metallicity gradient of the thick disc progenitor should not be negative, but either flat or even positive, to explain the current negative vertical metallicity gradient of the thick disc. Our study suggests that the thick disc was built-up in an inside-out and upside-down fashion, and older, smaller and thicker populations are more metal poor. In this case, star forming discs at different epochs of the thick disc formation are allowed to have different radial metallicity gradients, including a negative one, which helps to explain a variety of slopes observed in high redshift disc galaxies. This scenario helps to explain the positive slope of the metallicity-rotation velocity relation observed for the Galactic thick disc. On the other hand, radial mixing flattens the slope of an existing gradient. △ Less

Submitted 20 September, 2017; v1 submitted 5 June, 2017; originally announced June 2017.

Comments: 12 pages, 9 figures, 3 tables, accepted for publication in MNRAS

Diversity of dwarf galaxy IR-submm emission patterns: CLUES from hydrodynamical simulations

Authors: Isabel M. E. Santos-Santos, Rosa Domínguez-Tenreiro, Gian Luigi Granato, Chris B. Brook, Aura Obreja

Abstract: The spectral energy distributions (SEDs) of low-mass low-metallicity (dwarf) galaxies are a challenging piece of the puzzle of galaxy formation in the near Universe. These SEDs show some particular features in the submillimeter to far-infrared wavelength range compared to normal, larger, and metal-richer galaxies that cannot be explained by the current models. These can be summarized as: a broaden… ▽ More The spectral energy distributions (SEDs) of low-mass low-metallicity (dwarf) galaxies are a challenging piece of the puzzle of galaxy formation in the near Universe. These SEDs show some particular features in the submillimeter to far-infrared wavelength range compared to normal, larger, and metal-richer galaxies that cannot be explained by the current models. These can be summarized as: a broadening of the IR peak, which implies a warmer dust component; an excess of emission in the submm ($\sim$500 $μ$m), that causes a flattening of the submm/FIR slope; and a very low intensity of PAH emission features. With the aim of explaining these emission patterns, the SEDs of a sample of 27 simulated dwarf galaxies were calculated using the GRASIL-3D radiation transfer code. This code has the particularity that it separately treats the radiative transfer through dust grains within molecular clouds and within the cirrus, the dense and diffuse components of the gas phase, respectively. The simulated galaxies have stellar masses ranging from 10$^6$-10$^9$ M$_\odot$, and were obtained from a single simulation run within a Local Group environment with initial conditions from the CLUES project. We report a study of the IRAS, Spitzer and Herschel bands luminosities, and of the SFRs, dust, and gas (HI and H$_2$) mass contents. We find a satisfactory agreement with observational data, with GRASIL-3D naturally reproducing the specific spectral features mentioned above. We conclude that the GRASIL-3D two-component dust model gives a physical interpretation of the emission of dwarf galaxies: molecular clouds and cirrus represent the warm and cold dust components, respectively, needed to reproduce observational data. △ Less

Submitted 28 March, 2017; v1 submitted 9 March, 2017; originally announced March 2017.

Comments: 22 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics. Final version for journal

Journal ref: A&A 603, A4 (2017)

NIHAO XI: Formation of Ultra-Diffuse Galaxies by outflows

Authors: Arianna Di Cintio, Chris B. Brook, Aaron A. Dutton, Andrea V. Macciò, Aura C. Obreja, Avishai Dekel

Abstract: 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 h… ▽ More 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érsic 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}$. △ Less

Submitted 22 November, 2016; v1 submitted 3 August, 2016; originally announced August 2016.

Comments: matches accepted version, MNRAS Letter 2016-10-12

Testing Feedback-Modified Dark Matter Haloes with Galaxy Rotation Curves: Estimation of Halo Parameters and Consistency with $Λ$CDM

Authors: Harley Katz, Federico Lelli, Stacy S. McGaugh, Arianna Di Cintio, Chris B. Brook, James M. Schombert

Abstract: Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to th… ▽ More Cosmological $N$-body simulations predict dark matter (DM) haloes with steep central cusps (e.g. NFW, Navarro et al. 1996). This contradicts observations of gas kinematics in low-mass galaxies that imply the existence of shallow DM cores. Baryonic processes such as adiabatic contraction and gas outflows can, in principle, alter the initial DM density profile, yet their relative contributions to the halo transformation remain uncertain. Recent high resolution, cosmological hydrodynamic simulations (Di Cintio et al. 2014, DC14) predict that inner density profiles depend systematically on the ratio of stellar to DM mass (M$_*$/M$_{\text{halo}}$). Using a Markov Chain Monte Carlo approach, we test the NFW and the M$_*$/M$_{\text{halo}}$-dependent DC14 halo models against a sample of 147 galaxy rotation curves from the new {\it Spitzer} Photometry and Accurate Rotation Curves (SPARC) data set. These galaxies all have extended H{\small I} rotation curves from radio interferometry as well as accurate stellar mass density profiles from near-infrared photometry. The DC14 halo profile provides markedly better fits to the data compared to the NFW profile. Unlike NFW, the DC14 halo parameters found in our rotation curve fits naturally fall within two standard deviations of the mass-concentration relation predicted by $Λ$CDM and the stellar mass-halo mass relation inferred from abundance matching with few outliers. Halo profiles modified by baryonic processes are therefore more consistent with expectations from $Λ$ cold dark matter ($Λ$CDM) cosmology and provide better fits to galaxy rotation curves across a wide range of galaxy properties than do halo models that neglect baryonic physics. Our results offer a solution to the decade long cusp-core discrepancy. △ Less

Submitted 11 January, 2017; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: 23 Pages, 18 Figures, MNRAS Accepted

NIHAO IX: the role of gas inflows and outflows in driving the contraction and expansion of cold dark matter haloes

Authors: Aaron A. Dutton, Andrea V. Macciò, Avishai Dekel, Liang Wang, Gregory S. Stinson, Aura Obreja, Arianna Di Cintio, Chris B. Brook, Tobias Buck, Xi Kang

Abstract: We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with… ▽ More We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(Ω_b/Ω_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes. △ Less

Submitted 5 July, 2016; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: 20 pages, 15 figures, accepted to MNRAS

Impacts of a Flaring Star-forming Disc and Stellar Radial Mixing on the Vertical Metallicity Gradient

Authors: D. Kawata, R. J. J. Grand, B. K. Gibson, L. Casagrande, J. A. S. Hunt, C. B. Brook

Abstract: Using idealised N-body simulations of a Milky Way-sized disc galaxy, we qualitatively study how the metallicity distributions of the thin disc star particles are modified by the formation of the bar and spiral arm structures. The thin disc in our numerical experiments initially has a tight negative radial metallicity gradient and a constant vertical scale-height. We show that the radial mixing of… ▽ More Using idealised N-body simulations of a Milky Way-sized disc galaxy, we qualitatively study how the metallicity distributions of the thin disc star particles are modified by the formation of the bar and spiral arm structures. The thin disc in our numerical experiments initially has a tight negative radial metallicity gradient and a constant vertical scale-height. We show that the radial mixing of stars drives a positive vertical metallicity gradient in the thin disc. On the other hand, if the initial thin disc is flared, with vertical scale-height increasing with galactocentric radius, the metal poor stars originally in the outer disc become dominant in regions above the disc plane at every radii. This process can drive a negative vertical metallicity gradient, which is consistent with the current observed trend. This model mimics a scenario where the star-forming thin disc was flared in the outer region at earlier epochs. Our numerical experiment with an initial flared disc predicts that the negative vertical metallicity gradient of the mono-age relatively young thin disc population should be steeper in the inner disc, and the radial metallicity gradient of the mono-age population should be shallower at greater heights above the disc plane. We also predict that the metallicity distribution function of mono-age young thin disc populations above the disc plane would be more positively skewed in the inner disc compared to the outer disc. △ Less

Submitted 14 September, 2016; v1 submitted 25 April, 2016; originally announced April 2016.

Comments: 12 pages, 12 figures, 2 tables, accepted for publication in MNRAS

The different baryonic Tully-Fisher relations at low masses

Authors: C. B. Brook, I. Santos-Santos, G. Stinson

Abstract: We compare the Baryonic Tully-Fisher relation (BTFR) of simulations and observations of galaxies ranging from dwarfs to spirals, using various measures of rotational velocity Vrot. We explore the BTFR when measuring Vrot at the flat part of the rotation curve, Vflat, at the extent of HI gas, Vlast, and using 20% (W20) and 50% (W50) of the width of HI line profiles. We also compare with the maximum… ▽ More We compare the Baryonic Tully-Fisher relation (BTFR) of simulations and observations of galaxies ranging from dwarfs to spirals, using various measures of rotational velocity Vrot. We explore the BTFR when measuring Vrot at the flat part of the rotation curve, Vflat, at the extent of HI gas, Vlast, and using 20% (W20) and 50% (W50) of the width of HI line profiles. We also compare with the maximum circular velocity of the parent halo, Vmax, within dark matter only simulations. The different BTFRs increasingly diverge as galaxy mass decreases. Using Vlast one obtains a power law over four orders of magnitude in baryonic mass, with slope similar to the observed BTFR. Measuring Vflat gives similar results as Vlast when galaxies with rising rotation curves are excluded. However, higher rotation velocities would be found for low mass galaxies if the cold gas extended far enough for Vrot to reach a maximum. W20 gives a similar slope as Vlast but with slightly lower values of Vrot for low mass galaxies, although this may depend on the extent of the gas in your galaxy sample. W50 bends away from these other relations toward low velocities at low masses. By contrast, Vmax bends toward high velocities for low mass galaxies, as cold gas does not extend out to the radius at which halos reach Vmax. Our study highlights the need for careful comparisons between observations and models: one needs to be consistent about the particular method of measuring Vrot, and precise about the radius at which velocities are measured. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Comments: Accepted for publication MNRAS

Origin of the Metallicity Distribution in the Thick Disc

Authors: Maider S. Miranda, Kate Pilkington, Brad K. Gibson, Christopher B. Brook, Patricia Sánchez-Blázquez, Ivan Minchev, Christopher Gareth Few, Rory Smith, Rosa Domínguez-Tenreiro, Aura Obreja, Jeremy Bailin, Greg S. Stinson

Abstract: Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulat… ▽ More Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Comments: 15 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 587, A10 (2016)

The distribution of mass components in simulated disc galaxies

Authors: Isabel M. Santos-Santos, Chris B. Brook, Greg Stinson, Arianna Di Cintio, James Wadsley, Rosa Domínguez-Tenreiro, Stefan Gottlöber, Gustavo Yepes

Abstract: Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V… ▽ More Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter. △ Less

Submitted 13 October, 2015; v1 submitted 8 October, 2015; originally announced October 2015.

Comments: 8 pages, 7 pdf figures. Accepted for publication in MNRAS on 2015 October 06

NIHAO IV: Core creation and destruction in dark matter density profiles across cosmic time

Authors: Edouard Tollet, Andrea V. Macciò, Aaron A. Dutton, Greg S. Stinson, Liang Wang, Camilla Penzo, Thales A. Gutcke, Tobias Buck, Xi Kang, Chris Brook, Arianna Di Cintio, Ben W. Keller, James Wadsley

Abstract: We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter densi… ▽ More We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $α$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $α$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $α$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation. △ Less

Submitted 6 May, 2019; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: 11 pages, 13 figures. Corrected typo in table 2 (middle row) with respect to the version published in MNRAS

The Spanish Square Kilometre Array White Book

Authors: J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo , et al. (94 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on… ▽ More The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project. △ Less

Submitted 17 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Española de Astronomía

Signatures of Dark Matter Halo Expansion in Galaxy Populations

Authors: Chris B. Brook, Arianna Di Cintio

Abstract: Dark matter cores within galaxy haloes can be formed by energy feedback from star forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with M$_{\star}\sim10^{8.5}$M$_{\odot}$. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity func… ▽ More Dark matter cores within galaxy haloes can be formed by energy feedback from star forming regions: an energy balance suggests that the maximum core formation efficiency arises in galaxies with M$_{\star}\sim10^{8.5}$M$_{\odot}$. We show that a model population of galaxies, in which the density profile has been modified by such baryonic feedback, is able to explain the observed galaxy velocity function and Tully-Fisher relations significantly better than a model in which a universal cuspy density profile is assumed. Alternative models, namely warm or self-interacting dark matter, also provide a better match to these observed relations than a universal profile model does, but make different predictions for how halo density profiles vary with mass compared to the baryonic feedback case. We propose that different core formation mechanisms may be distinguished based on the imprint they leave on galaxy populations over a wide range of mass. Within the current observational data we find evidence of the expected signatures of the mass dependence of core formation generated by baryonic feedback. △ Less

Submitted 30 July, 2015; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 11 pages, 9 figures, accepted in MNRAS

A matter of measurement: rotation velocities and the velocity function of dwarf galaxies

Authors: Chris Brook, Francesco Shankar

Abstract: The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vro… ▽ More The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax. △ Less

Submitted 1 June, 2015; originally announced June 2015.

Comments: submitted to MNRAS. Comments very welcome

The Variation of Rotation Curve Shapes as a Signature of the Effects of Baryons on Dark Matter Density Profiles

Authors: Chris Brook

Abstract: Rotation curves of galaxies show a wide range of shapes, which can be paramaterized as scatter in Vrot(1kpc)/Vmax i.e.the ratio of the rotation velocity measured at 1kpc and the maximum measured rotation velocity. We examine whether the observed scatter can be accounted for by combining scatters in disc scale-lengths, the concentration-halo mass relation, and the M*-Mhalo relation. We use these sc… ▽ More Rotation curves of galaxies show a wide range of shapes, which can be paramaterized as scatter in Vrot(1kpc)/Vmax i.e.the ratio of the rotation velocity measured at 1kpc and the maximum measured rotation velocity. We examine whether the observed scatter can be accounted for by combining scatters in disc scale-lengths, the concentration-halo mass relation, and the M*-Mhalo relation. We use these scatters to create model galaxy populations; when housed within dark matter halos that have universal, NFW density profiles, the model does not match the lowest observed values of Vrot(1kpc)/Vmax and has too little scatter in Vrot(1kpc)/Vmax compared to observations. By contrast, a model using a mass dependent dark matter profile, where the inner slope is determined by the ratio of M*/Mhalo, produces galaxies with low values of Vrot(1kpc)/Vmax and a much larger scatter, both in agreement with observation. We conclude that the large observed scatter in Vrot(1kpc)/Vmax favours density profiles that are significantly affected by baryonic processes. Alternative dark matter core formation models such as SIDM may also account for the observed variation in rotation curve shapes, but these observations may provide important constraints in terms of core sizes, and whether they vary with halo mass and/or merger history. △ Less

Submitted 16 October, 2015; v1 submitted 31 May, 2015; originally announced June 2015.

Comments: MNRAS final version with updated treatment of alternative DM models

Expanded haloes, abundance matching and too-big-to-fail in the Local Group

Authors: Chris B. Brook, Arianna Di Cintio

Abstract: Observed kinematical data of 40 Local Group (LG) members are used to derive the dark matter halo mass of such galaxies. Haloes are selected from the theoretically expected LG mass function and two different density profiles are assumed, a standard universal cuspy model and a mass dependent profile which accounts for the effects of baryons in modifying the dark matter distribution within galaxies.… ▽ More Observed kinematical data of 40 Local Group (LG) members are used to derive the dark matter halo mass of such galaxies. Haloes are selected from the theoretically expected LG mass function and two different density profiles are assumed, a standard universal cuspy model and a mass dependent profile which accounts for the effects of baryons in modifying the dark matter distribution within galaxies. The resulting relations between stellar and halo mass are compared with expectations from abundance matching. Using a universal cuspy profile, the ensemble of LG galaxies is fit in relatively low mass haloes, leaving "dark" many massive haloes of \mhalo$\gtrsim$10$^{10}$\msun: this reflects the "too big to fail" problem and results in a \mstar-\mhalo\ relation that differs from abundance matching predictions. Moreover, the star formation efficiency of isolated LG galaxies increases with decreasing halo mass when adopting a cuspy model. By contrast, using the mass dependent density profile, dwarf galaxies with \mstar$\gtrsim$10$^{6}$\msun are assigned to more massive haloes, which have a central cored distribution of dark matter: the "too big to fail" problem is alleviated, the resultant \mstar-\mhalo\ relation follows abundance matching predictions down to the completeness limit of current surveys, and the star formation efficiency of isolated members decreases with decreasing halo mass, in agreement with theoretical expectations. Finally, the cusp/core space of LG galaxies is presented, providing a framework to understand the non-universality of their density profiles. △ Less

Submitted 15 April, 2015; v1 submitted 14 October, 2014; originally announced October 2014.

Comments: Accepted in MNRAS, 15 pages 7 figures. section 3.3 has been added after reports

A mass-dependent density profile for dark matter haloes including the influence of galaxy formation

Authors: Arianna Di Cintio, Chris B. Brook, Aaron A. Dutton, Andrea V. Macciò, Greg S. Stinson, Alexander Knebe

Abstract: We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range… ▽ More We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($γ$) and outer ($β$) slopes of dark matter density, and the sharpness of transition between the slopes($α$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile. △ Less

Submitted 23 April, 2014; originally announced April 2014.

Comments: 10 pages, 4 figures. Accepted for publication in MNRAS

The Main Sequence and the Fundamental Metallicity Relation in MaGICC Galaxies: Evolution and Scatter

Authors: A. Obreja, C. B. Brook, G. Stinson, R. Domínguez-Tenreiro, B. K. Gibson, L. Silva, G. L. Granato

Abstract: Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations… ▽ More Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence. △ Less

Submitted 31 March, 2014; originally announced April 2014.

Comments: 12 pages, 8 figures, submitted to MNRAS

The Role of Feedback in Shaping the Structure of the Interstellar Medium

Authors: A. P. Walker, B. K. Gibson, K. Pilkington, C. B. Brook, P. Dutta, S. Stanimirovic, G. S. Stinson, J. Bailin

Abstract: We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density)… ▽ More We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, 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 THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -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 = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI 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. △ Less

Submitted 3 March, 2014; originally announced March 2014.

Comments: Accepted for publication in MNRAS; 9 pages; 7 figures

GRASIL-3D: an Implemention of Dust Effects in the SEDs of Simulated Galaxies

Authors: R. Domínguez-Tenreiro, A. Obreja, G. L. Granato, A. Schurer, P. Alpresa, L. Silva, C. B. Brook, A. Serna

Abstract: We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as g… ▽ More We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands. △ Less

Submitted 21 February, 2014; v1 submitted 3 December, 2013; originally announced December 2013.

Comments: 24 pages, 13 Figures and 6 Tables. Accepted in MNRAS

The stellar-to-halo mass relation for Local Group galaxies

Authors: C. B. Brook, A. Di Cintio, A. Knebe, S. Gottlöber, Y. Hoffman, G. Yepes, S. Garrison-Kimmel

Abstract: We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local u… ▽ More We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local universe simulations (CLUES). For the stellar mass range 10$^7$<M*<10$^8$M$_\odot$, for which we likely have a complete census of observed galaxies, we find that the stellar mass-halo mass relation follows a power law with slope of 3.1, significantly steeper than most values in the literature. The steep relation between stellar and halo masses indicates that Local Group dwarf galaxies are hosted by dark matter halos with a small range of mass. Our methodology is robust down to the stellar mass to which the census of observed Local Group galaxies is complete, but the significant uncertainty in the currently measured slope of the stellar-to halo mass relation will decrease dramatically if the Local Group completeness limit was $10^{6.5}$M$\odot$ or below, highlighting the importance of pushing such limit to lower masses and larger volumes. △ Less

Submitted 1 February, 2014; v1 submitted 21 November, 2013; originally announced November 2013.

Comments: ApJL Accepted version. Comments to cbabrook@gmail.com

MaGICC-WDM: the effects of warm dark matter in hydrodynamical simulations of disc galaxy formation

Authors: Jakob Herpich, Gregory S. Stinson, Andrea V. Macciò, Chris Brook, James Wadsley, Hugh M. P. Couchman, Tom Quinn

Abstract: We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar prof… ▽ More We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation. △ Less

Submitted 27 November, 2013; v1 submitted 5 August, 2013; originally announced August 2013.

Comments: 13 pages, 9 figures, 2 tables; minor clarifications added in results section, conclusions unchanged; accepted for publication in MNRAS

Journal ref: 2014 MNRAS, 437, 293

The MaGICC Baryon Cycle: The Enrichment History of Simulated Disc Galaxies

Authors: C. B. Brook, G. Stinson, B. K. Gibson, S. Shen, A. V. Macciò, J. Wadsley, T. Quinn

Abstract: Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We de… ▽ More Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We derive the distributions of time between gas crossing the virial radius and being accreted to the star forming region (which allows mixing within the corona), as well as the time between gas being accreted to the star forming region and then forming stars (which allows 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.1Z 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 simulation (Mvir~2x10^{10}Msun, with M*~10^7Msun), by contrast, accretes primordial gas through the virial radius and onto the disc at all times. Much like the classical analytical solutions to the `G-dwarf problem', overproduction of low-metallicity stars is ameliorated by the inefficiency of star formation. Finally, gas outflow/metal removal rates from star forming regions as a function of galactic mass are presented. △ Less

Submitted 25 January, 2016; v1 submitted 24 June, 2013; originally announced June 2013.

Comments: published in MNRAS

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

Authors: Chris B. Brook, Maider S. Miranda, Brad. K. Gibson, Kate Pilkington, Greg S. Stinson

Abstract: 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… ▽ More 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 galaxy's progenitor, but also in locally collapsed sub-halos during the early phases of galaxy formation. The lowest metallicity stars in these sub-halos end up in a relatively extended distribution around the host, with these accreted stars possessing present-day galactocentric distances as great as ~40kpc. 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. △ Less

Submitted 8 July, 2013; v1 submitted 15 June, 2013; originally announced June 2013.

Comments: accepted for publication in PASA, as part of the "The Magellanic System: In Perspective" edition

The dependence of dark matter profiles on the stellar to halo mass ratio: a prediction for cusps vs cores

Authors: Arianna Di Cintio, Chris B. Brook, Andrea V. Maccio, Greg S. Stinson, Alexander Knebe, Aaron A. Dutton, James Wadsley

Abstract: We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks… ▽ More We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, α in ρr^α, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of α on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how α varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with α decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes α<-0.8, approaching the NFW profile. This novel prediction for the dependence of α on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies. △ Less

Submitted 7 October, 2013; v1 submitted 4 June, 2013; originally announced June 2013.

Comments: accepted version, MNRAS in press. results unchanged,expanded discussion in sec 3.1,3.2,conclusions

Constraining Sub-Grid Physics with High-Redshift Spatially-Resolved Metallicity Distributions

Authors: B. K. Gibson, K. Pilkington, C. B. Brook, G. S. Stinson, J. Bailin

Abstract: We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We mad… ▽ More We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We made use of a suite of L* discs, realised with and without `enhanced' feedback. All the simulations were produced using the smoothed particle hydrodynamics code Gasoline, and their in situ gas-phase metallicity gradients traced from redshift z~2 to the present-day. Present-day age-metallicity relations and metallicity distribution functions were derived for each system. The `enhanced' feedback models, which have been shown to be in agreement with a broad range of empirical scaling relations, distribute energy and re-cycled ISM material over large scales and predict the existence of relatively `flat' and temporally invariant abundance gradients. Enhanced feedback schemes reduce significantly the scatter in the local stellar age-metallicity relation and, especially, the [O/Fe]-[Fe/H] relation. The local [O/Fe] distribution functions for our L* discs show clear bimodality, with peaks at [O/Fe]=-0.05 and +0.05 (for stars with [Fe/H]>-1), consistent with our earlier work on dwarf discs. Our results with `enhanced' feedback are inconsistent with our earlier generation of simulations realised with `conservative' feedback. We conclude that spatially-resolved metallicity distributions, particularly at high-redshift, offer a unique and under-utilised constraint on the uncertain nature of stellar feedback processes. △ Less

Submitted 14 May, 2013; v1 submitted 10 April, 2013; originally announced April 2013.

Comments: Accepted for publication in A&A; 5 pages; 3 figures; inset/legend to figure fixed (misaligned symbol fixed!)

The MaGICC volume: reproducing statistical properties of high redshift galaxies

Authors: Rahul Kannan, Greg S. Stinson, Andrea V. Macciò, Chris Brook, Simone M. Weinmann, James Wadsley, H. M. P. Couchman

Abstract: We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compar… ▽ More We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high redshift galaxies ($z \ge 2$) including the stellar mass - halo mass ($M_\star - M_h$ ) relation, the Galaxy Stellar Mass Function (GSMF) at low masses ($M_\star \lt 5 \times 10^{10} M_\odot$ ) and the number density evolution of low mass galaxies. The poor match of $M_\star - M_h$ and the GSMF at high masses ($M_\star \ge 5 \times 10^{10} M_\odot$ ) indicates supernova feedback is insufficient to limit star formation in these haloes. At $z = 0$, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around $L_\star$ mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernovae feedback, plays a major role in regulating the properties of low mass galaxies at high redshift. △ Less

Submitted 5 November, 2013; v1 submitted 11 February, 2013; originally announced February 2013.

Comments: Simulation extended to z=0, conclusions unchanged. Accepted for publication in MNRAS

MaGICC Thick Disk I: Comparing a Simulated Disk Formed with Stellar Feedback to the Milky Way

Authors: G. S. Stinson, J. Bovy, H. -W. Rix, C. Brook, R. Roškar, J. J. Dalcanton, A. V. Macciò, J. Wadsley, H. M. P. Couchman, T. R. Quinn

Abstract: We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [α/Fe] abundance ranges, we find remarkable… ▽ More We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [α/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances. △ Less

Submitted 1 October, 2013; v1 submitted 22 January, 2013; originally announced January 2013.

Comments: 11 page, 10 figures, submitted. Comments welcome

A two-phase scenario for bulge assembly in LCDM cosmologies

Authors: A. Obreja, R. Domínguez-Tenreiro, C. Brook, F. J. Martínez-Serrano, M. Doménech-Moral, A. Serna, M. Mollá, G. Stinson

Abstract: We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback… ▽ More We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback which drives large scale outflows. In all cases, the marked knee-shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z=0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L* galaxies will generically have two bulge populations which can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories. △ Less

Submitted 16 November, 2012; originally announced November 2012.

Comments: 24 pages, 10 figures, accepted for publication in ApJ

Properties of simulated Milky Way-mass galaxies in loose group and field environments

Authors: C. G. Few, B. K. Gibson, S. Courty, L. Michel-Dansac, C. B. Brook, G. S. Stinson

Abstract: We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart fr… ▽ More We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart from the environments of the galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. Comparison of these two samples allow for systematic differences in the simulations to be identified. Metallicity gradients, disk scale lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link between these and environment of the galaxies is studied. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their velocity dispersion-age relations. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar assembly history is considered. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Comments: 16 pages, 11 figures, abstract abridged for arXiv. Accepted for publication in Astronomy & Astrophysics

Making Galaxies in a Cosmological Context: The Need for Early Stellar Feedback

Authors: G. Stinson, C. Brook, A. V. Macciò, J. Wadsley, T. R. Quinn, H. M. P. Couchman

Abstract: We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match th… ▽ More We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass--halo mass relationship, the galaxy forms too many stars before z=2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like a UV ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >10^6 K before cooling to 10^4 K. The pressure from this hot gas creates a more extended disk and prevents more star formation prior to z=1 than supernovae feedback alone. The resulting disk galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disk scaling relationships. The disk forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts. △ Less

Submitted 15 October, 2012; v1 submitted 31 July, 2012; originally announced August 2012.

Comments: 13 pages, 14 figures, accepted MNRAS, movies at http://www.mpia.de/~stinson/magicc

The Halo Shape and Evolution of Polar Disc Galaxies

Authors: O. N. Snaith, B. K. Gibson, C. B. Brook, A. Knebe, R. J. Thacker, T. R. Quinn, F. Governato, P. B. Tissera

Abstract: We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid f… ▽ More We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger. △ Less

Submitted 25 June, 2012; originally announced June 2012.

Comments: Accepted for publication in MNRAS; 14 pages; 14 figures

Thin disc, Thick Disc and Halo in a Simulated Galaxy

Authors: C. B. Brook, G. S. Stinson, B. K. Gibson, D. Kawata, E. L. House, M. S. Miranda, A. V. Macciò, K. Pilkington, R. Roškar, J. Wadsley, T. R. Quinn

Abstract: Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion (σU,V,W = 41, 3… ▽ More Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion (σU,V,W = 41, 31, 25 km/s), high [Fe/H], and low [O/Fe]. The thick disc stars are old (6.6<age<9.8 Gyrs), lag the thin disc by \sim21 km/s, possess higher velocity dispersion (σU,V,W = 49, 44, 35 km/s), relatively low [Fe/H] and high [O/Fe]. The halo component comprises less than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a velocity ellipsoid defined by (σU,V,W = 62, 46, 45 km/s) and is formed primarily in-situ during an early merger epoch. Gas-rich mergers during this epoch play a major role in fuelling the formation of the old disc stars (the thick disc). This is consistent with studies which show that cold accretion is the main source of a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength \sim1.7 kpc at z=1) and kinematically hot disc, primarily from gas accreted during the galaxy's merger epoch. Far from being a competing formation scenario, migration is crucial for reconciling the short, hot, discs which form at high redshift in ΛCDM, with the properties of the thick disc at z=0. The thick disc, as defined by its abundances maintains its relatively short scale-length at z = 0 (2.31 kpc) compared with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is imprinted the evolution of abundances such that the metal poor α-young population has a larger scale-length (4.07 kpc) than the more chemically evolved metal rich α-young population (2.74 kpc). △ Less

Submitted 4 June, 2012; originally announced June 2012.

Comments: Submitted to MNRAS. This version after helpful referee comments. Comments welcome to cbabrook@gmail.com

The Distribution of Metals in Cosmological Hydrodynamical Simulations of Dwarf Disk Galaxies

Authors: K. Pilkington, B. K. Gibson, C. B. Brook, F. Calura, G. S. Stinson, R. J. Thacker, L. Michel-Dansac, J. Bailin, H. M. P. Couchman, J. Wadsley, T. R. Quinn, A. Maccio

Abstract: We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass --… ▽ More We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted. △ Less

Submitted 21 May, 2012; originally announced May 2012.

Comments: Accepted for publication in MNRAS; 12 pages; 7 figures