Source author record

A. N. Taylor

A. N. Taylor 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

56works
3topics
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

56 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid's unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10's stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10's orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1).

preprint2023arXiv

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Euclid preparation. XVIII. The NISP photometric system

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μ$m range, to a 5$σ$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime.

preprint2022arXiv

Euclid: Constraining ensemble photometric redshift distributions with stacked spectroscopy

The ESA Euclid mission will produce photometric galaxy samples over 15000 square degrees of the sky that will be rich for clustering and weak lensing statistics. The accuracy of the cosmological constraints derived from these measurements will depend on the knowledge of the underlying redshift distributions based on photometric redshift calibrations. A new approach is proposed to use the stacked spectra from Euclid slitless spectroscopy to augment broad-band photometric information to constrain the redshift distribution with spectral energy distribution fitting. The high spectral resolution available in the stacked spectra complements the photometry and helps to break the colour-redshift degeneracy and constrain the redshift distribution of galaxy samples. We modelled the stacked spectra as a linear mixture of spectral templates. The mixture may be inverted to infer the underlying redshift distribution using constrained regression algorithms. We demonstrate the method on simulated Vera C. Rubin Observatory and Euclid mock survey data sets based on the Euclid Flagship mock galaxy catalogue. We assess the accuracy of the reconstruction by considering the inference of the baryon acoustic scale from angular two-point correlation function measurements. We selected mock photometric galaxy samples at redshift z>1 using the self-organising map algorithm. Considering the idealised case without dust attenuation, we find that the redshift distributions of these samples can be recovered with 0.5% accuracy on the baryon acoustic scale. The estimates are not significantly degraded by the spectroscopic measurement noise due to the large sample size. However, the error degrades to 2% when the dust attenuation model is left free. We find that the colour degeneracies introduced by attenuation limit the accuracy considering the wavelength coverage of Euclid near-infrared spectroscopy.

preprint2022arXiv

Euclid: Covariance of weak lensing pseudo-$C_\ell$ estimates. Calculation, comparison to simulations, and dependence on survey geometry

An accurate covariance matrix is essential for obtaining reliable cosmological results when using a Gaussian likelihood. In this paper we study the covariance of pseudo-$C_\ell$ estimates of tomographic cosmic shear power spectra. Using two existing publicly available codes in combination, we calculate the full covariance matrix, including mode-coupling contributions arising from both partial sky coverage and non-linear structure growth. For three different sky masks, we compare the theoretical covariance matrix to that estimated from publicly available N-body weak lensing simulations, finding good agreement. We find that as a more extreme sky cut is applied, a corresponding increase in both Gaussian off-diagonal covariance and non-Gaussian super-sample covariance is observed in both theory and simulations, in accordance with expectations. Studying the different contributions to the covariance in detail, we find that the Gaussian covariance dominates along the main diagonal and the closest off-diagonals, but further away from the main diagonal the super-sample covariance is dominant. Forming mock constraints in parameters describing matter clustering and dark energy, we find that neglecting non-Gaussian contributions to the covariance can lead to underestimating the true size of confidence regions by up to 70 per cent. The dominant non-Gaussian covariance component is the super-sample covariance, but neglecting the smaller connected non-Gaussian covariance can still lead to the underestimation of uncertainties by 10--20 per cent. A real cosmological analysis will require marginalisation over many nuisance parameters, which will decrease the relative importance of all cosmological contributions to the covariance, so these values should be taken as upper limits on the importance of each component.

preprint2022arXiv

Euclid: Forecast constraints on consistency tests of the $Λ$CDM model

The standard cosmological model is based on the fundamental assumptions of a spatially homogeneous and isotropic universe on large scales. An observational detection of a violation of these assumptions at any redshift would immediately indicate the presence of new physics. We quantify the ability of the Euclid mission, together with contemporary surveys, to improve the current sensitivity of null tests of the canonical cosmological constant $Λ$ and the cold dark matter (LCDM) model in the redshift range $0<z<1.8$. We considered both currently available data and simulated Euclid and external data products based on a LCDM fiducial model, an evolving dark energy model assuming the Chevallier-Polarski-Linder (CPL) parameterization or an inhomogeneous Lema\^ıtre-Tolman-Bondi model with a cosmological constant $Λ$ (LLTB), and carried out two separate but complementary analyses: a machine learning reconstruction of the null tests based on genetic algorithms, and a theory-agnostic parametric approach based on Taylor expansion and binning of the data, in order to avoid assumptions about any particular model. We find that in combination with external probes, Euclid can improve current constraints on null tests of the LCDM by approximately a factor of three when using the machine learning approach and by a further factor of two in the case of the parametric approach. However, we also find that in certain cases, the parametric approach may be biased against or missing some features of models far from LCDM. Our analysis highlights the importance of synergies between Euclid and other surveys. These synergies are crucial for providing tighter constraints over an extended redshift range for a plethora of different consistency tests of some of the main assumptions of the current cosmological paradigm.

preprint2022arXiv

KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using eight bands in the multipole range $76<\ell<1500$ for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-$C_{\ell}$ estimates by forward modelling the effects of the mask. Assuming a spatially flat $Λ$CDM cosmology, we constrained the structure growth parameter $S_8 = σ_8(Ω_{\rm m}/0.3)^{1/2} = 0.754_{-0.029}^{+0.027}$. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-$α$ forest and its cross-correlation with quasars, we tightened these constraints to $S_8 = 0.771^{+0.006}_{-0.032}$. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a $\sim 3σ$ tension with early-Universe constraints from cosmic microwave background experiments.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond.

preprint2021arXiv

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2021arXiv

Euclid: Forecast constraints on the cosmic distance duality relation with complementary external probes

In metric theories of gravity with photon number conservation, the luminosity and angular diameter distances are related via the Etherington relation, also known as the distance-duality relation (DDR). A violation of this relation would rule out the standard cosmological paradigm and point at the presence of new physics. We quantify the ability of Euclid, in combination with contemporary surveys, to improve the current constraints on deviations from the DDR in the redshift range $0<z<1.6$. We start by an analysis of the latest available data, improving previously reported constraints by a factor of 2.5. We then present a detailed analysis of simulated Euclid and external data products, using both standard parametric methods (relying on phenomenological descriptions of possible DDR violations) and a machine learning reconstruction using Genetic Algorithms. We find that for parametric methods Euclid can (in combination with external probes) improve current constraints by approximately a factor of six, while for non-parametric methods Euclid can improve current constraints by a factor of three. Our results highlight the importance of surveys like Euclid in accurately testing the pillars of the current cosmological paradigm and constraining physics beyond the standard cosmological model.

preprint2021arXiv

Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids

Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmology. In particular, we study the imprints of dynamic and geometric distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11000 voids and estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on $f/b$ and $D_M H$, where $f$ is the linear growth rate of density fluctuations, $b$ the galaxy bias, $D_M$ the comoving angular diameter distance, and $H$ the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of $f/b$ and 0.5% on $D_M H$ in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter $w$ for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts.

preprint2021arXiv

The Cosmological Trajectories Method: Modelling cosmic structure formation in the non-linear regime

We introduce a novel approach, the Cosmological Trajectories Method (CTM), to model nonlinear structure formation in the Universe by expanding gravitationally-induced particle trajectories around the Zel'dovich approximation. A new Beyond Zel'dovich approximation is presented, which expands the CTM to leading second-order in the gravitational interaction and allows for post-Born gravitational scattering. In the Beyond Zel'dovich approximation we derive the exact expression for the matter clustering power spectrum. This is calculated to leading order and is available in the CTM MODULE. We compare the Beyond Zel'dovich approximation power spectrum and correlation function to other methods including 1-loop Standard Perturbation Theory (SPT), 1-loop Lagrangian Perturbation Theory (LPT) and Convolution Lagrangian Perturbation Theory (CLPT). We find that the Beyond Zel'dovich approximation power spectrum performs well, matching simulations to within $\pm{10}\%$, on mildly non-linear scales, and at redshifts above $z=1$ it outperforms the Zel'dovich approximation. We also find that the Beyond Zel'dovich approximation models the BAO peak in the correlation function at $z=0$ more accurately, to within $\pm{5}\%$ of simulations, than the Zel'dovich approximation, SPT 1-loop and CLPT.

preprint2020arXiv

Euclid preparation: VIII. The Complete Calibration of the Colour-Redshift Relation survey: VLT/KMOS observations and data release

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

preprint2020arXiv

Euclid: Nonparametric point spread function field recovery through interpolation on a graph Laplacian

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy. Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario. Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations. Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

preprint2020arXiv

Euclid: The reduced shear approximation and magnification bias for Stage IV cosmic shear experiments

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

preprint2020arXiv

Euclid: The selection of quiescent and star-forming galaxies using observed colours

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

preprint2016arXiv

Cosmic Shear Bias and Calibration in Cosmic Shear Studies

With the advent of large-scale weak lensing surveys there is a need to understand how realistic, scale-dependent systematics bias cosmic shear and dark energy measurements, and how they can be removed. Here we describe how spatial variations in the amplitude and orientation of realistic image distortions convolve with the measured shear field, mixing the even-parity convergence and odd-parity modes, and bias the shear power spectrum. Many of these biases can be removed by calibration to external data, the survey itself, or by modelling in simulations. The uncertainty in the calibration must be marginalised over and we calculate how this propagates into parameter estimation, degrading the dark energy Figure-of-Merit. We find that noise-like biases affect dark energy measurements the most, while spikes in the bias power have the least impact, reflecting their correlation with the effect of cosmological parameters. We argue that in order to remove systematic biases in cosmic shear surveys and maintain statistical power effort should be put into improving the accuracy of the bias calibration rather than minimising the size of the bias. In general, this appears to be a weaker condition for bias removal. We also investigate how to minimise the size of the calibration set for a fixed reduction in the Figure-of-Merit. These results can be used to model the effect of biases and calibration on a cosmic shear survey accurately, assess their impact on the measurement of modified gravity and dark energy models, and to optimise surveys and calibration requirements.

preprint2015arXiv

On Scale-Dependent Cosmic Shear Systematic Effects

In this paper we investigate the impact that realistic scale-dependence systematic effects may have on cosmic shear tomography. We model spatially varying residual ellipticity and size variations in weak lensing measurements and propagate these through to predicted changes in the uncertainty and bias of cosmological parameters. We show that the survey strategy - whether it is regular or randomised - is an important factor in determining the impact of a systematic effect: a purely randomised survey strategy produces the smallest biases, at the expense of larger parameter uncertainties, and a very regularised survey strategy produces large biases, but unaffected uncertainties. However, by removing, or modelling, the affected scales (l-modes) in the regular cases the biases are reduced to negligible levels. We find that the integral of the systematic power spectrum is not a good metric for dark energy performance, and we advocate that systematic effects should be modelled accurately in real space, where they enter the measurement process, and their effect subsequently propagated into power spectrum contributions.

preprint2015arXiv

Path-integral Evidence

Here we present a Bayesian formalism for the goodness-of-fit that is the evidence for a fixed functional form over the evidence for all functions that are a general perturbation about this form. This is done under the assumption that the statistical properties of the data can be modelled by a multivariate Gaussian distribution. We use this to show how one can optimise an experiment to find evidence for a fixed function over perturbations about this function. We apply this formalism to an illustrative problem of measuring perturbations in the dark energy equation of state about a cosmological constant.

preprint2013arXiv

Forecasts of non-Gaussian parameter spaces using Box-Cox transformations

Forecasts of statistical constraints on model parameters using the Fisher matrix abound in many fields of astrophysics. The Fisher matrix formalism involves the assumption of Gaussianity in parameter space and hence fails to predict complex features of posterior probability distributions. Combining the standard Fisher matrix with Box-Cox transformations, we propose a novel method that accurately predicts arbitrary posterior shapes. The Box-Cox transformations are applied to parameter space to render it approximately multivariate Gaussian, performing the Fisher matrix calculation on the transformed parameters. We demonstrate that, after the Box-Cox parameters have been determined from an initial likelihood evaluation, the method correctly predicts changes in the posterior when varying various parameters of the experimental setup and the data analysis, with marginally higher computational cost than a standard Fisher matrix calculation. We apply the Box-Cox-Fisher formalism to forecast cosmological parameter constraints by future weak gravitational lensing surveys. The characteristic non-linear degeneracy between matter density parameter and normalisation of matter density fluctuations is reproduced for several cases, and the capabilities of breaking this degeneracy by weak lensing three-point statistics is investigated. Possible applications of Box-Cox transformations of posterior distributions are discussed, including the prospects for performing statistical data analysis steps in the transformed Gaussianised parameter space.

preprint2013arXiv

Image Analysis for Cosmology: Results from the GREAT10 Galaxy Challenge

In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.

preprint2013arXiv

Probing the accelerating Universe with radio weak lensing in the JVLA Sky Survey

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

preprint2011arXiv

Cosmological information in Gaussianised weak lensing signals

We investigate the information on cosmology contained in Gaussianised weak gravitational lensing convergence fields. Employing Box-Cox transformations to determine optimal transformations to Gaussianity, we develop analytical models for the transformed power spectrum, including effects of noise and smoothing. We find that optimised Box-Cox transformations perform substantially better than an offset logarithmic transformation in Gaussianising the convergence, but both yield very similar results for the signal-to-noise and parameter constraints. None of the transformations is capable of eliminating correlations of the power spectra between different angular frequencies, which we demonstrate to have a significant impact on the errors on cosmology. Analytic models of the Gaussianised power spectrum yield good fits to the simulations and produce unbiased parameter estimates in the majority of cases, where the exceptions can be traced back to the limitations in modelling the higher-order correlations of the original convergence. In the idealistic case, without galaxy shape noise, we find an increase in cumulative signal-to-noise by a factor of 2.6 for angular frequencies up to 1500, and a decrease in the area of the confidence region in the Omega_m-sigma_8 plane by a factor of 4.4 in terms of q-values for the best-performing transformation. When adding a realistic level of shape noise, all transformations perform poorly with little decorrelation of angular frequencies, a maximum increase in signal-to-noise of 34%, and even marginally degraded errors on cosmological parameters. We argue that, to find Gaussianising transformations of practical use, one will need to go beyond transformations of the one-point distribution of the convergence, extend the analysis deeper into the non-linear regime, and resort to an exploration of parameter space via simulations. (abridged)

preprint2011arXiv

Simulating the Effect of Non-Linear Mode-Coupling in Cosmological Parameter Estimation

Fisher Information Matrix methods are commonly used in cosmology to estimate the accuracy that cosmological parameters can be measured with a given experiment, and to optimise the design of experiments. However, the standard approach usually assumes both data and parameter estimates are Gaussian-distributed. Further, for survey forecasts and optimisation it is usually assumed the power-spectra covariance matrix is diagonal in Fourier-space. But in the low-redshift Universe, non-linear mode-coupling will tend to correlate small-scale power, moving information from lower to higher-order moments of the field. This movement of information will change the predictions of cosmological parameter accuracy. In this paper we quantify this loss of information by comparing naive Gaussian Fisher matrix forecasts with a Maximum Likelihood parameter estimation analysis of a suite of mock weak lensing catalogues derived from N-body simulations, based on the SUNGLASS pipeline, for a 2-D and tomographic shear analysis of a Euclid-like survey. In both cases we find that the 68% confidence area of the Omega_m-sigma_8 plane increases by a factor 5. However, the marginal errors increase by just 20 to 40%. We propose a new method to model the effects of nonlinear shear-power mode-coupling in the Fisher Matrix by approximating the shear-power distribution as a multivariate Gaussian with a covariance matrix derived from the mock weak lensing survey. We find that this approximation can reproduce the 68% confidence regions of the full Maximum Likelihood analysis in the Omega_m-sigma_8 plane to high accuracy for both 2-D and tomographic weak lensing surveys. Finally, we perform a multi-parameter analysis of Omega_m, sigma_8, h, n_s, w_0 and w_a to compare the Gaussian and non-linear mode-coupled Fisher matrix contours. (Abridged)

preprint2011arXiv

Spatial matter density mapping of the STAGES Abell A901/2 supercluster field with 3D lensing

We present weak lensing data from the HST/STAGES survey to study the three-dimensional spatial distribution of matter and galaxies in the Abell 901/902 supercluster complex. Our method improves over the existing 3D lensing mapping techniques by calibrating and removing redshift bias and accounting for the effects of the radial elongation of 3D structures. We also include the first detailed noise analysis of a 3D lensing map, showing that even with deep HST quality data, only the most massive structures, for example M200>~10^15 Msun/h at z~0.8, can be resolved in 3D with any reasonable redshift accuracy (Δz~0.15). We compare the lensing map to the stellar mass distribution and find luminous counterparts for all mass peaks detected with a peak significance >3σ. We see structures in and behind the z=0.165 foreground supercluster, finding structure directly behind the A901b cluster at z~0.6 and also behind the SW group at z~0.7. This 3D structure viewed in projection has no significant impact on recent mass estimates of A901b or the SW group components SWa and SWb.

preprint2010arXiv

Analytic Methods for Cosmological Likelihoods

We present general, analytic methods for Cosmological likelihood analysis and solve the "many-parameters" problem in Cosmology. Maxima are found by Newton's Method, while marginalization over nuisance parameters, and parameter errors and covariances are estimated by analytic marginalization of an arbitrary likelihood function with flat or Gaussian priors. We show that information about remaining parameters is preserved by marginalization. Marginalizing over all parameters, we find an analytic expression for the Bayesian evidence for model selection. We apply these methods to data described by Gaussian likelihoods with parameters in the mean and covariance. This method can speed up conventional likelihood analysis by orders of magnitude when combined with Monte-Carlo Markov Chain methods, while Bayesian model selection becomes effectively instantaneous.

preprint2010arXiv

Characterization of the Millimeter-Wave Polarization of Centaurus A with QUaD

Centaurus (Cen) A represents one of the best candidates for an isolated, compact, highly polarized source that is bright at typical cosmic microwave background (CMB) experiment frequencies. We present measurements of the 4 degree by 2 degree region centered on Cen A with QUaD, a CMB polarimeter whose absolute polarization angle is known to 0.5 degrees. Simulations are performed to assess the effect of misestimation of the instrumental parameters on the final measurement and systematic errors due to the field's background structure and temporal variability from Cen A's nuclear region are determined. The total (Q, U) of the inner lobe region is (1.00 +/- 0.07 (stat.) +/- 0.04 (sys.), -1.72 +/- 0.06 +/- 0.05) Jy at 100 GHz and (0.80 +/- 0.06 +/- 0.06, -1.40 +/- 0.07 +/- 0.08) Jy at 150 GHz, leading to polarization angles and total errors of -30.0 +/- 1.1 degrees and -29.1 +/- 1.7 degrees. These measurements will allow the use of Cen A as a polarized calibration source for future millimeter experiments.

preprint2010arXiv

On Mitigation of the Uncertainty in Nonlinear Matter Clustering for Cosmic Shear Tomography

We present a new method that deals with the uncertainty in matter-clustering in cosmic shear power spectrum analysis that arises mainly due to poorly understood nonlinear baryonic processes on small-scales. We show that the majority of information about dark energy physics contained in the shear power comes from these small-scales; removing these nonlinear scales from a cosmic shear analysis results in a 50% cut in the accuracy of measurements of dark energy parameters, marginalizing over all other parameters. In this paper we propose a method to recover the information on small-scales by allowing cosmic shear surveys to measure the nonlinear matter power spectrum themselves and marginalize over all possible power spectra using path integrals. Information is still recoverable in these nonlinear regimes from the geometric part of weak lensing. In this self-calibration regime we find we recover 90% of the information on dark energy. Including an informative prior, we find the nonlinear matter power spectrum needs to be accurately known to 1% down to k=50 h/Mpc to recover 99% of the dark energy information. This presents a significant theoretical challenge to understand baryonic effects on the scale of galaxy haloes. However self-calibration from weak lensing may also provide observational input to help constrain baryon physics.

preprint2010arXiv

Parameter Estimation from Improved Measurements of the CMB from QUaD

We evaluate the contribution of cosmic microwave background (CMB) polarization spectra to cosmological parameter constraints. We produce cosmological parameters using high-quality CMB polarization data from the ground-based QUaD experiment and demonstrate for the majority of parameters that there is significant improvement on the constraints obtained from satellite CMB polarization data. We split a multi-experiment CMB dataset into temperature and polarization subsets and show that the best-fit confidence regions for the LCDM 6-parameter cosmological model are consistent with each other, and that polarization data reduces the confidence regions on all parameters. We provide the best limits on parameters from QUaD EE/BB polarization data and we find best-fit parameters from the multi-experiment CMB dataset using the optimal pivot scale of k_p=0.013 Mpc-1 to be {omch2, ombh2, H_0, A_s, n_s, tau}= {0.113, 0.0224, 70.6, 2.29 times 10^-9, 0.960, 0.086}.

preprint2010arXiv

Path Integral Marginalization for Cosmology: Scale Dependent Galaxy Bias & Intrinsic Alignments

We present a path-integral likelihood formalism that extends parameterized likelihood analyses to include continuous functions. The method finds the maximum likelihood point in function-space, and marginalizes over all possible functions, under the assumption of a Gaussian-distributed function-space. We apply our method to the problem of removing unknown systematic functions in two topical problems for dark energy research : scale-dependent galaxy bias in redshift surveys; and galaxy intrinsic alignments in cosmic shear surveys. We find that scale-dependent galaxy bias will degrade information on cosmological parameters unless the fractional variance in the bias function is known to 10%. Measuring and removing intrinsic alignments from cosmic shear surveys with a flat-prior can reduce the dark energy Figure-of-Merit by 20%, however provided that the scale and redshift-dependence is known to better than 10% with a Gaussian-prior, the dark energy Figure-of-Merit can be enhanced by a factor of two with no extra assumptions.

preprint2010arXiv

SUNGLASS: A new weak lensing simulation pipeline

A new cosmic shear analysis pipeline SUNGLASS (Simulated UNiverses for Gravitational Lensing Analysis and Shear Surveys) is introduced. SUNGLASS is a pipeline that rapidly generates simulated universes for weak lensing and cosmic shear analysis. The pipeline forms suites of cosmological N-body simulations and performs tomographic cosmic shear analysis using line-of-sight integration through these simulations while saving the particle lightcone information. Galaxy shear and convergence catalogues with realistic 3D galaxy redshift distributions are produced for the purposes of testing weak lensing analysis techniques and generating covariance matrices for data analysis and cosmological parameter estimation. We present a suite of fast medium resolution simulations with shear and convergence maps for a generic 100 square degree survey out to a redshift of z = 1.5, with angular power spectra agreeing with the theory to better than a few percent accuracy up to l = 10^3 for all source redshifts up to z = 1.5 and wavenumbers up to l = 2000 for the source redshifts z > 1.1. At higher wavenumbers, there is a failure of the theoretical lensing power spectrum reflecting the known discrepancy of the Smith et al. (2003) fitting formula at high physical wavenumbers. A two-parameter Gaussian likelihood analysis of sigma_8 and Omega_m is also performed on the suite of simulations, demonstrating that the cosmological parameters are recovered from the simulations and the covariance matrices are stable for data analysis. We find no significant bias in the parameter estimation at the level of ~ 0.02. The SUNGLASS pipeline should be an invaluable tool in weak lensing analysis.

preprint2010arXiv

The QUaD Galactic Plane Survey 1: Maps And Analysis of Diffuse Emission

We present a survey of ~800 square degrees of the galactic plane observed with the QUaD telescope. The primary product of the survey are maps of Stokes I, Q and U parameters at 100 and 150 GHz, with spatial resolution 5 and 3.5 arcminutes respectively. Two regions are covered, spanning approximately 245-295 and 315-5 degrees in galactic longitude l, and -4<b<+4 degrees in galactic latitude b. At 0.02 degree square pixel size, the median sensitivity is 74 and 107 kJy/sr at 100 GHz and 150 GHz respectively in I, and 98 and 120 kJy/sr for Q and U. In total intensity, we find an average spectral index of 2.35+/-0.01 (stat) +/- 0.02 (sys) for |b|<1 degree, indicative of emission components other than thermal dust. A comparison to published dust, synchrotron and free-free models implies an excess of emission in the 100 GHz QUaD band, while better agreement is found at 150 GHz. A smaller excess is observed when comparing QUaD 100 GHz data to WMAP 5-year W band; in this case the excess is likely due to the wider bandwidth of QUaD. Combining the QUaD and WMAP data, a two-component spectral fit to the inner galactic plane (|b|<1 degree) yields mean spectral indices of -0.32 +/- 0.03 and 2.84 +/- 0.03; the former is interpreted as a combination of the spectral indices of synchrotron, free-free and dust, while the second is attributed largely to the thermal dust continuum. In the same galactic latitude range, the polarization data show a high degree of alignment perpendicular to the expected galactic magnetic field direction, and exhibit mean polarization fraction 1.38+/-0.08 (stat) +/-0.1 (sys) % at 100 GHz and 1.70 +/- 0.06 (stat) +/- 0.1 (sys) % at 150GHz. We find agreement in polarization fraction between QUaD 100 GHz and WMAP W band, the latter giving 1.1 +/- 0.4 %.

preprint2010arXiv

The QUaD Galactic Plane Survey II: A Compact Source Catalog

We present a catalog of compact sources derived from the QUaD Galactic Plane Survey. The survey covers ~800 square degrees of the inner galaxy (|b|<4 degrees) in Stokes I, Q and U parameters at 100 and 150 GHz, with angular resolution 5 and 3.5 arcminutes. 505 unique sources are identified in I, of which 239 are spatially matched between frequency bands, with 50 (216) detected at 100 (150) GHz alone; 182 sources are identified as ultracompact HII (UCHII) regions. Approximating the distribution of total intensity source fluxes as a power-law, we find a slope of $γ_{S,100}=-1.8\pm0.4$ at 100 GHz, and $γ_{S,150}=-2.2\pm0.4$ at 150 GHz. Similarly, the power-law index of the source two-point angular correlation function is $γ_{θ,100}=-1.21\pm0.04$ and $γ_{θ,150}=-1.25\pm0.04$. The total intensity spectral index distribution peaks at $α_{I}\sim0.25$, indicating that dust emission is not the only source of radiation produced by these objects between 100 and 150 GHz; free-free radiation is likely significant in the 100 GHz band. Four sources are detected in polarized intensity P, of which three have matching counterparts in I. Three of the polarized sources lie close to the galactic center, Sagittarius A*, Sagittarius B2 and the Galactic Radio Arc, while the fourth is RCW 49, a bright HII region. An extended polarized source, undetected by the source extraction algorithm on account of its $\sim0.5^{\circ}$ size, is identified visually, and is an isolated example of large-scale polarized emission oriented distinctly from the bulk galactic dust polarization.

preprint2009arXiv

Improved measurements of the temperature and polarization of the CMB from QUaD

We present an improved analysis of the final dataset from the QUaD experiment. Using an improved technique to remove ground contamination, we double the effective sky area and hence increase the precision of our CMB power spectrum measurements by ~30% versus that previously reported. In addition, we have improved our modeling of the instrument beams and have reduced our absolute calibration uncertainty from 5% to 3.5% in temperature. The robustness of our results is confirmed through extensive jackknife tests and by way of the agreement we find between our two fully independent analysis pipelines. For the standard 6-parameter LCDM model, the addition of QUaD data marginally improves the constraints on a number of cosmological parameters over those obtained from the WMAP experiment alone. The impact of QUaD data is significantly greater for a model extended to include either a running in the scalar spectral index, or a possible tensor component, or both. Adding both the QUaD data and the results from the ACBAR experiment, the uncertainty in the spectral index running is reduced by ~25% compared to WMAP alone, while the upper limit on the tensor-to-scalar ratio is reduced from r < 0.48 to r < 0.33 (95% c.l). This is the strongest limit on tensors to date from the CMB alone. We also use our polarization measurements to place constraints on parity violating interactions to the surface of last scattering, constraining the energy scale of Lorentz violating interactions to < 1.5 x 10^{-43} GeV (68% c.l.). Finally, we place a robust upper limit on the strength of the lensing B-mode signal. Assuming a single flat band power between l = 200 and l = 2000, we constrain the amplitude of B-modes to be < 0.57 micro-K^2 (95% c.l.).

preprint2009arXiv

Interacting Galaxies in the A901/902 Supercluster with STAGES

We present a study of galaxy mergers and the influence of environment in the Abell 901/902 supercluster at z~0.165. We use HST ACS F606W data from the STAGES survey, COMBO-17, Spitzer 24um, and XMM-Newton X-ray data. Our analysis utilizes both a visual classification system, and quantitative CAS parameters to identify systems which show evidence of a recent or ongoing merger of mass ratio >1/10. Our results are: (1) After visual classification and minimizing the contamination from false projection pairs, we find that the merger fraction f_merge is 0.023+/-0.007. The estimated fractions of likely major mergers, likely minor mergers, and ambiguous cases are 0.01+/-0.004, 0.006+/-0.003, and 0.007+/-0.003, respectively. (2) The mergers lie outside the cluster core of radius R < 0.25 Mpc: the lack of mergers in the core is likely due to the large galaxy velocity dispersion in the core. Mergers populate the region (0.25 Mpc < R <= 2 Mpc) between the core and outskirt. In this region, the estimated frequency of mergers is similar to those seen at typical group overdensities. This suggests ongoing growth of the clusters via accretion of group and field galaxies. (3) We compare our observed merger fraction with those reported in other clusters and groups out to z~0.4. Existing data points on the merger fraction for L<= L* galaxies in clusters allow for a range of evolutionary scenarios. (4) The fraction of mergers, which lie on the blue cloud is 80%+/-18% versus 34%+/-7% for non-interacting galaxies, implying that interacting galaxies are preferentially blue. (5) The average SFR, based on UV or UV+IR data, is enhanced by a factor of ~1.5 to 2 in mergers compared to non-interacting galaxies. However, mergers in the clusters contribute only a small fraction (between 10% and 15%) of the total SFR density.(Abridged)

preprint2009arXiv

QUaD: A High-Resolution Cosmic Microwave Background Polarimeter

We describe the QUaD experiment, a millimeter-wavelength polarimeter designed to observe the Cosmic Microwave Background (CMB) from a site at the South Pole. The experiment comprises a 2.64 m Cassegrain telescope equipped with a cryogenically cooled receiver containing an array of 62 polarization-sensitive bolometers. The focal plane contains pixels at two different frequency bands, 100 GHz and 150 GHz, with angular resolutions of 5 arcmin and 3.5 arcmin, respectively. The high angular resolution allows observation of CMB temperature and polarization anisotropies over a wide range of scales. The instrument commenced operation in early 2005 and collected science data during three successive Austral winter seasons of observation.

preprint2009arXiv

Second and third season QUaD CMB temperature and polarization power spectra

We report results from the second and third seasons of observation with the QUaD experiment. Angular power spectra of the Cosmic Microwave Background are derived for both temperature and polarization at both 100 GHz and 150 GHz, and as cross frequency spectra. All spectra are subjected to an extensive set of jackknife tests to probe for possible systematic contamination. For the implemented data cuts and processing technique such contamination is undetectable. We analyze the difference map formed between the 100 and 150 GHz bands and find no evidence of foreground contamination in polarization. The spectra are then combined to form a single set of results which are shown to be consistent with the prevailing LCDM model. The sensitivity of the polarization results is considerably better than that of any previous experiment -- for the first time multiple acoustic peaks are detected in the E-mode power spectrum at high significance.

preprint2008arXiv

Galaxy Evolution in Dense Environments: Properties of Interacting Galaxies in the A901/02 Supercluster

We present a study of galaxies in the STAGES survey of the Abell 901/902 supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps. We characterize galaxies with strong, externally-triggered morphological distortions and normal, relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We report here our preliminary results.

preprint2008arXiv

Morphological Transformations of Galaxies in the A901/02 Supercluster from STAGES

We present a study of galaxies in the Abell 901/902 Supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps, as part of the STAGES survey. We characterize galaxies with strong externally-triggered morphological distortions and normal relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We revisit the morphology density relation, which postulates a higher fraction of early type galaxies in dense environments, by considering separately galaxies with a low bulge-to-disk (B/D) ratio and a low gas content as these two parameters may not be correlated in clusters. We report here on our preliminary analysis.

preprint2008arXiv

Systematic effects on dark energy from 3D weak shear

We present an investigation into the potential effect of systematics inherent in multi-band wide field surveys on the dark energy equation of state determination for two 3D weak lensing methods. The weak lensing methods are a geometric shear-ratio method and 3D cosmic shear. The analysis here uses an extension of the Fisher matrix framework to jointly include photometric redshift systematics, shear distortion systematics and intrinsic alignments. We present results for DUNE and Pan-STARRS surveys. We show that assuming systematic parameters are fixed, but possibly biased, results in potentially large biases in dark energy parameters. We quantify any potential bias by defining a Bias Figure of Merit. We also show the effect on the dark energy Figure of Merit of marginalising over each systematic parameter individually. We find that the largest effect on the Figure of Merit comes from uncertainty in the photometric redshift systematic parameters. These can reduce the Figure of Merit by up to a factor of 2 to 4 in both 3D weak lensing methods, if no informative prior on the systematic parameters is applied. Shear distortion systematics have a smaller overall effect. Intrinsic alignment effects can reduce the Figure of Merit by up to a further factor of 2. This, however, is a worst case scenario. By including prior information on systematic parameters the Figure of Merit can be recovered to a large extent. We conclude that, as a rule of thumb, given a realistic current understanding of intrinsic alignments and photometric redshifts, then including all three primary systematic effects reduces the Figure of Merit by at most a factor of 2, but that in reality this factor should be much less. [abridged]

preprint1997arXiv

Design and analysis of redshift surveys

In this paper we consider methods of analysis and optimal design of redshift surveys. In the first part, we develop a formalism for analysing galaxy redshift surveys which are essentially two-dimensional, such as thin declination slices. The formalism is a power spectrum method, using spherical coordinates, allowing the distorting effects of galaxy peculiar velocities to be calculated to linear order on the assumption of statistical isotropy but without further approximation. In this paper, we calculate the measured two-dimensional power for a constant declination strip, widely used in redshift surveys. We present a likelihood method for estimating the three-dimensional real-space power spectrum and the redshift distortion simultaneously, and show that for thin surveys of reasonable depth, the large-scale 3D power cannot be measured with high accuracy. The redshift distortion may be estimated successfully, and with higher accuracy if the 3D power spectrum can be measured independently, for example from a large-scale sky-projected catalogue. In the second part, we show how a 3D survey design can be optimized to measure the power spectrum, considering whether areal coverage is more important than depth, and whether the survey should be sampled sparsely or not. We show quite generally that width is better than depth, and show how the optimal sparse-sampling fraction, f, depends on the power, P, to be measured. For a Schechter luminosity function, a simple optimization fP \simeq 500 h^{-3} Mpc^3 is found.

preprint1996arXiv

Nonlinear Cosmological Power Spectra in Real and Redshift--Space

We present an expression for the nonlinear evolution of the cosmological power spectrum based on following Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self- gravitating fields. We have extended this analysis into redshift-space and found a solution for the nonlinear, anisotropic redshift-space power spectrum in the limit of plane--parallel redshift distortions. The quadrupole-to- monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, $β$. The point of zero-crossing of the quadrupole, $k_0$, is found to obey a scaling relation. This model is found to be in agreement with $N$-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole--monopole ratio found in the merged QDOT+1.2 Jy IRAS redshift survey. We have estimated that the distortion parameter is constrained to be $β>0.5$ at the $95 \%$ level. The local primordial spectral slope is not well constrained, but analysis suggests $n \approx -2$ in the translinear regime. The zero-crossing scale of the quadrupole is $k_0=0.5 \pm 0.1 h/Mpc$ and from this we infer the amplitude of clustering is $σ_8=0.7 \pm 0.05$. We suggest that the success of this model is due to nonlinear redshift--space effects arising from infall onto caustics and is not dominated by virialised cluster cores.

preprint1994arXiv

A Spherical Harmonic Analysis of Redshift Space

We re-examine the effects of redshift space distortion in all-sky galaxy redshift surveys in the formalism of spherical harmonics. Within this framework we show how one can treat both the large-scale linear effects, and the small-scale nonlinear clustering, exactly to first order. The method also allows in principle a determination of the power spectrum of perturbations, requiring no assumptions beyond that of linear theory. The method therefore offers significant advantages over Fourier techniques when dealing with all-sky surveys. We apply our likelihood analysis to both simulated data, and real data, using the IRAS 1.2-Jy galaxy catalogue, for which we find a maximum likelihood $β\simeq 1.1\pm 0.3$, and a real-space fluctuation amplitude corresponding to $σ_{8,{\rm IRAS}} = 0.68\pm 0.05$.

preprint1993arXiv

Reconstruction Analysis I.: Redshift--space deprojection in the quasi--nonlinear regime

We present a new method for extracting the true 3-d velocity and density fields from the nonlinear redshift--space projected density field. The method is based on the nonlinear, nonlocal transformation of the density field. We assume a curl--free velocity field, although in the formulation presented here this assumption can be dropped. Toy models with special symmetry are used to test the intrinsic accuracy of the method compared to linear theory. Finding significant improvement we then test the method in more detail on N--body simulations, considering the problems associated with sampling the density field using sparse, magnitude--limited galaxy redshift surveys. A method of estimating the magnitude of shot--noise effects on the deprojected velocity field is derived. We find there is a good correlation between the true and deprojected density and velocity fields, even in the mildly nonlinear regime, provided that on large scales galaxies trace the mass. An initial application is made to the IRAS--QDOT all--sky redshift survey.