Showing 1–41 of 41 results for author: Simpson, C M

An ever-present $Gaia$ snail shell triggered by a dark matter wake

Authors: Robert J. J. Grand, Rüdiger Pakmor, Francesca Fragkoudi, Facundo A. Gómez, Wilma Trick, Christine M. Simpson, Freeke van de Voort, Rebekka Bieri

Abstract: We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation - called Superstars - to simulate a Milky Way-like galaxy with $\gtrsim10^8$ star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane… ▽ More We utilize a novel numerical technique to model star formation in cosmological simulations of galaxy formation - called Superstars - to simulate a Milky Way-like galaxy with $\gtrsim10^8$ star particles to study the formation and evolution of out-of-equilibrium stellar disc structures in a full cosmological setting. In the plane defined by the coordinate and velocity perpendicular to the mid-plane (vertical phase space, $\{Z,V_Z\}$), stars in Solar-like volumes at late times exhibit clear spirals qualitatively similar in shape and amplitude to the $Gaia$ ``Snail shell'' phase spiral. We show that the phase spiral forms at a look back time of $\sim 6$ Gyr during the pericentric passage of a $\sim10^{10}$ $\rm M_{\odot}$ satellite on a polar orbit. This satellite stimulates the formation of a resonant wake in the dark matter halo while losing mass at a rate of $\sim0.5$-$1$ dex per orbit loop. The peak magnitude of the wake-induced gravitational torque at the Solar radius is $\sim 8$ times that from the satellite, and triggers the formation of a disc warp that wraps up into a vertical phase spiral over time. As the wake decays, the phase spiral propagates several Gigayears to present-day and can be described as ``ever-present'' once stable disc evolution is established. These results suggest an alternative scenario to explain the $Gaia$ phase spiral which does not rely on a perturbation from bar buckling or a recent direct hit from a satellite. △ Less

Submitted 28 June, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: accepted for pub in mnras. 11 + 5 pages. For a high-cadence animated version of figure 1, see https://wwwmpa.mpa-garching.mpg.de/auriga/movies/multi_halo_6_sf64.mp4

Colour and infall time distributions of satellite galaxies in simulated Milky-Way analogs

Authors: Yue Pan, Christine M. Simpson, Andrey Kravtsov, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Viraj Manwadkar, Clarke J. Esmerian

Abstract: We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulati… ▽ More We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulations, which was originally designed to model Milky Way-like host galaxies, resemble the populations in the Exploration of Local VolumE Satellites (ELVES) Survey and the Satellites Around Galactic Analogs (SAGA) survey in their luminosity function in the luminosity range $-12 \lesssim M_V \lesssim -15$ and resemble ELVES in their quenched fraction and colour--magnitude distribution in the luminosity range $-12 \lesssim M_g \lesssim -15$. We find that satellites transition from blue colours to red colours at the luminosity range $-15 \lesssim M_g \lesssim -12$ in both the simulations and observations and we show that this shift is driven by environmental effects in the simulations. We demonstrate also that the colour distribution in both simulations and observations can be decomposed into two statistically distinct populations based on their morphological type or star-forming status that are statistically distinct. In the simulations, these two populations also have statistically distinct infall time distributions. The comparison presented here seems to indicate that the tension between the quenched fraction in SAGA and simulations is resolved by the improved target selection of ELVES, but there are still tensions in understanding the colours of faint galaxies, of which ELVES appears to have a significant population of faint blue satellites not recovered in Auriga. △ Less

Submitted 30 January, 2023; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 15 pages, 10 figures. Published in MNRAS on December 15, 2022

Machine Learning for Galactic Archaeology: A chemistry-based neural network method for identification of accreted disc stars

Authors: Thorold Tronrud, Patricia B. Tissera, Facundo A. Gómez, Robert J. J. Grand, Ruediger Pakmor, Federico Marinacci, Christine M. Simpson

Abstract: We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that th… ▽ More We develop a method ('Galactic Archaeology Neural Network', GANN) based on neural network models (NNMs) to identify accreted stars in galactic discs by only their chemical fingerprint and age, using a suite of simulated galaxies from the Auriga Project. We train the network on the target galaxy's own local environment defined by the stellar halo and the surviving satellites. We demonstrate that this approach allows the detection of accreted stars that are spatially mixed into the disc. Two performance measures are defined - recovery fraction of accreted stars, and the probability that a star with a positive (accreted) classification is a true-positive result, P(TP). As the NNM output is akin to an assigned probability, we are able to determine positivity based on flexible threshold values that can be adjusted easily to refine the selection of presumed-accreted stars. We find that GANN identifies accreted disc stars within simulated galaxies, with high recovery fraction and/or high P(TP). We also find that stars in Gaia-Enceladus-Sausage (GES) mass systems are over 50% recovered by our NNMs in the majority (18/24) of cases. Additionally, nearly every individual source of accreted stars is detected at 10% or more of its peak stellar mass in the disc. We also demonstrate that a conglomerated NNM, trained on the halo and satellite stars from all of the Auriga galaxies provides the most consistent results, and could prove to be an intriguing future approach as our observational capabilities expand. △ Less

Submitted 19 July, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 19 pages, 12 figures

How Cosmic Rays Mediate the Evolution of the Interstellar Medium

Authors: Christine M. Simpson, Rüdiger Pakmor, Christoph Pfrommer, Simon C. O. Glover, Rowan Smith

Abstract: We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or e… ▽ More We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or explosions in regions that have been cleared by prior SN, stellar winds, or radiation. We explore various mixed models by adjusting these fractions relative to each other. We find that in the simple system of a periodic stratified gas layer the ISM structure will evolve to one of two solutions: a "peak driving" state where warm gas is volume filling or a "thermal runaway" state where hot gas is volume filling. CR pressure and transport are important factors that strongly influence the solution state the ISM reaches and have the ability to flip the ISM between solutions. Observable signatures such as gamma ray emission and HI gas are explored. We find that gamma ray luminosity from pion decay is largely consistent with observations for a range of model parameters. The thickness of the HI gas layer may be too compact, however, this may be due to a large cold neutral fraction of midplane gas. The volume fraction of hot gas evolves to stable states in both solutions, but neither settles to a Milky Way-like configuration, suggesting that additional physics which is omitted here (e.g. a cosmological circum-galactic medium, radiation transport, or spectrally resolved and spatially varying CR transport) may be required. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 27 pages, 21 Figures, submitted to MNRAS, comments welcome

Formation and fate of low metallicity stars in IllustrisTNG50

Authors: Ruediger Pakmor, Christine M. Simpson, Freeke van de Voort, Lars Hernquist, Lieke van Son, Martyna Chruślińska, Rebekka Bieri, Selma E. de Mink, Volker Springel

Abstract: Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in… ▽ More Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: 15 pages, 11 figures, accepted by MNRAS, comments welcome

Simulating radio synchrotron emission in star-forming galaxies: small-scale magnetic dynamo and the origin of the far infrared-radio correlation

Authors: Christoph Pfrommer, Maria Werhahn, Rüdiger Pakmor, Philipp Girichidis, Christine M. Simpson

Abstract: In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shoc… ▽ More In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shock, which injects turbulence and drives a small-scale magnetic dynamo. As the shock propagates outwards and the associated turbulence decays, the large velocity shear between the supersonically rotating cool disc with respect to the (partially) pressure-supported hot circumgalactic medium excites Kelvin-Helmholtz surface and body modes. Those inject turbulence and drive multiple small-scale dynamos, which exponentially amplify magnetic fields. They grow in scale to reach equipartition with thermal and CR energies in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The radio-synchrotron emission is dominated by primary electrons, irradiates the magnetised disc, bulge, and bubble-shaped magnetically-loaded outflows of our simulated Milky Way-mass galaxy. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR and (iii) a varying radio intensity with galactic inclination. (abridged) △ Less

Submitted 5 July, 2022; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 37 pages, 26 figures, accepted for publication in MNRAS

Satellites Around Milky Way Analogs: Tension in the Number and Fraction of Quiescent Satellites Seen in Observations Versus Simulations

Authors: Ananthan Karunakaran, Kristine Spekkens, Kyle A. Oman, Christine M. Simpson, Azadeh Fattahi, David J. Sand, Paul Bennet, Denija Crnojević, Carlos S. Frenk, Facundo A. Gómez, Robert J. J. Grand, Michael G. Jones, Federico Marinacci, Burçin Mutlu-Pakdil, Julio F. Navarro, Dennis Zaritsky

Abstract: We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APO… ▽ More We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We use archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APOSTLE and Auriga. We compare our detection rates from the NUV and FUV bands to the SAGA-II H$α$ detections and find that they are broadly consistent with over $85\%$ of observed satellites detected in all three tracers. We apply the same spatial selection criteria used around SAGA-II hosts to select satellites around the MW-like hosts in APOSTLE and Auriga. We find very good overall agreement in the derived SFRs for the star-forming satellites as well as the number of star-forming satellites per host in observed and simulated samples. However, the number and fraction of quenched satellites in the SAGA-II sample are significantly lower than those in APOSTLE and Auriga below a stellar mass of $M_*\sim10^{8}\,M_{\odot}$, even when the SAGA-II incompleteness and interloper corrections are included. This discrepancy is robust with respect to the resolution of the simulations and persists when alternative star-formation tracers are employed. We posit that this disagreement is not readily explained by vagaries in the observed or simulated samples considered here, suggesting a genuine discrepancy that may inform the physics of satellite populations around MW analogs. △ Less

Submitted 19 July, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: Accepted for publication in ApJ Letters after minor changes to text

Linking the brightest stellar streams with the accretion history of Milky Way-like galaxies

Authors: Alex Vera-Casanova, Facundo A. Gómez, Antonela Monachesi, Ignacio Gargiulo, Diego Pallero, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Carlos S. Frenk, Gustavo Morales

Abstract: According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest… ▽ More According to the current galaxy formation paradigm, mergers and interactions play an important role in shaping present-day galaxies. The remnants of this merger activity can be used to constrain galaxy formation models. In this work we use a sample of thirty hydrodynamical simulations of Milky Way-mass halos, from the AURIGA project, to generate surface brightness maps and search for the brightest stream in each halo as a function of varying limiting magnitude. We find that none of the models shows signatures of stellar streams at $μ_{r}^{lim} \leq 25$ mag arcsec$^{-2}$. The stream detection increases significantly between 27 and 28 mag arcsec$^{-2}$. Nevertheless, even at 30 mag arcsec$^{-2}$, 13 percent of our models show no detectable streams. We study the properties of the brightest streams progenitors (BSPs). We find that BSPs are accreted within a broad range of infall times, from 1.6 to 10 Gyr ago, with only 25 percent accreted within the last 5 Gyrs; thus most BSPs correspond to relatively early accretion events. We also find that 37 percent of the BSPs survive to the present day. The median infall times for surviving and disrupted BSPs are 5.6 and 6.7 Gyr, respectively. We find a clear relation between infall time and infall mass of the BSPs, such that more massive progenitors tend to be accreted at later times. However, we find that the BSPs are not, in most cases, the dominant contributor to the accreted stellar halo of each galaxy. △ Less

Submitted 9 June, 2022; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS, 15 pages, 11 figures

Determining the full satellite population of a Milky Way-mass halo in a highly resolved cosmological hydrodynamic simulation

Authors: Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Ashley J. Kelly, Facundo A. Gómez, Adrian Jenkins, Volker Springel, Carlos S. Frenk, Simon D. M. White

Abstract: We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution $m_b =$ 800 $\rm M_{\odot}$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic vari… ▽ More We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magneto-hydrodynamic cosmological zoom-in simulation (baryonic mass resolution $m_b =$ 800 $\rm M_{\odot}$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are: i) robust to stochastic variations associated with the so-called ``Butterfly Effect''; and ii) well converged over 3.5 orders of magnitude in mass resolution. We find that there are approximately five times as many satellite galaxies at this high resolution compared to a standard ($m_b\sim 10^{4-5}\, \rm M_{\odot}$) resolution simulation of the same system. This is primarily because 2/3rds of the high resolution satellites do not form at standard resolution. A smaller fraction (1/6th) of the satellites present at high resolution form and disrupt at standard resolution; these objects are preferentially low-mass satellites on intermediate- to low-eccentricity orbits with impact parameters $\lesssim 30$ kpc. As a result, the radial distribution of satellites becomes substantially more centrally concentrated at higher resolution, in better agreement with recent observations of satellites around Milky Way-mass haloes. Finally, we show that our galaxy formation model successfully forms ultra-faint galaxies and reproduces the stellar velocity dispersion, half-light radii, and $V$-band luminosities of observed Milky Way and Local Group dwarf galaxies across 6 orders of magnitude in luminosity ($10^3$-$10^{9}$ $\rm L_{\odot}$). △ Less

Submitted 3 September, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: accepted for publication in MNRAS. Main changes include new figure 4 and figure 11

A tale of two populations: surviving and destroyed dwarf galaxies and the build up of the Milky Way's stellar halo

Authors: Azadeh Fattahi, Alis J. Deason, Carlos S. Frenk, Christine M. Simpson, Facundo A. Gomez, Robert J. J. Grand, Antonela Monachesi, Federico Marinacci, Ruediger Pakmor

Abstract: We use magneto-hydrodynamical simulations of Milky Way-mass haloes from the Auriga project to examine the properties of surviving and destroyed dwarf galaxies that are accreted by these haloes over cosmic time. We show that the combined luminosity function of surviving and destroyed dwarfs at infall is similar in the various Auriga haloes, and is dominated by the destroyed dwarfs. There is, howeve… ▽ More We use magneto-hydrodynamical simulations of Milky Way-mass haloes from the Auriga project to examine the properties of surviving and destroyed dwarf galaxies that are accreted by these haloes over cosmic time. We show that the combined luminosity function of surviving and destroyed dwarfs at infall is similar in the various Auriga haloes, and is dominated by the destroyed dwarfs. There is, however, a strong dependence on infall time: destroyed dwarfs have typically early infall times, $t_{infall}<6$ Gyr, whereas the majority of dwarfs accreted at $t_{infall}>10$ Gyr have survived to the present day. Because of their late infall the surviving satellites today had higher metallicites at infall than their destroyed counterparts of similar infall mass; the difference is even more pronounced for the present-day metallicites of satellites, many of which continue to form stars after infall. In agreement with previous work, we find that a small number of relatively massive destroyed dwarf galaxies dominate the mass of the stellar haloes. However, there is a significant radial dependence: while 90 per cent of the mass in the inner regions ($<\,20\,$kpc) is contributed, on average, by only 3 massive progenitors, the outer regions ($>\,100\,$kpc) typically have $\sim8$ main progenitors of relatively lower mass. Finally, we show that a few massive progenitors dominate the metallicity distribution of accreted stars, even at the metal poor end. Contrary to common assumptions in the literature, dwarf galaxies of mass $M_{*}<10^7 \, M_{\odot}$ make up less than 10 per cent of the accreted, metal poor stars ([Fe/H] $<\,-3$) in the inner $50\,$kpc. △ Less

Submitted 7 October, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: Replaced with the accepted version after minor modifications

The Edge of the Galaxy

Authors: Alis J. Deason, Azadeh Fattahi, Carlos S. Frenk, Robert J. J. Grand, Kyle A. Oman, Shea Garrison-Kimmel, Christine M. Simpson, Julio F. Navarro

Abstract: We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a s… ▽ More We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc. △ Less

Submitted 25 June, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: 14 pages, 11 figures. MNRAS in press

Magnetising the circumgalactic medium of disk galaxies

Authors: Ruediger Pakmor, Freeke van de Voort, Rebekka Bieri, Facundo A. Gomez, Robert J. J. Grand, Thomas Guillet, Federico Marinacci, Christoph Pfrommer, Christine M. Simpson, Volker Springel

Abstract: The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom s… ▽ More The circumgalactic medium (CGM) is one of the frontiers of galaxy formation and intimately connected to the galaxy via accretion of gas on to the galaxy and gaseous outflows from the galaxy. Here we analyse the magnetic field in the CGM of the Milky Way-like galaxies simulated as part of the \textsc{Auriga} project that constitutes a set of high resolution cosmological magnetohydrodynamical zoom simulations. We show that before $z=1$ the CGM becomes magnetised via galactic outflows that transport magnetised gas from the disk into the halo. At this time the magnetisation of the CGM closely follows its metal enrichment. We then show that at low redshift an in-situ turbulent dynamo that operates on a timescale of Gigayears further amplifies the magnetic field in the CGM and saturates before $z=0$. The magnetic field strength reaches a typical value of $0.1\,μG$ at the virial radius at $z=0$ and becomes mostly uniform within the virial radius. Its Faraday rotation signal is in excellent agreement with recent observations. For most of its evolution the magnetic field in the CGM is an unordered small scale field. Only strong coherent outflows at low redshift are able to order the magnetic field in parts of the CGM that are directly displaced by these outflows. △ Less

Submitted 7 September, 2020; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 14 pages, 15 figures, accepted by MNRAS, updated Fig. 15 in proofs with errorbars

The Milky Way total mass profile as inferred from Gaia DR2

Authors: Marius Cautun, Alejandro Benitez-Llambay, Alis J. Deason, Carlos S. Frenk, Azadeh Fattahi, Facundo A. Gómez, Robert J. J. Grand, Kyle A. Oman, Julio F. Navarro, Christine M. Simpson

Abstract: We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r<20 kpc. We provide an analytic expression that relates the baryonic d… ▽ More We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r<20 kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2 and 4 at radial distances of 20, 8 and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well. The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.97_{-0.19}^{+0.24}\times10^{12} M_\odot$, and concentration before baryon contraction of $9.4_{-2.6}^{+1.9}$, which lie close to the median halo mass--concentration relation predicted in $Λ$CDM. The inferred total mass, $M_{200}^{\rm total}=1.08_{-0.14}^{+0.20} \times 10^{12} M_\odot$, is in good agreement with recent measurements. The model gives a MW stellar mass of $5.04_{-0.52}^{+0.43}\times10^{10} M_\odot$ and infers that the DM density at the Solar position is $ρ_{\odot}^{\rm DM}=8.8_{-0.5}^{+0.5}\times10^{-3} M_\odot \rm{pc}^{-3}\equiv0.33_{-0.02}^{+0.02}~\rm{GeV}~\rm{cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo. △ Less

Submitted 16 April, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 21 pages, 14 figures, main results: figs 3, 7, 10 and 11; MNRAS accepted version, code and resources available at https://github.com/MariusCautun/Milky_Way_mass_profile

Simulating cosmological substructure in the solar neighbourhood

Authors: Christine M. Simpson, Ignacio Gargiulo, Facundo A. Gómez, Robert J. J. Grand, Nicolás Maffione, Andrew P. Cooper, Alis J. Deason, Carlos Frenk, John Helly, Federico Marinacci, Rüdiger Pakmor

Abstract: We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These str… ▽ More We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyrs, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the Aurigaia catalogues and show that the current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods. △ Less

Submitted 12 September, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS letters

The effects of dynamical substructure on Milky Way mass estimates from the high velocity tail of the local stellar halo

Authors: Robert J. J. Grand, Alis J. Deason, Simon D. M. White, Christine M. Simpson, Facundo A. Gómez, Federico Marinacci, Ruediger Pakmor

Abstract: We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase… ▽ More We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise and truncation of the high velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by $\sim 20 \%$, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of $1.29 ^{+0.37}_{-0.47} \times 10^{12}$ $\rm M_{\odot}$. △ Less

Submitted 3 June, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: Accepted to the Monthly Notices of the Royal Astronomical Society, Letters, with minor revision

The star formation histories of dwarf galaxies in Local Group cosmological simulations

Authors: Ruth Digby, Julio F. Navarro, Azadeh Fattahi, Christine M. Simpson, Kyle A. Oman, Facundo A. Gomez, Carlos S. Frenk, Robert J. J. Grand, Ruediger Pakmor

Abstract: We use the APOSTLE and Auriga cosmological simulations to study the star formation histories (SFHs) of field and satellite dwarf galaxies. Despite sizeable galaxy-to-galaxy scatter, the SFHs of APOSTLE and Auriga dwarfs exhibit robust average trends with galaxy stellar mass: faint field dwarfs ($10^5<M_{\rm star}/M_\odot<10^{6.5}$) have, on average, steadily declining SFHs, whereas brighter dwarfs… ▽ More We use the APOSTLE and Auriga cosmological simulations to study the star formation histories (SFHs) of field and satellite dwarf galaxies. Despite sizeable galaxy-to-galaxy scatter, the SFHs of APOSTLE and Auriga dwarfs exhibit robust average trends with galaxy stellar mass: faint field dwarfs ($10^5<M_{\rm star}/M_\odot<10^{6.5}$) have, on average, steadily declining SFHs, whereas brighter dwarfs ($10^{7.5}<M_{\rm star}/M_\odot<10^{9}$) show the opposite trend. Intermediate-mass dwarfs have roughly constant SFHs. Satellites exhibit similar average trends, but with substantially suppressed star formation in the most recent $\sim 5$ Gyr, likely as a result of gas loss due to tidal and ram-pressure stripping after entering the haloes of their primaries. These simple mass and environmental trends are in good agreement with the derived SFHs of Local Group (LG) dwarfs whose photometry reaches the oldest main sequence turnoff. SFHs of galaxies with less deep data show deviations from these trends, but this may be explained, at least in part, by the large galaxy-to-galaxy scatter, the limited sample size, and the large uncertainties of the inferred SFHs. Confirming the predicted mass and environmental trends will require deeper photometric data than currently available, especially for isolated dwarfs. △ Less

Submitted 12 March, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: 18 pages, 10 figures, submitted to MNRAS

Efficacy of early stellar feedback in low gas surface density environments

Authors: Rahul Kannan, Federico Marinacci, Christine M. Simpson, Simon C. O. Glover, Lars Hernquist

Abstract: We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($Σ_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$)… ▽ More We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($Σ_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$). We show that radiation fields decrease the star formation rate and therefore the total stellar mass formed by a factor of $\sim 2$. This increases the gas depletion timescale and brings the simulated Kennicutt-Schmidt relation closer to the observational estimates. Radiation feedback coupled with SNe is more efficient at driving outflows with the mass and energy loading increasing by a factor of $\sim 10$. This increase is mainly driven by the additional entrainment of medium density ($10^{-2} \leq n< 1 \ {\rm cm}^{-3}$), warm ($300 \ {\rm K}\leq T<8000 \ {\rm K}$) material. Therefore including radiation fields tends to launch colder, denser and higher mass and energy loaded outflows. This is because photoheating of the high density gas around a newly formed star over-pressurises the region, causing it to expand. This reduces the ambient density in which the SNe explode by a factor of $10-100$ which in turn increases their momentum output by a factor of $\sim 1.5-2.5$. Finally, we note that in these low gas surface density environments, radiation fields primarily impact the ISM via photoheating and radiation pressure has only a minimal role in regulating star formation. △ Less

Submitted 4 December, 2018; originally announced December 2018.

Comments: 17 pages, 18 figures, Submitted to MNRAS

The velocity anisotropy of the Milky Way satellite system

Authors: Alexander H. Riley, Azadeh Fattahi, Andrew B. Pace, Louis E. Strigari, Carlos S. Frenk, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Julio F. Navarro, Rüdiger Pakmor, Christine M. Simpson, Simon D. M. White

Abstract: We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $β$, as a function of Galactocentric distance and compare the MW sa… ▽ More We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $β$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $β\sim -2$ at $r\sim20$ kpc to $β\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $β(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $β(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system. △ Less

Submitted 5 August, 2019; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 16 pages, 8 figures. Key results are summarized in Figures 6 and 7. Accepted to MNRAS

No cores in dark matter-dominated dwarf galaxies with bursty star formation histories

Authors: Sownak Bose, Carlos S. Frenk, Adrian Jenkins, Azadeh Fattahi, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Julio F. Navarro, Kyle A. Oman, Ruediger Pakmor, Joop Schaye, Christine M. Simpson, Volker Springel

Abstract: Measurements of the rotation curves of dwarf galaxies are often interpreted as requiring a constant density core at the centre, at odds with the "cuspy" inner profiles predicted by $N$-body simulations of cold dark matter (CDM) haloes. It has been suggested that this conflict could be resolved by fluctuations in the inner gravitational potential caused by the periodic removal of gas following burs… ▽ More Measurements of the rotation curves of dwarf galaxies are often interpreted as requiring a constant density core at the centre, at odds with the "cuspy" inner profiles predicted by $N$-body simulations of cold dark matter (CDM) haloes. It has been suggested that this conflict could be resolved by fluctuations in the inner gravitational potential caused by the periodic removal of gas following bursts of star formation. Earlier work has suggested that core formation requires a bursty and extended star formation history (SFH). Here we investigate the structure of CDM haloes of dwarf galaxies ($M_{\rm DM} \sim 10^9-5\times10^{10}\,{\rm M}_\odot$) formed in the APOSTLE ('A Project of Simulating the Local Environment') and AURIGA cosmological hydrodynamic simulations. Our simulations have comparable or better resolution than others that make cores ($M_{\rm gas} \sim 10^4\,{\rm M}_\odot$, gravitational softening $\sim 150$ pc). Yet, we do not find evidence of core formation at {\it any} mass or any correlation between the inner slope of the DM density profile and temporal variations in the SFH. APOSTLE and AURIGA dwarfs display a similar diversity in their cumulative SFHs to available data for Local Group dwarfs. Dwarfs in both simulations are DM-dominated on all resolved scales at all times, likely limiting the ability of gas outflows to alter significantly the central density profiles of their haloes. We conclude that recurrent bursts of star formation are not sufficient to cause the formation of cores, and that other conditions must also be met for baryons to be able to modify the central DM cusp. △ Less

Submitted 30 May, 2019; v1 submitted 8 October, 2018; originally announced October 2018.

Comments: 15 pages, 9 figures, 2 tables. Updated version accepted for publication in MNRAS

Faraday rotation maps of disk galaxies

Authors: Ruediger Pakmor, Thomas Guillet, Christoph Pfrommer, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Christine M. Simpson, Volker Springel

Abstract: Faraday rotation is one of the most widely used observables to infer the strength and configuration of the magnetic field in the ionised gas of the Milky Way and nearby spiral galaxies. Here we compute synthetic Faraday rotation maps at $z=0$ for a set of disk galaxies from the Auriga high-resolution cosmological simulations, for different observer positions within and outside the galaxy. We find… ▽ More Faraday rotation is one of the most widely used observables to infer the strength and configuration of the magnetic field in the ionised gas of the Milky Way and nearby spiral galaxies. Here we compute synthetic Faraday rotation maps at $z=0$ for a set of disk galaxies from the Auriga high-resolution cosmological simulations, for different observer positions within and outside the galaxy. We find that the strength of the Faraday rotation of our simulated galaxies for a hypothetic observer at the solar circle is broadly consistent with the Faraday rotation seen for the Milky Way. The same holds for an observer outside the galaxy and the observed signal of the nearby spiral galaxy M51. However, we also find that the structure and angular power spectra of the synthetic all-sky Faraday rotation maps vary strongly with azimuthal position along the solar circle. We argue that this variation is a result of the structure of the magnetic field of the galaxy that is dominated by an azimuthal magnetic field ordered scales of several kpc, but has radial and vertical magnetic field components that are only ordered on scales of 1-2 kpc. Because the magnetic field strength decreases exponentially with height above the disk, the Faraday rotation for an observer at the solar circle is dominated by the local environment. This represents a severe obstacle for attempts to reconstruct the global magnetic field of the Milky Way from Faraday rotation maps alone without including additional observables. △ Less

Submitted 19 September, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Comments: 10 pages, 10 figures, accepted by MNRAS

Aurigaia: mock Gaia DR2 stellar catalogues from the Auriga cosmological simulations

Authors: Robert J. J. Grand, John Helly, Azadeh Fattahi, Marius Cautun, Shaun Cole, Andrew P. Cooper, Alis J. Deason, Carlos Frenk, Facundo A. Gómez, Jason A. S. Hunt, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Volker Springel, Dandan Xu

Abstract: We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $Λ$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at… ▽ More We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $Λ$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $Λ$CDM cosmological model. △ Less

Submitted 5 September, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS (1st of September, 2018) after improvements to several figures and some discussion. The mock data are available for access and download at http://icc.dur.ac.uk/data and documentation and supplementary material is available at http://auriga.h-its.org/gaiamock.html

The Auriga Stellar Haloes: Connecting stellar population properties with accretion and merging history

Authors: Antonela Monachesi, Facundo A. Gómez, Robert J. J. Grand, Christine M. Simpson, Guinevere Kauffmann, Sebastián Bustamante, Federico Marinacci, Rüdiger Pakmor, Volker Springel, Carlos S. Frenk, Simon D. M. White, Patricia B. Tissera

Abstract: We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and… ▽ More We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies. △ Less

Submitted 21 February, 2019; v1 submitted 20 April, 2018; originally announced April 2018.

Comments: Accepted to MNRAS. 30 pages, 19 figures

On the relevance of chaos for halo stars in the solar neighbourhood II

Authors: Nicolas P. Maffione, Facundo A. Gómez, Pablo M. Cincotta, Claudia M. Giordano, Robert J. J. Grand, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Volker Springel, Carlos S. Frenk

Abstract: In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of c… ▽ More In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of chaotic mixing using multicomponent Galactic potential models and solar neighbourhood-like volumes extracted from fully cosmological hydrodynamic simulations, thus naturally accounting for the gravitational potential associated with baryonic components, such as the bulge and disc. Despite the strong change in the global Galactic potentials with respect to those obtained in dark matter only simulations, our results confirm that a large fraction of halo particles evolving on chaotic orbits exhibit their chaotic behaviour after periods of time significantly larger than a Hubble time. In addition, significant diffusion in phase space is not observed on those particles that do exhibit chaotic behaviour within a Hubble time. △ Less

Submitted 15 May, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

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

The multiplicity and anisotropy of galactic satellite accretion

Authors: Shi Shao, Marius Cautun, Carlos S. Frenk, Robert J. J. Grand, Facundo A. Gómez, Federico Marinacci, Christine M. Simpson

Abstract: We study the incidence of group and filamentary dwarf galaxy accretion into Milky Way (MW) mass haloes using two types of hydrodynamical simulations: EAGLE, which resolves a large cosmological volume, and the AURIGA suite, which are very high resolution zoom-in simulations of individual MW-sized haloes. The present-day 11 most massive satellites are predominantly (75%) accreted in single events, 1… ▽ More We study the incidence of group and filamentary dwarf galaxy accretion into Milky Way (MW) mass haloes using two types of hydrodynamical simulations: EAGLE, which resolves a large cosmological volume, and the AURIGA suite, which are very high resolution zoom-in simulations of individual MW-sized haloes. The present-day 11 most massive satellites are predominantly (75%) accreted in single events, 14% in pairs and 6% in triplets, with higher group multiplicities being unlikely. Group accretion becomes more common for fainter satellites, with 60% of the top 50 satellites accreted singly, 12% in pairs, and 28% in richer groups. A group similar in stellar mass to the Large Magellanic Cloud (LMC) would bring on average 15 members with stellar mass larger than $10^4{~\rm M_\odot}$. Half of the top 11 satellites are accreted along the two richest filaments. The accretion of dwarf galaxies is highly anisotropic, taking place preferentially perpendicular to the halo minor axis, and, within this plane, preferentially along the halo major axis. The satellite entry points tend to be aligned with the present-day central galaxy disc and satellite plane, but to a lesser extent than with the halo shape. Dwarfs accreted in groups or along the richest filament have entry points that show an even larger degree of alignment with the host halo than the full satellite population. We also find that having most satellites accreted as a single group or along a single filament is unlikely to explain the MW disc of satellites. △ Less

Submitted 27 February, 2018; v1 submitted 14 December, 2017; originally announced December 2017.

Comments: 15 pages, 17 figures. Minor changes. MNRAS in press

The dependence of cosmic ray driven galactic winds on halo mass

Authors: Svenja Jacob, Rüdiger Pakmor, Christine M. Simpson, Volker Springel, Christoph Pfrommer

Abstract: Galactic winds regulate star formation in disk galaxies and help to enrich the circum-galactic medium. They are therefore crucial for galaxy formation, but their driving mechanism is still poorly understood. Recent studies have demonstrated that cosmic rays (CRs) can drive outflows if active CR transport is taken into account. Using hydrodynamical simulations of isolated galaxies with virial masse… ▽ More Galactic winds regulate star formation in disk galaxies and help to enrich the circum-galactic medium. They are therefore crucial for galaxy formation, but their driving mechanism is still poorly understood. Recent studies have demonstrated that cosmic rays (CRs) can drive outflows if active CR transport is taken into account. Using hydrodynamical simulations of isolated galaxies with virial masses between $10^{10}$ and $10^{13}\mathrm{~M_\odot}$, we study how the properties of CR-driven winds depend on halo mass. CRs are treated in a two-fluid approximation and their transport is modelled through isotropic or anisotropic diffusion. We find that CRs are only able to drive mass-loaded winds beyond the virial radius in haloes with masses below $10^{12}\mathrm{~M_\odot}$. For our lowest examined halo mass, the wind is roughly spherical and has velocities of $\sim20\mathrm{~km\;s^{-1}}$. With increasing halo mass, the wind becomes biconical and can reach ten times higher velocities. The mass loading factor drops rapidly with virial mass, a dependence that approximately follows a power-law with a slope between $-1$ and $-2$. This scaling is slightly steeper than observational inferences, and also steeper than commonly used prescriptions for wind feedback in cosmological simulations. The slope is quite robust to variations of the CR injection efficiency or the CR diffusion coefficient. In contrast to the mass loading, the energy loading shows no significant dependence on halo mass. While these scalings are close to successful heuristic models of wind feedback, the CR-driven winds in our present models are not yet powerful enough to fully account for the required feedback strength. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Simulating Gamma-ray Emission in Star-forming Galaxies

Authors: C. Pfrommer, R. Pakmor, C. M. Simpson, V. Springel

Abstract: Star forming galaxies emit GeV- and TeV-gamma rays that are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. To understand the emission, we have used the moving mesh code Arepo to perform magneto-hydrodynamical galaxy formation simulations with self-consistent CR physics. Our galaxy models exhibit a first burst of star formation that injects C… ▽ More Star forming galaxies emit GeV- and TeV-gamma rays that are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. To understand the emission, we have used the moving mesh code Arepo to perform magneto-hydrodynamical galaxy formation simulations with self-consistent CR physics. Our galaxy models exhibit a first burst of star formation that injects CRs at supernovae. Once CRs have sufficiently accumulated in our Milky-Way like galaxy, their buoyancy force overcomes the magnetic tension of the toroidal disk field. As field lines open up, they enable anisotropically diffusing CRs to escape into the halo and to accelerate a bubble-like, CR-dominated outflow. However, these bubbles are invisible in our simulated gamma-ray maps of hadronic pion-decay and secondary inverse-Compton emission because of low gas density in the outflows. By adopting a phenomenological relation between star formation rate (SFR) and far-infrared emission and assuming that gamma rays mainly originate from decaying pions, our simulated galaxies can reproduce the observed tight relation between far-infrared and gamma-ray emission, independent of whether we account for anisotropic CR diffusion. This demonstrates that uncertainties in modeling active CR transport processes only play a minor role in predicting gamma-ray emission from galaxies. We find that in starbursts, most of the CR energy is "calorimetrically" lost to hadronic interactions. In contrast, the gamma-ray emission deviates from this calorimetric property at low SFRs due to adiabatic losses, which cannot be identified in traditional one-zone models. △ Less

Submitted 15 September, 2017; originally announced September 2017.

Comments: Accepted by ApJL, 8 pages, 4 figures

Origin of chemically distinct discs in the Auriga cosmological simulations

Authors: Robert J. J. Grand, Sebastián Bustamante, Facundo A. Gómez, Daisuke Kawata, Federico Marinacci, Rüdiger Pakmor, Hans-Walter Rix, Christine M. Simpson, Martin Sparre, Volker Springel

Abstract: The stellar disk of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a cle… ▽ More The stellar disk of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the $[\rm α/Fe]$ -- $[\rm Fe/H]$ plane. We find two main pathways to creating this dichotomy which operate in different regions of the galaxies: a) an early ($z>1$) and intense high-$\rm[α/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[α/Fe]$ star formation; and b) an early phase of high-$\rm[α/Fe]$ star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-$\rm[α/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way. △ Less

Submitted 20 November, 2017; v1 submitted 25 August, 2017; originally announced August 2017.

Comments: 12 pages, MNRAS accepted

Quenching and ram pressure stripping of simulated Milky Way satellite galaxies

Authors: Christine M. Simpson, Robert J. J. Grand, Facundo A. Gómez, Federico Marinacci, Rüdiger Pakmor, Volker Springel, David J. R. Campbell, Carlos S. Frenk

Abstract: We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems… ▽ More We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems beyond 600 kpc are only strongly quenched below a stellar mass of $10^7$ M$_\odot$. Ram pressure stripping appears to be the dominant quenching mechanism and 50% of quenched systems cease star formation within 1 Gyr of first infall. We demonstrate that systems within a host galaxy's $R_{200}$ radius are comprised of two populations: (i) a first infall population that has entered the host halo within the past few Gyrs and (ii) a population of returning `backsplash' systems that have had a much more extended interaction with the host. Backsplash galaxies that do not return to the host galaxy by redshift zero exhibit quenching properties similar to galaxies within $R_{200}$ and are distinct from other external systems. The simulated quenching trend with stellar mass has some tension with observations, but our simulations are able reproduce the range of quenching times measured from resolved stellar populations of Local Group dwarf galaxies. △ Less

Submitted 20 March, 2018; v1 submitted 8 May, 2017; originally announced May 2017.

Comments: 22 pages, 17 figures, accepted for publication in MNRAS

Lessons from the Auriga discs: The hunt for the Milky Way's ex-situ disc is not yet over

Authors: Facundo A. Gómez, Robert J. J. Grand, Antonela Monachesi, Simon D. M. White, Sebastian Bustamante, Federico Marinacci, Rüdiger Pakmor, Christine M. Simpson, Volker Springel, Carlos S. Frenk

Abstract: We characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex-situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population… ▽ More We characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex-situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population, which is strongly dominated by stars born in-situ. One third of our models shows a significant ex-situ disc but this fraction would be larger if constraints on orbital circularity were relaxed. Most of the ex-situ material ($\gtrsim 50\%$) comes from single massive satellites ($> 6 \times 10^{10}~M_{\odot}$). These satellites are accreted with a wide range of infall times and inclination angles (up to $85^{\circ}$). Ex-situ discs are thicker, older and more metal-poor than their in-situ counterparts. They show a flat median age profile, which differs from the negative gradient observed in the in-situ component. As a result, the likelihood of identifying an ex-situ disc in samples of old stars on near-circular orbits increases towards the outskirts of the disc. We show three examples that, in addition to ex-situ discs, have a strongly rotating dark matter component. Interestingly, two of these ex-situ stellar discs show an orbital circularity distribution that is consistent with that of the in-situ disc. Thus, they would not be detected in typical kinematic studies. △ Less

Submitted 26 April, 2017; originally announced April 2017.

Comments: 13 pages, 11 figures, submitted to MNRAS

Magnetic field formation in the Milky Way-like disk galaxies of the Auriga project

Authors: Ruediger Pakmor, Facundo A. Gomez, Robert J. J. Grand, Federico Marinacci, Christine M. Simpson, Volker Springel, David J. R. Campbell, Carlos S. Frenk, Thomas Guillet, Christoph Pfrommer, Simon D. M. White

Abstract: The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogethe… ▽ More The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{μG}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution. △ Less

Submitted 2 May, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 16 pages, 16 figures, accepted by MNRAS, section 4.1 (turbulent dynamo) extended substantially, now includes magnetic and kinetic power spectra

The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test

Authors: Ji-hoon Kim, Oscar Agertz, Romain Teyssier, Michael J. Butler, Daniel Ceverino, Jun-Hwan Choi, Robert Feldmann, Ben W. Keller, Alessandro Lupi, Thomas Quinn, Yves Revaz, Spencer Wallace, Nickolay Y. Gnedin, Samuel N. Leitner, Sijing Shen, Britton D. Smith, Robert Thompson, Matthew J. Turk, Tom Abel, Kenza S. Arraki, Samantha M. Benincasa, Sukanya Chakrabarti, Colin DeGraf, Avishai Dekel, Nathan J. Goldbaum , et al. (18 additional authors not shown)

Abstract: Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st… ▽ More Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes. △ Less

Submitted 8 March, 2018; v1 submitted 10 October, 2016; originally announced October 2016.

Comments: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

Journal ref: Astrophysical Journal, 833, 202, 2016

GMC Evolution in a Barred Spiral Galaxy with Star Formation and Thermal Feedback

Authors: Yusuke Fujimoto, Greg L. Bryan, Elizabeth J. Tasker, Asao Habe, Christine M. Simpson

Abstract: We explore the impact of star formation and thermal stellar feedback on the giant molecular cloud (GMC) population forming in a M83-type barred spiral galaxy. We compare three high-resolution simulations (1.5 pc cell size) with different star formation/feedback models: one with no star formation, one with star formation but no feedback, and one with star formation and thermal energy injection. We… ▽ More We explore the impact of star formation and thermal stellar feedback on the giant molecular cloud (GMC) population forming in a M83-type barred spiral galaxy. We compare three high-resolution simulations (1.5 pc cell size) with different star formation/feedback models: one with no star formation, one with star formation but no feedback, and one with star formation and thermal energy injection. We analyze the resulting population of clouds, finding that we can identify the same population of massive, virialized clouds and transient, low-surface density clouds found in our previous work (that did not include star formation or feedback). Star formation and feedback can affect the mix of clouds we identify. In particular, star formation alone simply converts dense cloud gas into stars with only a small change to the cloud populations, principally resulting in a slight decrease in the transient population. Feedback, however, has a stronger impact: while it is not generally sufficient to entirely destroy the clouds, it does eject gas out of them, increasing the gas density in the inter-cloud region. This decreases the number of massive clouds, but substantially increases the transient cloud population. We also find that feedback tends to drive a net radial inflow of massive clouds, leading to an increase in the star formation rate in the bar region. We examine a number of possible reasons for this and conclude that it is possible that the drag force from the enhanced intercloud density could be responsible. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: 19 pages, 16 figures, accepted for publication in MNRAS

The role of cosmic ray pressure in accelerating galactic outflows

Authors: Christine M. Simpson, Ruediger Pakmor, Federico Marinacci, Christoph Pfrommer, Volker Springel, Simon C. O. Glover, Paul C. Clark, Rowan J. Smith

Abstract: We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are… ▽ More We study the formation of galactic outflows from supernova explosions (SNe) with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SNe placement and energy feedback, including cosmic rays (CR), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium. △ Less

Submitted 16 August, 2016; v1 submitted 7 June, 2016; originally announced June 2016.

Comments: 8 pages, 4 figures, accepted for publication in ApJL; movie of simulated gas densities can be found here: http://www.h-its.org/tap-images/galactic-outflows/

Galactic winds driven by isotropic and anisotropic cosmic ray diffusion in disk galaxies

Authors: Ruediger Pakmor, Christoph Pfrommer, Christine M. Simpson, Volker Springel

Abstract: The physics of cosmic rays (CR) is a promising candidate for explaining the driving of galactic winds and outflows. Recent galaxy formation simulations have demonstrated the need for active CR transport either in the form of diffusion or streaming to successfully launch winds in galaxies. However, due to computational limitations, most previous simulations have modeled CR transport isotropically.… ▽ More The physics of cosmic rays (CR) is a promising candidate for explaining the driving of galactic winds and outflows. Recent galaxy formation simulations have demonstrated the need for active CR transport either in the form of diffusion or streaming to successfully launch winds in galaxies. However, due to computational limitations, most previous simulations have modeled CR transport isotropically. Here, we discuss high resolution simulations of isolated disk galaxies in a $10^{11}\rm{M_\odot}$ halo with the moving mesh code {\sc Arepo} that include injection of CRs from supernovae, advective transport, CR cooling, and CR transport through isotropic or anisotropic diffusion. We show that either mode of diffusion leads to the formation of strong bipolar outflows. However, they develop significantly later in the simulation with anisotropic diffusion compared to the simulation with isotropic diffusion. Moreover, we find that isotropic diffusion allows most of the CRs to quickly diffuse out of the disk, while in the simulation with anisotropic diffusion, most CRs remain in the disk once the magnetic field becomes dominated by its azimuthal component, which occurs after $\sim 300\,{\rm Myrs}$. This has important consequences for the gas dynamics in the disk. In particular, we show that isotropic diffusion strongly suppresses the amplification of the magnetic field in the disk compared to anisotropic or no diffusion models. We therefore conclude that reliable simulations which include CR transport inevitably need to account for anisotropic diffusion. △ Less

Submitted 16 July, 2016; v1 submitted 2 May, 2016; originally announced May 2016.

Comments: 7 pages, 5 figures, accepted by ApJL

Semi-implicit anisotropic cosmic ray transport on an unstructured moving mesh

Authors: Ruediger Pakmor, Christoph Pfrommer, Christine M. Simpson, Rahul Kannan, Volker Springel

Abstract: In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport of cosmic rays in the form of diffusion or streaming not only affects the gas dynamics but also rearranges the magnetic fields themselves. The coupled dynamics o… ▽ More In the interstellar medium of galaxies and the intracluster gas of galaxy clusters, the charged particles making up cosmic rays are moving almost exclusively along (but not across) magnetic field lines. The resulting anisotropic transport of cosmic rays in the form of diffusion or streaming not only affects the gas dynamics but also rearranges the magnetic fields themselves. The coupled dynamics of magnetic fields and cosmic rays can thus impact the formation and evolution of galaxies and the thermal evolution of galaxy clusters in critical ways. Numerically studying these effects requires solvers for anisotropic diffusion that are accurate, efficient, and robust, requirements that have proven difficult to satisfy in practice. Here, we present an anisotropic diffusion solver on an unstructured moving mesh that is conservative, does not violate the entropy condition, allows for semi-implicit time integration with individual timesteps, and only requires solving a single linear system of equations per timestep. We apply our new scheme to a large number of test problems and show that it works as well or better than previous implementations. Finally, we demonstrate for a numerically demanding simulation of the formation of an isolated disk galaxy that our local time-stepping scheme reproduces the results obtained with global time-stepping at a fraction of the computational cost. △ Less

Submitted 16 July, 2016; v1 submitted 28 April, 2016; originally announced April 2016.

Comments: 16 pages, 16 figures, accepted by MNRAS

Simulating cosmic ray physics on a moving mesh

Authors: C. Pfrommer, R. Pakmor, K. Schaal, C. M. Simpson, V. Springel

Abstract: We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magneto-hydrodynamics (MHD) on an unstructured moving mesh, as realised in the massively parallel AREPO code for cosmological simulations. We account for diffusive shock acceleration of CRs at resolved shocks and at supernova remnants in the interstellar medium (ISM), and follow the advective CR transport within… ▽ More We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magneto-hydrodynamics (MHD) on an unstructured moving mesh, as realised in the massively parallel AREPO code for cosmological simulations. We account for diffusive shock acceleration of CRs at resolved shocks and at supernova remnants in the interstellar medium (ISM), and follow the advective CR transport within the magnetised plasma, as well as anisotropic diffusive transport of CRs along the local magnetic field. CR losses are included in terms of Coulomb and hadronic interactions with the thermal plasma. We demonstrate the accuracy of our formalism for CR acceleration at shocks through simulations of plane-parallel shock tubes that are compared to newly derived exact solutions of the Riemann shock tube problem with CR acceleration. We find that the increased compressibility of the post-shock plasma due to the produced CRs decreases the shock speed. However, CR acceleration at spherically expanding blast waves does not significantly break the self-similarity of the Sedov-Taylor solution; the resulting modifications can be approximated by a suitably adjusted, but constant adiabatic index. In first applications of the new CR formalism to simulations of isolated galaxies and cosmic structure formation, we find that CRs add an important pressure component to the ISM that increases the vertical scale height of disk galaxies, and thus reduces the star formation rate. Strong external structure formation shocks inject CRs into the gas, but the relative pressure of this component decreases towards halo centres as adiabatic compression favours the thermal over the CR pressure. △ Less

Submitted 24 November, 2016; v1 submitted 25 April, 2016; originally announced April 2016.

Comments: Accepted by MNRAS, 31 pages, 10 figures

Kinetic energy from supernova feedback in high-resolution galaxy simulations

Authors: Christine M. Simpson, Greg L. Bryan, Cameron Hummels, Jeremiah P. Ostriker

Abstract: We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our applicat… ▽ More We describe a new method for adding a prescribed amount of kinetic energy to simulated gas modeled on a cartesian grid by directly altering grid cells' mass and velocity in a distributed fashion. The method is explored in the context of supernova feedback in high-resolution ($\sim 10$ pc) hydrodynamic simulations of galaxy formation. Resolution-dependence is a primary consideration in our application of the method and simulations of isolated explosions (performed at different resolutions) motivate a resolution-dependent scaling for the injected fraction of kinetic energy that we apply in cosmological simulations of a $10^9$ Msun dwarf halo. We find that in high density media ($\gtrsim$ 50 cm$^{-3}$) with coarse resolution ($\gtrsim 4$ pc per cell), results are sensitive to the initial kinetic energy fraction due to early and rapid cooling. In our galaxy simulations, the deposition of small amounts of supernova energy in kinetic form (as little as 1%) has a dramatic impact on the evolution of the system, resulting in an order of magnitude suppression of stellar mass. The overall behavior of the galaxy in the two highest resolution simulations we perform appears to converge. We discuss the resulting distribution of stellar metallicities, an observable sensitive to galactic wind properties, and find that while the new method demonstrates increased agreement with observed systems, significant discrepancies remain, likely due to simplistic assumptions that neglect contributions from Type Ia supernovae and stellar winds. △ Less

Submitted 26 June, 2015; v1 submitted 14 October, 2014; originally announced October 2014.

Comments: 23 pages, 16 figures, accepted for publication in ApJ

Diffuse gas properties and stellar metallicities in cosmological simulations of disc galaxy formation

Authors: Federico Marinacci, Ruediger Pakmor, Volker Springel, Christine M. Simpson

Abstract: We analyse the properties of the circum-galactic medium and the metal content of the stars comprising the central galaxy in eight hydrodynamical `zoom-in' simulations of disc galaxy formation. We use these properties as a benchmark for our model of galaxy formation physics implemented in the moving-mesh code AREPO, which succeeds in forming quite realistic late-type spirals in the set of `Aquarius… ▽ More We analyse the properties of the circum-galactic medium and the metal content of the stars comprising the central galaxy in eight hydrodynamical `zoom-in' simulations of disc galaxy formation. We use these properties as a benchmark for our model of galaxy formation physics implemented in the moving-mesh code AREPO, which succeeds in forming quite realistic late-type spirals in the set of `Aquarius' initial conditions of Milky Way-sized haloes. Galactic winds significantly influence the morphology of the circum-galactic medium and induce bipolar features in the distribution of heavy elements. They also affect the thermodynamic properties of the circum-galactic gas by supplying an energy input that sustains its radiative losses. Although a significant fraction of the heavy elements are transferred from the central galaxy to the halo, and even beyond the virial radius, enough metals are retained by stars to yield a peak in their metallicity distributions at about $Z_{\odot}$. All our default runs overestimate the stellar [O/Fe] ratio, an effect that we demonstrate can be rectified by an increase of the adopted SN type Ia rate. Nevertheless, the models have difficulty in producing stellar metallicity gradients of the same strength as observed in the Milky Way. △ Less

Submitted 6 June, 2014; v1 submitted 19 March, 2014; originally announced March 2014.

Comments: 18 pages, 10 figures, 1 table, accepted for publication in MNRAS. Added 1 table, updated figures 1, 2, 4 and 6 and a few text changes to match the accepted version

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

Authors: The Enzo Collaboration, Greg L. Bryan, Michael L. Norman, Brian W. O'Shea, Tom Abel, John H. Wise, Matthew J. Turk, Daniel R. Reynolds, David C. Collins, Peng Wang, Samuel W. Skillman, Britton Smith, Robert P. Harkness, James Bordner, Ji-hoon Kim, Michael Kuhlen, Hao Xu, Nathan Goldbaum, Cameron Hummels, Alexei G. Kritsuk, Elizabeth Tasker, Stephen Skory, Christine M. Simpson, Oliver Hahn, Jeffrey S. Oishi , et al. (4 additional authors not shown)

Abstract: This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gra… ▽ More This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology. △ Less

Submitted 8 July, 2013; originally announced July 2013.

Comments: 60 pages, 22 figures, submitted to ApJS. The code is available at http://enzo-project.org

The effect of feedback and reionization on star formation in low-mass dwarf galaxy haloes

Authors: Christine M. Simpson, Greg L. Bryan, Kathryn V. Johnston, Britton D. Smith, Mordecai-Mark Mac Low, Sanjib Sharma, Jason Tumlinson

Abstract: We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range… ▽ More We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc. △ Less

Submitted 18 March, 2013; v1 submitted 5 November, 2012; originally announced November 2012.

Comments: 26 pages, 14 figures, accepted for publication in MNRAS; additional explanation added to sections 2,3 and 4; Figure 4 updated; results unchanged

Measuring Transverse Motions for Nearby Galaxy Clusters

Authors: Erika T. Hamden, Christine M. Simpson, Kathryn V. Johnston, Duane M. Lee

Abstract: Measuring the full three-dimensional motions of extra-galactic objects in the Universe presents a seemingly insurmountable challenge. In this paper we investigate the application of a technique to measure tangential motion that has previously only been applied nearby within the Local Group of galaxies, to clusters of galaxies far beyond its borders. We show that mapping the mean line-of-sight moti… ▽ More Measuring the full three-dimensional motions of extra-galactic objects in the Universe presents a seemingly insurmountable challenge. In this paper we investigate the application of a technique to measure tangential motion that has previously only been applied nearby within the Local Group of galaxies, to clusters of galaxies far beyond its borders. We show that mapping the mean line-of-sight motion throughout a galaxy cluster could in principle be used to detect the "perspective rotation" induced by the projection of the cluster's tangential motion into the line-of-sight. The signal will be most prominent for clusters of the largest angular extent, most symmetric intrinsic velocity distribution and surveyed with the largest number of pointings possible. We investigate the feasibility of detecting this signal using three different approaches: measuring line-of-sight motions of individual cluster members; taking spectra of intracluster gas; and mapping distortions of the Cosmic Microwave Background radiation. We conclude that future spectroscopic surveys of 1000's of members of nearby galaxy clusters hold the most promise of measuring cluster tangential motions using this technique. △ Less

Submitted 24 May, 2010; originally announced May 2010.

Comments: Accepted for publication in ApJ Letters