Source author record

Enrique Gaztanaga

Enrique Gaztanaga 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

51works
8topics
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

51 published item(s)

preprint2024arXiv

The Gravitational Lensing Imprints of DES Y3 Superstructures on the CMB: A Matched Filtering Approach

$ $Low density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence $κ$. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 data set, covering approximately 4,200 deg$^2$ of the sky. We use two distinct void-finding algorithms: a 2D void-finder which operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the MICE N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure $A_κ$, the amplitude of the observed lensing signal relative to the simulation template, obtaining $A_κ= 1.03 \pm 0.22$ ($4.6σ$ significance) for Voxel and $A_κ= 1.02 \pm 0.17$ ($5.9σ$ significance) for 2D voids, both consistent with $Λ$CDM expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure $A_κ= 0.87 \pm 0.15$ ($5.9σ$ significance). The leading source of noise in our measurements is Planck noise, implying that future data from the Atacama Cosmology Telescope (ACT), South Pole Telescope (SPT) and CMB-S4 will increase sensitivity and allow for more precise measurements.

preprint2022arXiv

A peek outside our Universe

According to general relativity (GR), a universe with a cosmological constant $Λ$, like ours, is trapped inside an event horizon, $r< \sqrt{3/Λ}$. What is outside? We show, using Israel (1967) junction conditions, that there could be a different universe outside. Our universe looks like a black hole for an outside observer. Outgoing radial null geodesics cannot escape our universe, but incoming photons can enter and leave an imprint on our CMB sky. We present a picture of such a fossil record from the analysis of CMB maps that agrees with the black hole universe predictions but challenges our understanding of the origin of the primordial universe.

preprint2022arXiv

Clustering with general photo-$z$ uncertainties: Application to Baryon Acoustic Oscillations

Photometric data can be analyzed using the three-dimensional correlation function $ξ_{\rm p}$ to extract cosmological information via e.g., measurement of the Baryon Acoustic Oscillations (BAO). Previous studies modeled $ξ_{\rm p} $ assuming a Gaussian photo-$z$ approximation. In this work we improve the modeling by incorporating realistic photo-$z$ distribution. We show that the position of the BAO scale in $ξ_{\rm p}$ is determined by the photo-$z$ distribution and the Jacobian of the transformation. The latter diverges at the transverse scale of the separation $s_\perp $, and it explains why $ξ_{\rm p } $ traces the underlying correlation function at $s_\perp $, rather than $s$, when the photo-$z$ uncertainty $ σ_z / (1+ z) \gtrsim 0.02$. We also obtain the Gaussian covariance for $ξ_{\rm p}$. Due to photo-$z$ mixing, the covariance of $ξ_{\mathrm{p}}$ shows strong off-diagonal elements. The high correlation of the data causes some issues to the data fitting. Nonetheless, we find that either it can be solved by suppressing the largest eigenvalues of the covariance or it is not directly related to the BAO. We test our BAO fitting pipeline using a set of mock catalogs. The data set is dedicated for Dark Energy Survey Year 3 (DES Y3) BAO analyses and includes realistic photo-$z$ distributions. The theory template is in good agreement with mock measurement. Based on the DES Y3 mocks, $ξ_{\rm p}$ statistic is forecast to constrain the BAO shift parameter $α$ to be $1.001 \pm 0.023$, which is well consistent with the corresponding constraint derived from the angular correlation function measurements. Thus $ξ_{\rm p}$ offers a competitive alternative for the photometric data analyses.

preprint2022arXiv

How the Big Bang Ends up Inside a Black Hole

The standard model of cosmology assumes that our Universe began 14 Gyrs (billion years) ago from a singular Big Bang creation. This can explain a vast range of different astrophysical data from a handful of free cosmological parameters. However, we have no direct evidence or fundamental understanding of some key assumptions: Inflation, Dark Matter and Dark Energy. Here we review the idea that cosmic expansion originates instead from gravitational collapse and bounce. The collapse generates a Black Hole (BH) of mass $ M \simeq 5 \times 10^{22} M_{\odot}$ that formed 25~Gyrs ago. As there is no pressure support, the cold collapse can continue inside in free fall until it reaches atomic nuclear saturation (GeV), when is halted by Quantum Mechanics, as two particles cannot occupy the same quantum state. The collapse then bounces like a core-collapse supernovae, producing the Big Bang expansion. Cosmic acceleration results from the BH event horizon. During collapse, perturbations exit the horizon to re-enter during expansion, giving rise to the observed universe without the need for Inflation or Dark Energy. Using Ockham's razor, this makes the BH Universe (BHU) model more compelling than the standard singular Big Bang creation.

preprint2022arXiv

Modeling Intrinsic Galaxy Alignment in the MICE Simulation

The intrinsic alignment (IA) of galaxies is potentially a major limitation in deriving cosmological constraints from weak lensing surveys. In order to investigate this effect we assign intrinsic shapes and orientations to galaxies in the light-cone output of the MICE simulation, spanning $\sim5000\,{\rm deg}^2$ and reaching redshift $z=1.4$. This assignment is based on a 'semi-analytic' IA model that uses photometric properties of galaxies as well as the spin and shape of their host halos. Advancing on previous work, we include more realistic distributions of galaxy shapes and a luminosity dependent galaxy-halo alignment. The IA model parameters are calibrated against COSMOS and BOSS LOWZ observations. The null detection of IA in observations of blue galaxies is accounted for by setting random orientations for these objects. We compare the two-point alignment statistics measured in the simulation against predictions from the analytical IA models NLA and TATT over a wide range of scales, redshifts and luminosities for red and blue galaxies separately. We find that both models fit the measurements well at scales above $8\,h^{-1}{\rm Mpc}$, while TATT outperforms NLA at smaller scales. The IA parameters derived from our fits are in broad agreement with various observational constraints from red galaxies. Lastly, we build a realistic source sample, mimicking DES Year 3 observations and use it to predict the IA contamination to the observed shear statistics. We find this prediction to be within the measurement uncertainty, which might be a consequence of the random alignment of blue galaxies in the simulation.

preprint2022arXiv

SOAR/Goodman Spectroscopic Assessment of Candidate Counterparts of the LIGO-Virgo Event GW190814

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond.

preprint2022arXiv

The Cosmological Constant as Event Horizon

General Relativity allows for a cosmological constant ($Λ$) which has inspired models of cosmic Inflation and Dark Energy. We show instead that $r_Λ= \sqrt{3/Λ}$ corresponds to an event horizon: a causal boundary term in the action. Our Universe is expanding inside its Schwarzschild radius $r_S=\rL=2GM$, which could have originated from a uniform free falling cloud of mass $M$ that collapsed as a Black Hole (BH) 25 Gyrs ago. Such a BH Universe allows for large-scale structure formation without the need of Inflation or Dark Energy.

preprint2022arXiv

The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim20\Gpchcube$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. \textsc{AbacusSummit} is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\Gpchcube$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-$N$-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k < 0.3\hMpc$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in \textsc{AbacusSummit} suite to extend the effective volume $\sim 20$ times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

preprint2022arXiv

The PAU Survey: Measurements of the 4000 Å spectral break with narrow-band photometry

The D4000 spectral break index is one of the most important features in the visible spectrum, as it is a proxy for stellar ages and is also used in galaxy classification. However, its direct measurement has always been reserved to spectroscopy. Here, we present a general method to directly measure the D4000 with narrow-band (NB) photometry; it has been validated using realistic simulations, and then evaluated with PAUS NBs, cross-matched with VIPERS spectra ($i_{\rm AB} < 22.5$, $0.562 < z < 0.967$). We also reconstruct the D4000 with the SED-fitting code CIGALE; the use of PAUS NBs instead of broad bands significantly improves the SED fitting results. For D4000$_{\rm n}$, the direct measurement has $\rm \langle SNR \rangle \sim 4$, but we find that for $i_{\rm AB}<21$ all direct D4000 measurements have $\rm SNR>3$. The CIGALE D4000$_{\rm n}$ has $\rm \langle SNR \rangle \sim 20$, but underestimates the error by $>$50\%. Furthermore, the direct method recreates well the D4000-SFR relation, as well as the D4000-mass relation for blue galaxies (for red galaxies, selection effects impact the results). On the other hand, CIGALE accurately classifies galaxies into red and blue populations. We conclude that the direct measurement of D4000 with narrow-band photometry is a promising tool to determine average properties of galaxy samples, with results compatible with spectroscopy.

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 Cosmological Constant as a Zero Action Boundary

The cosmological constant $Λ$ is usually interpreted as Dark Energy (DE) or modified gravity (MG). Here we propose instead that $Λ$ corresponds to a boundary term in the action of classical General Relativity. The action is zero for a perfect fluid solution and this fixes $Λ$ to the average density $ρ$ and pressure $p$ inside a primordial causal boundary: $Λ= 4πG <ρ+3p>$. This explains both why the observed value of $Λ$ is related to the matter density today and also why other contributions to $Λ$, such as DE or MG, do not produce cosmic expansion. Cosmic acceleration results from the repulsive boundary force that occurs when the expansion reaches the causal horizon. This universe is similar to the $Λ$CDM universe, except on the largest observable scales, where we expect departures from homogeneity/isotropy, such as CMB anomalies and variations in cosmological parameters indicated by recent observations.

preprint2021arXiv

The PAU Survey: Intrinsic alignments and clustering of narrow-band photometric galaxies

We present the first measurements of the projected clustering and intrinsic alignments (IA) of galaxies observed by the Physics of the Accelerating Universe Survey (PAUS). With photometry in 40 narrow optical passbands ($450\rm{nm}-850\rm{nm}$), the quality of photometric redshift estimation is $σ_{z} \sim 0.01(1 + z)$ for galaxies in the $19\,\rm{deg}^{2}$ Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) W3 field, allowing us to measure the projected 3D clustering and IA for flux-limited, faint galaxies ($i < 22.5$) out to $z\sim0.8$. To measure two-point statistics, we developed, and tested with mock photometric redshift samples, `cloned' random galaxy catalogues which can reproduce data selection functions in 3D and account for photometric redshift errors. In our fiducial colour-split analysis, we made robust null detections of IA for blue galaxies and tentative detections of radial alignments for red galaxies ($\sim1-3σ$), over scales of $0.1-18\,h^{-1}\rm{Mpc}$. The galaxy clustering correlation functions in the PAUS samples are comparable to their counterparts in a spectroscopic population from the Galaxy and Mass Assembly survey, modulo the impact of photometric redshift uncertainty which tends to flatten the blue galaxy correlation function, whilst steepening that of red galaxies. We investigate the sensitivity of our correlation function measurements to choices in the random catalogue creation and the galaxy pair-binning along the line of sight, in preparation for an optimised analysis over the full PAUS area.

preprint2021arXiv

The PAU survey: Ly$α$ intensity mapping forecast

In this work, we explore the application of intensity mapping to detect extended Ly$α$ emission from the IGM via cross-correlation of PAUS images with Ly$α$ forest data from eBOSS and DESI. Seven narrow-band (FWHM=13nm) PAUS filters have been considered, ranging from 455 to 515 nm in steps of 10 nm, which allows the observation of Ly$α$ emission in a range $2.7<z<3.3$. The cross-correlation is simulated first in an area of 100 deg$^2$ (PAUS projected coverage), and second in two hypothetical scenarios: a deeper PAUS (complete up to $i_{\rm AB}<24$ instead of $i_{\rm AB}<23$, observation time x6), and an extended PAUS coverage of 225 deg$^2$ (observation time x2.25). A hydrodynamic simulation of size 400 Mpc/h is used to simulate both extended Ly$α$ emission and absorption, while the foregrounds in PAUS images have been simulated using a lightcone mock catalogue. Using an optimistic estimation of uncorrelated PAUS noise, the total probability of a non-spurious detection is estimated to be 1.8\% and 4.5\% for PAUS-eBOSS and PAUS-DESI , from a run of 1000 simulated cross-correlations with different realisations of instrumental noise and quasar positions. The hypothetical PAUS scenarios increase this probability to 15.3\% (deeper PAUS) and 9.0\% (extended PAUS). With realistic correlated noise directly measured from PAUS images, these probabilities become negligible. Despite these negative results, some evidences suggest that this methodology may be more suitable to broad-band surveys.

preprint2020arXiv

Dark Energy Survey Identification of A Low-Mass Active Galactic Nucleus at Redshift 0.823 from Optical Variability

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$β$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{\rm 2-12\,keV}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$β$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$σ$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

preprint2020arXiv

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

preprint2020arXiv

The Curious Case of PHL 293B: A Long-Lived Transient in a Metal-Poor Blue Compact Dwarf Galaxy

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

preprint2020arXiv

The PAU Survey: An improved photo-$z$ sample in the COSMOS field

We present -- and make publicly available -- accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with $i_{\mathrm{AB}}\leq 23$. The redshifts are computed using a combination of narrow band photometry from PAUS, a survey with 40 narrow bands spaced at $100Å$ intervals covering the range from $4500Å$ to $8500Å$, and 26 broad, intermediate, and narrow bands covering the UV, visible and near infrared spectrum from the COSMOS2015 catalogue. We introduce a new method that models the spectral energy distributions (SEDs) as a linear combination of continuum and emission line templates and computes its Bayes evidence, integrating over the linear combinations. The correlation between the UV luminosity and the OII line is measured using the 66 available bands with the zCOSMOS spectroscopic sample, and used as a prior which constrains the relative flux between continuum and emission line templates. The flux ratios between the OII line and $\mathrm{H}_α$, $\mathrm{H}_β$ and $\mathrm{OIII}$ are similarly measured and used to generate the emission line templates. Comparing to public spectroscopic surveys via the quantity $Δ_z\equiv(z_{\mathrm{photo}}-z_{\mathrm{spec}})/(1+z_{\mathrm{spec}})$, we find the photometric redshifts to be more precise than previous estimates, with $σ_{68}(Δ_z) \approx (0.003, 0.009)$ for galaxies at magnitude $i_{\mathrm{AB}}\sim18$ and $i_{\mathrm{AB}}\sim23$, respectively, which is $3\times$ and $1.66\times$ tighter than COSMOS2015. Additionally, we find the redshifts to be very accurate on average, yielding a median of the $Δ_z$ distribution compatible with $|\mathrm{median}(Δ_z)|\leq0.001$ at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results.

preprint2020arXiv

The PAU Survey: Background light estimation with deep learning techniques

In any imaging survey, measuring accurately the astronomical background light is crucial to obtain good photometry. This paper introduces BKGnet, a deep neural network to predict the background and its associated error. BKGnet has been developed for data from the Physics of the Accelerating Universe Survey (PAUS), an imaging survey using a 40 narrow-band filter camera (PAUCam). Images obtained with PAUCam are affected by scattered light: an optical effect consisting of light multiply that deposits energy in specific detector regions contaminating the science measurements. Fortunately, scattered light is not a random effect, but it can be predicted and corrected for. We have found that BKGnet background predictions are very robust to distorting effects, while still being statistically accurate. On average, the use of BKGnet improves the photometric flux measurements by 7% and up to 20% at the bright end. BKGnet also removes a systematic trend in the background error estimation with magnitude in the i-band that is present with the current PAU data management method. With BKGnet, we reduce the photometric redshift outlier rate

preprint2020arXiv

The size of our causal Universe

A Universe with finite age also has a finite causal scale. Larger scales can not affect our local measurements or modeling, but far away locations could have different cosmological parameters. The size of our causal Universe depends on the details of inflation and is usually assumed to be larger than our observable Universe today. To account for causality, we propose a new boundary condition, that can be fulfill by fixing the cosmological constant (a free geometric parameter of gravity). This forces a cancellation of vacuum energy with the cosmological constant. As a consequence, the measured cosmic acceleration can not be explained by a simple cosmological constant or constant vacuum energy. We need some additional odd properties such as the existence of evolving dark energy (DE) with energy-density fine tuned to be twice that of dark matter today. We show here that we can instead explain cosmic acceleration without DE (or modified gravity) assuming that the causal scale is smaller than the observable Universe today. Such scale corresponds to half the sky at z=1 and 60 degrees at z=1100, which is consistent with the anomalous lack of correlations observed in the CMB. Late time cosmic acceleration could then be interpreted as the smoking gun of primordial Inflation.

preprint2019arXiv

Spectral Variability of a Sample of Extreme Variability Quasars and Implications for the MgII Broad-line Region

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed by more than 1.5 mag in the $g$ band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey (DES) epochs (separated by a few years in the quasar rest frame). The quasar sample covers a redshift range of $0.5 < z < 2.1$. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the $g$ band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on timescales of only months. Leveraging on the large dynamic range in continuum variability between the earlier SDSS and the new GMOS spectra, we explore the associated variations in the broad Mg II,$\lambda2798$ line, whose variability properties have not been well studied before. The broad Mg II flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the width (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad Mg II to estimate the black hole mass therefore introduces a luminosity-dependent bias.

preprint2018arXiv

The PAU Survey: A Forward Modeling Approach for Narrow-band Imaging

Weak gravitational lensing is a powerful probe of the dark sector, once measurement systematic errors can be controlled. In Refregier & Amara (2014), a calibration method based on forward modeling, called MCCL, was proposed. This relies on fast image simulations (e.g., UFig; Berge et al. 2013) that capture the key features of galaxy populations and measurement effects. The MCCL approach has been used in Herbel et al. (2017) to determine the redshift distribution of cosmological galaxy samples and, in the process, the authors derived a model for the galaxy population mainly based on broad-band photometry. Here, we test this model by forward modeling the 40 narrow-band photometry given by the novel PAU Survey (PAUS). For this purpose, we apply the same forced photometric pipeline on data and simulations using Source Extractor (Bertin & Arnouts 1996). The image simulation scheme performance is assessed at the image and at the catalogues level. We find good agreement for the distribution of pixel values, the magnitudes, in the magnitude-size relation and the interband correlations. A principal component analysis is then performed, in order to derive a global comparison of the narrow-band photometry between the data and the simulations. We use a `mixing' matrix to quantify the agreement between the observed and simulated sets of Principal Components (PCs). We find good agreement, especially for the first three most significant PCs. We also compare the coefficients of the PCs decomposition. While there are slight differences for some coefficients, we find that the distributions are in good agreement. Together, our results show that the galaxy population model derived from broad-band photometry is in good overall agreement with the PAUS data. This offers good prospect for incorporating spectral information to the galaxy model by adjusting it to the PAUS narrow-band data using forward modeling.

preprint2016arXiv

Cosmology from large scale galaxy clustering and galaxy-galaxy lensing with Dark Energy Survey Science Verification data

We present cosmological constraints from the Dark Energy Survey (DES) using a combined analysis of angular clustering of red galaxies and their cross-correlation with weak gravitational lensing of background galaxies. We use a 139 square degree contiguous patch of DES data from the Science Verification (SV) period of observations. Using large scale measurements, we constrain the matter density of the Universe as Omega_m = 0.31 +/- 0.09 and the clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/- 0.13 after marginalizing over seven nuisance parameters and three additional cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16} = 0.74 +/- 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 = 0.78 +/- 0.09 using two bins over the range 0.2 <z< 0.5. We study the robustness of the results under changes in the data vectors, modelling and systematics treatment, including photometric redshift and shear calibration uncertainties, and find consistency in the derived cosmological parameters. We show that our results are consistent with previous cosmological analyses from DES and other data sets and conclude with a joint analysis of DES angular clustering and galaxy-galaxy lensing with Planck CMB data, Baryon Accoustic Oscillations and Supernova type Ia measurements.

preprint2016arXiv

Host Galaxy Identification for Supernova Surveys

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.

preprint2016arXiv

Linear and non-linear bias: predictions vs. measurements

We study the linear and non-linear bias parameters which determine the mapping between the distributions of galaxies and the full matter density fields, comparing different measurements and predictions. Associating galaxies with dark matter haloes in the MICE Grand Challenge N-body simulation we directly measure the bias parameters by comparing the smoothed density fluctuations of haloes and matter in the same region at different positions as a function of smoothing scale. Alternatively we measure the bias parameters by matching the probability distributions of halo and matter density fluctuations, which can be applied to observations. These direct bias measurements are compared to corresponding measurements from two-point and different third-order correlations, as well as predictions from the peak-background model, which we presented in previous articles using the same data. We find an overall variation of the linear bias measurements and predictions of $\sim 5 \%$ with respect to results from two-point correlations for different halo samples with masses between $\sim 10^{12} - 10^{15}$ $h^{-1}M_\odot$ at the redshifts $z=0.0$ and $0.5$. Variations between the second- and third-order bias parameters from the different methods show larger variations, but with consistent trends in mass and redshift. The various bias measurements reveal a tight relation between the linear and the quadratic bias parameters, which is consistent with results from the literature based on simulations with different cosmologies. Such a universal relation might improve constraints on cosmological models, derived from second-order clustering statistics at small scales or higher-order clustering statistics.

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.

preprint2015arXiv

Combining spectroscopic and photometric surveys using angular cross-correlations II: Parameter constraints from different physical effects

Future spectroscopic and photometric surveys will measure accurate positions and shapes of an increasing number of galaxies. In the previous paper of this series we studied the effects of Redshift Space Distortions (RSD), baryon acoustic oscillations (BAO) and Weak gravitational Lensing (WL) using angular cross-correlation. Here, we provide a new forecast that explores the contribution of including different observables, physical effects (galaxy bias, WL, RSD, BAO) and approximations (non-linearities, Limber approximation, covariance between probes). The radial information is included by using the cross-correlation of separate narrow redshift bins. For the auto correlation the separation of galaxy pairs is mostly transverse, while the cross-correlations also includes a radial component. We study how this information adds to our figure of merit (FoM), which includes the dark energy equation of state $w(z)$ and the growth history, parameterized by $γ$. We show that the Limber approximation and galaxy bias are the most critical ingredients to the modelling of correlations. Adding WL increases our FoM by 4.8, RSD by 2.1 and BAO by 1.3. We also explore how overlapping surveys perform under the different assumption and for different figures of merit. Our qualitative conclusions depend on the survey choices and scales included, but we find some clear tendencies that highlight the importance of combining different probes and can be used to guide and optimise survey strategies.

preprint2015arXiv

Combining Spectroscopic and Photometric Surveys: Same or different sky?

This article looks at the combined constraints from a photometric and spectroscopic survey. These surveys will measure cosmology using weak lensing (WL), galaxy clustering, baryon acoustic oscillations (BAO) and redshift space distortions (RSD). We find, contrary to some findings in the recent literature, that overlapping surveys can give important benefits when measuring dark energy. We therefore try to clarify the status of this issue with a full forecast of two stage-IV surveys using a new approach to properly account for covariance between the different probes in the overlapping samples. The benefit of the overlapping survey can be traced back to two factors: additional observables and sample variance cancellation. Both needs to be taken into account and contribute equally when combining 3D power spectrum and 2D correlations for lensing. With an analytic example we also illustrate that for optimal constraints, one should minimize the (Pearson) correlation coefficient between cosmological and nuisance parameters and maximize the one among nuisance parameters (e.g. galaxy bias) in the two samples. This can be achieved by increasing the overlap between the spectroscopic and photometric surveys. We show how BAO, WL and RSD contribute to this benefit and also look at some other survey designs, such as photometric redshift errors and spectroscopic density.

preprint2015arXiv

Comparing halo bias from abundance and clustering

We model the abundance of haloes in the $\sim(3 \ \text{Gpc}/h)^3$ volume of the MICE Grand Challenge simulation by fitting the universal mass function with an improved Jack-Knife error covariance estimator that matches theory predictions. We present unifying relations between different fitting models and new predictions for linear ($b_1$) and non-linear ($c_2$ and $c_3$) halo clustering bias. Different mass function fits show strong variations in their performance when including the low mass range ($M_h \lesssim 3 \ 10^{12} \ M_{\odot}/h$) in the analysis. Together with fits from the literature we find an overall variation in the amplitudes of around $10$% in the low mass and up to $50$% in the high mass (galaxy cluster) range ($M_h > 10^{14} \ M_{\odot}/h$). These variations propagate into a $10$% change in $b_1$ predictions and a $50$% change in $c_2$ or $c_3$. Despite these strong variations we find universal relations between $b_1$ and $c_2$ or $c_3$ for which we provide simple fits. Excluding low mass haloes, different models fitted with reasonable goodness in this analysis, show percent level agreement in their $b_1$ predictions, but are systematically $5-10$% lower than the bias directly measured with two-point halo-mass clustering. This result confirms previous findings derived from smaller volumes (and smaller masses). Inaccuracies in the bias predictions lead to $5-10$% errors in growth measurements. They also affect any HOD fitting or (cluster) mass calibration from clustering measurements.

preprint2014arXiv

Are the halo occupation predictions consistent with large scale galaxy clustering?

We study how well we can reconstruct the 2-point clustering of galaxies on linear scales, as a function of mass and luminosity, using the halo occupation distribution (HOD) in several semi-analytical models (SAMs) of galaxy formation from the Millennium Simulation. We find that HOD with Friends of Friends groups can reproduce galaxy clustering better than gravitationally bound haloes. This indicates that Friends of Friends groups are more directly related to the clustering of these regions than the bound particles of the overdensities. In general we find that the reconstruction works at best to 5% accuracy: it underestimates the bias for bright galaxies. This translates to an overestimation of 50% in the halo mass when we use clustering to calibrate mass. We also found a degeneracy on the mass prediction from the clustering amplitude that affects all the masses. This effect is due to the clustering dependence on the host halo substructure, an indication of assembly bias. We show that the clustering of haloes of a given mass increases with the number of subhaloes, a result that only depends on the underlying matter distribution. As the number of galaxies increases with the number of subhaloes in SAMs, this results in a low bias for the HOD reconstruction. We expect this effect to apply to other models of galaxy formation, including the real universe, as long as the number of galaxies incresases with the number of subhaloes. We have also found that the reconstructions of galaxy bias from the HOD model fails for low mass haloes with M = 3-5x10^11 Msun/h. We find that this is because galaxy clustering is more strongly affected by assembly bias for these low masses.

preprint2014arXiv

Asymmetric galaxy correlation functions

We study the two-point cross-correlation function between two populations of galaxies: for instance a bright population and a faint population. We show that this cross-correlation is asymmetric under the exchange of the line-of-sight coordinate of the galaxies, i.e. that the correlation is different if the bright galaxy is in front of, or behind, the faint galaxy. We give an intuitive, quasi-Newtonian derivation of all the effects that contribute to such an asymmetry in large-scale structure: gravitational redshift, Doppler shift, lensing, light-cone, evolution and Alcock-Paczynski effects - interestingly, the gravitational redshift term is exactly canceled by some of the others, assuming geodesic motion. Most of these effects are captured by previous calculations of general relativistic corrections to the observed galaxy density fluctuation; the asymmetry arises from terms that are suppressed by the ratio H/k - H is the Hubble constant and k is the wavenumber - which are more readily observable than the terms suppressed by (H/k)^2. Some of the contributions to the asymmetry, however, arise from terms that are generally considered 'Newtonian' - the lensing and evolution - and thus represent a contaminant in the search for general relativistic corrections. We propose methods to disentangle these different contributions. A simple method reduces the contamination to a level of < 10% for redshifts z<1. We also clarify the relation to recent work on measuring gravitational redshifts by stacking clusters.

preprint2014arXiv

Combining spectroscopic and photometric surveys using angular cross-correlations I: Algorithm and modelling

Weak lensing (WL) clustering is studied using 2D (angular) coordinates, while redshift space distortions (RSD) and baryon acoustic oscillations (BAO) use 3D coordinates, which requires a model dependent conversion of angles and redshifts into comoving distances. This is the first paper of a series, which explore modelling multi-tracer galaxy clustering (of WL, BAO and RSD), using only angular (2D) cross-correlations in thin redshift bins. This involves evaluating many thousands cross-correlations, each a multidimensional integral, which is computationally demanding. We present a new algorithm that performs these calculations as matrix operations. Nearby narrow redshift bins are intrinsically correlated, which can be used to recover the full (radial) 3D information. We show that the Limber approximation does not work well for this task. In the exact calculation, both the clustering amplitude and the RSD effect increase when decreasing the redshift bin width. For narrow bins, the cross-correlations has a larger BAO peak than the auto-correlation because smaller scales are filtered out by the radial redshift separation. Moreover, the BAO peak shows a second (ghost) peak, shifted to smaller angles. We explore how WL, RSD and BAO contribute to the cross-correlations as a function of the redshift bin width and present a first exploration of non-linear effects and signal-to-noise ratio on these quantities. This illustrates that the new approach to clustering analysis provides new insights and is potentially viable in practice.

preprint2014arXiv

Measuring the growth of matter fluctuations with third-order galaxy correlations

Measurements of the linear growth factor $D$ at different redshifts $z$ are key to distinguish among cosmological models. One can estimate the derivative $dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy two-point correlation $ξ(z)$, but the degeneracy of its transverse (or projected) component with galaxy bias $b$, i.e. $ξ_{\perp}(z) \propto\ D^2(z) b^2(z)$, introduces large errors in the growth measurement. Here we present a comparison between two methods which break this degeneracy by combining second- and third-order statistics. One uses the shape of the reduced three-point correlation and the other a combination of third-order one- and two-point cumulants. These methods use the fact that, for Gaussian initial conditions and scales larger than $20$ $h^{-1}$Mpc, the reduced third-order matter correlations are independent of redshift (and therefore of the growth factor) while the third-order galaxy correlations depend on $b$. We use matter and halo catalogs from the MICE-GC simulation to test how well we can recover $b(z)$ and therefore $D(z)$ with these methods in 3D real space. We also present a new approach, which enables us to measure $D$ directly from the redshift evolution of second- and third-order galaxy correlations without the need of modelling matter correlations. For haloes with masses lower than $10^{14}$ $h^{-1}$M$_\odot$, we find $10%$ deviations between the different estimates of $D$, which are comparable to current observational errors. At higher masses we find larger differences that can probably be attributed to the breakdown of the bias model and non-Poissonian shot noise.

preprint2014arXiv

Redshift-space distortions from the cross-correlation of photometric populations

Several papers have recently highlighted the possibility of measuring redshift space distortions from angular auto-correlations of galaxies in photometric redshift bins. In this work we extend this idea to include as observables the cross-correlations between redshift bins, as an additional way of measuring radial information. We show that this extra information allows to reduce the recovered error in the growth rate index γby a factor of ~2. Although the final error in γdepends on the bias and the mean photometric accuracy of the galaxy sample, the improvement from adding cross-correlations is robust in different settings. Another factor of 2-3 improvement in the determination of γcan be achieved by considering two galaxy populations over the same photometric sky area but with different biases. This additional gain is shown to be much larger than the one from the same populations when observed over different areas of the sky (with twice the combined area). The total improvement of ~5 implies that a photometric survey such as the Dark Energy Survey should be able to recover γat the 5-10% from the angular clustering in linear scales of two different tracers. It can also constrain the evolution of f(z)xσ_8(z) in few bins beyond z~0.8-0.9 at the 10-15% level per-bin, compatible with recent constrains from lower-z spectroscopic surveys. We also show how further improvement can be achieved by reducing the photometric redshift error.

preprint2014arXiv

Subhaloes gone Notts: Subhaloes as tracers of the dark matter halo shape

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

preprint2014arXiv

Subhaloes gone Notts: the clustering properties of subhaloes

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

preprint2012arXiv

Recovering 3D clustering information with angular correlations

We study how to recover the full 3D clustering information of P(\vec{k},z), including redshift space distortions (RSD), from 2D tomography using the angular auto and cross spectra of different redshift bins C_\ell(z,z'). We focus on quasilinear scales where the minimum scale λ_{min} or corresponding maximum wavenumber k_{max}= 2π/λ_{min} is targeted to be between k_{max}={0.05-0.2} h/Mpc. For spectroscopic surveys, we find that we can recover the full 3D clustering information when the redshift bin width Δz used in the 2D tomography is similar to the targeted minimum scale, i.e. Δz ~ {0.6-0.8} λ_{min} H(z)/c which corresponds to Δz ~ 0.01-0.05 for z<1. This value of Δz is optimal in the sense that larger values of Δz lose information, while smaller values violate our minimum scale requirement. For a narrow-band photometric survey, with photo-z error σ_z=0.004, we find almost identical results to the spectroscopic survey because the photo-z error is smaller than the optimal bin width σ_z<Δz. For a typical broad-band photometric survey with σ_z=0.1, we have that σ_z>Δz and most radial information is intrinsically lost. The remaining information can be recovered from the 2D tomography if we use Δz ~ 2σ_z. While 3D and 2D analysis are shown here to be equivalent, the advantage of using angular positions and redshifts is that we do not need a fiducial cosmology to convert to 3D coordinates. This avoids assumptions and marginalization over the fiducial model. In addition, it becomes straight forward to combine RSD, clustering and weak lensing in 2D space.

preprint2011arXiv

Clustering of photometric luminous red galaxies I : Growth of Structure and Baryon Acoustic Feature

The possibility of measuring redshift space (RSD) distortions using photometric data have been recently highlighted. This effect complements and significantly alters the detectability of baryon acoustic oscillations (BAO) in photometric surveys. In this paper we present measurements of the angular correlation function of luminous red galaxies (LRGs) in the photometric catalog of the final data release (DR7) of the Sloan Digital Sky Survey II (SDSS). The sample compromise ~ 1.5 x 10^6 LRGs distributed in 0.45 < z < 0.65, with a characteristic photometric error of ~ 0.05. Our measured correlation centered at z=0.55 is in very good agreement with predictions from standard LCDM in a broad range of angular scales, $0.5^\circ < θ< 6^\circ$. We find that the growth of structure can indeed be robustly measured, with errors matching expectations. The velocity growth rate is recovered as $f σ_8 = 0.53 \pm 0.42$ when no prior is imposed on the growth factor and the background geometry follows a LCDM model with WMAP7+SNIa priors. This is compatible with the corresponding General Relativity (GR) prediction $f σ_8 = 0.45$ for our fiducial cosmology. If we adopt a parametrization such that $f=Ω^γ_m(z)$, with $γ\approx 0.55$ in GR, and combine our $fσ_8$ measurement with the corresponding ones from spectroscopic LRGs at lower redshifts we obtain $γ=0.54 \pm 0.17$. In addition we find evidence for the presence of the baryon acoustic feature matching the amplitude, location and shape of LCDM predictions. The photometric BAO feature is detected with 98 % confidence level at z=0.55.

preprint2011arXiv

Have Baryonic Acoustic Oscillations in the galaxy distribution really been measured?

Recent publications claim that there is no convincing evidence for measurements of the baryonic acoustic (BAO) feature in galaxy samples using either monopole or radial information. Different claims seem contradictory: data is either not consistent with the BAO model or data is consistent with both the BAO model and featureless models without BAO. We investigate this point with a set of 216 realistic mock galaxy catalogs extracted from MICE7680, one of the largest volume dark matter simulation run to date, with a volume of 1300 cubical gigaparsecs. Our mocks cover similar volume, densities and bias as the real galaxies and provide 216 realizations of the Lambda or w=-1 Cold Dark Matter (wCDM) BAO model. We find that only 20% of the mocks show a statistically significant (3 sigma) preference for the true (input) wCDM BAO model as compared to a featureless (non-physical) model without BAO. Thus the volume of current galaxy samples is not yet large enough to claim that the BAO feature has been detected. Does this mean that we can not locate the BAO position? Using a simple (non optimal) algorithm we show that in 50% (100%) of the mocks we can find the BAO position within 5% (20%) of the true value. These two findings are not in contradiction: the former is about model selection, the later is about parameter fitting within a model. We conclude that current monopole and radial BAO measurements can be used as standard rulers if we assume wCDM type of models.

preprint2011arXiv

Measuring Redshift-Space Distortions using Photometric Surveys

We outline how redshift-space distortions (RSD) can be measured from the angular correlation function w(θ), of galaxies selected from photometric surveys. The natural degeneracy between RSD and galaxy bias can be minimized by comparing results from bins with top-hat galaxy selection in redshift, and bins based on the radial position of galaxy pair centres. This comparison can also be used to test the accuracy of the photometric redshifts. The presence of RSD will be clearly detectable with the next generation of photometric redshift surveys. We show that the Dark Energy Survey (DES) will be able to measure f(z)σ_8(z) to a 1σ accuracy of (17 {\times} b)%, using galaxies drawn from a single narrow redshift slice centered at z = 1. Here b is the linear bias, and f is the logarithmic rate of change of the linear growth rate with respect to the scale factor. Extending to measurements of w(θ) for a series of bins of width 0.02(1 + z) over 0.5 < z < 1.4 will measure γ to a 1σ accuracy of 25%, given the model f = Ω_m(z)^γ, and assuming a linear bias model that evolves such that b = 0.5 + z (and fixing other cosmological parameters). The accuracy of our analytic predictions is confirmed using mock catalogs drawn from simulations conducted by the MICE collaboration.

preprint2011arXiv

Statistical Analysis of Galaxy Surveys-IV: An objective way to quantify the impact of superstructures on galaxy clustering statistics

For galaxy clustering to provide robust constraints on cosmological parameters and galaxy formation models, it is essential to make reliable estimates of the errors on clustering measurements. We present a new technique, based on a spatial Jackknife (JK) resampling, which provides an objective way to estimate errors on clustering statistics. Our approach allows us to set the appropriate size for the Jackknife subsamples. The method also provides a means to assess the impact of individual regions on the measured clustering, and thereby to establish whether or not a given galaxy catalogue is dominated by one or several large structures, preventing it to be considered as a "fair sample". We apply this methodology to the two- and three-point correlation functions measured from a volume limited sample of M* galaxies drawn from data release seven of the Sloan Digital Sky Survey (SDSS). The frequency of jackknife subsample outliers in the data is shown to be consistent with that seen in large N-body simulations of clustering in the cosmological constant plus cold dark matter cosmology. We also present a comparison of the three-point correlation function in SDSS and 2dFGRS using this approach and find consistent measurements between the two samples.

preprint2011arXiv

The Local Bias Model in the Large Scale Halo Distribution

We explore the biasing in the clustering statistics of halos as compared to dark matter (DM) in simulations. We look at the second and third order statistics at large scales of the (intermediate) MICEL1536 simulation and also measure directly the local bias relation h = f(δ) between DM fluctuations, δ, smoothed over a top-hat radius Rs at a point in the simulation and its corresponding tracer h (i.e. halos) at the same point. This local relation can be Taylor expanded to define a linear (b1) and non-linear (b2) bias parameters. The values of b1 and b2 in the simulation vary with Rs approaching a constant value around Rs > 30 - 60 Mpc/h. We use the local relation to predict the clustering of the tracer in terms of the one of DM. This prediction works very well (about percent level) for the halo 2-point correlation ξ(r_12) for r_12 > 15 Mpc/h, but only when we use the biasing values that we found at very large smoothing radius Rs > 30 - 60 Mpc/h. We find no effect from stochastic or next to leading order terms in the f(δ) expansion. But we do find some discrepancies in the 3-point function that needs further understanding. We also look at the clustering of the smoothed moments, the variance and skewness which are volume average correlations and therefore include clustering from smaller scales. In this case, we find that both next to leading order and discreetness corrections (to the local model) are needed at the 10 - 20% level. Shot-noise can be corrected with a term σe^2/n where σe^2 < 1, i.e., always smaller than the Poisson correction. We also compare these results with the peak-background split predictions from the measured halo mass function. We find 5-10% systematic (and similar statistical) errors in the mass estimation when we use the halo model biasing predictions to calibrate the mass.

preprint2009arXiv

Clustering of luminous red galaxies II: small scale redshift space distortions

This is the second paper of a series where we study the clustering of LRG galaxies in the latest spectroscopic SDSS data release, DR6, which has 75000 LRG galaxies covering over 1 $Gpc^3/h^3$ for $0.15<z<0.47$. Here we focus on modeling redshift space distortions in $\xisp$, the 2-point correlation in separate line-of-sight and perpendicular directions, at small scales and in the line-of-sight. We show that a simple Kaiser model for the anisotropic 2-point correlation function in redshift space, convolved with a distribution of random peculiar velocities with an exponential form, can describe well the correlation of LRG at all scales. We show that to describe with accuracy the so called "fingers-of-God" (FOG) elongations in the radial direction, it is necessary to model the scale dependence of both bias $b$ and the pairwise rms peculiar velocity $σ_{12}$ with the distance. We show how both quantities can be inferred from the $\xisp$ data. From $r \simeq 10$ Mpc/h to $r \simeq 1$ Mpc/h, both the bias and $σ_{12}$ are shown to increase by a factor of two: from $b=2$ to $b=4$ and from $σ_{12}=400$ to 800 Km/s. The later is in good agreement, within a 5 percent accuracy in the recovered velocities, with direct velocity measurements in dark matter simulations with $Ω_m=0.25$ and $σ_8$=0.85.

preprint2009arXiv

Clustering of Luminous Red Galaxies IV: Baryon Acoustic Peak in the Line-of-Sight Direction and a Direct Measurement of H(z)

We study the clustering of LRG galaxies in the latest spectroscopic SDSS data releases, DR6 and DR7, which sample over 1 Gpc^3/h^3 to z=0.47. The 2-point correlation function $\xisp$ is estimated as a function of perpendicular $σ$ and line-of-sight $π$ (radial) directions. We find a significant detection of a peak at $r\simeq 110$Mpc/h, which shows as a circular ring in the $σ-π$ plane. There is also significant evidence for a peak along the radial direction whose shape is consistent with its originating from the recombination-epoch baryon acoustic oscillations (BAO). A $\xisp$ model with no radial BAO peak is disfavored at $3.2σ$, whereas a model with no magnification bias is disfavored at $2σ$. The radial data enable, for the first time, a direct measurement of the Hubble parameter $H(z)$ as a function of redshift. This is independent from earlier BAO measurements which used the spherically averaged (monopole) correlation to constrain an integral of $H(z)$. Using the BAO peak position as a standard ruler in the radial direction, we find: $H(z=0.24)= 79.69 \pm 2.32 (\pm 1.29)$ km/s/Mpc for z=0.15-0.30 and $H(z=0.43)= 86.45 \pm 3.27 (\pm 1.69)$ km/s/Mpc for $z=0.40-0.47$. The first error is a model independent statistical estimation and the second accounts for systematics both in the measurements and in the model. For the full sample, $z=0.15-0.47$, we find $H(z=0.34)= 83.80 \pm 2.96 (\pm 1.59)$ km/s/Mpc.

preprint2009arXiv

Simulating the Universe with MICE: The abundance of massive clusters

We introduce a new set of large N-body runs, the MICE simulations, that provide a unique combination of very large cosmological volumes with good mass resolution. They follow the gravitational evolution of ~ 8.5 billion particles (2048^3) in volumes covering up to 450 (Gpc/h)^3. Our main goal is to accurately model and calibrate basic cosmological probes that will be used by upcoming astronomical surveys. Here we take advantage of the very large volumes of MICE to make a robust sampling of the high-mass tail of the halo mass function (MF). We discuss and avoid possible systematic effects in our study, and do a detailed analysis of different error estimators. We find that available fits to the local abundance of halos (Warren et al. (2006)) match well the abundance in MICE up to M ~ 10^{14}\Msun, but significantly deviate for larger masses, underestimating the mass function by 10% (30%) at M = 3.16 x 10^{14}\Msun (10^{15}\Msun). Similarly, the widely used Sheth & Tormen (1999) fit, if extrapolated to high redshift assuming universality, leads to an underestimation of the cluster abundance by 30%, 20% and 15% at z=0, 0.5, 1 for M ~ [7 - 2.5 - 0.8] x 10^{14}\Msun respectively ($ν= δ_c/σ~ 3$). We provide a re-calibration of the halo MF valid over 5 orders of magnitude in mass, 10^{10} < M/(\Msun) < 10^{15}, that accurately describes its redshift evolution up to z=1. We explore the impact of this re-calibration on the determination of dark-energy, and conclude that using available fits may systematically bias the estimate of w by as much as 50% for medium-depth (z <= 1) surveys. MICE halo catalogues are publicly available at http://www.ice.cat/mice

preprint2008arXiv

Statistical Analysis of Galaxy Surveys - I. Robust error estimation for 2-point clustering statistics

We present a test of different error estimators for 2-point clustering statistics, appropriate for present and future large galaxy redshift surveys. Using an ensemble of very large dark matter LambdaCDM N-body simulations, we compare internal error estimators (jackknife and bootstrap) to external ones (Monte-Carlo realizations). For 3-dimensional clustering statistics, we find that none of the internal error methods investigated are able to reproduce neither accurately nor robustly the errors of external estimators on 1 to 25 Mpc/h scales. The standard bootstrap overestimates the variance of xi(s) by ~40% on all scales probed, but recovers, in a robust fashion, the principal eigenvectors of the underlying covariance matrix. The jackknife returns the correct variance on large scales, but significantly overestimates it on smaller scales. This scale dependence in the jackknife affects the recovered eigenvectors, which tend to disagree on small scales with the external estimates. Our results have important implications for the use of galaxy clustering in placing constraints on cosmological parameters. For example, in a 2-parameter fit to the projected correlation function, we find that the standard bootstrap systematically overestimates the 95% confidence interval, while the jackknife method remains biased, but to a lesser extent. The scatter we find between realizations, for Gaussian statistics, implies that a 2-sigma confidence interval, as inferred from an internal estimator, could correspond in practice to anything from 1-sigma to 3-sigma. Finally, by an oversampling of sub-volumes, it is possible to obtain bootstrap variances and confidence intervals that agree with external error estimates, but it is not clear if this prescription will work for a general case.

preprint2001arXiv

Constraints on the Accuracy of Photometric Redshifts Derived from BLAST and Herschel/SPIRE Sub-mm Surveys

More than 150 galaxies have been detected in blank-field millimetre and sub-millimetre surveys. However the redshift distribution of sub-mm galaxies remains uncertain due to the difficulty in identifying their optical-IR counterparts, and subsequently obtaining their spectroscopic emission-line redshifts. In this paper we discuss results from a Monte-Carlo analysis of the accuracy with which one can determine redshifts from photometric measurements at sub-millimetre-FIR wavelengths. The analysis takes into account the dispersion in colours introduced by including galaxies with a distribution of SEDs, and by including photometric and absolute calibration errors associated with real observations. We present examples of the probability distribution of redshifts for individual galaxies detected in the future BLAST and Herschel/SPIRE surveys. We show that the combination of BLAST and 850um observations constrain the photometric redshifts with sufficient accuracy to pursue a program of spectroscopic follow-up with the 100m GBT.

preprint2001arXiv

The transition to nonlinearity and new constraints on biasing

We present two new dynamical tests of the biasing hypothesis. The first is based on the amplitude and the shape of the galaxy-galaxy correlation function, $ξ_g(r)$, where $r$ is the separation of the galaxy pair. The second test uses the mean relative peculiar velocity for galaxy pairs, $\vs(r)$. This quantity is a measure of the rate of growth of clustering and it is related to the two-point correlation function for the matter density fluctuations, $ξ(r)$. Under the assumption that galaxies trace the mass ($ξ_g = ξ$), the expected relative velocity can be calculated directly from the observed galaxy clustering. The above assumption can be tested by confronting the expected $\vs$ with direct measurements from velocity-distance surveys. Both our methods are checked against N-body experiments and then compared with the $ξ_g(r)$ and $\vs$ estimated from the {\sc APM} galaxy survey and the Mark III catalogue, respectively. Our results suggest that cosmological density parameter is low, $Ω_m \approx 0.3$, and that the {\sc APM} galaxies trace the mass at separations $r \ga 5 \Mlu$, where $h$ is the Hubble constant in units of 100 km s$^{-1}$Mpc. The present results agree with earlier studies, based on comparing higher order correlations in the {\sc APM} with weakly non-linear perturbation theory. Both approaches constrain the linear bias factor to be within 20% of unity. If the existence of the feature we identified in the {\sc APM} $ξ_g(r)$ -- the inflection point near $ξ_g = 1$ -- is confirmed by more accurate surveys, we may have discovered gravity's smoking gun: the long awaited ``shoulder'' in $ξ$, generated by gravitational dynamics and predicted by Gott and Rees 25 years ago.

preprint1997arXiv

The Space Density of Galaxy Peaks and the Linear Matter Power Spectrum

One way of recovering information about the initial conditions of the Universe is by measuring features of the cosmological density field which are preserved during gravitational evolution and galaxy formation. In this paper we study the total number density of peaks in a (galaxy) point distribution smoothed with a filter, evaluating its usefulness as a means of inferring the shape of the initial (matter) power spectrum. We find that in numerical simulations which start from Gaussian initial conditions, the peak density follows well that predicted by the theory of Gaussian density fields, even on scales where the clustering is mildly non-linear. For smaller filter scales, $r \simlt 4-6 \hmpc$, we see evidence of merging as the peak density decreases with time. On larger scales, the peak density is independent of time. One might also expect it to be fairly robust with respect to variations in biasing, i.e. the way galaxies trace mass fluctuations. We find that this is the case when we apply various biasing prescriptions to the matter distribution in simulations. If the initial conditions are Gaussian, it is possible to use the peak density measured from the evolved field to reconstruct the shape of the initial power spectrum. We describe a stable method for doing this and apply it to several biased and unbiased non-linear simulations. We are able to recover the slope of the linear matter power spectrum on scales $k \simlt 0.4 \hmpc^{-1}$. The reconstruction has the advantage of being independent of the cosmological parameters ($Ω$, $Λ$, $H_0$) and of the clustering normalisation ($σ_8$). The peak density and reconstructed power spectrum slope therefore promise to be powerful discriminators between popular cosmological scenarios.

preprint1994arXiv

High-order galaxy correlation functions in the APM Galaxy Survey

We estimate J-point galaxy averaged correlation functions $\wbar_J(θ)$ for $J=2,...,9$, in a sample of the APM Galaxy Survey with more than $1.3 \times 10^6$ galaxies and a depth $ \calD \sim 400 \Mpc$. The hierarchical amplitudes $s_J=\wbar_J/\wbar_2^{J-1}$ are roughly constant, up to $J=9$, between $0.5 \Mpc$ and $2 \Mpc$ and decrease slowly for larger scales. At scales larger than $7 \Mpc$ we find strong similarities between the statistical properties of the galaxy fluctuations and the theoretical properties of matter fluctuations evolving under the influence of gravity in an expanding universe on assumption that the initial fluctuations are small and Gaussian. This is most easily explained if at large scales there is no significant biasing between matter and galaxy fluctuations. The comparison of the skewness in the CfA and SSRS catalogues with comparable sub-samples of the APM indicates that the volume of a ``fair sample'' has to be much larger that the one in the combined CfA/SSRS catalogues.

preprint1993arXiv

The Three-point Function as a Probe of Models for Large-scale Structure

We analyze the consequences of models of structure formation for higher-order ($n$-point) galaxy correlation functions in the mildly non-linear regime. Several variations of the standard $Ω=1$ cold dark matter model with scale-invariant primordial perturbations have recently been introduced to obtain more power on large scales, 20 h$^{-1}$ Mpc, e.g., low-matter-density (non-zero cosmological constant) models, `tilted' primordial spectra, and scenarios with a mixture of cold and hot dark matter. They also include models with an effective scale-dependent bias, such as the cooperative galaxy formation scenario of Bower, etal. (1993). We show that higher-order galaxy correlation functions can provide a useful test of such models and can discriminate between models with true large-scale power in the density field and those where the galaxy power arises from scale-dependent bias: a bias with rapid scale-dependence leads to a dramatic decrease of the hierarchical amplitudes at large scales. Current observational constraints on the three-point amplitudes $Q_3$ and $S_3$ can place limits on the bias parameter(s) and appear to disfavor, but not yet rule out, the hypothesis that scale-dependent bias is responsible for the extra power observed on large scales.