Context. The determination of the relative frequency of active galactic nuclei (AGN) versus other spectral classes, for example, HII region-like (HII), transition objects (TRAN), passive (PAS), and ...retired (RET), in a complete set of galaxies in the local Universe is of primary importance to discriminate the source of ionization in the nuclear region of galaxies (e.g., supermassive black holes vs. young and old stars). Aims. Here we aim to provide a spectroscopic characterization of the nuclei of galaxies belonging to the Herschel Reference Survey (HRS), a volume and magnitude limited sample representative of the local Universe, which has become a benchmark for local and high-z studies, for semianalytical models and cosmological simulations. The comparison between the nuclear spectral classification and the one determined on the global galactic scale provides information about how galaxy properties change from the nuclear to the outer regions. Moreover, the extrapolation of the global star formation (SF) properties from the SDSS fiber spectroscopy compared to the one computed by Hα photometry can be useful for testing the method based on aperture correction for determining the global star formation rate for local galaxies. Methods. By collecting the existing nuclear spectroscopy available from the literature, complemented with new observations obtained using the Loiano 1.52 m telescope, we analyze the 322 nuclear spectra of HRS galaxies; their integrated spectroscopy is available from the literature as well. Results. Using two diagnostic diagrams (the BPT and the WHAN) we provide a nuclear and an integrated spectral classification for the HRS galaxies. The BPT and the WHAN methods for nuclei consistently give a frequency of 53–64% HII, around 21–27% AGNs (including TRAN), and 15–20% of PAS (including RET), whereas for integrated spectra they give 69–84% HII, 4–11% of AGNs and 12–20% PAS. Solely among late-type galaxies (LTGs) do the nuclear percentages become 67–77% HII, 22–27% AGNs (including TRAN), and only 1–7% of PAS. For the integrated spectra these frequencies become: 80–85% HII, 9–11% AGNs and 4–9% PAS. Conclusions. We find that the fraction of HII region-like spectra is strongly anticorrelated with the stellar mass. On the contrary the frequency of AGNs increases significantly with stellar mass, such that at M* > 1010.0 M⊙~ 66% of the LTGs are AGNs or TRAN. Moreover there is not a significant dependence of the frequency of AGNs as a function of environment: AGNs+TRAN above 109.0 M⊙ are consistent with ~30% irrespective of their membership to the Virgo cluster, suggesting that the AGNs population is not sensitive to the environment. Finally, extrapolation of the global SF properties from the nuclear spectroscopy including aperture corrections leads to underestimates with respect to values derived from direct integrated Hα photometry.
Aims.
We construct the molecular mass function using the bivariate
K
-band-mass function (BMF) of the
Herschel
Reference Survey (HRS), which is a volume-limited sample that has already been widely ...studied at the entire electromagnetic spectrum.
Methods.
The molecular mass function was derived from the
K
-band and the gas mass cumulative distribution using a copula method, which is described in detail in our previous papers.
Results.
The H
2
mass is relatively strongly correlated with the
K
-band luminosity because of the tight relation between the stellar mass and the molecular gas mass within the sample with a scatter, which is likely due to those galaxies which have lost their molecular content because of environmental effects or because of a larger gas consumption due to past star formation processes. The derived H
2
MF samples the molecular mass range from ∼4 × 10
6
M
⊙
to ∼10
10
M
⊙
, and when compared with theoretical models, it agrees well with the theoretical predictions at the lower end of the mass values; whereas at masses larger than 10
10
M
⊙
, the HRS sample may miss galaxies with a large content of molecular hydrogen and the outcomes are not conclusive. The value of the local density of the molecular gas mass inferred from our analysis is ∼1.5 × 10
7
M
⊙
Mpc
−3
, and it is compared with the results at larger redshifts, confirming the lack of strong evolution for the molecular mass density between
z
= 0 and
z
= 4.
Conclusions.
This is the first molecular mass function that has been derived on a complete sample in the local Universe, which can be used as a reliable calibration at redshift
z
= 0 for models aiming to predict the evolution of the molecular mass density.
We use H Delta *a and far-ultraviolet (FUV, 1539 A) Galaxy Evolution Explorer (GALEX) data for a large sample of nearby objects to study the high-mass (m>= 2 M ) star formation activity of normal ...late-type galaxies. The data are corrected for dust attenuation using the most accurate techniques available at present, namely the Balmer decrement for H Delta *a data and the total far-infrared to FUV flux ratio for GALEX data. The sample shows a highly dispersed distribution in the H Delta *a to FUV flux ratio (log f(H Delta *a)/f(FUV) = 1.10 +/- 0.34 A) indicating that two of the most commonly used star formation tracers give star formation rates (SFRs) with uncertainties up to a factor of 2-3. The high dispersion is partly due to the presence of active galactic nuclei, where the UV and the H Delta *a emission can be contaminated by nuclear activity, highly inclined galaxies, for which the applied extinction corrections are probably inaccurate, or starburst galaxies, where the stationarity in the star formation history required for transforming H Delta *a and UV luminosities into SFRs is not satisfied. Excluding these objects, normal late-type galaxies have log f(H Delta *a)/f(FUV) = 0.94 +/- 0.16 A, which corresponds to an uncertainty of ~50% on the SFR. The H Delta *a to FUV flux ratio of the observed galaxies increases with their total stellar mass. If limited to normal star-forming galaxies, however, this relationship reduces to a weak trend that might be totally removed using different extinction correction recipes. In these objects, the H Delta *a to FUV flux ratio seems also barely related to the FUV - H color, the H-band effective surface brightness, the total star formation activity, and the gas fraction. The data are consistent with a Kroupa and Salpeter initial mass function (IMF) in the high-mass stellar range (m > 2 M ) and imply, for a Salpeter IMF, that the variations of the slope Delta *g cannot exceed 0.25, from Delta *g = 2.35 for massive galaxies to Delta *g = 2.60 in low luminosity systems. We show however that these observed trends, if real, can be due to the different micro-history of star formation in massive galaxies with respect to dwarf systems.
Although it accounts only for a small fraction of the baryonic mass, dust has a profound impact on the physical processes at play in galaxies. Thus, to understand the evolution of galaxies, it is ...essential not only to characterize dust properties per se, but also in relation to global galaxy properties. To do so, we derive the dust properties of galaxies in a volume limited, K-band selected sample, the Herschel Reference Survey (HRS). We gather infrared photometric data from 8 μm to 500 μm from Spitzer, WISE, IRAS, and Herschel for all of the HRS galaxies. Draine & Li (2007, ApJ, 663, 866) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the parameters of the Draine & Li models and that a strong constraint on the 20−60 μm range is mandatory to estimate the relative contribution of the photo-dissociation regions to the infrared spectral energy distribution (SED). The SED models tend to systematically underestimate the observed 500 μm flux densities, especially for low-mass systems. We provide the output parameters for all of the galaxies, i.e., the minimum intensity of the interstellar radiation field, the fraction of polycyclic aromatic hydrocarbon (PAH), the relative contribution of PDR and evolved stellar population to the dust heating, the dust mass, and the infrared luminosity. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the main integrated properties of galaxies, such as stellar mass, star formation rate, infraredluminosity, metallicity, Hα and H-band surface brightness, and the far-ultraviolet attenuation. A good correlation between the fraction of PAH and the metallicity is found, implying a weakening of the PAH emission in galaxies with low metallicities and, thus, low stellar masses. The intensity of the diffuse interstellar radiation field and the H-band and Hα surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star-forming regions and the diffuse evolved population. We use these results to provide a new set of infrared templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z = 0 reference for cosmological studies. For the same purpose, we place our sample on the SFR − M∗ diagram. The templates are compared to the most popular infrared SED libraries, enlightening a large discrepancy between all of them in the 20−100 μm range.
Abstract
The Herschel Multi-tiered Extragalactic Survey (HerMES) is a legacy programme designed to map a set of nested fields totalling ∼380 deg2. Fields range in size from 0.01 to ∼20 deg2, using ...the Herschel-Spectral and Photometric Imaging Receiver (SPIRE) (at 250, 350 and 500 μm) and the Herschel-Photodetector Array Camera and Spectrometer (PACS) (at 100 and 160 μm), with an additional wider component of 270 deg2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the reprocessed optical and ultraviolet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multiwavelength understanding of galaxy formation and evolution.
The survey will detect of the order of 100 000 galaxies at 5σ in some of the best-studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to facilitate redshift determination, rapidly identify unusual objects and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include the total infrared emission of galaxies, the evolution of the luminosity function, the clustering properties of dusty galaxies and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques.
This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.
We exploit the deep and extended far-IR data sets (at 70, 100 and 160 μm) of the Herschel Guaranteed Time Observation (GTO) PACS Evolutionary Probe (PEP) Survey, in combination with the Herschel ...Multi-tiered Extragalactic Survey data at 250, 350 and 500 μm, to derive the evolution of the rest-frame 35-, 60-, 90- and total infrared (IR) luminosity functions (LFs) up to z ∼ 4. We detect very strong luminosity evolution for the total IR LF (L
IR ∝ (1 + z)3.55 ± 0.10 up to z ∼ 2, and ∝ (1 + z)1.62 ± 0.51 at 2 < z 4) combined with a density evolution (∝(1 + z)−0.57 ± 0.22 up to z ∼ 1 and ∝ (1 + z)−3.92 ± 0.34 at 1 < z 4). In agreement with previous findings, the IR luminosity density (ρIR) increases steeply to z ∼ 1, then flattens between z ∼ 1 and z ∼ 3 to decrease at z 3. Galaxies with different spectral energy distributions, masses and specific star formation rates (SFRs) evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to ρIR. Galaxies occupying the well-established SFR-stellar mass main sequence (MS) are found to dominate both the total IR LF and ρIR at all redshifts, with the contribution from off-MS sources (≥0.6 dex above MS) being nearly constant (∼20 per cent of the total ρIR) and showing no significant signs of increase with increasing z over the whole 0.8 < z < 2.2 range. Sources with mass in the range 10 ≤ log(M/M) ≤ 11 are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z ( 2) LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the Super Massive Black Holes (SMBH) grows and is obscured by dust (possibly triggered by a major merging event), is followed by an AGN-dominated phase, then evolving towards a local elliptical. On the other hand, moderately star-forming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.
The radiation from stars heats dust grains in the diffuse interstellar medium and in star-forming regions in galaxies. Modelling this interaction provides information on dust in galaxies, a vital ...ingredient for their evolution. It is not straightforward to identify the stellar populations heating the dust, and to link attenuation to emission on a sub-galactic scale. Radiative transfer models are able to simulate this dust-starlight interaction in a realistic, three-dimensional setting. We investigate the dust heating mechanisms on a local and global galactic scale, using the Andromeda galaxy (M31) as our laboratory. We have performed a series of panchromatic radiative transfer simulations of Andromeda with our code SKIRT. The high inclination angle of M31 complicates the 3D modelling and causes projection effects. However, the observed morphology and flux density are reproduced fairly well from UV to sub-millimeter wavelengths. Our model reveals a realistic attenuation curve, compatible with previous, observational estimates. We find that the dust in M31 is mainly (91% of the absorbed luminosity) heated by the evolved stellar populations. The bright bulge produces a strong radiation field and induces non-local heating up to the main star-forming ring at 10 kpc. The relative contribution of unevolved stellar populations to the dust heating varies strongly with wavelength and with galactocentric distance. The dust heating fraction of unevolved stellar populations correlates strongly with NUV-r colour and specific star formation rate. These two related parameters are promising probes for the dust heating sources at a local scale.
The Next Generation Virgo Cluster Survey (NGVS) is a program that uses the 1 deg super(2) MegaCam instrument on the Canada-France-Hawaii Telescope to carry out a comprehensive optical imaging survey ...of the Virgo cluster, from its core to its virial radius-covering a total area of 104 deg super(2)-in the u*griz bandpasses. Thanks to a dedicated data acquisition strategy and processing pipeline, the NGVS reaches a point-source depth of g approx = 25.9 mag (10sigma) and a surface brightness limit of mu sub(g) ~ 29 mag arcsec super(-2) (2sigma above the mean sky level), thus superseding all previous optical studies of this benchmark galaxy cluster. In this paper, we give an overview of the technical aspects of the survey, such as areal coverage, field placement, choice of filters, limiting magnitudes, observing strategies, data processing and calibration pipelines, survey timeline, and data products. We also describe the primary scientific topics of the NGVS, which include: the galaxy luminosity and mass functions; the color-magnitude relation; galaxy scaling relations; compact stellar systems; galactic nuclei; the extragalactic distance scale; the large-scale environment of the cluster and its relationship to the Local Supercluster; diffuse light and the intracluster medium; galaxy interactions and evolutionary processes; and extragalactic star clusters. In addition, we describe a number of ancillary programs dealing with "foreground" and "background" science topics, including the study of high-inclination trans-Neptunian objects; the structure of the Galactic halo in the direction of the Virgo Overdensity and Sagittarius Stream; the measurement of cosmic shear, galaxy-galaxy, and cluster lensing; and the identification of distant galaxy clusters, and strong-lensing events.
Aims. Our knowledge of the cosmic mass assembly relies on measurements of star formation rates (SFRs) and stellar masses (Mstar), of galaxies as a function of redshift. These parameters must be ...estimated in a consistent way with a good knowledge of systematics before studying their correlation and the variation of the specific SFR. Constraining these fundamental properties of galaxies across the Universe is of utmost importance if we want to understand galaxy formation and evolution. Methods. We seek to derive SFRs and stellar masses in distant galaxies and to quantify the main uncertainties affecting their measurement. We explore the impact of the assumptions made in their derivation with standard calibrations or through a fitting process, as well as the impact of the available data, focusing on the role of infrared emission originating from dust. Results. We build a sample of galaxies with z > 1, all observed from the ultraviolet to the infrared in their rest frame. The data are fitted with the code CIGALE, which is also used to build and analyse a catalogue of mock galaxies. Models with different star formation histories are introduced: an exponentially decreasing or increasing SFR and a more complex one coupling a decreasing SFR with a younger burst of constant star formation. We define different sets of data, with or without a good sampling of the ultraviolet range, near-infrared, and thermal infrared data. Variations of the metallicity are also investigated. The impact of these different cases on the determination of stellar mass and SFR are analysed. Conclusions. Exponentially decreasing models with a redshift formation of the stellar population zf ≃ 8 cannot fit the data correctly. All the other models fit the data correctly at the price of unrealistically young ages when the age of the single stellar population is taken to be a free parameter, especially for the exponentially decreasing models. The best fits are obtained with two stellar populations. As long as one measurement of the dust emission continuum is available, SFR are robustly estimated whatever the chosen model is, including standard recipes. The stellar mass measurement is more subject to uncertainty, depending on the chosen model and the presence of near-infrared data, with an impact on the SFR-Mstar scatter plot. Conversely, when thermal infrared data from dust emission are missing, the uncertainty on SFR measurements largely exceeds that of stellar mass. Among all physical properties investigated here, the stellar ages are found to be the most difficult to constrain and this uncertainty acts as a second parameter in SFR measurements and as the most important parameter for stellar mass measurements.
We present new Herschel-SPIRE imaging spectroscopy (194-671 mu m) of the bright starburst galaxy M82. Covering the CO ladder from J = 4 arrow right 3 to J = 13 arrow right 12, spectra were obtained ...at multiple positions for a fully sampled ~3 x 3 arcmin map, including a longer exposure at the central position. We present measurements of super(12)CO, super(13)CO, CI, NII, HCN, and HCO+ in emission, along with OH+, H sub(2)O+, and HF in absorption and H sub(2)O in both emission and absorption, with discussion. We use a radiative transfer code and Bayesian likelihood analysis to model the temperature, density, column density, and filling factor of multiple components of molecular gas traced by super(12)CO and super(13)CO, adding further evidence to the high-J lines tracing a much warmer (~500 K), less massive component than the low-/lines. The addition of super(13)CO (and CI) is new and indicates that CI may be tracing different gas than super(12)CO. No temperature/density gradients can be inferred from the map, indicating that the single-pointing spectrum is descriptive of the bulk properties of the galaxy. At such a high temperature, cooling is dominated by molecular hydrogen. Photon-dominated region (PDR) models require higher densities than those indicated by our Bayesian likelihood analysis in order to explain the high-J line ratios, though cosmic-ray-enhanced PDR models can do a better job reproducing the emission at lower densities. Shocks and turbulent heating are likely required to explain the bright high-J emission.