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.
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.
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 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 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 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.
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.
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 present a simple model for galaxy attenuation by fitting skirt radiative transfer calculations for ∼10 000 eagle galaxies at redshifts z = 2 − 0. Our model adapts the two-component screen ...model of Charlot & Fall, parametrizing the optical depth and slope of the interstellar medium screen using the average dust surface density, Σdust. We recover relatively tight relations between these parameters for the eagle sample, but also provide the scatter in these parameters owing to the morphological variation and orientation of galaxies. We also find that these relations are nearly independent of redshift in the eagle model. By pairing our model with an empirical prescription for birth clouds below the resolution scale of the simulation, we reproduce the observed relation between attenuation slope and optical depth for the first time in a cosmological simulation. We demonstrate that this result is remarkably independent of the attenuation properties assumed for the birth cloud screen, merely requiring a boosted attenuation for infant stars. We present this model with a view to interpreting observations, as well as processing semi-analytic models and other hydrodynamic simulations.