Search

We introduce the public version of the cosmological magnetohydrodynamical moving-mesh simulation code Arepo. This version contains a finite-volume magnetohydrodynamics algorithm on an unstructured, dynamic Voronoi tessellation coupled to a tree-particle-mesh algorithm for the Poisson equation either on a Newtonian or cosmologically expanding spacetime. Time-integration is performed adopting local timestep constraints for each cell individually, solving the fluxes only across active interfaces, and calculating gravitational forces only between active particles, using an operator-splitting approach. This allows simulations with high dynamic range to be performed efficiently. Arepo is a massively distributed-memory parallel code, using the Message Passing Interface (MPI) communication standard and employing a dynamical work-load and memory balancing scheme to allow optimal use of multi-node parallel computers. The employed parallelization algorithms of Arepo are deterministic and produce binary-identical results when re-run on the same machine and with the same number of MPI ranks. A simple primordial cooling and star formation model is included as an example of sub-resolution models commonly used in simulations of galaxy formation. Arepo also contains a suite of computationally inexpensive test problems, ranging from idealized tests for automated code verification to scaled-down versions of cosmological galaxy formation simulations, and is extensively documented in order to assist adoption of the code by new scientific users., 39 pages, 15 figures, accepted in ApJS, https://www.arepo-code.org, repository: https://gitlab.mpcdf.mpg.de/vrs/arepo

We use the Illustris simulation to study the distributions of ellipticities and kinematic misalignments of galactic projections, as well as the intrinsic shapes and rotation of the simulated galaxies. Our results for the projections of galaxies display clear trends of an overall increase of kinematic misalignment and a slight decrease of ellipticity for fast rotators with increasing masses, while revealing no clear distinction between slow rotators of different mass. It is also found that the number of very slow rotators with large ellipticities is much larger than found in observations. The intrinsic properties of the galaxies are then analysed. The results for the intrinsic shapes of the galaxies are mostly consistent with previous results inferred from observational data. The distributions of intrinsic misalignment of the galaxies suggest that some of the galaxies produced by Illustris have significant rotation around their medium axes. Further analysis reveals that most of these galaxies display signs of non-equilibrium. We then study the evolution of the intrinsic misalignments and shapes of three specific Illustris galaxies, which we consider as typical ones, along the main progenitor line of their merger trees, revealing how mergers influence the intrinsic shapes and kinematics: the spin axis in general stays close to the shortest axis, and tends to quickly relax to such an equilibrium state within a few dynamical times of the galaxy after major perturbations; triaxiality and intrinsic flatness in general decrease with time, however, sometimes increases occur that are clearly seen to correlate with major merger events., 16 pages, 18 figures

We present a new cosmological zoom-in simulation, where the zoom region consists of two halos with virial mass M_v~5x10^{12}M_{sun} and a ~Mpc long cosmic filament connecting them at z~2. Using this simulation, we study the evolution of the intergalactic medium in between these two halos at unprecedented resolution. At 5>z>3, the two halos are found to lie in a large intergalactic sheet, or "pancake", consisting of multiple co-planar dense filaments along which nearly all halos with M_v>10^9M_{sun} are located. This sheet collapses at z~5 from the merger of two smaller sheets. The strong shock generated by this merger leads to thermal instabilities in the post-shock region, and to a shattering of the sheet resulting in ~2x10^4K and densities of n>~10^{-3}cm^{-3}, which are pressure confined in a hot medium with T~10^6K and n>~10^{-5}cm^{-3}. When the sheet is viewed face on, these cold clouds have neutral hydrogen column densities of N_{HI}>10^{17.2}cm^{-2}, making them detectable as Lyman limit systems, though they lie well outside the virial radius of any halo and even well outside the dense filaments. Their chemical composition is pristine, having zero metalicity, similar to several recently observed systems. Since these systems form far from any galaxies, these results are robust to galaxy formation physics, resulting purely from the collapse of large scale structure and radiative cooling, provided sufficient spatial resolution is available., Comment: 7 pages, 5 figures. Resubmitted to ApJL following first referee report. Movie generalizing figure 2 available with ancillary material. Additional comments welcome!

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto \sigma^b I^c$ with $R$, $\sigma$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log \sigma\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=\delta$, where $\delta \sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5\sigma^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations., Comment: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars ( 1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.

We present a new cosmological, magnetohydrodynamical simulation for galaxy formation: TNG50, the third and final installment of the IllustrisTNG project. TNG50 evolves 2x2160^3 dark-matter particles and gas cells in a volume 50 comoving Mpc across. It hence reaches a numerical resolution typical of zoom-in simulations, with a baryonic element mass of 8.5x10^4 Msun and an average cell size of 70-140 parsecs in the star-forming regions of galaxies. Simultaneously, TNG50 samples ~700 (6,500) galaxies with stellar masses above 10^10 (10^8) Msun at z=1. Here we investigate the structural and kinematical evolution of star-forming galaxies across cosmic time (0 < z < 6). We quantify their sizes, disk heights, 3D shapes, and degree of rotational vs. dispersion-supported motions as traced by rest-frame V-band light (i.e. roughly stellar mass) and by Halpha light (i.e. star-forming and dense gas). The unprecedented resolution of TNG50 enables us to model galaxies with sub-kpc half-light radii and with 2-3, z1.5). Despite large velocity dispersion at high redshift, cold and dense gas in galaxies predominantly arranges in disky or elongated shapes at all times and masses; these gaseous components exhibit rotationally-dominated motions far exceeding the collisionless stellar bodies., MNRAS. Highlights: Figures 9, 14, 15. See companion paper by Nelson et al. 2019b. Visuals at http://www.tng-project.org

Cosmological hydrodynamical simulations of galaxy formation are a powerful theoretical tool, and enable us to directly calculate the observable signatures resulting from the complex process of cosmic structure formation. Here we present early results from the ongoing TNG50 run, an unprecedented ‘next generation’ cosmological, magnetohydrodynamical simulation—the third and final volume of the IllustrisTNG project, with over 20 billion resolution elements and capturing spatial scales of \(\sim \)100 parsecs. It incorporates the comprehensive TNG model for galaxy formation physics, and we here describe the simulation scope, novel achievements, and early investigations on resolved galactic and halo structural properties.

Understanding how baryonic processes shape the intracluster medium (ICM) is of critical importance to the next generation of galaxy cluster surveys. However, most models of structure formation neglect potentially important physical processes, like anisotropic thermal conduction (ATC). In this letter, we explore the impact of ATC on the prevalence of cool-cores (CCs) using 12 pairs of magnetohydrodynamical galaxy cluster simulations, simulated using the IllustrisTNG model with and without ATC. Although the impact of ATC varies from cluster to cluster and with CC criterion, its inclusion produces a systematic shift to larger CC fractions at z = 0 for all CC criteria considered. Additionally, the inclusion of ATC yields a flatter CC fraction redshift evolution, easing the tension with the observed evolution. With ATC included, the energy required for the central black hole to achieve self-regulation is reduced and the gas fraction in the cluster core increases, resulting in larger CC fractions. ATC makes the ICM unstable to perturbations and the increased efficiency of AGN feedback suggests that its inclusion results in a greater level of mixing in the ICM. Therefore, ATC is potentially an important physical process in reproducing the thermal structure of the ICM., Comment: 5 pages, 4 figures, submitted to MNRAS, comments welcome

Galaxies are surrounded by massive gas reservoirs (i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated Milky Way-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L⋆ galaxies is extremely diverse: column densities of commonly observed species span ∼3 − 4 dex and their covering fractions range from ${\sim } 5$ to $90{{\ \rm per\ cent}}$. Despite this diversity, we identify the following correlations: 1) the covering fractions (CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H i, C iv, and Si ii anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H i, C iv, and Si ii positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L⋆ galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.

Using the IllustrisTNG simulations, we investigate the connection between galaxy morphology and star formation in central galaxies with stellar masses in the range $10^9-10^{11.5}~\mathrm{M}_{\odot}$. We quantify galaxy morphology by a kinematical decomposition of the stellar component into a spheroidal and a disc component (spheroid-to-total ratio, S/T) and by the concentration of the stellar mass density profile ($C_{82}$). S/T is correlated with stellar mass and star-formation activity, while $C_{82}$ correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and $C_{82}$. Low and high mass galaxies are dominated by random stellar motion, while only intermediate-mass galaxies ($M_{\star}\approx10^{10}-10^{10.5}~\mathrm{M}_{\odot}$) are dominated by ordered rotation. Whereas higher-mass galaxies are typical spheroids with high concentrations, lower-mass galaxies have low concentration, pointing to different formation channels. Although we find a correlation between S/T and star-formation activity, in the TNG model galaxies do not necessarily change their morphology when they transition through the green valley or when they cease their star formation, this depending on galaxy stellar mass and morphological estimator. Instead, the morphology (S/T and $C_{82}$) is generally set during the star-forming phase of galaxies. The apparent correlation between S/T and star formation arises because earlier-forming galaxies had, on average, a higher S/T at a given stellar mass. Furthermore, we show that mergers drive in-situ bulge formation in intermediate-mass galaxies and are responsible for the recent spheroidal mass assembly in the massive galaxies with $M_{\star}>10^{11}~\mathrm{M}_{\odot}$. In particular, these massive galaxies assemble about half of the spheroidal mass while star-forming and the other half through mergers while quiescent., 20 pages with 14 figures (+ appendix). Accepted by MNRAS

Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the Inter-Galactic Medium (IGM). About $90 \%$ of all material has been ejected and re-accreted in galactic winds at least once. The vast majority of smoothly accreted gas enters into a galactic fountain that extends to a median galactocentric distance of $\sim 20$ kpc with a median recycling timescale of $\sim 500$ Myr. We demonstrate that, in most cases, galactic fountains acquire angular momentum via mixing of low-angular momentum, wind-recycled gas with high-angular momentum gas in the Circum-Galactic Medium (CGM). Prograde mergers boost this activity by helping to align the disc and CGM rotation axes, whereas retrograde mergers cause the fountain to lose angular momentum. Fountain flows that promote angular momentum growth are conducive to smooth evolution on tracks quasi-parallel to the disc sequence of the stellar mass-specific angular momentum plane, whereas retrograde minor mergers, major mergers and bar-driven secular evolution move galaxies towards the bulge-sequence. Finally, we demonstrate that fountain flows act to flatten and narrow the radial metallicity gradient and metallicity dispersion of disc stars, respectively. Thus, the evolution of galactic fountains depends strongly on the cosmological merger history and is crucial for the chemo-dynamical evolution of Milky Way-sized disc galaxies., Comment: MNRAS; accepted. Changes include: clarifications and additional discussion, particularly of the fountain flow statistics shown in Fig. 7; and a resolution study that shows the results are converged for two resolution levels (~10^5 and ~10^4 Msun per baryonic element)

The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyse the distribution and evolution of metals within the IllustrisTNG simulation suite with a focus on the gas-phase mass-metallicity relation (MZR). We find that the IllustrisTNG model broadly reproduces the slope and normalization evolution of the MZR across the redshift range 0 < z < 2 and mass range 109 < M*/M☉ < 1010.5. We make predictions for the high-redshift (2 < z < 10) metal content of galaxies which is described by a gradual decline in the normalization of the metallicity with an average high- redshift (z > 2) evolution fit by d log(Z)/dz ≈ -0.064. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals (̃85 per cent at z = 0) live outside of the ISM, either in an extended gas disc, the circumgalactic medium, or outside the halo. Nevertheless, the redshift evolution in the simulated MZR normalization is driven by the higher gas fractions of high-redshift galaxies, not by changes to the metal retention efficiency. The scatter in the simulated MZR contains a clear correlation with the gas-mass or star formation rate of the system, in agreement with the observed fundamental metallicity relation. The scatter in the MZR is driven by a competition between periods of enrichment- and accretion-dominated metallicity evolution. We expect that while the normalization of the MZR declines with redshift, the slope of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the gas fraction dependence of `regulator' style models allows them to simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR.

It is well established that the properties of supermassive black holes and their host galaxies are correlated through scaling relations. While hydrodynamical cosmological simulations have begun to account for the co-evolution of BHs and galaxies, they typically have neglected the BH spin, even though it may play an important role in modulating the growth and feedback of BHs. Here we introduce a new sub-grid model for the BH spin evolution in the moving-mesh code {\small AREPO} in order to improve the physical faithfulness of the BH modelling in galaxy formation simulations. We account for several different channels of spin evolution, in particular gas accretion through a Shakura-Sunyaev $\alpha$-disc, chaotic accretion, and BH mergers. For BH feedback, we extend the IllustrisTNG model, which considers two different BH feedback modes, a thermal quasar mode for high accretion states and a kinetic mode for low Eddington ratios, with a self-consistent accounting of spin-dependent radiative efficiencies and thus feedback strength. We find that BHs with mass $M_{\rm{bh}}\lesssim 10^{8}\, {\rm M}_\odot$ reach high spin values as they typically evolve in the coherent gas accretion regime. On the other hand, BHs with mass $M_{\rm{bh}}\gtrsim 10^{8}\, {\rm M}_\odot$ have lower spins as BH mergers become more frequent, and their accretion discs fragment due to self-gravity, inducing chaotic accretion. We also explore the hypothesis that the transition between the quasar and kinetic feedback modes is mediated by the accretion mode of the BH disc itself, i.e.~the kinetic feedback mode is activated when the disc enters the self-gravity regime. We find excellent agreement between the galaxy and BH populations for this approach and the fiducial TNG model with no spin evolution. Furthermore, our new approach alleviates a tension in the galaxy morphology -- colour relation of the original TNG model., Comment: 22 pages, 13 figure, Submitted to MNRAS

We present the full public release of all data from the TNG50, TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~1.1 PB, including 2,000 full volume snapshots and ~110,000 high time-resolution subbox snapshots. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets., TNG50 joins TNG100 and TNG300 in public release (1 Feb 2021) at http://www.tng-project.org/data

We analyze the IllustrisTNG simulations to study the mass, volume fraction and phase distribution of gaseous baryons embedded in the knots, filaments, sheets and voids of the Cosmic Web from redshift $z=8$ to redshift $z=0$. We find that filaments host more star-forming gas than knots, and that filaments also have a higher relative mass fraction of gas in this phase than knots. We also show that the cool, diffuse Intergalactic Medium (IGM; $T, Comment: Accepted by MNRAS

We measure power spectrum response functions in the presence of baryonic physical processes using separate universe simulations with the IllustrisTNG galaxy formation model. The response functions describe how the small-scale power spectrum reacts to long-wavelength perturbations and they can be efficiently measured with the separate universe technique by absorbing the effects of the long modes into a modified cosmology. Specifically, we focus on the total first-order matter power spectrum response to an isotropic density fluctuation $R_1(k,z)$, which is fully determined by the logarithmic derivative of the nonlinear matter power spectrum ${\rm dln}P_m(k,z)/{\rm dln}k$ and the growth-only response function $G_1(k,z)$. We find that $G_1(k,z)$ is not affected by the baryonic physical processes in the simulations at redshifts $z < 3$ and on all scales probed ($k \lesssim 15h/{\rm Mpc}$, i.e. length scales $\gtrsim 0.4 {\rm Mpc}/h$). In practice, this implies that the power spectrum fully specifies the baryonic dependence of its response function. Assuming an idealized lensing survey setup, we evaluate numerically the baryonic impact on the squeezed-lensing bispectrum and the lensing super-sample power spectrum covariance, which are given in terms of responses. Our results show that these higher-order lensing statistics can display varying levels of sensitivity to baryonic effects compared to the power spectrum, with the squeezed-bispectrum being the least sensitive. We also show that ignoring baryonic effects on lensing covariances slightly overestimates the error budget (and is therefore conservative from the point of view of parameter error bars) and likely has negligible impact on parameter biases in inference analyses., 15 pages, 6 figures, 1 table; comments welcomed! v2 matches version published in MNRAS

We use the IllustrisTNG (TNG) simulations to explore the galaxy-halo connection as inferred from state-of-the-art cosmological, magnetohydrodynamical simulations. With the high mass resolution and large volume achieved by combining the 100 Mpc (TNG100) and 300 Mpc (TNG300) volumes, we establish the mean occupancy of central and satellite galaxies and their dependence on the properties of the dark matter haloes hosting them. We derive best-fitting HOD parameters from TNG100 and TNG300 for target galaxy number densities of $\bar{n}_g = 0.032\,h^3$Mpc$^{-3}$ and $\bar{n}_g = 0.016\,h^3$Mpc$^{-3}$, respectively, corresponding to a minimum galaxy stellar mass of $M_\star\sim1.9\times10^9\,{\rm M}_\odot$ and $M_\star\sim3.5\times10^9\,{\rm M}_\odot$, respectively, in hosts more massive than $10^{11}\,{\rm M}_\odot$. Consistent with previous work, we find that haloes located in dense environments, with low concentrations, later formation times, and high angular momenta are richest in their satellite population. At low mass, highly-concentrated haloes and those located in overdense regions are more likely to contain a central galaxy. The degree of environmental dependence is sensitive to the definition adopted for the physical boundary of the host halo. We examine the extent to which correlations between galaxy occupancy and halo properties are independent and demonstrate that HODs predicted by halo mass and present-day concentration capture the qualitative dependence on the remaining halo properties. At fixed halo mass, concentration is a strong predictor of the stellar mass of the central galaxy, which may play a defining role in the fate of the satellite population. The radial distribution of satellite galaxies, which exhibits a universal form across a wide range of host halo mass, is described accurately by the best-fit NFW density profile of their host haloes., Comment: 20 pages, 13 figures, 1 table, comments welcome

Feedback by active galactic nuclei (AGNs) is essential for regulating the fast radiative cooling of low-entropy gas at the centers of galaxy clusters and for reducing star formation rates of central ellipticals. The details of self-regulation depend critically on the unknown contents of AGN-inflated bubbles. Observations of the Sunyaev-Zeldovich (SZ) signal of AGN bubbles provide us with the ability to directly measure the lobe electron pressure given a bubble morphology. Here we compute the SZ signal of jet-inflated bubbles in three-dimensional magnetohydrodynamical simulations of the galaxy cluster MS0735.6+7421 with the Arepo code, and compare our synthetic SZ results to inferences obtained with popular modelling approaches. We find that cutting out ellipsoidal bubbles from a double-beta pressure profile only matches the inner bubble edges in the simulations and fails to account for the emission of the shock-enhanced pressure cocoon outside the bubbles. This additional contribution significantly worsens the accuracy of the cut-out method for jets with small inclinations with respect to the line of sight. Also, the kinetic SZ effect of the bubbles, a previously neglected contribution, becomes relevant at these smaller inclinations due to entrainment and mixing of the intracluster medium with low-density jet material. Fortunately, the different signs of the kinetic SZ signal in opposite lobes allow modelling this effect. We present an approximate method to determine the jet inclination, which combines jet power and lifetime estimates, the stand-off distance between jet head and bow shock, and the kinetic SZ effect, thereby helping to correctly infer the bubble contents., 7 pages, 4 figures, submitted to ApJ letters. Comments welcome!

We explore orbit properties of 35 prolate-triaxial galaxies selected from the Illustris cosmological hydrodynamic simulation. We present a detailed study of their orbit families, and also analyse relations between the relative abundance of the orbit families and the spin parameter, triaxiality, the ratio of the angular momentum and the baryon fraction. We find that box orbits dominate the orbit structure for most prolate-triaxial galaxies, especially in the central region. The fraction of irregular orbits in the prolate-triaxial galaxies is small, and it increases with galaxy radius. Both the x-tube and z-tube orbits are important in prolate-triaxial galaxies, especially in the outer regions. The fraction of box orbits for prolate-triaxial galaxies ($0.7, Comment: 12 pages, 13 Figures; accepted to MNRAS

We study the stellar angular momentum of thousands of galaxies in the Illustris cosmological simulation, which captures gravitational and gas dynamics within galaxies, as well as feedback from stars and black holes. We find that the angular momentum of the simulated galaxies matches observations well, and in particular two distinct relations are found for late-type versus early-type galaxies. The relation for late-type galaxies corresponds to the value expected from full conservation of the specific angular momentum generated by cosmological tidal torques. The relation for early-type galaxies corresponds to retention of only ~30% of that, but we find that those early-type galaxies with low angular momentum at z=0 nevertheless reside at high redshift on the late-type relation. Some of them abruptly lose angular momentum during major mergers. To gain further insight, we explore the scaling relations in simulations where the galaxy formation physics is modified with respect to the fiducial model. We find that galactic winds with high mass-loading factors are essential for obtaining the high angular momentum relation typical for late-type galaxies, while AGN feedback largely operates in the opposite direction. Hence, feedback controls the stellar angular momentum of galaxies, and appears to be instrumental for establishing the Hubble sequence., Comment: Updated to match accepted ApJL version (minor modifications). 7 pages, 4 figures

We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an unstructured representation of the magnetic vector potential, making the numerical method simple and computationally efficient. The scheme is implemented in the moving mesh code Arepo. We demonstrate the performance of the approach with simulations of driven MHD turbulence, a magnetized disc galaxy, and a cosmological volume with primordial magnetic field. We compare the outcomes of these experiments to those obtained with a previously implemented Powell divergence-cleaning scheme. While CT and the Powell technique yield similar results in idealized test problems, some differences are seen in situations more representative of astrophysical flows. In the turbulence simulations, the Powell cleaning scheme artificially grows the mean magnetic field, while CT maintains this conserved quantity of ideal MHD. In the disc simulation, CT gives slower magnetic field growth rate and saturates to equipartition between the turbulent kinetic energy and magnetic energy, whereas Powell cleaning produces a dynamically dominant magnetic field. Such difference has been observed in adaptive-mesh refinement codes with CT and smoothed-particle hydrodynamics codes with divergence-cleaning. In the cosmological simulation, both approaches give similar magnetic amplification, but Powell exhibits more cell-level noise. CT methods in general are more accurate than divergence-cleaning techniques, and, when coupled to a moving mesh can exploit the advantages of automatic spatial/temporal adaptivity and reduced advection errors, allowing for improved astrophysical MHD simulations., Comment: 12 pages, 8 figures, mnras accepted

We investigate the impact of f(R) modified gravity on the internal properties of Milky Way sized dark matter halos in a set of cosmological zoom simulations of seven halos from the Aquarius suite, carried out with our code MG-GADGET in the Hu & Sawicki f(R) model. Also, we calculate the fifth forces in ideal NFW-halos as well as in our cosmological simulations and compare them against analytic model predictions for the fifth force inside spherical objects. We find that these theoretical predictions match the forces in the ideal halos very well, whereas their applicability is somewhat limited for realistic cosmological halos. Our simulations show that f(R) gravity significantly affects the dark matter density profile of Milky Way sized objects as well as their circular velocities. In unscreened regions, the velocity dispersions are increased by up to 40% with respect to LCDM for viable f(R) models. This difference is larger than reported in previous works. The Solar circle is fully screened in $f_{R0} = -10^{-6}$ models for Milky Way sized halos, while this location is unscreened for slightly less massive objects. Within the scope of our limited halo sample size, we do not find a clear dependence of the concentration parameter of dark matter halos on $f_{R0}$., 13 pages, 9 figures, published in MNRAS

Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding algorithms due to radiative cooling and star-forming processes, prompting us to develop a number of methodology improvements. We find in Illustris a total shock surface area which is about 1.4 times larger at the present epoch compared to non-radiative runs, and an energy dissipation rate at shocks which is higher by a factor of around 7. Remarkably, shocks with Mach numbers above and below $\mathcal{M}\approx10$ contribute about equally to the total dissipation across cosmic time. This is in sharp contrast to non-radiative simulations, and we demonstrate that a large part of the difference arises due to strong black hole radio-mode feedback in Illustris. We also provide an overview of the large diversity of shock morphologies, which includes complex networks of halo-internal shocks, shocks on to cosmic sheets, feedback shocks due to black holes and galactic winds, as well as ubiquitous accretion shocks. In high redshift systems more massive than $10^{12}\,\mathrm{M}_\odot$ we discover the existence of a double accretion shock pattern in haloes. They are created when gas streams along filaments without being shocked at the outer accretion shock, but then forms a second, roughly spherical accretion shock further inside., Comment: 26 pages, 15 figures, 1 movie, published in mnras

Vertically extended, high velocity dispersion stellar distributions appear to be a ubiquitous feature of disc galaxies, and both internal and external mechanisms have been proposed to be the major driver of their formation. However, it is unclear to what extent each mechanism can generate such a distribution, which is likely to depend on the assembly history of the galaxy. To this end, we perform 16 high resolution cosmological-zoom simulations of Milky Way-sized galaxies using the state-of-the-art cosmological magneto-hydrodynamical code \textlcsc{AREPO}, and analyse the evolution of the vertical kinematics of the stellar disc in connection with various heating mechanisms. We find that the bar is the dominant heating mechanism in most cases, whereas spiral arms, radial migration, and adiabatic heating from mid-plane density growth are all sub-dominant. The strongest source, though less prevalent than bars, originates from external perturbations from satellites/sub-halos of masses log$_{10} (M/\rm M_{\odot}) \gtrsim 10$. However, in many simulations the orbits of newborn star particles become cooler with time, such that they dominate the shape of the age-velocity dispersion relation and overall vertical disc structure unless a strong external perturbation takes place., 24 pages, 15 figures (new figs. 6&7). Accepted in MNRAS. Much improved over originally submitted version

We study how X-rays from stellar binary systems and the hot intracluster medium (ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel implementation of gas cooling in the moving-mesh hydrodynamics code \textsc{arepo}. X-rays from stellar binaries do not affect cooling at all as their emission spectrum is too hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple well with gas in the temperature range $10^4 - 10^6$ K. Idealised simulations show that the hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters because of the high virial temperature ($> 10^7$K), making the interaction between the gas and incident photons very ineffective. Satellite galaxies in cluster environments, on the other hand, experience a high radiation flux due to the emission from the host halo. Low mass satellites ($< 10^{12}\rm{M_\odot}$) in particular have virial temperatures that are exactly in the regime where the effect of the radiation field is maximal. Idealised simulations of satellite galaxies including only the effect of host halo radiation (no ram pressure stripping or tidal effects) fields show a drastic reduction in the amount of cool gas formed ($\sim 40\%$) on a short timescale of about $0.5$ Gyrs. A galaxy merger simulation including all the other environmental quenching mechanisms, shows about $20\%$ reduction in the stellar mass of the satellite and about $\sim 30\%$ reduction in star formation rate after $1$ Gyr due to the host hot halo radiation field. These results indicate that the hot halo radiation fields potentially play an important role in quenching galaxies in cluster environments., 14 pages, 12 figures, submitted to MNRAS

We study how optical galaxy morphology depends on mass and star formation rate (SFR) in the Illustris Simulation. To do so, we measure automated galaxy structures in 10808 simulated galaxies at z=0 with stellar masses 10^9.7 < M_*/M_sun < 10^12.3. We add observational realism to idealized synthetic images and measure non-parametric statistics in rest-frame optical and near-IR images from four directions. We find that Illustris creates a morphologically diverse galaxy population, occupying the observed bulge strength locus and reproducing median morphology trends versus stellar mass, SFR, and compactness. Morphology correlates realistically with rotation, following classification schemes put forth by kinematic surveys. Type fractions as a function of environment agree roughly with data. These results imply that connections among mass, star formation, and galaxy structure arise naturally from models matching global star formation and halo occupation functions when simulated with accurate methods. This raises a question of how to construct experiments on galaxy surveys to better distinguish between models. We predict that at fixed halo mass near 10^12 M_sun, disc-dominated galaxies have higher stellar mass than bulge-dominated ones, a possible consequence of the Illustris feedback model. While Illustris galaxies at M_* ~ 10^11 M_sun have a reasonable size distribution, those at M_* ~ 10^10 M_sun have half-light radii larger than observed by a factor of two. Furthermore, at M_* ~ 10^10.5-10^11 M_sun, a relevant fraction of Illustris galaxies have distinct "ring-like" features, such that the bright pixels have an unusually wide spatial extent., 26 pages, 15 figures, MNRAS accepted version

We characterize the photometric and kinematic properties of simulated early-type galaxy (ETG) stellar halos, and compare them to observations. We select a sample of ~1200 ETGs in the TNG100 and TNG50 simulations, spanning a stellar mass range of $10^{10.3}-10^{12}M_{\odot}$ and within the range of (g-r) colour and lambda-ellipticity diagram populated by observed ETGs. We determine photometric parameters, intrinsic shapes, and kinematic observables in their extended stellar halos. We study the variation in kinematics from center to halo and connect it to a change in the intrinsic shape of the galaxies. We find that the simulated galaxy sample reproduces the diversity of kinematic properties observed in ETG halos. Simulated fast rotators (FRs) divide almost evenly in one third having flat lambda profiles and high halo rotational support, a third with gently decreasing profiles, and another third with low halo rotation. Slow rotators (SRs) tend to have increased rotation in the outskirts, with half of them exceeding lambda=0.2. For $M_{*}>10^{11.5}M_{\odot}$ halo rotation is unimportant. A similar variety of properties is found for the stellar halo intrinsic shapes. Rotational support and shape are deeply related: the kinematic transition to lower rotational support is accompanied by a change towards rounder intrinsic shape. Triaxiality in the halos of FRs increases outwards and with stellar mass. Simulated SRs have relatively constant triaxiality profiles. Simulated stellar halos show a large variety of structural properties, with quantitative but no clear qualitative differences between FRs and SRs. At the same stellar mass, stellar halo properties show a gradual transition and significant overlap between the two families, despite the clear bimodality in the central regions. This is in agreement with observations of extended photometry and kinematics. [abridged], accepted for publication in A&A, 25 pages, 22 figures

Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M_* >10e11M_sun) than in Milky Way like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ~5-10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical timescale for star-formation to be quenched ~2 Gyr in groups but more than ~5 Gyr for satellites around Milky Way like primaries. The gas contents we infer are consistent with z=0 observations of HI gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time., 5 pages, 3 figures. Accepted for publication in MNRAS Letters

Publisher's version (útgefin grein), We use the IllustrisTNG (TNG) cosmological simulations to provide theoretical expectations for the dark matter mass fractions (DMFs) and circular velocity profiles of galaxies. TNG predicts flat circular velocity curves for z = 0 Milky Way (MW)-like galaxies beyond a few kpc from the galaxy centre, in better agreement with observational constraints than its predecessor, Illustris. TNG also predicts an enhancement of the dark matter mass within the 3D stellar half-mass radius (⁠rhalf⁠; M200c=1010−1013M⊙⁠, z ≤ 2) compared to its dark matter only and Illustris counterparts. This enhancement leads TNG present-day galaxies to be dominated by dark matter within their inner regions, with fDM(, MRL is supported by a COFUND/Durham Junior Research Fellowship under EU grant 609412. MRL acknowledges support by a Grant of Excellence from the Icelandic Research Fund (grant number 173929-051). VS, RP, and RW acknowledge support through the European Research Council under ERC-StG grant EXAGAL-308037, and would like to thank the Klaus Tschira Foundation. The Flatiron Institute is supported by the Simons Foundation.

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., 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 have finished two new, extremely large hydrodynamical simulations of galaxy formation that significantly advance the state of the art in cosmology. Together with accompanying dark matter only runs, we call them ‘IllustrisTNG’, the next generation Illustris simulations. Our largest and most ambitious calculation follows a cosmological volume 300 megaparsecs on a side and self-consistently solves the equations of magnetohydrodynamics and self-gravity coupled to the fundamental physical processes driving galaxy formation. We have employed AREPO, a sophisticated moving-mesh code developed by our team over the past 7 years and equipped with an improved, multi-purpose galaxy formation physics model. The simulated universe contains tens of thousands of galaxies encompassing a variety of environments, mass scales and evolutionary stages. The groundbreaking volume of TNG enables us to sample statistically significant sets of rare astrophysical objects like rich galaxy clusters, and to study galaxy formation and the spatial clustering of matter over a very large range of spatial scales. Here we report some early results on the matter and galaxy clustering found in the simulations. The two-point galaxy correlation function of our largest simulation agrees extremely well with the best available observational constraints from the Sloan Digital Sky Survey, both as a function of galaxy stellar mass and color. The predicted impact of baryonic physics on the matter power spectrum is sizeable and needs to be taken into account in precision studies of cosmology. Interestingly, this impact appears to be fairly robust to the details of the modelling of supermassive black holes, provided this reproduces the scaling properties of the intracluster medium of galaxy clusters.

The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the Arepo code and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} \leq M_{\rm 200c}/M_{\rm sun} \leq 10^{15}$) at recent times ($z \leq 1$). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light. The exact ICL fraction depends sensitively on the definition of a central galaxy's mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5\times 10^{14}M_{\rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matter's ( $-3.5 < ��_{\rm 3D} < -3$). Total halo mass is a very good predictor of stellar mass, and vice versa: at $z=0$, the 3D stellar mass measured within 30 kpc scales as $\propto (M_{\rm 500c})^{0.49}$ with a $\sim 0.12$ dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of TNG less massive galaxies ($< 10^{11}M_{\rm sun}$ in stars) is in good agreement with the measured galaxy stellar mass functions at recent epochs. The 3D sizes of massive galaxies fall too on a tight ($\sim$0.16 dex scatter) power-law relation with halo mass, with $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 500c})^{0.53}$. Even more fundamentally, halo mass alone is a good predictor for the whole stellar mass profiles beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo., Accepted by MNRAS, updated to match published version. Highlights: Figures 5, 9, 11. The IllustrisTNG website can be found at http://www.tng-project.org/

The distribution of elements in galaxies provides a wealth of information about their production sites and their subsequent mixing into the interstellar medium. Here we investigate the distribution of elements within stars in the IllustrisTNG simulations. In particular, we analyze the abundance ratios of magnesium and europium in Milky Way-like galaxies from the TNG100 simulation (stellar masses ${\log} (M_\star / {\rm M}_\odot) \sim 9.7 - 11.2$). As abundances of magnesium and europium for individual stars in the Milky Way are observed across a variety of spatial locations and metallicities, comparison with the stellar abundances in our more than $850$ Milky Way-like galaxies provides stringent constraints on our chemical evolutionary methods. To this end we use the magnesium to iron ratio as a proxy for the effects of our SNII and SNIa metal return prescription, and a means to compare our simulated abundances to a wide variety of galactic observations. The europium to iron ratio tracks the rare ejecta from neutron star -- neutron star mergers, the assumed primary site of europium production in our models, which in turn is a sensitive probe of the effects of metal diffusion within the gas in our simulations. We find that europium abundances in Milky Way-like galaxies show no correlation with assembly history, present day galactic properties, and average galactic stellar population age. In general, we reproduce the europium to iron spread at low metallicities observed in the Milky Way, with the level of enhancement being sensitive to gas properties during redshifts $z \approx 2-4$. We show that while the overall normalization of [Eu/Fe] is susceptible to resolution and post-processing assumptions, the relatively large spread of [Eu/Fe] at low [Fe/H] when compared to that at high [Fe/H] is very robust., Comment: 18 pages, 14 figures, accepted to MNRAS

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

We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multi-mode feedback. In this paper we give a comprehensive description of the physical and numerical advances which form the core of the IllustrisTNG (The Next Generation) framework. We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content, and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black hole driven kinetic feedback at low accretion rates, magnetohydrodynamics, and improvements to the numerical scheme. Using a suite of (25 Mpc $h^{-1}$)$^3$ cosmological boxes we assess the outcome of the new model at our fiducial resolution. The presence of a self-consistently amplified magnetic field is shown to have an important impact on the stellar content of $10^{12} M_{\rm sun}$ haloes and above. Finally, we demonstrate that the new galactic winds promise to solve key problems identified in Illustris in matching observational constraints and affecting the stellar content and sizes of the low mass end of the galaxy population., 31 pages, 16 figures; published in MNRAS (Figs. 6 and 7 revised + additional minor changes to match accepted version)

We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about $10^{8.5}\,{\rm M}_\odot$, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass-stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function., 16 pages, 11 figures, submitted to MNRAS, the IllustrisTNG project website: www.tng-project.org, comments welcome

We study the chaotic-like behavior of cosmological simulations by quantifying how minute perturbations grow over time and manifest as macroscopic differences in galaxy properties. When we run pairs of 'shadow' simulations that are identical except for random minute initial displacements to particle positions (e.g. of order 1e-7pc), the results diverge from each other at the individual galaxy level (while the statistical properties of the ensemble of galaxies are unchanged). After cosmological times, the global properties of pairs of 'shadow' galaxies that are matched between the simulations differ from each other generally at a level of ~2-25%, depending on the considered physical quantity. We perform these experiments using cosmological volumes of (25-50Mpc/h)^3 evolved either purely with dark matter, or with baryons and star-formation but no feedback, or using the full feedback model of the IllustrisTNG project. The runs cover four resolution levels spanning a factor of 512 in mass. We find that without feedback the differences between shadow galaxies generally become smaller as the resolution increases, but with the IllustrisTNG model the results are mostly converging towards a 'floor'. This hints at the role of feedback in setting the chaotic properties of galaxy formation. Importantly, we compare the macroscopic differences between shadow galaxies to the overall scatter in various galaxy scaling relations, and conclude that for the star formation-mass and the Tully-Fisher relations the butterfly effect in our simulations contributes significantly to the overall scatter. We find that our results are robust to whether random numbers are used in the sub-grid models or not. We discuss the implications for galaxy formation theory in general and for cosmological simulations in particular., Comment: Key figures: 10 & 11. Accepted for publication in ApJ. This final version includes a new verification of the conclusions in simulations that completely avoid the usage of random numbers

In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We include these non-ideal terms for two MHD techniques: the Powell 8-wave formalism and a constrained transport scheme, which evolves the cell-centred magnetic vector potential. We test our implementation against problems of increasing complexity, such as one- and two-dimensional diffusion problems, and the evolution of progressive and stationary Alfv\'en waves. On these test problems, our implementation recovers the analytic solutions to second-order accuracy. As first applications, we investigate the tearing instability in magnetized plasmas and the gravitational collapse of a rotating magnetized gas cloud. In both systems, resistivity plays a key role. In the former case, it allows for the development of the tearing instability through reconnection of the magnetic field lines. In the latter, the adopted (constant) value of ohmic resistivity has an impact on both the gas distribution around the emerging protostar and the mass loading of magnetically driven outflows. Our new non-ideal MHD implementation opens up the possibility to study magneto-hydrodynamical systems on a moving mesh beyond the ideal MHD approximation., Comment: 18 pages, 11 figures, accepted for publication in MNRAS. Revisions to match the accepted version

The distribution of metals in the intra-cluster medium encodes important information about the enrichment history and formation of galaxy clusters. Here we explore the metal content of clusters in IllustrisTNG - a new suite of galaxy formation simulations building on the Illustris project. Our cluster sample contains 20 objects in TNG100 - a ~(100 Mpc)^3 volume simulation with 2x1820^3 resolution elements, and 370 objects in TNG300 - a ~(300 Mpc)^3 volume simulation with 2x2500^3 resolution elements. The z=0 metallicity profiles agree with observations, and the enrichment history is consistent with observational data going beyond z~1, showing nearly no metallicity evolution. The abundance profiles vary only minimally within the cluster samples, especially in the outskirts with a relative scatter of ~15%. The average metallicity profile flattens towards the center, where we find a logarithmic slope of -0.1 compared to -0.5 in the outskirts. Cool core clusters have more centrally peaked metallicity profiles (~0.8 solar) compared to non-cool core systems (~0.5 solar), similar to observational trends. Si/Fe and O/Fe radial profiles follow positive gradients. The outer abundance profiles do not evolve below z~2, whereas the inner profiles flatten towards z=0. More than ~80% of the metals in the intra-cluster medium have been accreted from the proto-cluster environment, which has been enriched to ~0.1 solar already at z~2. We conclude that the intra-cluster metal distribution is uniform among our cluster sample, nearly time-invariant in the outskirts for more than 10 Gyr, and forms through a universal enrichment history., Comment: 21 pages, 16 figures. MNRAS accepted. The IllustrisTNG project website can be found at http://www.tng-project.org/

The fundamental metallicity relation (FMR) is a postulated correlation between galaxy stellar mass, star formation rate (SFR), and gas-phase metallicity. At its core, this relation posits that offsets from the mass-metallicity relation (MZR) at a fixed stellar mass are correlated with galactic SFR. In this Letter, we use hydrodynamical simulations to quantify the time-scales over which populations of galaxies oscillate about the average SFR and metallicity values at fixed stellarmass.We find that Illustris and IllustrisTNG predict that galaxy offsets from the star formation main sequence and MZR oscillate over similar time-scales, are often anticorrelated in their evolution, evolve with the halo dynamical time, and produce a pronounced FMR. Our models indicate that galaxies oscillate about equilibrium SFR and metallicity values - set by the galaxy's stellar mass - and that SFR and metallicity offsets evolve in an anticorrelated fashion. This anticorrelated variability of the metallicity and SFR offsets drives the existence of the FMR in our models. In contrast to Illustris and IllustrisTNG, we speculate that the SFR and metallicity evolution tracks may become decoupled in galaxy formation models dominated by feedback-driven globally bursty SFR histories, which could weaken the FMR residual correlation strength. This opens the possibility of discriminating between bursty and non-bursty feedback models based on the strength and persistence of the FMR - especially at high redshift., United States. National Aeronautics and Space Administration (HST-HF2-51341.001-A), United States. Department of Energy. Computational Science Graduate Fellowship Program (Grant DEFG02-97ER25308)

Using the astrometry from the ESA's Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significant components with high radial velocity anisotropy, beta > 0.6. However, only in one third of the hosts do the high-beta stars contribute significantly to the accreted stellar halo overall, similar to what is observed in the Milky Way. For this particular subset we reveal the origin of the dominant stellar halo component with high metallicity, [Fe/H]~-1, and high orbital anisotropy, beta>0.8, by tracing their stars back to the epoch of accretion. It appears that, typically, these stars come from a single dwarf galaxy with a stellar mass of order of 10^9-10^10 Msol that merged around 6-10 Gyr ago, causing a sharp increase in the halo mass. Our study therefore establishes a firm link between the excess of radially anisotropic stellar debris in the Milky Way halo and an ancient head-on collision between the young Milky Way and a massive dwarf galaxy, replaced with the accepted version (very minor changes)

We use the IllustrisTNG simulations to investigate the evolution of the mass-metallicity relation (MZR) for star-forming cluster galaxies as a function of the formation history of their cluster host. The simulations predict an enhancement in the gas-phase metallicities of star-forming cluster galaxies (10^9< M_star, 5 pages, 4 figures, accepted for publication in MNRAS Letters

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

Feedback processes by active galactic nuclei in the centres of galaxy clusters appear to prevent large-scale cooling flows and impede star formation. However, the detailed heating mechanism remains uncertain. One promising heating scenario invokes the dissipation of Alfv\'en waves that are generated by streaming cosmic rays (CRs). In order to study this idea, we use three-dimensional magneto-hydrodynamical simulations with the AREPO code that follow the evolution of jet-inflated bubbles that are filled with CRs in a turbulent cluster atmosphere. We find that a single injection event produces the CR distribution and heating rate required for a successful CR heating model. As a bubble rises buoyantly, cluster magnetic fields drape around the leading interface and are amplified to strengths that balance the ram pressure. Together with helical magnetic fields in the bubble, this initially confines the CRs and suppresses the formation of interface instabilities. But as the bubble continues to rise, bubble-scale eddies significantly amplify radial magnetic filaments in its wake and enable CR transport from the bubble to the cooling intracluster medium. By varying the jet parameters, we obtain a rich and diverse set of jet and bubble morphologies ranging from Fanaroff-Riley type I-like (FRI) to FRII-like jets. We identify jet energy as the leading order parameter (keeping the ambient density profiles fixed), whereas jet luminosity is primarily responsible for setting the Mach numbers of shocks around FRII-like sources. Our simulations also produce FRI-like jets that inflate bubbles without detectable shocks and show morphologies consistent with cluster observations., Comment: 25 pages, 20 figures, accepted by MNRAS

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

A method for improving the accuracy of hydrodynamical codes that use a moving Voronoi mesh is described. Our scheme is based on a new regularization scheme that constrains the mesh to be centroidal to high precision while still allowing the cells to move approximately with the local fluid velocity, thereby retaining the quasi-Lagrangian nature of the approach. Our regularization technique significantly reduces mesh noise that is attributed to changes in mesh topology and deviations from mesh regularity. We demonstrate the advantages of our method on various test problems, and note in particular improvements obtained in handling shear instabilities, mixing, and in angular momentum conservation. Calculations of adiabatic jets in which shear excites Kelvin Helmholtz instability show reduction of mesh noise and entropy generation. In contrast, simulations of the collapse and formation of an isolated disc galaxy are nearly unaffected, showing that numerical errors due to the choice of regularization do not impact the outcome in this case., 9 pages, 14 figures, MNRAS submitted

We simulate the phase-space distribution of stellar mass in nine massive Λ cold dark matter galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 ≈ 7.5–33 × 1014 M⊙). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ≳90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ≳10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a ‘cD envelope'. These clouds typically originate from tidal stripping at z ≲ 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc2 of a z ∼ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80–95 per cent of stellar mass below a SB of μV ∼ 26.5 mag arcsec−2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ∼30 per cent.