Search

Scientific workflows have become integral tools in broad scientific computing use cases. Science discovery is increasingly dependent on workflows to orchestrate large and complex scientific experiments that range from execution of a cloud-based data preprocessing pipeline to multi-facility instrument-to-edge-to-HPC computational workflows. Given the changing landscape of scientific computing and the evolving needs of emerging scientific applications, it is paramount that the development of novel scientific workflows and system functionalities seek to increase the efficiency, resilience, and pervasiveness of existing systems and applications. Specifically, the proliferation of machine learning/artificial intelligence (ML/AI) workflows, need for processing large scale datasets produced by instruments at the edge, intensification of near real-time data processing, support for long-term experiment campaigns, and emergence of quantum computing as an adjunct to HPC, have significantly changed the functional and operational requirements of workflow systems. Workflow systems now need to, for example, support data streams from the edge-to-cloud-to-HPC enable the management of many small-sized files, allow data reduction while ensuring high accuracy, orchestrate distributed services (workflows, instruments, data movement, provenance, publication, etc.) across computing and user facilities, among others. Further, to accelerate science, it is also necessary that these systems implement specifications/standards and APIs for seamless (horizontal and vertical) integration between systems and applications, as well as enabling the publication of workflows and their associated products according to the FAIR principles. This document reports on discussions and findings from the 2022 international edition of the Workflows Community Summit that took place on November 29 and 30, 2022.

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., Comment: 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

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., Comment: 15 pages, 10 figures. Published in MNRAS on December 15, 2022

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., Comment: 19 pages, 12 figures

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., Comment: 27 pages, 21 Figures, submitted to MNRAS, comments welcome

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$., Comment: 15 pages, 11 figures, accepted by MNRAS, comments welcome

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), Comment: 37 pages, 26 figures, accepted for publication in MNRAS

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$\alpha$ 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., Comment: Accepted for publication in ApJ Letters after minor changes to text

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 $\mu_{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., Comment: Accepted for publication in MNRAS, 15 pages, 11 figures

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}$)., Comment: accepted for publication in MNRAS. Main changes include new figure 4 and figure 11

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., Comment: Replaced with the accepted version after minor modifications

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., Comment: 14 pages, 11 figures. MNRAS in press

We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sampled and its density evaluated, then so can the complex distribution. Our first application to gravitational waves uses an autoregressive flow, conditioned on detector strain data, to map a multivariate standard normal distribution into the posterior distribution over system parameters. We train the model on artificial strain data consisting of IMRPhenomPv2 waveforms drawn from a five-parameter $(m_1, m_2, \phi_0, t_c, d_L)$ prior and stationary Gaussian noise realizations with a fixed power spectral density. This gives performance comparable to current best deep-learning approaches to gravitational-wave parameter estimation. We then build a more powerful latent variable model by incorporating autoregressive flows within the variational autoencoder framework. This model has performance comparable to Markov chain Monte Carlo and, in particular, successfully models the multimodal $\phi_0$ posterior. Finally, we train the autoregressive latent variable model on an expanded parameter space, including also aligned spins $(\chi_{1z}, \chi_{2z})$ and binary inclination $\theta_{JN}$, and show that all parameters and degeneracies are well-recovered. In all cases, sampling is extremely fast, requiring less than two seconds to draw $10^4$ posterior samples., Comment: 14 pages, 7 figures

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\,\mu 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., Comment: 14 pages, 15 figures, accepted by MNRAS, updated Fig. 15 in proofs with errorbars

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 $\Lambda$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 $\rho_{\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., Comment: 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

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}$., Comment: Accepted to the Monthly Notices of the Royal Astronomical Society, Letters, with minor revision

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., Comment: 6 pages, 4 figures, accepted for publication in MNRAS letters

We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin., Comment: 7 pages, 3 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

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

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 ($\Sigma_{\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., Comment: 17 pages, 18 figures, Submitted to MNRAS

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, $\beta$, 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 $\beta\sim -2$ at $r\sim20$ kpc to $\beta\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 $\beta(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 $\beta(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., Comment: 16 pages, 8 figures. Key results are summarized in Figures 6 and 7. Accepted to MNRAS

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., Comment: 15 pages, 9 figures, 2 tables. Updated version accepted for publication in MNRAS

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., Comment: 10 pages, 10 figures, accepted by MNRAS

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 $\Lambda$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 $\Lambda$CDM cosmological model., Comment: 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

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., Comment: Accepted to MNRAS. 30 pages, 19 figures

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., Comment: 17 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

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., Comment: 15 pages, 17 figures. Minor changes. MNRAS in press

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., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

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 \alpha/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[\alpha/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[\alpha/Fe]$ star formation; and b) an early phase of high-$\rm[\alpha/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[\alpha/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way., Comment: 12 pages, MNRAS accepted

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., Comment: 22 pages, 17 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 11 figures, submitted to MNRAS

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{\mu 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., Comment: 16 pages, 16 figures, accepted by MNRAS, section 4.1 (turbulent dynamo) extended substantially, now includes magnetic and kinetic power spectra

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., Comment: 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/

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., Comment: 19 pages, 16 figures, accepted for publication in MNRAS

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., Comment: 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/

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., Comment: 7 pages, 5 figures, accepted by ApJL

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., Comment: 16 pages, 16 figures, accepted by MNRAS

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., Comment: 23 pages, 16 figures, accepted for publication in ApJ