Quantifying cosmic variance Driver, Simon P.; Robotham, Aaron S. G.
Monthly notices of the Royal Astronomical Society,
October 2010, Letnik:
407, Številka:
4
Journal Article
Recenzirano
Odprti dostop
We determine an expression for the cosmic variance of any ‘normal’ galaxy survey based on examination of M*± 1 mag galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) data cube. We ...find that cosmic variance will depend on a number of factors principally: total survey volume, survey aspect ratio and whether the area surveyed is contiguous or comprising independent sightlines. As a rule of thumb cosmic variance falls below 10 per cent once a volume of 107 h−30.7 Mpc3 is surveyed for a single contiguous region with a 1:1 aspect ratio. Cosmic variance will be lower for higher aspect ratios and/or non-contiguous surveys. Extrapolating outside our test region we infer that cosmic variance in the entire SDSS DR7 main survey region is ∼7 per cent to z < 0.1. The equation obtained from the SDSS DR7 region can be generalized to estimate the cosmic variance for any density measurement determined from normal galaxies (e.g. luminosity densities, stellar mass densities and cosmic star formation rates) within the volume range 103–107 h−30.7 Mpc3. We apply our equation to show that two sightlines are required to ensure that cosmic variance is <10 per cent in any ASKAP galaxy survey (divided into Δz∼ 0.1 intervals, i.e. ∼1 Gyr intervals for z < 0.5). Likewise 10 MeerKAT sightlines will be required to meet the same conditions. GAMA, VVDS and zCOSMOS all suffer less than 10 per cent cosmic variance (∼3–8 per cent) in Δz intervals of 0.1, 0.25 and 0.5, respectively. Finally we show that cosmic variance is potentially at the 50–70 per cent level, or greater, in the Hubble Space Telescope (HST) Ultra Deep Field depending on assumptions as to the evolution of clustering. 100 or 10 independent sightlines will be required to reduce cosmic variance to a manageable level (<10 per cent) for HST ACS or HST WFC3 surveys, respectively (in Δz∼ 1 intervals). Cosmic variance is therefore a significant factor in the z > 6 HST studies currently underway.
We present a new, open source, free, semi-analytic model (SAM) of galaxy formation, SHARK, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ...ways of modelling them. We introduce the philosophy behind SHARK and provide an overview of the physical processes included in the model. SHARK is written in C++11 and has been parallelized with OpenMP. In the released version (V1.1), we implement several different models for gas cooling, active galactic nuclei, stellar and photo-ionization feedback, and star formation (SF). We demonstrate the basic performance of SHARK using the Planck Collaboration et al. (2016) cosmology SURFS simulations, by comparing against a large set of observations, including: the stellar mass function (SMF) and stellar-halo mass relation at z = 0-4; the cosmic evolution of the star formation rate density (SFRD), stellar mass, atomic and molecular hydrogen; local gas scaling relations; and structural galaxy properties, finding excellent agreement. Significant improvements over previous SAMs are seen in the mass-size relation for discs/bulges, the gas-stellar mass and stellar mass-metallicity relations. To illustrate the power of SHARK in exploring the systematic effects of the galaxy formation modelling, we quantify how the scatter of the SF main sequence and the gas scaling relations changes with the adopted SF law, and the effect of the starbursts H2 depletion time-scale on the SFRD and Ω _H_2. We compare SHARK with other SAMs and the hydrodynamical simulation EAGLE, and find that SAMs have a much higher halo baryon fractions due to large amounts of intra-halo gas, which in the case of EAGLE is in the intergalactic medium.
Abstract
We introduce ProFound, a source finding and image analysis package. ProFound provides methods to detect sources in noisy images, generate segmentation maps identifying the pixels belonging ...to each source, and measure statistics like flux, size, and ellipticity. These inputs are key requirements of ProFit, our recently released galaxy profiling package, where the design aim is that these two software packages will be used in unison to semi-automatically profile large samples of galaxies. The key novel feature introduced in ProFound is that all photometry is executed on dilated segmentation maps that fully contain the identifiable flux, rather than using more traditional circular or ellipse-based photometry. Also, to be less sensitive to pathological segmentation issues, the de-blending is made across saddle points in flux. We apply ProFound in a number of simulated and real-world cases, and demonstrate that it behaves reasonably given its stated design goals. In particular, it offers good initial parameter estimation for ProFit, and also segmentation maps that follow the sometimes complex geometry of resolved sources, whilst capturing nearly all of the flux. A number of bulge–disc decomposition projects are already making use of the ProFound and ProFit pipeline, and adoption is being encouraged by publicly releasing the software for the open source r data analysis platform under an LGPL-3 license on GitHub (github.com/asgr/ProFound).
We present a simple method of generating approximate synthetic halo catalogues: halogen. This method uses a combination of second-order Lagrangian Perturbation Theory (2lpt) in order to generate the ...large-scale matter distribution, analytical mass functions to generate halo masses, and a single-parameter stochastic model for halo bias to position haloes. halogen represents a simplification of similar recently published methods. Our method is constrained to recover the two-point function at intermediate (10 h
−1 Mpc < r < 50 h
−1 Mpc) scales, which we show is successful to within 2 per cent. Larger scales (∼100 h
−1 Mpc) are reproduced to within 15 per cent. We compare several other statistics (e.g. power spectrum, point distribution function, redshift space distortions) with results from N-body simulations to determine the validity of our method for different purposes. One of the benefits of halogen is its flexibility, and we demonstrate this by showing how it can be adapted to varying cosmologies and simulation specifications. A driving motivation for the development of such approximate schemes is the need to compute covariance matrices and study the systematic errors for large galaxy surveys, which requires thousands of simulated realizations. We discuss the applicability of our method in this context, and conclude that it is well suited to mass production of appropriate halo catalogues. The code is publicly available at https://github.com/savila/halogen.
We present single-Sérsic two-dimensional (2D) model fits to 167 600 galaxies modelled independently in the ugrizYJHK bandpasses using reprocessed Sloan Digital Sky Survey Data Release Seven (SDSS ...DR7) and UKIRT Infrared Deep Sky Survey Large Area Survey imaging data available from the Galaxy And Mass Assembly (GAMA) data base. In order to facilitate this study we developed Structural Investigation of Galaxies via Model Analysis (sigma), an r wrapper around several contemporary astronomy software packages including source extractor, psf extractor and galfit 3. sigma produces realistic 2D model fits to galaxies, employing automatic adaptive background subtraction and empirical point spread function measurements on the fly for each galaxy in GAMA. Using these results, we define a common coverage area across the three GAMA regions containing 138 269 galaxies. We provide Sérsic magnitudes truncated at 10r
e which show good agreement with SDSS Petrosian and GAMA photometry for low Sérsic index systems (n < 4), and much improved photometry for high Sérsic index systems (n > 4), recovering as much as Δm= 0.5 mag in the r band. We employ a K-band Sérsic index/u−r colour relation to delineate the massive (n > ∼2) early-type galaxies (ETGs) from the late-type galaxies (LTGs). The mean Sérsic index of these ETGs shows a smooth variation with wavelength, increasing by 30 per cent from g through K. LTGs exhibit a more extreme change in Sérsic index, increasing by 52 per cent across the same range. In addition, ETGs and LTGs exhibit a 38 and 25 per cent decrease, respectively, in half-light radius from g through K. These trends are shown to arise due to the effects of dust attenuation and stellar population/metallicity gradients within galaxy populations.
We present ProFit, a new code for Bayesian two-dimensional photometric galaxy profile modelling. ProFit consists of a low-level c++ library (libprofit), accessible via a command-line interface and ...documented API, along with high-level r (ProFit) and python (PyProFit) interfaces (available at github.com/ICRAR/libprofit, github.com/ICRAR/ProFit, and github.com/ICRAR/pyprofit, respectively). r ProFit is also available pre-built from cran; however, this version will be slightly behind the latest GitHub version. libprofit offers fast and accurate two-dimensional integration for a useful number of profiles, including Sersic, Core-Sersic, broken-exponential, Ferrer, Moffat, empirical King, point-source, and sky, with a simple mechanism for adding new profiles. We show detailed comparisons between libprofit and galfit. Libprofit is both faster and more accurate than galfit at integrating the ubiquitous Sersic profile for the most common values of the Sersic index n (0.5 < n < 8). The high-level fitting code ProFit is tested on a sample of galaxies with both SDSS and deeper KiDS imaging. We find good agreement in the fit parameters, with larger scatter in best-fitting parameters from fitting images from different sources (SDSS versus KiDS) than from using different codes (ProFit versus galfit). A large suite of Monte Carlo-simulated images are used to assess prospects for automated bulge-disc decomposition with ProFit on SDSS, KiDS, and future LSST imaging. We find that the biggest increases in fit quality come from moving from SDSS- to KiDS-quality data, with less significant gains moving from KiDS to LSST.
We determine the low-redshift field galaxy stellar mass function (GSMF) using an area of 143 deg2 from the first three years of the Galaxy And Mass Assembly (GAMA) survey. The magnitude limits of ...this redshift survey are r < 19.4 mag over two-thirds and 19.8 mag over one-third of the area. The GSMF is determined from a sample of 5210 galaxies using a density-corrected maximum volume method. This efficiently overcomes the issue of fluctuations in the number density versus redshift. With H
0= 70 km s−1 Mpc−1, the GSMF is well described between 108 and 1011.5 M⊙ using a double Schechter function with
,
, α1=−0.35,
and α2=−1.47. This result is more robust to uncertainties in the flow-model corrected redshifts than from the shallower Sloan Digital Sky Survey main sample (r < 17.8 mag). The upturn in the GSMF is also seen directly in the i-band and K-band galaxy luminosity functions. Accurately measuring the GSMF below 108 M⊙ is possible within the GAMA survey volume but as expected requires deeper imaging data to address the contribution from low surface-brightness galaxies.
We report an expanded sample of visual morphological classifications from the Galaxy and Mass Assembly survey phase two, which now includes 7556 objects (previously 3727 in phase one). We define a ...local (z < 0.06) sample and classify galaxies into E, S0-Sa, SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr, and ‘little blue spheroid’ types. Using these updated classifications, we derive stellar mass function fits to individual galaxy populations divided both by morphological class and more general spheroid- or disc-dominated categories with a lower mass limit of log(M
*/M⊙) = 8 (one dex below earlier morphological mass function determinations). We find that all individual morphological classes and the combined spheroid-/bulge-dominated classes are well described by single Schechter stellar mass function forms. We find that the total stellar mass densities for individual galaxy populations and for the entire galaxy population are bounded within our stellar mass limits and derive an estimated total stellar mass density of ρ* = 2.5 × 108 M⊙ Mpc−3 h
0.7, which corresponds to an approximately 4 per cent fraction of baryons found in stars. The mass contributions to this total stellar mass density by galaxies that are dominated by spheroidal components (E and S0-Sa classes) and by disc components (Sab-Scd and Sd-Irr classes) are approximately 70 and 30 per cent, respectively.
ABSTRACT We combine wide and deep galaxy number-count data from the Galaxy And Mass Assembly, COSMOS/G10, Hubble Space Telescope (HST) Early Release Science, HST UVUDF, and various near-, mid-, and ...far-IR data sets from ESO, Spitzer, and Herschel. The combined data range from the far UV (0.15 m) to far-IR (500 m), and in all cases the contribution to the integrated galaxy light (IGL) of successively fainter galaxies converges. Using a simple spline fit, we derive the IGL and the extrapolated IGL in all bands. We argue that undetected low-surface-brightness galaxies and intracluster/group light are modest, and that our extrapolated-IGL measurements are an accurate representation of the extragalactic background light (EBL). Our data agree with most earlier IGL estimates and with direct measurements in the far IR, but disagree strongly with direct estimates in the optical. Close agreement between our results and recent very high-energy experiments (H.E.S.S. and MAGIC) suggests that there may be an additional foreground affecting the direct estimates. The most likely culprit could be the adopted model of zodiacal light. Finally we use a modified version of the two-component model to integrate the EBL and obtain measurements of the cosmic optical background (COB) and cosmic infrared background of nW m−2 sr−1 and nW m−2 sr−1 respectively (48%:52%). Over the next decade, upcoming space missions such as Euclid and the Wide Field Infrared Space Telescope will have the capacity to reduce the COB error to <1%, at which point comparisons to the very high-energy data could have the potential to provide a direct detection and measurement of the reionization field.
ABSTRACT
We combine the shark semi-analytic model of galaxy formation with the prospect software tool for spectral energy distribution (SED) generation to study the multiwavelength emission of ...galaxies from the far-ultraviolet (FUV) to the far-infrared (FIR) at 0 ≤ z ≤ 10. We produce a physical model for the attenuation of galaxies across cosmic time by combining a local Universe empirical relation to compute the dust mass of galaxies from their gas metallicity and mass, attenuation curves derived from radiative transfer calculations of galaxies in the eagle hydrodynamic simulation suite, and the properties of shark galaxies. We are able to produce a wide range of galaxies, from the z = 8 star-forming galaxies with almost no extinction, z = 2 submillimetre galaxies, down to the normal star-forming and red-sequence galaxies at z = 0. Quantitatively, we find that shark reproduces the observed (i) z = 0 FUV-to-FIR, (ii) 0 ≤ z ≤ 3 rest-frame K-band, and (iii) 0 ≤ z ≤ 10 rest-frame FUV luminosity functions, (iv) z ≤ 8 UV slopes, (v) the FUV-to-FIR number counts (including the widely disputed 850 μm), (vi) redshift distribution of bright $850\, \mu$m galaxies, and (vii) the integrated cosmic SED from z = 0 to 1 to an unprecedented level. This is achieved without the need to invoke changes in the stellar initial mass function, dust-to-metal mass ratio, or metal enrichment time-scales. Our model predicts star formation in galaxy discs to dominate in the FUV-to-optical, while bulges dominate at the NIR at all redshifts. The FIR sees a strong evolution in which discs dominate at z ≤ 1 and starbursts (triggered by both galaxy mergers and disc instabilities, in an even mix) dominate at higher redshifts, even out to z = 10.
PDF
