Zarija Lukić

Haloes gone MAD: The Halo-Finder Comparison Project

We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends, spherical-overdensity and phase-space-based algorithms. We

The Halo Mass Function: High-Redshift Evolution and Universality

We study the formation of dark matter halos in the concordance ΛCDM model over a wide range of redshifts, from -->z = 20 to the present. Our primary focus is the halo mass function, a key probe of cosmology. By performing a large suite of nested-box N-body simulations with careful convergence and error controls (60 simulations with box sizes from 4 to 256 h−1 Mpc), we determine the mass function and its evolution with excellent statistical and systematic errors, reaching a few percent over most of the considered redshift and mass range. Across the studied redshifts, the halo mass is probed over 6 orders of magnitude (107-1013.5 h−1 M☉). Historically, there has been considerable variation in the high-redshift mass function as obtained by different groups. We have made a concerted effort to identify and correct possible systematic errors in computing the mass function at high redshift and to explain the discrepancies between some of the previous results. We discuss convergence criteria for the required force resolution, simulation box size, halo mass range, initial and final redshifts, and time stepping. Because of conservative cuts on the mass range probed by individual boxes, our results are relatively insensitive to simulation volume, the remaining sensitivity being consistent with extended Press-Schechter theory. Previously obtained mass function fits near -->z = 0, when scaled by linear theory, are in good agreement with our results at all redshifts, although a mild redshift dependence consistent with that found by Reed et al. may exist at low redshifts. Overall, our results are consistent with a universal form for the mass function at high redshifts.

MASS FUNCTION PREDICTIONS BEYOND ΛCDM

The statistics of dark matter halos is an essential component of precision cosmology. The mass distribution of halos, as specified by the halo mass function, is a key input for several cosmological probes. The sizes of N-body simulations are now such that, for the most part, results need no longer be statistics-limited, but are still subject to various systematic uncertainties. Discrepancies in the results of simulation campaigns for the halo mass function remain in excess of statistical uncertainties and of roughly the same size as the error limits set by near-future observations; we investigate and discuss some of the reasons for these differences. Quantifying error sources and compensating for them as appropriate, we carry out a high-statistics study of dark matter halos from 67 N-body simulations to investigate the mass function and its evolution for a reference ?CDM cosmology and for a set of wCDM cosmologies. For the reference ?CDM cosmology (close to WMAP5), we quantify the breaking of universality in the form of the mass function as a function of redshift, finding an evolution of as much as 10% away from the universal form between redshifts z = 0 and z = 2. For cosmologies very close to this reference we provide a fitting formula to our results for the (evolving) ?CDM mass function over a mass range of 6 ? 1011-3 ? 1015 M ? to an estimated accuracy of about 2%. The set of wCDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a ?CDM cosmology and others with w ?1) are described by the fitting formula for the reference ?CDM case at an accuracy level of 10%, but with clear systematic deviations. We argue that, as a consequence, fitting formulae based on a universal form for the mass function may have limited utility in high-precision cosmological applications.

The cosmic code comparison project

Current and upcoming cosmological observations allow us to probe structures on smaller and smaller scales, entering highly nonlinear regimes. In order to obtain theoretical predictions in these regimes, large cosmological simulations have to be carried out. The promised high accuracy from observations makes the simulation task very demanding: the simulations have to be at least as accurate as the observations. This requirement can only be fulfilled by carrying out an extensive code verification program. The first step of such a program is the comparison of different cosmology codes including gravitational interactions only. In this paper, we extend a recently carried out code comparison project to include five more simulation codes. We restrict our analysis to a small cosmological volume which allows us to investigate properties of halos. For the matter power spectrum and the mass function, the previous results hold, with the codes agreeing at the 10% level over wide dynamic ranges. We extend our analysis to the comparison of halo profiles and investigate the halo count as a function of local density. We introduce and discuss ParaView as a flexible analysis tool for cosmological simulations, the use of which immensely simplifies the code comparison task.

The Structure of Halos: Implications for Group and Cluster Cosmology

The dark matter halo mass function is a key repository of cosmological information over a wide range of mass scales, from individual galaxies to galaxy clusters. N-body simulations have established that the friends-of-friends (FOF) mass function has a universal form to a surprising level of accuracy (10%). The high-mass tail of the mass function is exponentially sensitive to the amplitude of the initial density perturbations, the mean matter density parameter, Ω m , and to the dark energy controlled late-time evolution of the density field. Observed group and cluster masses, however, are usually stated in terms of a spherical overdensity (SO) mass which does not map simply to the FOF mass. Additionally, the widely used halo models of structure formation—and halo occupancy distribution descriptions of galaxies within halos—are often constructed exploiting the universal form of the FOF mass function. This again raises the question of whether FOF halos can be simply related to the notion of a spherical overdensity mass. By employing results from Monte Carlo realizations of ideal Navarro-Frenk-White (NFW) halos and N-body simulations, we study the relationship between the two definitions of halo mass. We find that the vast majority of halos (80%-85%) in the mass-range 1012.5-1015.5 h –1 M ☉ indeed allow for an accurate mapping between the two definitions (~5%), but only if the halo concentrations are known. Nonisolated halos fall into two broad classes: those with complex substructure that are poor fits to NFW profiles and those bridged by the (isodensity-based) FOF algorithm. A closer investigation of the bridged halos reveals that the fraction of these halos and their satellite mass distribution is cosmology dependent. We provide a preliminary discussion of the theoretical and observational ramifications of these results.

Nyx: A Massively Parallel AMR Code for Computational Cosmology

We present a new N-body and gas dynamics code, called Nyx, for large-scale cosmological simulations. Nyx follows the temporal evolution of a system of discrete dark matter particles gravitationally coupled to an inviscid ideal fluid in an expanding universe. The gas is advanced in an Eulerian framework with block-structured adaptive mesh refinement; a particle-mesh scheme using the same grid hierarchy is used to solve for self-gravity and advance the particles. Computational results demonstrating the validation of Nyx on standard cosmological test problems, and the scaling behavior of Nyx to 50,000 cores, are presented.

PARTICLE MESH SIMULATIONS OF THE Lyα FOREST AND THE SIGNATURE OF BARYON ACOUSTIC OSCILLATIONS IN THE INTERGALACTIC MEDIUM

We present a set of ultra-large particle-mesh simulations of the Lyα forest targeted at understanding the imprint of baryon acoustic oscillations in the inter-galactic medium. We use nine dark matter only simulations which can, for the first time, simultaneously resolve the Jeans scale of the intergalactic gas while covering the large volumes required to adequately sample the acoustic feature. Mock absorption spectra are generated using the fluctuating Gunn-Peterson approximation which have approximately correct flux probability density functions and small-scale power spectra. On larger scales, there is clear evidence in the redshift-space correlation function for an acoustic feature, which matches a linear theory template with constant bias. These spectra, which we make publicly available, can be used to test pipelines, plan future experiments, and model various physical effects. As an illustration, we discuss the basic properties of the acoustic signal in the forest, the scaling of errors with noise and source number density, modified statistics to treat mean flux evolution and mis-estimation, and non-gravitational sources such as fluctuations in the photoionizing background and temperature fluctuations due to He II reionization.

HACC: Simulating sky surveys on state-of-the-art supercomputing architectures

Abstract Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind the two major components of the ‘Dark Universe’, dark energy and dark matter. In addition, the surveys also probe primordial perturbations and carry out fundamental measurements, such as determining the sum of neutrino masses. Large-scale simulations of structure formation in the Universe play a critical role in the interpretation of the data and extraction of the physics of interest. Just as survey instruments continue to grow in size and complexity, so do the supercomputers that enable these simulations. Here we report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently developed and evolving cosmology N-body code framework, designed to run efficiently on diverse computing architectures and to scale to millions of cores and beyond. HACC can run on all current supercomputer architectures and supports a variety of programming models and algorithms. It has been demonstrated at scale on Cell- and GPU-accelerated systems, standard multi-core node clusters, and Blue Gene systems. HACC’s design allows for ease of portability, and at the same time, high levels of sustained performance on the fastest supercomputers available. We present a description of the design philosophy of HACC, the underlying algorithms and code structure, and outline implementation details for several specific architectures. We show selected accuracy and performance results from some of the largest high resolution cosmological simulations so far performed, including benchmarks evolving more than 3.6 trillion particles.

Capturing halos at high redshifts

We study the evolution of the mass function of dark matter halos in the concordance ΛCDM model at high redshift. We employ overlapping (multiple realization) numerical simulations to cover a wide range of halo masses, 107-1015 h-1 M☉, with redshift coverage beginning at z = 20. The Press-Schechter mass function is significantly discrepant from the simulation results at high redshifts. Of the more recently proposed mass functions, our results are in best agreement with those of Warren and coworkers. The statistics of the simulations—along with good control over systematics—allow for fits accurate to the level of 20% at all redshifts. We provide a concise discussion of various issues in defining and computing the halo mass function and how these are addressed in our simulations.

Imaging Fukushima Daiichi reactors with muons

A study of imaging the Fukushima Daiichi reactors with cosmic-ray muons to assess the damage to the reactors is presented. Muon scattering imaging has high sensitivity for detecting uranium fuel and debris even through thick concrete walls and a reactor pressure vessel. Technical demonstrations using a reactor mockup, detector radiation test at Fukushima Daiichi, and simulation studies have been carried out. These studies establish feasibility for the reactor imaging. A few months of measurement will reveal the spatial distribution of the reactor fuel. The muon scattering technique would be the best and probably the only way for Fukushima Daiichi to make this determination in the near future.