-
Constraining Reionization with the $z \sim 5-6$ Lyman-$α$ Forest Power Spectrum: the Outlook after Planck
Authors:
J. Oñorbe,
J. F. Hennawi,
Z. Lukić,
M. Walther
Abstract:
The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparabl…
▽ More
The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparable to the Hubble time, memory of reionization heating is retained. Therefore, a late ending reionization has significant implications for the structure of the $z\sim5-6$ Lyman-$α$ (ly$α$) forest. Using state-of-the-art hydrodynamical simulations that allow us to vary the timing of reionization and its associated heat injection, we argue that extant thermal signatures from reionization can be detected via the Ly$α$ forest power spectrum at $5< z<6$. This arises because the small-scale cutoff in the power depends not only the the IGM temperature at these epochs, but is also particularly sensitive to the pressure smoothing scale set by the IGM full thermal history. Comparing our different reionization models with existing measurements of the Ly$α$ forest flux power spectrum at $z=5.0-5.4$, we find that models satisfying Planck's $τ_e$, constraint favor a moderate amount of heat injection consistent with galaxies driving reionization, but disfavoring quasar-driven scenarios. We explore the impact of different reionization histories and heating models on the shape of the power spectrum, and find that they can produce similar effects, but argue that this degeneracy can be broken with high enough quality data. We study the feasibility of measuring the flux power spectrum at $z\simeq 6$ using mock quasar spectra and conclude that a sample of $\sim10$ high-resolution spectra with an attainable signal-to-noise ratio will allow distinguishing between different reionization scenarios.
△ Less
Submitted 27 April, 2019; v1 submitted 24 March, 2017;
originally announced March 2017.
-
Implications of $z \sim 6$ Quasar Proximity Zones for the Epoch of Reionization and Quasar Lifetimes
Authors:
Anna-Christina Eilers,
Frederick B. Davies,
Joseph F. Hennawi,
J. Xavier Prochaska,
Zarija Lukić,
Chiara Mazzucchelli
Abstract:
We study quasar proximity zones in the redshift range $5.77 \leq z \leq 6.54$ by homogeneously analyzing $34$ medium resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift, and argued that this provides evidence for a r…
▽ More
We study quasar proximity zones in the redshift range $5.77 \leq z \leq 6.54$ by homogeneously analyzing $34$ medium resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift, and argued that this provides evidence for a rapidly evolving intergalactic medium (IGM) neutral fraction during reionization, we measure a much shallower trend $\propto(1+z)^{-1.44}$. We compare our measured proximity zone sizes to predictions from hydrodynamical simulations post-processed with one-dimensional radiative transfer, and find good agreement between observations and theory irrespective of the ionization state of the ambient IGM. This insensitivity to IGM ionization state has been previously noted, and results from the fact that the definition of proximity zone size as the first drop of the smoothed quasar spectrum below the $10\%$ flux transmission level probes locations where the ionizing radiation from the quasar is an order of magnitude larger than the expected ultraviolet ionizing background that sets the neutral fraction of the IGM. Our analysis also uncovered three objects with exceptionally small proximity zones (two have $R_p < 1$proper Mpc), which constitute outliers from the observed distribution and are challenging to explain with our radiative transfer simulations. We consider various explanations for their origin, such as strong absorption line systems associated with the quasar or patchy reionization, but find that the most compelling scenario is that these quasars have been shining for $\lesssim 10^5$yr.
△ Less
Submitted 10 May, 2017; v1 submitted 7 March, 2017;
originally announced March 2017.
-
Measuring alignments between galaxies and the cosmic web at $z \sim 2-3$ using IGM tomography
Authors:
Alex Krolewski,
Khee-Gan Lee,
Zarija Lukic,
Martin White
Abstract:
Many galaxy formation models predict alignments between galaxy spin and the cosmic web (i.e. the directions of filaments and sheets), leading to intrinsic alignment between galaxies that creates a systematic error in weak lensing measurements. These effects are often predicted to be stronger at high-redshifts ($z\gtrsim1$) that are inaccessible to massive galaxy surveys on foreseeable instrumentat…
▽ More
Many galaxy formation models predict alignments between galaxy spin and the cosmic web (i.e. the directions of filaments and sheets), leading to intrinsic alignment between galaxies that creates a systematic error in weak lensing measurements. These effects are often predicted to be stronger at high-redshifts ($z\gtrsim1$) that are inaccessible to massive galaxy surveys on foreseeable instrumentation, but IGM tomography of the Ly$α$ forest from closely-spaced quasars and galaxies is starting to measure the $z\sim2-3$ cosmic web with the requisite fidelity. Using mock surveys from hydrodynamical simulations, we examine the utility of this technique, in conjunction with coeval galaxy samples, to measure alignment between galaxies and the cosmic web at $z\sim2.5$. We show that IGM tomography surveys with $\lesssim5$ $h^{-1}$ Mpc sightline spacing can accurately recover the eigenvectors of the tidal tensor, which we use to define the directions of the cosmic web. For galaxy spins and shapes, we use a model parametrized by the alignment strength, $Δ\langle\cosθ\rangle$, with respect to the tidal tensor eigenvectors from the underlying density field, and also consider observational effects such as errors in the galaxy position angle, inclination, and redshift. Measurements using the upcoming $\sim1\,\mathrm{deg}^2$ CLAMATO tomographic survey and 600 coeval zCOSMOS-Deep galaxies should place $3σ$ limits on extreme alignment models with $Δ\langle\cosθ\rangle\sim0.1$, but much larger surveys encompassing $>10,000$ galaxies, such as Subaru PFS, will be required to constrain models with $Δ\langle\cosθ\rangle\sim0.03$. These measurements will constrain models of galaxy-cosmic web alignment and test tidal torque theory at $z\sim2$, improving our understanding of the redshift dependence of galaxy-cosmic web alignment and the physics of intrinsic alignments.
△ Less
Submitted 30 November, 2016;
originally announced December 2016.
-
The DESI Experiment Part II: Instrument Design
Authors:
DESI Collaboration,
Amir Aghamousa,
Jessica Aguilar,
Steve Ahlen,
Shadab Alam,
Lori E. Allen,
Carlos Allende Prieto,
James Annis,
Stephen Bailey,
Christophe Balland,
Otger Ballester,
Charles Baltay,
Lucas Beaufore,
Chris Bebek,
Timothy C. Beers,
Eric F. Bell,
José Luis Bernal,
Robert Besuner,
Florian Beutler,
Chris Blake,
Hannes Bleuler,
Michael Blomqvist,
Robert Blum,
Adam S. Bolton,
Cesar Briceno
, et al. (268 additional authors not shown)
Abstract:
DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from…
▽ More
DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.
△ Less
Submitted 13 December, 2016; v1 submitted 31 October, 2016;
originally announced November 2016.
-
The DESI Experiment Part I: Science,Targeting, and Survey Design
Authors:
DESI Collaboration,
Amir Aghamousa,
Jessica Aguilar,
Steve Ahlen,
Shadab Alam,
Lori E. Allen,
Carlos Allende Prieto,
James Annis,
Stephen Bailey,
Christophe Balland,
Otger Ballester,
Charles Baltay,
Lucas Beaufore,
Chris Bebek,
Timothy C. Beers,
Eric F. Bell,
José Luis Bernal,
Robert Besuner,
Florian Beutler,
Chris Blake,
Hannes Bleuler,
Michael Blomqvist,
Robert Blum,
Adam S. Bolton,
Cesar Briceno
, et al. (268 additional authors not shown)
Abstract:
DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure…
▽ More
DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.
△ Less
Submitted 13 December, 2016; v1 submitted 31 October, 2016;
originally announced November 2016.
-
Self-Consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations
Authors:
Jose Oñorbe,
Joseph F Hennawi,
Zarija Lukić
Abstract:
The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rate…
▽ More
The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. We show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier z ~ 15 than they should. This problem arises because at z > 6, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyman-alpha forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. We make our new UVB photoionization and photoheating rates publicly available for use in future simulations.
△ Less
Submitted 16 March, 2017; v1 submitted 14 July, 2016;
originally announced July 2016.
-
ASCR/HEP Exascale Requirements Review Report
Authors:
Salman Habib,
Robert Roser,
Richard Gerber,
Katie Antypas,
Katherine Riley,
Tim Williams,
Jack Wells,
Tjerk Straatsma,
A. Almgren,
J. Amundson,
S. Bailey,
D. Bard,
K. Bloom,
B. Bockelman,
A. Borgland,
J. Borrill,
R. Boughezal,
R. Brower,
B. Cowan,
H. Finkel,
N. Frontiere,
S. Fuess,
L. Ge,
N. Gnedin,
S. Gottlieb
, et al. (29 additional authors not shown)
Abstract:
This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 ti…
▽ More
This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.
△ Less
Submitted 31 March, 2016; v1 submitted 30 March, 2016;
originally announced March 2016.
-
Modeling the Lyman-alpha Forest in Collisionless Simulations
Authors:
Daniele Sorini,
José Oñorbe,
Zarija Lukić,
Joseph F. Hennawi
Abstract:
Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we pres…
▽ More
Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present Iteratively Matched Statistics (IMS), a novel method to accurately model the Lyman-alpha forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyman-alpha forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can perfectly reproduce line-of-sight observables, such as the PDF and power spectrum, and accurately reproduce the 3D flux power spectrum (5-20%). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic "mock" skies for Lyman-alpha forest surveys.
△ Less
Submitted 25 February, 2016;
originally announced February 2016.
-
HACC: Simulating Sky Surveys on State-of-the-Art Supercomputing Architectures
Authors:
Salman Habib,
Adrian Pope,
Hal Finkel,
Nicholas Frontiere,
Katrin Heitmann,
David Daniel,
Patricia Fasel,
Vitali Morozov,
George Zagaris,
Tom Peterka,
Venkatram Vishwanath,
Zarija Lukic,
Saba Sehrish,
Wei-keng Liao
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…
▽ More
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.
△ Less
Submitted 8 October, 2014;
originally announced October 2014.
-
The Lyman-$α$ forest in optically-thin hydrodynamical simulations
Authors:
Zarija Lukić,
Casey Stark,
Peter Nugent,
Martin White,
Avery Meiksin,
Ann Almgren
Abstract:
We study the statistics of the Lyman-$α$ forest in a flat LCDM cosmology with the N-body + Eulerian hydrodynamics code Nyx. We produce a suite of simulations, covering the observationally relevant redshift range $2 \leq z \leq 4$. We find that a grid resolution of 20 kpc/h is required to produce one percent convergence of Lyman-$α$ flux statistics, up to k = 10 h/Mpc. In addition to establishing r…
▽ More
We study the statistics of the Lyman-$α$ forest in a flat LCDM cosmology with the N-body + Eulerian hydrodynamics code Nyx. We produce a suite of simulations, covering the observationally relevant redshift range $2 \leq z \leq 4$. We find that a grid resolution of 20 kpc/h is required to produce one percent convergence of Lyman-$α$ flux statistics, up to k = 10 h/Mpc. In addition to establishing resolution requirements, we study the effects of missing modes in these simulations, and find that box sizes of L > 40 Mpc/h are needed to suppress numerical errors to a sub-percent level. Our optically-thin simulations with the ionizing background prescription of Haardt & Madau (2012) reproduce an IGM equation of state with $T_0 \approx 10^4 K$ and $γ\approx 1.55$ at z=2, with a mean transmitted flux close to the observed values. When using the ionizing background prescription of Faucher-Giguere et al. (2009), the mean flux is 10-15 per cent below observed values at z=2, and a factor of 2 too small at z = 4. We show the effects of the common practice of rescaling optical depths to the observed mean flux and how it affects convergence rates. We also investigate the common practice of `splicing' results from a number of different simulations to estimate the 1D flux power spectrum and show it is accurate at the 10 percent level. Finally, we find that collisional heating of the gas from dark matter particles is negligible in modern cosmological simulations.
△ Less
Submitted 23 December, 2015; v1 submitted 24 June, 2014;
originally announced June 2014.
-
A new method for imaging nuclear threats using cosmic ray muons
Authors:
C. L. Morris,
Jeffrey Bacon,
Konstantin Borozdin,
Haruo Miyadera,
John Perry,
Evan Rose,
Scott Watson,
Timothy White,
Derek Aberle,
J. Andrew Green,
George G. McDuff,
Zarija Lukić,
Edward C. Milner
Abstract:
Muon tomography is a technique that uses cosmic ray muons to generate three dimensional images of volumes using information contained in the Coulomb scattering of the muons. Advantages of this technique are the ability of cosmic rays to penetrate significant overburden and the absence of any additional dose delivered to subjects under study above the natural cosmic ray flux. Disadvantages include…
▽ More
Muon tomography is a technique that uses cosmic ray muons to generate three dimensional images of volumes using information contained in the Coulomb scattering of the muons. Advantages of this technique are the ability of cosmic rays to penetrate significant overburden and the absence of any additional dose delivered to subjects under study above the natural cosmic ray flux. Disadvantages include the relatively long exposure times and poor position resolution and complex algorithms needed for reconstruction. Here we demonstrate a new method for obtaining improved position resolution and statistical precision for objects with spherical symmetry.
△ Less
Submitted 4 June, 2013; v1 submitted 3 June, 2013;
originally announced June 2013.
-
Structure Finding in Cosmological Simulations: The State of Affairs
Authors:
Alexander Knebe,
Frazer R. Pearce,
Hanni Lux,
Yago Ascasibar,
Peter Behroozi,
Javier Casado,
Christine Corbett Moran,
Juerg Diemand,
Klaus Dolag,
Rosa Dominguez-Tenreiro,
Pascal Elahi,
Bridget Falck,
Stefan Gottloeber,
Jiaxin Han,
Anatoly Klypin,
Zarija Lukic,
Michal Maciejewski,
Cameron K. McBride,
Manuel E. Merchan,
Stuart I. Muldrew,
Mark Neyrinck,
Julian Onions,
Susana Planelles,
Doug Potter,
Vicent Quilis
, et al. (10 additional authors not shown)
Abstract:
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 remain…
▽ More
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.
△ Less
Submitted 26 July, 2013; v1 submitted 2 April, 2013;
originally announced April 2013.
-
Nyx: A Massively Parallel AMR Code for Computational Cosmology
Authors:
Ann Almgren,
John Bell,
Mike Lijewski,
Zarija Lukić,
Ethan Van Andel
Abstract:
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 (AMR); a particle-mesh (PM) scheme using the same…
▽ More
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 (AMR); a particle-mesh (PM) 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.
△ Less
Submitted 18 January, 2013;
originally announced January 2013.
-
The Universe at Extreme Scale: Multi-Petaflop Sky Simulation on the BG/Q
Authors:
Salman Habib,
Vitali Morozov,
Hal Finkel,
Adrian Pope,
Katrin Heitmann,
Kalyan Kumaran,
Tom Peterka,
Joe Insley,
David Daniel,
Patricia Fasel,
Nicholas Frontiere,
Zarija Lukic
Abstract:
Remarkable observational advances have established a compelling cross-validated model of the Universe. Yet, two key pillars of this model -- dark matter and dark energy -- remain mysterious. Sky surveys that map billions of galaxies to explore the `Dark Universe', demand a corresponding extreme-scale simulation capability; the HACC (Hybrid/Hardware Accelerated Cosmology Code) framework has been de…
▽ More
Remarkable observational advances have established a compelling cross-validated model of the Universe. Yet, two key pillars of this model -- dark matter and dark energy -- remain mysterious. Sky surveys that map billions of galaxies to explore the `Dark Universe', demand a corresponding extreme-scale simulation capability; the HACC (Hybrid/Hardware Accelerated Cosmology Code) framework has been designed to deliver this level of performance now, and into the future. With its novel algorithmic structure, HACC allows flexible tuning across diverse architectures, including accelerated and multi-core systems.
On the IBM BG/Q, HACC attains unprecedented scalable performance -- currently 13.94 PFlops at 69.2% of peak and 90% parallel efficiency on 1,572,864 cores with an equal number of MPI ranks, and a concurrency of 6.3 million. This level of performance was achieved at extreme problem sizes, including a benchmark run with more than 3.6 trillion particles, significantly larger than any cosmological simulation yet performed.
△ Less
Submitted 19 November, 2012;
originally announced November 2012.
-
Obtaining material identification with cosmic ray radiography
Authors:
Christopher Morris,
Konstantin Borozdin,
Jeffrey Bacon,
Elliott Chen,
Zarija Lukić,
Edward Milner,
Haruo Miyadera,
John Perry,
Dave Schwellenbach,
Derek Aberle,
Wendi Dreesen,
J. Andrew Green,
George G. McDuff,
Kanetada Nagamine,
Michael Sossong,
Candace Spore,
Nathan Toleman
Abstract:
The passage of muons through matter is mostly affected by their Coulomb interactions with electrons and nuclei. The muon interactions with electrons lead to continuous energy loss and stopping of muons, while their scattering off nuclei lead to angular 'diffusion'. By measuring both the number of stopped muons and angular changes in muon trajectories we can estimate density and identify materials.…
▽ More
The passage of muons through matter is mostly affected by their Coulomb interactions with electrons and nuclei. The muon interactions with electrons lead to continuous energy loss and stopping of muons, while their scattering off nuclei lead to angular 'diffusion'. By measuring both the number of stopped muons and angular changes in muon trajectories we can estimate density and identify materials. Here we demonstrate the material identification using data taken at Los Alamos with the Mini Muon Tracker.
△ Less
Submitted 22 October, 2012;
originally announced October 2012.
-
Cosmic Ray Radiography of the Damaged Cores of the Fukushima Reactors
Authors:
Konstantin Borozdin,
Steven Greene,
Zarija Lukić,
Edward Cas Milner,
Haruo Miyadera,
Christopher Morris,
John Perry
Abstract:
The passage of muons through matter is dominated by the Coulomb interaction with electrons and nuclei. The interaction with the electrons leads to continuous energy loss and stopping of the muons. The interaction with nuclei leads to angle diffusion. Two muon imaging methods that use flux attenuation and multiple Coulomb scattering of cosmic-ray muons are being studied as tools for diagnosing the…
▽ More
The passage of muons through matter is dominated by the Coulomb interaction with electrons and nuclei. The interaction with the electrons leads to continuous energy loss and stopping of the muons. The interaction with nuclei leads to angle diffusion. Two muon imaging methods that use flux attenuation and multiple Coulomb scattering of cosmic-ray muons are being studied as tools for diagnosing the damaged cores of the Fukushima reactors. Here we compare these two methods. We conclude that the scattering method can provide detailed information about the core. Attenuation has low contrast and little sensitivity to the core.
△ Less
Submitted 12 September, 2012;
originally announced September 2012.
-
Haloes gone MAD: The Halo-Finder Comparison Project
Authors:
Alexander Knebe,
Steffen R. Knollmann,
Stuart I. Muldrew,
Frazer R. Pearce,
Miguel Angel Aragon-Calvo,
Yago Ascasibar,
Peter S. Behroozi,
Daniel Ceverino,
Stephane Colombi,
Juerg Diemand,
Klaus Dolag,
Bridget L. Falck,
Patricia Fasel,
Jeff Gardner,
Stefan Gottloeber,
Chung-Hsing Hsu,
Francesca Iannuzzi,
Anatoly Klypin,
Zarija Lukic,
Michal Maciejewski,
Cameron McBride,
Mark C. Neyrinck,
Susana Planelles,
Doug Potter,
Vicent Quilis
, et al. (12 additional authors not shown)
Abstract:
[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 develope…
▽ More
[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).
△ Less
Submitted 5 April, 2011;
originally announced April 2011.
-
Mass Function Predictions Beyond LCDM
Authors:
Suman Bhattacharya,
Katrin Heitmann,
Martin White,
Zarija Lukić,
Christian Wagner,
Salman Habib
Abstract:
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. We investigate and discuss some of the reasons for these differences. Quantifying error sources and com…
▽ More
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. 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 $w$CDM 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\cdot 10^{11}-3\cdot 10^{15}$ M$_{\odot}$ to an estimated accuracy of about 2%. The set of $w$CDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a $Λ$CDM cosmology and others with $w\simeq-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.
△ Less
Submitted 18 March, 2011; v1 submitted 12 May, 2010;
originally announced May 2010.
-
Galaxy Clusters as a probe of early dark energy
Authors:
Ujjaini Alam,
Zarija Lukić,
Suman Bhattacharya
Abstract:
We study a class of early dark energy (EDE) models, in which, unlike in standard dark energy models, a substantial amount of dark energy exists in the matter-dominated era. We self-consistently include dark energy perturbations, and show that these models may be successfully constrained using future observations of galaxy clusters, in particular the redshift abundance, and the Sunyaev-Zel'dovich (…
▽ More
We study a class of early dark energy (EDE) models, in which, unlike in standard dark energy models, a substantial amount of dark energy exists in the matter-dominated era. We self-consistently include dark energy perturbations, and show that these models may be successfully constrained using future observations of galaxy clusters, in particular the redshift abundance, and the Sunyaev-Zel'dovich (SZ) power spectrum. We make predictions for EDE models, as well as LCDM for incoming X-ray (eROSITA) and microwave (South Pole Telescope) observations. We show that galaxy clusters' mass function and the SZ power spectrum will put strong constraints both on the equation of state of \de today and the redshift at which EDE transits to present-day LCDM like behavior for these models, thus providing complementary information to the geometric probes of dark energy. Not including perturbations in EDE models leads to those models being practically indistinguishable from LCDM. An MCMC analysis of future galaxy cluster surveys provides constraints for EDE parameters that are competitive with and complementary to background expansion observations such as supernovae.
△ Less
Submitted 15 December, 2010; v1 submitted 3 April, 2010;
originally announced April 2010.
-
Particle mesh simulations of the Lyman-alpha forest and the signature of Baryon Acoustic Oscillations in the intergalactic medium
Authors:
Martin White,
Adrian Pope,
Jordan Carlson,
Katrin Heitmann,
Salman Habib,
Patricia Fasel,
David Daniel,
Zarija Lukic
Abstract:
We present a set of ultra-large particle-mesh simulations of the LyA forest targeted at understanding the imprint of baryon acoustic oscillations (BAO) in the inter-galactic medium. We use 9 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…
▽ More
We present a set of ultra-large particle-mesh simulations of the LyA forest targeted at understanding the imprint of baryon acoustic oscillations (BAO) in the inter-galactic medium. We use 9 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 (PDFs) 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 misestimation, and non-gravitational sources such as fluctuations in the photo-ionizing background and temperature fluctuations due to HeII reionization.
△ Less
Submitted 26 February, 2010; v1 submitted 27 November, 2009;
originally announced November 2009.
-
The Structure of Halos: Implications for Group and Cluster Cosmology
Authors:
Zarija Lukić,
Darren Reed,
Salman Habib,
Katrin Heitmann
Abstract:
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 t…
▽ More
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, Omega_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 10^{12.5}-10^{15.5} M_sun/h 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. (abridged)
△ Less
Submitted 25 March, 2008;
originally announced March 2008.
-
The Cosmic Code Comparison Project
Authors:
Katrin Heitmann,
Zarija Lukic,
Patricia Fasel,
Salman Habib,
Michael S. Warren,
Martin White,
James Ahrens,
Lee Ankeny,
Ryan Armstrong,
Brian O'Shea,
Paul M. Ricker,
Volker Springel,
Joachim Stadel,
Hy Trac
Abstract:
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 make the simulation task very demanding: the simulations have to be at least as accurate as the…
▽ More
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 make 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 validation program. The first step of such a program is the comparison of different cosmology codes including gravitation 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.
△ Less
Submitted 8 June, 2007;
originally announced June 2007.
-
The Halo Mass Function: High-Redshift Evolution and Universality
Authors:
Zarija Lukic,
Katrin Heitmann,
Salman Habib,
Sergei Bashinsky,
Paul M. Ricker
Abstract:
We study the formation of dark matter halos in the concordance LCDM 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 Mpc/h, we determine the mass function and its…
▽ More
We study the formation of dark matter halos in the concordance LCDM 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 Mpc/h, 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 (10^7 - 10^13.5 M_sun/h). 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 redshift, 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 and collaborators exists at low redshifts.
△ Less
Submitted 13 January, 2008; v1 submitted 14 February, 2007;
originally announced February 2007.
-
Capturing Halos at High Redshifts
Authors:
Katrin Heitmann,
Zarija Lukic,
Salman Habib,
Paul M. Ricker
Abstract:
We study the evolution of the mass function of dark matter halos in the concordance LCDM model at high redshift. We employ overlapping (multiple-realization) numerical simulations to cover a wide range of halo masses, 10^7-10^{15}h^{-1}M_\odot, with redshift coverage beginning at z=20. The Press-Schechter mass function is significantly discrepant from the simulation results at high redshifts. Of…
▽ More
We study the evolution of the mass function of dark matter halos in the concordance LCDM model at high redshift. We employ overlapping (multiple-realization) numerical simulations to cover a wide range of halo masses, 10^7-10^{15}h^{-1}M_\odot, 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 Warren et al. (2005). 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.
△ Less
Submitted 11 January, 2006;
originally announced January 2006.