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 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
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.
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
We analyse the metallicity histories of ∼4500 galaxies from the GAMA survey at z < 0.06 modelled by the SED-fitting code ProSpect using an evolving metallicity implementation. These ...metallicity histories, in combination with the associated star formation histories, allow us to analyse the inferred gas-phase mass–metallicity relation. Furthermore, we extract the mass–metallicity relation at a sequence of epochs in cosmic history, to track the evolving mass–metallicity relation with time. Through comparison with observations of gas-phase metallicity over a large range of redshifts, we show that, remarkably, our forensic SED analysis has produced an evolving mass–metallicity relationship that is consistent with observations at all epochs. We additionally analyse the three-dimensional mass–metallicity–SFR space, showing that galaxies occupy a clearly defined plane. This plane is shown to be subtly evolving, displaying an increased tilt with time caused by general enrichment, and also the slowing down of star formation with cosmic time. This evolution is most apparent at lookback times greater than 7 Gyr. The trends in metallicity recovered in this work highlight that the evolving metallicity implementation used within the SED-fitting code ProSpect produces reasonable metallicity results over the history of a galaxy. This is expected to provide a significant improvement to the accuracy of the SED-fitting outputs.
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 reinvestigate the correlation between black hole mass and bulge concentration. With an increased galaxy sample (totaling 27) and updated estimates of galaxy distances, black hole masses, and ...Sersic indices n--a measure of concentration--we perform a least-squares regression analysis to obtain a relation suitable for the purpose of predicting black hole masses in other galaxies. In addition to the linear relation, log M sub(bh) = super((7.81 c 0.08)) + super((2.69 c 0.28)) log(n/3) with e sub(intrinsic) = 0.31 super(+) sub(-) super(0) sub(0) super(.) sub(.) super(0) sub(0) super(9) sub(7) dex, we investigated the possibility of a higher order M sub(bh)-n relation, finding the second-order term in the best-fitting quadratic relation to be inconsistent with a value of zero at greater than the 99.99% confidence level. The optimal relation is given by log M sub(bh) = (7.98 c 0.09) + (3.70 c 0.46) log(n/3) - super((3.10 c 0.84)) log(n/3) super(2), with e sub(intrinsic) = 0.18 super(+) sub(-) super(0) sub(0) super(.) sub(.) super(0) sub(0) super(7) sub(6) dex and a total absolute scatter of 0.31 dex. When the quadratic relation is extrapolated, it predicts black holes with masses of 610 super(3) M sub( )in n = 0.5 dwarf elliptical galaxies, compared to 610 super(5) M sub( )from the linear relation, and an upper bound on the largest black hole masses in the local universe equal to 1.2 super(+) sub(-) super(2) sub(0) super(.) sub(.) super(6) sub(4) x 10 super(9) M sub( ). In addition, we show that the nuclear star clusters at the centers of low-luminosity elliptical galaxies follow an extrapolation of the same quadratic relation, strengthening suggestions for a possible evolutionary link between supermassive black holes and nuclear star clusters. Moreover, we speculate that the merger of two such nucleated galaxies, accompanied by the merger and runaway collision of their central star clusters, could result in the late-time formation of some supermassive black holes. Finally, we predict the existence of, and provide equations for, an M sub(bh)-k sub(0) relation, in which k sub(0) is the (extrapolated) central surface brightness of a bulge.
ABSTRACT
We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for ∼7000 galaxies from the GAMA survey at z < 0.06, in order to extract their star formation ...histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies and stellar masses. Our analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago. In other morphological types, the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 11 Gyr ago, whereas the least massive galaxies reached this stage as late as 4 Gyr ago (the well-known effect of ‘galaxy downsizing’). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with independent, higher redshift observations of metal densities in the Universe.