Source author record

Francisco Prada

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

104works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

104 published item(s)

preprint2025arXiv

Galaxy-Multiplet Clustering from DESI DR2

We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

preprint2023arXiv

The Target-selection Pipeline for the Dark Energy Spectroscopic Instrument

In 2021 May, the Dark Energy Spectroscopic Instrument (DESI) began a 5 yr survey of approximately 50 million total extragalactic and Galactic targets. The primary DESI dark-time targets are emission line galaxies (ELGs), luminous red galaxies (LRGs) and quasars (QSOs). In bright time, DESI will focus on two surveys known as the Bright Galaxy Survey (BGS) and the Milky Way Survey (MWS). DESI also observes a selection of "secondary" targets for bespoke science goals. This paper gives an overview of the publicly available pipeline (desitarget) used to process targets for DESI observations. Highlights include details of the different DESI survey targeting phases, the targeting ID (TARGETID) used to define unique targets, the bitmasks used to indicate a particular type of target, the data model and structure of DESI targeting files, and examples of how to access and use the desitarget code base. This paper will also describe "supporting" DESI target classes, such as standard stars, sky locations, and random catalogs that mimic the angular selection function of DESI targets. The DESI target selection pipeline is complex and sizable; this paper attempts to summarize the most salient information required to understand and work with DESI targeting data.

preprint2022arXiv

Fast full N-body simulations of generic modified gravity: conformal coupling models

We present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for three classes of MG models with conformal coupling terms: the $f(R)$ gravity, symmetron and coupled quintessence models. Derived from the parallel particle-mesh code GLAM, MG-GLAM incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For $f(R)$ gravity, we have included new variants to diversify the model behaviour, and we have tailored the relaxation algorithms to these to maintain high computational efficiency. In a companion paper, we describe versions of this code developed for derivative coupling MG models, including the Vainshtein- and K-mouflage-type models. MG-GLAM can model the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous speedup of a factor of more than a hundred compared with earlier $N$-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. MG-GLAM is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and the production of accurate emulators for ongoing and future galaxy surveys.

preprint2022arXiv

Towards an accurate model of small-scale redshift-space distortions in modified gravity

The coming generation of galaxy surveys will provide measurements of galaxy clustering with unprecedented accuracy and data size, which will allow us to test cosmological models at much higher precision than achievable previously. This means that we must have more accurate theoretical predictions to compare with future observational data. As a first step towards more accurate modelling of the redshift space distortions (RSD) of small-scale galaxy clustering in modified gravity (MG) cosmologies, we investigate the validity of the so-called Skew-T (ST) probability distribution function (PDF) of halo pairwise peculiar velocities in these models. We show that, combined with the streaming model of RSD, the ST PDF substantially improves the small-scale predictions by incorporating skewness and kurtosis, for both $Λ$CDM and two leading MG models: $f(R)$ gravity and the DGP braneworld model. The ST model reproduces the velocity PDF and redshift-space halo clustering measured from MG $N$-body simulations very well down to $\sim 5 \, h^{-1}\mathrm{Mpc}$. In particular, we investigate the enhancements of halo pairwise velocity moments with respect to $Λ$CDM for a larger range of MG variants than previous works, and present simple explanations to the behaviours observed. By performing a simple Fisher analysis, we find a significnat increase in constraining power to detect modifications of General Relativity by introducing small-scale information in the RSD analyses.

preprint2021arXiv

Fast full $N$-body simulations of generic modified gravity: derivative coupling models

We present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological simulations for two broad classes of MG models with derivative coupling terms -- the Vainshtein- and Kmouflage-type models -- which respectively features the Vainshtein and Kmouflage screening mechanism. Derived from the parallel particle-mesh code GLAM, MG-GLAM incorporates an efficient multigrid relaxation technique to solve the characteristic nonlinear partial differential equations of these models. For Kmouflage, we have proposed a new algorithm for the relaxation solver, and run the first simulations of the model to understand its cosmological behaviour. In a companion paper, we describe versions of this code developed for conformally-coupled MG models, including several variants of $f(R)$ gravity, the symmetron model and coupled quintessence. Altogether, MG-GLAM has so far implemented the prototypes for most MG models of interest, and is broad and versatile. The code is highly optimised, with a tremendous (over two orders of magnitude) speedup when comparing its running time with earlier $N$-body codes, while still giving accurate predictions of the matter power spectrum and dark matter halo abundance. MG-GLAM is ideal for the generation of large numbers of MG simulations that can be used in the construction of mock galaxy catalogues and accurate emulators for ongoing and future galaxy surveys.

preprint2021arXiv

Final Targeting Strategy for the SDSS-IV APOGEE-2N Survey

APOGEE-2 is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemo-dynamical mapping of the Milky Way Galaxy. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs" competitively awarded to SDSS-IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5-year expansion of the survey, known as the Bright Time Extension, made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The Bright Time Extension permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new datasets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, Santana et al. (submitted), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

preprint2021arXiv

Installation of the Dark Energy Spectroscopic Instrument at the Mayall 4-meter telescope

The Dark Energy Spectroscopic Instrument (DESI) is a Stage IV ground-based dark energy experiment that will measure the expansion history of the Universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasars over 14000 square degrees will be measured during the life of the experiment. We describe the installation of the major elements of the instrument at the Mayall 4m telescope, completed in late 2019. The previous prime focus corrector, spider vanes, and upper rings were removed from the Mayall's Serrurier truss and replaced with the newly-constructed DESI ring, vanes, cage, hexapod, and optical corrector. The new corrector was optically aligned with the primary mirror using a laser tracker system. The DESI focal plane system was integrated to the corrector, with each of its ten 500-fiber-positioner petal segments installed using custom installation hardware and the laser tracker. Ten DESI spectrographs with 30 cryostats were installed in a newly assembled clean room in the Large Coude Room. The ten cables carrying 5000 optical fibers from the positioners in the focal plane were routed down the telescope through cable wraps at the declination and hour angle axes, and their integral slitheads were integrated with the ten spectrographs. The fiber view camera assembly was installed to the Mayall's primary mirror cell. Servers for the instrument control system replaced existing computer equipment. The fully integrated instrument has been commissioned and is ready to start its operations phase.

preprint2021arXiv

Performance of Kitt Peak's Mayall 4-meter Telescope During DESI Commissioning

In preparation for the Dark Energy Spectroscopic Instrument (DESI), a new top end was installed on the Mayall 4-meter telescope at Kitt Peak National Observatory. The refurbished telescope and the DESI instrument were successfully commissioned on sky between 2019 October and 2020 March. Here we describe the pointing, tracking and imaging performance of the Mayall telescope equipped with its new DESI prime focus corrector, as measured by six guider cameras sampling the outer edge of DESI's focal plane. Analyzing ~500,000 guider images acquired during commissioning, we find a median delivered image FWHM of 1.1 arcseconds (in the r-band at 650 nm), with the distribution extending to a best-case value of ~0.6 arcseconds. The point spread function is well characterized by a Moffat profile with a power-law index of $β$ ~ 3.5 and little dependence of $β$ on FWHM. The shape and size of the PSF delivered by the new corrector at a field angle of 1.57 degrees are very similar to those measured with the old Mayall corrector on axis. We also find that the Mayall achieves excellent pointing accuracy (several arcseconds RMS) and minimal open-loop tracking drift (< 1 milliarcsecond per second), improvements on the telecope's pre-DESI performance. In the future, employing DESI's active focus adjustment capabilities will likely further improve the Mayall/DESI delivered image quality.

preprint2021arXiv

The DESI Sky Continuum Monitor System

The Dark Energy Spectroscopic Instrument (DESI) is an ongoing spectroscopic survey to measure the dark energy equation of state to unprecedented precision. We describe the DESI Sky Continuum Monitor System, which tracks the night sky brightness as part of a system that dynamically adjusts the spectroscopic exposure time to produce more uniform data quality and to maximize observing efficiency. The DESI dynamic exposure time calculator (ETC) will combine sky brightness measurements from the Sky Monitor with data from the guider system to calculate the exposure time to achieve uniform signal-to-noise ratio (SNR) in the spectra under various observing conditions. The DESI design includes 20 sky fibers, and these are split between two identical Sky Monitor units to provide redundancy. Each Sky Monitor unit uses an SBIG STXL-6303e CCD camera and supports an eight-position filter wheel. Both units have been completed and delivered to the Mayall Telescope at the Kitt Peak National Observatory. Commissioning results show that the Sky Monitor delivers the required performance necessary for the ETC.

preprint2020arXiv

Baryon Acoustic Oscillations in the projected cross-correlation function between the eBOSS DR16 quasars and photometric galaxies from the DESI Legacy Imaging Surveys

We search for the Baryon Acoustic Oscillations in the projected cross-correlation function binned into transverse comoving radius between the SDSS-IV DR16 eBOSS quasars and a dense photometric sample of galaxies selected from the DESI Legacy Imaging Surveys. We estimate the density of the photometric sample of galaxies in this redshift range to be about 2900 deg$^{-2}$, which is deeper than the official DESI ELG selection, and the density of the spectroscopic sample is about 20 deg$^{-2}$. In order to mitigate the systematics related to the use of different imaging surveys close to the detection limit, we use a neural network approach that accounts for complex dependencies between the imaging attributes and the observed galaxy density. We find that we are limited by the depth of the imaging surveys which affects the density and purity of the photometric sample and its overlap in redshift with the quasar sample, which thus affects the performance of the method. When cross-correlating the photometric galaxies with quasars in $0.6 \leq z \leq 1.2$, the cross-correlation function can provide better constraints on the comoving angular distance, $D_{\rm M}$ (6\% precision) compared to the constraint on the spherically-averaged distance $D_{\rm V}$ (9\% precision) obtained from the auto-correlation. Although not yet competitive, this technique will benefit from the arrival of deeper photometric data from upcoming surveys which will enable it to go beyond the current limitations we have identified in this work.

preprint2020arXiv

Building a digital twin of a luminous red galaxy spectroscopic survey: galaxy properties and clustering covariance

Luminous red galaxies (LRGs) are one of the key tracers of the large-scale structure of the Universe used by galaxy surveys. Hence, it is important to make accurate predictions for their properties and clustering, including the errors on these statistics. Here, we describe a novel technique which uses the semi-analytical model of galaxy formation {\sc Galform}, embedded in the high-resolution $N$-body Planck-Millennium simulation, to populate a thousand halo catalogues generated using the Parallel-PM $N$-body {\sc glam} code. Our hybrid scheme allows us to make clustering predictions on scales that cannot be modelled in the original $N$-body simulation. LRGs are selected in the redshift range $z=0.6-1$ from the {\sc Galform} output using similar colour-magnitude cuts in the $r$, $z$ and $W1$ bands to those that will be applied in the Dark Energy Spectroscopic Instrument (DESI) survey. We find that the LRG-halo connection is non-trivial, leading to the prediction of a non-standard halo occupation distribution; in particular, the occupation of central galaxies does not reach unity for the most massive haloes, and drops with increasing mass. The {\sc glam} catalogues reproduce the abundance and clustering of the LRGs predicted by {\sc Galform}, and show good agreement with recent measurements of the clustering of DESI-like LRGs using photometric redshifts. We use the \glam{} mocks to compute the covariance matrices for the two-point correlation function and power spectrum of the LRGs and their background dark matter density field, revealing important differences. We also make predictions for the linear-growth rate and the baryon acoustic oscillations distances at $z=0.6$, $0.74$ and $0.93$. All DESI-like LRG catalogues are made publicly available.

preprint2020arXiv

Imaging Systematics and Clustering of DESI Main Targets

We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG) and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits $r_0 = 7.78 \pm 0.26$ $h^{-1}$Mpc, $γ= 1.98 \pm 0.02$ for LRGs and $r_0 = 5.45 \pm 0.1$ $h^{-1}$Mpc, $γ= 1.54 \pm 0.01$ for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale dependent biases in the weakly nonlinear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer $dN/dz$ and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.

preprint2020arXiv

The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: GLAM-QPM mock galaxy catalogs for the Emission Line Galaxy Sample

We present 2000 mock galaxy catalogs for the analysis of baryon acoustic oscillations in the Emission Line Galaxy (ELG) sample of the Extended Baryon Oscillation Spectroscopic Survey Data Release 16 (eBOSS DR16). Each mock catalog has a number density of $6.7 \times 10^{-4} h^3 \rm Mpc^{-3}$, covering a redshift range from 0.6 to 1.1. The mocks are calibrated to small-scale eBOSS ELG clustering measurements at scales of around 10 $h^{-1}$Mpc. The mock catalogs are generated using a combination of GaLAxy Mocks (GLAM) simulations and the Quick Particle-Mesh (QPM) method. GLAM simulations are used to generate the density field, which is then assigned dark matter halos using the QPM method. Halos are populated with galaxies using a halo occupation distribution (HOD). The resulting mocks match the survey geometry and selection function of the data, and have slightly higher number density which allows room for systematic analysis. The large-scale clustering of mocks at the baryon acoustic oscillation (BAO) scale is consistent with data and we present the correlation matrix of the mocks.

preprint2020arXiv

The Local Universe from Calar Alto (LUCA)

LUCA (Local Universe from Calar Alto) was conceived as a new generation science program for the Calar Alto Observatory. It proposed the construction of a large Integral Field Unit (IFU) spectrograph with six thousands optical fibers (IFU-6000) at the CAHA 3.5-m telescope to map our universe neighborhood in 3D with an unprecedented spatial resolution. Two galaxy samples were defined to map the local universe: (i) 102 galaxies in the Local Volume, out to 11 Mpc, and (ii) 218 galaxies in the Virgo cluster. A complementary project was developed to map the three largest galaxies in the local universe, M31, M33, and M101 with the Schmidt telescope. In this white paper we describe the LUCA Project science justification and survey strategy, as well as a technical description of the IFU-6000 instrument, its performance, and design.

preprint2020arXiv

TXS 2116$-$077: A Gamma-ray Emitting Relativistic Jet Hosted in a Galaxy Merger

What triggers collimated relativistic outflows or jets, from the centers of galaxies remains a fundamental question in astrophysics. The merging of two galaxies has been proposed to realize the conditions to successfully launch and drive such jets into the intergalactic medium. However, evidences for the operation of this mechanism are scarce. Here we report the first unambiguous detection of an ongoing merger of a narrow-line Seyfert 1 galaxy, TXS 2116$-$077, hosting a closely aligned, $γ$-ray emitting relativistic jet with a Seyfert 2 galaxy at a separation of $\sim$12 kpc, using the observations taken with 8.2 m Subaru telescope. Our subsequent followup observations with 10.4 m Gran Telescopio Canarias, 4.2 m William Herschel Telescope, and Chandra X-ray observatory have provided what is likely to be the first glimpse of the merging environment hosting a closely aligned relativistic jet. Our finding that the jet is considerably younger than the merger demonstrates that jet activity can be triggered by galaxy mergers and that $γ$-ray detected narrow-line Seyfert 1 galaxies represent the beginning phase of that activity. These results also highlight the crucial role of mergers in shaping the fate of galaxies in their cosmological evolution and are consistent with recent studies focused on the host galaxy imaging of this enigmatic class of active galactic nuclei.

preprint2020arXiv

White Paper on MAAT@GTC

MAAT is proposed as a visitor mirror-slicer optical system that will allow the OSIRIS spectrograph on the 10.4-m Gran telescopio CANARIAS (GTC) the capability to perform Integral Field Spectroscopy (IFS) over a seeing-limited FoV 14.20''x10'' with a slice width of 0.303''. MAAT@GTC will enhance the resolution power of OSIRIS by 1.6 times as compared to its 0.6'' wide long-slit. All the eleven OSIRIS grisms and volume-phase holographic gratings will be available to provide broad spectral coverage with moderate resolution (R=600 up to 4100) in the 3600 - 10000 Å wavelength range. MAAT unique observing capabilities will broaden its use to the needs of the GTC community to unveil the nature of most striking phenomena in the universe well beyond time-domain astronomy. The GTC equipped with OSIRIS+MAAT will also play a fundamental role in synergy with other facilities, some of them operating on the northern ORM at La Palma. This White Paper presents the different aspects of MAAT@GTC - including scientific and technical specifications, outstanding science cases, and an outline of the instrument concept.

preprint2019arXiv

Suppressing cosmic variance with paired-and-fixed cosmological simulations: average properties and covariances of dark matter clustering statistics

Making cosmological inferences from the observed galaxy clustering requires accurate predictions for the mean clustering statistics and their covariances. Those are affected by cosmic variance -- the statistical noise due to the finite number of harmonics. The cosmic variance can be suppressed by fixing the amplitudes of the harmonics instead of drawing them from a Gaussian distribution predicted by the inflation models. Initial realizations also can be generated in pairs with 180 degrees flipped phases to further reduce the variance. Here, we compare the consequences of using paired-and-fixed vs Gaussian initial conditions on the average dark matter clustering and covariance matrices predicted from N-body simulations. As in previous studies, we find no measurable differences between paired-and-fixed and Gaussian simulations for the average density distribution function, power spectrum and bispectrum. Yet, the covariances from paired-and-fixed simulations are suppressed in a complicated scale- and redshift-dependent way. The situation is particularly problematic on the scales of Baryon Acoustic Oscillations where the covariance matrix of the power spectrum is lower by only 20% compared to the Gaussian realizations, implying that there is not much of a reduction of the cosmic variance. The non-trivial suppression, combined with the fact that paired-and-fixed covariances are noisier than from Gaussian simulations, suggests that there is no path towards obtaining accurate covariance matrices from paired-and-fixed simulations. Because the covariances are crucial for the observational estimates of galaxy clustering statistics and cosmological parameters, paired-and-fixed simulations, though useful for some applications, cannot be used for the production of mock galaxy catalogs.

preprint2016arXiv

Clustering properties of $g$-selected galaxies at $z\sim0.8$

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

preprint2016arXiv

Galaxy Three-Point Correlation Functions and Halo/Subhalo Models

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

preprint2016arXiv

Gravitational effects on the Higgs field within the Solar System

The Higgs mechanism predicts, apart from the existence of a new scalar boson, the presence of a constant Higgs field that permeates all of space. The vacuum expectation value (VEV) of this field is affected by quantum corrections which are mainly generated by the self-interactions and couplings of the Higgs field to gauge bosons and heavy quarks. In this work we show that gravity can affect, in a non-trivial way, these quantum corrections through the finite parts of the one-loop contributions to the effective potential. In particular, we consider the corrections generated by the Standard Model Higgs self-interactions in slowly-varying weak gravitational backgrounds. The obtained results amount to the existence of non-negligible inhomogeneities in the Higgs VEV. Such inhomogeneities translate into spatial variations of the particle masses, and in particular of the proton-to-electron mass ratio. We find that these Higgs perturbations in our Solar System are controlled by the Eddington parameter, and are absent in pure General Relativity. Yet, they may be present in modified gravity theories. This predicted effect may be constrained by atomic clocks or high-resolution spectroscopic measurements, which could allow to improve current limits on modifications of Einstein's gravity.

preprint2016arXiv

Higgs effective potential in a perturbed Robertson-Walker background

We calculate the one-loop effective potential of a scalar field in a Robertson-Walker background with scalar metric perturbations. A complete set of orthonormal solutions of the perturbed equations is obtained by using the adiabatic approximation for comoving observers. After analyzing the problem of renormalization in inhomogeneous backgrounds, we get the explicit contribution of metric perturbations to the effective potential. We apply these results to the Standard Model Higgs field and evaluate the effects of metric perturbations on the Higgs mass and on its vacuum expectation value. Space-time variations are found, which are proportional to the gravitational slip parameter, with a typical amplitude of the order of $Δϕ/ϕ\simeq 10^{-11}$ on cosmological scales. We also discuss possible astrophysical signatures in the Solar System and in the Milky Way that could open new possibilities to explore the symmetry breaking sector of the electroweak interactions.

preprint2016arXiv

Modelling galaxy clustering: halo occupation distribution versus subhalo matching

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $χ^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

preprint2016arXiv

MultiDark simulations: the story of dark matter halo concentrations and density profiles

Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations.

preprint2016arXiv

MultiDarkLens Simulations: weak lensing light-cones and data base presentation

In this paper we present a large database of weak lensing light cones constructed using different snapshots from the Big MultiDark simulation (BigMDPL). The ray-tracing through different multiple planes has been performed with the GLAMER code accounting both for single source redshifts and for sources distributed along the cosmic time. This first paper presents weak lensing forecasts and results according to the geometry of the VIPERS-W1 and VIPERS-W4 field of view. Additional fields will be available on our database and new ones can be run upon request. Our database also contains some tools for lensing analysis. In this paper we present results for convergence power spectra, one point and high order weak lensing statistics useful for forecasts and for cosmological studies. Covariance matrices have also been computed for the different realisations of the W1 and W4 fields. In addition we compute also galaxy-shear and projected density contrasts for different halo masses at two lens redshifts according to the CFHTLS source redshift distribution both using stacking and cross-correlation techniques, finding very good agreement.

preprint2016arXiv

The angular power spectrum of the diffuse gamma-ray emission as measured by the Fermi Large Area Telescope and constraints on its Dark Matter interpretation

The isotropic gamma-ray background arises from the contribution of unresolved sources, including members of confirmed source classes and proposed gamma-ray emitters such as the radiation induced by dark matter annihilation and decay. Clues about the properties of the contributing sources are imprinted in the anisotropy characteristics of the gamma-ray background. We use 81 months of Pass 7 Reprocessed data from the Fermi Large Area Telescope to perform a measurement of the anisotropy angular power spectrum of the gamma-ray background. We analyze energies between 0.5 and 500 GeV, extending the range considered in the previous measurement based on 22 months of data. We also compute, for the first time, the cross-correlation angular power spectrum between different energy bins. We find that the derived angular spectra are compatible with being Poissonian, i.e. constant in multipole. Moreover, the energy dependence of the anisotropy suggests that the signal is due to two populations of sources, contributing, respectively, below and above 2 GeV. Finally, using data from state-of-the-art numerical simulations to model the dark matter distribution, we constrain the contribution from dark matter annihilation and decay in Galactic and extragalactic structures to the measured anisotropy. These constraints are competitive with those that can be derived from the average intensity of the isotropic gamma-ray background. Data are available at https://www-glast.stanford.edu/pub_data/552.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourier-space

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z < 0.75$. We divide this dataset into three (overlapping) redshift bins with the effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as $\sim 1000$ MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately. We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and $2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high ($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and $0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of $8σ$ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within $Λ$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: combining correlated Gaussian posterior distributions

The cosmological information contained in anisotropic galaxy clustering measurements can often be compressed into a small number of parameters whose posterior distribution is well described by a Gaussian. We present a general methodology to combine these estimates into a single set of consensus constraints that encode the total information of the individual measurements, taking into account the full covariance between the different methods. We illustrate this technique by applying it to combine the results obtained from different clustering analyses, including measurements of the signature of baryon acoustic oscillations (BAO) and redshift-space distortions (RSD), based on a set of mock catalogues of the final SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our results show that the region of the parameter space allowed by the consensus constraints is smaller than that of the individual methods, highlighting the importance of performing multiple analyses on galaxy surveys even when the measurements are highly correlated. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The methodology presented here is used in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the configuration-space clustering wedges

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $Λ$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $Ω_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_ν < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=Ω_{\rm m}(z)^γ$ and obtain $γ= 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $γ=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $fσ_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $Λ$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: double-probe measurements from BOSS galaxy clustering \& Planck data -- towards an analysis without informative priors

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $Ω_b h^2$, $n_s$, $log(A_s)$, $Ω_k$, $H(z)$, $D_A(z)$, $f(z)σ_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{Ω_m H_0^2}\,r(z_*)$, and $l_a=πr(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $Ω_b h^2$ and $n_s$. Using our double-probe results, we obtain $Ω_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $σ_8=0.806\pm0.014$ assuming $Λ$CDM; and $Ω_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $Λ$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $Σm_ν<0.10/0.22$ (68\%/95\%) assuming $Λ$CDM and $Σm_ν<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

preprint2016arXiv

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: BAO measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

[abridged] We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains contains 361 762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$; the CMASS sample consists of 777 202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. We extract the BAO peak position from the monopole power spectrum moment, $α_0$, and from the $μ^2$ moment, $α_2$, where $μ$ is the cosine of the angle to the line-of-sight. The $μ^2$-moment provides equivalent information to that available in the quadrupole but is simpler to analyse. After applying a reconstruction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO peak position in the monopole and $μ^2$-moment, which are related to radial and angular shifts in scale. We report $H(z_{\rm LOWZ})r_s(z_d)=(11.60\pm0.60)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.66\pm0.16$ with a cross-correlation coefficient of $r_{HD_A}=0.41$, for the LOWZ sample; and $H(z_{\rm CMASS})r_s(z_d)=(14.56\pm0.37)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.42\pm0.13$ with a cross-correlation coefficient of $r_{HD_A}=0.47$, for the CMASS sample. We combine these results with the measurements of the BAO peak position in the monopole and quadrupole correlation function of the same dataset \citep[][companion paper]{Cuestaetal2015} and report the consensus values: $H(z_{\rm LOWZ})r_s(z_d)=(11.63\pm0.69)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.67\pm0.15$ with $r_{HD_A}=0.35$ for the LOWZ sample; $H(z_{\rm CMASS})r_s(z_d)=(14.67\pm0.42)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.47\pm0.12$ with $r_{HD_A}=0.52$ for the CMASS sample.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: mock galaxy catalogues for the BOSS Final Data Release

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $σ$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $Λ$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$σ$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$σ$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies

We measure and analyse the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm lowz}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm cmass}=0.57$. From the power spectrum monopole and quadrupole moments around the LOS, we measure the growth of structure parameter $f$ times the amplitude of dark matter density fluctuations $σ_8$ by modeling the Redshift-Space Distortion signal. When the geometrical Alcock-Paczynski effect is also constrained from the same data, we find joint constraints on $fσ_8$, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon $D_A(z)/r_s(z_d)$. We find $f(z_{\rm lowz})σ_8(z_{\rm lowz})=0.394\pm0.062$, $D_A(z_{\rm lowz})/r_s(z_d)=6.35\pm0.19$, $H(z_{\rm lowz})r_s(z_d)=(11.41\pm 0.56)\,{10^3\rm km}s^{-1}$ for the LOWZ sample, and $f(z_{\rm cmass})σ_8(z_{\rm cmass})=0.444\pm0.038$, $D_A(z_{\rm cmass})/r_s(z_d)=9.42\pm0.15$, $H(z_{\rm cmass})r_s(z_d)=(13.92 \pm 0.44)\, {10^3\rm km}s^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 measurements. Assuming the Hubble parameter and angular distance parameter are fixed at fiducial $Λ$CDM values, we find $f(z_{\rm lowz})σ_8(z_{\rm lowz})=0.485\pm0.044$ and $f(z_{\rm cmass})σ_8(z_{\rm cmass})=0.436\pm0.022$ for the LOWZ and CMASS samples, respectively.

preprint2016arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: single-probe measurements from CMASS anisotropic galaxy clustering

With the largest spectroscopic galaxy survey volume drawn from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), we can extract cosmological constraints from the measurements of redshift and geometric distortions at quasi-linear scales (e.g. above 50 $h^{-1}$Mpc). We analyze the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS galaxy sample, at the effective redshift $z=0.59$, to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)σ_8(z)$, and the physical matter density $Ω_mh^2$. We obtain robust measurements by including a polynomial as the model for the systematic errors, and find it works very well against the systematic effects, e.g., ones induced by stars and seeing. We provide accurate measurements $\{D_A(0.59)r_{s,fid}/r_s$ $\rm Mpc$, $H(0.59)r_s/r_{s,fid}$ $km s^{-1} Mpc^{-1}$, $f(0.59)σ_8(0.59)$, $Ω_m h^2\}$ = $\{1427\pm26$, $97.3\pm3.3$, $0.488 \pm 0.060$, $0.135\pm0.016\}$, where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$ Mpc is the sound scale of the fiducial cosmology used in this study. The parameters which are not well constrained by our galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets, e.g., cosmic microwave background (CMB), are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e., CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. The uncertainty on the dark energy equation of state parameter, $w$, from CMB+CMASS is about 8 per cent. The uncertainty on the curvature fraction, $Ω_k$, is 0.3 per cent. We do not find deviation from flat $Λ$CDM.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2016arXiv

The evolution of the [OII], Hβ and [OIII] emission-line luminosity functions over the last nine billions years

Emission line galaxies are one of the main tracers of the large-scale structure to be targeted by the next-generation dark energy surveys. To provide a better understanding of the properties and statistics of these galaxies, we have collected spectroscopic data from the VVDS and DEEP2 deep surveys and estimated the galaxy luminosity functions (LFs) of three distinct emission lines, [OII] ($0.5 < z < 1.3$), Hβ ($0.3 < z < 0.8$) and [OIII] ($0.3 < z < 0.8$). Our measurements are based on 35,639 emission line galaxies and cover a volume of $\sim10^7$Mpc$^3$. We present the first measurement of the Hβ LF at these redshifts. We have also compiled LFs from the literature that were based on independent data or covered different redshift ranges, and we fit the entire set over the whole redshift range with analytic Schechter and Saunders models, assuming a natural redshift dependence of the parameters. We find that the characteristic luminosity ($L_*$) and density ($ϕ_*$) of all LFs increase with redshift. Using the Schechter model over the redshift ranges considered, we find that, for [OII] emitters, the characteristic luminosity $L_*(z=0.5)=3.2\times10^{41}$ erg/s increases by a factor of $2.7 \pm 0.2$ from z=0.5 to 1.3; for Hβ emitters $L_*(z=0.3)=1.3\times10^{41}$ erg/s increases by a factor of $2.0 \pm 0.2$ from z=0.3 to 0.8; and for [OIII] emitters $L_*(z=0.3)=7.3\times10^{41}$ erg/s increases by a factor of $3.5 \pm 0.4$ from z=0.3 to 0.8.

preprint2016arXiv

The extended Baryon Oscillation Spectroscopic Survey (eBOSS): a cosmological forecast

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$α$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$α$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $fσ_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $σ(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

preprint2016arXiv

The XMM Cluster Survey: The Halo Occupation Number of BOSS galaxies in X-ray clusters

We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-selected galaxy clusters drawn from the XMM Cluster Survey (XCS). Halo masses are determined for each galaxy cluster based on X-ray temperature measurements, and range between ${\rm log_{10}} (M_{180}/M_{\odot}) = 13-15$. Our directly measured HODs are consistent with the HOD-model fits inferred via the galaxy-clustering analyses of Parejko et al. for the BOSS LOWZ sample and White et al. for the BOSS CMASS sample. Under the simplifying assumption that the other parameters that describe the HOD hold the values measured by these authors, we have determined a best-fit alpha-index of 0.91$\pm$0.08 and $1.27^{+0.03}_{-0.04}$ for the CMASS and LOWZ HOD, respectively. These alpha-index values are consistent with those measured by White et al. and Parejko et al. In summary, our study provides independent support for the HOD models assumed during the development of the BOSS mock-galaxy catalogues that have subsequently been used to derive BOSS cosmological constraints.

preprint2015arXiv

A Phenomenological Model for the Intracluster Medium that matches X-ray and Sunyaev-Zel'dovich observations

Cosmological hydrodynamical simulations of galaxy clusters are still challenged to produce a model for the intracluster medium that matches all aspects of current X-ray and Sunyaev-Zel'dovich observations. To facilitate such comparisons with future simulations and to enable realistic cluster population studies for modeling e.g., non-thermal emission processes, we construct a phenomenological model for the intracluster medium that is based on a representative sample of observed X-ray clusters. We create a mock galaxy cluster catalog based on the large collisionless N-body simulation MultiDark, by assigning our gas density model to each dark matter cluster halo. Our clusters are classified as cool-core and non cool-core according to a dynamical disturbance parameter. We demonstrate that our gas model matches the various observed Sunyaev-Zel'dovich and X-ray scaling relations as well as the X-ray luminosity function, thus enabling to build a reliable mock catalog for present surveys and forecasts for future experiments. In a companion paper, we apply our catalogs to calculate non-thermal radio and gamma-ray emission of galaxy clusters. We make our cosmologically complete multi-frequency mock catalogs for the (non-)thermal cluster emission at different redshifts publicly and freely available online through the MultiDark database (www.multidark.org).

preprint2015arXiv

Bayesian redshift-space distortions correction from galaxy redshift surveys

We present a Bayesian reconstruction method which maps a galaxy distribution from redshift-space to real-space inferring the distances of the individual galaxies. The method is based on sampling density fields assuming a lognormal prior with a likelihood given by the negative binomial distribution function modelling stochastic bias. We assume a deterministic bias given by a power law relating the dark matter density field to the expected halo or galaxy field. Coherent redshift-space distortions are corrected in a Gibbs-sampling procedure by moving the galaxies from redshift-space to real-space according to the peculiar motions derived from the recovered density field using linear theory with the option to include tidal field corrections from second order Lagrangian perturbation theory. The virialised distortions are corrected by sampling candidate real-space positions (being in the neighbourhood of the observations along the line of sight), which are compatible with the bulk flow corrected redshift-space position adding a random dispersion term in high density collapsed regions. The latter are defined according to the eigenvalues of the Hessian. This approach presents an alternative method to estimate the distances to galaxies using the three dimensional spatial information, and assuming isotropy. Hence the number of applications is very broad. In this work we show the potential of this method to constrain the growth rate up to $k$ ~ 0.3 $h$ Mpc$^{-1}$. Furthermore it could be useful to correct for photo-metric redshift errors, and to obtain improved BAO reconstructions.

preprint2015arXiv

Building a Better Understanding of the High Redshift BOSS Galaxies as Tools for Cosmology

We explore the bluer star-forming population of the Sloan Digital Sky Survey (SDSS) III/BOSS CMASS DR11 galaxies at $z>0.55$ to quantify their differences, in terms of redshift-space distortions and large-scale bias, with respect to the luminous red galaxy sample. We perform a qualitative analysis to understand the significance of these differences and whether we can model and reproduce them in mock catalogs. Specifically, we measure galaxy clustering in CMASS on small and intermediate scales ($r\lesssim 50\,h^{-1}$Mpc) by computing the two-point correlation function $-$ both projected and redshift-space $-$ of these galaxies, and a new statistic, $Σ(π)$, able to provide robust information about redshift-space distortions and large-scale bias. We interpret our clustering measurements by adopting a Halo Occupation Distribution (HOD) scheme that maps them onto high-resolution N-body cosmological simulations to produce suitable mock galaxy catalogs. The traditional HOD prescription can be applied to the red and the blue samples, independently, but this approach is unphysical since it allows the same mock galaxies to be either red or blue. To overcome this failure, we modify the standard formulation and infer the red and the blue mock catalogs directly from the full one, so that they are complementary and non-overlapping. This separation is performed by matching the observed CMASS red and blue galaxy fractions and produces reliable and accurate models.

preprint2015arXiv

Constraint on the time variation of the fine-structure constant with the SDSS-III/BOSS DR12 quasar sample

From the Sloan Digital Sky Survey (SDSS) Data Release 12, which covers the full Baryonic Oscillation Spectroscopic Survey (BOSS) footprint, we investigate the possible variation of the fine-structure constant over cosmological time-scales. We analyse the largest quasar sample considered so far in the literature, which contains 13175 spectra (10363 from SDSS-III/BOSS DR12 + 2812 from SDSS-II DR7) with redshift $z<\,$1. We apply the emission-line method on the [O III] doublet (4960, 5008 A) and obtain $Δα/α= \left(0.9 \pm 1.8\right)\times10^{-5}$ for the relative variation of the fine-structure constant. We also investigate the possible sources of systematics: misidentification of the lines, sky OH lines, H$\,β$ and broad line contamination, Gaussian and Voigt fitting profiles, optimal wavelength range for the Gaussian fits, chosen polynomial order for the continuum spectrum, signal-to-noise ratio and good quality of the fits. The uncertainty of the measurement is dominated by the sky subtraction. The results presented in this work, being systematics limited, have sufficient statistics to constrain robustly the variation of the fine-structure constant in redshift bins ($Δz\approx$ 0.06) over the last 7.9 Gyr. In addition, we study the [Ne III] doublet (3870, 3969 A) present in 462 quasar spectra and discuss the systematic effects on using these emission lines to constrain the fine-structure constant variation. Better constraints on $Δα/α $ ($<$10$^{-6}$) using the emission-line method would be possible with high-resolution spectroscopy and large galaxy/qso surveys.

preprint2015arXiv

Cosmological implications of baryon acoustic oscillation (BAO) measurements

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $Ω_m=0.301 \pm 0.008$ and curvature $Ω_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $χ^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

preprint2015arXiv

Galaxy clusters and groups in the ALHAMBRA Survey

We present a catalogue of 348 galaxy clusters and groups with $0.2<z<1.2$ selected in the 2.78 $deg^2$ ALHAMBRA Survey. The high precision of our photometric redshifts, close to $1\%$, and the wide spread of the seven ALHAMBRA pointings ensure that this catalogue has better mass sensitivity and is less affected by cosmic variance than comparable samples. The detection has been carried out with the Bayesian Cluster Finder (BCF), whose performance has been checked in ALHAMBRA-like light-cone mock catalogues. Great care has been taken to ensure that the observable properties of the mocks photometry accurately correspond to those of real catalogues. From our simulations, we expect to detect galaxy clusters and groups with both $70\%$ completeness and purity down to dark matter halo masses of $M_h\sim3\times10^{13}\rm M_{\odot}$ for $z<0.85$. Cluster redshifts are expected to be recovered with $\sim0.6\%$ precision for $z<1$. We also expect to measure cluster masses with $σ_{M_h|M^*_{CL}}\sim0.25-0.35\, dex$ precision down to $\sim3\times10^{13}\rm M_{\odot}$, masses which are $50\%$ smaller than those reached by similar work. We have compared these detections with previous optical, spectroscopic and X-rays work, finding an excellent agreement with the rates reported from the simulations. We have also explored the overall properties of these detections such as the presence of a colour-magnitude relation, the evolution of the photometric blue fraction and the clustering of these sources in the different ALHAMBRA fields. Despite the small numbers, we observe tentative evidence that, for a fixed stellar mass, the environment is playing a crucial role at lower redshifts (z$<$0.5).

preprint2015arXiv

Halo mass distribution reconstruction across the cosmic web

We study the relation between halo mass and its environment from a probabilistic perspective. We find that halo mass depends not only on local dark matter density, but also on non-local quantities such as the cosmic web environment and the halo-exclusion effect. Given these accurate relations, we have developed the HADRON-code (Halo mAss Distribution ReconstructiON), a technique which permits us to assign halo masses to a distribution of haloes in three-dimensional space. This can be applied to the fast production of mock galaxy catalogues, by assigning halo masses, and reproducing accurately the bias for different mass cuts. The resulting clustering of the halo populations agree well with that drawn from the BigMultiDark $N$-body simulation: the power spectra are within 1-$σ$ up to scales of $k=0.2\,h\,{\rm Mpc}^{-1}$, when using augmented Lagrangian perturbation theory based mock catalogues. Only the most massive haloes show a larger deviation. For these, we find evidence of the halo-exclusion effect. A clear improvement is achieved when assigning the highest masses to haloes with a minimum distance separation. We also compute the 2- and 3-point correlation functions, and find an excellent agreement with $N$-body results. Our work represents a quantitative application of the cosmic web classification. It can have further interesting applications in the multi-tracer analysis of the large-scale structure for future galaxy surveys.

preprint2015arXiv

nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

preprint2015arXiv

On the Physics of Radio Halos in Galaxy Clusters: Scaling Relations and Luminosity Functions

The underlying physics of giant and mini radio halos in galaxy clusters is still an open question. We find that mini halos (such as in Perseus and Ophiuchus) can be explained by radio-emitting electrons that are generated in hadronic cosmic ray (CR) interactions with protons of the intracluster medium. By contrast, the hadronic model either fails to explain the extended emission of giant radio halos (as in Coma at low frequencies) or would require a flat CR profile, which can be realized through outward streaming and diffusion of CRs (in Coma and A2163 at 1.4 GHz). We suggest that a second, leptonic component could be responsible for the missing flux in the outer parts of giant halos within a new hybrid scenario and we describe its possible observational consequences. To study the hadronic emission component of the radio halo population statistically, we use a cosmological mock galaxy cluster catalog built from the MultiDark simulation. Because of the properties of CR streaming and the different scalings of the X-ray luminosity (L_X) and the Sunyaev-Zel'dovich flux (Y) with gas density, our model can simultaneously reproduce the observed bimodality of radio-loud and radio-quiet clusters at the same L_X as well as the unimodal distribution of radio-halo luminosity versus Y; thereby suggesting a physical solution to this apparent contradiction. We predict radio halo emission down to the mass scale of galaxy groups, which highlights the unique prospects for low-frequency radio surveys (such as the LOFAR Tier 1 survey) to increase the number of detected radio halos by at least an order of magnitude.

preprint2015arXiv

SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 12: galaxy target selection and large scale structure catalogues

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2015arXiv

The SDSS-IV extended Baryonic Oscillation Spectroscopic Survey: Luminous Red Galaxy Target Selection

We describe the algorithm used to select the Luminous Red Galaxy (LRG) sample for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV) using photometric data from both the SDSS and the Wide-Field Infrared Survey Explorer (WISE). LRG targets are required to meet a set of color selection criteria and have z-band and i-band MODEL magnitudes z < 19.95 and 19.9 < i < 21.8, respectively. Our algorithm selects roughly 50 LRG targets per square degree, the great majority of which lie in the redshift range 0.6 < z < 1.0 (median redshift 0.71). We demonstrate that our methods are highly effective at eliminating stellar contamination and lower-redshift galaxies. We perform a number of tests using spectroscopic data from SDSS-III/BOSS to determine the redshift reliability of our target selection and its ability to meet the science requirements of eBOSS. The SDSS spectra are of high enough signal-to-noise ratio that at least 89% of the target sample yields secure redshift measurements. We also present tests of the uniformity and homogeneity of the sample, demonstrating that it should be clean enough for studies of the large-scale structure of the universe at higher redshifts than SDSS-III/BOSS LRGs reached.

preprint2014arXiv

Collision-free motion planning for fiber positioner robots: discretization of velocity profiles

The next generation of large-scale spectroscopic survey experiments such as DESI, will use thousands of fiber positioner robots packed on a focal plate. In order to maximize the observing time with this robotic system we need to move in parallel the fiber-ends of all positioners from the previous to the next target coordinates. Direct trajectories are not feasible due to collision risks that could undeniably damage the robots and impact the survey operation and performance. We have previously developed a motion planning method based on a novel decentralized navigation function for collision-free coordination of fiber positioners. The navigation function takes into account the configuration of positioners as well as their envelope constraints. The motion planning scheme has linear complexity and short motion duration (~2.5 seconds with the maximum speed of 30 rpm for the positioner), which is independent of the number of positioners. These two key advantages of the decentralization designate the method as a promising solution for the collision-free motion-planning problem in the next-generation of fiber-fed spectrographs. In a framework where a centralized computer communicates with the positioner robots, communication overhead can be reduced significantly by using velocity profiles consisting of a few bits only. We present here the discretization of velocity profiles to ensure the feasibility of a real-time coordination for a large number of positioners. The modified motion planning method that generates piecewise linearized position profiles guarantees collision-free trajectories for all the robots. The velocity profiles fit few bits at the expense of higher computational costs.

preprint2014arXiv

Constraining the halo bispectrum in real and redshift space from perturbation theory and non-linear stochastic bias

We present a method to produce mock galaxy catalogues with efficient perturbation theory schemes, which match the number density, power spectra and bispectra in real and in redshift space from N-body simulations. The essential contribution of this work is the way in which we constrain the bias parameters in the PATCHY-code. In addition of aiming at reproducing the two-point statistics, we seek the set of bias parameters, which constrain the univariate halo probability distribution function (PDF) encoding higher-order correlation functions. We demonstrate that halo catalogues based on the same underlying dark matter field with a fix halo number density, and accurately matching the power spectrum (within 2%), can lead to very different bispectra depending on the adopted halo bias model. A model ignoring the shape of the halo PDF can lead to deviations up to factors of 2. The catalogues obtained additionally constraining the shape of the halo PDF can significantly lower the discrepancy in the three-point statistics, yielding closely unbiased bispectra both in real and in redshift space; which are in general compatible with those corresponding to an N-body simulation within 10% (deviating at most up to 20%). Our calculations show that the constant linear bias of ~2 for Luminous Red Galaxy (LRG) like galaxies seen in the power spectrum, mainly comes from sampling halos in high density peaks, choosing a high density threshold rather than from a factor multiplying the dark matter density field. Our method contributes towards an efficient modelling of the halo/galaxy distribution required to estimate uncertainties in the clustering measurements from galaxy redshift surveys. We have also demonstrated that it represents a powerful tool to test various bias models.

preprint2014arXiv

Developing Micro DC-Brushless Motor Driver and Position Control for Fiber Positioners

In the large-scale, Dark Energy Spectroscopic Instrument (DESI), thousands of fiber positioners will be used. Those are robotic positioners, with two axis, and having the size of a pen. They are tightly packed on the focal plane of the telescope. Dedicated micro-robots have been developed and they use 4mm brushless DC motors. To simplify the implementation and reduce the space occupancy, each actuator will integrate its own electronic control board. This board will be used to communicate with the central trajectory generator, manage low level control tasks and motor current feeding. In this context, we present a solution for a highly compact electronic. This electronic is composed of two layers. The first is the power stage that can drive simultaneously two brushless motors. The second one consists of a fast microcontroller and deals with different control tasks: communication, acquisition of the hall sensor signals, commutation of the motors phases, and performing position and current regulation. A set of diagnostic functions are also implemented to detect failure in the motors or the sensors, and to sense abnormal load change that may be the result of two robots colliding.

preprint2014arXiv

Hunting down systematics in baryon acoustic oscillations after cosmic high noon

Future dark energy experiments will require better and more accurate theoretical predictions for the baryonic acoustic oscillations (BAO) signature in the spectrum of cosmological perturbations. Here, we use large N-body simulations of the \LambdaCDM Planck cosmology to study any possible systematic shifts and damping in BAO due to the impact of nonlinear gravitational growth of structure, scale dependent and non-local bias, and redshift-space distortions. The effect of cosmic variance is largely reduced by dividing the tracer power spectrum by that from a BAO-free simulation starting with the same phases. This permits us to study with unprecedented accuracy (better than 0.02% for dark matter and 0.07% for low-bias halos) small shifts of the pristine BAO wavenumbers towards larger k, and non-linear damping of BAO wiggles in the power spectrum of dark matter and halo populations in the redshift range z=0-1. For dark matter, we provide an accurate parametrization of the evolution of αas a function of the linear growth factor D(z). For halo samples, with bias ranging from 1.2 to 2.8, we measure a typical BAO shift of ~0.25%, observed in real-space, which does not show an appreciable evolution with redshift within the uncertainties. Moreover, we report a constant shift as a function of halo bias. We find a different evolution of the damping of the acoustic feature in all halo samples as compared to dark matter with haloes suffering less damping, and also find some weak dependence on bias. A larger BAO shift and damping is measured in redshift-space which can be well explained by linear theory due to redshift-space distortions. A clear modulation in phase with the acoustic scale is observed in the scale-dependent halo bias due to the presence of the baryonic acoustic oscillations.

preprint2014arXiv

Tests of redshift-space distortions models in configuration space for the analysis of the BOSS final data release

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. In preparation for analysis of redshift-space distortions from the Baryon Oscillation Spectroscopic Survey (BOSS) final data release we compare a number of analytic and phenomenological `streaming' models, specified in configuration space, to mock catalogs derived in different ways from several N-body simulations. The galaxies in each mock catalog have properties similar to those of the higher redshift galaxies measured by BOSS but differ in the details of how small-scale velocities and halo occupancy are determined. We find that all of the analytic models fit the simulations over a limited range of scales while failing at small scales. We discuss which models are most robust and on which scales they return reliable estimates of the rate of growth of structure: we find that models based on some form of resummation can fit our N-body data for BOSS-like galaxies above $30\,h^{-1}$Mpc well enough to return unbiased parameter estimates.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $Λ$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,σ$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

preprint2014arXiv

The flattening of the concentration-mass relation towards low halo masses and its implications for the annihilation signal boost

In the standard Cold Dark Matter (CDM) theory for understanding the formation of structure in the universe, there exists a tight connection between the properties of dark matter (DM) haloes, and their formation epochs. Such relation can be expressed in terms of a single key parameter, namely the halo concentration. In this work, we examine the median concentration-mass relation, c(M), at present time, over more than 20 orders of magnitude in halo mass, i.e., from tiny Earth-mass microhalos up to galaxy clusters. The c(M) model proposed by Prada et al. (2012), which links the halo concentration with the r.m.s. amplitude of matter linear fluctuations, describes remarkably well all the available N-body simulation data down to ~10^(-6) Msun microhalos. A clear fattening of the halo concentration-mass relation towards smaller masses is observed, that excludes the commonly adopted power-law c(M) models, and stands as a natural prediction for the CDM paradigm. We provide a parametrization for the c(M) relation that works accurately for all halo masses. This feature in the c(M) relation at low masses has decisive consequences e.g. for gamma-ray DM searches, as it implies more modest boosts of the DM annihilation flux due to substructure, i.e., ~35 for galaxy clusters and ~15 for galaxies like our own, as compared to those huge values adopted in the literature that rely on such power-law c(M) extrapolations. We provide a parametrization of the boosts that can be safely used for dwarfs to galaxy cluster-size halos.

preprint2014arXiv

The SDSS-III Baryonic Oscillation Spectroscopic Survey: Constraints on the Integrated Sachs Wolfe effect

In the context of the study of the Integrated Sachs Wolfe effect (ISW), we construct a template of the projected density distribution up to $z\simeq 0.7$ by using the Luminous Galaxies (LGs) from the Sloan Digital Sky Survey DR8. We use a photo-z catalogue trained with more than a hundred thousand galaxies from BOSS in the SDSS DR8 imaging area. We consider two different LG samples whose selection matches that of SDSS-III/BOSS: the LOWZ sample ($z\in [0.15,0.5]$) and the CMASS sample ($z\in[0.4,0.7]$). When building the LG density maps we use the information from star density, survey footprint, seeing conditions, sky emission, dust extinction and airmass to explore the impact of these artifacts on the two LG samples. In agreement with previous studies, we find that the CMASS sample is particularly sensitive to Galactic stars, which dominate the contribution to the auto-angular power spectrum below $\ell=7$. Other potential systematics affect mostly the low multipole range ($\ell\in[2,7]$), but leave fluctuations on smaller scales practically unchanged. The resulting power spectra in the multipole range $\ell\in[2,100]$ for the LOWZ, CMASS and LOWZ+CMASS samples are compatible with linear $Λ$CDM expectations and constant bias values of $b=1.98 \pm 0.11$, $2.08\pm0.14$ and $1.88\pm 0.11$, respectively, with no traces of non-Gaussianity: $f_{\rm NL}^{\rm local}=59\pm 75$ at 95% confidence level for the full LOWZ+CMASS sample in the range $\ell\in[4,100]$. After cross-correlating WMAP-9yr data with the LOWZ+CMASS LG density field, the ISW signal is detected at the level of 1.62--1.69$\,σ$. While this result is in close agreement with predictions from Monte Carlo simulations in the concordance $Λ$CDM model, it cannot rule out by itself an Einstein-de Sitter scenario, and has a moderately low signal compared to previous studies conducted on subsets of this LG sample.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2013arXiv

Collision Avoidance in Next-generation Fiber Positioner Robotic System for Large Survey Spectrograph

Some of the next generation massive spectroscopic survey projects, such as DESI and PFS, plan to use thousands of fiber positioner robots packed at a focal plane to quickly move in parallel the fiber-ends from the previous to the next target points. The most direct trajectories are prone to collision that could damage the robots and impact the survey operation. We thus present here a motion planning method based on a novel decentralized navigation function for collision-free coordination of fiber positioners. The navigation function takes into account the configuration of positioners as well as the actuator constraints. We provide details for the proof of convergence and collision avoidance. Decentralization results in linear complexity for the motion planning as well as dependency of motion duration with respect to the number of positioners. Therefore the coordination method is scalable for large-scale spectrograph robots. The short in-motion duration of positioner robots (~2.5 seconds using typical actuator constraints), will thus allow the time dedicated for observation to be maximized.

preprint2013arXiv

Constraints on WIMP Annihilation for Contracted Dark Matter in the Inner Galaxy with the Fermi-LAT

We derive constraints on parameters of generic dark matter candidates by comparing theoretical predictions with the gamma-ray emission observed by the Fermi-LAT from the region around the Galactic Center. Our analysis is conservative since it simply requires that the expected dark matter signal does not exceed the observed emission. The constraints obtained in the likely case that the collapse of baryons to the Galactic Center is accompanied by the contraction of the dark matter are strong. In particular, we find that for bb and τ+τ- or W+W- dark matter annihilation channels, the upper limits on the annihilation cross section imply that the thermal cross section is excluded for a Weakly Interacting Massive Particle (WIMP) mass smaller than about 700 and 500 GeV, respectively. For the μ+ μ- channel, where the effect of the inverse Compton scattering is important, depending on models of the Galactic magnetic field the exclusion of the thermal cross-section is for a WIMP mass smaller than about 150 to 400 GeV. The upper limits on the annihilation cross section of dark matter particles obtained are two orders of magnitude stronger than without contraction. In the latter case our results are compatible with the upper limits from the Galactic halo analysis reported by the Fermi-LAT collaboration for the case in which the same conservative approach without modeling of the astrophysical background is employed.

preprint2013arXiv

Halo Abundance Matching: accuracy and conditions for numerical convergence

Accurate predictions of the abundance and clustering of dark matter haloes play a key role in testing the standard cosmological model. Here, we investigate the accuracy of one of the leading methods of connecting the simulated dark matter haloes with observed galaxies -- the Halo Abundance Matching (HAM) technique. We show how to choose the optimal values of the mass and force resolution in large-volume N-body simulations so that they provide accurate estimates for correlation functions and circular velocities for haloes and their subhaloes -- crucial ingredients of the HAM method. At the 10% accuracy, results converge for 50 particles for haloes and 150 particles for progenitors subhaloes. In order to achieve this level of accuracy a number of conditions should be satisfied. The force resolution for the smallest resolved (sub)haloes should be in the range (0.1-0.3)rs, where rs is the scale radius of (sub)haloes. The number of particles for progenitors of subhaloes should be 150. We also demonstrate that the two-body scattering plays a minor role for the accuracy of N-body simulations thanks to the relatively small number of crossing-times of dark matter in haloes, and the limited force resolution of cosmological simulations.

preprint2013arXiv

Modelling Baryon Acoustic Oscillations with Perturbation Theory and Stochastic Halo Biasing

In this work we investigate the generation of mock halo catalogues based on perturbation theory and nonlinear stochastic biasing with the novel PATCHY-code. In particular, we use Augmented Lagrangian Perturbation Theory (ALPT) to generate a dark matter density field on a mesh starting from Gaussian fluctuations and to compute the peculiar velocity field. ALPT is based on a combination of second order LPT (2LPT) on large scales and the spherical collapse model on smaller scales. We account for the systematic deviation of perturbative approaches from N-body simulations together with halo biasing adopting an exponential bias model. We then account for stochastic biasing by defining three regimes: a low, an intermediate and a high density regime, using a Poisson distribution in the intermediate regime and the negative binomial distribution to model over-dispersion in the high density regime. Since we focus in this study on massive halos, we suppress the generation of halos in the low density regime. The various nonlinear and stochastic biasing parameters, and density thresholds (five) are calibrated with the large BigMultiDark N-body simulation to match the power spectrum of the corresponding halo population. Our mock catalogues show power spectra, both in real- and redshift-space, which are compatible with N-body simulations within about 2% up to k ~ 1 h Mpc^-1 at z = 0.577 for a sample of halos with the typical BOSS CMASS galaxy number density. The corresponding correlation functions are compatible down to a few Mpc. We also find that neglecting over-dispersion in high density regions produces power spectra with deviations of 10% at k ~ 0.4 h Mpc^-1. These results indicate the need to account for an accurate statistical description of the galaxy clustering for precise studies of large-scale surveys.

preprint2013arXiv

On the Physics of Radio Halos in Galaxy Clusters: Scaling Relations and Luminosity Functions

The underlying physics of giant radio halos and mini halos in galaxy clusters is still an open question, which becomes more pressing with the growing number of detections. In this paper, we explore the possibility that radio-emitting electrons are generated in hadronic cosmic ray (CR) proton interactions with ambient thermal protons of the intra-cluster medium. Our CR model derives from cosmological hydrodynamical simulations of cluster formation and additionally accounts for CR transport in the form of CR streaming and diffusion. This opens the possibility of changing the radio halo luminosity by more than an order of magnitude on a dynamical time scale. We build a mock galaxy cluster catalog from the large MultiDark N-body LCDM simulation by adopting a phenomenological gas density model for each cluster based on X-ray measurements that matches Sunyaev-Zel'dovich (SZ) and X-ray scaling relations and luminosity function. Using magnetic field strength estimates from Faraday rotation measure studies, our model successfully reproduces the observed surface brightness profiles of giant radio halos (Coma, A2163) as well as radio mini-halos (Perseus, Ophiuchus), while obeying upper limits on the gamma-ray emission in these clusters. Our model is also able to simultaneously reproduce the observed bimodality of radio-loud and radio-quiet clusters at the same L_X as well as the unimodal distribution of radio-halo luminosity versus the SZ flux Y; thereby suggesting a physical solution to this apparent contradiction. For a plausible fraction of 10% radio-loud clusters, our model matches the NVSS radio-halo luminosity function. Constructing an analytical radio-halo luminosity function, we demonstrate the unique prospects for low-frequency radio surveys (such as the LOFAR Tier 1 survey) to detect ~3500 radio halos back to redshift two and to probe the underlying physics of radio halos. [abridged]

preprint2013arXiv

Stochastic bias of colour-selected BAO tracers by joint clustering-weak lensing analysis

The baryon acoustic oscillation (BAO) feature in the two-point correlation function of galaxies supplies a standard ruler to probe the expansion history of the Universe. We study here several galaxy selection schemes, aiming at building an emission-line galaxy (ELG) sample in the redshift range $0.6<z<1.7$, that would be suitable for future BAO studies, providing a highly biased galaxy sample. We analyse the angular galaxy clustering of galaxy selections at the redshifts 0.5, 0.7, 0.8, 1 and 1.2 and we combine this analysis with a halo occupation distribution (HOD) model to derive the properties of the haloes these galaxies inhabit, in particular the galaxy bias on large scales. We also perform a weak lensing analysis (aperture statistics) to extract the galaxy bias and the cross-correlation coefficient and compare to the HOD prediction. We apply this analysis on a data set composed of the photometry of the deep co-addition on Sloan Digital Sky Survey (SDSS) Stripe 82 (225 deg$^2$), of Canda-France-Hawai Telescope/Stripe 82 deep \emph{i}-band weak lensing survey and of the {\it Wide-Field Infrared Survey Explorer}infrared photometric band W1. The analysis on the SDSS-III/constant mass galaxies selection at $z=0.5$ is in agreement with previous studies on the tracer, moreover we measure its cross-correlation coefficient $r=1.16\pm0.35$. For the higher redshift bins, we confirm the trends that the brightest galaxy populations selected are strongly biased ($b>1.5$), but we are limited by current data sets depth to derive precise values of the galaxy bias. A survey using such tracers of the mass field will guarantee a high significance detection of the BAO.

preprint2013arXiv

The clustering of galaxies at z~0.5 in the SDSS-III Data Release 9 BOSS-CMASS sample: a test for the LCDM cosmology

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological constraints from the full shape of the clustering wedges

We explore the cosmological implications of the clustering wedges, xi_perp(s) and xi_para(s), of the CMASS Data Release 9 (DR9) sample of the Baryon Oscillation Spectroscopic Survey (BOSS). These clustering wedges are defined by averaging the full two-dimensional correlation function, xi(mu,s), over the ranges 0<mu<0.5 and 0.5<mu<1, respectively. These measurements allow us to constrain the parameter combinations D_A(z)/r_s(z_d)=9.03 +- 0.21 and cz/(r_s(z_d)H(z)) = 12.14 +- 0.43 at the mean redsfhit of the sample, z=0.57. We combine the information from the clustering wedges with recent measurements of CMB, BAO and type Ia supernovae to obtain constraints on the cosmological parameters of the standard LCDM model and a number of potential extensions. The information encoded in the clustering wedges is most useful when the dark energy equation of state is allowed to deviate from its standard LCDM value. The combination of all datasets shows no evidence of a deviation from a constant dark energy equation of state, in which case we find w_DE = -1.013 +- 0.064, in complete agreement with a cosmological constant. We explore potential deviations from general relativity by constraining the growth rate f(z)=d ln D(a)/ d ln a, in which case the combination of the CMASS clustering wedges with CMB data implies f(z=0.57)=0.719 +- 0.094, in accordance with the predictions of GR. Our results clearly illustrate the additional constraining power of anisotropic clustering measurements with respect to that of angle-averaged quantities.

preprint2013arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Measuring D_A and H at z=0.57 from the Baryon Acoustic Peak in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of the angular diameter distance to and Hubble parameter at z=0.57 from the measurement of the baryon acoustic peak in the correlation of galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey. Our analysis is based on a sample from Data Release 9 of 264,283 galaxies over 3275 square degrees in the redshift range 0.43<z<0.70. We use two different methods to provide robust measurement of the acoustic peak position across and along the line of sight in order to measure the cosmological distance scale. We find D_A(0.57) = 1408 +/- 45 Mpc and H(0.57) = 92.9 +/- 7.8 km/s/Mpc for our fiducial value of the sound horizon. These results from the anisotropic fitting are fully consistent with the analysis of the spherically averaged acoustic peak position presented in Anderson et al, 2012. Our distance measurements are a close match to the predictions of the standard cosmological model featuring a cosmological constant and zero spatial curvature.

preprint2013arXiv

The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Measuring H(z) and D_A(z) at z = 0.57 with Clustering Wedges

We analyze the 2D correlation function of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to measure cosmic expansion H and the angular diameter distance D_A at a mean redshift of <z> = 0.57. We apply, for the first time, a new correlation function technique called clustering wedges. Using a physically motivated model, the anisotropic baryonic acoustic feature in the galaxy sample is detected at a significance level of 4.7 sigma compared to a featureless model. The baryonic acoustic feature is used to obtain model independent constraints cz/H/r_s = 12.28 +- 0.82 (6.7 per-cent accuracy) and D_A/r_s = 9.05 +- 0.27 (3.0 per-cent) with a correlation coefficient of -0.5, where r_s is the sound horizon scale at the end of the baryonic drag era. We conduct thorough tests on the data and 600 simulated realizations, finding robustness of the results regardless of the details of the analysis method. Combining with r_s constraints from the Cosmic Microwave Background we obtain H(0.57) = 90.8 +- 6.2 kms-1Mpc-1 and D_A(0.57) = 1386 +- 45 Mpc. We use simulations to forecast results of the final BOSS CMASS data set. We apply the reconstruction technique on the simulations demonstrating that the sharpening of the anisotropic baryonic acoustic feature should improve the detection as well as tighten constraints of H and D_A by 30 per-cent on average.

preprint2012arXiv

Acoustic scale from the angular power spectra of SDSS-III DR8 photometric luminous galaxies

We measure the acoustic scale from the angular power spectra of the Sloan Digital Sky Survey III (SDSS-III) Data Release 8 imaging catalog that includes 872,921 galaxies over ~ 10,000 deg^2 between 0.45<z<0.65. The extensive spectroscopic training set of the Baryon Oscillation Spectroscopic Survey (BOSS) luminous galaxies allows precise estimates of the true redshift distributions of galaxies in our imaging catalog. Utilizing the redshift distribution information, we build templates and fit to the power spectra of the data, which are measured in our companion paper, Ho et al. 2011, to derive the location of Baryon acoustic oscillations (BAO) while marginalizing over many free parameters to exclude nearly all of the non-BAO signal. We derive the ratio of the angular diameter distance to the sound horizon scale D_A/r_s= 9.212 + 0.416 -0.404 at z=0.54, and therefore, D_A= 1411+- 65 Mpc at z=0.54; the result is fairly independent of assumptions on the underlying cosmology. Our measurement of angular diameter distance D_A is 1.4 σhigher than what is expected for the concordance LCDM (Komatsu et al. 2011), in accordance to the trend of other spectroscopic BAO measurements for z >~ 0.35. We report constraints on cosmological parameters from our measurement in combination with the WMAP7 data and the previous spectroscopic BAO measurements of SDSS (Percival et al. 2010) and WiggleZ (Blake et al. 2011). We refer to our companion papers (Ho et al. 2011; de Putter et al. 2011) for investigations on information of the full power spectrum.

preprint2012arXiv

Characterization of Dark-Matter-induced anisotropies in the diffuse gamma-ray background

The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing template all-sky gamma-ray maps for the emission produced in the galactic halo and its substructures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art N-body simulations of cosmic structure formation, namely Millennium-II for extragalactic (sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of extrapolation to account for (sub)structures that are below the resolution limit of the simulations, allowing us to estimate the total emission all the way down to the minimal self-bound halo mass. We describe in detail the features appearing in the APS of our template maps and we estimate the effect of various uncertainties such as the value of the minimal halo mass, the fraction of substructures hosted in a halo and the shape of the DM density profile. Our results indicate that the fluctuation APS of the DM-induced emission is of the same order as the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the comparison with the measured anisotropy. We also quantify the uncertainties affecting our results, finding "theoretical error bands" spanning more than two orders of magnitude and dominated (for a given particle physics model) by our lack of knowledge of the abundance of low-mass (sub)halos.

preprint2012arXiv

Clustering of Sloan Digital Sky Survey III Photometric Luminous Galaxies: The Measurement, Systematics and Cosmological Implications

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find Ω_Λ= 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $Ω_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

preprint2012arXiv

Quantifying the coherent outflows of galaxies around voids in the SDSS DR7

We report the detection, with a high level of confidence, of coherent outflows around voids found in the seventh data release of the Sloan Digital Sky Survey (SDSS DR7). In particular, we developed a robust <|cos theta|> statistical test to quantify the strength of redshift-space distortions (RSD) associated with extended coherent velocity fields. We consistently find that the vector that joins void centers with galaxies that lie in shells around them is more likely to be perpendicular to the line-of-sight than parallel to it. This effect is clear evidence for the existence of outflows in the vicinity of voids. We show that the RSD exist for a wide range of void radius and shell thickness, but they are more evident in the largest voids in our sample. For instance, we find that the <|cos theta|> for galaxies located in shells within 2 h^-1 Mpc from the edge of voids larger than 15 h^-1 Mpc deviates 3.81sigma from uniformity. The measurements presented in this work provide useful information to constrain cosmological parameters, in particular Omega_m and Sigma_8.

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

preprint2012arXiv

The Clustering of Galaxies in SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Constraints on Primordial Non-Gaussianity

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Analysis of potential systematics

We analyze the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of 3,275 deg^2. We identify, characterize, and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the `CMASS' sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120h^-1Mpc or k < 0.01hMpc^-1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. We use mock galaxy catalogs that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on correlation function measurements and that this systematic error is negligible for the spherically averaged correlation function. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature (Anderson et al. 2012), constrain cosmological models using the full shape of the correlation function (Sanchez et al. 2012), and measure the rate of structure growth (Reid et al. 2012). (abridged)

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: the low redshift sample

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

Towards a panchromatic picture of galaxy evolution during the reionization epoch

There are thousands of confirmed detections of star forming galaxies at high redshift (z > 4). These observations rely primarily on the detection of the spectral Lyman Break and the Lyman-alpha emission line. Theoretical modelling of these sources helps to interpret the observations in the framework of the standard cosmological paradigm. We present results from the High-z MareNostrum Project, aimed at constructing a panchromatic picture of the high redshift galaxy evolution that will improve our understanding of young star forming galaxies. Our simulation successfully reproduces the observational constraints from Lyman Break Galaxies and Lyman-alpha emitters at 5 < z < 7 . Based on this model we make predictions on the expected Far Infrared (FIR) emission that should be observed for LAEs. These predictions will help to settle down the question on the dust content of massive high-z galaxies, an issue that will be feasible to probe observationally with the Atacama Large Millimetre Array (ALMA).

preprint2011arXiv

Ameliorating Systematic Uncertainties in the Angular Clustering of Galaxies: A Study using SDSS-III

We investigate the effects of potential sources of systematic error on the angular and photometric redshift, z_phot, distributions of a sample of redshift 0.4 < z < 0.7 massive galaxies whose selection matches that of the Baryon Oscillation Spectroscopic Survey (BOSS) constant mass sample. Utilizing over 112,778 BOSS spectra as a training sample, we produce a photometric redshift catalog for the galaxies in the SDSS DR8 imaging area that, after masking, covers nearly one quarter of the sky (9,913 square degrees). We investigate fluctuations in the number density of objects in this sample as a function of Galactic extinction, seeing, stellar density, sky background, airmass, photometric offset, and North/South Galactic hemisphere. We find that the presence of stars of comparable magnitudes to our galaxies (which are not traditionally masked) effectively remove area. Failing to correct for such stars can produce systematic errors on the measured angular auto-correlation function, w, that are larger than its statistical uncertainty. We describe how one can effectively mask for the presence of the stars, without removing any galaxies from the sample, and minimize the systematic error. Additionally, we apply two separate methods that can be used to correct the systematic errors imparted by any parameter that can be turned into a map on the sky. We find that failing to properly account for varying sky background introduces a systematic error on w. We measure w, in four z_phot slices of width 0.05 between 0.45 < z_phot < 0.65 and find that the measurements, after correcting for the systematic effects of stars and sky background, are generally consistent with a generic LambdaCDM model, at scales up to 60 degrees. At scales greater than 3 degrees and z_phot > 0.5, the magnitude of the corrections we apply are greater than the statistical uncertainty in w.

preprint2011arXiv

Dark halo response and the stellar initial mass function in early-type and late-type galaxies

We investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

preprint2011arXiv

Dark Matter implications of the Fermi-LAT measurement of anisotropies in the diffuse gamma-ray background: status report

For the first time, the Fermi-LAT measured the angular power spectrum (APS) of anisotropies in the diffuse gamma-ray background. The data is found to be broadly compatible with a model with contributions from the point sources in the 1-year catalog, the galactic diffuse background, and the extragalactic isotropic emission; however deviations are present at both large and small angular scales. In this study, we complement the model with a contribution from Dark Matter (DM) whose distribution is modeled exploiting the results of the most recent N-body simulations, considering the contribution of extragalactic halos and subhalos (from Millennium-II) and of galactic substructures (from Aquarius). With the use of the Fermi Science Tools, these simulations serve as templates to produce mock gamma-ray count maps for DM gamma-ray emission, both in the case of an annihilating and a decaying DM candidate. The APS will then be computed and compared with the Fermi-LAT results to derive constraints on the DM particle physics properties. The possible systematic due to an imperfect model of the galactic foreground is also studied and taken into account properly. The present paper reports on the status of the project.

preprint2011arXiv

Fiber Assignment in Next-generation Wide-field Spectrographs

We present an optimized algorithm for assigning fibers to targets in next-generation fiber-fed multi-object spectrographs. The method, that we named draining algorithm, ensures that the maximum number of targets in a given target field is observed in the first few tiles. Using randomly distributed targets and mock galaxy catalogs we have estimated that the gain provided by the draining algorithm as compared to a random assignment can be as much as 2% for the first tiles. This would imply for a survey like BigBOSS saving for observation several hundred thousand objects or, alternatively, reducing the covered area in ~350 sq. deg. An important advantage of this method is that the fiber collision problem can be solved easily and in an optimal way. We also discuss additional optimizations of the fiber positioning process. In particular, we show that allowing for rotation of the focal plane can improve the efficiency of the process in ~3.5-4.5% even if only small adjustments are permitted (up to 2 deg). For instruments that allow large rotations of the focal plane the expected gain increases to ~5-6%. These results, therefore, strongly support focal plane rotation in future spectrographs, as far as the efficiency of the fiber positioning process is concerned. Finally, we discuss on the implications of our optimizations and provide some basic hints for an optimal survey strategy based on the number of targets per positioner.

preprint2011arXiv

Halo concentrations in the standard LCDM cosmology

We study the concentration of dark matter halos and its evolution in N-body simulations of the standard LCDM cosmology. The results presented in this paper are based on 4 large N-body simulations with about 10 billion particles each: the Millennium-I and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simulation is introduced in this paper. This suite of simulations with high mass resolution over a large volume allows us to compute with unprecedented accuracy the concentration over a large range of scales (about six orders of magnitude in mass), which constitutes the state-of-the-art of our current knowledge on this basic property of dark matter halos in the LCDM cosmology. We find that there is consistency among the different simulation data sets. We confirm a novel feature for halo concentrations at high redshifts: a flattening and upturn with increasing mass. The concentration c(M,z) as a function of mass and the redshift and for different cosmological parameters shows a remarkably complex pattern. However, when expressed in terms of the linear rms fluctuation of the density field sigma(M,z), the halo concentration c(sigma) shows a nearly-universal simple U-shaped behaviour with a minimum at a well defined scale at sigma=0.71. Yet, some small dependences with redshift and cosmology still remain. At the high-mass end (sigma < 1) the median halo kinematic profiles show large signatures of infall and highly radial orbits. This c-sigma(M,z) relation can be accurately parametrized and provides an analytical model for the dependence of concentration on halo mass. When applied to galaxy clusters, our estimates of concentrations are substantially larger -- by a factor up to 1.5 -- than previous results from smaller simulations, and are in much better agreement with results of observations. (abridged)

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The MultiDark Database: Release of the Bolshoi and MultiDark Cosmological Simulations

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

preprint2010arXiv

A Fiber Positioner Robot for the Gran Telescopio Canarias

Fiber-fed spectrographs dedicated to observing massive portions of the sky are increasingly being more demanded within the astronomical community. For all the fiber-fed instruments, the primordial and common problem is the positioning of the fiber ends, which must match the position of the objects of a target field on the sky. Amongst the different approaches found in the state of the art, actuator arrays are one of the best. Indeed, an actuator array is able to position all the fiber heads simultaneously, thus making the reconfiguration time extremely short and the instrument efficiency high. The SIDE group (see http://side.iaa.es) at the Instituto de Astrofísica de Andalucía, together with the industrial company AVS and the University of Barcelona, has been developing an actuator suitable for a large and scalable array. A real-scale prototype has been built and tested in order to validate its innovative design concept, as well as to verify the fulfillment of the mechanical requirements. The present article describes both the concept design and the test procedures and conditions. The main results are shown and a full justification of the validity of the proposed concept is provided.

preprint2010arXiv

The Kinematic Connection Between Galaxies and Dark Matter Haloes

Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt \ne V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.

preprint2010arXiv

The search of axion-like-particles with Fermi and Cherenkov telescopes

Axion Like Particles (ALPs), postulated to solve the strong-CP problem, are predicted to couple with photons in the presence of magnetic fields, which may lead to a significant change in the observed spectra of gamma-ray sources such as AGNs. Here we simultaneously consider in the same framework both the photon/axion mixing that takes place in the gamma-ray source and that one expected to occur in the intergalactic magnetic fields. We show that photon/axion mixing could explain recent puzzles regarding the observed spectra of distant gamma-ray sources as well as the recently published lower limit to the EBL intensity. We finally summarize the different signatures expected and discuss the best strategy to search for ALPs with the Fermi satellite and current Cherenkov telescopes like CANGAROO, HESS, MAGIC and VERITAS.

preprint2009arXiv

Modeling gamma-ray burst observations by Fermi and MAGIC including attenuation due to diffuse background light

Gamma rays from extragalactic sources are attenuated by pair-production interactions with diffuse photons of the extragalactic background light (EBL). Gamma-ray bursts (GRBs) are a source of high-redshift photons above 10 GeV, and could be therefore useful as a probe of the evolving UV background radiation. In this paper, we develop a simple phenomenological model for the number and redshift distribution of gamma-ray bursts that can be seen at GeV energies with the Fermi satellite and MAGIC atmospheric Cherenkov telescope. We estimate the observed number of gamma rays per year, and show how this result is modified by considering interactions with different realizations of the evolving EBL. We also discuss the bright Fermi GRB 080916C in the context of this model. We find that the LAT on Fermi can be expected to see a small number of photons above 10 GeV each year from distant GRBs. Annual results for ground-based instruments like MAGIC are highly variable due to the low duty cycle and sky coverage of the telescope. However, successfully viewing a bright or intermediate GRB from the ground could provide hundreds of photons from high redshift, which would almost certainly be extremely useful in constraining both GRB physics and the high-redshift EBL.

preprint2009arXiv

The orientation and kinematics of inner tidal tails around dwarf galaxies orbiting the Milky Way

Using high-resolution collisionless N-body simulations we study the properties of tidal tails formed in the immediate vicinity of a two-component dwarf galaxy evolving in a static potential of the Milky Way (MW). The stellar component of the dwarf is initially in the form of a disk and the galaxy is placed on an eccentric orbit motivated by CDM-based cosmological simulations. We measure the orientation, density and velocity distribution of the stars in the tails. Due to the geometry of the orbit, in the vicinity of the dwarf, where the tails are densest and therefore most likely to be detectable, they are typically oriented towards the MW and not along the orbit. We report on an interesting phenomenon of `tidal tail flipping': on the way from the pericentre to the apocentre the old tails following the orbit are dissolved and new ones pointing towards the MW are formed over a short timescale. We also find a tight linear relation between the velocity of stars in the tidal tails and their distance from the dwarf. Using mock data sets we demonstrate that if dwarf spheroidal (dSph) galaxies in the vicinity of the MW are tidally affected their kinematic samples are very likely contaminated by tidally stripped stars which tend to artificially inflate the measured velocity dispersion. The effect is stronger for dwarfs on their way from the peri- to the apocentre due to the formation of new tidal tails after pericentre. Realistic mass estimates of dSph galaxies thus require removal of these stars from kinematic samples.

preprint2009arXiv

The origin of the diffuse non-thermal X-ray and radio emission in the Ophiuchus cluster of galaxies

We present high resolution 240 and 607 MHz GMRT radio observations, complemented with 74 MHz archival VLA radio observations of the Ophiuchus cluster of galaxies, whose radio mini-halo has been recently detected at 1400 MHz. We also present archival Chandra and XMM-Newton data of the Ophiuchus cluster. Our observations do not show significant radio emission from the mini-halo, hence we present upper limits to the integrated, diffuse non-thermal radio emission of the core of the Ophiuchus cluster. The XMM-Newton observations can be well explained by a two-temperature thermal model with temperatures of ~=1.8 keV and ~=9.0 keV, respectively, which confirms previous results that suggest that the innermost central region of the Ophiuchus cluster is a cooling core. We also used the XMM-Newton data to set up an upper limit to the (non-thermal) X-ray emission from the cluster. The combination of available radio and X-ray data has strong implications for the currently proposed models of the spectral energy distribution (SED) from the Ophiuchus cluster. In particular, a synchrotron+IC model is in agreement with the currently available data, if the average magnetic field is in the range (0.02-0.3) microG. A pure WIMP annihilation scenario can in principle reproduce both radio and X-ray emission, but at the expense of postulating very large boost factors from dark matter substructures, jointly with extremely low values of the average magnetic field. Finally, a scenario where synchrotron and inverse Compton emission arise from PeV electron-positron pairs (via interactions with the CMB), can be ruled out, as it predicts a non-thermal soft X-ray emission that largely exceeds the thermal Bremsstrahlung measured by INTEGRAL.

preprint2009arXiv

The statistics of voids as a tool to constrain cosmological parameters: sigma_8 and Omega_m h

We present a general analytical formalism to calculate accurately several statistics related to underdense regions in the Universe. The statistics are computed for dark matter halo and galaxy distributions both in real space and redshift space at any redshift. Using this formalism, we found that void statistics for galaxy distributions can be obtained, to a very good approximation, assuming galaxies to have the same clustering properties as halos above a certain mass. We deducted a relationship between this mass and that of halos with the same accumulated number density as the galaxies. We also found that the dependence of void statistics on redshift is small. For instance, the number of voids larger than 13 Mpc/h (defined to not contain galaxies brighter than M_r=-20.4 +5logh change less than 20% between z=1 and z=0. However, the dependence of void statistics on sigma_8 and Omega_m h is considerably larger, making them appropriate to develop tests to measure these parameters. We have shown how to efficiently construct several of these tests and discussed in detail the treatment of several observational effects. The formalism presented here along with the observed statistics extracted from current and future large galaxy redshift surveys will provide an independent measurement of the relevant cosmological parameters. Combining these measurements with those found using other methods will contribute to reduce their uncertainties.

preprint2007arXiv

Detailed theoretical predictions of the outskirts of dark matter halos

In the present work we describe the formalism necessary to derive the properties of dark matter halos beyond two virial radius using the spherical collapse model (without shell crossing), and provide the framework for the theoretical prediction presented in Prada et al. (2005). We show in detail how to obtain within this model the probability distribution for the spherically-averaged enclosed density at any radii P(delta,r). Using this probability distribution, we compute the most probable and mean density profiles, which turns out to differ considerably from each other. We also show how to obtain the typical profile, as well as the probability distribution and mean profile for the spherically averaged radial velocity. Two probability distributions are obtained: a first one is derived using a simple assumption, that is, if Q is the virial radius in Lagrangian coordinates, then the enclosed linear contrast delta_l(q,Q) must satisfy the condition that delta_l(q=Q) = delta_vir, where delta_vir is the linear density contrast within the virial radius Rvir at the moment of virialization. Then we introduce an additional constraint to obtain a more accurate P(delta,r) which reproduces to a higher degree of precision the distribution of the spherically averaged enclosed density found in the simulations. This new constraint is delta_l(q,Q) < delta_vir for all q > Q, which means that there are no radii larger than Rvir where the density contrast is larger than that used to define the virial radius. Finally, we compare in detail our theoretical predictions for the probability distributions with the results found in the simulations.

preprint2005arXiv

The ALHAMBRA Survey: For a systematic Study of Cosmic Evolution

ALHAMBRA is a project to gather data to sample a fraction of the Universe with enough precision to follow the evolution of its content and properties with z, a kind of Cosmic Tomography. It is defined as a 4 square degrees photometric survey with 20 contiguous, equal width, medium band filters covering from 3500 to 9700 A, plus the JHKs NIR bands. It was optimized to get (for a fixed observing time) the maximum number of objects with accurate classification and redshift and to be sensitive to relatively faint emission features. The observations will be carried out with the 3.5m telescope in Calar Alto (Spain). We intend to reach the limit AB=25 in all the optical filters from the bluest to 8300 A, and from 24.7 to 23.4 for the remainder. The expected limit in the NIR is fixed at Ks=20, H=21, J=22. The homogeneous and contiguous spectral coverage will result in several hundred thousand objects with accurate SED identification and z-values. The accuracy of the survey will allow to make detailed studies, including galactic, without the need for any further follow-up. Indeed, it will provide exciting targets for large telescopes, the GTC in particular.

preprint2003arXiv

Observing the dark matter density profile of isolated galaxies

Using the Sloan Digital Sky Survey (SDSS), we probe the halo mass distribution by studying the velocities of satellites orbiting isolated galaxies. In a subsample that covers 2500 sq. degrees on the sky, we detect about 3000 satellites with absolute blue magnitudes going down to M_B = -14; most of the satellites have M_B=-16 to -18, comparable to the magnitudes of M32 and the Magellanic Clouds. After a careful, model-independent removal of interlopers, we find that the line-of-sight velocity dispersion of satellites declines with distance to the primary. For an L* galaxy the r.m.s. line-of-sight velocity changes from ~120 km/s at 20 kpc to ~60 km/s at 350 kpc. This decline agrees remarkably well with theoretical expectations, as all modern cosmological models predict that the density of dark matter in the peripheral parts of galaxies declines as rho_DM propto r^{-3}. Thus, for the first time we find direct observational evidence of the density decline predicted by cosmological models; we also note that this result contradicts alternative theories of gravity such as MOND. We also find that the velocity dispersion of satellites within 100 kpc scales with the absolute magnitude of the central galaxy as sigma propto L^{0.3}; this is very close to the Tully-Fisher relation for normal spiral galaxies.

preprint2003arXiv

The near-IR properties and continuum shapes of high redshift quasars from the Sloan Digital Sky Survey

We present J-H-K' photometry for a sample of 45 high redshift quasars found by the Sloan Digital Sky Survey. The sample was originally selected on the basis of optical colors and spans a redshift range from 3.6 to 5.03. Our photometry reflects the rest-frame SED longward of Ly alpha for all redshifts. The results show that the near-IR colors of high redshift quasars are quite uniform. We have modelled the continuum shape of the quasars (from just beyond Ly alpha to ~4000 A) with a power law of the form f_nu \propto nu^alpha, and find <alpha > =-0.57 with a scatter of 0.33. This value is similar to what is found for lower redshift quasars over the same restframe wavelength range, and we conclude that there is hardly any evolution in the continuum properties of optically selected quasars up to redshift 5. The spectral indices found by combining near-IR with optical photometry are in general consistent but slightly flatter than what is found for the same quasars using the optical spectra and photometry alone, showing that the continuum region used to determine the spectral indices can somewhat influence the results.

preprint1995arXiv

New Parametrizations of Non-Gaussian Line-of-sight Velocity Distribution

A five-parameter fitting formula for the line-of-sight stellar velocity distributions of steady state systems is proposed. It can faithfully reproduce velocity distributions of theoretical models ranging from nearly Gaussian profiles to strongly skewed or mildly double-peaked profiles. In contrast to van der Marel and Franx (1993) and Kuijken and Merrifield (1993), the line profiles are required to have neither multi-peaks nor negative velocity wings. Information of the profile is mostly specified by five physically meaningful and nearly orthogonal fitting parameters.