Abstract
We report results from a large molecular line survey of luminous infrared galaxies (LIRGs; L
IR ≳1011L⊙) in the local Universe (z ≤ 0.1), conducted during the last decade with the James ...Clerk Maxwell Telescope and the IRAM 30-m telescope. This work presents the CO and 13CO line data for 36 galaxies, further augmented by multi-J total CO line luminosities available for other infrared (IR) bright galaxies from the literature. This yields a combined sample of N = 70 galaxies with the star formation (SF) powered fraction of their IR luminosities spanning and a wide range of morphologies. Simple comparisons of their available CO spectral line energy distributions (SLEDs) with local ones, as well as radiative transfer models, discern a surprisingly wide range of average interstellar medium (ISM) conditions, with most of the surprises found in the high-excitation regime. These take the form of global CO SLEDs dominated by a very warm (T
kin ≳100 K) and dense (n ≥ 104 cm−3) gas phase, involving galaxy-sized (∼(few) × 109 M⊙) gas mass reservoirs under conditions that are typically found only for ∼(1-3) per cent of mass per typical SF molecular cloud in the Galaxy. Furthermore, some of the highest excitation CO SLEDs are found in ultraluminous infrared galaxies (ULIRGs; L
IR ≥ 1012 L⊙) and surpass even those found solely in compact SF-powered hot spots in Galactic molecular clouds. Strong supersonic turbulence and high cosmic ray energy densities rather than far-ultraviolet/optical photons or supernova remnant induced shocks from individual SF sites can globally warm the large amounts of dense gas found in these merger-driven starbursts and easily power their extraordinary CO line excitation. This exciting possibility can now be systematically investigated with Herschel and the Atacama Large Milimeter Array (ALMA). As expected for an IR-selected (and thus SF rate selected) galaxy sample, only few 'cold' CO SLEDs are found, and for fewer still a cold low/moderate-density and gravitationally bound state (i.e. Galactic type) emerges as the most likely one. The rest remain compatible with a warm and gravitationally unbound low-density phase often found in ULIRGs. Such degeneracies, prominent when only the low-J SLED segment (J = 1-0, 2-1 and 3-2) is available, advise against using its CO line ratios and the so-called X
co = M(H2)/L
co(1-0) factor as SF mode indicators, a practice that may have led to the misclassification of the ISM environments of IR-selected gas-rich discs in the distant Universe. Finally, we expect that the wide range of ISM conditions found among LIRGs will strongly impact the X
co factor, an issue we examine in detail in Paper II.
We report the results of a pilot study with the Expanded Very Large Array (EVLA) of 12CO J= 1-0 emission from four submillimetre-selected galaxies at z= 2.2-2.5, each with an existing detection of ...12CO J= 3-2, one of which comprises two distinct spatial components. Using the EVLA's most compact configuration, we detect strong, broad medians: 990 km s−1 full width at zero intensity; 540 km s−1 full width at half-maximum (FWHM)J= 1-0 line emission from all of our targets - coincident in position and velocity with their J= 3-2 emission. The median linewidth ratio, σ1-0/σ3-2= 1.15 ± 0.06, suggests that the J= 1-0 is more spatially extended than the J= 3-2 emission, a situation confirmed by our maps which reveal velocity structure in several cases and typical sizes of ∼16 kpc FWHM. The median brightness temperature (T
b) ratio is r
3-2/1-0= 0.55 ± 0.05, consistent with local galaxies with L
IR > 1011 L⊙, noting that our value may be biased high because of the J= 3-2 based sample selection. Naively, this suggests gas masses roughly two times higher than estimates made using higher J transitions of CO, with the discrepancy due entirely to the difference in assumed T
b ratio. We also estimate molecular gas masses using the 12CO J= 1-0 line and the observed global T
b ratios, assuming standard underlying T
b ratios for the non-star-forming and star-forming gas phases as well as a limiting star formation efficiency for the latter in all systems, i.e. without calling upon X
CO (≡α). Using this new method, we find a median molecular gas mass of (2.5 ± 0.8) × 1010 M⊙, with a plausible range stretching up to three times higher. Even larger masses cannot be ruled out, but are not favoured by dynamical constraints: the median dynamical mass within R∼ 7 kpc for our sample is (2.3 ± 1.4) × 1011 M⊙ or ∼6 times more massive than UV-selected galaxies at this epoch. We examine the Schmidt-Kennicutt (S-K) relation for all the distant galaxy populations for which CO J= 1-0 or J= 2-1 data are available, finding small systematic differences between galaxy populations. These have previously been interpreted as evidence for different modes of star formation, but we argue that these differences are to be expected, given the still considerable uncertainties, certainly when considering the probable excitation biases due to the molecular lines used, and the possibility of sustained S-K offsets during the evolution of individual gas-rich systems. Finally, we discuss the morass of degeneracies surrounding molecular gas mass estimates, the possibilities for breaking them, and the future prospects for imaging and studying cold, quiescent molecular gas at high redshifts.
Aims. The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the ...DustPedia sample. Methods. This study is performed by exploiting the multi-wavelength DustPedia database, from ultraviolet (UV) to sub-millimeter bands, in addition to molecular (12CO) and atomic (Hi) gas maps and metallicity abundance information available in the literature. We fitted the surface-brightness profiles of the tracers of dust and stars, the mass surface-density profiles of dust, stars, molecular gas, and total gas, and the SFR surface-density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (αCO) per galaxy by using dust- and gas-mass profiles. Results. Although each galaxy has its own peculiar behavior, we identified a common trend of the exponential scale-lengths versus wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec2 radius decrease from UV to 70 μm, from 0.4 to 0.2, and then increase back up to ~0.3 at 500 microns. The main result is that, on average, the dust-mass surface-density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 μm surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 μm scale-length tend to increase from earlier to later types, the scale-length at 70 μm tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 μm that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types. Our αCO determinations are in the range (0.3−9) M⊙ pc-2 (K km s-1)-1, almost invariant by using a fixed dust-to-gas ratio mass (DGR) or a DGR depending on metallicity gradient.
Aims.
The purpose of this work is the characterization of the main scaling relations between all of the interstellar medium (ISM) components, namely dust, atomic, molecular, and total gas, and ...gas-phase metallicity, as well as other galaxy properties, such as stellar mass (
M
star
) and galaxy morphology, for late-type galaxies in the Local Universe.
Methods.
This study was performed by extracting late-type galaxies from the entire DustPedia sample and by exploiting the large and homogeneous dataset available thanks to the DustPedia project. The sample consists of 436 galaxies with morphological stage spanning from
T
= 1−10,
M
star
from 6 × 10
7
to 3 × 10
11
M
⊙
, star formation rate from 6 × 10
−4
to 60
M
⊙
yr
−1
, and oxygen abundance from 12 + log(O/H) = 8−9.5. Molecular and atomic gas data were collected from the literature and properly homogenized. All the masses involved in our analysis refer to the values within the optical disks of galaxies. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent CO-to-H
2
conversion factor (
X
CO
). The analysis was performed by means of the survival analysis technique, in order to properly take into account the presence of both detection and nondetection in the data.
Results.
We confirm that the dust mass correlates very well with the total gas mass, and find –for the first time– that the dust mass correlates better with the atomic gas mass than with the molecular one. We characterize important mass ratios such as the gas fraction, the molecular-to-atomic gas mass ratio, the dust-to-total gas mass ratio (DGR), and the dust-to-stellar mass ratio, and study how they relate to each other, to galaxy morphology, and to gas-phase metallicity. Only the assumption of a metallicity-dependent
X
CO
reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. The DGR, the gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for galaxies rich in molecular gas this trend breaks down. To our knowledge, this trend has never been found before, and provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies. We discuss several scenarios related to this finding.
Conclusions.
The DustPedia database of late-type galaxies is an extraordinary tool for the study of the ISM scaling relations, thanks to its homogeneous collection of data for the different ISM components. The database is made publicly available to the whole community.
Aims. We aim to study the fraction of stellar radiation absorbed by dust, fabs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia ...sample, including almost all large (optical diameter larger than 1′), nearby (v ≤ 3000 km s−1) galaxies observed with the Herschel Space Observatory. Methods. For each object, we modelled the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the Code Investigating GALaxy Evolution (CIGALE). The value of fabs was obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. Results. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of fabs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of fabs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of fabs with stellar mass and specific star-formation rate. We confirm that the local Universe is missing the high fabs, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light.
ABSTRACT
The dust mass absorption coefficient, κd is the conversion function used to infer physical dust masses from observations of dust emission. However, it is notoriously poorly constrained, and ...it is highly uncertain how it varies, either between or within galaxies. Here we present the results of a proof-of-concept study, using the DustPedia data for two nearby face-on spiral galaxies M 74 (NGC 628) and M 83 (NGC 5236), to create the first ever maps of κd in galaxies. We determine κd using an empirical method that exploits the fact that the dust-to-metals ratio of the interstellar medium is constrained by direct measurements of the depletion of gas-phase metals. We apply this method pixel-by-pixel within M 74 and M 83, to create maps of κd. We also demonstrate a novel method of producing metallicity maps for galaxies with irregularly sampled measurements, using the machine learning technique of Gaussian process regression. We find strong evidence for significant variation in κd. We find values of κd at 500 $\mu$m spanning the range 0.11–0.25 ${\rm m^{2}\, kg^{-1}}$ in M 74, and 0.15–0.80 ${\rm m^{2}\, kg^{-1}}$ in M 83. Surprisingly, we find that κd shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method – only the adoption of unphysical or highly unusual assumptions would be able to suppress it.
Aims. Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with stellar radiation (originating from different stellar populations) for 814 ...galaxies in the nearby Universe, all observed by the Herschel Space Observatory. Methods. We take advantage of the widely used fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. For comparison purposes, an estimation of the dust properties is provided by approximating the emission at far-infrared and sub-millimeter wavelengths with a modified blackbody. Using the DustPedia photometry we determine the physical properties of the galaxies, such as the dust and stellar mass, the star-formation rate, the bolometric luminosity, the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (≤200 Myr) stellar populations. Results. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ∼77% of the total unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ∼24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ∼90%. Young stars progressively contribute more for “typical” spiral galaxies and they become the dominant source of dust heating for Sm-type and irregular galaxies, with ∼60% of their luminosity contributing to that purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate.
The Herschel Virgo Cluster Survey (HeViCS) is the deepest, confusion-limited survey of the Virgo Cluster at far-infrared (FIR) wavelengths. The entire survey at full depth covers ∼55 deg2 in five ...bands (100-500 μm), encompassing the areas around the central dominant elliptical galaxies (M87, M86 and M49) and extends as far as the NW cloud, the W cloud and the Southern extension. The survey extends beyond this region with lower sensitivity so that the total area covered is 84 deg2. In this paper we describe the data, the data acquisition techniques and present the detection rates of the optically selected Virgo Cluster Catalogue (VCC). We detect 254 (34 per cent) of 750 VCC galaxies found within the survey boundary in at least one band and 171 galaxies are detected in all five bands. For the remainder of the galaxies we have measured strict upper limits for their FIR emission. The population of detected galaxies contains early as well as late types although the latter dominate the detection statistics. We have modelled 168 galaxies, showing no evidence of a strong synchrotron component in their FIR spectra, using a single-temperature modified blackbody spectrum with a fixed emissivity index (β = 2). A study of the χ2 distribution indicates that this model is not appropriate in all cases, and this is supported by the FIR colours which indicate a spread in β = 1-2. Statistical comparison of the dust mass and temperature distributions from 140 galaxies with χ2
d.o.f. = 3 < 7.8 (95 per cent confidence level) shows that late types have typically colder, more massive dust reservoirs; the early-type dust masses have a mean of log〈M〉/M = 6.3 ± 0.3, while for late types log〈M〉/M = 7.1 ± 0.1. The late-type dust temperatures have a mean of 〈T〉 = 19.4 ± 0.2 K, while for the early types, 〈T〉 = 21.1 ± 0.8 K. Late-type galaxies in the cluster exhibit slightly lower dust masses than those in the field, but the cluster environment seems to have little effect on the bulk dust properties of early types. In future papers we will focus more on the scientific analysis of the catalogue (e.g. measuring FIR luminosity functions, dust mass functions and resolved gas and dust properties).
Aims.
The general consensus is that late-type galaxies undergo intense star-formation, activity while early-type galaxies are mostly inactive. We question this general rule and investigate the ...existence of star-forming early-type and quiescent late-type galaxies in the local Universe. By computing the physical properties of these galaxies and by using information on their structural properties as well as the density of their local environment, we seek to understand the differences from their ‘typical’ counterparts.
Methods.
We made use of the multi-wavelength photometric data (from the ultraviolet to the sub-millimetre), for 2209 morphologically classified galaxies in the Galaxy And Mass Assembly survey. Furthermore, we separated the galaxies into subsets of star-forming and quiescent based on their dominant ionising process, making use of established criteria based on the
W
H
α
width and the N
II
/H
α
ratio. Taking advantage of the spectral energy distribution fitting code
CIGALE
, we derived galaxy properties, such as the stellar mass, dust mass, and star-formation rate, and we also estimated the unattenuated and the dust-absorbed stellar emission, for both the young (≤200 Myr) and old (> 200 Myr) stellar populations.
Results.
We find that about 47% of E/S0 galaxies in our sample show ongoing star-formation activity and 8% of late-type galaxies are quiescent. The star-forming elliptical galaxies, together with the little blue spheroids, constitute a population that follows the star-forming main sequence of spiral galaxies very well. The fraction of the luminosity originating from young stars in the star-forming early-type galaxies is quite substantial (∼25%) and similar to that of the star-forming late-type galaxies. The stellar luminosity absorbed by the dust (and used to heat the dust grains) is highest in star-forming E/S0 galaxies (an average of 35%) followed by star-forming Sa-Scd galaxies (27%) with this fraction becoming significantly smaller for their quiescent analogues (6% and 16%, for E/S0 and Sa-Scd, respectively). Star-forming and quiescent E/S0 galaxies donate quite different fractions of their young stellar luminosities to heat up the dust grains (74% and 36%, respectively), while these fractions are very similar for star-forming and quiescent Sa-Scd galaxies (59% and 60%, respectively). Investigating possible differences between star-forming and quiescent galaxies, we find that the intrinsic (unattenuated) shape of the SED of the star-forming galaxies is, on average, very similar for all morphological types. Concerning their structural parameters, quiescent galaxies tend to show larger values of the
r
-band Sérsic index and larger effective radii (compared to star-forming galaxies). Finally, we find that star-forming galaxies preferably reside in lower density environments compared to the quiescent ones, which exhibit a higher percentage of sources being members of groups.