Source author record

Volker Springel

Volker Springel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

207works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

207 published item(s)

preprint2026arXiv

Rapid emergence of overmassive black holes in the early Universe

The origin of supermassive black holes (SMBHs) remains a long-standing problem in astrophysics. Recent JWST observations reveal an unexpectedly abundant population of overmassive black holes at z>4-6, where the BH masses lie far above local scaling relations and not reproduced by current cosmological models. How such overmassive black holes form and rapidly grow within young galaxies has remained unclear. Here we present fully cosmological radiation-hydrodynamic simulations that, for the first time, self-consistently follow the birth, early growth, and emergent observable signatures of SMBHs in proto-cluster environments. We find that heavy seeds of order $10^6 M_\text{sun}$ naturally form, exceeding typical theoretical expectations by an order of magnitude. These seeds rapidly develop dense, optically thick disks whose strong electron scattering produces broad H$α$ emission comparable to that seen in little red dots (LRDs). Sustained super-Eddington accretion then drives fast growth to $\sim 3 \times 10^7 ~M_\text{sun}$ by $z \sim 8$. These results provide a unified physical scenario in which LRDs correspond to a short-lived, enshrouded phase of heavy-seed formation, naturally evolving into the overmassive quasars detected by JWST and ultimately the progenitors of today's SMBHs.

preprint2025arXiv

NOEMA$^{\rm 3D}$: A first kpc resolution study of a $z\sim1.5$ main sequence barred galaxy channeling gas into a growing bulge

We present a very deep CO(3-2) observation of a massive, gas-rich, main sequence, barred spiral galaxy at $z\approx1.52$. Our data were taken with the IRAM-NOEMA interferometer for a 12-antenna equivalent on-source integration time of $\sim$ 50 hours. We fit the major axis kinematics using forward modelling of a rotating disk, and then subtract the two-dimensional beam convolved best-fit model revealing signatures of planar non-circular motions in the residuals. The inferred in-plane radial velocities are remarkably large, of the order of $\approx60$ km/s. Direct comparisons with a high-resolution, simulated, gas-rich, barred galaxy, obtained with the moving mesh code AREPO and the TNG sub-grid model, show that the observed non-circular gas flows can be explained as radial flows driven by the central bar, with an inferred net inflow rate of the order of the SFR. Given the recent evidence for a higher-than-expected fraction of barred disk galaxies at cosmic noon, our results suggest that rapid gas inflows due to bars could be important evolutionary drivers for the dominant population of star-forming galaxies at the peak epoch of star and galaxy formation.

preprint2022arXiv

Disc instability and bar formation: view from the IllustrisTNG simulations

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

preprint2022arXiv

Early-type galaxy density profiles from IllustrisTNG: III. Effects on outer kinematic structure

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 $γ_{\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 $γ_{\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.

preprint2022arXiv

Formation and fate of low metallicity stars in IllustrisTNG50

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

preprint2022arXiv

Galactic angular momentum in the IllustrisTNG simulation -- I. Connection to morphology, halo spin, and black hole mass

We use the TNG100 simulation of the IllustrisTNG project to investigate the stellar specific angular momenta ($j_{\ast}$) of $\sim$12,000 central galaxies at $z=0$ in a full cosmological context, with stellar masses ($M_{\ast}$) ranging from $10^{9}$ to $10^{12} \, {\rm M}_{\odot}$. We find that the $j_{\ast}$-$M_{\ast}$ relations for early-type and late-type galaxies in IllustrisTNG are in good overall agreement with observations, and that these galaxy types typically `retain' $\sim$10-20 and $\sim$50-60 per cent of their host haloes' specific angular momenta, respectively, with some dependence on the methodology used to measure galaxy morphology. We present results for kinematic as well as visual-like morphological measurements of the simulated galaxies. Next, we explore the scatter in the $j_{\ast}$-$M_{\ast}$ relation with respect to the spin of the dark matter halo and the mass of the supermassive black hole (BH) at the galactic centre. We find that galaxies residing in faster spinning haloes, as well as those hosting less massive BHs, tend to have a higher specific angular momentum. We also find that, at fixed galaxy or halo mass, halo spin and BH mass are anticorrelated with each other, probably as a consequence of more efficient gas flow toward the galactic centre in slowly rotating systems. Finally, we show that halo spin plays an important role in determining galaxy sizes - larger discs form at the centres of faster-rotating haloes - although the trend breaks down for massive galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, roughly the mass scale at which a galaxy's stellar mass becomes dominated by accreted stars.

preprint2022arXiv

LYRA II: Cosmological dwarf galaxy formation with inhomogeneous Population III enrichment

We present the simulation of a $2\times10^9 M_\odot$ halo mass cosmological dwarf galaxy run to $z=0$ at 4 solar mass gas resolution with resolved supernova feedback. We compare three simple subgrid implementations for the inhomogeneous chemical enrichment from Population III stars and compare them to constraints from Local Group dwarf galaxies. The employed model, LYRA, is a novel high resolution galaxy formation model built for the moving mesh code AREPO, which is marked by a resolved multi-phase interstellar medium, single stars and individual supernova events. The resulting reionization relic is characterized by a short ($<1.5$ Gyr) star formation history that is repeatedly brought to a standstill by violent bursts of feedback. Star formation is reignited for a short duration due to a merger at $z\approx4$ and then again at $z\approx0.2-0$ after sustained gas accretion. Our model $z=0$ galaxy matches the stellar mass, size, stellar kinematics and metallicity relations of Local Group dwarf galaxies well. The dark matter profile does not exhibit a core in any version of the model. We show that the host halo masses of Population III stars affect the assembly history of dwarf galaxies. This manifests itself through the initial gaseous collapse in the progenitor halos, affecting the central density of the stellar component and through the accretion of luminous substructure.

preprint2022arXiv

Quenching in Cosmic Sheets: Tracing the Impact of Large Scale Structure Collapse on the Evolution of Dwarf Galaxies

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.

preprint2022arXiv

Stellar migration in the Auriga simulations

We study the presence and importance of stellar migration in the evolution of 17 Milky-Way like disk galaxies with stellar mass $10 < \textrm{log}(M_{*}/{\rm M}_\odot) < 11$ from the Auriga suite of zoom-in cosmological hydrodynamical simulations. We compare the birth radii of the stars to their radii at $z=0$ for each system and present mean values of the strength of stellar migration as a function of radius and stellar age which vary between 1-4 kpc. We also investigate the effect of migration on age and metallicity radial profiles in the disks. We find several cases of age gradient flattening due to migration, but significant changes to metallicity profiles only for older stellar populations and disks that develop a strong bar. Furthermore, we study stellar migration from the perspective of the change of the galactocentric radius ($ΔR$) and orbital guiding centre radius ($ΔR_g$) of stellar particles between given time intervals. We find that stars migrate approximately as a diffusion process only in the outer parts of the disks and for particular galaxies that have a weak bar. Strongly barred galaxies in our sample show larger stellar migration but its timestep evolution is slower-than-diffusion. Finally, we give parametrisations that encapsulate the dependence of the strength of the radial migration as a function of time and radius, for incorporation into (semi-)analytic models of galaxy evolution.

preprint2022arXiv

Sunyaev-Zel'dovich effect and X-ray scaling relations of galaxies, groups and clusters in the IllustrisTNG simulations

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.

preprint2022arXiv

The effects of AGN feedback on the structural and dynamical properties of Milky Way-mass galaxies in cosmological simulations

Feedback from active galactic nuclei (AGN) has become established as a fundamental process in the evolution of the most massive galaxies. Its impact on Milky Way (MW)-mass systems, however, remains comparatively unexplored. In this work, we use the Auriga simulations to probe the impact of AGN feedback on the dynamical and structural properties of galaxies, focussing on the bar, bulge, and disc. We analyse three galaxies -- two strongly and one unbarred/weakly barred -- using three setups: (i) the fiducial Auriga model, which includes both radio and quasar mode feedback, (ii) a setup with no radio mode, and (iii) one with neither the radio nor the quasar mode. When removing the radio mode, gas in the circumgalactic medium cools more efficiently and subsequently settles in an extended disc, with little effect on the inner disc. Contrary to previous studies, we find that although the removal of the quasar mode results in more massive central components, these are in the form of compact discs, rather than spheroidal bulges. Therefore, galaxies without quasar mode feedback are more baryon-dominated and thus prone to forming stronger and shorter bars, which reveals an anti-correlation between the ejective nature of AGN feedback and bar strength. Hence, we report that the effect of AGN feedback (i.e. ejective or preventive) can significantly alter the dynamical properties of MW-like galaxies. Therefore, the observed dynamical and structural properties of MW-mass galaxies can be used as additional constraints for calibrating the efficiency of AGN feedback models.

preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

preprint2022arXiv

The THESAN project: predictions for multi-tracer line intensity mapping in the epoch of reionization

Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

preprint2022arXiv

The THESAN project: properties of the intergalactic medium and its connection to Reionization-era galaxies

The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extra-galactic astronomy. We investigate the IGM properties and their connection to galaxies at $z\geq5.5$ under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel THESAN radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic `late' reionization histories, match available constraints on global IGM properties and reproduce the recently-observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-$α$ transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterising the spectral shape and separation of Lyman-$α$ transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, THESAN reproduces the modulation of the Lyman-$α$ flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation-hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, THESAN produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.

preprint2022arXiv

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

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.

preprint2021arXiv

EAGLE-Auriga: effects of different subgrid models on the baryon cycle around Milky Way-mass galaxies

Modern hydrodynamical simulations reproduce many properties of the real universe. These simulations model various physical processes, but many of these are included using `subgrid models' due to resolution limits. Although different subgrid models have been successful in modelling the effects of supernovae (SNe) feedback on galactic properties, it remains unclear if, and by how much, these differing implementations affect observable halo gas properties. In this work, we use `zoom-in' cosmological initial conditions of two volumes selected to resemble the Local Group (LG) evolved with both the Auriga and EAGLE galaxy formation models. While the subgrid physics models in both simulations reproduce realistic stellar components of $L^\star$ galaxies, they exhibit different gas properties. Namely, Auriga predicts that the Milky Way (MW) is almost baryonically closed, whereas EAGLE suggests that only half of the expected baryons reside within the halo. Furthermore, EAGLE predicts that this baryon deficiency extends to the LG, ($r \leq 1 \mathrm{~Mpc}$). The baryon deficiency in EAGLE is likely due to SNe feedback at high redshift, which generates halo-wide outflows, with high covering fractions and radial velocities, which both eject baryons and significantly impede cosmic gas accretion. Conversely, in Auriga, gas accretion is almost unaffected by feedback. These differences appear to be the result of the different energy injection methods from SNe to gas. Our results suggest that both quasar absorption lines and fast radio burst dispersion measures could constrain these two regimes with future observations.

preprint2021arXiv

Hot and counter-rotating star-forming disk galaxies in IllustrisTNG and their real-world counterparts

A key feature of a large population of low-mass, late-type disk galaxies are star-forming disks with exponential light distributions. They are typically also associated with thin and flat morphologies, blue colors, and dynamically cold stars moving along circular orbits within co-planar thin gas disks. However, the latter features do not necessarily always imply the former, in fact, a variety of different kinematic configurations do exist. In this work, we use the cosmological hydrodynamical IllustrisTNG Simulation to study the nature and origin of dynamically hot, sometimes even counter-rotating, star-forming disk galaxies in the lower stellar mass range (between $5\times 10^9\,\mathrm{M_{\odot}}$ and $2\times 10^{10}\,\mathrm{M_{\odot}}$). We find that being dynamically hot arises in most cases as an induced transient state, for example due to galaxy interactions and merger activities, rather than as an age-dependent evolutionary phase of star-forming disk galaxies. The dynamically hot but still actively star-forming disks show a common feature of hosting kinematically misaligned gas and stellar disks, and centrally concentrated on-going star formation. The former is often accompanied by disturbed gas morphologies, while the latter is reflected in low gas and stellar spins in comparison to their dynamically cold, normal disk counterparts. Interestingly, observed galaxies from MaNGA with kinematic misalignment between gas and stars show remarkably similar general properties as the IllustrisTNG galaxies, and therefore are plausible real-world counterparts. In turn, this allows us to make predictions for the stellar orbits and gas properties of these misaligned galaxies.

preprint2021arXiv

LYRA I: Simulating the multi-phase ISM of a dwarf galaxy with variable energy supernovae from individual stars

We introduce the LYRA project, a new high resolution galaxy formation model built within the framework of the cosmological hydro-dynamical moving mesh code AREPO. The model resolves the multi-phase interstellar medium down to 10 K. It forms individual stars sampled from the initial mass function (IMF), and tracks their lifetimes and death pathways individually. Single supernova (SN) blast waves with variable energy are followed within the hydrodynamic calculation to interact with the surrounding interstellar medium (ISM). In this paper, we present the methods and apply the model to a $10^{10} M_{\odot}$ isolated halo. We demonstrate that the majority of supernovae are Sedov-resolved at our fiducial gas mass resolution of $4 M_{\odot}$. We show that our SN feedback prescription self-consistently produces a hot phase within the ISM that drives significant outflows, reduces the gas density and suppresses star formation. Clustered SN play a major role in enhancing the effectiveness of feedback, because the majority of explosions occur in low density material. Accounting for variable SN energy allows the feedback to respond directly to stellar evolution. We show that the ISM is sensitive to the spatially distributed energy deposition. It strongly affects the outflow behaviour, reducing the mass loading by a factor of 2-3, thus allowing the galaxy to retain a higher fraction of mass and metals. LYRA makes it possible to use a comprehensive multi-physics ISM model directly in cosmological (zoom) simulations of dwarf and higher mass galaxies.

preprint2021arXiv

Magnetogenesis around the first galaxies: the impact of different field seeding processes on galaxy formation

We study the evolution of magnetic fields generated by charge segregation ahead of ionization fronts during the Epoch of Reionization, and their effects on galaxy formation. We compare this magnetic seeding process with the Biermann battery, injection from supernovae, and an imposed seed field at redshift $z\gtrsim127$. Using a suite of self-consistent cosmological and zoom-in simulations based on the Auriga galaxy-formation model, we determine that all mechanisms produce galactic magnetic fields that equally affect galaxy formation, and are nearly indistinguishable at $z\lesssim1.5$. The former is compatible with observed values, while the latter is correlated with the gas metallicity below a seed-dependent redshift. Low-density gas and haloes below a seed-dependent mass threshold retain memory of the initial magnetic field. We produce synthetic Faraday rotation measure maps, showing that they have the potential to constrain the seeding process, although current observations are not yet sensitive enough. Our results imply that the ad-hoc assumption of a primordial seed field - widely used in galaxy formation simulations but of uncertain physical origin - can be replaced by physically-motivated mechanisms for magnetogenesis with negligible impact on galactic properties. Additionally, magnetic fields generated ahead of ionization fronts appear very similar but weaker than those produced by the Biermann battery. Hence, in a realistic scenario where both mechanisms are active, the former will be negligible compared to the latter. Finally, our results highlight that the high-redshift Universe is a fruitful testing ground for our understanding of magnetic fields generation.

preprint2021arXiv

Submillimetre galaxies in cosmological hydrodynamical simulations -- an opportunity for constraining feedback models

Submillimetre galaxies (SMGs) have long posed a challenge for theorists, and self-consistently reproducing the properties of the SMG population in a large-volume cosmological hydrodynamical simulation has not yet been achieved. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects based on the simulated galaxies' star formation rates (SFRs) and dust masses and compare the predicted number counts with observations. We find that the predicted SMG number counts based on IllustrisTNG are significantly less than observed (more than 1 dex at $S_{850} \gtrsim 4$ mJy). The simulation from the original Illustris project yields more SMGs than IllustrisTNG: the predicted counts are consistent with those observed at both $S_{850} \lesssim 5$ mJy and $S_{850} \gtrsim 9$ mJy and only a factor of $\sim 2$ lower than observed at intermediate flux densities. The redshift distribution of SMGs with $S_{850} > 3$ mJy in IllustrisTNG is consistent with the observed distribution, whereas the Illustris redshift distribution peaks at significantly lower redshift (1.5 vs. 2.8). We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, \text{M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and SFRs than in Illustris owing to differences in the sub-grid models for stellar or/and active galactic nucleus (AGN) feedback between the two simulations (we unfortunately cannot isolate the specific cause(s) post hoc). Our results demonstrate that because our method enables predicting SMG number counts in post-processing with a negligible computational expense, SMGs can provide useful constraints for tuning sub-grid models in future large-volume cosmological simulations.

preprint2021arXiv

The IllustrisTNG Simulations: Public Data Release

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.

preprint2020arXiv

Connecting turbulent velocities and magnetic fields in galaxy cluster simulations with active galactic nuclei jets

The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small-scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relationship between large-scale gas kinematics and magnetic fields through non-radiative magnetohydrodynamical simulations of the creation, evolution and disruption of AGN jet-inflated lobes in an isolated Perseus-like galaxy cluster, with and without pre-existing turbulence. In particular, we connect cluster velocity measurements with mock RM maps to highlight their underlying physical connection, which opens up the possibility of comparing turbulence levels in two different observables. For single jet outbursts, we find only a local impact on the velocity field, i.e. the associated increase in velocity dispersion is not volume-filling. Furthermore, in a setup with pre-existing turbulence, this increase in velocity dispersion is largely hidden. We use mock X-ray observations to show that at arcmin resolution, the velocity dispersion is therefore dominated by existing large-scale turbulence and is only minimally altered by the presence of a jet. For the velocity structure of central gas uplifted by buoyantly rising lobes, we find fast, coherent outflows with low velocity dispersion. Our results highlight that projected velocity distributions show complex structures which pose challenges for the interpretation of observations.

preprint2020arXiv

Early-type galaxy density profiles from IllustrisTNG: I. Galaxy correlations and the impact of baryons

We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langleγ^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $σ_{γ^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $γ^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $γ^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $γ^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $γ^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.

preprint2020arXiv

Ejective and preventative: the IllustrisTNG black hole feedback and its effects on the thermodynamics of the gas within and around galaxies

Supermassive black holes (SMBHs) which reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star-formation in massive galaxies. Here we study the $10^{9\textrm{--}12.5}\,\mathrm{M_\odot}$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at \zeq{0}, and show how the three components -- SMBH, galaxy, and circumgalactic medium (CGM) -- are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion BH feedback mode, realised in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star-formation quenching at stellar masses $\gtrsim10^{10.5}\mathrm{M_\odot}$ but also to a change in thermodynamic properties of the (\emph{non}-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\textrm{--}100\,\mathrm{Myr}$ to $1\textrm{--}10\,\mathrm{Gyr}$, effectively ceasing ongoing star-formation and inhibiting radiative cooling and future gas accretion. In practice, the same AGN feedback channel is simultaneously `ejective' and `preventative' and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.

preprint2020arXiv

High redshift JWST predictions from IllustrisTNG: II. Galaxy line and continuum spectral indices and dust attenuation curves

We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${\rm H}_α$ and ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions are consistent with present observations at $z\lesssim 3$ with $\lesssim 0.1\,{\rm dex}$ differences in luminosities. However, the predicted ${\rm H}_α$ luminosity function is $\sim 0.3\,{\rm dex}$ dimmer than the observed one at $z\simeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$Å(D4000) and the UV continuum slope $β$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{\rm UV}$ versus $β$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.

preprint2020arXiv

Magnetising the circumgalactic medium of disk galaxies

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

preprint2020arXiv

Neutron star mergers and rare core-collapse supernovae as sources of r-process enrichment in simulated galaxies

We use cosmological, magnetohydrodynamical simulations of Milky Way-mass galaxies from the Auriga project to study their enrichment with rapid neutron capture (r-process) elements. We implement a variety of enrichment models from both binary neutron star mergers and rare core-collapse supernovae. We focus on the abundances of (extremely) metal-poor stars, most of which were formed during the first ~Gyr of the Universe in external galaxies and later accreted onto the main galaxy. We find that the majority of metal-poor stars are r-process enriched in all our enrichment models. Neutron star merger models result in a median r-process abundance ratio which increases with metallicity, whereas the median trend in rare core-collapse supernova models is approximately flat. The scatter in r-process abundance increases for models with longer delay times or lower rates of r-process producing events. Our results are nearly perfectly converged, in part due to the mixing of gas between mesh cells in the simulations. Additionally, different Milky Way-mass galaxies show only small variation in their respective r-process abundance ratios. Current (sparse and potentially biased) observations of metal-poor stars in the Milky Way seem to prefer rare core-collapse supernovae over neutron star mergers as the dominant source of r-process elements at low metallicity, but we discuss possible caveats to our models. Dwarf galaxies which experience a single r-process event early in their history show highly enhanced r-process abundances at low metallicity, which is seen both in observations and in our simulations. We also find that the elements produced in a single event are mixed with ~10^8 Msun of gas relatively quickly, distributing the r-process elements over a large region.

preprint2020arXiv

Powering galactic super-winds with small-scale AGN winds

We present a new implementation for active galactic nucleus (AGN) feedback through small-scale, ultra-fast winds in the moving-mesh hydrodynamic code AREPO. The wind is injected by prescribing mass, momentum and energy fluxes across a spherical boundary centred on a supermassive black hole according to available constraints for accretion disc winds. After sweeping-up a mass equal to their own, small-scale winds thermalise, powering energy-driven outflows with dynamics, structure and cooling properties in excellent agreement with those of analytic wind solutions. Momentum-driven solutions do not easily occur, because the Compton cooling radius is usually much smaller than the free-expansion radius of the small-scale winds. Through various convergence tests, we demonstrate that our implementation yields wind solutions which are well converged down to the typical resolution achieved in cosmological simulations. We test our model in hydrodynamic simulations of isolated Milky Way - mass galaxies. Above a critical AGN luminosity, initially spherical, small-scale winds power bipolar, energy-driven super-winds that break out of the galactic nucleus, flowing at speeds $> 1000 \rm \, km \, s^{-1}$ out to $\sim 10 \, \rm kpc$. These energy-driven outflows result in moderate, but long-term, reduction in star formation, which becomes more pronounced for higher AGN luminosities and faster small-scale winds. Suppression of star formation proceeds through a rapid mode that involves the removal of the highest-density, nuclear gas and through a slower mode that effectively halts halo gas accretion. Our new implementation makes it possible to model AGN-driven winds in a physically meaningful and validated way in simulations of galaxy evolution, the interstellar medium and black hole accretion flows.

preprint2020arXiv

Redshift Evolution of the Fundamental Plane Relation in the IllustrisTNG Simulation

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 σ^b I^c$ with $R$, $σ$, 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 σ\,+\,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=δ$, where $δ\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σ^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.

preprint2020arXiv

Resolving small-scale cold circumgalactic gas in TNG50

We use the high-resolution TNG50 cosmological magnetohydrodynamical simulation to explore the properties and origin of cold circumgalactic medium (CGM) gas around massive galaxies (M* > 10^11 Msun) at intermediate redshift (z~0.5). We discover a significant abundance of small-scale, cold gas structure in the CGM of 'red and dead' elliptical systems, as traced by neutral HI and MgII. Halos can host tens of thousands of discrete absorbing cloudlets, with sizes of order a kpc or smaller. With a Lagrangian tracer analysis, we show that cold clouds form due to strong drho/rho >> 1 gas density perturbations which stimulate thermal instability. These local overdensities trigger rapid cooling from the hot virialized background medium at ~10^7 K to radiatively inefficient ~10^4 K clouds, which act as cosmologically long-lived, 'stimulated cooling' seeds in a regime where the global halo does not satisfy the classic tcool/tff < 10 criterion. Furthermore, these small clouds are dominated by magnetic rather than thermal pressure, with plasma beta << 1, suggesting that magnetic fields may play an important role. The number and total mass of cold clouds both increase with resolution, and the ~8x10^4 Msun cell mass of TNG50 enables the ~few hundred pc, small-scale CGM structure we observe to form. Finally, we make a preliminary comparison against observations from the COS-LRG, LRG-RDR, COS-Halos, and SDSS LRG surveys. We broadly find that our recent, high-resolution cosmological simulations produce sufficiently high covering fractions of extended, cold gas as observed to surround massive galaxies.

preprint2020arXiv

Separate Universe Simulations with IllustrisTNG: baryonic effects on power spectrum responses and higher-order statistics

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.

preprint2020arXiv

The Arepo public code release

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.

preprint2020arXiv

The effects of cosmic rays on the formation of Milky Way-mass galaxies in a cosmological context

We investigate the impact of cosmic rays (CR) and different modes of CR transport on the properties of Milky Way-mass galaxies in cosmological magneto-hydrodynamical simulations in the context of the AURIGA project. We systematically study how advection, anisotropic diffusion and additional Alfvén-wave cooling affect the galactic disc and the circum-galactic medium (CGM). Global properties such as stellar mass and star formation rate vary little between simulations with and without various CR transport physics, whereas structural properties such as disc sizes, CGM densities or temperatures can be strongly affected. In our simulations, CRs affect the accretion of gas onto galaxies by modifying the CGM flow structure. This alters the angular momentum distribution which manifests itself as a difference in stellar and gaseous disc size. The strength of this effect depends on the CR transport model: CR advection results in the most compact discs while the Alfvén-wave model resembles more the AURIGA model. The advection and diffusion models exhibit large ($r\sim50$ kpc) CR pressure-dominated gas haloes causing a smoother and partly cooler CGM. The additional CR pressure smoothes small-scale density peaks and compensates for the missing thermal pressure support at lower CGM temperatures. In contrast, the Alfvén-wave model is only CR pressure dominated at the disc-halo interface and only in this model the gamma-ray emission from hadronic interactions agrees with observations. In contrast to previous findings, we conclude that details of CR transport are critical for accurately predicting the impact of CR feedback on galaxy formation.

preprint2020arXiv

The fate of disk galaxies in IllustrisTNG clusters

We study the stellar morphological evolution of disc galaxies within clusters in the TNG50 and TNG100 runs from the IllustrisTNG simulation suite. We select satellites of masses $10^{9.7} \leq M_{*,z=0}/\text{M}_{\odot} \leq 10^{11.6}$ residing in clusters of masses $10^{14} \lesssim M_{\text{200c,z=0}}/\text{M}_{\odot} \leq 10^{14.6}$ at $z=0$ and that were discs at accretion according to a kinematic morphology indicator (the circularity fraction). These are traced from the time of accretion to $z=0$ and compared to a control sample of central galaxies mass-matched at accretion. Most cluster discs become non-discy by $z=0$, in stark contrast with the control discs, of which a significant fraction remains discy over the same timescales. Cluster discs become non-discy accompanied by gas removal and star formation quenching, loss of dark matter and little growth or a loss of stellar mass. In contrast, control discs transform while also losing gas mass and quenching, but growing significantly in dark matter and stellar mass. Most cluster satellites change morphologies on similar timescales regardless of stellar mass, in $\sim0.5-4$ Gyr after accretion. Cluster discs that experienced more numerous and closer pericentric passages show the largest change in morphology. Morphological change in all cases requires the presence of a gravitational perturbation to drive stellar orbits to non-discy configurations, along with gas removal/heating to prevent replenishment of the disc through continued star-formation. For cluster discs, the perturbation is impulsive tidal shocking at pericentres and not tidal stripping of outer disc stellar material, whereas for control discs, a combination of mergers and AGN feedback appears to be the key driving force behind morphological transformations.

preprint2020arXiv

The Hestia project: simulations of the Local Group

We present the Hestia simulation suite: High-resolutions Environmental Simulations of The Immediate Area, a set of cosmological simulations of the Local Group. Initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to accurately simulate the local cosmography. Halo pairs that resemble the Local Group are found in low resolutions constrained, dark matter only simulations, and selected for higher resolution magneto hydrodynamic simulation using the Arepo code. Baryonic physics follows the Auriga model of galaxy formation. The simulations contain a high resolution region of 3-5 Mpc in radius from the Local Group midpoint embedded in the correct cosmographic landscape. Within this region a simulated Local Group consisting of a Milky Way and Andromeda like galaxy forms, whose description is in excellent agreement with observations. The simulated Local Group galaxies resemble the Milky Way and Andromeda in terms of their halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution and in some cases, the presence of a Magellanic cloud like object. Because these simulations properly model the Local Group in their cosmographic context, they provide a testing ground for questions where environment is thought to play an important role.

preprint2020arXiv

Universal structure of dark matter haloes over a mass range of 20 orders of magnitude

Cosmological models in which dark matter consists of cold elementary particles predict that the dark halo population should extend to masses many orders of magnitude below those at which galaxies can form. Here we report a cosmological simulation of the formation of present-day haloes over the full range of observed halo masses (20 orders of magnitude) when dark matter is assumed to be in the form of weakly interacting massive particles of mass approximately 100 gigaelectronvolts. The simulation has a full dynamic range of 30 orders of magnitude in mass and resolves the internal structure of hundreds of Earth-mass haloes in as much detail as it does for hundreds of rich galaxy clusters. We find that halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae. Halo mass and concentration are tightly related in a way that depends on cosmology and on the nature of the dark matter. For a fixed mass, the concentration is independent of the local environment for haloes less massive than those of typical galaxies. Haloes over the mass range of 10^3 to 10^11 solar masses contribute about equally (per logarithmic interval) to the luminosity produced by dark matter annihilation, which we find to be smaller than all previous estimates by factors ranging up to one thousand.

preprint2020arXiv

X-ray Signatures of Black Hole Feedback: Hot Galactic Atmospheres in IllustrisTNG and X-ray Observations

Hot gaseous atmospheres that permeate galaxies and extend far beyond their stellar distribution, where they are commonly referred to as the circumgalactic medium (CGM), imprint important information about feedback processes powered by the stellar populations of galaxies and their central supermassive black holes (SMBH). In this work we study the properties of this hot X-ray emitting medium using the IllustrisTNG cosmological simulations. We analyse their mock X-ray spectra, obtained from the diffuse and metal-enriched gas in TNG100 and TNG50, and compare the results with X-ray observations of nearby early-type galaxies. The simulations reproduce the observed X-ray luminosities ($L_{\rm X}$) and temperature ($T_{\rm X})$ at small ($<R_{\rm e}$) and intermediate ($<5R_{\rm e}$) radii reasonably well. We find that the X-ray properties of lower mass galaxies depend on their star formation rates. In particular, in the magnitude range where the star-forming and quenched populations overlap, $M_{\rm K}\sim-24$ $ (M_*\sim10^{10.7}M_\odot)$, we find that the X-ray luminosities of star-forming galaxies are on average about an order of magnitude higher than those of their quenched counterparts. We show that this diversity in $L_{\rm X}$ is a direct manifestation of the quenching mechanism in the simulations, where the galaxies are quenched due to gas expulsion driven by SMBH kinetic feedback. The observed dichotomy in $L_{\rm X}$ is thus an important observable prediction for the SMBH feedback-based quenching mechanisms implemented in state-of-the-art cosmological simulations. While the current X-ray observations of star forming galaxies are broadly consistent with the predictions of the simulations, the observed samples are small and more decisive tests are expected from the sensitive all-sky X-ray survey with eROSITA.

preprint2019arXiv

High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes

We present a new family of very high order accurate direct Arbitrary-Lagrangian-Eulerian (ALE) Finite Volume (FV) and Discontinuous Galerkin (DG) schemes for the solution of nonlinear hyperbolic PDE systems on moving 2D Voronoi meshes that are regenerated at each time step and which explicitly allow topology changes in time. The Voronoi tessellations are obtained from a set of generator points that move with the local fluid velocity. We employ an AREPO-type approach, which rapidly rebuilds a new high quality mesh rearranging the element shapes and neighbors in order to guarantee a robust mesh evolution even for vortex flows and very long simulation times. The old and new Voronoi elements associated to the same generator are connected to construct closed space--time control volumes, whose bottom and top faces may be polygons with a different number of sides. We also incorporate degenerate space--time sliver elements, needed to fill the space--time holes that arise because of topology changes. The final ALE FV-DG scheme is obtained by a redesign of the fully discrete direct ALE schemes of Boscheri and Dumbser, extended here to moving Voronoi meshes and space--time sliver elements. Our new numerical scheme is based on the integration over arbitrary shaped closed space--time control volumes combined with a fully-discrete space--time conservation formulation of the governing PDE system. In this way the discrete solution is conservative and satisfies the GCL by construction. Numerical convergence studies as well as a large set of benchmarks for hydrodynamics and magnetohydrodynamics (MHD) demonstrate the accuracy and robustness of the proposed method. Our numerical results clearly show that the new combination of very high order schemes with regenerated meshes with topology changes lead to substantial improvements compared to direct ALE methods on conforming meshes.

preprint2019arXiv

High redshift JWST predictions from IllustrisTNG: Dust modelling and galaxy luminosity functions

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $Λ$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$α$ versus star formation rate relation, and the H$α$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.

preprint2019arXiv

Satellites of Satellites: The Case for Carina and Fornax

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

preprint2016arXiv

A fully cosmological model of a Monoceros-like ring

We study the vertical structure of a stellar disk obtained from a fully cosmological high-resolution hydrodynamical simulation of the formation of a Milky Way-like galaxy. At the present day, the disk's mean vertical height shows a well-defined and strong pattern, with amplitudes as large as 3 kpc in its outer regions. This pattern is the result of a satellite - host halo - disk interaction and reproduces, qualitatively, many of the observable properties of the Monoceros Ring. In particular we find disk material at the distance of Monoceros ($R \sim$ 12-16 kpc, galactocentric) extending far above the mid plane ( 30$^{\circ}$, $\langle Z \rangle \sim$ 1-2 kpc) in both hemispheres, as well as well-defined arcs of disk material at heliocentric distances $\gtrsim 5$ kpc. The pattern was first excited $\approx 3$ Gyr ago as an $m=1$ mode that later winds up into a leading spiral pattern. Interestingly, the main driver behind this perturbation is a low-mass low-velocity fly-by encounter. The satellite has total mass, pericentre distance and pericentric velocity of $\sim 5\%$ of the host, $\sim 80$ kpc, and 215 km/s, respectively. The satellite is not massive enough to directly perturb the galactic disk but we show that the density field of the host dark matter halo responds to this interaction resulting in a strong amplification of the perturbative effects. This subsequently causes the onset and development of the Monoceros-like feature.

preprint2016arXiv

A moving mesh unstaggered constrained transport scheme for magnetohydrodynamics

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.

preprint2016arXiv

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

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

preprint2016arXiv

Large-scale mass distribution in the Illustris simulation

Observations at low redshifts thus far fail to account for all of the baryons expected in the Universe according to cosmological constraints. A large fraction of the baryons presumably resides in a thin and warm-hot medium between the galaxies, where they are difficult to observe due to their low densities and high temperatures. Cosmological simulations of structure formation can be used to verify this picture and provide quantitative predictions for the distribution of mass in different large-scale structure components. Here we study the distribution of baryons and dark matter at different epochs using data from the Illustris simulation. We identify regions of different dark matter density with the primary constituents of large-scale structure, allowing us to measure mass and volume of haloes, filaments and voids. At redshift zero, we find that 49 % of the dark matter and 23 % of the baryons are within haloes more massive than the resolution limit of $2\times 10^8$ M$_\odot$. The filaments of the cosmic web host a further 45 % of the dark matter and 46 % of the baryons. The remaining 31 % of the baryons reside in voids. The majority of these baryons have been transported there through active galactic nuclei feedback. We note that the feedback model of Illustris is too strong for heavy haloes, therefore it is likely that we are overestimating this amount. Categorizing the baryons according to their density and temperature, we find that 17.8 % of them are in a condensed state, 21.6 % are present as cold, diffuse gas, and 53.9 % are found in the state of a warm-hot intergalactic medium.

preprint2016arXiv

Magnetic Field Amplification During the Common Envelope Phase

During the common envelope (CE) phase, a giant star in a binary system overflows its Roche lobe and unstable mass transfer leads to a spiral-in of the companion, resulting in a close binary system or in a merger of the stellar cores. Dynamo processes during the CE phase have been proposed as a mechanism to generate magnetic fields that are important for forming magnetic white dwarfs (MWDs) and for shaping planetary nebulae. Here, we present the first magnetohydrodynamics simulations of the dynamical spiral-in during a CE phase. We find that magnetic fields are strongly amplified in the accretion stream around the $1M_\odot$ companion as it spirals into the envelope of a $2M_\odot$ RG. This leads to field strengths of 10 to 100 kG throughout the envelope after 120 d. The magnetic field amplification is consistent with being driven by the magnetorotational instability. The field strengths reached in our simulation make the magnetic field interesting for diagnostic purposes, but they are dynamically irrelevant. They are also too small to explain the formation of the highest fields found in MWDs, but may be relevant for luminous red novae, and detecting magnetic fields in these events would support the scenario as proposed here.

preprint2016arXiv

Matter power spectrum and the challenge of percent accuracy

Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day $N$-body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used $N$-body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at $k\leq1$ $h\,\rm Mpc^{-1}$ and to within three percent at $k\leq10$ $h\,\rm Mpc^{-1}$. We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for $k\leq 2$ $h\,\rm Mpc^{-1}$. In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than $L=0.5$ $h^{-1}\rm Gpc$ to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of $M_{\rm p}=10^{9}$ $h^{-1}\rm M_{\odot}$ is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.

preprint2016arXiv

On the stellar halo metallicity profile of Milky Way-like galaxies in the Auriga simulations

A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high resolution cosmological hydrodynamical simulations to study the [Fe/H] distribution of galactic haloes. We find that kinematically selected stellar haloes, including both in-situ and accreted particles, have an oblate [Fe/H] distribution. Spherical [Fe/H] radial profiles show strong negative gradients within 100 kpc, in agreement with previous numerical results. However, the projected median [Fe/H] profiles along the galactic disc minor axis, typically obtained in observations, are significantly flatter. The median [Fe/H] values at a given radius are larger for the spherical profiles than for the minor axis profiles by as much as 0.4 dex within the inner 50 kpc. Similar results are obtained if only the accreted stellar component is considered indicating that the differences between spherical and minor axis profiles are not purely driven by heated disc star particles formed in situ. Our study highlights the importance of performing careful comparisons between models and observations of halo [Fe/H] distributions.

preprint2016arXiv

Rotation curve fitting and its fatal attraction to cores in realistically simulated galaxy observations

We study the role of systematic effects in observational studies of the cusp-core problem under the minimum disc approximation using a suite of high-resolution (25-pc softening length) hydrodynamical simulations of dwarf galaxies. We mimic realistic kinematic observations and fit the mock rotation curves with two analytic models commonly used to differentiate cores from cusps in the dark matter distribution. We find that the cored pseudo-isothermal sphere (ISO) model is strongly favoured by the reduced $χ^2_ν$ of the fits in spite of the fact that our simulations contain cuspy Navarro-Frenk-White profiles (NFW). We show that even idealized measurements of the gas circular motions can lead to the incorrect answer if velocity underestimates induced by pressure support, with a typical size of order $\sim$5 km s$^{-1}$ in the central kiloparsec, are neglected. Increasing the spatial resolution of the mock observations leads to more misleading results because the inner region, where the effect of pressure support is most significant, is better sampled. Fits to observations with a spatial resolution of 100 pc (2 arcsec at 10 Mpc) favour the ISO model in 78-90 per cent of the cases, while at 800-pc resolution, 40-78 per cent of the galaxies indicate the fictitious presence of a dark matter core. The coefficients of our best-fit models agree well with those reported in observational studies; therefore, we conclude that NFW haloes can not be ruled out reliably from this type of analysis.

preprint2016arXiv

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

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

preprint2016arXiv

Shock finding on a moving-mesh: II. Hydrodynamic shocks in the Illustris universe

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.

preprint2016arXiv

Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

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.

preprint2016arXiv

Simulations of ram-pressure stripping in galaxy-cluster interactions

Observationally, the quenching of star-forming galaxies appears to depend both on their mass and environment. The exact cause of the environmental dependence is still poorly understood, yet semi-analytic models (SAMs) of galaxy formation need to parameterise it to reproduce observations of galaxy properties. In this work, we use hydrodynamical simulations to investigate the quenching of disk galaxies through ram-pressure stripping (RPS) as they fall into galaxy clusters with the goal of characterising the importance of this effect for the reddening of disk galaxies. Our set-up employs a live model of a galaxy cluster that interacts with infalling disk galaxies on different orbits. We use the moving-mesh code AREPO, augmented with a special refinement strategy to yield high resolution around the galaxy on its way through the cluster in a computationally efficient way. Our direct simulations differ substantially from stripping models employed in current SAMs, which in most cases overpredict the mass loss from RPS. Furthermore, after pericentre passage, as soon as ram pressure becomes weaker, gas that remains bound to the galaxy is redistributed to the outer parts, an effect that is not captured in simplified treatements of RPS. Star formation in our model galaxies is quenched mainly because the hot gas halo is stripped, depriving the galaxy of its gas supply. The cold gas disk is only stripped completely in extreme cases, leading to full quenching and significant reddening on timescale of ~200 Myr. On the other hand, galaxies experiencing only mild ram pressure actually show an enhanced star formation rate that is consistent with observations and are quenched on timescales > 1 Gyr. Stripped gas in the wake is mixed efficiently with intracluster gas already a few tens of kpc behind the disk, and this gas is free of residual star formation.

preprint2016arXiv

Spiral-induced velocity and metallicity patterns in a cosmological zoom simulation of a Milky Way-sized galaxy

We use a high resolution cosmological zoom simulation of a Milky Way-sized halo to study the observable features in velocity and metallicity space associated with the dynamical influence of spiral arms. For the first time, we demonstrate that spiral arms, that form in a disc in a fully cosmological environment with realistic galaxy formation physics, drive large-scale systematic streaming motions. In particular, on the trailing edge of the spiral arms the peculiar galacto-centric radial and azimuthal velocity field is directed radially outward and azimuthally backward, whereas it is radially inward and azimuthally forward on the leading edge. Owing to the negative radial metallicity gradient, this systematic motion drives, at a given radius, an azimuthal variation in the residual metallicity that is characterised by a metal rich trailing edge and a metal poor leading edge. We show that these signatures are theoretically observable in external galaxies with Integral Field Unit instruments such as VLT/MUSE, and if detected, would provide evidence for large-scale systematic radial migration driven by spiral arms.

preprint2016arXiv

The role of cosmic ray pressure in accelerating galactic outflows

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

preprint2016arXiv

Vertical disc heating in Milky Way-sized galaxies in a cosmological context

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.

preprint2016arXiv

Zoomed cosmological simulations of Milky Way sized halos in f(R)-gravity

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}$.

preprint2016arXiv

Zooming in on major mergers: dense, starbursting gas in cosmological simulations

We introduce the `Illustris zoom simulation project', which allows the study of selected galaxies forming in the $Λ$CDM cosmology with a 40 times better mass resolution than in the parent large-scale hydrodynamical Illustris simulation. We here focus on the starburst properties of the gas in four cosmological simulations of major mergers. The galaxies in our high-resolution zoom runs exhibit a bursty mode of star formation with gas consumption timescales 10 times shorter than for the normal star formation mode. The strong bursts are only present in the simulations with the highest resolution, hinting that a too low resolution is the reason why the original Illustris simulation showed a dearth of starburst galaxies. Very pronounced bursts of star formation occur in two out of four major mergers we study. The high star formation rates, the short gas consumption timescales and the morphology of these systems strongly resemble observed nuclear starbursts. This is the first time that a sample of major mergers is studied through self-consistent cosmological hydrodynamical simulations instead of using isolated galaxy models setup on a collision course. We also study the orbits of the colliding galaxies and find that the starbursting gas preferentially appears in head-on mergers with very high collision velocities. Encounters with large impact parameters do typically not lead to the formation of starbursting gas.

preprint2015arXiv

Accurately simulating anisotropic thermal conduction on a moving mesh

We present a novel implementation of an extremum preserving anisotropic diffusion solver for thermal conduction on the unstructured moving Voronoi mesh of the AREPO code. The method relies on splitting the one-sided facet fluxes into normal and oblique components, with the oblique fluxes being limited such that the total flux is both locally conservative and extremum preserving. The approach makes use of harmonic averaging points and a simple, robust interpolation scheme that works well for strong heterogeneous and anisotropic diffusion problems. Moreover, the required discretisation stencil is small. Efficient fully implicit and semi-implicit time integration schemes are also implemented. We perform several numerical tests that evaluate the stability and accuracy of the scheme, including applications such as point explosions with heat conduction and calculations of convective instabilities in conducting plasmas. The new implementation is suitable for studying important astrophysical phenomena, such as the conductive heat transport in galaxy clusters, the evolution of supernova remnants, or the distribution of heat from blackhole-driven jets into the intracluster medium.

preprint2015arXiv

An analysis of the evolving comoving number density of galaxies in hydrodynamical simulations

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

preprint2015arXiv

Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement

Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second order finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.

preprint2015arXiv

Characterizing the Pressure Smoothing Scale of the Intergalactic Medium

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Ly-alpha forest is degenerate with the 3-D structure of the IGM on ~100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T~$10^4$ K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected Jeans filtering cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Ly-alpha flux, $F_\mathrm{real}$, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The $F_\mathrm{real}$ power spectrum is accurately described by a simple fitting function with cutoff at $λ_F$, allowing us to rigorously quantify the filtering scale for the first time: we find $λ_F$ = 79 kpc (comoving) at z=3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Ly-alpha forest absorption in close quasar pairs, recently proposed as a method to measure the filtering scale. Our results enable one to quantify the filtering scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude $T_0$ and slope $γ$ of the temperature-density relation $T = T_0Δ^{γ-1}$ should be augmented with a third filtering scale parameter $λ_F$.

preprint2015arXiv

Galactic Angular Momentum in the Illustris Simulation: Feedback and the Hubble Sequence

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.

preprint2015arXiv

Galaxy formation in the Planck Cosmology - I. Matching the observed evolution of star formation rates, colours and stellar masses

We have updated the Munich galaxy formation model to the Planck first-year cosmology, while modifying the treatment of baryonic processes to reproduce recent data on the abundance and passive fractions of galaxies from z= 3 down to z=0. Matching these more extensive and more precise observational results requires us to delay the reincorporation of wind ejecta, to lower the surface density threshold for turning cold gas into stars, to eliminate ram-pressure stripping in haloes less massive than ~10^14 Msun, and to modify our model for radio mode feedback. These changes cure the most obvious failings of our previous models, namely the overly early formation of low-mass galaxies and the overly large fraction of them that are passive at late times. The new model is calibrated to reproduce the observed evolution both of the stellar mass function and of the distribution of star formation rate at each stellar mass. Massive galaxies (M>10^11 [Msun]) assemble most of their mass before z=1 and are predominantly old and passive at z=0, while lower mass galaxies assemble later and, for M<10^9.5 (Msun), are still predominantly blue and star forming at z=0. This phenomenological but physically based model allows the observations to be interpreted in terms of the efficiency of the various processes that control the formation and evolution of galaxies as a function of their stellar mass, gas content, environment and time.

preprint2015arXiv

Galaxy Formation with Local Photoionization Feedback -II. Effect of X-Ray Emission from Binaries and Hot Gas

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.

preprint2015arXiv

Galaxy Morphology and Star Formation in the Illustris Simulation at z=0

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.

preprint2015arXiv

How well can cold-dark-matter substructures account for the observed radio flux-ratio anomalies?

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

preprint2015arXiv

Hydrodynamic moving-mesh simulations of the common envelope phase in binary stellar systems

The common envelope (CE) phase is an important stage in binary stellar evolution. It is needed to explain many close binary stellar systems, such as cataclysmic variables, Type Ia supernova progenitors, or X-ray binaries. To form the resulting close binary, the initial orbit has to shrink, thereby transferring energy to the primary giant's envelope that is hence ejected. The details of this interaction, however, are still not understood. Here, we present new hydrodynamic simulations of the dynamical spiral-in forming a CE system. We apply the moving-mesh code AREPO to follow the interaction of a $1M_\odot$ compact star with a $2M_\odot$ red giant possessing a $0.4M_\odot$ core. The nearly Lagrangian scheme combines advantages of smoothed particle hydrodynamics and traditional grid-based hydrodynamic codes and allows us to capture also small flow features at high spatial resolution. Our simulations reproduce the initial transfer of energy and angular momentum from the binary core to the envelope by spiral shocks seen in previous studies, but after about 20 orbits a new phenomenon is observed. Large-scale flow instabilities are triggered by shear flows between adjacent shock layers. These indicate the onset of turbulent convection in the common envelope, thus altering the transport of energy on longer time scales. At the end of our simulation, only 8% of the envelope mass is ejected. The failure to unbind the envelope completely may be caused by processes on thermal time scales or unresolved microphysics.

preprint2015arXiv

Hydrogen Reionization in the Illustris universe

Hydrodynamical simulations of galaxy formation such as the Illustris simulations have progressed to a state where they approximately reproduce the observed stellar mass function from high to low redshift. This in principle allows self-consistent models of reionization that exploit the accurate representation of the diffuse gas distribution together with the realistic growth of galaxies provided by these simulations, within a representative cosmological volume. In this work, we apply and compare two radiative transfer algorithms implemented in a GPU-accelerated code to the $106.5\,{\rm Mpc}$ wide volume of Illustris in postprocessing in order to investigate the reionization transition predicted by this model. We find that the first generation of galaxies formed by Illustris is just about able to reionize the universe by redshift $z\sim 7$, provided quite optimistic assumptions about the escape fraction and the resolution limitations are made. Our most optimistic model finds an optical depth of $τ\simeq 0.065$, which is in very good agreement with recent Planck 2015 determinations. Furthermore, we show that moment-based approaches for radiative transfer with the M1 closure give broadly consistent results with our angular-resolved radiative transfer scheme. In our favoured fiducial model, 20\% of the hydrogen is reionized by redshift $z=9.20$, and this rapidly climbs to 80\% by redshift $z=6.92$. It then takes until $z=6.24$ before 99\% of the hydrogen is ionized. On average, reionization proceeds `inside-out' in our models, with a size distribution of reionized bubbles that progressively features regions of ever larger size while the abundance of small bubbles stays fairly constant.

preprint2015arXiv

Improving the convergence properties of the moving-mesh code AREPO

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

preprint2015arXiv

Modeling Galactic Conformity with the Color-Halo Age Relation in the Illustris Simulation

Comparisons between observational surveys and galaxy formation models find that the mass of dark matter haloes can largely explain galaxies' stellar mass. However, it remains uncertain whether additional environmental variables, generally referred to as assembly bias, are necessary to explain other galaxy properties. We use the Illustris Simulation to investigate the role of assembly bias in producing galactic conformity by considering 18,000 galaxies with $M_{stellar}$ > $2 \times 10^9$ $M_{\odot}$. We find a significant signal of galactic conformity: out to distances of about 10 Mpc, the mean red fraction of galaxies around redder galaxies is higher than around bluer galaxies at fixed stellar mass. Dark matter haloes exhibit an analogous conformity signal, in which the fraction of haloes formed at earlier times (old haloes) is higher around old haloes than around younger ones at fixed halo mass. A plausible interpretation of galactic conformity can be given as a combination of the halo conformity signal with the galaxy color-halo age relation: at fixed stellar mass, particularly toward the low-mass end, Illustris' galaxy colors correlate with halo age, with the reddest galaxies (often satellites) being preferentially found in the oldest haloes. In fact, we can explain the galactic conformity effect with a simple semi-empirical model, by assigning stellar mass based on halo mass (abundance matching) and by assigning galaxy color based on halo age (age matching). We investigate other interpretations for the galactic conformity, particularly its dependence on the isolation criterion and on the central-satellite information. Regarding comparison to observations, we conclude that the adopted selection/isolation criteria, projection effects, and stacking techniques can have a significant impact on the measured amplitude of the conformity signal.

preprint2015arXiv

Planet-disc interaction on a freely moving mesh

General-purpose, moving-mesh schemes for hydrodynamics have opened the possibility of combining the accuracy of grid-based numerical methods with the flexibility and automatic resolution adaptivity of particle-based methods. Due to their supersonic nature, Keplerian accretion discs are in principle a very attractive system for applying such freely moving mesh techniques. However, the high degree of symmetry of simple accretion disc models can be difficult to capture accurately by these methods, due to the generation of geometric grid noise and associated numerical diffusion, which is absent in polar grids. To explore these and other issues, in this work we study the idealized problem of two-dimensional planet-disc interaction with the moving-mesh code AREPO. We explore the hydrodynamic evolution of discs with planets through a series of numerical experiments that vary the planet mass, the disc viscosity and the mesh resolution, and compare the resulting surface density, vortensity field and tidal torque with results from the literature. We find that the performance of the moving-mesh code in this problem is in accordance with published results, showing good consistency with grid codes written in polar coordinates. We also conclude that grid noise and mesh distortions do not introduce excessive numerical diffusion. Finally, we show how the moving-mesh approach can naturally increase resolution in regions of high densityaround planets and planetary wakes, while retaining the background flow at low resolution. This provides an alternative to the difficult task of implementing adaptive mesh refinement in conventional polar-coordinate codes.

preprint2015arXiv

Recoiling black holes: prospects for detection and implications of spin alignment

Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.

preprint2015arXiv

Reducing noise in moving-grid codes with strongly-centroidal Lloyd mesh regularization

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.

preprint2015arXiv

Semi-analytic galaxy formation in coupled dark energy cosmologies

Among the possible alternatives to the standard cosmological model ($Λ$CDM), coupled Dark Energy models postulate that Dark Energy (DE), seen as a dynamical scalar field, may interact with Dark Matter (DM), giving rise to a "fifth-force", felt by DM particles only. In this paper, we study the impact of these cosmologies on the statistical properties of galaxy populations by combining high-resolution numerical simulations with semi-analytic models (SAM) of galaxy formation and evolution. New features have been implemented in the reference SAM in order to have it run self-consistently and calibrated on these cosmological simulations. They include an appropriate modification of the mass temperature relation and of the baryon fraction in DM haloes, due to the different virial scalings and to the gravitational bias, respectively. Our results show that the predictions of our coupled-DE SAM do not differ significantly from theoretical predictions obtained with standard SAMs applied to a reference $Λ$CDM simulation, implying that the statistical properties of galaxies provide only a weak probe for these alternative cosmological models. On the other hand, we show that both galaxy bias and the galaxy pairwise velocity distribution are sensitive to coupled DE models: this implies that these probes might be successfully applied to disentangle among quintessence, $f(R)$-Gravity and coupled DE models.

preprint2015arXiv

Shock finding on a moving-mesh: I. Shock statistics in non-radiative cosmological simulations

Cosmological shock waves play an important role in hierarchical structure formation by dissipating and thermalizing kinetic energy of gas flows, thereby heating the universe. Furthermore, identifying shocks in hydrodynamical simulations and measuring their Mach number accurately is critical for calculating the production of non-thermal particle components through diffusive shock acceleration. However, shocks are often significantly broadened in numerical simulations, making it challenging to implement an accurate shock finder. We here introduce a refined methodology for detecting shocks in the moving-mesh code AREPO, and show that results for shock statistics can be sensitive to implementation details. We put special emphasis on filtering against spurious shock detections due to tangential discontinuities and contacts. Both of them are omnipresent in cosmological simulations, for example in the form of shear-induced Kelvin-Helmholtz instabilities and cold fronts. As an initial application of our new implementation, we analyse shock statistics in non-radiative cosmological simulations of dark matter and baryons. We find that the bulk of energy dissipation at redshift zero occurs in shocks with Mach numbers around ${\cal M}\approx2.7$. Furthermore, almost $40\%$ of the thermalization is contributed by shocks in the warm hot intergalactic medium (WHIM), whereas $\approx60\%$ occurs in clusters, groups and smaller halos. Compared to previous studies, these findings revise the characterization of the most important shocks towards higher Mach numbers and lower density structures. Our results also suggest that regions with densities above and below $δ_b=100$ should be roughly equally important for the energetics of cosmic ray acceleration through large-scale structure shocks.

preprint2015arXiv

Star forming filaments in warm dark matter models

We performed a hydrodynamical cosmological simulation of the formation of a Milky Way-like galaxy in a warm dark matter (WDM) cosmology. Smooth and dense filaments, several co-moving mega parsec long, form generically above z 2 in this model. Atomic line cooling allows gas in the centres of these filaments to cool to the base of the cooling function, resulting in a very striking pattern of extended Lyman-limit systems (LLSs). Observations of the correlation function of LLSs might hence provide useful limits on the nature of the dark matter. We argue that the self-shielding of filaments may lead to a thermal instability resulting in star formation. We implement a sub-grid model for this, and find that filaments rather than haloes dominate star formation until z 6. Reionisation decreases the gas density in filaments, and the more usual star formation in haloes dominates below z 6, although star formation in filaments continues until z=2. Fifteen per cent of the stars of the z=0 galaxy formed in filaments. At higher redshift, these stars give galaxies a stringy appearance, which, if observed, might be a strong indication that the dark matter is warm.

preprint2015arXiv

Stellar orbit evolution in close circumstellar disc encounters

The formation and early evolution of circumstellar discs often occurs within dense, newborn stellar clusters. For the first time, we apply the moving-mesh code AREPO, to circumstellar discs in 3-D, focusing on disc-disc interactions that result from stellar fly-bys. Although a small fraction of stars are expected to undergo close approaches, the outcomes of the most violent encounters might leave an imprint on the discs and host stars that will influence both their orbits and their ability to form planets. We first construct well-behaved 3-D models of self-gravitating discs, and then create a suite of numerical experiments of parabolic encounters, exploring the effects of pericenter separation r_p, disc orientation and disc-star mass ratio (M_d/M_*) on the orbital evolution of the host stars. Close encounters (2r_p<~ disc radius) can truncate discs on very short time scales. If discs are massive, close encounters facilitate enough orbital angular momentum extraction to induce stellar capture. We find that for realistic primordial disc masses M_d<~0.1M_*, non-colliding encounters induce minor orbital changes, which is consistent with analytic calculations of encounters in the linear regime. The same disc masses produce entirely different results for grazing/colliding encounters. In the latter case, rapidly cooling discs lose orbital energy by radiating away the energy excess of the shock-heated gas, thus causing capture of the host stars into a bound orbit. In rare cases, a tight binary with a circumbinary disc forms as a result of this encounter.

preprint2015arXiv

Surface photometry of BCGs and intracluster stars in Lambda-CDM

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.

preprint2015arXiv

The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

preprint2015arXiv

The Illustris Simulation: Public Data Release

We present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)^3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial "full" baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume "subboxes". The total data volume is ~265 TB, including ~800 full volume snapshots and ~30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.

preprint2015arXiv

The Illustris simulation: the evolving population of black holes across cosmic time

We study the properties of black holes and their host galaxies across cosmic time in the Illustris simulation. Illustris is a large scale cosmological hydrodynamical simulation which resolves a (106.5 Mpc)^3 volume with more than 12 billion resolution elements and includes state-of-the-art physical models relevant for galaxy formation. We find that the black hole mass density for redshifts z = 0 - 5 and the black hole mass function at z = 0 predicted by Illustris are in very good agreement with the most recent observational constraints. We show that the bolometric and hard X-ray luminosity functions of AGN at z = 0 and 1 reproduce observational data very well over the full dynamic range probed. Unless the bolometric corrections are largely underestimated, this requires radiative efficiencies to be on average low, epsilon_r <= 0.1, noting however that in our model radiative efficiencies are degenerate with black hole feedback efficiencies. Cosmic downsizing of the AGN population is in broad agreement with the findings from X-ray surveys, but we predict a larger number density of faint AGN at high redshifts than currently inferred. We also study black hole -- host galaxy scaling relations as a function of galaxy morphology, colour and specific star formation rate. We find that black holes and galaxies co-evolve at the massive end, but for low mass, blue and star-forming galaxies there is no tight relation with either their central black hole masses or the nuclear AGN activity.

preprint2015arXiv

The impact of feedback on cosmological gas accretion

We investigate how the way galaxies acquire their gas across cosmic time in cosmological hydrodynamic simulations is modified by a comprehensive physical model for baryonic feedback processes. To do so, we compare two simulations -- with and without feedback -- both evolved with the moving mesh code AREPO. The feedback runs implement the full physics model of the Illustris simulation project, including star formation driven galactic winds and energetic feedback from supermassive blackholes. We explore: (a) the accretion rate of material contributing to the net growth of galaxies and originating directly from the intergalactic medium, finding that feedback strongly suppresses the raw, as well as the net, inflow of this "smooth mode" gas at all redshifts, regardless of the temperature history of newly acquired gas. (b) At the virial radius the temperature and radial flux of inflowing gas is largely unaffected at z=2. However, the spherical covering fraction of inflowing gas at 0.25 rvir decreases substantially, from more than 80% to less than 50%, while the rates of both inflow and outflow increase, indicative of recycling across this boundary. (c) The fractional contribution of smooth accretion to the total accretion rate is lower in the simulation with feedback, by roughly a factor of two across all redshifts. Moreover, the smooth component of gas with a cold temperature history, is entirely suppressed in the feedback run at z<1. (d) The amount of time taken by gas to cross from the virial radius to the galaxy -- the "halo transit time" -- increases in the presence of feedback by a factor of ~2-3, and is notably independent of halo mass. We discuss the possible implications of this invariance for theoretical models of hot halo gas cooling.

preprint2015arXiv

The Impact of Galactic Feedback on the Circumgalactic Medium

Galactic feedback strongly affects the way galactic environments are enriched. We examine this connection by performing a suite of cosmological hydrodynamic simulations, exploring a range of parameters based on the galaxy formation model developed in Vogelsberger et al. 2013 (henceforth V13). We examine the effects of AGN feedback, wind mass loading, wind specific energy, and wind metal-loading on the properties of the circumgalactic medium (CGM) of galaxies with $M_\text{halo} > 10^{11} M_\odot$. Note that while the V13 model was tuned to match observations including the stellar mass function, no explicit tuning was done for the CGM. The wind energy per unit outflow mass has the most significant effect on the CGM enrichment. High energy winds launch metals far beyond the virial radius. AGN feedback also has a significant effect, but only at $z < 3$. We compare to high redshift HI and CIV observations. All our simulations produce the observed number of Damped Lyman-$α$ Absorbers. At lower column density, several of our simulations produce enough Lyman Limit Systems (LLS) $100$ kpc from the galaxy, but in all cases the LLS abundance falls off with distance faster than observations, with too few LLS at $200$ kpc. Further, in all models the CIV abundance drops off too sharply with distance, with too little CIV $100$-$200$ kpc from the galaxy. Higher energy wind models produce more extended CIV but also produce less stars, in tension with star-formation rate density observations. This highlights the fact that circumgalactic observations are a strong constraint on galactic feedback models.

preprint2015arXiv

The merger rate of galaxies in the Illustris Simulation: a comparison with observations and semi-empirical models

We have constructed merger trees for galaxies in the Illustris Simulation by directly tracking the baryonic content of subhalos. These merger trees are used to calculate the galaxy-galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy-galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from `orphaned' galaxies by allowing some objects to `skip' a snapshot when finding a descendant, and by only considering mergers which show a well-defined `infall' moment. Adopting these definitions, we obtain well-converged predictions for the galaxy-galaxy merger rate with the following main features, which are qualitatively similar to the halo-halo merger rate except for the last one: a strong correlation with redshift that evolves as $\sim (1+z)^{2.4-2.8}$, a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than $\sim 2 \times 10^{11} \, {\rm M_{\odot}}$. These trends are consistent with observational constraints for medium-sized galaxies ($M_{\ast} \gtrsim 10^{10} \, {\rm M_{\odot}}$), but in tension with some recent observations of the close pair fraction for massive galaxies ($M_{\ast} \gtrsim 10^{11} \, {\rm M_{\odot}}$), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy-galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.

preprint2015arXiv

The stability of stellar disks in Milky-Way sized dark matter halos

We employ an improved methodology to insert live stellar disks into high-resolution dark matter simulations of Milky Way sized halos, allowing us to investigate the fate of thin stellar disks in the tumultuous environment of cold dark matter structures. We study a set of eight different halos, drawn from the Aquarius simulation project, in which stellar disks are adiabatically grown with a prescribed structure, and then allowed to self-consistently evolve. The initial velocity distribution is set-up in very good equilibrium with the help of the GALIC code. We find that the residual triaxiality of the halos leads to significant disk tumbling, qualitatively confirming earlier work. We show that the disk turning motion is unaffected by structural properties of the galaxies such as the presence or absence of a bulge or bar. In typical Milky Way sized dark matter halos, we expect an average turning of the disks by about 40 degrees between z=1 and z=0, over the coarse of 6 Gyr. We also investigate the impact of the disks on substructures, and conversely, the disk heating rate caused by the dark matter halo substructures. The presence of disks reduces the central subhalo abundance by a about a factor of two, due to an increased evaporation rate by gravitational shocks from disk passages. We find that substructures are important for heating the outer parts of stellar disks but do not appear to significantly affect their inner parts.

preprint2015arXiv

Zooming in on accretion - I. The structure of halo gas

We study the properties of gas in and around 10^12 solar mass halos at z=2 using a suite of high-resolution cosmological hydrodynamic 'zoom' simulations. We quantify the thermal and dynamical structure of these gaseous reservoirs in terms of their mean radial distributions and angular variability along different sightlines. With each halo simulated at three levels of increasing resolution, the highest reaching a baryon mass resolution of ~10,000 solar masses, we study the interaction of filamentary inflow and the quasi-static hot halo atmosphere. We highlight the discrepancy between the spatial resolution available in the halo gas as opposed to within the galaxy itself, and find that stream morphologies become increasingly complex at higher resolution, with large coherent flows revealing density and temperature structure at progressively smaller scales. Moreover, multiple gas components co-exist at the same radius within the halo, making radially averaged analyses misleading. This is particularly true where the hot, quasi-static, high entropy halo atmosphere interacts with cold, rapidly inflowing, low entropy accretion. We investigate the process of gas virialization and identify different regimes for the heating of gas as it accretes from the intergalactic medium. Haloes at this mass have a well-defined virial shock, associated with a sharp jump in temperature and entropy at ~1.25 r_vir. The presence, radius, and radial width of this boundary feature, however, vary not only from halo to halo, but also as a function of angular direction, covering roughly ~85% of the 4pi sphere. Our findings are relevant for the proper interpretation of observations pertaining to the circumgalactic medium, including evidence for large amounts of cold gas surrounding massive haloes at intermediate redshifts.

preprint2014arXiv

An iterative method for the construction of N-body galaxy models in collisionless equilibrium

(Abridged) We describe a new iterative approach for the realization of equilibrium N-body systems for given density distributions. Our method uses elements of Schwarzschild's technique and of the made-to-measure method, but is based on a different principle. Starting with some initial assignment of particle velocities, we calculate their orbits using a flexible tree-based force algorithm in the stationary potential of the target density distribution. The difference of the time-averaged density response produced by these orbits with respect to the initial density configuration is characterized through a merit function, and a stationary solution of the collisionless Boltzmann equation is found by minimizing this merit function directly by iteratively adjusting the initial velocities. Because the distribution function is in general not unique for a given density structure, we augment the merit function with additional constraints that single out a desired target solution. We do this for broad classes of axisymmetric density distributions by numerically solving the Jeans equations to obtain the second velocity moments, which are then imposed as further constraints in the optimization process. The velocity adjustment is carried out with a stochastic process in which new velocities are randomly drawn from an approximate solution of the distribution function, but are kept only when they improve the fit. Our method converges rapidly and is flexible enough to allow the construction of solutions with third integrals of motion, including disk galaxies with realistic velocity distributions. We demonstrate that the new method reproduces analytic distribution functions where they are known (such as the Hernquist sphere), and that it yields very stable N-body realizations for a variety of examples of compound galaxy models, considerably improving on widely used moment-based approaches.

preprint2014arXiv

CO-dark gas and molecular filaments in Milky Way type galaxies

We use the moving mesh code AREPO coupled to a time-dependent chemical network to investigate the formation and destruction of molecular gas in simulated spiral galaxies. This allows us to determine the characteristics of the gas that is not traced by CO emission. Our extremely high resolution AREPO simulations allow us to capture the chemical evolution of the disc, without recourse to a parameterised `clumping factor'. We calculate H2 and CO column densities through our simulated disc galaxies, and estimate the CO emission and CO-H2 conversion factor. We find that in conditions akin to those in the local interstellar medium, around 42% of the total molecular mass should be in CO-dark regions, in reasonable agreement with observational estimates. This fraction is almost insensitive to the CO integrated intensity threshold used to discriminate between CO-bright and CO-dark gas, as long as this threshold is less than 10 K km/s. The CO-dark molecular gas primarily resides in extremely long (>100 pc) filaments that are stretched between spiral arms by galactic shear. Only the centres of these filaments are bright in CO, suggesting that filamentary molecular clouds observed in the Milky Way may only be small parts of much larger structures. The CO-dark molecular gas mainly exists in a partially molecular phase which accounts for a significant fraction of the total disc mass budget. The dark gas fraction is higher in simulations with higher ambient UV fields or lower surface densities, implying that external galaxies with these conditions might have a greater proportion of dark gas.

preprint2014arXiv

Cosmic Degeneracies I: Joint N-body Simulations of Modified Gravity and Massive Neutrinos

We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $Λ$CDM cosmological scenario - namely an $f(R)$ theory of Modified Gravity (MG) and a cosmological background of massive neutrinos - with the aim to investigate their possible observational degeneracies. We focus on three basic statistics of the large-scale matter distribution, more specifically the nonlinear matter power spectrum, the halo mass function, and the halo bias, for which we determine the deviation with respect to the fiducial $Λ$CDM cosmology in the context of both separate and combined simulations of $f(R)$ MG and massive neutrinos scenarios. Our results show that while these two extended models separately determine very prominent and potentially detectable features in all the three statistics, when we allow them to be simultaneously at work these features are strongly suppressed, resulting in much weaker deviations from the standard model's predictions. In particular, when an $f(R)$ gravity model with $f_{R0}=-1\times 10^{-4}$ is combined with a total neutrino mass of $Σ_{i}m_{ν_{i}}=0.4$ eV, the resulting matter power spectrum, halo mass function, and bias at z=0 are found to be consistent with the standard model's predictions at the 10%, 20%, and 5% accuracy levels, respectively. Therefore, our results imply an intrinsic theoretical limit to the effective discriminating power of present and future observational data sets with respect to these widely considered extensions of the standard cosmological scenario in the absence of independent measurements of the neutrino masses from laboratory experiments, even though the high-redshift evolution might still allow to partially break the degeneracy [Abridged].

preprint2014arXiv

Damped Lyman-alpha absorbers as a probe of stellar feedback

We examine the abundance, clustering and metallicity of Damped Lyman-alpha Absorbers (DLAs) in a suite of hydrodynamic cosmological simulations using the moving mesh code AREPO. We incorporate models of supernova and AGN feedback, as well as molecular hydrogen formation. We compare our simulations to the column density distribution function at $z=3$, the total DLA abundance at $z=2-4$, the measured DLA bias at $z=2.3$ and the DLA metallicity distribution at $z=2-4$. Our preferred models produce populations of DLAs in good agreement with most of these observations. The exception is the DLA abundance at $z < 3$, which we show requires stronger feedback in $10^{11-12} \, h^{-1} M_\odot$ mass halos. While the DLA population probes a wide range of halo masses, we find the cross-section is dominated by halos of mass $10^{10} - 10^{11} \, h^{-1} M_\odot$ and virial velocities $50 - 100 \;\mathrm{km/s}$. The simulated DLA population has a linear theory bias of $1.7$, whereas the observations require $2.17 \pm 0.2$. We show that non-linear growth increases the bias in our simulations to $2.3$ at $k=1\; \mathrm{Mpc/}h$, the smallest scale observed. The scale-dependence of the bias is, however, very different in the simulations compared against the observations. We show that, of the observations we consider, the DLA abundance and column density function provide the strongest constraints on the feedback model.

preprint2014arXiv

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

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

preprint2014arXiv

Formation of massive protostars in atomic cooling haloes

We present the highest-resolution three-dimensional simulation to date of the collapse of an atomic cooling halo in the early Universe. We use the moving-mesh code arepo with the primordial chemistry module introduced in Greif (2014), which evolves the chemical and thermal rate equations for over more than 20 orders of magnitude in density. Molecular hydrogen cooling is suppressed by a strong Lyman-Werner background, which facilitates the near-isothermal collapse of the gas at a temperature of about $10^4\,$K. Once the central gas cloud becomes optically thick to continuum emission, it settles into a Keplerian disc around the primary protostar. The initial mass of the protostar is about $0.1\,{\rm M}_\odot$, which is an order of magnitude higher than in minihaloes that cool via molecular hydrogen. The high accretion rate and efficient cooling of the gas catalyse the fragmentation of the disc into a small protostellar system with 5-10 members. After about 12 yr, strong gravitational interactions disrupt the disc and temporarily eject the primary protostar from the centre of the cloud. By the end of the simulation, a secondary clump has collapsed at a distance of $\simeq 150\,$au from the primary clump. If this clump undergoes a similar evolution as the first, the central gas cloud may evolve into a wide binary system. High accretion rates of both the primary and secondary clumps suggest that fragmentation is not a significant barrier for forming at least one massive black hole seed.

preprint2014arXiv

Galaxy mergers on a moving mesh: a comparison with smoothed-particle hydrodynamics

Galaxy mergers have been investigated for decades using smoothed particle hydrodynamics (SPH), but recent work highlighting inaccuracies inherent in the traditional SPH technique calls into question the reliability of previous studies. We explore this issue by comparing a suite of Gadget-3 SPH simulations of idealised (i.e., non-cosmological) isolated discs and galaxy mergers with otherwise identical calculations performed using the moving-mesh code Arepo. When black hole (BH) accretion and active galactic nucleus (AGN) feedback are not included, the star formation histories (SFHs) obtained from the two codes agree well. When BHs are included, the code- and resolution-dependent variations in the SFHs are more significant, but the agreement is still good, and the stellar mass formed over the course of a simulation is robust to variations in the numerical method. During a merger, the gas morphology and phase structure are initially similar prior to the starburst phase. However, once a hot gaseous halo has formed from shock heating and AGN feedback (when included), the agreement is less good. In particular, during the post-starburst phase, the SPH simulations feature more prominent hot gaseous haloes and spurious clumps, whereas with Arepo, gas clumps and filaments are less apparent and the hot halo gas can cool more efficiently. We discuss the origin of these differences and explain why the SPH technique yields trustworthy results for some applications (such as the idealised isolated disc and galaxy merger simulations presented here) but not others (e.g., gas flows onto galaxies in cosmological hydrodynamical simulations).

preprint2014arXiv

Halo Mass and Assembly History Exposed in the Faint Outskirts: the Stellar and Dark Matter Haloes of Illustris Galaxies

We use the Illustris Simulations to gain insight into the build-up of the outer, low-surface brightness regions which surround galaxies. We characterize the stellar haloes by means of the logarithmic slope of the spherically-averaged stellar density profiles, alphaSTARS at z=0, and we relate these slopes to the properties of the underlying Dark-Matter (DM) haloes, their central galaxies, and their assembly histories. We analyze a sample of ~5,000 galaxies resolved with more than 5x10^4 particles each, and spanning a variety of morphologies and halo masses (3x10^11 < Mvir < 10^14 Msun). We find a strong trend between stellar halo slope and total halo mass, where more massive objects have shallower stellar haloes than the less massive ones (-5.5 \pm 0.5 < alphaSTARS <-3.5 \pm 0.2 in the studied mass range). At fixed halo mass, we show that disk-like, blue, young, and more massive galaxies are surrounded by significantly steeper stellar haloes than elliptical, red, older, and less massive galaxies. Overall, the stellar density profiles fall off much more steeply than the underlying DM, and no clear trend holds between stellar slope and DM halo concentration. However, DM haloes which formed more recently, or which accreted larger fractions of stellar mass from infalling satellites, exhibit shallower stellar haloes than their older analogs with similar masses, by up to Delta(alphaSTARS) ~ 0.5-0.7. Our findings, combined with the most recent measurements of the strikingly different stellar power-law indexes for M31 and the Milky Way, appear to favour a massive M31, and a Milky Way characterized by a much quieter accretion history over the past 10 Gyrs than its companion.

preprint2014arXiv

High performance computing and numerical modelling

Numerical methods play an ever more important role in astrophysics. This is especially true in theoretical works, but of course, even in purely observational projects, data analysis without massive use of computational methods has become unthinkable. The key utility of computer simulations comes from their ability to solve complex systems of equations that are either intractable with analytic techniques or only amenable to highly approximative treatments. Simulations are best viewed as a powerful complement to analytic reasoning, and as the method of choice to model systems that feature enormous physical complexity such as star formation in evolving galaxies, the topic of this 43rd Saas Fee Advanced Course. The organizers asked me to lecture about high performance computing and numerical modelling in this winter school, and to specifically cover the basics of numerically treating gravity and hydrodynamics in the context of galaxy evolution. This is still a vast field, and I necessarily had to select a subset of the relevant material. The written notes presented here quite closely follow the lectures as held in Villars-sur-Ollon, which were meant to provide a general overview about some of the most pertinent techniques that may be relevant for students working on numerical models of galaxy evolution and star formation. The discussion is hence often at an introductory level, giving precedence to a presentation of the main numerical concepts rather than to a mathematically detailed exposition of the techniques.

preprint2014arXiv

Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe

We introduce the Illustris Project, a series of large-scale hydrodynamical simulations of galaxy formation. The highest resolution simulation, Illustris-1, covers a volume of $(106.5\,{\rm Mpc})^3$, has a dark mass resolution of ${6.26 \times 10^{6}\,{\rm M}_\odot}$, and an initial baryonic matter mass resolution of ${1.26 \times 10^{6}\,{\rm M}_\odot}$. At $z=0$ gravitational forces are softened on scales of $710\,{\rm pc}$, and the smallest hydrodynamical gas cells have an extent of $48\,{\rm pc}$. We follow the dynamical evolution of $2\times 1820^3$ resolution elements and in addition passively evolve $1820^3$ Monte Carlo tracer particles reaching a total particle count of more than $18$ billion. The galaxy formation model includes: primordial and metal-line cooling with self-shielding corrections, stellar evolution, stellar feedback, gas recycling, chemical enrichment, supermassive black hole growth, and feedback from active galactic nuclei. At $z=0$ our simulation volume contains about $40,000$ well-resolved galaxies covering a diverse range of morphologies and colours including early-type, late-type and irregular galaxies. The simulation reproduces reasonably well the cosmic star formation rate density, the galaxy luminosity function, and baryon conversion efficiency at $z=0$. It also qualitatively captures the impact of galaxy environment on the red fractions of galaxies. The internal velocity structure of selected well-resolved disk galaxies obeys the stellar and baryonic Tully-Fisher relation together with flat circular velocity curves. In the well-resolved regime the simulation reproduces the observed mix of early-type and late-type galaxies. Our model predicts a halo mass dependent impact of baryonic effects on the halo mass function and the masses of haloes caused by feedback from supernova and active galactic nuclei.

preprint2014arXiv

Introducing the Illustris Project: the evolution of galaxy populations across cosmic time

We present an overview of galaxy evolution across cosmic time in the Illustris Simulation. Illustris is an N-body/hydrodynamical simulation that evolves 2*1820^3 resolution elements in a (106.5Mpc)^3 box from cosmological initial conditions down to z=0 using the AREPO moving-mesh code. The simulation uses a state-of-the-art set of physical models for galaxy formation that was tuned to reproduce the z=0 stellar mass function and the history of the cosmic star-formation rate density. We find that Illustris successfully reproduces a plethora of observations of galaxy populations at various redshifts, for which no tuning was performed, and provide predictions for future observations. In particular, we discuss (a) the buildup of galactic mass, showing stellar mass functions and the relations between stellar mass and halo mass from z=7 to z=0, (b) galaxy number density profiles around massive central galaxies out to z=4, (c) the gas and total baryon content of both galaxies and their halos for different redshifts, and as a function of mass and radius, and (d) the evolution of galaxy specific star-formation rates up to z=8. In addition, we (i) present a qualitative analysis of galaxy morphologies from z=5 to z=0, for the stellar as well as the gaseous components, and their appearance in HST mock observations, (ii) follow galaxies selected at z=2 to their z=0 descendants, and quantify their growth and merger histories, and (iii) track massive z=0 galaxies to high redshift and study their joint evolution in star-formation activity and compactness. We conclude with a discussion of several disagreements with observations, and lay out possible directions for future research.

preprint2014arXiv

Properties of galaxies reproduced by a hydrodynamic simulation

Previous simulations of the growth of cosmic structures have broadly reproduced the 'cosmic web' of galaxies that we see in the Universe, but failed to create a mixed population of elliptical and spiral galaxies due to numerical inaccuracies and incomplete physical models. Moreover, because of computational constraints, they were unable to track the small scale evolution of gas and stars to the present epoch within a representative portion of the Universe. Here we report a simulation that starts 12 million years after the Big Bang, and traces 13 billion years of cosmic evolution with 12 billion resolution elements in a volume of $(106.5\,{\rm Mpc})^3$. It yields a reasonable population of ellipticals and spirals, reproduces the distribution of galaxies in clusters and statistics of hydrogen on large scales, and at the same time the metal and hydrogen content of galaxies on small scales.

preprint2014arXiv

Scaling relations and mass bias in hydrodynamical f(R) gravity simulations of galaxy clusters

We investigate the impact of chameleon-type f(R) gravity models on the properties of galaxy clusters and groups. Our f(R) simulations follow for the first time also the hydrodynamics of the intracluster and intragroup medium. This allows us to assess how f(R) gravity alters the X-ray scaling relations of clusters and how hydrostatic and dynamical mass estimates are biased when modifications of gravity are ignored in their determination. We find that velocity dispersions and intracluster medium temperatures are both increased by up to 1/3 in f(R) gravity in low-mass halos, while the difference disappears in massive objects. The mass scale of the transition depends on the background value f_R0 of the scalar degree of freedom. These changes in temperature and velocity dispersion alter the mass-temperature and X-ray luminosity-temperature scaling relations and bias dynamical and hydrostatic mass estimates that do not explicitly account for modified gravity towards higher values. Recently, a relative enhancement of X-ray compared to weak lensing masses was found by the Planck Collaboration (2013). We demonstrate that an explanation for this offset may be provided by modified gravity and the associated bias effects, which interestingly are of the required size. Finally, we find that the abundance of subhalos at fixed cluster mass is only weakly affected by f(R) gravity.

preprint2014arXiv

Statistical properties of dark matter mini-haloes at z >= 15

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.

preprint2014arXiv

Stellar feedback by radiation pressure and photoionization

The relative impact of radiation pressure and photoionization feedback from young stars on surrounding gas is studied with hydrodynamic radiative transfer (RT) simulations. The calculations focus on the single-scattering (direct radiation pressure) and optically thick regime, and adopt a moment-based RT-method implemented in the moving-mesh code AREPO. The source luminosity, gas density profile and initial temperature are varied. At typical temperatures and densities of molecular clouds, radiation pressure drives velocities of order ~20 km/s over 1-5 Myr; enough to unbind the smaller clouds. However, these estimates ignore the effects of photoionization that naturally occur concurrently. When radiation pressure and photoionization act together, the latter is substantially more efficient, inducing velocities comparable to the sound speed of the hot ionized medium (10-15 km/s) on timescales far shorter than required for accumulating similar momentum with radiation pressure. This mismatch allows photoionization to dominate the feedback as the heating and expansion of gas lowers the central densities, further diminishing the impact of radiation pressure. Our results indicate that a proper treatment of the impact of young stars on the interstellar medium needs to primarily account for their ionization power whereas direct radiation pressure appears to be a secondary effect. This conclusion may change if extreme boosts of the radiation pressure by photon trapping are assumed.

preprint2014arXiv

Synthetic Galaxy Images and Spectra from the Illustris Simulation

We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift z = 0 galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.

preprint2014arXiv

The colors of satellite galaxies in the Illustris Simulation

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.

preprint2014arXiv

The Lyman-alpha forest in f(R) modified gravity

In this work, we analyze the Lyman-$α$ forest in cosmological hydrodynamical simulations of chameleon-type f(R) gravity with the goal to assess whether the impact of such models is detectable in absorption line statistics. We carry out a set of hydrodynamical simulations with the cosmological simulation code MG-GADGET, including star formation and cooling effects, and create synthetic Lyman-$α$ absorption spectra from the simulation outputs. We statistically compare simulations with f(R) and ordinary general relativity, focusing on flux probability distribution functions (PDFs) and flux power-spectra, an analysis of the column density and line width distributions, as well as the matter power spectrum. We find that the influence of f(R) gravity on the Lyman-$α$ forest is rather small. Even models with strong modifications of gravity, like $|\bar{f}_{R0}| = 10^{-4}$, do not change the statistical Lyman-$α$ properties by more than 10%. The column density and line width distributions are hardly affected at all. It is therefore not possible to get competitive constraints on the background field $f_R$ using current observational data. An improved understanding of systematics in the observations and a more accurate modeling of the baryonic/radiative physics would be required to allow this in the future. The impact of f(R) on the matter power spectrum in our results is consistent with previous works.

preprint2014arXiv

The star formation main sequence and stellar mass assembly of galaxies in the Illustris simulation

Understanding the physical processes that drive star formation is a key challenge for galaxy formation models. In this article we study the tight correlation between the star formation rate (SFR) and stellar mass of galaxies at a given redshift, how halo growth influences star formation, and star formation histories of individual galaxies. We study these topics using Illustris, a state-of-the-art cosmological hydrodynamical simulation of galaxy formation. Illustris reproduces the observed relation (the star formation main sequence; SFMS) between SFR and stellar mass at redshifts z=0 and z=4, but at intermediate redshifts of z~2, the simulated SFMS has a significantly lower normalisation than reported by observations. The scatter in the relation is consistent with the observed scatter. However, the fraction of outliers above the SFR-stellar mass relation in Illustris is less than that observed. Galaxies with halo masses of ~10^{12} solar masses dominate the SFR density of the Universe, in agreement with the results of abundance matching. Furthermore, more-massive galaxies tend to form the bulk of their stars at high redshift, which indicates that `downsizing' occurs in Illustris. We also studied the star formation histories of individual galaxies, including the use of a principal component analysis decomposition. We find that for fixed stellar mass, galaxies that form earlier have more-massive black holes at z=0, indicating that star formation and black hole growth are tightly linked processes in Illustris. While many of the properties of normal star-forming galaxies are well-reproduced in the Illustris simulation, forming a realistic population of starbursts will likely require higher resolution and probably a more sophisticated treatment of star formation and feedback from stars and black holes.

preprint2013arXiv

A model for cosmological simulations of galaxy formation physics

We present a new comprehensive model of the physics of galaxy formation designed for large-scale hydrodynamical simulations of structure formation using the moving mesh code AREPO. Our model includes primordial and metal line cooling with self-shielding corrections, stellar evolution and feedback processes, gas recycling, chemical enrichment, a novel subgrid model for the metal loading of outflows, black hole (BH) seeding, BH growth and merging procedures, quasar- and radio-mode feedback, and a prescription for radiative electro-magnetic (EM) feedback from active galactic nuclei (AGN). The metal mass loading of outflows can be adjusted independently of the wind mass loading. This is required to simultaneously reproduce the stellar mass content of low mass haloes and their gas oxygen abundances. Radiative EM AGN feedback is implemented assuming an average spectral energy distribution and a luminosity-dependent scaling of obscuration effects. This form of feedback suppresses star formation more efficiently than continuous thermal quasar-mode feedback alone, but is less efficient than mechanical radio-mode feedback in regulating star formation in massive haloes. We contrast simulation predictions for different variants of our galaxy formation model with key observations. Our best match model reproduces, among other things, the cosmic star formation history, the stellar mass function, the stellar mass - halo mass relation, g-, r-, i-, z-band SDSS galaxy luminosity functions, and the Tully-Fisher relation. We can achieve this success only if we invoke very strong forms of stellar and AGN feedback such that star formation is adequately reduced in both low and high mass systems. In particular, the strength of radio-mode feedback needs to be increased significantly compared to previous studies to suppress efficient cooling in massive, metal-enriched haloes.

preprint2013arXiv

A physical model for cosmological simulations of galaxy formation: multi-epoch validation

We present a multi-epoch analysis of the galaxy populations formed within the cosmological hydrodynamical simulations presented in Vogelsberger et al. (2013). These simulations explore the performance of a recently implemented feedback model which includes primordial and metal line radiative cooling with self-shielding corrections; stellar evolution with associated mass loss and chemical enrichment; feedback by stellar winds; black hole seeding, growth and merging; and AGN quasar- and radio-mode heating with a phenomenological prescription for AGN electro-magnetic feedback. We illustrate the impact of the model parameter choices on the resulting simulated galaxy population properties at high and intermediate redshifts. We demonstrate that our scheme is capable of producing galaxy populations that broadly reproduce the observed galaxy stellar mass function extending from redshift z=0 to z=3. We also characterise the evolving galactic B-band luminosity function, stellar mass to halo mass ratio, star formation main sequence, Tully-Fisher relation, and gas-phase mass-metallicity relation and confront them against recent observational estimates. This detailed comparison allows us to validate elements of our feedback model, while also identifying areas of tension that will be addressed in future work.

preprint2013arXiv

Exploring the non-linear density field in the Millennium simulations with tessellations - I. The probability distribution function

We use the Delaunay Tessellation Field Estimator (DTFE) to study the one-point density distribution functions of the Millennium (MS) and Millennium-II (MS-II) simulations. The DTFE technique is based directly on the particle positions, without requiring any type of smoothing or analysis grid, thereby providing high sensitivity to all non-linear structures resolved by the simulations. In order to identify the detailed origin of the shape of the one-point density probability distribution function (PDF), we decompose the simulation particles according to the mass of their host FoF halos, and examine the contributions of different halo mass ranges to the global density PDF. We model the one-point distribution of the FoF halos in each halo mass bin with a set of Monte Carlo realizations of idealized NFW dark matter halos, finding that this reproduces the measurements from the N-body simulations reasonably well, except for a small excess present in simulation results. This excess increases with increasing halo mass. We show that its origin lies in substructure, which becomes progressively more abundant and better resolved in more massive dark matter halos. We demonstrate that the high density tail of the one-point distribution function in less massive halos is severely affected by the gravitational softening length and the mass resolution. In particular, we find these two parameters to be more important for an accurate measurement of the density PDF than the simulated volume. Combining our results from individual halo mass bins we find that the part of the one-point density PDF originating from collapsed halos can nevertheless be quite well described by a simple superposition of a set of NFW halos with the expected cosmological abundance over the resolved mass range. The transition region to the low-density unbound material is however not well captured by such an analytic halo model.

preprint2013arXiv

Following the flow: tracer particles in astrophysical fluid simulations

We present two numerical schemes for passive tracer particles in the hydrodynamical moving-mesh code AREPO, and compare their performance for various problems, from simple setups to cosmological simulations. The purpose of tracer particles is to allow the flow to be followed in a Lagrangian way, tracing the evolution of the fluid with time, and allowing the thermodynamical history of individual fluid parcels to be recorded. We find that the commonly-used `velocity field tracers', which are advected using the fluid velocity field, do not in general follow the mass flow correctly, and explain why this is the case. This method can result in orders-of-magnitude biases in simulations of driven turbulence and in cosmological simulations, rendering the velocity field tracers inappropriate for following these flows. We then discuss a novel implementation of `Monte Carlo tracers', which are moved along with fluid cells, and are exchanged probabilistically between them following the mass flux. This method reproduces the mass distribution of the fluid correctly. The main limitation of this approach is that it is more diffusive than the fluid itself. Nonetheless, we show that this novel approach is more reliable than what has been employed previously and demonstrate that it is appropriate for following hydrodynamical flows in mesh-based codes. The Monte Carlo tracers can also naturally be transferred between fluid cells and other types of particles, such as stellar particles, so that the mass flow in cosmological simulations can be followed in its entirety.

preprint2013arXiv

Galaxy formation on the largest scales: The impact of astrophysics on the BAO peak

We investigate the effects of galaxy formation on the baryonic acoustic oscillations (BAO) peak by applying semi-analytic modelling techniques to the Millennium-XXL, a $3 \times 10^{11}$ particle N-body simulation of similar volume to the future EUCLID survey. Our approach explicitly incorporates the effects of tidal fields and stochasticity on halo formation, as well as the presence of velocity bias, spatially correlated merger histories, and the connection of all these with the observable and physical properties of galaxies. We measure significant deviations in the shape of the BAO peak from the expectations of a linear bias model built on top of the nonlinear dark matter distribution. We find that the galaxy correlation function shows an excess close to the maximum of the BAO peak ($r\sim110 Mpc/h$) and a deficit at $r\sim90 Mpc/h$. Depending on the redshift, selection criteria and number density of the galaxy samples, these bias distortions can be up to 5% in amplitude. They are, however, largely absorbed by marginalization over nuisance parameters in current analytical modelling of the BAO peak in configuration space, in particular into the parameter that controls the broadening due to nonlinear evolution. As a result, the galaxy formation effects detected here are unlikely to bias the high-precision measurements planned by the upcoming generation of wide-field galaxy surveys.

preprint2013arXiv

Modified Gravity-GADGET: A new code for cosmological hydrodynamical simulations of modified gravity models

We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation solver on an adaptive mesh to efficiently solve for perturbations in the scalar degree of freedom of the modified gravity model. As this new algorithm is implemented as a module for the P-Gadget3 code, it can at the same time follow the baryonic physics included in P-Gadget3, such as hydrodynamics, radiative cooling and star formation. We demonstrate that the code works reliably by applying it to simple test problems that can be solved analytically, as well as by comparing cosmological simulations to results from the literature. Using the new code, we perform the first non-radiative and radiative cosmological hydrodynamical simulations of an f(R)-gravity model. We also discuss the impact of AGN feedback on the matter power spectrum, as well as degeneracies between the influence of baryonic processes and modifications of gravity.

preprint2013arXiv

Moving mesh cosmology: tracing cosmological gas accretion

We investigate the nature of gas accretion onto haloes and galaxies at z=2 using cosmological hydrodynamic simulations run with the moving mesh code AREPO. Implementing a Monte Carlo tracer particle scheme to determine the origin and thermodynamic history of accreting gas, we make quantitative comparisons to an otherwise identical simulation run with the smoothed particle hydrodynamics (SPH) code GADGET-3. Contrasting these two numerical approaches, we find significant physical differences in the thermodynamic history of accreted gas in haloes above 10^10.5 solar masses. In agreement with previous work, GADGET simulations show a cold fraction near unity for galaxies forming in massive haloes, implying that only a small percentage of accreted gas heats to an appreciable fraction of the virial temperature during accretion. The same galaxies in AREPO show a much lower cold fraction, <20% in haloes above 10^11 solar masses. This results from a hot gas accretion rate which, at this same halo mass, is an order of magnitude larger than with GADGET, while the cold accretion rate is also lower. These discrepancies increase for more massive systems, and we explain both as due to numerical inaccuracies in the standard formulation of SPH. We also observe that the relatively sharp transition from cold to hot mode dominated accretion, at a halo mass of ~10^11, is a consequence of comparing past gas temperatures to a constant threshold value independent of virial temperature. Examining the spatial distribution of accreting gas, we find that gas filaments in GADGET tend to remain collimated and flow coherently to small radii, or artificially fragment and form a large number of purely numerical "blobs". Similar gas streams in AREPO show increased heating and disruption at 0.25-0.5 virial radii and contribute to the hot gas accretion rate in a manner distinct from classical cooling flows.

preprint2013arXiv

On the operation of the chemothermal instability in primordial star-forming clouds

We investigate the operation of the chemothermal instability in primordial star-forming clouds with a suite of three-dimensional, moving-mesh simulations. In line with previous studies, we find that the gas at the centre of high-redshift minihaloes becomes chemothermally unstable as three-body reactions convert the atomic hydrogen into a fully molecular gas. The competition between the increasing rate at which the gas cools and the increasing optical depth to H2 line emission creates a characteristic dip in the cooling time over the free-fall time on a scale of 100 au. As a result, the free-fall time decreases to below the sound-crossing time, and the cloud may become gravitationally unstable and fragment on a scale of a few tens of au during the initial free-fall phase. In three of the nine haloes investigated, secondary clumps condense out of the parent cloud, which will likely collapse in their own right before they are accreted by the primary clump. In the other haloes, fragmentation at such an early stage is less likely. However, given that previous simulations have shown that the infall velocity decreases substantially once the gas becomes rotationally supported, the amount of time available for perturbations to develop may be much greater than is evident from the limited period of time simulated here.

preprint2013arXiv

Physical properties of simulated galaxy populations at z=2 - I. Effect of metal-line cooling and feedback from star formation and AGN

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on metal-line cooling and feedback from star formation and active galactic nuclei (AGN). We find that if the sub-grid feedback from star formation is implemented kinetically, the feedback is only efficient if the initial wind velocity exceeds a critical value. This critical velocity increases with galaxy mass and also if metal-line cooling is included. This suggests that radiative losses quench the winds if their initial velocity is too low. If the feedback is efficient, then the star formation rate is inversely proportional to the amount of energy injected per unit stellar mass formed (which is proportional to the initial mass loading for a fixed wind velocity). This can be understood if the star formation is self-regulating, i.e. if the star formation rate (and thus the gas fraction) increase until the outflow rate balances the inflow rate. Feedback from AGN is efficient at high masses, while increasing the initial wind velocity with gas pressure or halo mass allows one to generate galaxy-wide outflows at all masses. Matching the observed galaxy mass function requires efficient feedback. In particular, the predicted faint-end slope is too steep unless we resort to highly mass loaded winds for low-mass objects. Such efficient feedback from low-mass galaxies (M_* << 10^10 Msun) also reduces the discrepancy with the observed specific star formation rates, which are higher than predicted unless the feedback transitions from highly efficient to inefficient just below the observed stellar mass range.

preprint2013arXiv

Physical properties of simulated galaxy populations at z=2 - II. Effects of cosmology, reionization and ISM physics

We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on the assumed star formation law, the equation of state imposed on the unresolved interstellar medium, the stellar initial mass function, the reionization history, and the assumed cosmology. This work complements that of Paper I, where we studied the effects of varying models for galactic winds driven by star formation and AGN. The normalisation of the matter power spectrum strongly affects the galaxy mass function, but has a relatively small effect on the physical properties of galaxies residing in haloes of a fixed mass. Reionization suppresses the stellar masses and gas fractions of low-mass galaxies, but by z = 2 the results are insensitive to the timing of reionization. The stellar initial mass function mainly determines the physical properties of galaxies through its effect on the efficiency of the feedback, while changes in the recycled mass and metal fractions play a smaller role. If we use a recipe for star formation that reproduces the observed star formation law independently of the assumed equation of state of the unresolved ISM, then the latter is unimportant. The star formation law, i.e. the gas consumption time scale as a function of surface density, determines the mass of dense, star-forming gas in galaxies, but affects neither the star formation rate nor the stellar mass. This can be understood in terms of self-regulation: the gas fraction adjusts until the outflow rate balances the inflow rate.

preprint2013arXiv

Semi-analytic galaxy formation in f(R)-gravity cosmologies

Modifications of the equations of general relativity at large distances offer one possibility to explain the observed properties of our Universe without invoking a cosmological constant. Numerous proposals for such modified gravity cosmologies exist, but often their consequences for structure formation in the non-linear sector are not yet accurately known. In this work, we employ high-resolution numerical simulations of f(R)-gravity models coupled with a semi-analytic model (SAM) for galaxy formation to obtain detailed predictions for the evolution of galaxy properties. The f(R)-gravity models imply the existence of a `fifth-force', which is however locally suppressed, preserving the successes of general relativity on solar system scales. We show that dark matter haloes in f(R)-gravity models are characterized by a modified virial scaling with respect to the LCDM scenario, reflecting a higher dark matter velocity dispersion at a given mass. This effect is taken into account in the SAM by an appropriate modification of the mass--temperature relation. We find that the statistical properties predicted for galaxies (such as the stellar mass function and the cosmic star formation rate) in f(R)-gravity show generally only very small differences relative to LCDM, smaller than the dispersion between the results of different SAM models, which can be viewed as a measure of their systematic uncertainty. We also demonstrate that galaxy bias is not able to disentangle between f(R)-gravity and the standard cosmological scenario. However, f(R)-gravity imprints modifications in the linear growth rate of cosmic structures at large scale, which can be recovered from the statistical properties of large galaxy samples.

preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

preprint2013arXiv

The formation of disc galaxies in high resolution moving-mesh cosmological simulations

We present cosmological hydrodynamical simulations of eight Milky Way-sized haloes that have been previously studied with dark matter only in the Aquarius project. For the first time, we employ the moving-mesh code AREPO in zoom simulations combined with a comprehensive model for galaxy formation physics designed for large0 cosmological simulations. Our simulations form in most of the eight haloes strongly disc-dominated systems with realistic rotation curves, close to exponential surface density profiles, a stellar-mass to halo-mass ratio that matches expectations from abundance matching techniques, and galaxy sizes and ages consistent with expectations from large galaxy surveys in the local Universe. There is no evidence for any dark matter core formation in our simulations, even so they include repeated baryonic outflows by supernova-driven winds and black hole quasar feedback. For one of our haloes, the object studied in the recent `Aquila' code comparison project, we carried out a resolution study with our techniques, covering a dynamic range of 64 in mass resolution. Without any change in our feedback parameters, the final galaxy properties are reassuringly similar, in contrast to other modelling techniques used in the field that are inherently resolution dependent. This success in producing realistic disc galaxies is reached, in the context of our interstellar medium treatment, without resorting to a high density threshold for star formation, a low star formation efficiency, or early stellar feedback, factors deemed crucial for disc formation by other recent numerical studies.

preprint2013arXiv

The Mass Profile and Accretion History of Cold Dark Matter Halos

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(ρ_{crit}(z)), resembles that of the enclosed density profile, M(<ρ>), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<ρ>) and M(ρ_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.

preprint2013arXiv

The Mass-Concentration-Redshift Relation of Cold Dark Matter Halos

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $ν(M,z)=δ_{\rm crit}(z)/σ(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.

preprint2013arXiv

The satellites of the Milky Way - Insights from semi-analytic modelling in a LambdaCDM cosmology

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

preprint2012arXiv

Formation and evolution of primordial protostellar systems

We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R_sun, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here.

preprint2012arXiv

Gas stripping and mixing in galaxy clusters: A numerical comparison study

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.

preprint2012arXiv

Moving mesh cosmology: characteristics of galaxies and haloes

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.

preprint2012arXiv

Moving mesh cosmology: numerical techniques and global statistics

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.

preprint2012arXiv

Moving Mesh Cosmology: Properties of Gas Disks

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.

preprint2012arXiv

Moving mesh cosmology: properties of neutral hydrogen in absorption

We examine the distribution of neutral hydrogen in cosmological simulations carried out with the new moving-mesh code AREPO and compare it with the corresponding GADGET simulations based on the smoothed particle hydrodynamics (SPH) technique. The two codes use identical gravity solvers and baryonic physics implementations, but very different methods for solving the Euler equations, allowing us to assess how numerical effects associated with the hydro-solver impact the results of simulations. Here we focus on an analysis of the neutral gas, as detected in quasar absorption lines. We find that the high column density regime probed by Damped Lyman-alpha (DLA) and Lyman Limit Systems (LLS) exhibits significant differences between the codes. GADGET produces spurious artefacts in large halos in the form of gaseous clumps, boosting the LLS cross-section. Furthermore, it forms halos with denser central baryonic cores than AREPO, which leads to a substantially greater DLA cross-section from smaller halos. AREPO thus produces a significantly lower cumulative abundance of DLAs, which is intriguingly in much closer agreement with observations. The column density function, however, is not altered enough to significantly reduce the discrepancy with the observed value. For the low column density gas probed by the Lyman-alpha forest, the codes differ only at the level of a few percent, suggesting that this regime is quite well described by both methods, a fact that is reassuring for the many Lyman-alpha studies carried out with SPH thus far. While the residual differences are smaller than the errors on current Lyman-alpha forest data, we note that this will likely change for future precision experiments.

preprint2012arXiv

Moving mesh cosmology: the hydrodynamics of galaxy formation

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.

preprint2012arXiv

Multi-Dimensional, Compressible Viscous Flow on a Moving Voronoi Mesh

Numerous formulations of finite volume schemes for the Euler and Navier-Stokes equations exist, but in the majority of cases they have been developed for structured and stationary meshes. In many applications, more flexible mesh geometries that can dynamically adjust to the problem at hand and move with the flow in a (quasi) Lagrangian fashion would, however, be highly desirable, as this can allow a significant reduction of advection errors and an accurate realization of curved and moving boundary conditions. Here we describe a novel formulation of viscous continuum hydrodynamics that solves the equations of motion on a Voronoi mesh created by a set of mesh-generating points. The points can move in an arbitrary manner, but the most natural motion is that given by the fluid velocity itself, such that the mesh dynamically adjusts to the flow. Owing to the mathematical properties of the Voronoi tessellation, pathological mesh-twisting effects are avoided. Our implementation considers the full Navier-Stokes equations and has been realized in the AREPO code both in 2D and 3D. We propose a new approach to compute accurate viscous fluxes for a dynamic Voronoi mesh, and use this to formulate a finite volume solver of the Navier-Stokes equations. Through a number of test problems, including circular Couette flow and flow past a cylindrical obstacle, we show that our new scheme combines good accuracy with geometric flexibility, and hence promises to be competitive with other highly refined Eulerian methods. This will in particular allow astrophysical applications of the AREPO code where physical viscosity is important, such as in the hot plasma in galaxy clusters, or for viscous accretion disk models.

preprint2012arXiv

Semi-analytic galaxy formation in early dark energy cosmologies

We study the impact of early dark energy (EDE) cosmologies on galaxy properties by coupling high-resolution numerical simulations with semi-analytic modeling (SAM) of galaxy formation and evolution. EDE models are characterized by a non-vanishing high-redshift contribution of dark energy, producing an earlier growth of structures and a modification of large-scale structure evolution. They can be viewed as typical representatives of non-standard dark energy models in which only the expansion history is modified, and hence the impact on galaxy formation is indirect. We show that in EDE cosmologies the predicted space density of galaxies is enhanced at all scales with respect to the standard LCDM scenario, and the corresponding cosmic star formation history and stellar mass density is increased at high-redshift. We compare these results with a set of theoretical predictions obtained with alternative SAMs applied to our reference LCDM simulation, yielding a rough measure of the systematic uncertainty of the models. We find that the modifications in galaxy properties induced by EDE cosmologies are of the same order of magnitude as intra-SAM variations for a standard LCDM realization (unless rather extreme EDE models are considered), suggesting that is difficult to use such predictions alone to disentangle between different cosmological scenarios. However, when independent information on the underlying properties of host dark matter haloes is included, the SAM predictions on galaxy bias may provide important clues on the expansion history and the equation-of-state evolution.

preprint2012arXiv

Shallow Dark Matter Cusps in Galaxy Clusters

We study the evolution of the stellar and dark matter components in a galaxy cluster of $10^{15} \, \rm{M_{\odot}}$ from $z=3$ to the present epoch using the high-resolution collisionless simulations of Ruszkowski & Springel (2009). At $z=3$ the dominant progenitor halos were populated with spherical model galaxies with and without accounting for adiabatic contraction. We apply a weighting scheme which allows us to change the relative amount of dark and stellar material assigned to each simulation particle in order to produce luminous properties which agree better with abundance matching arguments and observed bulge sizes at $z=3$. This permits the study of the effect of initial compactness on the evolution of the mass-size relation. We find that for more compact initial stellar distributions the size of the final Brightest Cluster Galaxy grows with mass according to $r\propto M^{2}$, whereas for more extended initial distributions, $r\propto M$. Our results show that collisionless mergers in a cosmological context can reduce the strength of inner dark matter cusps with changes in logarithmic slope of 0.3 to 0.5 at fixed radius. Shallow cusps such as those found recently in several strong lensing clusters thus do not necessarily conflict with CDM, but may rather reflect on the initial structure of the progenitor galaxies, which was shaped at high redshift by their formation process.

preprint2012arXiv

Shaping the galaxy stellar mass function with supernova- and AGN-driven winds

Cosmological hydrodynamical simulations of galaxy formation in representative regions of the Universe typically need to resort to subresolution models to follow some of the feedback processes crucial for galaxy formation. Here, we show that an energy-driven outflow model in which the wind velocity decreases and the wind mass loading increases in low-mass galaxies, as suggested by observations, can produce a good match to the low-mass end of the observed galaxy stellar mass function. The high-mass end can be recovered simultaneously if feedback from active galactic nuclei (AGN) and a correction for diffuse stellar light plausibly missed in observations are included. At the same time, our model is in good agreement with the stellar mass functions at redshifts z=1 and z=2, and with the observed redshift evolution of the cosmic star formation rate density. In addition, it accurately reproduces the observed gas to stellar mass ratios and specific star formation rates of galaxies as a function of their stellar mass. This agreement with a diverse set of data marks significant progress in hydrodynamically modelling the formation of a representative galaxy population. It also suggests that the mass flux in real galactic winds should strongly increase towards low-mass galaxies. Without this assumption, an overproduction of galaxies at the faint-end of the galaxy luminosity function seems inevitable in our models.

preprint2012arXiv

Subsonic turbulence in smoothed particle hydrodynamics and moving-mesh simulations

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]

preprint2012arXiv

The Lyman-alpha forest in a blazar-heated Universe

It has been realised only recently that TeV emission from blazars can significantly heat the intergalactic medium (IGM) by pair-producing high-energy electrons and positrons, which in turn excite vigorous plasma instabilities, leading to a local dissipation of the pairs' kinetic energy. In this work, we use cosmological hydrodynamical simulations to model the impact of this blazar heating on the Lyman-alpha forest at redshifts z~2-3. We find that blazar heating produces an inverted temperature-density relation in the IGM and naturally resolves many of the problems present in previous simulations of the forest that included photoheating alone. In particular, our simulations with blazar heating simultaneously reproduce the observed effective optical depth and temperature as a function of redshift, the observed probability distribution functions of the transmitted flux, and the observed flux power spectra, over the full redshift range 2<z<3 analysed here. Additionally, by deblending the Lyman-alpha forest into a sum of thermally broadened individual lines, we find superb agreement with the observed lower cutoff of the line-width distribution and abundances of neutral hydrogen column densities. Using the most recent constraints on the cosmic ultraviolet (UV) background, this excellent agreement with observations does not require rescaling the amplitude of the UV background; a procedure that was routinely used in the past to match the observed level of transmitted flux. We also show that our blazar-heated model matches the data better than standard simulations even when such a rescaling is allowed. This concordance between Lyman-alpha data and simulations, which are based on the most recent cosmological parameters, suggests that the inclusion of blazar heating alleviates previous tensions on constraints for sigma_8 derived from Lyman-alpha measurements and other cosmological data. [abridged]

preprint2011arXiv

Cold flows and the first quasars

Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred 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 cosmological hydrodynamic simulation (MassiveBlack) 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 ~10^9 solar mass 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 M_halo > 10^{12} solar masses. 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.

preprint2011arXiv

Dark matter halo occupation: environment and clustering

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.

preprint2011arXiv

Early Black Holes in Cosmological Simulations: Luminosity Functions and Clustering Behaviour

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.

preprint2011arXiv

Formation history, structure and dynamics of discs and spheroids in simulated Milky Way mass galaxies

We study the stellar discs and spheroids in eight simulations of galaxy formation within Milky Way-mass haloes in a Lambda Cold Dark Matter cosmology. A first paper in this series concentrated on disc properties. Here, we extend this analysis to study how the formation history, structure and dynamics of discs and spheroids relate to the assembly history and structure of their haloes. We find that discs are generally young, with stars spanning a wide range in stellar age: the youngest stars define thin discs and have near-circular orbits, while the oldest stars form thicker discs which rotate ~2 times slower than the thin components, and have 2-3 times larger velocity dispersions. Unlike the discs, spheroids form early and on short time-scales, and are dominated by velocity dispersion. We find great variety in their structure. The inner regions are bar- or bulge-like, while the extended outer haloes are rich in complex non-equilibrium structures such as stellar streams, shells and clumps. Our discs have very high in-situ fractions, i.e. most of their stars formed in the disc itself. Nevertheless, there is a non-negligible contribution (~15 percent) from satellites that are accreted on nearly coplanar orbits. The inner regions of spheroids also have relatively high in-situ fractions, but 65-85 percent of their outer stellar population is accreted. We analyse the circular velocities, rotation velocities and velocity dispersions of our discs and spheroids, both for gas and stars, showing that the dynamical structure is complex as a result of the non-trivial interplay between cooling and SN heating.

preprint2011arXiv

From dwarf spheroidals to cDs: Simulating the galaxy population in a LCDM cosmology

We apply updated semi-analytic galaxy formation models simultaneously to the stored halo/subhalo merger trees of the Millennium and Millennium-II simulations. These differ by a factor of 125 in mass resolution, allowing explicit testing of resolution effects on predicted galaxy properties. We have revised the treatments of the transition between the rapid infall and cooling flow regimes of gas accretion, of the sizes of bulges and of gaseous and stellar disks, of supernova feedback, of the transition between central and satellite status as galaxies fall into larger systems, and of gas and star stripping once they become satellites. Plausible values of efficiency and scaling parameters yield an excellent fit not only to the observed abundance of low-redshift galaxies over 5 orders of magnitude in stellar mass and 9 magnitudes in luminosity, but also to the observed abundance of Milky Way satellites. This suggests that reionisation effects may not be needed to solve the "missing satellite" problem except, perhaps, for the faintest objects. The same model matches the observed large-scale clustering of galaxies as a function of stellar mass and colour. The fit remains excellent down to ~30kpc for massive galaxies. For M* < 6 x 10^10Msun, however, the model overpredicts clustering at scales below 1 Mpc, suggesting that the sigma_8 adopted in the simulations (0.9) is too high. Galaxy distributions within rich clusters agree between the simulations and match those observed, but only if galaxies without dark matter subhalos (so-called orphans) are included. Our model predicts a larger passive fraction among low-mass galaxies than is observed, as well as an overabundance of ~10^10Msun galaxies beyond z~0.6, reflecting deficiencies in the way star-formation rates are modelled.

preprint2011arXiv

Gravitational recoils of supermassive black holes in hydrodynamical simulations of gas rich galaxies

We study the evolution of gravitationally recoiled supermassive black holes (BHs) in massive gas-rich galaxies by means of high-resolution hydrodynamical simulations. We find that the presence of a massive gaseous disc allows recoiled BHs to return to the centre on a much shorter timescale than for purely stellar discs. Also, BH accretion and feedback can strongly modify the orbit of recoiled BHs and hence their return timescale, besides affecting the distribution of gas and stars in the galactic centre. However, the dynamical interaction of kicked BHs with the surrounding medium is in general complex and can facilitate both a fast return to the centre as well as a significant delay. The Bondi-Hoyle-Lyttleton accretion rates of the recoiling BHs in our simulated galaxies are favourably high for the detection of off-centred AGN if kicked within gas-rich discs -- up to a few per cent of the Eddington accretion rate -- and are highly variable on timescales of a few 10^7 yrs. In major merger simulations of gas-rich galaxies, we find that gravitational recoils increase the scatter in the BH mass -- host galaxy relationships compared to simulations without kicks, with the BH mass being more sensitive to recoil kicks than the bulge mass. A generic result of our numerical models is that the clumpy massive discs suggested by recent high-redshift observations, as well as the remnants of gas-rich mergers, exhibit a gravitational potential that falls steeply in the central regions, due to the dissipative concentration of baryons. As a result, supermassive BHs should only rarely be able to escape from massive galaxies at high redshifts, which is the epoch where the bulk of BH recoils is expected to occur.[Abridged]

preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] 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 (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows 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. Via 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).

preprint2011arXiv

Hydrodynamic simulations on a moving Voronoi mesh

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.

preprint2011arXiv

Magnetohydrodynamics on an unstructured moving grid

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.

preprint2011arXiv

Properties of fossil groups in cosmological simulations and galaxy formation models

It has been a long-standing question whether fossil groups are just sampling the tail of the distribution of ordinary groups, or whether they are a physically distinct class of objects, characterized by an unusual and special formation history. To study this question, we here investigate fossil groups identified in the hydrodynamical simulations of the GIMIC project, which consists of resimulations of five regions in the Millennium Simulation (MS) that are characterized by different large-scale densities, ranging from a deep void to a proto-cluster region. For comparison, we also consider semi-analytic models built on top of the MS, as well as a conditional luminosity function approach. We identify galaxies in the GIMIC simulations as groups of stars and use a spectral synthesis code to derive their optical properties. The X-ray luminosity of the groups is estimated in terms of the thermal bremsstrahlung emission of the gas in the host halos, neglecting metallicity effects. We focus on comparing the properties of fossil groups in the theoretical models and observational results, highlighting the differences between them, and trying to identify possible dependencies on environment for which our approach is particularly well set-up. We find that the optical fossil fraction in all of our theoretical models declines with increasing halo mass, and there is no clear environmental dependence. Combining the optical and X-ray selection criteria for fossil groups, the halo mass dependence of the fossil groups seen in optical vanishes. Over the GIMIC halo mass range we resolve best, 9.0\times1012 \sim 4.0\times1013 h-1 M, the central galaxies in the fossil groups show similar properties as those in ordinary groups, in terms of age, metallicity, color, concentration, and mass-to-light ratio. [abridged]

preprint2011arXiv

Simulations of galaxy formation with radiative transfer: Hydrogen reionization and radiative feedback

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-$α$ 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.

preprint2011arXiv

Simulations on a Moving Mesh: The Clustered Formation of Population III Protostars

The cosmic dark ages ended a few hundred million years after the Big Bang, when the first stars began to fill the universe with new light. It has generally been argued that these stars formed in isolation and were extremely massive - perhaps 100 times as massive as the Sun. In a recent study, Clark and collaborators showed that this picture requires revision. They demonstrated that the accretion disks that build up around Population III stars are strongly susceptible to fragmentation and that the first stars should therefore form in clusters rather than in isolation. We here use a series of high-resolution hydrodynamical simulations performed with the moving mesh code AREPO to follow up on this proposal and to study the influence of environmental parameters on the level of fragmentation. We model the collapse of five independent minihalos from cosmological initial conditions, through the runaway condensation of their central gas clouds, to the formation of the first protostar, and beyond for a further 1000 years. During this latter accretion phase, we represent the optically thick regions of protostars by sink particles. Gas accumulates rapidly in the circumstellar disk around the first protostar, fragmenting vigorously to produce a small group of protostars. After an initial burst, gravitational instability recurs periodically, forming additional protostars with masses ranging from ~ 0.1 to 10 M_sun. Although the shape, multiplicity, and normalization of the protostellar mass function depend on the details of the sink-particle algorithm, fragmentation into protostars with diverse masses occurs in all cases, confirming earlier reports of Population III stars forming in clusters. Depending on the efficiency of later accretion and merging, Population III stars may enter the main sequence in clusters and with much more diverse masses than are commonly assumed.

preprint2011arXiv

Smoothed Particle Hydrodynamics in Astrophysics

This review discusses Smoothed Particle Hydrodynamics (SPH) in the astrophysical context, with a focus on inviscid gas dynamics. The particle-based SPH technique allows an intuitive and simple formulation of hydrodynamics that has excellent conservation properties and can be coupled to self-gravity easily and highly accurately. The Lagrangian character of SPH allows it to automatically adjust its resolution to the clumping of matter, a property that makes the scheme ideal for many applications in astrophysics, where often a large dynamic range in density is encountered. We discuss the derivation of the basic SPH equations in their modern formulation, and give an overview about extensions of SPH developed to treat physics such as radiative transfer, thermal conduction, relativistic dynamics or magnetic fields. We also briefly describe some of the most important applications areas of SPH in astrophysical research. Finally, we provide a critical discussion of the accuracy of SPH for different hydrodynamical problems, including measurements of its convergence rate for important classes of problems.

preprint2011arXiv

Terapixel imaging of cosmological simulations

The increasing size of cosmological simulations has led to the need for new visualization techniques. We focus on Smoothed Particle Hydrodynamical (SPH) simulations run with the GADGET code and describe methods for visually accessing the entire simulation at full resolution. The simulation snapshots are rastered and processed on supercomputers into images that are ready to be accessed through a web interface (GigaPan). This allows any scientist with a web-browser to interactively explore simulation datasets in both in spatial and temporal dimensions, datasets which in their native format can be hundreds of terabytes in size or more. We present two examples, the first a static terapixel image of the MassiveBlack simulation, a P-GADGET SPH simulation with 65 billion particles, and the second an interactively zoomable animation of a different simulation with more than one thousand frames, each a gigapixel in size. Both are available for public access through the GigaPan web interface. We also make our imaging software publicly available.

preprint2011arXiv

The cosmic X-ray and gamma-ray background from dark matter annihilation

(Abridged) The extragalactic background light (EBL) observed at multiple wavelengths is a promising tool to probe the nature of dark matter since it might contain a significant contribution from gamma-rays produced promptly by dark matter annihilation. Additionally, the electrons and positrons produced in the annihilation give energy to the CMB photons to populate the EBL with X-rays and gamma-rays. We here create full-sky maps of the radiation from both of these contributions using the high-resolution Millennium-II simulation. We use upper limits on the contributions of unknown sources to the EBL to constrain the intrinsic properties of dark matter using a model-independent approach that can be employed as a template to test different particle physics models (including those with a Sommerfeld enhancement). These upper limits are based on observations spanning eight orders of magnitude in energy (from soft X-rays measured by CHANDRA to gamma-rays measured by Fermi), and on expectations for the contributions from blazars and star forming galaxies. To exemplify this approach, we analyze a set of benchmark Sommerfeld-enhanced models that give the correct dark matter abundance, satisfy CMB constraints, and fit the cosmic ray spectra measured by PAMELA and Fermi without any contribution from local subhalos. We find that these models are in conflict with the EBL constraints unless the contribution of unresolved subhalos is small and the annihilation signal dominates the EBL. We conclude that provided the collisionless cold dark matter paradigm is accurate, even for conservative estimates of the contribution from unresolved subhalos and astrophysical backgrounds, the EBL is at least as sensitive a probe of these types of scenarios as the CMB. Our results disfavor an explanation of the positron excess measured by PAMELA based only on dark matter annihilation in the smooth Galactic halo.

preprint2011arXiv

The Delay of Population III Star Formation by Supersonic Streaming Velocities

It has recently been demonstrated that coherent relative streaming velocities of order 30 km / s between dark matter and gas permeated the universe on scales below a few Mpc directly after recombination. We here use a series of high-resolution moving-mesh calculations to show that these supersonic motions significantly influence the virialization of the gas in minihalos, and delay the formation of the first stars. As the gas streams into minihalos with bulk velocities around 1 km / s at z ~ 20, the additional momentum and energy input reduces the gas fractions and central densities of the halos, increasing the typical virial mass required for efficient cooling by a factor of three, and delaying Population III star formation by dz ~ 4. Since the distribution of the magnitude of the streaming velocities is narrowly peaked around a non-negligible value, this effect is important in most regions of the universe. As a consequence, the increased minimum halo mass implies a reduction of the absolute number of minihalos that can be expected to cool and form Population III stars by up to an order of magnitude. We further find that the streaming velocities increase the turbulent velocity dispersion of the minihalo gas, which could affect its ability to fragment and hence alter the mass function of the first stars.

preprint2011arXiv

The Density and Pseudo-Phase-Space Density Profiles of CDM halos

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, γ=dlnρ/dlnr ~ r^α$. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=ρ/σ^3 ~ r^-χ. 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 α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 αand χ, 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 αor χ, are determined by evolutionary history, environment or initial conditions.

preprint2011arXiv

The population of Milky Way satellites in the LambdaCDM cosmology

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.

preprint2011arXiv

The Shape of Dark Matter Haloes in the Aquarius Simulations: Evolution and Memory

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.

preprint2010arXiv

A novel approach for accurate radiative transfer in cosmological hydrodynamic simulations

We present a numerical implementation of radiative transfer based on an explicitly photon-conserving advection scheme, where radiative fluxes over the cell interfaces of a structured or unstructured mesh are calculated with a second-order reconstruction of the intensity field. The approach employs a direct discretisation of the radiative transfer equation in Boltzmann form with adjustable angular resolution that in principle works equally well in the optically thin and optically thick regimes. In our most general formulation of the scheme, the local radiation field is decomposed into a linear sum of directional bins of equal solid-angle, tessellating the unit sphere. Each of these "cone-fields" is transported independently, with constant intensity as a function of direction within the cone. Photons propagate at the speed of light (or optionally using a reduced speed of light approximation to allow larger timesteps), yielding a fully time-dependent solution of the radiative transfer equation that can naturally cope with an arbitrary number of sources, as well as with scattering. The method casts sharp shadows, subject to the limitations induced by the adopted angular resolution. If the number of point sources is small and scattering is unimportant, our implementation can alternatively treat each source exactly in angular space, producing shadows whose sharpness is only limited by the grid resolution. A third hybrid alternative is to treat only a small number of the locally most luminous point sources explicitly, with the rest of the radiation intensity followed in a radiative diffusion approximation. We have implemented the method in the moving-mesh code {\small AREPO}, where it is coupled to the hydrodynamics in an operator splitting approach that subcycles the radiative transfer alternatingly with the hydrodynamical evolution steps.

preprint2010arXiv

A Scaling Relation of the Evolving Tidal Fields in a LCDM Cosmology

We report the finding of a scaling relation among the cosmic-web anisotropy parameter $A$, the linear density rms fluctuation sigma(r) and the linear growth factor D(z). Using the tidal field derived from the Millennium Simulation on 512^{3} grids at z=0, 2, 5 and 127, we calculate the largest eigenvalues $λ$ of the local tidal tensor at each grid resolution and measure its distance-averaged two-point correlation function, xi_{lambda}, as a function of the cosines of polar angles cos(theta) in the local principal axis frame. We show that xi_{lambda} is quite anisotropic, increasing toward the directions of minimal matter compression, and that the anisotropy of xi_{lambda} increases as the redshift, z, decreases and as the upper distance cutoff r_{c} decreases. Fitting the numerical results to an analytic fitting model ξ_λ(\cosθ)\propto (1+A\cos^{n}θ)^{-1}, it is found that the best fit value of A, dubbed the cosmic-web anisotropy parameter, varies systematically with sigma(r_{c}) and D(z), allowing us to determine the simple empiral scaling relation A(r_{c},z)=0.8, D^{0.76}(z), sigma (r_{c})$.

preprint2010arXiv

Bound and unbound substructures in Galaxy-scale Dark Matter haloes

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.

preprint2010arXiv

Cross-Correlation between Damped Ly-alpha Systems and Ly-break Galaxies in Cosmological SPH Simulations

We calculate the cross-correlation function (CCF) between damped Ly-a systems (DLAs) and Lyman break galaxies (LBGs) using cosmological hydrodynamic simulations at z=3. We compute the CCF with two different methods. First, we assume that there is one DLA in each dark matter halo if its DLA cross section is non-zero. In our second approach we weight the pair-count by the DLA cross section of each halo, yielding a cross-section-weighted CCF. We also compute the angular CCF for direct comparison with observations. Finally, we calculate the auto-correlation functions of LBGs and DLAs, and their bias against the dark matter distribution. For these different approaches, we consistently find that there is good agreement between our simulations and observational measurements by Cooke et al. and Adelberger et al. Our results thus confirm that the spatial distribution of LBGs and DLAs can be well described within the framework of the concordance Lambda CDM model. We find that the correlation strengths of LBGs and DLAs are consistent with the actual observations, and in the case of LBGs it is higher than would be predicted by low-mass galaxy merger models.

preprint2010arXiv

Detecting neutral hydrogen in emission at redshift z ~ 1

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.

preprint2010arXiv

Extragalactic gamma-ray background radiation from dark matter annihilation

If dark matter is composed of neutralinos, one of the most exciting prospects for its detection lies in observations of the gamma-ray radiation created in pair annihilations between neutralinos, a process that may contribute significantly to the extragalactic gamma-ray background (EGB) radiation. We here use the high-resolution Millennium-II simulation of cosmic structure formation to produce the first full-sky maps of the expected radiation coming from extragalactic dark matter structures. Our map making procedure takes into account the total gamma-ray luminosity from all haloes and their subhaloes, and includes corrections for unresolved components of the emission as well as an extrapolation to the damping scale limit of neutralinos. Our analysis also includes 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 and may be detectable in forthcoming observations by the recently launched FERMI satellite. 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 difference maps at different energies can reveal cosmic large-scale structure at low and intermediate redshifts. If the intrinsic emission spectrum is characterized by an emission peak, cosmological tomography with gamma ray annihilation radiation is in principle possible.

preprint2010arXiv

Feedback and the Structure of Simulated Galaxies at redshift z=2

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.

preprint2010arXiv

Gas expulsion by quasar-driven winds as a solution to the over-cooling problem in galaxy groups and clusters

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').

preprint2010arXiv

Genus statistics using the Delaunay tessellation field estimation method: (I) tests with the Millennium Simulation and the SDSS DR7

We study the topology of cosmic large-scale structure through the genus statistics, using galaxy catalogues generated from the Millennium Simulation and observational data from the latest Sloan Digital Sky Survey Data Release (SDSS DR7). We introduce a new method for constructing galaxy density fields and for measuring the genus statistics of its isodensity surfaces. It is based on a Delaunay tessellation field estimation (DTFE) technique that allows the definition of a piece-wise continuous density field and the exact computation of the topology of its polygonal isodensity contours, without introducing any free numerical parameter. Besides this new approach, we also employ the traditional approaches of smoothing the galaxy distribution with a Gaussian of fixed width, or by adaptively smoothing with a kernel that encloses a constant number of neighboring galaxies. Our results show that the Delaunay-based method extracts the largest amount of topological information. Unlike the traditional approach for genus statistics, it is able to discriminate between the different theoretical galaxy catalogues analyzed here, both in real space and in redshift space, even though they are based on the same underlying simulation model. In particular, the DTFE approach detects with high confidence a discrepancy of one of the semi-analytic models studied here compared with the SDSS data, while the other models are found to be consistent.

preprint2010arXiv

Hydrodynamical N-body simulations of coupled dark energy cosmologies

If the accelerated expansion of the Universe at the present epoch is driven by a dark energy scalar field, there may well be a non-trivial coupling between the dark energy and the cold dark matter (CDM) fluid. Such interactions give rise to new features in cosmological structure growth, like an additional long-range attractive force between CDM particles, or variations of the dark matter particle mass with time. We have implemented these effects in the N-body code GADGET-2 and present results of a series of high-resolution N-body simulations where the dark energy component is directly interacting with the cold dark matter. As a consequence of the new physics, CDM and baryon distributions evolve differently both in the linear and in the nonlinear regime of structure formation. Already on large scales a linear bias develops between these two components, which is further enhanced by the nonlinear evolution. We also find, in contrast with previous work, that the density profiles of CDM halos are less concentrated in coupled dark energy cosmologies compared with LCDM, and that this feature does not depend on the initial conditions setup, but is a specific consequence of the extra physics induced by the coupling. Also, the baryon fraction in halos in the coupled models is significantly reduced below the universal baryon fraction. These features alleviate tensions between observations and the LCDM model on small scales. Our methodology is ideally suited to explore the predictions of coupled dark energy models in the fully non-linear regime, which can provide powerful constraints for the viable parameter space of such scenarios.

preprint2010arXiv

Intracluster stars in simulations with AGN feedback

We use a set of high-resolution hydrodynamical simulations of clusters of galaxies to study the build-up of the intracluster light (ICL), an interesting and likely significant component of their total stellar mass. Our sample of groups and clusters includes AGN feedback and is of high enough resolution to accurately resolve galaxy populations down to the smallest galaxies that are expected to significantly contribute to the stellar mass budget. We describe and test four different methods to identify the ICL in simulations, thereby allowing us to assess the reliability of the measurements. For all of the methods, we consistently find a very significant ICL stellar fraction (~45%) which exceeds the values typically inferred from observations. However, we show that this result is robust with respect to numerical resolution and integration accuracy, remarkably insensitive to changes in the star formation model, and almost independent of halo mass. It is also almost invariant when black hole growth is included, even though AGN feedback successfully prevents excessive overcooling in clusters and leads to a drastically improved agreement of the simulated cluster galaxy population with observations. In particular, the luminosities of central galaxies and the ages of their stellar populations are much more realistic when including AGN. In the light of these findings, it appears challenging to construct a simulation model that simultaneously matches the cluster galaxy population and at the same time produces a low ICL component. We find that intracluster stars are preferentially stripped in a cluster's densest region from massive galaxies that fall into the cluster at z>1. Surprisingly, some of the intracluster stars also form in the intracluster medium inside cold gas clouds that are stripped out of infalling galaxies.

preprint2010arXiv

Lyman-alpha Emitters and Lyman-break Galaxies at z=3-6 in Cosmological SPH Simulations

We study the properties of Lyman-alpha emitters (LAEs) and Lyman-break galaxies (LBGs) at z=3-6 using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed Ly-a and rest-frame UV luminosity functions (LFs) of LAEs: (i) the "escape fraction" scenario, in which the "effective" escape fraction (including the IGM attenuation) of Ly-a photons is f_Lya ~0.1 (0.15) at z=3 (6), and (ii) the "stochastic" scenario, in which the fraction of LAEs that are turned on at z=3 (6) is \Cstoc ~0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with a number of different observations suggest that the stochastic scenario is preferred over the escape fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass hosted by LAEs are all smaller in the stochastic scenario than in the escape fraction scenario. In our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure formation. However, its evolution is largely compensated by a beginning decline in the specific star formation rate, resulting in little evolution of the rest-frame UV LF from z=6 to 3. The rest-frame UV LF of both LAEs and LBGs at z=3 & 6 can be described well by the stochastic scenario provided the extinction is moderate, E(B-V) ~0.15, for both populations, although our simulation might be overpredicting the number of bright LBGs at z=6. We also discuss the correlation function and bias of LAEs. The Ly-a LFs at z=6 in a field-of-view of 0.2 deg^2 show a significantly larger scatter owing to cosmic variance relative to that in a 1 deg^2 field, and the scatter seen in the current observational estimates of the Ly-a LF can be accounted for by cosmic variance.

preprint2010arXiv

Mapping extragalactic dark matter structures through gamma-rays

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.

preprint2010arXiv

Quasar Clustering in Cosmological Hydrodynamic Simulations: Evidence for mergers

We examine the clustering properties of a population of quasars drawn from fully hydrodynamic cosmological simulations that directly follow black hole growth. We find that the black hole correlation function is best described by two distinct components: contributions from BH pairs occupying the same dark matter halo ('1-halo term') which dominate at scales below 300 kpc/h, and contributions from BHs occupying separate halos ('2-halo term') which dominate at larger scales. From the 2-halo BH term we find a typical host halo mass for faint-end quasars (those probed in our simulation volumes) ranging from 10^11 to a few 10^12 solar masses from z=5 to z=1 respectively (consistent with the mean halo host mass). The BH correlation function shows a luminosity dependence as a function of redshift, though weak enough to be consistent with observational constraints. At small scales, the high resolution of our simulations allows us to probe the 1-halo clustering in detail, finding that the 1-halo term follows an approximate power law, lacking the characteristic decrease in slope at small scales found in 1-halo terms for galaxies and dark matter. We show that this difference is a direct result of a boost in the small-scale quasar bias caused by galaxies hosting multiple quasars (1-subhalo term) following a merger event, typically between a large central subgroup and a smaller, satellite subgroup hosting a relatively small black hole. We show that our predicted small-scale excess caused by such mergers is in good agreement with both the slope and amplitude indicated by recent small-scale measurements. Finally, we note the excess to be a strong function of halo mass, such that the observed excess is well matched by the multiple black holes of intermediate mass (10^7-10^8 solar masses) found in hosts of 4-8*10^11 solar masses, a range well probed by our simulations.

preprint2010arXiv

Secondary Infall and the Pseudo-Phase-Space Density Profiles of Cold Dark Matter Halos

We use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual particles, Q_i; profiles based on the ellipsoidal radius dictated by the triaxial structure of the halo, Q_i(r'); and by carefully removing substructures in order to focus on the profile of the smooth halo, Q^s. The resulting Q_i^s(r') profiles follow closely a power law near the center, but show a clear upturn from this trend near the virial radius, r_{200}. The location and magnitude of the deviations are in excellent agreement with the predictions from Bertschinger's spherical secondary-infall similarity solution. In this model, Q \propto r^{-1.875} in the inner, virialized regions, but departures from a power-law occur near r_{200} because of the proximity of this radius to the location of the first shell crossing - the shock radius in the case of a collisional fluid. Particles there have not yet fully virialized, and so Q departs from the inner power-law profile. Our results imply that the power-law nature of $Q$ profiles only applies to the inner regions and cannot be used to predict accurately the structure of CDM halos beyond their characteristic scale radius.

preprint2010arXiv

The case for AGN feedback in galaxy groups

[Abridged] The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine its effects on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the universe. To explore the effects of BH feedback on groups, we analyse two high resolution cosmological hydro simulations from the OverWhelmingly Large Simulations project. While both include galactic winds driven by supernovae, only one includes feedback from BHs. We compare the properties of the simulated groups to a wide range of observational data, including hot gas radial profiles and gas mass fractions (fgas), luminosity-mass-temperature (L-M-T) scaling relations, K-band luminosity of the group and its central brightest galaxy (CBG), SFRs and ages of the CBG, and gas/stellar metallicities. Both runs yield entropy profiles similar to the data, while the run without AGN feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from BHs significantly reduces fgas for groups with masses less than ~10^14 Msun, yielding fgas-T and L-T relations that are in agreement with the data. The run without AGN feedback suffers from the well known overcooling problem; the resulting K-band luminosities are much larger than observed. By contrast, the run that includes BH feedback yields K-band luminosities and CBG SFRs and ages in agreement with current estimates. Both runs yield very similar gas-phase metallicities that match X-ray data, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that BH feedback is a crucial ingredient in the formation of massive galaxies.

preprint2010arXiv

The effect of neutrinos on the matter distribution as probed by the Intergalactic Medium

We present a suite of full hydrodynamical cosmological simulations that quantitatively address the impact of neutrinos on the (mildly non-linear) spatial distribution of matter and in particular on the neutral hydrogen distribution in the Intergalactic Medium (IGM), which is responsible for the intervening Lyman-alpha absorption in quasar spectra. The free-streaming of neutrinos results in a (non-linear) scale-dependent suppression of power spectrum of the total matter distribution at scales probed by Lyman-alpha forest data which is larger than the linear theory prediction by about 25% and strongly redshift dependent. By extracting a set of realistic mock quasar spectra, we quantify the effect of neutrinos on the flux probability distribution function and flux power spectrum. The differences in the matter power spectra translate into a ~2.5% (5%) difference in the flux power spectrum for neutrino masses with Sigma m_ν = 0.3 eV (0.6 eV). This rather small effect is difficult to detect from present Lyman-alpha forest data and nearly perfectly degenerate with the overall amplitude of the matter power spectrum as characterised by sigma_8. If the results of the numerical simulations are normalized to have the same sigma_8 in the initial conditions, then neutrinos produce a smaller suppression in the flux power of about 3% (5%) for Sigma m_ν = 0.6$ eV (1.2 eV) when compared to a simulation without neutrinos. We present constraints on neutrino masses using the Sloan Digital Sky Survey flux power spectrum alone and find an upper limit of Sigma m_ν < 0.9$ eV (2 sigma C.L.), comparable to constraints obtained from the cosmic microwave background data or other large scale structure probes.

preprint2010arXiv

The Stellar Mass Components of Galaxies: Comparing Semi-Analytical Models with Observation

We compare the stellar masses of central and satellite galaxies predicted by three independent semianalytical models with observational results obtained from a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. In particular, we compare the stellar mass functions of centrals and satellites, the relation between total stellar mass and halo mass, and the conditional stellar mass functions, which specify the average number of galaxies of stellar mass M_* that reside in a halo of mass M_h. The semi-analytical models only predict the correct stellar masses of central galaxies within a limited mass range and all models fail to reproduce the sharp decline of stellar mass with decreasing halo mass observed at the low mass end. In addition, all models over-predict the number of satellite galaxies by roughly a factor of two. The predicted stellar mass in satellite galaxies can be made to match the data by assuming that a significant fraction of satellite galaxies are tidally stripped and disrupted, giving rise to a population of intra-cluster stars in their host halos. However, the amount of intra-cluster stars thus predicted is too large compared to observation. This suggests that current galaxy formation models still have serious problems in modeling star formation in low-mass halos.

preprint2010arXiv

There's no place like home? Statistics of Milky Way-mass dark matter halos

We present an analysis of the distribution of structural properties for Milky Way-mass halos in the Millennium-II Simulation (MS-II). This simulation of structure formation within the standard LCDM cosmology contains thousands of Milky Way-mass halos and has sufficient resolution to properly resolve many subhalos per host. It thus provides a major improvement in the statistical power available to explore the distribution of internal structure for halos of this mass. In addition, the MS-II contains lower resolution versions of the Aquarius Project halos, allowing us to compare our results to simulations of six halos at a much higher resolution. We study the distributions of mass assembly histories, of subhalo mass functions and accretion times, and of merger and stripping histories for subhalos capable of impacting disks at the centers of halos. We show that subhalo abundances are not well-described by Poisson statistics at low mass, but rather are dominated by intrinsic scatter. Using the masses of subhalos at infall and the abundance-matching assumption, there is less than a 10% chance that a Milky Way halo with M_vir =10^12 M_sun will host two galaxies as bright as the Magellanic Clouds. This probability rises to ~25% for a halo with M_vir=2.5 x 10^12 M_sun. The statistics relevant for disk heating are very sensitive to the mass range that is considered relevant. Mergers with infall mass : redshift zero virial mass greater than 1:30 could well impact a central galactic disk and are a near inevitability since z=2, whereas only half of all halos have had a merger with infall mass : redshift zero virial mass greater than 1:10 over this same period.

preprint2009arXiv

Caustics in growing Cold Dark Matter Haloes

We simulate the growth of isolated dark matter haloes from self-similar and spherically symmetric initial conditions. Our N-body code integrates the geodesic deviation equation in order to track the streams and caustics associated with individual simulation particles. The radial orbit instability causes our haloes to develop major-to-minor axis ratios approaching 10 to 1 in their inner regions. They grow similarly in time and have similar density profiles to the spherical similarity solution, but their detailed structure is very different. The higher dimensionality of the orbits causes their stream and caustic densities to drop much more rapidly than in the similarity solution. This results in a corresponding increase in the number of streams at each point. At 1% of the turnaround radius (corresponding roughly to the Sun's position in the Milky Way) we find of order 10^6 streams in our simulations, as compared to 10^2 in the similarity solution. The number of caustics in the inner halo increases by a factor of several, because a typical orbit has six turning points rather than one, but caustic densities drop by a much larger factor. This reduces the caustic contribution to the annihilation radiation. For the region between 1% and 50% of the turnaround radius, this is 4% of the total in our simulated haloes, as compared to 6.5% in the similarity solution. Caustics contribute much less at smaller radii. These numbers assume a 100 GeV c^-2 neutralino with present-day velocity dispersion 0.03 cm s^-1, but reducing the dispersion by ten orders of magnitude only doubles the caustic luminosity. We conclude that caustics will be unobservable in the inner parts of haloes. Only the outermost caustic might potentially be detectable.

preprint2009arXiv

E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh

Hydrodynamic cosmological simulations at present usually employ either the Lagrangian SPH technique, or Eulerian hydrodynamics on a Cartesian mesh with adaptive mesh refinement. Both of these methods have disadvantages that negatively impact their accuracy in certain situations. We here propose a novel scheme which largely 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 unsplit Godunov scheme with an exact Riemann solver. The mesh-generating points can in principle be moved arbitrarily. If they are chosen to be stationary, the scheme is equivalent to an ordinary Eulerian method with second order accuracy. If they instead move with the velocity of the local flow, one obtains a Lagrangian formulation of hydrodynamics that does not suffer from the mesh distortion limitations inherent in other mesh-based Lagrangian schemes. In this mode, our new method is fully Galilean-invariant, unlike ordinary Eulerian codes, a property that is of significant importance for cosmological simulations. In addition, the new scheme can adjust its spatial resolution automatically and continuously, and hence inherits the principal advantage of SPH for simulations of cosmological structure growth. The high accuracy of Eulerian methods in the treatment of shocks is retained, while the treatment of contact discontinuities improves. We discuss how this approach is implemented in our new parallel code AREPO, both in 2D and 3D. We use a suite of test problems to examine the performance of the new code and argue that it provides an attractive and competitive alternative to current SPH and Eulerian techniques. (abridged)

preprint2009arXiv

Faint-end Quasar Luminosity Functions from Cosmological Hydrodynamic Simulations

We investigate the predictions for the faint-end quasar luminosity function (QLF) and its evolution using fully cosmological hydrodynamic simulations which self-consistently follow star formation, black hole growth and associated feedback processes. We find remarkably good agreement between predicted and observed faint end of the optical and X-ray QLFs (the bright end is not accessible in our simulated volumes) at z < 2. At higher redshifts our simulations tend to overestimate the QLF at the faintest luminosities. We show that although the low (high) luminosity ranges of the faint-end QLF are dominated by low (high) mass black holes, a wide range of black hole masses still contributes to any given luminosity range. This is consistent with the complex lightcurves of black holes resulting from the detailed hydrodynamics followed in the simulations. Consistent with the results on the QLFs, we find good agreement for the evolution of the comoving number density (in optical, soft and hard X-ray bands) of AGN for luminosities above 10^43 erg/s. However, the luminosity density evolution from the simulation appears to imply a peak at higher redshift than constrained from hard X-ray data (but not in optical). Our predicted excess at the faintest fluxes at z >= 2 does not lead to an overestimate to the total X-ray background and its contribution is at most a factor of two larger than the unresolved fraction of the 2-8 keV background. Even though this could be explained by some yet undetected, perhaps heavily obscured faint quasar population, we show that our predictions for the faint sources at high redshifts (which are dominated by the low mass black holes) in the simulations are likely affected by resolution effects.

preprint2009arXiv

Galaxies-Intergalactic Medium Interaction Calculation --I. Galaxy formation as a function of large-scale environment

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

preprint2009arXiv

Growing the first bright quasars in cosmological simulations of structure formation

We employ cosmological hydrodynamical simulations to study the growth of massive black holes (BHs) at high redshifts subject to BH merger recoils from gravitational wave emission. We select the most massive dark matter halo at z=6 from the Millennium simulation, and resimulate its formation at much higher resolution including gas physics and a model for BH seeding, growth and feedback. Assuming that the initial BH seeds are relatively massive, of the order of 10^5 Msun, and that seeding occurs around z~15 in dark matter haloes of mass 10^9-10^10 Msun, we find that it is possible to build up supermassive BHs (SMBHs) by z=6 that assemble most of their mass during extended Eddington-limited accretion periods. The properties of the simulated SMBHs are consistent with observations of z=6 quasars in terms of the estimated BH masses and bolometric luminosities, the amount of star formation occurring within the host halo, and the presence of highly enriched gas in the innermost regions of the host galaxy. After a peak in the BH accretion rate at z=6, the most massive BH has become sufficiently massive for the growth to enter into a much slower phase of feedback-regulated accretion. We explore the full range of expected recoils and radiative efficiencies, and also consider models with spinning BHs. In the most `pessimistic' case where BH spins are initially high, we find that the growth of the SMBHs can be potentially hampered if they grow mostly in isolation and experience only a small number of mergers. Whereas BH kicks can expel a substantial fraction of low mass BHs, they do not significantly affect the build up of the SMBHs. On the contrary, a large number of BH mergers has beneficial consequences for the growth of the SMBHs by considerably reducing their spin. [Abridged]

preprint2009arXiv

Lensed CMB temperature and polarization maps from the Millennium Simulation

We have constructed the first all-sky CMB temperature and polarization lensed maps based on a high-resolution cosmological N-body simulation, the Millennium Simulation (MS). We have exploited the lensing potential map obtained using a map-making procedure (Carbone et al. 2008) which integrates along the line-of-sight the MS dark matter distribution by stacking and randomizing the simulation boxes up to $z = 127$, and which semi-analytically supplies the large-scale power in the angular lensing potential that is not correctly sampled by the N-body simulation. The lensed sky has been obtained by properly modifying the latest version of the LensPix code (Lewis 2005) to account for the MS structures. We have also produced all-sky lensed maps of the so-called $ψ_E$ and $ψ_B$ potentials, which are directly related to the electric and magnetic types of polarization. The angular power spectra of the simulated lensed temperature and polarization maps agree well with semi-analytic estimates up to $l < 2500$, while on smaller scales we find a slight excess of power which we interpret as being due to non-linear clustering in the MS. We also observe how non-linear lensing power in the polarised CMB is transferred to large angular scales by suitably misaligned modes in the CMB and the lensing potential. This work is relevant in view of the future CMB probes, as a way to analyse the lensed sky and disentangle the contribution from primordial gravitational waves.

preprint2009arXiv

Modeling the cosmological co-evolution of supermassive black holes and galaxies: II. The clustering of quasars and their dark environment

We use semi-analytic modeling on top of the Millennium simulation to study the joint formation of galaxies and their embedded supermassive black holes. Our goal is to test scenarios in which black hole accretion and quasar activity are triggered by galaxy mergers, and to constrain different models for the lightcurves associated with individual quasar events. In the present work we focus on studying the spatial distribution of simulated quasars. At all luminosities, we find that the simulated quasar two-point correlation function is fit well by a single power-law in the range 0.5 < r < 20 h^{-1} Mpc, but its normalization is a strong function of redshift. When we select only quasars with luminosities within the range typically accessible by today's quasar surveys, their clustering strength depends only weakly on luminosity, in agreement with observations. This holds independently of the assumed lightcurve model, since bright quasars are black holes accreting close to the Eddington limit, and are hosted by dark matter haloes with a narrow mass range of a few 10^12 h^{-1} M_sun. Therefore the clustering of bright quasars cannot be used to disentangle lightcurve models, but such a discrimination would become possible if the observational samples can be pushed to significantly fainter limits. Overall, our clustering results for the simulated quasar population agree rather well with observations, lending support to the conjecture that galaxy mergers could be the main physical process responsible for triggering black hole accretion and quasar activity.

preprint2009arXiv

On merger bias and the clustering of quasars

We use the large catalogues of haloes available for the Millennium Simulation to test whether recently merged haloes exhibit stronger large-scale clustering than other haloes of the same mass. This effect could help to understand the very strong clustering of quasars at high redshift. However, we find no statistically significant excess bias for recently merged haloes over the redshift range 2 < z < 5, with the most massive haloes showing an excess of at most ~5%. We also consider galaxies extracted from a semianalytic model built on the Millennium Simulation. At fixed stellar mass, we find an excess bias of ~ 20-30% for recently merged objects, decreasing with increasing stellar mass. The fact that recently-merged galaxies are found in systematically more massive haloes than other galaxies of the same stellar mass accounts for about half of this signal, and perhaps more for high-mass galaxies. The weak merger bias of massive systems suggests that objects of merger-driven nature, such as quasars, do not cluster significantly differently than other objects of the same characteristic mass. We discuss the implications of these results for the interpretation of clustering data with respect to quasar duty cycles, visibility times, and evolution in the black hole-host mass relation.

preprint2009arXiv

Simulations of AGN feedback in galaxy clusters and groups

There is compelling evidence that black holes (BHs) in cluster centers vigorously interact with their surroundings, indicating that any realistic model of cluster formation needs to account for these processes. Here we use high-resolution cosmological simulations of a large cluster and group sample to study how BHs affect their host systems. We focus on two specific properties, the halo gas fraction and the X-ray luminosity-temperature scaling relation, both of which are notoriously difficult to reproduce in self-consistent hydrodynamical simulations. We show that BH feedback can solve both of these issues, bringing them in excellent agreement with observations, without alluding to the `cooling only' solution that produces unphysically bright central galaxies.

preprint2009arXiv

Substructure depletion in the Milky Way halo by the disk

We employ numerical simulations and simple analytical estimates to argue that dark matter substructures orbiting in the inner regions of the Galaxy can be efficiently destroyed by disk shocking, a dynamical process known to affect globular star clusters. We carry out a set of fiducial high-resolution collisionless simulations in which we adiabatically grow a disk, allowing us to examine the impact of the disk on the substructure abundance. We also track the orbits of dark matter satellites in the high-resolution Aquarius simulations and analytically estimate the cumulative halo and disk shocking effect. Our calculations indicate that the presence of a disk with only 10% of the total Milky Way mass can significantly alter the mass function of substructures in the inner parts of halos. This has important implications especially for the relatively small number of satellites seen within ~30 kpc of the Milky Way center, where disk shocking is expected to reduce the substructure abundance by a factor of ~2 at 10^9 M$_{\odot}$ and ~3 at 10^7 M$_{\odot}$. The most massive subhalos with 10^10 M$_{\odot}$ survive even in the presence of the disk. This suggests that there is no inner missing satellite problem, and calls into question whether these substructures can produce transient features in disks, like multi-armed spiral patterns. Also, the depletion of dark matter substructures through shocking on the baryonic structures of the disk and central bulge may aggravate the problem to fully account for the observed flux anomalies in gravitational lens systems, and significantly reduces the dark matter annihilation signal expected from nearby substructures in the inner halo.

preprint2009arXiv

The Diversity and Similarity of Simulated Cold Dark Matter Halos

We study the mass, velocity dispersion, and anisotropy profiles of $Λ$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $γ\equiv -$d$\lnρ/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $ρ\propto r^{-1.2}$ are clearly ruled out. The radial dependence of $γ$ is well approximated by a power-law, $γ\propto r^α$ (the Einasto profile). The shape parameter, $α$, varies slightly but significantly from halo to halo, implying that the mass profiles of $Λ$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $ρ/σ^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

preprint2009arXiv

The Formation and Survival of Discs in a Lambda-CDM Universe

We study the formation of galaxies in a Lambda-CDM Universe using high resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at z = 0, none of our galaxies contain a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The z = 0 spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at z = 0 nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at z > ~2, regardless of their z = 0 morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic Lambda-CDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.

preprint2009arXiv

The impact of feedback on the low redshift Intergalactic Medium

We analyse the evolution of the properties of the low-redshift Intergalactic Medium (IGM) using high-resolution hydrodynamic simulations that include a detailed chemical evolution model. We focus on the effects that two different forms of energy feedback, strong galactic winds driven by supernova explosion and Active Galactic Nuclei (AGN) powered by gas accretion onto super-massive black holes (BHs), have on the thermo- and chemo-dynamical properties of of the low redshift IGM. We find that feedback associated to winds (W) and BHs leave distinct signatures in both the chemical and thermal history of the baryons, especially at redshift z<3 [..] We present results for the enrichment in terms of mass and metallicity distributions for the WHIM phase, both as a function of density and temperature. Finally, we compute the evolution of the relative abundances between different heavy elements, namely Oxygen, Carbon and Iron. While both C/O and O/Fe evolve differently at high redshifts for different feedback models, their values are similar at z=0 [..]. The sensitivity of WHIM properties on the implemented feedback scheme could be important both for discriminating between different feedback physics and for detecting the WHIM with future far-UV and X-ray telescopes.

preprint2009arXiv

The role of dry mergers for the formation and evolution of brightest cluster galaxies

Using a resimulation technique, we perform high-resolution cosmological simulations of dry mergers in a massive galaxy cluster identified in the Millennium Run. Our initial conditions include well resolved compound galaxy models consisting of dark matter halos and stellar bulges that are used to replace the most massive cluster progenitor halos at redshift z=3. By construction, our galaxy models obey the stellar mass-size relation initially. We demonstrate that the brightest cluster galaxy (BCG) evolves away from the Kormendy relation as defined by the smaller mass galaxies (i.e., the relation bends). As a result of the comparatively large number of mergers the BCG experiences, its total mass-to-light ratio becomes significantly higher than in typical elliptical galaxies. We also show that the mixing processes between dark matter and stars lead to a small but numerically robust tilt in the fundamental plane and that the BCG lies on the tilted plane. Our model is consistent with the observed steepening of the logarithmic mass-to-light gradient as a function of the stellar mass. Surprisingly, we find only tentative weak distortion in the Faber-Jackson relation that depends on the aperture size, unlike expected based on studies of isolated merger simulations. This may be due to differences in the distribution of galaxy orbits, which is given in our approach directly by the cosmological context while it has to be assumed in isolated merger simulations, and the fact that the BCG is located deep in the cluster potential well.

preprint2009arXiv

The spatial distribution of X-ray selected AGN in the Chandra deep fields: a theoretical perspective

We study the spatial distribution of X-ray selected AGN in the framework of hierarchical co-evolution of supermassive black holes and their host galaxies and dark matter haloes. To this end, we have applied the model developed by Croton et al.(2006), De Lucia & Blaizot(2007) and Marulli et al.(2008) to the output of the Millennium Run and obtained hundreds of realizations of past light-cones from which we have extracted realistic mock AGN catalogues that mimic the Chandra deep fields. We find that the model AGN number counts are in fair agreement with observations, except at fluxes <1e-15 erg/cm^2/s. The spatial two-point correlation function predicted by the model is well described by a power-law relation out to 20 Mpc/h, in close agreement with observations. Our model matches the correlation length r_0 of AGN in the Chandra Deep Field North but underestimates it in the Chandra Deep Field South. When fixing the slope to gamma = 1.4, as in Gilli et al. (2005), the statistical significance of the mismatch is 2-2.5 sigma, suggesting that the predicted cosmic variance, which dominates the error budget, may not account for the different correlation length of the AGN in the two fields. While our results are robust to changes in the model prescriptions for the AGN lightcurves, the luminosity dependence of the clustering is sensitive to the different lightcurve models adopted. However, irrespective of the model considered, the luminosity dependence of the AGN clustering in our mock fields seems to be weaker than in the real Chandra fields. The significance of this mismatch needs to be confirmed using larger datasets.

preprint2008arXiv

Effects of Supernova Feedback on the Formation of Galaxies

We study the effects of Supernova (SN) feedback on the formation of galaxies using hydrodynamical simulations in a Lambda-CDM cosmology. We use an extended version of the code GADGET-2 which includes chemical enrichment and energy feedback by Type II and Type Ia SN, metal-dependent cooling and a multiphase model for the gas component. We focus on the effects of SN feedback on the star formation process, galaxy morphology, evolution of the specific angular momentum and chemical properties. We find that SN feedback plays a fundamental role in galaxy evolution, producing a self-regulated cycle for star formation, preventing the early consumption of gas and allowing disks to form at late times. The SN feedback model is able to reproduce the expected dependence on virial mass, with less massive systems being more strongly affected.

preprint2008arXiv

Galaxy morphology, kinematics and clustering in a hydrodynamic simulation of a LambdaCDM universe

We explore galaxy properties and their link with environment and clustering using a population of ~1000 galaxies formed in a high resolution hydrodynamic simulation of the Lambda CDM cosmology. At the redshift we concentrate on, z=1, the spatial resolution is 1.4 proper kpc/h and Milky-way sized disk galaxies contain ~10^5 particles within their virial radii. We include supermassive black hole accretion and feedback as well as a multiphase model for star formation. We find that a number of familiar qualitative relationships hold approximately between galaxy properties, for example, galaxies lie between two broad extremes of type, where ``late'' types tend to be smaller in size, have lower circular velocities, younger stars, higher star formation rates, larger disk to bulge ratios and lower Sersic indices than ``early types''. As in previous studies the stellar component of disk galaxies is not as rotationally supported as in observations. Bulges contain too much of the stellar mass, although disks do have scale lengths compatible with observations. The addition of black hole physics to the simulations does not appear to have an impact on the angular momentum results, nor do we find that it is affected in an identical simulation with significantly lower mass resolution. Despite this, we can profitably use the rank order of either disk to total ratio, Sersic index, or galaxy age to separate galaxies into morphological classes and examine the density-morphology relation and morphology dependence of clustering. We find that while at redshift z=0, the well known preponderance of early types in dense environments is seen, at z=1 the density-morphology relation becomes flatter and late type galaxies are even seen to have a higher clustering amplitude than early types (abridged).

preprint2008arXiv

Quantifying the Cosmic Web I: The large-scale halo ellipticity-ellipticity and ellipticity-direction correlations

The formation of dark matter halos tends to occur anisotropically along the filaments of the Cosmic Web, which induces both ellipticity-ellipticity (EE) correlations between the shapes of halos, as well as ellipticity-direction (ED) cross-correlations between halo shapes and the directions to neighboring halos. We analyze the halo catalogue and the semi-analytic galaxy catalogue of the recent Millennium Run Simulation to measure the EE and ED correlations numerically at four different redshifts (z=0, 0.5, 1 and 2). For the EE correlations, we find that (i) the major-axis correlation is strongest while the intermediate-axis correlation is weakest; (ii) the signal is significant at distances out to 10 Mpc/h; (iii) the signal decreases as z decreases; (iv) and its behavior depends strongly on the halo mass scale, with larger masses showing stronger correlations at large distances. For the ED correlations, we find that (i) the correlations are much stronger than the EE correlations, and are significant even out to distances of 50 Mpc/h; (ii) the signal also decreases as z decreases; (iii) and it increases with halo mass at all distances. We also provide empirical fitting functions for the EE and ED correlations. The EE correlations are found to scale linearly with the linear density correlation function, xi(r). While the ED cross-correlation is found to scale as xi^{1/2}(r) at large distances beyond 10 Mpc/h. The best-fit values of the fitting parameters for the EE and the ED correlations are all determined through chi^{2}-statistics. Our results may be useful for quantifying the filamentary distribution of dark matter halos over a wide range of scales.

preprint2008arXiv

The Aquarius Project: the subhalos of galactic halos

We have performed the largest ever particle simulation of a Milky Way-sized dark matter halo, and present the most comprehensive convergence study for an individual dark matter halo carried out thus far. We have also simulated a sample of 6 ultra-highly resolved Milky-way sized halos, allowing us to estimate the halo-to-halo scatter in substructure statistics. In our largest simulation, we resolve nearly 300,000 gravitationally bound subhalos within the virialized region of the halo. Simulations of the same object differing in mass resolution by factors up to 1800 accurately reproduce the largest subhalos with the same mass, maximum circular velocity and position, and yield good convergence for the abundance and internal properties of dark matter substructures. We detect up to four generations of subhalos within subhalos, but contrary to recent claims, we find less substructure in subhalos than in the main halo when regions of equal mean overdensity are compared. The overall substructure mass fraction is much lower in subhalos than in the main halo. Extrapolating the main halo's subhalo mass spectrum down to an Earth mass, we predict the mass fraction in substructure to be well below 3% within 100 kpc, and to be below 0.1% within the Solar Circle. The inner density profiles of subhalos show no sign of converging to a fixed asymptotic slope and are well fit by gently curving profiles of Einasto form. The mean concentrations of isolated halos are accurately described by the fitting formula of Neto et al. down to maximum circular velocities of 1.5 km/s, an extrapolation over some 5 orders of magnitude in mass. However, at equal maximum circular velocity, subhalos are more concentrated than field halos, with a characteristic density that is typically ~2.6 times larger and increases towards the halo centre.

preprint2008arXiv

The impact of Early Dark Energy on non-linear structure formation

We study non-linear structure formation in high-resolution simulations of Early Dark Energy (EDE) cosmologies and compare their evolution with the standard LCDM model. Extensions of the spherical top-hat collapse model predict that the virial overdensity and linear threshold density for collapse should be modified in EDE model, yielding significant modifications in the expected halo mass function. Here we present numerical simulations that directly test these expectations. Interestingly, we find that the Sheth & Tormen formalism for estimating the abundance of dark matter halos continues to work very well in its standard form for the Early Dark Energy cosmologies, contrary to analytic predictions. The residuals are even slightly smaller than for LCDM. We also study the virial relationship between mass and dark matter velocity dispersion in different dark energy cosmologies, finding excellent agreement with the normalization for Lambda as calibrated by Evrard et al.(2008). The earlier growth of structure in EDE models relative to LCDM produces large differences in the mass functions at high redshift. This could be measured directly by counting groups as a function of the line-of-sight velocity dispersion, skirting the ambiguous problem of assigning a mass to the halo. Using dark matter substructures as a proxy for member galaxies, we demonstrate that even with 3-5 members sufficiently accurate measurements of the halo velocity dispersion function are possible. Finally, we determine the concentration-mass relationship for our EDE cosmologies. Consistent with the earlier formation time, the EDE halos show higher concentrations at a given halo mass. We find that the magnitude of the difference in concentration is well described by the prescription of Eke et al.(2001) for estimating halo concentrations.

preprint2007arXiv

The statistics of LCDM Halo Concentrations

We use the Millennium Simulation (MS) to study the statistics of LCDM halo concentrations at z = 0. Our results confirm that the average halo concentration declines monotonically with mass; a power-law fits well the concentration-mass relation for over 3 decades in mass, up to the most massive objects to form in a LCDM universe (~ 10^15 h^-1 Msol). This is in clear disagreement with the predictions of the model proposed by Bullock et al. for these rare objects, and agrees better with the original predictions of Navarro, Frenk, & White. The large volume surveyed, together with the unprecedented numerical resolution of the MS, allow us to estimate with confidence the distribution of concentrations and, consequently, the abundance of systems with unusual properties. About one in a hundred cluster haloes (M200 >~ 3x10^14 h^-1 Msol) have concentrations exceeding c200 = 7.5, a result that may be used to interpret the likelihood of unusually strong massive gravitational lenses, such as Abell 1689, in the LCDM cosmogony. A similar fraction (1 in 100) of galaxy-sized haloes (M200 ~ 10^12 h^-1 Msol) have c200 < 4.5, an important constraint on models that attempt to reconcile the rotation curves of low surface-brightness galaxies by appealing to haloes of unexpectedly low concentration. We find that halo concentrations are independent of spin once haloes manifestly out of equilibrium are removed from the sample. Compared to their relaxed brethren, the concentrations of out-of-equilibrium haloes tend to be lower and to have more scatter, while their spins tend to be higher... (continue)

preprint2006arXiv

Cosmological structure formation shocks and cosmic rays in hydrodynamical simulations

Cosmological shock waves during structure formation not only play a decisive role for the thermalization of gas in virializing structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. We discuss a novel numerical treatment of the physics of cosmic rays in combination with a formalism for identifying and measuring the shock strength on-the-fly during a smoothed particle hydrodynamics simulation. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. Using this formalism, we study the history of the thermalization process in high-resolution hydrodynamic simulations of the Lambda cold dark matter model. Collapsed cosmological structures are surrounded by shocks with high Mach numbers up to 1000, but they play only a minor role in the energy balance of thermalization. However, this finding has important consequences for our understanding of the spatial distribution of CRs in the large-scale structure. In high resolution simulations of galaxy clusters, we find a low contribution of the averaged CR pressure, due to the small acceleration efficiency of lower Mach numbers of flow shocks inside halos and the softer adiabatic index of CRs. However, within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading there to a lower effective adiabatic index and thus to an enhanced compressibility of the central intracluster medium. This effect increases the central density and pressure of the cluster and thus the resulting X-ray emission and the central Sunyaev-Zel'dovich flux decrement. The integrated Sunyaev-Zel'dovich effect, however, is only slightly changed.

preprint2006arXiv

The large-scale structure of the Universe

Research over the past 25 years has led to the view that the rich tapestry of present-day cosmic structure arose during the first instants of creation, where weak ripples were imposed on the otherwise uniform and rapidly expanding primordial soup. Over 14 billion years of evolution, these ripples have been amplified to enormous proportions by gravitational forces, producing ever-growing concentrations of dark matter in which ordinary gases cool, condense and fragment to make galaxies. This process can be faithfully mimicked in large computer simulations, and tested by a variety of observations that probe the history of the Universe starting from just 400,000 yr after the Big Bang.

preprint2006arXiv

The role of AGN feedback and gas viscosity in hydrodynamical simulations of galaxy clusters

We study the imprints of AGN feedback and physical viscosity on the properties of galaxy clusters using hydrodynamical simulation models carried out with the TreeSPH code GADGET-2. Besides self-gravity of dark matter and baryons, our approach includes radiative cooling and heating processes of the gas component and a multiphase model for star formation and SNe feedback. Additionally, we introduce a prescription for physical viscosity in GADGET-2, based on a SPH discretization of the Navier-Stokes and general heat transfer equations. Adopting the Braginskii parameterization for the shear viscosity coefficient, we explore how gas viscosity influences the properties of AGN-driven bubbles. We also introduce a novel, self-consistent AGN feedback model where we simultaneously follow the growth and energy release of massive black holes embedded in a cluster environment. We assume that black holes accreting at low rates with respect to the Eddington limit are in a radiatively inefficient regime, and that most of the feedback energy will appear in a mechanical form. Thus, we introduce AGN-driven bubbles into the ICM with properties, such as radius and energy content, that are directly linked to the black hole physics. This model leads to a self-regulated mechanism for the black hole growth and overcomes the cooling flow problem in host halos, ranging from the scale of groups to that of massive clusters. (Abridged)

preprint2005arXiv

Constrained Simulations of the Magnetic Field in the Local Universe and the Propagation of UHECRs

We use simulations of LSS formation to study the build-up of magnetic fields (MFs) in the ICM. Our basic assumption is that cosmological MFs grow in a MHD amplification process driven by structure formation out of a seed MF present at high z. Our LCDM initial conditions for the density fluctuations have been statistically constrained by the observed galaxies, based on the IRAS 1.2-Jy all-sky redshift survey. As a result, prominent galaxy clusters in our simulation coincide closely with their real counterparts. We find excellent agreement between RMs of our simulated clusters and observational data. The improved resolution compared to previous work also allows us to study the MF in large-scale filaments, sheets and voids. By tracing the propagation of UHE protons in the simulated MF we construct full-sky maps of expected deflection angles of protons with arrival energies E=1e20eV and 4e19eV, respectively. Strong deflections are only produced if UHE protons cross clusters, however covering only a small area on the sky. Multiple crossings of sheets and filaments over larger distances may give rise to noticeable deflections, depending on the model adopted for the magnetic seed field. Based on our results we argue that over a large fraction of the sky the deflections are likely to remain smaller than the present experimental angular sensitivity. Therefore, we conclude that forthcoming air shower experiments should be able to locate sources of UHE protons and shed more light on the nature of cosmological MFs.

preprint2005arXiv

The age dependence of halo clustering

We use a very large simulation of the concordance LCDM cosmogony to study the clustering of dark matter haloes. For haloes less massive than about 1e13Msun/h, the amplitude of the two-point correlation function on large scales depends strongly on halo formation time. Haloes that assembled at high redshift are substantially more clustered than those that assembled more recently. The effect is a smooth function of halo formation time and its amplitude increases with decreasing halo mass. At 1e11 Msun/h the ``oldest'' 10% of haloes are more than 5 times more strongly correlated than the ``youngest'' 10%. This unexpected result is incompatible with the standard excursion set theory for structure growth, and it contradicts a fundamental assumption of the halo occupation distribution models often used to study galaxy clustering, namely that the galaxy content of a halo of given mass is statistically independent of its larger scale environment.

preprint2005arXiv

The Evolution of the M-sigma Relation

(Abridged) We examine the evolution of the black hole mass - stellar velocity dispersion (M-sigma) relation over cosmic time using simulations of galaxy mergers that include feedback from supermassive black hole growth. We consider mergers of galaxies varying the properties of the progenitors to match those expected at redshifts z=0-6. We find that the slope of the resulting M-sigma relation is the same at all redshifts considered. For the same feedback efficiency that reproduces the observed amplitude of the M-sigma relation at z=0, there is a weak redshift-dependence to the normalization that results from an increasing velocity dispersion for a given galactic stellar mass. We develop a formalism to connect redshift evolution in the M-sigma relation to the scatter in the local relation at z=0. We show that the scatter in the local relation places severe constraints on the redshift evolution of both the normalization and slope of the M-sigma relation. Furthermore, we demonstrate that cosmic downsizing introduces a black hole mass-dependent dispersion in the M-sigma relation and that the skewness of the distribution about the locally observed M-sigma relation is sensitive to redshift evolution in the normalization and slope. In principle, these various diagnostics provide a method for differentiating between theories for producing the M-sigma relation. In agreement with existing constraints, our simulations imply that hierarchical structure formation should produce the relation with small intrinsic scatter.

preprint2004arXiv

Detecting Sunyaev-Zel'dovich clusters with PLANCK: I. Construction of all-sky thermal and kinetic SZ-maps

All-sky thermal and kinetic Sunyaev-Zel'dovich (SZ) maps are presented for assessing how well the PLANCK-mission can find and characterise clusters of galaxies, especially in the presence of primary anisotropies of the cosmic microwave background (CMB) and various galactic and ecliptic foregrounds. The maps have been constructed from numerical simulations of structure formation in a standard LCDM cosmology and contain all clusters out to redshifts of z = 1.46 with masses exceeding 5e13 M_solar/h. By construction, the maps properly account for the evolution of cosmic structure, the halo-halo correlation function, the evolving mass function, halo substructure and adiabatic gas physics. The velocities in the kinetic map correspond to the actual density environment at the cluster positions. We characterise the SZ-cluster sample by measuring the distribution of angular sizes, the integrated thermal and kinetic Comptonisations, the source counts in the three relevant PLANCK-channels, and give the angular power-spectra of the SZ-sky. While our results are broadly consistent with simple estimates based on scaling relations and spherically symmetric cluster models, some significant differences are seen which may affect the number of cluster detectable by PLANCK.

preprint2004arXiv

Massive galaxies at redshift 2 in cosmological hydrodynamic simulations

We study the properties of galaxies at z=2 in a Lambda cold dark matter universe, using two different types of hydrodynamic simulation methods -- Eulerian TVD and smoothed particle hydrodynamics (SPH) -- and a spectrophotometric analysis in the U_n, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) and Steidel et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 1e-2 h^3 Mpc^-3 for E(B-V)=0.15, which is roughly twice that of the number density found by Erb et al. (2004) for the UV bright sample. This suggests that roughly half of the massive galaxies with M*>10^{10} Msun/h at z=2 are UV bright population, and the other half is bright in the infra-red wavelengths. The most massive galaxies at z=2 have stellar masses >= 10^{11-12} Msun. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, together with significant contribution by starbursts reaching up to 1000 Msun/yr which lie on top of the continuous component. TVD simulations indicate a more sporadic star formation history than the SPH simulations. Our results do not imply that hierarchical galaxy formation fails to account for the observed massive galaxies at z>=1. The global star formation rate density in our simulations peaks at z>=5, a much higher redshift than predicted by the semianalytic models. This star formation history suggests early build-up of the stellar mass density, and predicts that 70 (50, 30)% of the total stellar mass at z=0 had already been formed by z=1 (2, 3). Upcoming observations by Spitzer and Swift might help to better constrain the star formation history at high redshift.

preprint2004arXiv

Massive galaxies in cosmological simulations: UV-selected sample at redshift z=2

We study the properties of galaxies at z=2 in a Lambda CDM universe, using two different types of hydrodynamic simulation methods (Eulerian TVD and SPH) and a spectrophotometric analysis in the Un, G, R filter set. The simulated galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that the number density of simulated galaxies brighter than R<25.5 at z=2 is about 2e-2 h^3/Mpc^3, roughly one order of magnitude larger than that of Lyman break galaxies at z=3. The most massive galaxies at z=2 have stellar masses >~1e11 Msun, and their observed-frame G-R colors lie in the range 0.0<G-R<1.0. They typically have been continuously forming stars with a rate exceeding 30 Msun/yr over a few Gyrs from z=10 to z=2, although the TVD simulation indicates a more sporadic star formation history than the SPH simulations. Of order half of their stellar mass was already assembled by z~4. The reddest massive galaxies at z=2 with G-R >= 1.0 and Mstar>1e10 Msun/h finished the build-up of their stellar mass by z~3. Interestingly, our study suggests that the majority of the most massive galaxies at z=2 should be detectable at rest-frame UV wavelengths, contrary to some recent claims made on the basis of near-IR studies of galaxies at the same epoch, provided the median extinction is less than E(B-V)<0.3. However, our results also suggest that the fraction of stellar mass contained in galaxies that pass the color-selection criteria could be as low as 50% of the total stellar mass in the Universe at z=2. Our simulations suggest that the missing stellar mass is contained in fainter (R>25.5) and intrinsically redder galaxies. Our results do not suggest that hierarchical galaxy formation fails to account for the massive galaxies at z>=1. (abridged)

preprint2003arXiv

Early Structure Formation and Reionization in a Warm Dark Matter Cosmology

We study first structure formation in Lambda-dominated universes using large cosmological N-body/SPH simulations. We consider a standard LCDM model and a LWDM model in which the mass of the dark matter particles is taken to be m_X=10 keV. The linear power spectrum for the LWDM model has a characteristic cut-off at a wavenumber k=200 /Mpc, suppressing the formation of low mass (< 10^6 Msun) nonlinear objects early on. The absence of low mass halos in the WDM model makes the formation of primordial gas clouds with molecular hydrogen very inefficient at high redshifts. The first star-forming gas clouds form at z~21 in the WDM model, considerably later than in the CDM counterpart, and the abundance of these gas clouds differs by an order of magnitude between the two models. We carry out radiative transfer calculations by embedding massive Population III stars in the gas clouds. We show that the volume fraction of ionized gas rises up close to 100% by z=18 in the CDM case, whereas that of the WDM model remains extremely small at a level of a few percent. Thus the WDM model with m_X=10 keV is strongly inconsistent with the observed high optical depth by the WMAP satellite.

preprint2002arXiv

Cosmological SPH simulations: A hybrid multi-phase model for star formation

We present a model for star formation and supernova feedback that describes the multi-phase structure of star forming gas on scales that are typically not resolved in cosmological simulations. Our approach includes radiative heating and cooling, the growth of cold clouds embedded in an ambient hot medium, star formation in these clouds, feedback from supernovae in the form of thermal heating and cloud evaporation, galactic winds and outflows, and metal enrichment. Implemented using SPH, our scheme is a significantly modified and extended version of the grid-based method of Yepes et al. (1997), and enables us to achieve high dynamic range in simulations of structure formation. We discuss properties of the feedback model in detail and show that it predicts a self-regulated, quiescent mode of star formation, which, in particular, stabilises the star forming gaseous layers of disk galaxies. The parameterisation of this mode can be reduced to a single free quantity which determines the overall timescale for star formation. We fix this parameter to match the observed rates of star formation in local disk galaxies. When normalised in this manner, cosmological simulations nevertheless overproduce the observed cosmic abundance of stellar material. We are thus motivated to extend our feedback model to include galactic winds associated with star formation. Using small-scale simulations of individual star-forming disk galaxies, we show that these winds produce either galactic fountains or outflows, depending on the depth of the gravitational potential. Moreover, outflows from galaxies in these simulations drive chemical enrichment of the intergalactic medium, in principle accounting for the presence of metals in the Lyman alpha forest. (abridged)

preprint2002arXiv

High redshift galaxies and the Lyman-alpha forest in a CDM universe

We use a cosmological hydrodynamic simulation of a cold dark matter universe to investigate theoretically the relationship between high redshift galaxies and the Lyman=alpha forest at redshift z=3. Galaxies in the simulation are surrounded by halos of hot gas, which nevertheless contain enough neutral hydrogen to cause a Ly-alpha flux decrement, its strength increasing with galaxy mass. A comparison with recent observational data by Adelberger et. al on the Ly-alpha forest around galaxies reveals that actual galaxies may have systematically less Ly-alpha absorption within 1 Mpc of them than our simulated galaxies. In order to investigate this possibility, we add several simple prescriptions for galaxy feedback on the IGM to the evolved simulation. These include the effect of photoionizing background radiation coming from galactic sources, galactic winds whose only effect is to deposit thermal energy into the IGM, and another, kinetic model for winds, which are assumed to evacuate cavities in the IGM around galaxies. We find that only the latter is able to produce a large effect, enough to match the tentative observational data, given the energy available from star formation in the simulated galaxies. Another intriguing possibility is that a selection effect is responsible, so that galaxies with low Ly-alpha absorption are preferentially included in the sample. This is also viable, but predicts very different galaxy properties (including clustering) than the other scenarios.

preprint2000arXiv

Weakly Self-Interacting Dark Matter and the Structure of Dark Halos

We study the formation of dark halos in a $Λ$CDM universe under the assumption that Cold Dark Matter particles have a finite cross-section for elastic collisions. We compare evolution when CDM mean free paths are comparable to halo sizes with the collisionless and fluid limits. We show that a few collisions per particle per Hubble time at halo centre can substantially affect the central density profile. Cross-sections an order of magnitude larger produce sufficient relaxation for rich clusters to develop core radii in the range 100-200 $h^{-1}$kpc. The structural evolution of halos is a competition between collisional relaxation caused by individual particle interactions and violent relaxation resulting from the infall and merging processes by which clusters grow. Although our simulations concentrate on systems of cluster size, we can scale our results to address the halo structure expected for dwarf galaxies. We find that collision cross-sections sufficiently large to significantly modify the cores of such galaxies produce cluster cores which are too large and/or too round to be consistent with observation. Thus the simplest model for self-interacting dark matter is unable to improve fits to published dwarf galaxy rotation curves without violating other observational constraints.

preprint1998arXiv

Tidal tails in CDM cosmologies

We study the formation of tidal tails in pairs of merging disk galaxies with structural properties motivated by current theories of cold dark matter (CDM) cosmologies. In a recent study, Dubinski, Mihos & Hernquist (1996) showed that the formation of prominent tidal tails can be strongly suppressed by massive and extended dark haloes. For the large halo-to-disk mass ratio expected in CDM cosmologies their sequence of models failed to produce strong tails like those observed in many well-known pairs of interacting galaxies. In order to test whether this effect can constrain the viability of CDM cosmologies, we construct N-body models of disk galaxies with structural properties derived in analogy to the analytical work of Mo, Mao & White (1998). With a series of self-consistent collisionless simulations of galaxy-galaxy mergers we demonstrate that even the disks of very massive dark haloes have no problems developing long tidal tails, provided the halo spin parameter is large enough. We show that the halo-to-disk mass ratio is a poor indicator for the ability to produce tails. Instead, the relative size of disk and halo, or alternatively, the ratio of circular velocity to local escape speed at the half mass radius of the disk are more useful criteria. This result holds in all CDM cosmologies. The length of tidal tails is thus unlikely to provide useful constraints on such models.