Does galaxy evolution proceed through the green valley via multiple pathways or as a single population? Motivated by recent results highlighting radically different evolutionary pathways between ...early- and late-type galaxies, we present results from a simple Bayesian approach to this problem wherein we model the star formation history (SFH) of a galaxy with two parameters, t, τ and compare the predicted and observed optical and near-ultraviolet colours. We use a novel method to investigate the morphological differences between the most probable SFHs for both disc-like and smooth-like populations of galaxies, by using a sample of 126 316 galaxies (0.01 < z < 0.25) with probabilistic estimates of morphology from Galaxy Zoo. We find a clear difference between the quenching time-scales preferred by smooth- and disc-like galaxies, with three possible routes through the green valley dominated by smooth- (rapid time-scales, attributed to major mergers), intermediate- (intermediate time-scales, attributed to minor mergers and galaxy interactions) and disc-like (slow time-scales, attributed to secular evolution) galaxies. We hypothesize that morphological changes occur in systems which have undergone quenching with an exponential time-scale τ < 1.5 Gyr, in order for the evolution of galaxies in the green valley to match the ratio of smooth to disc galaxies observed in the red sequence. These rapid time-scales are instrumental in the formation of the red sequence at earlier times; however, we find that galaxies currently passing through the green valley typically do so at intermediate time-scales.
†
When combined with infrared observations with the Spitzer telescope (3 to 160 μm), the Herschel
Space Observatory now fully samples the thermal dust emission up to 500 μm and enables us to better ...estimate the total infrared-submm energy budget (L
TIR) of nearby galaxies. We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved spectral energy distribution (SED) modelling of these objects using the Draine & Li dust models and investigate the influence of the addition of Spectral and Photometric Imaging Receiver (SPIRE) measurements in the estimation of L
TIR. We find that using data up to 250 μm leads to local L
TIR values consistent with those obtained with a complete coverage (up to 500 μm) within ±10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the 70-160 μm band. The 24-70 μm fraction decreases with increasing metallicity. The 160-1100 μmsubmillimetre band can account for up to 25 per cent of the L
TIR in metal-rich galaxies. We investigate the correlation between the total infrared (TIR) surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three Photodetector Array Camera and Spectrometer (PACS) bands can be used as reliable monochromatic estimators of the L
TIR, the 100 μm band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250 μm and L
TIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated L
TIR of nearby galaxies as well as L
TIR of z ∼ 1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.
This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilize the positions, velocities and colours of ...galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealized clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R
200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 1014M⊙, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of 2. Below 1014M⊙, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.
The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We ...study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters with precise spectroscopic redshifts. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples, respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated–true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.
This paper is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilize the positions, velocities and colours of ...galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the rms errors in log M
200c delivered by the different methods (0.18–1.08 dex, i.e. a factor of ∼1.5–12), with abundance-matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or overestimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness- or abundance-matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high-mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilize the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.
Abstract
Taking advantage of the unprecedented combination of sensitivity and angular resolution afforded by the Herschel Space Observatory at far-infrared and submillimetre wavelengths, we aim to ...characterize the physical properties of cold dust within nearby galaxies, as well as the associated uncertainties, namely the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) key programme, we perform two-temperature fits of the Spitzer and Herschel photometric data (from 24 to 500 μm), with a warm and a cold component, both globally and in each resolution element. We compare the results obtained from different analysis strategies. At global scale, we observe a range of values of the modified blackbody fit parameters βc (0.8-2.5) and T
c (19.1-25.1 K). We compute maps of our modelling parameters with βc fixed or treated as a free parameter to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe steeper temperature gradients as a function of radius than when it is allowed to vary. When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities. Gathering the parameters obtained in each resolution element in a T
c-βc diagram underlines an anticorrelation between the two parameters. It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar-ends (coinciding with molecular gas density enhancements usually found in these locations). We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (β < 2.0) lead to significantly lower dust masses compared to the β < 2.0 case, with dust masses lower by up to 50 per cent if βc = 1.5, for instance. The working resolution affects our total dust mass estimates: masses increase from global fits to spatially resolved fits.
Interstellar dust and starlight are modeled for the galaxies of the project "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel." The galaxies were observed by the Infrared Array ...Camera and the Multiband Imaging Photometer for Spitzer on Spitzer Space Telescope, and the Photodetector Array Camera and Spectrometer and the Spectral and Photometric Imaging Receiver on Herschel Space Observatory. With data from 3.6 to 500 m, dust models are strongly constrained. Using a physical dust model, for each pixel in each galaxy we estimate (1) dust surface density, (2) dust mass fraction in polycyclic aromatic hydrocarbons (PAHs), (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in subregions with high starlight intensity. The dust models successfully reproduce the observed global and resolved spectral energy distributions. With the angular resolution of Herschel, we obtain well-resolved maps (available online) for the dust properties. As in previous studies, we find the PAH fraction to be an increasing function of metallicity, with a threshold oxygen abundance Z/Z 0.1, but we find the data to be fitted best with increasing linearly with above a threshold value of 0.15(O/H) . We obtain total dust masses for each galaxy by summing the dust mass over the individual map pixels; these "resolved" dust masses are consistent with the masses inferred from a model fit to the global photometry. The global dust-to-gas ratios obtained from this study are found to correlate with galaxy metallicities. Systems with Z/Z 0.5 have most of their refractory elements locked up in dust, whereas in systems with Z/Z 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that is suppressed in regions with unusually warm dust with . With knowledge of one long-wavelength flux density ratio (e.g., f160/f500), the minimum starlight intensity heating the dust ( ) can be estimated to within ∼50%, despite a variation in of more than two orders of magnitude. For the adopted dust model, dust masses can be estimated to within ∼0.2 dex accuracy using the f160/f500 flux ratio and the integrated dust luminosity, and to ∼0.07 dex accuracy using the 500 m luminosity alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions for estimating starlight heating intensity and dust mass may be useful for studies of high-redshift galaxies.
Using free-free emission measured in the Ka band (26-40 GHz) for 10 star-forming regions in the nearby galaxy NGC 6946, including its starbursting nucleus, we compare a number of star formation rate ...(SFR) diagnostics that are typically considered to be unaffected by interstellar extinction. These diagnostics include non-thermal radio (i.e., 1.4 GHz), total infrared (IR; 8-1000 Delta *mm), and warm dust (i.e., 24 Delta *mm) emission, along with hybrid indicators that attempt to account for obscured and unobscured emission from star-forming regions including H Delta *a + 24 Delta *mm and UV + IR measurements. The assumption is made that the 33 GHz free-free emission provides the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24 Delta *mm, H Delta *a + 24 Delta *mm, and UV + IR SFR calibrations are in good agreement with the 33 GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimates the nuclear SFR by a factor of ~2 relative to the 33 GHz free-free SFR. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the 33 GHz free-free SFRs for the nucleus and underestimate the SFRs from the extranuclear star-forming regions by an average factor of ~2 and ~4-5 before and after subtracting local background emission, respectively. This result likely arises from the cosmic-ray (CR) electrons decaying within the starburst region with negligible escape, whereas the transient nature of star formation in the young extranuclear star-forming complexes allows for CR electrons to diffuse significantly further than dust-heating photons, resulting in an underestimate of the true SFR. Finally, we find that the SFRs estimated using the total 33 GHz flux density appear to agree well with those estimated using free-free emission due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33 GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.
Abstract
With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical ...uncertainties, the control of systematics in cluster masses is crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an ‘unrelaxed’ state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by ∼ 10 per cent at 1014 and ≳ 20 per cent for ≲ 1013.5. The use of cluster samples with different levels of substructure can therefore bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.