Search

Dwarf galaxies are thought to quench primarily due to environmental processes most typically occurring in galaxy groups and clusters or around single, massive galaxies. However, at earlier epochs, ($5 < z < 2$), the collapse of large scale structure (forming Zel'dovich sheets and subsequently filaments of the cosmic web) can produce volume-filling accretion shocks which elevate large swaths of the intergalactic medium (IGM) in these structures to a hot ($T>10^6$ K) phase. We study the impact of such an event on the evolution of central dwarf galaxies ($5.5 < \log M_* < 8.5$) in the field using a spatially large, high resolution cosmological zoom simulation which covers the cosmic web environment between two protoclusters. We find that the shock-heated sheet acts as an environmental quencher much like clusters and filaments at lower redshift, creating a population of quenched, central dwarf galaxies. Even massive dwarfs which do not quench are affected by the shock, with reductions to their sSFR and gas accretion. This process can potentially explain the presence of isolated quenched dwarf galaxies, and represents an avenue of pre-processing, via which quenched satellites of bound systems quench before infall., 15 pages, 10 figures. Submitted to MNRAS

Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{\mathrm{eff}}$ (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass $M_{\ast}\geqslant 10^{11} \mathrm{M_{\odot}}$ in TNG100-1. The mock observations reproduce the key outer ($1.5R_{\mathrm{eff}}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $\gamma_{\mathrm{outer}}$ and the kurtosis $h_{4}$'s gradient. We find that $h_{4}$ is uncorrelated with stellar orbital anisotropy beyond $R_{\mathrm{eff}}$; instead we find that the variations in $\gamma_{\mathrm{outer}}$ and outer $h_{4}$ (a good proxy for $h_{4}$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities., Comment: Accepted for publication in MNRAS. 19 pages, 11 figures. Moderate revision since v1

The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies and thus provide a powerful approach to predict and understand the IA signal. Here we employ the novel, large-volume hydrodynamical simulation MTNG740, a product of the MillenniumTNG (MTNG) project, to study the IA of galaxies. We measure the projected correlation functions between the intrinsic shape/shear of galaxies and various tracers of large-scale structure, $w_{+g},\ w_{+m},\ w_{++}$ over the radial range $r_{\rm p} \in [0.02 , 200]\,h^{-1}{\rm Mpc}$ and at redshifts $z=0.0$, $0.5$ and $1.0$. We detect significant signal-to-noise IA signals with the density field for both elliptical and spiral galaxies. We also find significant intrinsic shear-shear correlations for ellipticals. We further examine correlations of the intrinsic shape of galaxies with the local tidal field. Here we find a significant IA signal for elliptical galaxies assuming a linear model. We also detect a weak IA signal for spiral galaxies under a quadratic tidal torquing model. Lastly, we measure the alignment between central galaxies and their host dark-matter halos, finding small to moderate misalignments between their principal axes that decline with halo mass., Comment: 16 pages, 14 figures

The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 Mpc/h of the Milky Way. Our analysis relies on constrained realization simulations and forward modelling of the 2M++ galaxy catalogue. We find that the angular distribution of neutrinos is anti-correlated with the projected matter density, due to the capture and deflection of neutrinos by massive structures along the line of sight. Of relevance to tritium capture experiments, we find that the gravitational clustering effect of the large-scale structure on the local number density of neutrinos is more important than that of the Milky Way for neutrino masses less than 0.1 eV. Nevertheless, we predict that the density of relic neutrinos is close to the cosmic average, with a suppression or enhancement over the mean of (-0.3%, +7%, +27%) for masses of (0.01, 0.05, 0.1) eV. This implies no more than a marginal increase in the event rate for tritium capture experiments like PTOLEMY. We also predict that the CNB and CMB rest frames coincide for 0.01 eV neutrinos, but that neutrino velocities are significantly perturbed for masses larger than 0.05 eV. Regardless of mass, we find that the angle between the neutrino dipole and the ecliptic plane is small, implying a near-maximal annual modulation in the bulk velocity. Along with this paper, we publicly release our simulation data, comprising more than 100 simulations for six different neutrino masses., Comment: 26 pages, 8 figures, submitted to JCAP

We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics hydrodynamical simulations with corresponding dark-matter-only runs we quantify the impact of baryonic physics on the most important weak lensing statistics. Likewise, we predict the impact of massive neutrinos reliably far into the non-linear regime. We also demonstrate that the "fixed & paired" variance suppression technique increases the statistical robustness of the simulation predictions on large scales not only for time slices but also for continuously output lightcone data. We find that both baryonic and neutrino effects substantially impact weak lensing shear measurements, with the latter dominating over the former on large angular scales. Thus, both effects must explicitly be included to obtain sufficiently accurate predictions for stage IV lensing surveys. Reassuringly, our results agree accurately with other simulation results where available, supporting the promise of simulation modelling for precision cosmology far into the non-linear regime., 14 pages, 10 figures, comments welcome

We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $\Lambda$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stellar disky morphology, with a stellar mass in the range of $M_* = 10^{10.5 - 11.2}$ Msun, and within a MW-like Mpc-scale environment at z=0. These are resolved with baryonic (dark matter) mass resolution of $8.5\times10^4$ Msun ($4.5\times10^5$ Msun) and $\sim150$ pc of average spatial resolution in the star-forming regions: we therefore expand by many factors (2 orders of magnitude) the sample size of cosmologically-simulated analogs with similar ($\times 10$ better) numerical resolution. The majority of TNG50 MW/M31 analogs at $z=0$ exhibit a bar, 60 per cent are star-forming, the sample includes 3 Local Group (LG)-like systems, and a number of galaxies host one or more satellites as massive as e.g. the Magellanic Clouds. Even within such a relatively narrow selection, TNG50 reveals a great diversity in galaxy and halo properties, as well as in past histories. Within the TNG50 sample, it is possible to identify several simulated galaxies whose integral and structural properties are consistent, one or more at a time, with those measured for the Galaxy and Andromeda. With this paper, we document and release a series of broadly applicable data products that build upon the IllustrisTNG public release and aim to facilitate easy access and analysis by public users. These include datacubes across snapshots ($0 \le z \le 7$) for each TNG50 MW/M31-like galaxy, and a series of value-added catalogs that will be continually expanded to provide a convenient and up to date community resource., Comment: Submitted to MNRAS. Feedback from the community encouraged. This is also a data-release paper: see visuals and documentation at https://www.tng-project.org/data/milkyway+andromeda/. See other TNG50 MW/M31-based papers also on astro-ph today: Ramesh+ and Sotillo-Ramos+

Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magneto-hydrodynamical simulations of entire galaxies with masses ranging from $10^{10}$ to $10^{12}\,\mathrm{M}_\odot$ and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from $0.1\,\mathrm{GeV}~c^{-1}$ to $100\,\mathrm{TeV}~c^{-1}$. We find that hadronic cooling of CRs dominates over Alfv\'{e}n cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of four in dwarf galaxies. All simulated low-mass halos ($M=10^{10}$, $10^{11}$, and $3\times10^{11}\,\mathrm{M}_\odot$) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with $M=10^{12}\,\mathrm{M}_\odot$ does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of $\langle{D}\rangle_{e_\mathrm{cr}}\sim10^{28}\,\mathrm{cm^2\,s^{-1}}$ in the disc up to $3\times10^{29}\,\mathrm{cm^2\,s^{-1}}$ in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models., Comment: 20 pages, submitted to MNRAS, comments welcome

In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, $\tau$, and the Compton-y parameter, $Y$, receive large line-of-sight contributions which can be removed effectively by applying a Compensated Aperture Photometry (CAP) filter. In contrast with other $\tau$ probes (e.g., X-rays and Fast Radio Bursts), the kSZ-inferred $\tau$ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line-of-sight. Additionally, we perform fits to the $Y-M$ and $\tau-M$ scaling relations and report best-fit parameters adopting the smoothly broken power law (SBPL) formalism. We note that subgrid physics modeling can broaden the error bar on these by 30\% for intermediate-mass halos ($\sim$$10^{13} \, {\rm M}_{\odot}$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration, and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than halos in real observations, which can bias the inferred SZ profiles by $\sim$20\% for $L_\ast$-galaxies., Comment: 14 pages, 6 figures

We analyze the internal cross-section structure and dynamics of inter-halo cosmic-web filaments that connect to massive high-$z$ galaxies. Our analysis is based on a high-resolution \texttt{AREPO} cosmological simulation zooming in on a volume encompassing three filaments of $\sim 1 \, \Mpc$ length feeding three massive haloes of $\sim 10^{12}\,\msun$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a strong cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \kpc$, within which the gas is in virial equilibrium in the potential well set by an isothermal dark-matter filament. The cross-section of the filament can be divided into three radial zones. In the outer ``thermal'' (\textbf{T}) zone, $r \geq 0.65 \, r_{\rm shock}$, the inward gravity and ram-pressure forces from the inflowing gas are over-balanced by outward thermal pressure forces, causing the inflowing gas to decelerate and the shock to expand outwards. In the intermediate ``vortex'' (\textbf{V}) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outward force is provided mainly by centrifugal forces associated with these vortices and, to a lesser extent, by global rotation and thermal pressure, while shear and turbulent forces that are also associated with the vortex structure add to the force inward. The inner ``stream'' (\textbf{S}) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core defining the cold streams that feed the galaxies, with $T\sim 3 \times 10^4 \, \Kdegree$ and $n_{\rm H}\sim 0.01 \cmc$. The core is formed by an isobaric cooling flow that results from post-shock cooling on a free-fall time, and is associated with a decrease in the outward forces, though the gas within the core exhibits both inflows and outflows. [abridged], 30 pages, 20 figures, submitted to MNRAS

We investigate the properties of mergers that are comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. It was previously shown that these mergers occur over a wide range of times ($6-10\,$Gyr ago). We find that the merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^810^6\,\rm{M}_{\odot}$), but most of these do not follow the merger to low orbital energies. Between $0-1$ of these satellites may survive to $z=0$, but with no clear signatures of their past association. We show that the fraction of stars originating from GSE-like mergers is reduced for lower metallicities (reaching a minimum around $\text{[Fe/H]} = -2$), and also within $5\,$kpc of the galactic centre. Whilst these central regions are dominated by $\textit{in-situ}$ stars, the $\textit{ex-situ}$ fraction trends towards a 100 per cent asymptote when considering the most metal-poor stars ($\text{[Fe/H]}\ll-2.5$). Considering this, its near proximity, and its small volume on the sky, the Galactic centre lends itself as a prime environment in the search for the stars from the earliest galaxies, whilst avoiding contamination from GSE stars., Comment: 20 pages, 14 figures, submitted to MNRAS

Magnetic fields are one of most concealed components of the universe. They are observed as part of the intergalactic medium and on galaxy cluster scales, however their origin and evolution is unclear. In this work we use the IllustrisTNG simulation to investigate the effects of magnetic fields in cosmic voids, the least dense regions of the universe. We find that, under the hypothesis of a uniform primordial magnetic field, the voids still reflect the primordial properties of the fields. On the other hand, the galaxies in their interior acquire weaker magnetic fields than galaxies in denser environments., Comment: 3 pages, 2 figures, accepted in Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia

We present the new public version of the KETJU supermassive black hole (SMBH) dynamics module, as implemented into GADGET-4. KETJU adds a small region around each SMBH where the dynamics of the SMBHs and stellar particles are integrated using an algorithmically regularised integrator instead of the leapfrog integrator with gravitational softening used by GADGET-4. This enables modelling SMBHs as point particles even during close interactions with stellar particles or other SMBHs, effectively removing the spatial resolution limitation caused by gravitational softening. KETJU also includes post-Newtonian corrections, which allows following the dynamics of SMBH binaries to sub-parsec scales and down to tens of Schwarzschild radii. Systems with multiple SMBHs are also supported, with the code also including the leading non-linear cross terms that appear in the post-Newtonian equations for such systems. We present tests of the code showing that it correctly captures, at sufficient mass resolution, the sinking driven by dynamical friction and binary hardening driven by stellar scattering. We also present an example application demonstrating how the code can be applied to study the dynamics of SMBHs in mergers of multiple galaxies and the effect they have on the properties of the surrounding galaxy. We expect that the presented KETJU SMBH dynamics module can also be straightforwardly incorporated into other codes similar to GADGET-4, which would allow coupling small-scale SMBH dynamics to the rich variety of galactic physics models that exist in the literature., Comment: 21 pages, 19 figures. Code available from https://www.mv.helsinki.fi/home/phjohans/ketju

Observable thermodynamical properties of the intracluster medium (ICM) reflect the complex interplay between AGN feedback and the gravitational collapse of haloes. Using the large volume TNG300 simulation of the IllustrisTNG project we provide predictions for X-ray and Sunyaev-Zel'dovich (SZ) scaling relations for a sample of over 30,000 haloes that cover a wide mass range from galaxies to massive galaxy clusters ($M_{\rm 500crit}$ $\in [10^{12}$ M$_{\odot} - 2\times 10^{15}$ M$_{\odot}$]). We produce mock X-ray observations of simulated haloes using methods that are consistent with observational techniques. Thus, we investigate the scaling relations between the soft-band X-ray luminosity, spectroscopic temperature, gas mass fraction, $Y_{\rm X}$ and $Y_{\rm SZ}$ as a function of halo mass, and we find broad agreement between IllustrisTNG and the observed relations. Our results highlight the scatter and bias introduced by estimated masses, and thus the importance of converting simulated ICM properties to the observable space when comparing simulations to current X-ray observations. The wide range of halo masses in our sample provides new insights into the shape of the X-ray and SZ scaling relations across three orders of magnitude in mass. Our findings show strong evidence for a break in $z=0$ scaling relations. We introduce a smoothly broken power law model which robustly captures the location of this break, the width of the transition region around the break, as well as the slope dependence on halo mass. Our results inform the next generation of subgrid black hole feedback models and provide predictions for ongoing and future observational surveys., Comment: 29 pages, 15 figures, 7 tables, 2 appendices; submitted to MNRAS. This is a companion paper to Pop et al. 2022: "Unifying Sunyaev-Zel'dovich and X-ray predictions from clusters to galaxy groups: the impact of X-ray mass estimates on the Y-M scaling relation"

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_*, 15 pages, 11 figures, accepted by MNRAS, comments welcome

Early galaxies were the radiation source for reionization, with the photoheating feedback from the reionization process expected to reduce the efficiency of star formation in low mass haloes. Hence, to fully understand reionization and galaxy formation, we must study their impact on each other. The THESAN project has so far aimed to study the impact of galaxy formation physics on reionization, but here we present the new THESAN simulations with a factor 50 higher resolution ($m_{\rm b} \approx 10^4$~M$_\odot$) that aim to self-consistently study the back-reaction of reionization on galaxies. By resolving haloes with virial temperatures $T_{\rm vir} < 10^4$~K, we are able to demonstrate that simplistic, spatially-uniform, reionization models are not sufficient to study early galaxy evolution. Comparing the self-consistent THESAN model (employing fully coupled radiation hydrodynamics) to a uniform UV background, we are able to show that galaxies in THESAN are predicted to be larger in physical extent (by a factor $\sim 2$), less metal enriched (by $\sim 0.2$~dex), and less abundant (by a factor $\sim 10$ at $M_{\rm 1500}~=~-10$) by $z=5$. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation $z \gtrsim 3$ must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties., Submitted to MNRAS

Luminous red galaxies (LRGs) and blue star-forming emission-line galaxies (ELGs) are key tracers of large-scale structure used by cosmological surveys. Theoretical predictions for such data are often done via simplistic models for the galaxy-halo connection. In this work, we use the large, high-fidelity hydrodynamical simulation of the MillenniumTNG project (MTNG) to inform a new phenomenological approach for obtaining an accurate and flexible galaxy-halo model on small scales. Our aim is to study LRGs and ELGs at two distinct epochs, $z = 1$ and $z = 0$, and recover their clustering down to very small scales, $r \sim 0.1 \ {\rm Mpc}/h$, i.e. the one-halo regime, while a companion paper extends this to a two-halo model for larger distances. The occupation statistics of ELGs in MTNG inform us that: (1) the satellite occupations exhibit a slightly super-Poisson distribution, contrary to commonly made assumptions, and (2) that haloes containing at least one ELG satellite are twice as likely to host a central ELG. We propose simple recipes for modeling these effects, each of which calls for the addition of a single free parameter to simpler halo occupation models. To construct a reliable satellite population model, we explore the LRG and ELG satellite radial and velocity distributions and compare them with those of subhalos and particles in the simulation. We find that ELGs are anisotropically distributed within halos, which together with our occupation results provides strong evidence for cooperative galaxy formation (manifesting itself as one-halo galaxy conformity); i.e.~galaxies with similar properties form in close proximity to each other. Our refined galaxy-halo model represents a useful improvement of commonly used analysis tools and thus can be of help to increase the constraining power of large-scale structure surveys., Comment: 14 pages, 10 figures; submitted to MNRAS

Approximate methods to populate dark matter halos with galaxies are of great utility to large galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, $z=1$ and $z=0$, in the large volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to $r\sim 0.1{\rm Mpc}/h$, while here we focus on modeling galaxy assembly bias in the two-halo regime, $r\gtrsim 1{\rm Mpc}/h$. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales ($r\sim 20{\rm Mpc}/h$), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10-15\% discrepancy present in the standard halo model prescription is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties ("shear" being the best-performing one) rather than intrinsic ones (e.g., concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and "shear") is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher-order statistics such as nearest-neighbor counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters., Comment: 16 pages, 8 figures; submitted to MNRAS

One of the main limitations in precision cluster cosmology arises from systematic errors and uncertainties in estimating cluster masses. Using the Mock-X pipeline, we produce synthetic X-ray images and derive cluster and galaxy group X-ray properties for a sample of over 30,000 simulated galaxy groups and clusters with $M_{\rm 500crit}$ between $10^{12}$ and $2\times 10^{15}$ M$_{\odot}$ in IllustrisTNG. We explore the similarities and differences between IllustrisTNG predictions of the Sunyaev-Zel'dovich and X-ray scaling relations with mass. We find a median hydrostatic mass bias $b = 0.125 \pm 0.003$ for $M_{\rm 500crit}$ $>10^{13}$ M$_{\odot}$. The bias increases to $b = 0.17 \pm 0.004$ when masses are derived from synthetic X-ray observations. We model how different underlying assumptions about the dependence of $Y_{\rm X}$ on halo mass can generate biases in the observed $Y_{\rm SZ} - M_{Y_{\rm X}}$ scaling relation. In particular, the simplifying assumption that $Y_{\rm X} - M_{\rm tot}$ is self-similar at all mass scales largely hides the break in $Y_{\rm SZ} - M_{\rm tot}$ and overestimates $Y_{\rm SZ}$ at galaxy and groups scales. We show that calibrating the $Y_{\rm X}-$mass proxy using a new model for a smoothly broken power law reproduces the true underlying $Y_{\rm SZ} - M_{\rm tot}$ scaling relation with high accuracy. Moreover, $M_{Y_{\rm X}}$ estimates calibrated with this method lead to $Y_{\rm SZ} - M_{Y_{\rm X}}$ predictions that are not biased by the presence of lower mass clusters or galaxy groups in the sample. Finally, we show that our smoothly broken power law model provides a robust way to derive the $Y_{\rm X}-$mass proxy, significantly reducing the level of mass bias for clusters, groups, and galaxies., 12 pages, 5 figures; submitted to MNRAS. This is a companion paper to Pop et al. 2022: "Sunyaev-Zel'dovich effect and X-ray scaling relations of galaxies, groups and clusters in the IllustrisTNG simulations"

Modern redshift surveys are tasked with mapping out the galaxy distribution over enormous distance scales. Existing hydrodynamical simulations, however, do not reach the volumes needed to match upcoming surveys. We present results for the clustering of galaxies using a new, large volume hydrodynamical simulation as part of the MillenniumTNG (MTNG) project. With a computational volume that is $\approx15$ times larger than the next largest such simulation currently available, we show that MTNG is able to accurately reproduce the observed clustering of galaxies as a function of stellar mass. When separated by colour, there are some discrepancies with respect to the observed population, which can be attributed to the quenching of satellite galaxies in our model. We combine MTNG galaxies with those generated using a semi-analytic model to emulate the sample selection of luminous red galaxies (LRGs) and emission line galaxies (ELGs), and show that although the bias of these populations is approximately (but not exactly) constant on scales larger than $\approx10$ Mpc, there is significant scale-dependent bias on smaller scales. The amplitude of this effect varies between the two galaxy types, and also between the semi-analytic model and MTNG. We show that this is related to the distribution of haloes hosting LRGs and ELGs. Using mock SDSS-like catalogues generated on MTNG lightcones, we demonstrate the existence of prominent baryonic acoustic features in the large-scale galaxy clustering. We also demonstrate the presence of realistic redshift space distortions in our mocks, finding excellent agreement with the multipoles of the redshift-space clustering measured in SDSS data., Comment: 15 pages, 11 figures. Accepted for publication in MNRAS

The early release science results from $\textit{JWST}$ have yielded an unexpected abundance of high-redshift luminous galaxies that seems to be in tension with current theories of galaxy formation. However, it is currently difficult to draw definitive conclusions form these results as the sources have not yet been spectroscopically confirmed. It is in any case important to establish baseline predictions from current state-of-the-art galaxy formation models that can be compared and contrasted with these new measurements. In this work, we use the new large-volume ($L_\mathrm{box}\sim 740 \, \mathrm{cMpc}$) hydrodynamic simulation of the MillenniumTNG project, suitably scaled to match results from higher resolution - smaller volume simulations, to make predictions for the high-redshift ($z\gtrsim8$) galaxy population and compare them to recent $\textit{JWST}$ observations. We show that the simulated galaxy population is broadly consistent with observations until $z\sim10$. From $z\approx10-12$, the observations indicate a preference for a galaxy population that is largely dust-free, but is still consistent with the simulations. Beyond $z\gtrsim12$, however, our simulation results underpredict the abundance of luminous galaxies and their star-formation rates by almost an order of magnitude. This indicates either an incomplete understanding of the new $\textit{JWST}$ data or a need for more sophisticated galaxy formation models that account for additional physical processes such as Population~III stars, variable stellar initial mass functions, or even deviations from the standard $\Lambda$CDM model. We emphasise that any new process invoked to explain this tension should only significantly influence the galaxy population beyond $z\gtrsim10$, while leaving the successful galaxy formation predictions of the fiducial model intact below this redshift., Comment: Accepted for publication in MNRAS -- Part of the initial set of papers introducing the MillenniumTNG project. Visit www.mtng-project.org for more details

Early-type galaxies (ETGs) possess total radial density profiles that are largely described by singular isothermal spheres, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{eff}$ that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on the formation of these non-canonical kinematic features through mock ETGs with well-studied density profiles in IllustrisTNG. We create 2D maps mimicking integral-field-unit (IFU) observations to extract projected stellar kinematic features for 207 ETGs in the TNG100-1 box with stellar mass $M_{\ast} > 10^{11} M_{\odot}$. The mock observations reproduce the key outer ($1.5R_{eff}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $γ_{outer}$ and the kurtosis $h_4$'s gradient. We find that $h_4$ is uncorrelated with stellar orbital anisotropy beyond $R_{eff}$; instead we find that the variations in $γ_{outer}$ and outer $h_4$ (which dominates the $h_4$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the dominant driver of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates these two properties., Submitted to MNRAS. 19 pages, 11 figures

The pattern speed with which galactic bars rotate is intimately linked to the amount of dark matter in the inner regions of their host galaxies. In particular, dark matter haloes act to slow down bars via torques exerted through dynamical friction. Observational studies of barred galaxies tend to find that bars rotate fast, while hydrodynamical cosmological simulations of galaxy formation and evolution in the $��$CDM framework have previously found that bars slow down excessively. This has led to a growing tension between fast bars and the $��$CDM cosmological paradigm. In this study we revisit this issue, using the Auriga suite of high resolution, magneto-hydrodynamical cosmological zoom-in simulations of galaxy formation and evolution in the $��$CDM framework, finding that bars remain fast down to $z=0$. In Auriga, bars form in galaxies that have higher stellar-to-dark matter ratios and are more baryon-dominated than in previous cosmological simulations; this suggests that in order for bars to remain fast, massive spiral galaxies must lie above the commonly used abundance matching relation. While this reduces the aforementioned tension between the rotation speed of bars and $��$CDM, it accentuates the recently reported discrepancy between the dynamically inferred stellar-to-dark matter ratios of massive spirals and those inferred from abundance matching. Our results highlight the potential of using bar dynamics to constrain models of galaxy formation and evolution., 5 pages + 4 pages of appendices; 7 figures; Accepted for publication in A&A Letters

Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBHs) have to be assembled within the first 800 million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing question whether and how such high black hole accretion rates can be achieved and sustained at the centers of early galaxies. Here we use our new MassiveBlack cosmological hydrodynamic simulation covering a volume (0.75 Gpc)3 appropriate for studying the rare first quasars to show that steady high density cold gas flows responsible for assembling the first galaxies produce the high gas densities that lead to sustained critical accretion rates and hence rapid growth commensurate with the existence of ~109 M ☉ black holes as early as z ~ 7. We find that under these conditions quasar feedback is not effective at stopping the cold gas from penetrating the central regions and hence cannot quench the accretion until the host galaxy reaches . This cold-flow-driven scenario for the formation of quasars implies that they should be ubiquitous in galaxies in the early universe and that major (proto)galaxy mergers are not a requirement for efficient fuel supply and growth, particularly for the earliest SMBHs.

The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyze 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 $02$) evolution fit by $\mathrm{d\;log(Z)}/\mathrm{dz} \approx - 0.064$. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals ($\sim$ 85 per cent at $z=0$) live outside of the ISM, either in an extended gas disk, 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 strength of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the "regulator" style models are best suited for simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR., Comment: 21 pages, 13 figures, submitted to MNRAS; The IllustrisTNG project website can be found at http://www.tng-project.org/

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications., Comment: 21 pages, 7 figures, to appear in Proceedings of the SPPEXA symposium, Lecture Notes in Computational Science and Engineering (LNCSE), Springer

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM 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 to 33 x 10^14 Msol). 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 Mpc^2 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 mu_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., Accepted by MNRAS, conclusions and figures unchanged, minor revision to discussion; 18 pages, 16 figures (appendix 2 pages, 3 figures)

Understanding the formation of the first objects in the universe critically depends on knowing whether the properties of small dark matter structures at high-redshift (z > 15) are different from their more massive lower-redshift counterparts. To clarify this point, we performed a high-resolution N-body simulation of a cosmological volume 1 Mpc/h comoving on a side, reaching the highest mass resolution to date in this regime. We make precision measurements of various physical properties that characterize dark matter haloes (such as the virial ratio, spin parameter, shape, and formation times, etc.) for the high-redshift (z > 15) dark matter mini-haloes we find in our simulation, and compare them to literature results and a moderate-resolution comparison run within a cube of side-length 100 Mpc/h. We find that dark matter haloes at high-redshift have a log-normal distribution of the dimensionless spin parameter centered around �� $\sim$ 0.03, similar to their more massive counterparts. They tend to have a small ratio of the length of the shortest axis to the longest axis (sphericity), and are highly prolate. In fact, haloes of given mass that formed recently are the least spherical, have the highest virial ratios, and have the highest spins. Interestingly, the formation times of our mini-halos depend only very weakly on mass, in contrast to more massive objects. This is expected from the slope of the linear power spectrum of density perturbations at this scale, but despite this difference, dark matter structures at high-redshift share many properties with their much more massive counterparts observed at later times., 17 pages. Accepted for publication in MNRAS

We examine predictions for the quasar luminosity functions (QLF) and quasar clustering at high redshift (z > 4.75) using MassiveBlack, our new hydrodynamic cosmological simulation which includes a self-consistent model for black hole growth and feedback. We show that the model reproduces the Sloan QLF within observational constraints at z >= 5. We find that the high-z QLF is consistent with a redshift-independent occupation distribution of BHs among dark matter halos (which we provide) such that the evolution of the QLF follows that of the halo mass function. The sole exception is the bright-end at z=6 and 7, where BHs in high-mass halos tend to be unusually bright due to extended periods of Eddington growth caused by high density cold flows into the halo center. We further use these luminosity functions to make predictions for the number density of quasars in upcoming surveys, predicting there should be ~119+-28 (~87+-28) quasars detectable in the F125W band of the WIDE (DEEP) fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) from z=5-6, ~19+-7 (~18+-9) from z=6-7, and ~1.7+-1.5 (~1.5+-1.5) from z=7-8. We also investigate quasar clustering, finding that the correlation length is fully consistent with current constraints for Sloan quasars (r_0~17 h^{-1} Mpc at z=4 for quasars above m_i = 20.2), and grows slowly with redshift up to z=6 (r_0~22 h^{-1} Mpc). Finally, we note that the quasar clustering strength depends weakly on luminosity for low L_BH, but gets stronger at higher L_BH as the BHs are found in higher mass halos., 8 pages, 3 figures

The ambient hot intra-halo gas in clusters of galaxies is constantly fed and stirred by in-falling galaxies, a process that can be studied in detail with cosmological hydro- dynamical simulations. However, different numerical methods yield discrepant predic- tions for crucial hydrodynamical processes, leading for example to different entropy profiles in clusters of galaxies. In particular, the widely used Lagrangian smoothed particle hydrodynamics (SPH) scheme is suspected to strongly damp fluid instabili- ties and turbulence, which are both crucial to establish the thermodynamic structure of clusters. In this study, we test to which extent our recently developed Voronoi particle hydrodynamics (VPH) scheme yields different results for the stripping of gas out of infalling galaxies, and for the bulk gas properties of cluster. We consider both the evolution of isolated galaxy models that are exposed to a stream of intra cluster medium or are dropped into cluster models, as well as non-radiative cosmological sim- ulations of cluster formation. We also compare our particle-based method with results obtained with a fundamentally different discretisation approach as implemented in the moving-mesh code AREPO. We find that VPH leads to noticeably faster stripping of gas out of galaxies than SPH, in better agreement with the mesh-code than with SPH. We show that despite the fact that VPH in its present form is not as accurate as the moving mesh code in our investigated cases, its improved accuracy of gradient estimates makes VPH an attractive alternative to SPH., Comment: accepted by MNRAS, 24 pages, 25 figures

Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, \gamma =dln\rho /dlnr ~ r^\alpha $. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=\rho /\sigma ^3 ~ r^-\chi . We show that this is expected from dynamical equilibrium considerations, since Jeans' equations predict that the pseudo-phase-space density profiles of Einasto halos should resemble power laws over a wide range of radii. For the values of \alpha typical of CDM halos, the inner Q profiles of equilibrium halos deviate significantly from a power law only very close to the center, and simulations of extremely high-resolution would be needed to detect such deviations unambiguously. We use an ensemble of halos drawn from the Millennium-II simulation to study which of these two alternatives describe best the mass profile of CDM halos. Our analysis indicates that at the resolution of the best available simulations, both Einasto and power-law PPSD profiles (with adjustable exponents \alpha and \chi, respectively) provide equally acceptable fits to the simulations. A full account of the structure of CDM halos requires understanding how the shape parameters that characterize departures from self-similarity, like \alpha or \chi, are determined by evolutionary history, environment or initial conditions., Comment: 9 pages, 7 figures, submitted to MNRAS

We use a large dark matter simulation of a LambdaCDM model to investigate the clustering and environmental dependence of the number of substructures in a halo. Focusing on redshift z=1, we find that the halo occupation distribution is sensitive at the tens of percent level to the surrounding density and to a lesser extent to asymmetry of the surrounding density distribution. We compute the autocorrelation function of halos as a function of occupation, building on the finding of Wechsler et al. (2006) and Gao and White (2007) that halos (at fixed mass) with more substructure are more clustered. We compute the relative bias as a function of occupation number at fixed mass, finding a strong relationship. At fixed mass, halos in the top 5% of occupation can have an autocorrelation function ~ 1.5-2 times higher than the mean. We also compute the bias as a function of halo mass, for fixed halo occupation. We find that for group and cluster sized halos, when the number of subhalos is held fixed, there is a strong anticorrelation between bias and halo mass. Such a relationship represents an additional challenge to the halo model., Comment: 13 pages, 10 figures, MNRAS submitted

At the heart of any method for computational fluid dynamics lies the question of how the simulated fluid should be discretized. Traditionally, a fixed Eulerian mesh is often employed for this purpose, which in modern schemes may also be adaptively refined during a calculation. Particle-based methods on the other hand discretize the mass instead of the volume, yielding an approximately Lagrangian approach. It is also possible to achieve Lagrangian behavior in mesh-based methods if the mesh is allowed to move with the flow. However, such approaches have often been fraught with substantial problems related to the development of irregularity in the mesh topology. Here we describe a novel scheme that eliminates these weaknesses. It is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points. The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite volume approach, based on a second-order Godunov scheme with an exact Riemann solver. A particularly powerful feature of the approach is that the mesh-generating points can in principle be moved arbitrarily. If they are given the velocity of the local flow, a highly accurate Lagrangian formulation of continuum hydrodynamics is obtained that is free of mesh distortion problems, while it is at the same time fully Galilean-invariant, unlike ordinary Eulerian codes. We describe the formulation and implementation of our new Voronoi-based hydrodynamics, and we discuss a number of illustrative test problems that highlight its performance in practical applications., Comment: 35 pages, 10 figures; invited review for the volume "Tessellations in the Sciences: Virtues, Techniques and Applications of Geometric Tilings", eds. R. van de Weijgaert, G. Vegter, J. Ritzerveld and V. Icke, Springer (accepted)

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations., Comment: The paper contains 23 pages, 13 figures. Accepted to MNRAS (2011, June 28)

We present the first hydrodynamical simulations of structure formation using the new moving mesh code AREPO and compare the results with GADGET simulations based on a traditional smoothed particle hydrodynamics (SPH) technique. The two codes share the same Tree-PM gravity solver and include identical sub-resolution physics, but employ different methods to solve the equations of hydrodynamics. This allows us to assess the impact of hydro-solver uncertainties on the results of cosmological studies of galaxy formation. We focus on predictions for global baryon statistics, such as the cosmic star formation rate density, after we introduce our simulation suite and numerical methods. Properties of individual galaxies and haloes are examined by Keres et al. (2011), while a third paper by Sijacki et al. (2011) uses idealised simulations to analyse the differences between the hydrodynamical schemes. We find that the global baryon statistics differ significantly between the two simulation approaches. AREPO shows higher star formation rates at late times, lower mean temperatures, and different gas mass fractions in characteristic phases of the intergalactic medium, in particular a reduced amount of hot gas. Although both codes use the same implementation of cooling, more gas cools out of haloes in AREPO compared with GADGET towards low redshifts. We show that this is caused by a higher heating rate with SPH in the outer parts of haloes, owing to viscous dissipation of SPH's inherent sonic velocity noise and SPH's efficient damping of subsonic turbulence injected in the halo infall region, and because of a higher efficiency of gas stripping in AREPO. As a result of such differences, AREPO leads also to more disk-like morphologies compared to GADGET. Our results indicate that inaccuracies in hydrodynamic solvers can lead to comparatively large systematic differences., Comment: 35 pages, 28 figures, MNRAS accepted. Added discussion about generic issues with SPH. Volume-rendering movies and high-resolution images can be found at http://www.cfa.harvard.edu/itc/research/movingmeshcosmology

Highly supersonic, compressible turbulence is thought to be of tantamount importance for star formation processes in the interstellar medium. Likewise, cosmic structure formation is expected to give rise to subsonic turbulence in the intergalactic medium, which may substantially modify the thermodynamic structure of gas in virialized dark matter halos and affect small-scale mixing processes in the gas. Numerical simulations have played a key role in characterizing the properties of astrophysical turbulence, but thus far systematic code comparisons have been restricted to the supersonic regime, leaving it unclear whether subsonic turbulence is faithfully represented by the numerical techniques commonly employed in astrophysics. Here we focus on comparing the accuracy of smoothed particle hydrodynamics (SPH) and our new moving-mesh technique AREPO in simulations of driven subsonic turbulence. To make contact with previous results, we also analyze simulations of transsonic and highly supersonic turbulence. We find that the widely employed standard formulation of SPH yields problematic results in the subsonic regime. Instead of building up a Kolmogorov-like turbulent cascade, large-scale eddies are quickly damped close to the driving scale and decay into small-scale velocity noise. Reduced viscosity settings improve the situation, but the shape of the dissipation range differs compared with expectations for a Kolmogorov cascade. In contrast, our moving-mesh technique does yield power-law scaling laws for the power spectra of velocity, vorticity and density, consistent with expectations for fully developed isotropic turbulence. We show that large errors in SPH's gradient estimate and the associated subsonic velocity noise are ultimately responsible for producing inaccurate results in the subsonic regime. In contrast, SPH's performance is much better for supersonic turbulence. [Abridged], Comment: 22 pages, 20 figures, accepted in MNRAS. Includes a rebuttal to arXiv:1111.1255 of D. Price and significant revisions to address referee comments. Conclusions of original submission unchanged

We discuss cosmological hydrodynamic simulations of galaxy formation performed with the new moving-mesh code AREPO, which promises higher accuracy compared with the traditional SPH technique that has been widely employed for this problem. We use an identical set of physics in corresponding simulations carried out with the well-tested SPH code GADGET, adopting also the same high-resolution gravity solver. We are thus able to compare both simulation sets on an object-by-object basis, allowing us to cleanly isolate the impact of different hydrodynamical methods on galaxy and halo properties. In accompanying papers, we focus on an analysis of the global baryonic statistics predicted by the simulation codes, (Vogelsberger et al. 2011) and complementary idealized simulations that highlight the differences between the hydrodynamical schemes (Sijacki et al. 2011). Here we investigate their influence on the baryonic properties of simulated galaxies and their surrounding haloes. We find that AREPO leads to significantly higher star formation rates for galaxies in massive haloes and to more extended gaseous disks in galaxies, which also feature a thinner and smoother morphology than their GADGET counterparts. Consequently, galaxies formed in AREPO have larger sizes and higher specific angular momentum than their SPH correspondents. The more efficient cooling flows in AREPO yield higher densities and lower entropies in halo centers (and the opposite trend in halo outskirts) leading to higher star formation rates of massive galaxies. While both codes agree to acceptable accuracy on a number of baryonic properties of cosmic structures, our results clearly demonstrate that galaxy formation simulations greatly benefit from the use of more accurate hydrodynamical techniques such as AREPO., Comment: 24 pages, 14 figures. MNRAS accepted (in press). Expanded introduction and discussion of galactic angular momentum; added references; one additional figure. Movies and high-resolution images can be found at http://www.cfa.harvard.edu/itc/research/movingmeshcosmology/

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society © 2011 RAS © 2011 The authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends, spherical-overdensity and phase-space-based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allow halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo-finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo, and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Through a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase-space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high-resolution cosmological volume, we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity and peak of the rotation curve). We further suggest to utilize the peak of the rotation curve, vmax, as a proxy for mass, given the arbitrariness in defining a proper halo edge, We are greatly indebted to the ASTROSIM network of the European Science Foundation (Science Meeting 2910) for financially supporting the workshop ‘Haloes going MAD’ held in Miraflores de la Sierra near Madrid in 2010 May where all of this work was initiated.AK is supported by the Spanish Ministerio de Ciencia e Innovación (MICINN) in Spain through the Ramón y Cajal programme as well as the grants AYA 2009-13875-C03-02, AYA2009- 12792-C03-03, CSD2009-00064 and CAM S2009/ESP-1496. SRK acknowledges the support by the MICINN under the Consolider-Ingenio, SyeC project CSD-2007-00050. SP and VQ have also been supported by the MICINN (grants AYA2010-21322-C03-01 and CONSOLIDER2007-00050) and the Generalitat Valenciana (grant PROMETEO-2009-103). MZ is supported by NSF grant AST-0708087. GY acknowledges financial support from MICINN (Spain) under project AYA 2009-13875-C03-02 and the ASTROMADRID project S2009/ESP-1496 financed by Comunidad de Madrid

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse., Comment: 16 pages, 14 figures. Accepted for publication in MNRAS. Minor changes to match published version

We compare the structural properties of galaxies formed in cosmological simulations using the smoothed particle hydrodynamics (SPH) code GADGET with those using the moving-mesh code AREPO. Both codes employ identical gravity solvers and the same sub-resolution physics but use very different methods to track the hydrodynamic evolution of gas. This permits us to isolate the effects of the hydro solver on the formation and evolution of galactic gas disks in GADGET and AREPO haloes with comparable numerical resolution. In a matching sample of GADGET and AREPO haloes we fit simulated gas disks with exponential profiles. We find that the cold gas disks formed using the moving mesh approach have systematically larger disk scale lengths and higher specific angular momenta than their GADGET counterparts across a wide range in halo masses. For low mass galaxies differences between the properties of the simulated galaxy disks are caused by an insufficient number of resolution elements which lead to the artificial angular momentum transfer in our SPH calculation. We however find that galactic disks formed in massive halos, resolved with 10^6 particles/cells, are still systematically smaller in the GADGET run by a factor of ~2. The reasons for this are: 1) The excessive heating of haloes close to the cooling radius due to spurious dissipation of the subsonic turbulence in GADGET; and 2) The efficient delivery of low angular momentum gaseous blobs to the bottom of the potential well. While this large population of gaseous blobs in GADGET originates from the filaments which are pressure confined and fragment due to the SPH surface tension while infalling into hot halo atmospheres, it is essentially absent in the moving mesh calculation, clearly indicating numerical rather than physical origin of the blob material., Comment: 17 pages, 14 figures, MNRAS accepted. Movies and high-resolution images can be found at http://www.cfa.harvard.edu/itc/research/movingmeshcosmology

We present a detailed comparison between the well-known SPH code GADGET and the new moving-mesh code AREPO on a number of hydrodynamical test problems. Through a variety of numerical experiments we establish a clear link between test problems and systematic numerical effects seen in cosmological simulations of galaxy formation. Our tests demonstrate deficiencies of the SPH method in several sectors. These accuracy problems not only manifest themselves in idealized hydrodynamical tests, but also propagate to more realistic simulation setups of galaxy formation, ultimately affecting gas properties in the full cosmological framework, as highlighted in papers by Vogelsberger et al. (2011) and Keres et al. (2011). We find that an inadequate treatment of fluid instabilities in GADGET suppresses entropy generation by mixing, underestimates vorticity generation in curved shocks and prevents efficient gas stripping from infalling substructures. In idealized tests of inside-out disk formation, the convergence rate of gas disk sizes is much slower in GADGET due to spurious angular momentum transport. In simulations where we follow the interaction between a forming central disk and orbiting substructures in a halo, the final disk morphology is strikingly different. In AREPO, gas from infalling substructures is readily depleted and incorporated into the host halo atmosphere, facilitating the formation of an extended central disk. Conversely, gaseous sub-clumps are more coherent in GADGET simulations, morphologically transforming the disk as they impact it. The numerical artefacts of the SPH solver are particularly severe for poorly resolved flows, and thus inevitably affect cosmological simulations due to their hierarchical nature. Our numerical experiments clearly demonstrate that AREPO delivers a physically more reliable solution., Comment: 32 pages, 27 figures, MNRAS accepted. Movies and high-resolution images can be found at http://www.cfa.harvard.edu/itc/research/movingmeshcosmology/

Magnetic fields play an important role in astrophysics on a wide variety of scales, ranging from the Sun and compact objects to galaxies and galaxy clusters. Here we discuss a novel implementation of ideal magnetohydrodynamics (MHD) in the moving mesh code AREPO which combines many of the advantages of Eulerian and Lagrangian methods in a single computational technique. The employed grid is defined as the Voronoi tessellation of a set of mesh-generating points which can move along with the flow, yielding an automatic adaptivity of the mesh and a substantial reduction of advection errors. Our scheme solves the MHD Riemann problem in the rest frame of the Voronoi interfaces using the HLLD Riemann solver. To satisfy the divergence constraint of the magnetic field in multiple dimensions, the Dedner divergence cleaning method is applied. In a set of standard test problems we show that the new code produces accurate results, and that the divergence of the magnetic field is kept sufficiently small to closely preserve the correct physical solution. We also apply the code to two first application problems, namely supersonic MHD turbulence and the spherical collapse of a magnetized cloud. We verify that the code is able to handle both problems well, demonstrating the applicability of this MHD version of AREPO to a wide range of problems in astrophysics., Comment: 11 pages, 9 figures, accepted by MNRAS

We study the properties of simulated high-redshift galaxies using cosmological N-body/gasdynamical runs from the OverWhelmingly Large Simulations (OWLS) project. The runs contrast several feedback implementations of varying effectiveness: from no-feedback, to supernova-driven winds to powerful AGN-driven outflows. These different feedback models result in large variations in the abundance and structural properties of bright galaxies at z=2. We find that feedback affects the baryonic mass of a galaxy much more severely than its spin, which is on average roughly half that of its surrounding dark matter halo in our runs. Feedback induces strong correlations between angular momentum content and galaxy mass that leave their imprint on galaxy scaling relations and morphologies. Encouragingly, we find that galaxy disks are common in moderate-feedback runs, making up typically ~50% of all galaxies at the centers of haloes with virial mass exceeding 1e11 M_sun. The size, stellar masses, and circular speeds of simulated galaxies formed in such runs have properties that straddle those of large star-forming disks and of compact early-type galaxies at z=2. Once the detailed abundance and structural properties of these rare objects are well established it may be possible to use them to gauge the overall efficacy of feedback in the formation of high redshift galaxies., Comment: 16 pages, 12 figures. Accepted for publication in MNRAS. Minor changes to match published version

If dark matter is composed of neutralinos, the gamma-ray radiation produced in their annihilation offers an attractive possibility for dark matter detection. This process may contribute significantly to the extragalactic gamma-ray background (EGB) radiation, which is being measured by the FERMI satellite with unprecedented sensitivity. Using the high-resolution Millennium-II simulation of cosmic structure formation we have produced the first full-sky maps of the expected contribution of dark matter annihilation to the EGB radiation. Our maps include a proper normalization of the signal according to a specific supersymmetric model based on minimal supergravity. The new simulated maps allow a study of the angular power spectrum of the gamma-ray background from dark matter annihilation, which has distinctive features associated with the nature of the annihilation process. Our results are in broad agreement with analytic models for the gamma-ray background, but they also include higher-order correlations not readily accessible in analytic calculations and, in addition, provide detailed spectral information for each pixel. In particular, we find that color maps combining different energies can reveal the cosmic large-scale structure at low and intermediate redshifts., Comment: 7 pages, 5 figures, 2009 Fermi Symposium, eConf Proceedings C091122

Galaxy groups are not scaled down versions of massive galaxy clusters - the hot gas in groups (known as the intragroup medium, IGrM for short) is, on average, less dense than the intracluster medium, implying that one or more non-gravitational processes (e.g., radiative cooling, star formation, and/or feedback) has had a relatively larger effect on groups. In the present study, we compare a number of cosmological hydrodynamic simulations that form part of the OverWhelmingly Large Simulations project to isolate and quantify the effects of cooling and feedback from supernovae (SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing Lagrangian thermal histories of the gas in the different runs, which were all started from identical initial conditions. While radiative cooling, star formation, and SN feedback are all necessary ingredients, only runs that also include AGN feedback are able to successfully reproduce the optical and X-ray properties of groups and low-mass clusters. We isolate how, when, and exactly what gas is heated by AGN. Interestingly, we find that the gas that constitutes the present-day IGrM is that which was not strongly heated by AGN. Instead, the low median density/high median entropy of the gas in present-day groups is achieved by the ejection of lower entropy gas from low-mass progenitor galaxies at high redshift (primarily 2 < z < 4). This corresponds to the epoch when supermassive black holes accreted most of their mass, typically at a rate that is close to the Eddington limit (i.e., when the black holes are in a `quasar mode')., 22 pages, 14 figures. MNRAS, accepted with minor revisions. For a version with high-resolution figures, see http://www.ast.cam.ac.uk/~mccarthy/owls_paper2.pdf

We carry out hydrodynamical simulations of galaxy formation that simultaneously follow radiative transfer of hydrogen-ionizing photons, based on the optically-thin variable Eddington tensor approximation as implemented in the {\small GADGET} code. We consider only star-forming galaxies as sources and examine to what extent they can yield a reasonable reionization history and thermal state of the intergalactic medium at redshifts around $z\sim 3$. This serves as an important benchmark for our self-consistent methodology to simulate galaxy formation and reionization, and for future improvements through accounting of other sources and other wavelength ranges. We find that star formation alone is sufficient for reionizing the Universe by redshift $z\sim6$. For a suitable choice of the escape fraction and the heating efficiency, our models are approximately able to account at the same time for the one-point function and the power spectrum of the Lyman-$\alpha$ forest. The radiation field has an important impact on the star formation rate density in our simulations and significantly lowers the gaseous and stellar fractions in low-mass dark matter halos. Our results thus directly demonstrate the importance of radiative feedback for galaxy formation. The spatial and temporal importance of this effect can be studied accurately with the modelling technique explored here, allowing more faithful simulations of galaxy formation., Comment: 14 pages, 17 figures, published in MNRAS

We analyse the coarse-grained phase-space structure of the six Galaxy-scale dark matter haloes of the Aquarius Project using a state-of-the-art 6D substructure finder. Within r_50, we find that about 35% of the mass is in identifiable substructures, predominantly tidal streams, but including about 14% in self-bound subhaloes. The slope of the differential substructure mass function is close to -2, which should be compared to around -1.9 for the population of self-bound subhaloes. Near r_50 about 60% of the mass is in substructures, with about 30% in self-bound subhaloes. The inner 35 kpc of the highest resolution simulation has only 0.5% of its mass in self-bound subhaloes, but 3.3% in detected substructure, again primarily tidal streams. The densest tidal streams near the solar position have a 3-D mass density about 1% of the local mean, and populate the high velocity tail of the velocity distribution., Comment: Submitted to MNRAS on 12/10/2010, 11 pages, 10 figures

We use a large N-body simulation to examine the detectability of HI in emission at redshift z ~ 1, and the constraints imposed by current observations on the neutral hydrogen mass function of galaxies at this epoch. We consider three different models for populating dark matter halos with HI, designed to encompass uncertainties at this redshift. These models are consistent with recent observations of the detection of HI in emission at z ~ 0.8. Whilst detection of 21 cm emission from individual halos requires extremely long integrations with existing radio interferometers, such as the Giant Meter Radio Telescope (GMRT), we show that the stacked 21 cm signal from a large number of halos can be easily detected. However, the stacking procedure requires accurate redshifts of galaxies. We show that radio observations of the field of the DEEP2 spectroscopic galaxy redshift survey should allow detection of the HI mass function at the 5-12 sigma level in the mass range 10^(11.4) M_sun/h < M_halo < 10^(12.5)M_sun/h, with a moderate amount of observation time. Assuming a larger noise level that corresponds to an upper bound for the expected noise for the GMRT, the detection significance for the HI mass function is still at the 1.7-3 sigma level. We find that optically undetected satellite galaxies enhance the HI emission profile of the parent halo, leading to broader wings as well as a higher peak signal in the stacked profile of a large number of halos. We show that it is in principle possible to discern the contribution of undetected satellites to the total HI signal, even though cosmic variance limitation make this challenging for some of our models., 14 pages, 9 figures, Submitted To MNRAS

We use the OWLS (OverWhelmingly Large Simulations) set of cosmological Nbody/gasdynamical simulations to study the properties of simulated galaxies at z=2. We focus on the effect of supernova feedback from evolving stars on the baryonic mass and angular momentum content of galaxies that assemble at the center of 1e11-3e12 h^{-1}M_\odot halos. Our main finding is that the mass and angular momentum of such galaxies are strongly coupled, in a way that is approximately independent of feedback: varying the feedback implementation leads, in a given halo, to large variations in galaxy mass but leaves the galaxy mass-angular momentum correlation largely unaltered. In particular, the ratio between the angular momentum of a galaxy and that of its surrounding halo (j_d=J_gal/J_vir) correlates closely with the galaxy mass (expressed in units of the virial mass of the halo; m_d=M_gal/M_vir). This correlation differs substantially from the m_d=j_d assumption commonly adopted in semianalytic models of galaxy formation. We use these results to infer the sizes of disk galaxies at z=2 expected in the LCDM scenario and to interpret recent observations of extended disks at z~2 by the SINS collaboration, Comment: Conference Proceeding: "Galaxy Evolution: Emerging Insights and Future Challenges", November 11-14, 2008. The University of Texas at Austin. 8 pages, 7 figures

Deep Chandra exposures of the hot galaxy cluster 1E0657-56 have revealed that the cluster is observed shortly after the first core-passage of a massive infalling subcluster, which is preceded by a prominent bow shock with Mach number M~3. The inferred shock velocity of ~4700 km/s has been commonly interpreted as the velocity of the `bullet' subcluster itself. This velocity is unexpectedly high in the LCDM cosmology, which may require non-trivial modifications in the dark sector to be accommodated if taken at face value. Here we present explicit hydrodynamical toy models of galaxy cluster mergers which very well reproduce the observed dynamical state of 1E0657-56 and the mass models inferred from gravitational lensing observations. However, despite a shock speed of 4500 km/s, the subcluster's mass centroid is moving only with 2600 km/s in the rest frame of the system. The difference arises in part due to a gravitationally induced inflow velocity of the gas ahead of the shock, which amounts to ~1100 km/s for our assumed 10:1 mass ratio of the merger. A second effect is that the shock front moves faster than the subcluster itself. A generic LCDM collision model, where a bullet subcluster with concentration c=7.2 merges with a parent cluster with concentration c=3 on a zero-energy orbit, reproduces all the main observational features seen in 1E0657-56 with good accuracy, suggesting that 1E0657-56 is well in line with expectations from standard cosmological models. In theories with an additional 5th-force in the dark sector, the subcluster can be accelerated beyond the velocity reached in LCDM, and the spatial offset between the X-ray peak and the mass centroid of the subcluster can be significantly enlarged. (abridged), Comment: 15 pages, submitted, 14 figures, animation available at http://www.mpa-garching.mpg.de/~volker/bullet