Quantifying cosmic variance Driver, Simon P.; Robotham, Aaron S. G.
Monthly notices of the Royal Astronomical Society,
October 2010, Volume:
407, Issue:
4
Journal Article
Peer reviewed
Open access
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.
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 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.
In this paper, we demonstrate a new method for fitting galaxy profiles which makes use of the full multiwavelength data provided by modern large optical-near-infrared imaging surveys. We present a ...new version of galapagos, which utilizes a recently developed multiwavelength version of galfit, and enables the automated measurement of wavelength-dependent Sérsic profile parameters for very large samples of galaxies. Our new technique is extensively tested to assess the reliability of both pieces of software, galfit and galapagos on both real ugrizY JHK imaging data from the Galaxy And Mass Assembly survey and simulated data made to the same specifications. We find that fitting galaxy light profiles with multiwavelength data increases the stability and accuracy of the measured parameters, and hence produces more complete and meaningful multiwavelength photometry than has been available previously. The improvement is particularly significant for magnitudes in low-S/N bands and for structural parameters like half-light radius r
e and Sérsic index n for which a prior is used by constraining these parameters to a polynomial as a function of wavelength. This allows the fitting routines to push the magnitude of galaxies for which sensible values can be derived to fainter limits. The technique utilizes a smooth transition of galaxy parameters with wavelength, creating more physically meaningful transitions than single-band fitting and allows accurate interpolation between passbands, perfect for derivation of rest-frame values.
Abstract
Our nearest large cosmological neighbour, the Andromeda galaxy (M31), is a dynamical system, and an accurate measurement of its total mass is central to our understanding of its assembly ...history, the life-cycles of its satellite galaxies, and its role in shaping the Local Group environment. Here, we apply a novel approach to determine the dynamical mass of M31 using high-velocity Planetary Nebulae, establishing a hierarchical Bayesian model united with a scheme to capture potential outliers and marginalize over tracers unknown distances. With this, we derive the escape velocity run of M31 as a function of galactocentric distance, with both parametric and non-parametric approaches. We determine the escape velocity of M31 to be 470 ± 40 km s−1 at a galactocentric distance of 15 kpc, and also, derive the total potential of M31, estimating the virial mass and radius of the galaxy to be 0.8 ± 0.1 × 1012 M⊙ and 240 ± 10 kpc, respectively. Our M31 mass is on the low side of the measured range, this supports the lower expected mass of the M31–Milky Way system from the timing and momentum arguments, satisfying the H i constraint on circular velocity between 10 ≲ R/ kpc < 35, and agreeing with the stellar mass Tully–Fisher relation. To place these results in a broader context, we compare them to the key predictions of the ΛCDM cosmological paradigm, including the stellar-mass–halo-mass and the dark matter halo concentration–virial mass correlation, and finding it to be an outlier to this relation.
Abstract
We use the energy-balance code magphys to determine stellar and dust masses, and dust corrected star formation rates for over 200 000 GAMA galaxies, 170 000 G10-COSMOS galaxies, and 200 000 ...3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous data set spanning a broad range in stellar-mass (108–1012 M⊙), dust-mass (106–109 M⊙), and star formation rates (0.01–100 M⊙yr−1), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous data set using consistent mass and star formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our CSFH we precisely reproduce the SMD
with an interstellar medium replenishment factor of 0.50 ± 0.07, consistent with our choice of Chabrier initial mass function plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust-mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065–0.004 units of dust mass is also formed. (2) Over the history of the Universe approximately 90–95 per cent of all dust formed has been destroyed and/or ejected.
PDF
