We present results of the detailed dust energy balance study for the seven large edge-on galaxies in the HEROES sample using 3D radiative transfer (RT) modelling. Based on available optical and near-infrared observations of the HEROES galaxies, we derive the 3D distribution of stars and dust in these galaxies. For the sake of uniformity, we apply the same technique to retrieve galaxy properties for the entire sample: we use a stellar model consisting of a S\'ersic bulge and three double-exponential discs (a superthin disc for a young stellar population and thin and thick discs for old populations). For the dust component, we adopt a double-exponential disc with the new THEMIS dust-grain model. We fit oligochromatic radiative transfer (RT) models to the optical and near-infrared images with the fitting algorithm FitSKIRT and do panchromatic simulations with the SKIRT code at wavelengths ranging from ultraviolet to submillimeter. We confirm the previously stated dust energy balance problem in galaxies: for the HEROES galaxies, the dust emission derived from our RT calculations underestimates the real observations by a factor 1.5-4 for all galaxies except NGC 973 and NGC 5907 (apparently, the latter galaxy has a more complex geometry than we used). The comparison between our RT simulations and the observations at mid-infrared-submillimeter wavelengths shows that most of our galaxies exhibit complex dust morphologies (possible spiral arms, star-forming regions, more extended dust structure in the radial and vertical directions). We suggest that, in agreement with the results from Saftly et al. (2015), the large- and small-scale structure is the most probable explanation for the dust energy balance problem.

We investigate the dust energy balance for the edge-on galaxy IC 2531, one of the seven galaxies in the HEROES sample. We perform a state-of-the-art radiative transfer modelling based, for the first time, on a set of optical and near-infrared galaxy images. We show that taking into account near-infrared imaging in the modelling significantly improves the constraints on the retrieved parameters of the dust content. We confirm the result from previous studies that including a young stellar population in the modelling is important for explaining the observed stellar energy distribution. However, the discrepancy between the observed and modelled thermal emission at far-infrared wavelengths, the so-called dust energy balance problem, is still present: the model underestimates the observed fluxes by a factor of about two. We compare two different dust models, and find that dust parameters and thus the spectral energy distribution in the infrared domain are sensitive to the adopted dust model. In general, the THEMIS model reproduces the observed emission in the infrared wavelength domain better than the popular Zubko et al. BARE-GR-S model. Our study of IC 2531 is a pilot case for detailed and uniform radiative transfer modelling of the entire HEROES sample, which will shed more light on the strength and origins of the dust energy balance problem.

Context. Edge-on galaxies can offer important insights in galaxy evolution as they are the only systems where the distribution of the different components can be studied both radially and vertically. The HEROES project was designed to investigate the interplay between the gas, dust, stars and dark matter (DM) in a sample of 7 massive edge-on spiral galaxies. Aims. In this second HEROES paper we present an analysis of the atomic gas content of 6 out of 7 galaxies in our sample. The remaining galaxy was recently analysed according to the same strategy. The primary aim of this work is to constrain the surface density distribution, the rotation curve and the geometry of the gas disks in a homogeneous way. In addition we identify peculiar features and signs of recent interactions. Methods. We construct detailed tilted-ring models of the atomic gas disks based on new GMRT 21-cm observations of NGC 973 and UGC 4277 and re-reduced archival HI data of NGC 5907, NGC 5529, IC 2531 and NGC 4217. Potential degeneracies between different models are resolved by requiring a good agreement with the data in various representations of the data cubes. Results. From our modelling we find that all but one galaxy are warped along the major axis. In addition, we identify warps along the line of sight in three galaxies. A flaring gas layer is required to reproduce the data only for one galaxy, but (moderate) flares cannot be ruled for the other galaxies either. A coplanar ring-like structure is detected outside the main disk of NGC 4217, which we suggest could be the remnant of a recent minor merger event. We also find evidence for a radial inflow of 15 +- 5 km/s in the disk of NGC 5529, which might be related to the ongoing interaction with two nearby companions. (Abridged)

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

To assess the effects of the cluster environment on the different components of the interstellar medium, we analyse the FIR-submm properties of a sample of star-forming dwarf (SFD) galaxies detected by the Herschel Virgo Cluster Survey (HeViCS). We determine dust masses and dust temperatures by fitting a modified black body (MBB) function to the spectral energy distributions (SEDs). Stellar and gas masses, star formation rates (SFRs), and metallicities are obtained from the analysis of a set of ancillary data. Dust is detected in 49 out of 140 optically identified dwarfs covered by the HeViCS field; considering only dwarfs brighter than $m_B$ = 18 mag, this gives a detection rate of 43%. After evaluating different emissivity indices, we find that the FIR-submm SEDs are best-fit by $\beta$=1.5, with a median dust temperature $T_d$ = 22.4 K. Assuming $\beta$=1.5, 67% of the 23 galaxies detected in all five Herschel bands show emission at 500 $\mu$m in excess of the MBB model. The excess is inversely correlated with SFR and stellar masses. To study the variations in the global properties of our sample due to environmental effects, we compare the Virgo SFDs to other Herschel surveys, such as KINGFISH, the Dwarf Galaxy Survey (DGS), and the HeViCS bright galaxy catalogue (BGC). We explore the relations between stellar mass and HI fraction, specific SFR, dust fraction, gas-to-dust ratio over a wide range of stellar masses. Highly HI-deficient Virgo dwarf galaxies are mostly characterised by quenched star formation activity and lower dust fractions giving hints for dust stripping in cluster dwarfs. However, we find that the fraction of dust removed has to be less than that of the HI component. Since the Virgo SFDs are likely to be crossing the cluster for the first time, a longer timescale might be necessary to strip the more centrally concentrated dust distribution.

We present three independent catalogs of point-sources extracted from SPIRE images at 250, 350, and 500 micron as a part of the Herschel Virgo Cluster Survey (HeViCS). The source positions are determined by estimating the likelihood to be a real source for each peak on the maps and the flux densities are estimated using the sourceExtractorTimeline, a timeline-based point source fitter. Afterwards, each source is subtracted from the maps, removing a Gaussian function in every position with the full width half maximum equal to that estimated in sourceExtractorTimeline. This procedure improves the robustness of our algorithm in terms of source identification. The HeViCS catalogs contain about 52000, 42200, and 18700 sources selected at 250, 350, and 500 micron above 3sigma and are ~ 75%, 62%, and 50% complete at flux densities of 20 mJy at 250, 350, 500 micron, respectively. We then measured source number counts at 250, 350, and 500 micron and we also cross-correlated the catalogs with the Sloan Digital Sky Survey to investigate the redshift distribution of the nearby sources. From this cross-correlation, we select ~2000 sources with reliable fluxes and a high signal-to-noise ratio, finding an average redshift z~0.3+/-0.22. The number counts at 250, 350, and 500 micron show an increase in the slope below 200 mJy, indicating a strong evolution in number of density for galaxies at these fluxes. In general, models tend to overpredict the counts at brighter flux densities, underlying the importance of studying the Rayleigh-Jeans part of the spectral energy distribution to refine the theoretical recipes of the models. Our iterative method for source identification allowed the detection of a family of 500 micron sources that are not foreground objects belonging to Virgo and not found in other catalogs.

The imprints of dust-starlight interactions are visible in scaling relations between stellar characteristics, star formation parameters and dust properties. We aim to examine dust scaling relations on a sub-kpc resolution in the Andromeda galaxy (M31) by comparing the properties on a local and global scale to other galaxies of the local universe. New Herschel observations are combined with available data from GALEX, SDSS, WISE and Spitzer to construct a dataset covering UV to submm wavelengths. We work at the resolution of the SPIRE $500\; \mu$m beam, with pixels corresponding to physical regions of 137 x 608 pc in the galaxy's disk. A panchromatic spectral energy distribution was modelled for each pixel and several dust scaling relations are investigated. We find, on a sub-kpc scale, strong correlations between $M_d/M_\star$ and NUV-r, and between $M_d/M_\star$ and $\mu_\star$ (the stellar mass surface density). Striking similarities with corresponding relations based on integrated galaxies are found. We decompose M31 in four macro-regions based on their FIR morphology; the bulge, inner disk, star forming ring and the outer disk. All regions closely follow the galaxy-scale average trends. The specific star formation characteristics we derive for these macro-regions give strong hints of an inside-out formation of the bulge-disk geometry, as well as an internal downsizing process. However, within each macro-region, a great diversity in individual micro-regions is found. Furthermore, we confirm that dust in the bulge of M31 is heated only by the old stellar populations. In general, the local dust scaling relations indicate that the dust content in M31 is maintained by a subtle interplay of past and present star formation. The similarity with galaxy-based relations strongly suggests that they are in situ correlations, with underlying processes that must be local in nature. (Abriged)

In this paper we present a catalogue of Giant Molecular Clouds (GMCs) in the Andromeda (M31) galaxy extracted from the Hershel Exploitation of Local Galaxy Andromeda (HELGA) dataset. GMCs are identified from the Herschel maps using a hierarchical source extraction algorithm. We present the results of this new catalogue and characterise the spatial distribution and spectral energy properties of its clouds based on the radial dust/gas properties found by Smith et al (2012). 236 GMCs in the mass range 10^4-10^7 M_sol are identified, their cumulative mass distribution is found to be proportional to M^-1.45 in agreement with earlier studies. The GMCs appear to follow the same cloud mass to L_CO correlation observed in the Milky Way. However, comparison between this catalogue and interferometry studies also shows that the GMCs are substructured below the Herschel resolution limit suggesting that we are observing associations of GMCs. Following Gordon et al. (2006), we study the spatial structure of M31 by splitting the observed structure into a set of spiral arms and offset rings. We fit radii of 10.5 and 15.5 kpc to the two most prominent rings. We then fit a logarithmic spiral with a pitch angle of 8.9 deg to the GMCs not associated with either ring. Lastly, we comment upon the effects of deprojection on our results and investigate the effect different models for M31's inclination will have upon the projection of an unperturbed spiral arm system.

We investigate the connection between dust and gas in the nearby edge-on spiral galaxy NGC 891. High resolution Herschel PACS and SPIRE 70, 100, 160, 250, 350, and 500 $\mu$m images are combined with JCMT SCUBA 850 $\mu$m observations to trace the far-infrared/submillimetre spectral energy distribution (SED). Maps of the HI 21 cm line and CO(J=3-2) emission trace the atomic and molecular hydrogen gas, respectively. We fit one-component modified blackbody models to the integrated SED, finding a global dust mass of 8.5$\times$10$^{7}$ M$_{\odot}$ and an average temperature of 23$\pm$2 K. We also fit the pixel-by-pixel SEDs to produce maps of the dust mass and temperature. The dust mass distribution correlates with the total stellar population as traced by the 3.6 $\mu$m emission. The derived dust temperature, which ranges from approximately 17 to 24 K, is found to correlate with the 24 $\mu$m emission. Allowing the dust emissivity index to vary, we find an average value of $\beta$ = 1.9$\pm$0.3. We confirm an inverse relation between the dust emissivity spectral index and dust temperature, but do not observe any variation of this relationship with vertical height from the mid-plane of the disk. A comparison of the dust properties with the gaseous components of the ISM reveals strong spatial correlations between the surface mass densities of dust and the molecular hydrogen and total gas surface densities. Observed asymmetries in the dust temperature, and the H$_{2}$-to-dust and total gas-to-dust ratios hint that an enhancement in the star formation rate may be the result of larger quantities of molecular gas available to fuel star formation in the NE compared to the SW. Whilst the asymmetry likely arises from dust obscuration due to the geometry of the line-of-sight projection of the spiral arms, we cannot exclude an enhancement in the star formation rate in the NE side of the disk.

We cross-correlate the Planck Catalogue of Compact Sources (PCCS) with the fully sampled 84 deg2 Herschel Virgo Cluster Survey (HeViCS) fields. We search for and identify the 857 and 545 GHz PCCS sources in the HeViCS fields by studying their FIR/submm and optical counterparts. We find 84 and 48 compact Planck sources in the HeViCS fields at 857 and 545 GHz, respectively. Almost all sources correspond to individual bright Virgo Cluster galaxies. The vast majority of the Planck detected galaxies are late-type spirals, with the Sc class dominating the numbers, while early-type galaxies are virtually absent from the sample, especially at 545 GHz. We compare the HeViCS SPIRE flux densities for the detected galaxies with the four different PCCS flux density estimators and find an excellent correlation with the aperture photometry flux densities, even at the highest flux density levels. We find only seven PCCS sources in the HeViCS fields without a nearby galaxy as obvious counterpart, and conclude that all of these are dominated by Galactic cirrus features or are spurious detections. No Planck sources in the HeViCS fields seem to be associated to high-redshift proto-clusters of dusty galaxies or strongly lensed submm sources. Finally, our study is the first empirical confirmation of the simulation-based estimated completeness of the PCCS, and provides a strong support of the internal PCCS validation procedure.

Herschel FIR observations are used to construct Virgo cluster galaxy luminosity functions and to show that the cluster lacks the very bright and the numerous faint sources detected in field galaxy surveys. The far-infrared SEDs are fitted to obtain dust masses and temperatures and the dust mass function. The cluster is over dense in dust by about a factor of 100 compared to the field. The same emissivity (beta) temperature relation applies for different galaxies as that found for different regions of M31. We use optical and HI data to show that Virgo is over dense in stars and atomic gas by about a factor of 100 and 20 respectively. Metallicity values are used to measure the mass of metals in the gas phase. The mean metallicity is about 0.7 solar and 50% of the metals are in the dust. For the cluster as a whole the mass density of stars in galaxies is 8 times that of the gas and the gas mass density is 130 times that of the metals. We use our data to consider the chemical evolution of the individual galaxies, inferring that the measured variations in effective yield are due to galaxies having different ages, being affected to varying degrees by gas loss. Four galaxy scaling relations are considered: mass-metallicity, mass-velocity, mass-star formation rate and mass-radius - we suggest that initial galaxy mass is the prime driver of a galaxy's ultimate destiny. Finally, we use X-ray observations and galaxy dynamics to assess the dark and baryonic matter content compared to the cosmological model.

We use dust scaling relations to investigate the hypothesis that Virgo cluster transition-type dwarfs are infalling star-forming field galaxies, which is argued based on their optical features (e.g. disks, spiral arms, bars) and kinematic properties similar to late-type galaxies. After their infall, environmental effects gradually transform them into early-type galaxies through the removal of their interstellar medium and quenching of all star formation activity. In this paper, we aim to verify whether this hypothesis holds using far-infrared diagnostics based on Herschel observations of the Virgo cluster taken as part of the Herschel Virgo Cluster Survey (HeViCS). We select transition-type objects in the nearest cluster, Virgo, based on spectral diagnostics indicative for their residual or ongoing star formation. We detect dust Md ~ 10^{5-6} Msun in 36% of the transition-type dwarfs located on the high end of the stellar mass distribution. This suggests that the dust reservoirs present in non-detections fall just below the Herschel detection limit (< 1.1x10^5 Msun). Dust scaling relations support the hypothesis of a transformation between infalling late-type galaxies to quiescent low-mass spheroids governed by environmental effects, with dust-to-stellar mass fractions for transition-type dwarfs in between values characteristic for late-type objects and the lower dust fractions observed in early-type galaxies. Several transition-type dwarfs demonstrate blue central cores, hinting at the radially outside-in removal of gas and quenching of star formation activity. The fact that dust is also confined to the inner regions suggests that metals are stripped in the outer regions along with the gas. In the scenario of most dust being stripped from the galaxy along with the gas, we argue that... (abridged)

Context. Edge-on spiral galaxies with prominent dust lanes provide us with an excellent opportunity to study the distribution and properties of the dust within them. The HEROES project was set up to observe a sample of seven large edge-on galaxies across various wavelengths for this investigation. Aims. Within this first paper, we present the Herschel observations and perform a qualitative and quantitative analysis on them, and we derive some global properties of the far infrared and submillimetre emission. Methods. We determine horizontal and vertical profiles from the Herschel observations of the galaxies in the sample and describe the morphology. Modified black-body fits to the global fluxes, measured using aperture photometry, result in dust temperatures and dust masses. The latter values are compared to those that are derived from radiative transfer models taken from the literature. Results. On the whole, our Herschel flux measurements agree well with archival values. We find that the exponential horizontal dust distribution model often used in the literature generally provides a good description of the observed horizontal profiles. Three out of the seven galaxies show signatures of extended vertical emission at 100 and 160 {\mu}m at the 5{\sigma} level, but in two of these it is probably due to deviations from an exactly edge-on orientation. Only for NGC 4013, a galaxy in which vertically extended dust has already been detected in optical images, we can detect vertically extended dust, and the derived scaleheight agrees with the value estimated through radiative transfer modelling. Our analysis hints at a correlation between the dust scaleheight and its degree of clumpiness, which we infer from the difference between the dust masses as calculated from modelling of optical data and from fitting the spectral energy distribution of Herschel datapoints.

We present a detailed study of how the Star Formation Rate (SFR) relates to the interstellar medium (ISM) of M31 at ~140pc scales. The SFR is calculated using the far-ultraviolet and 24um emission, corrected for the old stellar population in M31. We find a global value for the SFR of 0.25+/-0.05Msun/yr and compare this with the SFR found using the total far-infrared (FIR) luminosity. There is general agreement in regions where young stars dominate the dust heating. Atomic hydrogen (HI) and molecular gas (traced by carbon monoxide, CO) or the dust mass is used to trace the total gas in the ISM. We show that the global surface densities of SFR and gas mass place M31 amongst a set of low-SFR galaxies in the plot of Kennicutt (1998b). The relationship between SFR and gas surface density is tested in six radial annuli across M31, assuming a power law relationship with index, N. The star formation law using total gas traced by HI and CO gives a global index of N=2.03+/-0.04, with a significant variation with radius; the highest values are observed in the 10kpc ring. We suggest that this slope is due to HI turning molecular at ~10Msun/pc2. When looking at H2 regions, we measure a higher mean SFR suggesting a better spatial correlation between H2 and SF. We find N~0.6 with consistent results throughout the disk - this is at the low end of values found in previous work and argues against a superlinear SF law on small scales.

We present Herschel observations of the Fornax cluster at 100, 160, 250, 350 and 500u with a spatial resolution of 7 - 36 arc sec (10" = 1 kpc at d_Fornax=17.9 Mpc). We define a sample of 11 bright galaxies, selected at 500u, directly comparable with our past work on Virgo. We find good agreement with previous observations made by IRAS and Planck. The FIR luminosity density is higher (factor of three) in Fornax compared to Virgo. The 100u (42.5-122.5u) luminosity is two orders of magnitude larger in Fornax than in the local field as measured by IRAS. Using stellar (L_{0.4-2.5}) and FIR (L_{100-500}) luminosities we estimate a mean optical depth of tau=0.4+/-0.1 - the same value as Virgo. For 10 of the 11 galaxies (NGC1399 excepted) we fit a modified blackbody curve (beta=2.0) to the SEDs to derive dust masses and temperatures of 10^{6.54-8.35} M_0 and T=14.6-24.2K respectively, comparable to Virgo. The derived stars-to-gas(atomic) and gas(atomic)-to-dust ratios vary from 1.1-67.6 and 9.8-436.5 respectively, again consistent with Virgo. Fornax is a mass overdensity in stars and dust of about 120 compared to the local field (30 for Virgo). Fornax and Virgo are both a factor of 6 lower over densities in gas(atomic) than in stars and dust indicating loss of gas, but not dust and stars, in the cluster environment. As the brightest FIR source in either Fornax and Virgo, NGC1365 is detected by Planck. The Planck data fit the PACS/SPIRE SED out to 1382u with no evidence of other sources of emission ('spinning dust', free-free, synchrotron). At the opposite end of the scale NGC1399 is detected only at 500$\mu$m with the emission probably arising from the nuclear radio source rather than inter-stellar dust.

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 $\sim$55 sq. deg. in 5 bands (100-500 \micron), encompassing the areas around the central dominant elliptical galaxies (M87, M86 & 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 sq. deg. 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%) 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 ($\beta = 2$). A study of the $\chi^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 $\beta$=1--2. Statistical comparison of the dust mass and temperature distributions from 140 galaxies with $\chi^2_{dof=3} < 7.8$ (95% confidence level) shows that late-types have typically colder, more massive dust reservoirs; the early-type dust masses have a mean of ${\rm log}(<M> / M_{\sun}) = 6.3 \pm 0.3 $, while for late-types ${\rm log}(<M> / M_{\sun}) =7.1 \pm 0.1$... (abridged)

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

We investigate the dust-to-gas mass ratio and the environmental effects on the various components of the interstellar medium for a spatially resolved sample of Virgo spirals. We have used the IRAM-30m telescope to map over their full extent NGC 4189, NGC 4298, NGC 4388, and NGC 4299 in the 12CO(1-0) and the 12CO(2-1) lines. We observed the same lines in selected regions of NGC 4351, NGC 4294, and NGC 4424. The CO observations are combined with Herschel maps in 5 bands between 100-500 {\mu}m from the HeViCS survey, and with HI data from the VIVA survey, to obtain spatially resolved dust and gas distributions. We studied the environmental dependencies by adding to our sample eight galaxies with 12CO(1-0) maps from the literature. We estimate the integrated mass of molecular hydrogen for the galaxies observed in the CO lines. We find molecular-to-total gas mass fractions between 0.04 \leq fmol \leq 0.65, with the lowest values for the dimmest galaxy in the B-band. The integrated dust-to-gas ratio ranges between 0.011 and 0.004. For the 12 mapped galaxies we derive the radial distributions of the atomic gas, molecular gas, and dust. We also study the effect of different CO-to-H2 conversion factors. Both the molecular gas and the dust distributions show steeper radial profiles for HI-deficient galaxies and the average dust-to-gas ratio for these galaxies increases or stays radially constant. On scales of \sim 3 kpc, we find a strong correlation between the molecular gas and the 250 micron surface brightness that is tighter than average for non-deficient galaxies. The correlation becomes linear if we consider the total gas surface mass density. However, the inclusion of atomic hydrogen does not improve the statistical significance of the correlation. The environment can modify the distributions of molecules and dust within a galaxy, although these components are more tightly bound than the atomic gas.

We present an analysis of the dust and gas in Andromeda, using Herschel images sampling the entire far-infrared peak. We fit a modified-blackbody model to ~4000 quasi-independent pixels with spatial resolution of ~140pc and find that a variable dust-emissivity index (beta) is required to fit the data. We find no significant long-wavelength excess above this model suggesting there is no cold dust component. We show that the gas-to-dust ratio varies radially, increasing from ~20 in the center to ~70 in the star-forming ring at 10kpc, consistent with the metallicity gradient. In the 10kpc ring the average beta is ~1.9, in good agreement with values determined for the Milky Way (MW). However, in contrast to the MW, we find significant radial variations in beta, which increases from 1.9 at 10kpc to ~2.5 at a radius of 3.1kpc and then decreases to 1.7 in the center. The dust temperature is fairly constant in the 10kpc ring (ranging from 17-20K), but increases strongly in the bulge to ~30K. Within 3.1kpc we find the dust temperature is highly correlated with the 3.6 micron flux, suggesting the general stellar population in the bulge is the dominant source of dust heating there. At larger radii, there is a weak correlation between the star formation rate and dust temperature. We find no evidence for 'dark gas' in M31 in contrast to recent results for the MW. Finally, we obtained an estimate of the CO X-factor by minimising the dispersion in the gas-to-dust ratio, obtaining a value of (1.9+/-0.4)x10^20 cm^-2 [K kms^-1]^-1.

Using the far-infrared emission, as observed by the Herschel Virgo Cluster Survey (HeViCS), and the integrated HI and CO brightness, we infer the dust and total gas mass for a magnitude limited sample of 35 metal rich spiral galaxies in Virgo. The CO flux correlates tightly and linearly with far-infrared fluxes observed by Herschel. Molecules in these galaxies are more closely related to cold dust rather than to dust heated by star formation or to optical/NIR brightness. We show that dust mass establishes a stronger correlation with the total gas mass than with the atomic or molecular component alone. The dust-to-gas ratio increases as the HI deficiency increases, but in highly HI deficient galaxies it stays constant. Dust is in fact less affected than atomic gas by weak cluster interactions, which remove most of the HI gas from outer and high latitudes regions. Highly disturbed galaxies, in a dense cluster environment, can instead loose a considerable fraction of gas and dust from the inner regions of the disk keeping constant the dust-to-gas ratio. There is evidence that the molecular phase is also quenched. This quencing becomes evident by considering the molecular gas mass per unit stellar mass. Its amplitude, if confirmed by future studies, highlights that molecules are missing in Virgo HI deficient spirals, but to a somewhat lesser extent than dust.

We present Herschel dust continuum, James Clerk Maxwell Telescope CO(3-2) observations and a search for [CII] 158 micron and [OI] 63 micron spectral line emission for the brightest early-type dwarf satellite of Andromeda, NGC 205. While direct gas measurements (Mgas ~ 1.5e+6 Msun, HI + CO(1-0)) have proven to be inconsistent with theoretical predictions of the current gas reservoir in NGC 205 (> 1e+7 Msun), we revise the missing interstellar medium mass problem based on new gas mass estimates (CO(3-2), [CII], [OI]) and indirect measurements of the interstellar medium content through dust continuum emission. Based on Herschel observations, covering a wide wavelength range from 70 to 500 micron, we are able to probe the entire dust content in NGC 205 (Mdust ~ 1.1-1.8e+4 Msun at Tdust ~ 18-22 K) and rule out the presence of a massive cold dust component (Mdust ~ 5e+5 Msun, Tdust ~ 12 K), which was suggested based on millimeter observations from the inner 18.4 arcsec. Assuming a reasonable gas-to-dust ratio of ~ 400, the dust mass in NGC 205 translates into a gas mass Mgas ~ 4-7e+6 Msun. The non-detection of [OI] and the low L_[CII]-to-L_CO(1-0) line intensity ratio (~ 1850) imply that the molecular gas phase is well traced by CO molecules in NGC 205. We estimate an atomic gas mass of 1.5e+4 Msun associated with the [CII] emitting PDR regions in NGC 205. From the partial CO(3-2) map of the northern region in NGC 205, we derive a molecular gas mass of M_H2 ~ 1.3e+5 Msun. [abridged]

We present Herschel observations of 62 Early-Type Galaxies (ETGs), including 39 galaxies morphologically classified as S0+S0a and 23 galaxies classified as ellipticals using SPIRE at 250, 350 and 500 microns (and PACS 100 and 160 microns for 19 sources) as part of the volume-limited Herschel Reference Survey. We detect dust emission in 24% of the ellipticals and 62% of the S0s. The mean temperature of the dust is 23.9+/-0.8 K, warmer than that found for late-type galaxies in the Virgo Cluster. Including the non-detections, the mean dust mass is log(Mdust) = 5.9+/-0.1 and 5.2+/-0.1 Msun for the S0s and elliptical galaxies respectively. The mean dust-to-stellar mass is log(Mdust/Mstar) = -4.4+/-0.1 (S0s) and -5.8+/-0.1 (ellipticals). Virtually all the galaxies lie close to the red sequence yet the large number of detections of cool dust, the gas-to-dust ratios and the ratios of far-infrared to radio emission all suggest that many ETGs contain a cool interstellar medium similar to that in late-type galaxies. The mean dust-to-stellar mass ratio for S0s is approximatly a factor of ten less than for early-type spirals and the sizes of the dust sources in the S0s are also much smaller. We show that the difference cannot be explained by either the different bulge-to-disk ratios or environmental effects such as ram-pressure stripping. The wide range in the dust-to-stellar mass ratio for ETGs and the lack of a correlation between dust mass and optical luminosity suggest that much of the dust in the ETGs detected by Herschel has been acquired as the result of gravitational interactions; these interactions are unlikely to have had a major effect on the stellar masses of the ETGs. The Herschel observations tentatively suggest that in the most massive ETGs, the mass of the interstellar medium is unconnected to the evolution of the stellar populations.

It has often been suggested that an alternative to the standard CO/21-cm method for estimating the mass of the interstellar medium (ISM) in a galaxy might be to estimate the mass of the ISM from the continuum dust emission. In this paper, we investigate the potential of this technique using Herschel observations of ten galaxies in the Herschel Reference Survey and in the Herschel Virgo Cluster Survey. We show that the emission detected by Herschel is mostly from dust that has a temperature and emissivity index similar to that of dust in the local ISM in our galaxy, with the temperature generally increasing towards the centre of each galaxy. We calibrate the dust method using the CO and 21-cm observations to provide an independent estimate of the mass of hydrogen in each galaxy, solving the problem of the uncertain `X factor' for the molecular gas by minimizing the dispersion in the ratio of the masses estimated using the two methods. With the calibration for the dust method and the estimate of the X-factor produced in this way, the dispersion in the ratio of the two gas masses is 30%, which gives an upper limit on the fundamental accuracy of the dust method. The calibration we obtain for the dust method is very similar to an independent Herschel measurement for M31 and to the calibration for the Milky Way from Planck measurements.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data that cover the complete survey area (four 4 x 4 deg2 regions). We use these data to measure and compare the global far infrared properties of 78 optically bright galaxies that are selected at 500 \mum and detected in all five far-infrared bands. We show that our measurements and calibration are broadly consistent with previous data obtained by IRAS, ISO, Spitzer and Planck. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 \mum cluster luminosity distributions. These luminosity distributions are not power laws, but peaked, with small numbers of both faint and bright galaxies. We measure a cluster 100-500 micron far-infrared luminosity density of 1.6(7.0) \pm 0.2 x 10^9 Lsun/Mpc3. This compares to a cluster 0.4-2.5 \mum optical luminosity density of 5.0(20.0) x 10^9 Lsun/Mpc3, some 3.2(2.9) times larger than the far-infrared. A typical photon originates from an optical depth of 0.4\pm0.1. Most of our sample galaxies are well fitted by a single modified blackbody (beta=2), leading to a mean dust mass of log Mdust = 7.31 Msun and temperature of 20.0 K. We also derive both stellar and atomic hydrogen masses from which we calculate mean values for the stars:gas(atomic) and gas(atomic): dust mass ratios of 15.1 and 58.2 respectively. Using our derived dust, atomic gas and stellar masses we estimate cluster mass densities of 8.6(27.8) x 10^6, 4.6(13.9) x 10^8, 7.8(29.7) x 10^9 Msun/Mpc3, respectively for dust, atomic gas and stars. These values are higher than those derived for field galaxies by factors of 39(126), 6(18) and 34(129) respectively. In the above luminosity/mass densities are given using the whole sample with values in brackets using just those galaxies that lie between 17 and 23 Mpc. We provide a data table of flux densities in all the Herschel bands for all 78 bright Virgo cluster galaxies.

We present a detailed study of the dust energy balance in the Sombrero galaxy M104. From a full radiative transfer analysis, including scattering, absorption and thermal re-emission, we construct models that can reproduce images at optical/near-infrared wavelengths, the observed stellar SED and the minor axis extinction profiles in the V and R_C band. A standard model, that contains only an old stellar population to heat the dust, underestimates the observations of dust emission at infrared wavelengths by a factor of ~ 3. Supplementing this basic model with a young stellar component of low star formation activity in both the inner disk (SFR ~ 0.21 Msun/yr) and dust ring (SFR ~ 0.05 Msun/yr), we are capable of solving the discrepancy in the dust energy budget of the Sombrero galaxy at wavelengths shortwards of 100 \mum. To account for the increased FIR/submm emission beyond 100 \mum, we propose a additional dust component distributed in quiescent clumps. This model with a clumpy dust structure predicts three-quarters of the total dust content (~ 2.8 x 10^7 Msun) to reside in compact dust clouds with no associated embedded sources. Although the assumption of a clumpy dust structure in the Sombrero galaxy is supported by high-resolution optical data, we cannot rule out the possibility that dust grains with a higher dust emissivity account for part of the discrepancy in the energy budget at submm wavelengths.

We present an updated version of SKIRT, a 3D Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, 2D, LTE model to a fully 3D model that includes NLTE dust emission and a clumpy structure of the dusty ISM. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

Using Herschel data from the Open Time Key Project the Herschel Virgo Cluster Survey (HeViCS), we investigated the relationship between the metallicity gradients expressed by metal abundances in the gas phase as traced by the chemical composition of HII regions, and in the solid phase, as traced by the dust-to-gas mass ratio. We derived the radial gradient of the dust-to-gas mass ratio for all galaxies observed by HeViCS whose metallicity gradients are available in the literature. They are all late type Sbc galaxies, namely NGC4254, NGC4303, NGC4321, and NGC4501. We examined different dependencies on metallicity of the CO-to-H$_2$ conversion factor (\xco), used to transform the $^{12}$CO observations into the amount of molecular hydrogen. We found that in these galaxies the dust-to-gas mass ratio radial profile is extremely sensitive to choice of the \xco\ value, since the molecular gas is the dominant component in the inner parts. We found that for three galaxies of our sample, namely NGC4254, NGC4321, and NGC4501, the slopes of the oxygen and of the dust-to-gas radial gradients agree up to $\sim$0.6-0.7R$_{25}$ using \xco\ values in the range 1/3-1/2 Galactic \xco. For NGC4303 a lower value of \xco$\sim0.1\times$ 10$^{20}$ is necessary. We suggest that such low \xco\ values might be due to a metallicity dependence of \xco (from close to linear for NGC4254, NGC4321, and NGC4501 to superlinear for NGC4303), especially in the radial regions R$_G<$0.6-0.7R$_{25}$ where the molecular gas dominates. On the other hand, the outer regions, where the atomic gas component is dominant, are less affected by the choice of \xco, and thus we cannot put constraints on its value.

The GALEX Ultraviolet Virgo Cluster Survey (GUViCS) is a complete blind survey of the Virgo cluster covering about 40 sq. deg. in the far UV (FUV, lambda_eff=1539A, Delta-lambda=442A) and about 120 sq. deg. in the near UV (NUV, lambda_eff=2316A, Delta-lambda=1060A). The goal of the survey is to study the ultraviolet (UV) properties of galaxies in a rich cluster environment, spanning a wide luminosity range from giants to dwarfs, and regardless of prior knowledge of their star formation activity. The UV data will be combined with those in other bands (optical: NGVS; far-infrared - submm: HeViCS; HI: ALFALFA) and with our multizone chemo-spectrophotometric models of galaxy evolution to make a complete and exhaustive study of the effects of the environment on the evolution of galaxies in high density regions. We present here the scientific objectives of the survey, describing the observing strategy and briefly discussing different data reduction techniques. Using UV data already in-hand for the central 12 sq. deg. we determine the FUV and NUV luminosity functions of the Virgo cluster core for all cluster members and separately for early- and late-type galaxies and compare it to the one obtained in the field and other nearby clusters (Coma, A1367). This analysis shows that the FUV and NUV luminosity functions of the core of the Virgo clusters are flatter (alpha about -1.1) than those determined in Coma and A1367. We discuss the possible origin of this difference

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We observed M87 with PACS and SPIRE as part of the Herschel Virgo Cluster Survey (HeViCS). We compare the new Herschel data with a synchrotron model based on infrared, submm and radio data to investigate the origin of the far-infrared emission. We find that both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87.

Interstellar dust grains efficiently absorb and scatter UV and optical radiation in galaxies, and therefore can significantly affect the apparent structure of spiral galaxies. We discuss the effect of dust attenuation on the observed structural properties of bulges and discs. We also present some first results on modelling the dust content of edge-on spiral galaxies using both optical and Herschel far-infrared data. Both of these results demonstrate the complex interplay of dust and star light in spiral galaxies.

We describe the Herschel Virgo Cluster Survey (HeViCS) and the first data obtained as part of the Science Demonstration Phase (SDP). The data cover a central 4x4 sq deg region of the cluster. We use SPIRE and PACS photometry data to produce 100, 160, 250, 350 and 500 micron luminosity functions (LFs) for optically bright galaxies that are selected at 500 micron and detected in all bands. We compare these LFs with those previously derived using IRAS, BLAST and Herschel-ATLAS data. The Virgo Cluster LFs do not have the large numbers of faint galaxies or examples of very luminous galaxies seen previously in surveys covering less dense environments.

By combining Herschel-SPIRE observations obtained as part of the Herschel Virgo Cluster Survey with 21 cm HI data from the literature, we investigate the role of the cluster environment on the dust content of Virgo spiral galaxies.We show for the first time that the extent of the dust disk is significantly reduced in HI-deficient galaxies, following remarkably well the observed 'truncation' of the HI disk. The ratio of the submillimetre-to- optical diameter correlates with the HI-deficiency, suggesting that the cluster environment is able to strip dust as well as gas. These results provide important insights not only into the evolution of cluster galaxies but also into the metal enrichment of the intra-cluster medium.

Passive early-type galaxies (ETGs) provide an ideal laboratory for studying the interplay between dust formation around evolved stars and its subsequent destruction in a hot gas. Using Spitzer-IRS and Herschel data we compare the dust production rate in the envelopes of evolved AGB stars with a constraint on the total dust mass. Early-type galaxies which appear to be truly passively evolving are not detected by Herschel. We thus derive a distance independent upper limit to the dust grain survival time in the hostile environment of ETGs of < 46 +/- 25 Myr for amorphous silicate grains. This implies that ETGs which are detected at far-infrared wavelengths have acquired a cool dusty medium via interaction. Given likely time-scales for ram-pressure stripping, this also implies that only galaxies with dust in a cool (atomic) medium can release dust into the intra-cluster medium.

We present a resolved dust analysis of three of the largest angular size spiral galaxies, NGC 4501 and NGC 4567/8, in the Herschel Virgo Cluster Survey (HeViCS) Science Demonstration field. Herschel has unprecedented spatial resolution at far-infrared wavelengths and with the PACS and SPIRE instruments samples both sides of the peak in the far infrared spectral energy distribution (SED).We present maps of dust temperature, dust mass, and gas-to-dust ratio, produced by fitting modified black bodies to the SED for each pixel. We find that the distribution of dust temperature in both systems is in the range ~19 - 22 K and peaks away from the centres of the galaxies. The distribution of dust mass in both systems is symmetrical and exhibits a single peak coincident with the galaxy centres. This Letter provides a first insight into the future analysis possible with a large sample of resolved galaxies to be observed by Herschel.

We present the dust properties of a small sample of Virgo cluster dwarf galaxies drawn from the science demonstration phase data set of the Herschel Virgo Cluster Survey. These galaxies have low metallicities (7.8 < 12 + log(O/H) < 8.3) and star-formation rates < 10^{-1} M_{sun}/yr. We measure the spectral energy distribution (SED) from 100 to 500 um and derive dust temperatures and dust masses. The SEDs are fitted by a cool component of temperature T < 20 K, implying dust masses around 10^{5} M_{sun} and dust-to-gas ratios D within the range 10^{-3}-10^{-2}. The completion of the full survey will yield a larger set of galaxies, which will provide more stringent constraints on the dust content of star-forming dwarf galaxies.

The origin of the far-infrared emission from the nearby radio galaxy M87 remains a matter of debate. Some studies find evidence of a far-infrared excess due to thermal dust emission, whereas others propose that the far-infrared emission can be explained by synchrotron emission without the need for an additional dust emission component. We present Herschel PACS and SPIRE observations of M87, taken as part of the science demonstration phase observations of the Herschel Virgo Cluster Survey. We compare these data with a synchrotron model based on mid-infrared, far-infrared, submm and radio data from the literature to investigate the origin of the far-infrared emission. Both the integrated SED and the Herschel surface brightness maps are adequately explained by synchrotron emission. At odds with previous claims, we find no evidence of a diffuse dust component in M87, which is not unexpected in the harsh X-ray environment of this radio galaxy sitting at the core of the Virgo Cluster.

We use the Science Demonstration Phase data of the Herschel Virgo Cluster Survey to search for dust emission of early-type dwarf galaxies in the central regions of the Virgo Cluster as an alternative way of identifying the interstellar medium.We present the first possible far-infrared detection of cluster early-type dwarf galaxies: VCC781 and VCC951 are detected at the 10 sigma level in the SPIRE 250 micron image. Both detected galaxies have dust masses of the order of 10^5 Msun and average dust temperatures ~20K. The detection rate (less than 1%) is quite high compared to the 1.7% detection rate for Hi emission, considering that dwarfs in the central regions are more Hi deficient. We conclude that the removal of interstellar dust from dwarf galaxies resulting from ram pressure stripping, harassment, or tidal effects must be as efficient as the removal of interstellar gas.

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.