28 research outputs found
Degeneracies between Modified Gravity and Baryonic Physics
In order to determine the observable signatures of modified gravity theories,
it is important to consider the effect of baryonic physics. We use a modified
version of the ISIS code to run cosmological hydrodynamic simulations to study
degeneracies between modified gravity and radiative hydrodynamical processes.
Of these, one was the standard CDM model and four were variations of
the Symmetron model. For each model we ran three variations of baryonic
processes: non-radiative hydrodynamics; cooling and star formation; and
cooling, star formation, and supernova feedback. We construct stacked gas
density, temperature, and dark matter density profiles of the halos in the
simulations, and study the differences between them. We find that both
radiative variations of the models show degeneracies between their processes
and at least two of the three parameters defining the Symmetron model.Comment: 9 pages, 4 figures, matches version accepted to A&
Baryon acoustic oscillations reconstruction using convolutional neural networks
We propose a new scheme to reconstruct the baryon acoustic oscillations (BAO) signal, which contains key cosmological information, based on deep convolutional neural networks (CNN). Trained with almost no fine tuning, the network can recover large-scale modes accurately in the test set: the correlation coefficient between the true and reconstructed initial conditions reaches 90% at k †0.2âhMpcâ1, which can lead to significant improvements of the BAO signal-to-noise ratio down to k â 0.4âhMpcâ1. Since this new scheme is based on the configuration-space density field in sub-boxes, it is local and less affected by survey boundaries than the standard reconstruction method, as our tests confirm. We find that the network trained in one cosmology is able to reconstruct BAO peaks in the others, i.e. recovering information lost to non-linearity independent of cosmology. The accuracy of recovered BAO peak positions is far less than that caused by the difference in the cosmology models for training and testing, suggesting that different models can be distinguished efficiently in our scheme. It is very promising that Our scheme provides a different new way to extract the cosmological information from the ongoing and future large galaxy surveys
Characterizing the contaminating distance distribution for Bayesian supernova cosmology
Measurements of the equation of state of dark energy from surveys of
thousands of Type Ia Supernovae (SNe Ia) will be limited by spectroscopic
follow-up and must therefore rely on photometric identification, increasing the
chance that the sample is contaminated by Core Collapse Supernovae (CC SNe).
Bayesian methods for supernova cosmology can remove contamination bias while
maintaining high statistical precision but are sensitive to the choice of
parameterization of the contaminating distance distribution. We use simulations
to investigate the form of the contaminating distribution and its dependence on
the absolute magnitudes, light curve shapes, colors, extinction, and redshifts
of core collapse supernovae. We find that the CC luminosity function dominates
the distance distribution function, but its shape is increasingly distorted as
the redshift increases and more CC SNe fall below the survey magnitude limit.
The shapes and colors of the CC light curves generally shift the distance
distribution, and their effect on the CC distances is correlated. We compare
the simulated distances to the first year results of the SDSS-II SN survey and
find that the SDSS distance distributions can be reproduced with simulated CC
SNe that are ~1 mag fainter than the standard Richardson et al. (2002)
luminosity functions, which do not produce a good fit. To exploit the full
power of the Bayesian parameter estimation method, parameterization of the
contaminating distribution should be guided by the current knowledge of the CC
luminosity functions, coupled with the effects of the survey selection and
magnitude-limit, and allow for systematic shifts caused by the parameters of
the distance fit.Comment: 17 pages, 5 figures; accepted for publication in the Astrophysical
Journa
Photometric Supernova Cosmology with BEAMS and SDSS-II
Supernova cosmology without spectroscopic confirmation is an exciting new
frontier which we address here with the Bayesian Estimation Applied to Multiple
Species (BEAMS) algorithm and the full three years of data from the Sloan
Digital Sky Survey II Supernova Survey (SDSS-II SN). BEAMS is a Bayesian
framework for using data from multiple species in statistical inference when
one has the probability that each data point belongs to a given species,
corresponding in this context to different types of supernovae with their
probabilities derived from their multi-band lightcurves. We run the BEAMS
algorithm on both Gaussian and more realistic SNANA simulations with of order
10^4 supernovae, testing the algorithm against various pitfalls one might
expect in the new and somewhat uncharted territory of photometric supernova
cosmology. We compare the performance of BEAMS to that of both mock
spectroscopic surveys and photometric samples which have been cut using typical
selection criteria. The latter typically are either biased due to contamination
or have significantly larger contours in the cosmological parameters due to
small data-sets. We then apply BEAMS to the 792 SDSS-II photometric supernovae
with host spectroscopic redshifts. In this case, BEAMS reduces the area of the
(\Omega_m,\Omega_\Lambda) contours by a factor of three relative to the case
where only spectroscopically confirmed data are used (297 supernovae). In the
case of flatness, the constraints obtained on the matter density applying BEAMS
to the photometric SDSS-II data are \Omega_m(BEAMS)=0.194\pm0.07. This
illustrates the potential power of BEAMS for future large photometric supernova
surveys such as LSST.Comment: 25 pages, 15 figures, submitted to Ap
ORIGAMI: Delineating Halos using Phase-Space Folds
We present the ORIGAMI method of identifying structures, particularly halos,
in cosmological N-body simulations. Structure formation can be thought of as
the folding of an initially flat three-dimensional manifold in six-dimensional
phase space. ORIGAMI finds the outer folds that delineate these structures.
Halo particles are identified as those that have undergone shell-crossing along
3 orthogonal axes, providing a dynamical definition of halo regions that is
independent of density. ORIGAMI also identifies other morphological structures:
particles that have undergone shell-crossing along 2, 1, or 0 orthogonal axes
correspond to filaments, walls, and voids respectively. We compare this method
to a standard Friends-of-Friends halo-finding algorithm and find that ORIGAMI
halos are somewhat larger, more diffuse, and less spherical, though the global
properties of ORIGAMI halos are in good agreement with other modern
halo-finding algorithms.Comment: 11 pages, 14 figures; matches version accepted to Ap
Straightening the Density-Displacement Relation with a Logarithmic Transform
We investigate the use of a logarithmic density variable in estimating the
Lagrangian displacement field, motivated by the success of a logarithmic
transformation in restoring information to the matter power spectrum. The
logarithmic relation is an extension of the linear relation, motivated by the
continuity equation, in which the density field is assumed to be proportional
to the divergence of the displacement field; we compare the linear and
logarithmic relations by measuring both of these fields directly in a
cosmological N-body simulation. The relative success of the logarithmic and
linear relations depends on the scale at which the density field is smoothed.
Thus we explore several ways of measuring the density field, including
Cloud-In-Cell smoothing, adaptive smoothing, and the (scale-independent)
Delaunay tessellation, and we use both a Fourier space and a geometrical
tessellation approach to measuring the divergence. We find that the relation
between the divergence of the displacement field and the density is
significantly tighter with a logarithmic density variable, especially at low
redshifts and for very small (~2 Mpc/h) smoothing scales. We find that the
grid-based methods are more reliable than the tessellation-based method of
calculating both the density and the divergence fields, though in both cases
the logarithmic relation works better in the appropriate regime, which
corresponds to nonlinear scales for the grid-based methods and low densities
for the tessellation-based method.Comment: 6 pages, 3 figures, accepted to Ap
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).Comment: 27 interesting pages, 20 beautiful figures, and 4 informative tables
accepted for publication in MNRAS. The high-resolution version of the paper
as well as all the test cases and analysis can be found at the web site
http://popia.ft.uam.es/HaloesGoingMA