Showing papers by "Oliver Hahn published in 2013"

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

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

The AGORA high-resolution galaxy simulations comparison project

Abstract: We introduce the Assembling Galaxies Of Resolved Anatomy (AGORA) project, a comprehensive numerical study of well-resolved galaxies within the ΛCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of eight galaxies with halo masses M_(vir) ≃ 10^(10), 10^(11), 10^(12), and 10^(13) M_☉ at z = 0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit and validated against observations to verify that the solutions are robust—i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The initial conditions for the AGORA galaxies as well as simulation outputs at various epochs will be made publicly available to the community. The proof-of-concept dark-matter-only test of the formation of a galactic halo with a z = 0 mass of M_(vir) ≃ 1.7 × 10^(11) M_☉ by nine different versions of the participating codes is also presented to validate the infrastructure of the project.

The AGORA High-Resolution Galaxy Simulations Comparison Project

Abstract: We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

The warm dark matter halo mass function below the cut-off scale

Abstract: Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts.Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of lowmass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

Halo-to-halo Similarity and Scatter in the Velocity Distribution of Dark Matter

Abstract: We examine the velocity distribution function (VDF) in dark matter halos from Milky Way to cluster mass scales. We identify an empirical model for the VDF with a wider peak and a steeper tail than a Maxwell-Boltzmann distribution, and discuss physical explanations. We quantify sources of scatter in the VDF of cosmological halos and their implication for direct detection of dark matter. Given modern simulations and observations, we find that the most significant uncertainty in the VDF of the Milky Way arises from the unknown radial position of the solar system relative to the dark matter halo scale radius.

A new approach to simulating collisionless dark matter fluids

Abstract: Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.

Rhapsody. I. Structural Properties and Formation History from a Statistical Sample of Re-simulated Cluster-size Halos

Abstract: We present the rst results from the Rhapsody cluster re-simulation project: a sample of 96 \zoom-in" simulations of dark matter halos of 10 14:8 0:05 h 1 M , selected from a 1 h 3 Gpc 3 volume. This simulation suite is the rst to resolve this many halos with 5 10 6 particles per halo in the cluster mass regime, allowing us to statistically characterize the distribution of and correlation between halo properties at xed mass. We focus on the properties of the main halos and how they are aected by formation history, which we track back to z = 12, over ve decades in mass. We give particular attention to the impact of the formation history on the density proles of the halos. We nd that the deviations from the Navarro{Frenk{White (NFW) model and the Einasto model depend on formation time. Late-forming halos tend to have considerable deviations from both models, partly due to the presence of massive subhalos, while early-forming halos deviate less but still signicantly from the NFW model and are better described by the Einasto model. We nd that the halo shapes depend only moderately on formation time. Departure from spherical symmetry impacts the density proles through the anisotropic distribution of massive subhalos. Further evidence of the impact of subhalos is provided by analyzing the phase-space structure. A detailed analysis of the properties of the subhalo population in Rhapsody is presented in a companion paper.

Virial scaling of galaxies in clusters: bright to faint is cool to hot

Abstract: By combining galaxy tracers from high-resolution N-body and hydrodynamical simulations, we present a consistent picture of the behaviour of galaxy velocities in massive clusters. In haloes above ~ 10^14 Msun, the brightest satellite galaxies are slightly cooler compared to the dark matter, while fainter satellites are hotter. Within the virial radius of a cluster, the mean velocity dispersion based on the 100 brightest galaxies is a factor of 1.065 +/- 0.005 (stat) +/- 0.027 (sys) higher than that of the dark matter (corresponding to a ~10-15 per cent bias in the dynamical mass estimate) while that based on only the five brightest galaxies is 0.868 +/- 0.039 (stat) +/- 0.035 (sys). These trends are approximately independent of redshift. The velocity structure is sensitive to the modelling of galaxies in clusters, indicative of the complex interplay of tidal stripping, dynamical friction, and merging. Velocity dispersions derived from instantaneous subhalo properties are larger than those employing either peak subhalo properties or hydrodynamical galaxy tracers. The latter two methods are consistent, implying that stacked spectroscopic analysis of cluster samples should, after correction for projection, show a trend towards slightly higher velocities when fainter galaxies are included, with an unbiased measure of dark matter velocity dispersion coming from approximately 30 galaxies per cluster. We show evidence that the velocity distribution function of bright galaxies near the cluster centre has a low-velocity tail due to strong dynamical friction.

How closely do baryons follow dark matter on large scales

Abstract: We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe.

Rhapsody. ii. subhalo properties and the impact of tidal stripping from a statistical sample of cluster-size halos

Abstract: We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the RHAPSODY simulations introduced in Wu et al. We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.

Population iii star formation in large cosmological volumes. i. halo temporal and physical environment

Abstract: We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population?III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population?III star forming halos to be constructed, which will lead to Population?III star formation simulations that more accurately reflect the diversity of Population?III stars, both in time and halo mass. This model shows that Population?III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z ~ 10, and likely beyond, though Population?III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z = 10. A catalog of more than 40,000 candidate Population?III forming halos were identified, with formation times temporally ranging from z = 30 to z = 10, and ranging in mass from 2.3 ? 105 M ? to 1.2 ? 1010 M ?. At early times, the environment that Population?III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population?III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population?III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population?III star formation by z ~ 10.

Population III Star Formation In Large Cosmological Simulations I. Halo Temporal and Physical Environment

Abstract: We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population III star forming halos to be constructed, which will lead to Population III star formation simulations that more accurately reflect the diversity of Population III stars, both in time and halo mass. This model shows that Population III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z~10, and likely beyond, though Population III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z=10. A catalog of more than 40,000 candidate Population III forming halos were identified, with formation times temporally ranging from z=30 to z=10, and ranging in mass from 2.3x10^5 M_sun to 1.2x10^10 M_sun. At early times, the environment that Population III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population III star formation by z~10.

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

Abstract: We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.