Source author record

Ilse de Looze

Ilse de Looze appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2026arXiv

REBELS-MOSFIRE: Weak CIII] Emission is Typical Among Extremely UV-bright, Massive Galaxies at $z\sim7$

We present Keck/MOSFIRE H-band spectroscopic measurements covering the [CIII]1907, CIII]1909 doublet for a sample of 8 z~7 spectroscopically-confirmed star-forming galaxies drawn from the Reionization Era Bright Emission Line Survey (REBELS). This REBELS-MOSFIRE sample is notable for its bright median UV luminosity (Muv=-22.5 AB) and large median stellar mass (log(Mstar/Msun)=9.2). Although three sources show tentative evidence of a CIII] detection, we obtain no confident detections for any of the 8 REBELS-MOSFIRE sources. The median [CIII]1907+CIII]1909 3-sigma upper limit in equivalent width (EW) for the REBELS-MOSFIRE sample is 6.5 AA, and a stack of their H-band MOSFIRE spectra yields a non-detection with an associated 3-sigma upper limit of 2.6 AA. These upper limits fall significantly below the CIII] EW measured in a composite spectrum of representative z~7 star-forming galaxies, as well as those measured for notable early star-forming galaxies such as GN-z11, GHZ2, GS-z12, and RXCJ2248-ID. The lack of strong CIII] emission can be understood within the context of the stellar populations of the REBELS galaxies, as well as the ionization conditions and gas-phase metallicity implied by rest-frame optical spectroscopic properties ([OIII]+Hb EWs, and [OIII]5007/[OII]3727 and [NeIII]3869/[OII]3727 line ratios). The REBELS-MOSFIRE sample represents the higher-mass, higher-metallicity, lower-excitation tail of the z~7 galaxy population, whose ionizing properties must be fully characterized to constrain the role of star-forming galaxies during cosmic reionization.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Dust masses for a large sample of core-collapse supernovae from optical emission line asymmetries: dust formation on 30-year timescales

Modelling the red-blue asymmetries seen in the broad emission lines of core-collapse supernovae (CCSNe) is a powerful technique to quantify total dust mass formed in the ejecta at late times ($>5$ years after outburst) when ejecta dust temperatures become too low to be detected by mid-IR instruments. Following our success in using the Monte Carlo radiative transfer code DAMOCLES to measure the dust mass evolution in SN~1987A and other CCSNe, we present the most comprehensive sample of dust mass measurements yet made with DAMOCLES, for CCSNe aged between four and sixty years after outburst. Our sample comprises of multi-epoch late-time optical spectra taken with the Gemini GMOS and VLT X-Shooter spectrographs, supplemented by archival spectra. For the fourteen CCSNe that we have modelled, we confirm a dust mass growth with time that can be fit by a sigmoid curve which is found to saturate beyond an age of $\sim30$ years, at a mass of 0.23$^{+0.17}_{-0.12}$ M$_\odot$. An expanded sample including dust masses found in the literature for a further eleven CCSNe and six CCSN remnants, the dust mass at saturation is found to be 0.42$^{+0.09}_{-0.05}$~M$_\odot$. Uncertainty limits for our dust masses were determined from a Bayesian analysis using the affine invariant Markov Chain Monte Carlo ensemble sampler emcee with DAMOCLES. The best-fitting line profile models for our sample all required grain radii between 0.1 and 0.5 $μ$m. Our results are consistent with CCSNe forming enough dust in their ejecta to significantly contribute to the dust budget of the Universe.

preprint2022arXiv

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates.

preprint2022arXiv

The First Maps of $κ_{d}$ -- the Dust Mass Absorption Coefficient -- in Nearby Galaxies, with DustPedia

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 M74 (NGC 628) and M83 (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 M74 and M83, 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 $μ$m spanning the range 0.11-0.25 ${\rm m^{2}\,kg^{-1}}$ in M74, and 0.15-0.80 ${\rm m^{2}\,kg^{-1}}$ in M83. 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.

preprint2021arXiv

AlFoCS + F3D II: unexpectedly low gas-to-dust ratios in the Fornax galaxy cluster

We combine observations from ALMA, ATCA, MUSE, andHerschel to study gas-to-dust ratios in 15 Fornax cluster galaxies detected in the FIR/sub-mm by Herschel and observed by ALMA as part of the ALMA Fornax Cluster Survey (AlFoCS). The sample spans a stellar mass range of 8.3 $\leq$ log (M$_*$ / M$_\odot$) $\leq$ 11.16, and a variety of morphological types. We use gas-phase metallicities derived from MUSE observations (from the Fornax3D survey) to study these ratios as a function of metallicity, and to study dust-to-metal ratios, in a sub-sample of nine galaxies. We find that gas-to-dust ratios in Fornax galaxies are systematically lower than those in field galaxies at fixed stellar mass/metallicity. This implies that a relatively large fraction of the metals in these Fornax systems is locked up in dust, which is possibly due to altered chemical evolution as a result of the dense environment. The low ratios are not only driven by HI deficiencies, but H$_2$-to-dust ratios are also significantly decreased. This is different in the Virgo cluster, where low gas-to-dust ratios inside the virial radius are driven by low HI-to-dust ratios, while H$_2$-to-dust ratios are increased. Resolved observations of NGC1436 show a radial increase in H$_2$-to-dust ratio, and show that low ratios are present throughout the disc. We propose various explanations for the low H$_2$-to-dust ratios in the Fornax cluster, including the more efficient stripping of H$_2$ compared to dust, more efficient enrichment of dust in the star formation process, and altered ISM physics in the cluster environment.

preprint2021arXiv

The Case for Thermalization as a Contributor to the [CII] Deficit

The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission from 18 Local Volume galaxies has been split based on conditions within the interstellar medium where it originated. This is completed using the Key Insights in Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) surveys and the wide-range of wavelength information, from UV to far-infrared emission lines, available for a selection of star-forming regions within these samples. By comparing these subdivided [CII] emissions to isolated infrared emission and other properties, we find that the thermalization (collisional de-excitation) of the [CII] line in HII regions plays a significant role in the deficit observed in our sample.

preprint2020arXiv

A Complete Catalogue of Dusty Supernova Remnants

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (360 degrees in longitude and b = +/- 1 deg ) at 70 - 500 micron with Herschel. We detect dust signatures in 39 SNRs out of 190, made up of 13 core-collapse supernovae (CCSNe), including 4 Pulsar Wind Nebulae (PWNe), and 2 Type Ia SNe. A further 24 FIR detected SNRs have unknown types. We confirm the FIR detection of ejecta dust within G350.1-0.3, adding to the known sample of ~10 SNRs containing ejecta dust. We discover dust features at the location of a radio core at the centre of G351.2+0.1, indicating FIR emission coincident with a possible Crab-like compact object, with dust temperature and mass of Td = 45.8 K and Md = 0.18 solar mass, similar to the PWN G54.1+0.3. We show that the detection rate is higher among young SNRs. We produce dust temperature maps of 11 SNRs and mass maps of those with distance estimates, finding dust at temperatures 15 < Td < 40 K. If the dust is heated by shock interactions the shocked gas must be relatively cool and/or have a low density to explain the observed low grain temperatures.

preprint2020arXiv

AlFoCS + Fornax3D: resolved star formation in the Fornax cluster with ALMA and MUSE

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the 'uncertainty principle for star formation' (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (<10 Myr) in this galaxy. Galaxy-averaged depletion times are compared with other parameters such as stellar mass and cluster-centric distance. We find that the star formation relation in the Fornax cluster is close to those from Kennicutt (1998) and Bigiel et al. (2008}), but overlaps mostly with the shortest depletion times predicted by Bigiel et al. (2008). This slight decrease in depletion time is mostly driven by dwarf galaxies with disturbed molecular gas reservoirs close to the virial radius. In FCC90, a dwarf galaxy with a molecular gas tail, we find that depletion times are a factor >~10 higher in its tail than in its stellar body.

preprint2020arXiv

ALMA Observations of the Molecular Clouds in NGC 625

We present the highest resolution (1") $^{12}$CO(1-0) observations of molecular gas in the dwarf starburst galaxy NGC 625, to date, obtained with ALMA. Molecular gas is distributed in discrete clouds within an area of $0.4$ kpc$^2$ and does not have well-ordered large-scale motions. We measure a total molecular mass in NGC 625 of $5.3\times 10^6 M_\odot$, assuming a Milky Way CO-to-H$_2$ conversion factor. We use the CPROPS package to identify molecular clouds and measure their properties. The 19 resolved CO clouds have a median radius of 20 pc, a median line width 2.5 km s$^{-1}$, and a median surface density of 169 $M_\odot$pc$^{-2}$. Larson scaling relations suggest that molecular clouds in NGC 625 are mostly in virial equilibrium. Comparison of our high-resolution CO observations with a star formation rate map, inferred from ancillary optical observations observations, suggests that about 40% of the molecular clouds coincide with the brightest HII regions. These bright HII regions have a range of molecular gas depletion timescales, all within a factor of $\sim 3$ of the global depletion time in NGC 625 of 106-134 Myr. The highest surface density molecular clouds towards the southwest of the galaxy, in a region we call the Butterfly, do not show strong star formation activity and suggest a depletion time scale longer than 5 Gyr.

preprint2020arXiv

High-resolution, 3D radiative transfer modelling III. The DustPedia barred galaxies

Context: Dust in late-type galaxies in the local Universe is responsible for absorbing approximately one third of the energy emitted by stars. It is often assumed that dust heating is mainly attributable to the absorption of UV and optical photons emitted by the youngest (<= 100 Myr) stars. Consequently, thermal re-emission by dust at FIR wavelengths is often linked to the star-formation activity of a galaxy. However, several studies argue that the contribution to dust heating by much older stars might be more significant. Advances in radiation transfer (RT) simulations finally allow us to actually quantify the heating mechanisms of diffuse dust by the stellar radiation field. Aims: As one of the main goals in the DustPedia project, we have constructed detailed 3D stellar and dust RT models for nearby galaxies. We analyse the contribution of the different stellar populations to the dust heating in four face-on barred galaxies: NGC1365, M83, M95, and M100. We aim to quantify the fraction directly related to young stars, both globally and on local scales, and to assess the influence of the bar on the heating fraction. Results: We derive global attenuation laws for each galaxy and confirm that galaxies of high sSFR have shallower attenuation curves and weaker UV bumps. On average, 36.5% of the bolometric luminosity is absorbed by dust. We report a clear effect of the bar structure on the radial profiles of the dust-heating fraction by the young stars, and the dust temperature. We find that the young stars are the main contributors to the dust heating, donating, on average ~59% of their luminosity to this purpose throughout the galaxy. This dust-heating fraction drops to ~53% in the bar region and ~38% in the bulge region where the old stars are the dominant contributors to the dust heating. We also find a strong link between the heating fraction by the young stars and the sSFR.

preprint2020arXiv

High-resolution, 3D radiative transfer modelling. II. The early-type spiral galaxy M81

Interstellar dust absorbs stellar light very efficiently and thus shapes the energetic output of galaxies. Studying the impact of different stellar populations on the dust heating remains hard because it requires decoupling the relative geometry of stars and dust, and involves complex processes as scattering and non-local dust heating. We aim to constrain the relative distribution of dust and stellar populations in the spiral galaxy M81 and create a realistic model of the radiation field that describes the observations. Investigating the dust-starlight interaction on local scales, we want to quantify the contribution of young and old stellar populations to the dust heating. We aim to standardise the setup and model selection of such inverse radiative transfer simulations so this can be used for comparable modelling of other nearby galaxies. We present a semi-automated radiative transfer modelling pipeline that implements the necessary steps such as the geometric model construction and the normalisation of the components through an optimisation routine. We use the Monte Carlo radiative transfer code SKIRT to calculate a self-consistent, panchromatic model of the interstellar radiation field. By looking at different stellar populations independently, we can quantify to what extent different stellar age populations contribute to the dust heating. Our method takes into account the effects of non-local heating. We obtain a realistic 3D radiative transfer model of the face-on galaxy M81. We find that only 50.2\% of the dust heating can be attributed to young stellar populations. We confirm a tight correlation between the specific star formation rate and the heating fraction by young stellar populations, both in sky projection and in 3D, also found for radiative transfer models of M31 and M51. We conclude that... (abridged)

preprint2020arXiv

Nonparametric galaxy morphology from UV to submm wavelengths

We present the first nonparametric morphological analysis of a set of spiral galaxies from UV to submm wavelengths. Our study is based on high-quality multi-wavelength imaging for nine well-resolved spiral galaxies from the DustPedia database, combined with nonparametric morphology indicators calculated in a consistent way using the {\tt{StatMorph}} package. We measure the half-light radius, the concentration index, the asymmetry index, the smoothness index, the Gini coefficient and the $M_{20}$ indicator in various wavebands from UV to submm wavelengths, as well as in stellar mass, dust mass and star formation rate maps. We find that the interstellar dust in galaxies is distributed in a more extended, less centrally concentrated, more asymmetric, and more clumpy way than the stars. This is particularly evident when comparing morphological indicators based on the stellar mass and dust mass maps. This should serve as a warning sign against treating the dust in galaxies as a simple smooth component. We argue that the nonparametric galaxy morphology of galaxies from UV to submm wavelengths is an interesting test for cosmological hydrodynamics simulations.

preprint2020arXiv

Reproducing the Universe: a comparison between the EAGLE simulations and the nearby DustPedia galaxy sample

We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models.

preprint2019arXiv

The ISM scaling relations in DustPedia late-type galaxies: A benchmark study for the Local Universe

The purpose of this work is the characterization of the main scaling relations between all the ISM components (dust, atomic/molecular/total gas), gas-phase metallicity, and other galaxy properties, such as Mstar and galaxy morphology, for late-type galaxies in the Local Universe. This study is 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 from T = 1 to 10, Mstar from 6 x 10^7 to 3 x 10^11 Msun, SFR from 6 x 10^(-4) to 60 Msun/yr, and 12 + log(O/H) from 8 to 9.5. The scaling relations involving the molecular gas are studied by assuming both a constant and a metallicity-dependent XCO. The analysis has been performed by means of the survival analysis technique. 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 the molecular one. We characterize important mass ratios such as gas fraction, molecular-to-atomic gas mass ratio, dust-to-total gas mass ratio (DGR), and dust-to-stellar mass ratio. Only the assumption of a metallicity-dependent XCO reproduces the expected decrease of the DGR with increasing morphological stage and decreasing gas-phase metallicity, with a slope of about 1. DGR, gas-phase metallicity, and the dust-to-stellar mass ratio are, for our galaxy sample, directly linked to the galaxy morphology. The molecular-to-atomic gas mass ratio and the DGR show a positive correlation for low molecular gas fractions, but for molecular gas rich galaxies this trend breaks down. This trend has never been found previously, to our knowledge. It provides new constraints for theoretical models of galaxy evolution and a reference for high-redshift studies.

preprint2019arXiv

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample.

preprint2018arXiv

A Catalogue of Galactic Supernova Remnants in the far-infrared: revealing ejecta dust in pulsar wind nebulae

We search for far-infrared (FIR) counterparts of known supernova remnants (SNRs) in the Galactic plane (10 deg <| l |< 60 deg) at 70-500 micron using the Herschel Infrared Galactic Plane Survey (Hi-GAL). Of 71 sources studied, we find that 29 (41 per cent) SNRs have a clear FIR detection of dust emission associated with the SNR. Dust from 8 of these is in the central region, and 4 indicate pulsar wind nebulae (PWNe) heated ejecta dust. A further 23 have dust emission in the outer shell structures which is potentially related to swept up material. Many Galactic SNe have dust signatures but we are biased towards detecting ejecta dust in young remnants and those with a heating source (shock or PWN). We estimate the dust temperature and mass contained within three PWNe, G11.2-0.3, G21.5-0.9, and G29.7-0.3 using modified blackbody fits. To more rigorously analyse the dust properties at various temperatures and dust emissivity index beta, we use point process mapping (PPMAP). We find significant quantities of cool dust (at 20-40 K) with dust masses of Md = 0.34 +/- 0.14 solar mass, Md = 0.29 +/- 0.08 solar mass, and Md = 0.51 +/- 0.13 solar mass for G11.2-0.3, G21.5-0.9, and G29.7-0.3 respectively. We derive the dust emissivity index for the PWN ejecta dust in G21.5-0.3 to be beta = 1.4 +/- 0.5 compared to dust in the surrounding medium where beta = 1.8 +/- 0.1.

preprint2016arXiv

Far-reaching Dust Distribution in Galaxy Disks

In most studies of dust in galaxies, dust is only detected from its emission to approximately the optical radius of the galaxy. By combining the signal of 110 spiral galaxies observed as part of the Herschel Reference Survey, we are able to improve our sensitivity by an order-of-magnitude over that for a single object. Here we report the direct detection of dust from its emission that extends out to at least twice the optical radius. We find that the distribution of dust is consistent with an exponential at all radii with a gradient of ~-1.7 dex R$_{25}^{-1}$. Our dust temperature declines linearly from ~25 K in the centre to 15 K at R$_{25}$ from where it remains constant out to ~2.0 R$_{25}$. The surface-density of dust declines with radius at a similar rate to the surface-density of stars but more slowly than the surface-density of the star-formation rate. Studies based on dust extinction and reddening of high-redshift quasars have concluded that there are substantial amounts of dust in intergalactic space. By combining our results with the number counts and angular correlation function from the SDSS, we show that with Milky Way type dust we can explain the reddening of the quasars by the dust within galactic disks alone. Given the uncertainties in the properties of any intergalactic dust, we cannot rule out its existence, but our results show that statistical investigations of the dust in galactic halos that use the reddening of high-redshift objects must take account of the dust in galactic disks.

preprint2016arXiv

HERschel Observations of Edge-on Spirals (HEROES). III. Dust energy balance study of IC 2531

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.

preprint2016arXiv

The Herschel Virgo Cluster Survey. XIX. Physical properties of low luminosity FIR sources at $z <$ 0.5

The Star formation rate (SFR) is a crucial parameter to investigate galaxy evolution. At low redshift the cosmic SFR density declines smoothly, and massive active galaxies become passive, reducing their star formation activity. This implies that the bulk of the SFR density at low redshift is mainly driven by low mass objects. We investigate the properties of a sample of low luminosity Far-Infrared (FIR) sources selected at 250 microns from Pappalardo et al. (2015). We have collected data from Ultraviolet to FIR to perform a multi-wavelengths analysis. The main goal is to investigate the correlation between SFR, stellar mass, and dust mass for a galaxy population with a wide range in dust content and stellar mass, including the low mass regime that most probably dominates the SFR density at low z. We define a main sample of ~800 sources with full Spectral Energy Distribution (SED) coverage between 0.15 < lambda < 500 microns and an extended sample with ~5000 sources in which we remove the constraints on the Ultraviolet and Near-Infrared bands. We analyze both samples with two different SED fitting methods: MAGPHYS and CIGALE. In the SFR versus stellar mass plane our samples occupy a region included between local spirals and higher redshift star forming galaxies. The galaxies subsample with the higher masses (M* > 3e10 Msol) does not lie on the main sequence, but shows a small offset, as a consequence of the decreased star formation. Low mass galaxies (M* < 1e10 Msol) settle in the main sequence with SFR and stellar mass consistent with local spirals. Deep Herschel data allow the identification of a mixed galaxy population, with galaxies still in an assembly phase, or galaxies at the beginning of their passive evolution. We find that the dust luminosity is the parameter that discriminates these two galaxy populations.

preprint2016arXiv

The interstellar medium in Andromeda's dwarf spheroidal galaxies: I. Content and origin of the interstellar dust

Dwarf spheroidal galaxies are among the most numerous galaxy population in the Universe, but their main formation and evolution channels are still not well understood. The three dwarf spheroidal satellites (NGC147, NGC185, and NGC205) of the Andromeda galaxy are characterised by very different interstellar medium (ISM) properties, which might suggest them being at different galaxy evolutionary stages. While the dust content of NGC205 has been studied in detail by De Looze et al. (2012), we present new Herschel dust continuum observations of NGC147 and NGC185. The non-detection of NGC147 in Herschel SPIRE maps puts a strong constraint on its dust mass (< 128 Msun). For NGC185, we derive a total dust mass M_d = 5.1 x 10^3 Msun, which is a factor of ~2-3 higher than that derived from ISO and Spitzer observations and confirms the need for longer wavelength observations to trace more massive cold dust reservoirs. We, furthermore, estimate the dust production by asymptotic giant branch (AGB) stars and supernovae (SNe). For NGC147, the upper limit on the dust mass is consistent with expectations of the material injected by the evolved stellar population. In NGC185 and NGC205, the observed dust content is one order of magnitude higher compared to the estimated dust production by AGBs and SNe. Efficient grain growth, and potentially longer dust survival times (3-6 Gyr) are required to account for their current dust content. Our study confirms the importance of grain growth in the gas phase to account for the current dust reservoir in galaxies.

preprint2016arXiv

The interstellar medium in Andromeda's dwarf spheroidal galaxies: II. Multi-phase gas content and ISM conditions

We make an inventory of the interstellar medium material in three low-metallicity dwarf spheroidal galaxies of the Local Group (NGC147, NGC185 and NGC205). Ancillary HI, CO, Spitzer IRS spectra, Hα and X-ray observations are combined to trace the atomic, cold and warm molecular, ionised and hot gas phases. We present new Nobeyama CO(1-0) observations and Herschel SPIRE FTS [CI] observations of NGC205 to revise its molecular gas content. We derive total gas masses of M_gas = 1.9-5.5x10^5 Msun for NGC185 and M_gas = 8.6-25.0x10^5 Msun for NGC205. Non-detections combine to an upper limit on the gas mass of M_gas =< 0.3-2.2x10^5 Msun for NGC147. The observed gas reservoirs are significantly lower compared to the expected gas masses based on a simple closed-box model that accounts for the gas mass returned by planetary nebulae and supernovae. The gas-to-dust mass ratios GDR~37-107 and GDR~48-139 are also considerably lower compared to the expected GDR~370 and GDR~520 for the low metal abundances in NGC 185 (0.36 Zsun) and NGC205 (0.25 Zsun), respectively. To simultaneously account for the gas deficiency and low gas-to-dust ratios, we require an efficient removal of a large gas fraction and a longer dust survival time (~1.6 Gyr). We believe that efficient galactic winds (combined with heating of gas to sufficiently high temperatures in order for it to escape from the galaxy) and/or environmental interactions with neighbouring galaxies are responsible for the gas removal from NGC147, NGC185 and NGC205.

preprint2016arXiv

The Spatially Resolved [CII] Cooling Line Deficit in Galaxies

We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas.

preprint2015arXiv

Dust energy balance study of two edge-on spiral galaxies in the Herschel-ATLAS survey

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

preprint2015arXiv

Far-infrared observations of an unbiased sample of gamma-ray burst host galaxies

Gamma-ray bursts (GRBs) are the most energetic phenomena in the Universe; believed to result from the collapse and subsequent explosion of massive stars. Even though it has profound consequences for our understanding of their nature and selection biases, little is known about the dust properties of the galaxies hosting GRBs. We present analysis of the far-infrared properties of an unbiased sample of 20 \textit{BeppoSAX} and \textit{Swift} GRB host galaxies (at an average redshift of $z\,=\,3.1$) located in the {\it Herschel} Astrophysical Terahertz Large Area Survey, the {\it Herschel} Virgo Cluster Survey, the {\it Herschel} Fornax Cluster Survey, the {\it Herschel} Stripe 82 Survey and the {\it Herschel} Multi-tiered Extragalactic Survey, totalling $880$ deg$^2$, or $\sim 3$\% of the sky in total. Our sample selection is serendipitous, based only on whether the X-ray position of a GRB lies within a large-scale {\it Herschel} survey -- therefore our sample can be considered completely unbiased. Using deep data at wavelengths of 100\,--\,500$\,μ$m, we tentatively detected 1 out of 20 GRB hosts located in these fields. We constrain their dust masses and star formation rates (SFRs), and discuss these in the context of recent measurements of submillimetre galaxies and ultraluminous infrared galaxies. The average far-infrared flux of our sample gives an upper limit on SFR of $<114\,{\rm M}\odot\,\mbox{yr}^{-1}$. The detection rate of GRB hosts is consistent with that predicted assuming that GRBs trace the cosmic SFR density in an unbiased way, i.e. that the fraction of GRB hosts with $\mbox{SFR}>500\,{\rm M}\odot\,\mbox{yr}^{-1}$ is consistent with the contribution of such luminous galaxies to the cosmic star formation density.

preprint2014arXiv

An extremely low gas-to-dust ratio in the dust-lane lenticular galaxy NGC 5485

Evidence is mounting that a significant fraction of the early-type galaxy population contains substantial reservoirs of cold interstellar gas and dust. We investigate the gas and dust in NGC 5485, an early-type galaxy with a prominent minor-axis dust lane. Using new Herschel PACS and SPIRE imaging data, we detect 3.8 x 10^6 Msun of cool interstellar dust in NGC 5485, which is in stark contrast with the non-detection of the galaxy in sensitive HI and CO observations from the ATLAS3D consortium. The resulting gas-to-dust ratio upper limit is Mgas/Md < 14.5, almost an order of magnitude lower than the canonical value for the Milky Way. We scrutinize the reliability of the dust, atomic gas and molecular gas mass estimates, but these do not show systematic uncertainties that can explain the extreme gas-to-dust ratio. Also a warm or hot ionized gas medium does not offer an explanation. A possible scenario could be that NGC 5485 merged with an SMC-type metal-poor galaxy with a substantial CO-dark molecular gas component and that the bulk of atomic gas was lost during the interaction, but it remains to be investigated whether such a scenario is possible.

preprint2014arXiv

High-resolution, 3D radiative transfer modeling : I. The grand-design spiral galaxy M51

Context: Dust reprocesses about half of the stellar radiation in galaxies. The thermal re-emission by dust of absorbed energy is considered driven merely by young stars and, consequently, often applied to trace the star formation rate in galaxies. Recent studies have argued that the old stellar population might anticipate a non-negligible fraction of the radiative dust heating. Aims: In this work, we aim to analyze the contribution of young (< 100 Myr) and old (~ 10 Gyr) stellar populations to radiative dust heating processes in the nearby grand-design spiral galaxy M51 using radiative transfer modeling. High-resolution 3D radiative transfer (RT) models are required to describe the complex morphologies of asymmetric spiral arms and clumpy star-forming regions and model the propagation of light through a dusty medium. Methods: In this paper, we present a new technique developed to model the radiative transfer effects in nearby face-on galaxies. We construct a high-resolution 3D radiative transfer model with the Monte-Carlo code SKIRT accounting for the absorption, scattering and non-local thermal equilibrium (NLTE) emission of dust in M51. The 3D distribution of stars is derived from the 2D morphology observed in the IRAC 3.6 μm, GALEX FUV, Hα and MIPS 24 μm wavebands, assuming an exponential vertical distribution with an appropriate scale height. The dust geometry is constrained through the far-ultraviolet (FUV) attenuation, which is derived from the observed total-infrared-to-far-ultraviolet luminosity ratio. The stellar luminosity, star formation rate and dust mass have been scaled to reproduce the observed stellar spectral energy distribution (SED), FUV attenuation and infrared SED. (abridged)

preprint2014arXiv

Spatially resolved physical conditions of molecular gas and potential star formation tracers in M83, revealed by the Herschel SPIRE FTS

Since the launch of the Herschel Space Observatory, our understanding about the photo-dissociation regions (PDR) has taken a step forward. In the bandwidth of the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) on board Herschel, ten CO rotational transitions, including J=4-3 to J=13-12, and three fine structure lines, including [CI] 609, [CI] 370, and [NII] 250 micron, are covered. In this paper, we present our findings from the FTS observations at the nuclear region of M83, based on the spatially resolved physical parameters derived from the CO spectral line energy distribution (SLED) map and the comparisons with the dust properties and star-formation tracers. We discuss (1) the potential of using [NII] 250 and [CI] 370 micron as star-formation tracers; (2) the reliability of tracing molecular gas with CO; (3) the excitation mechanisms of warm CO; (4) the possibility of studying stellar feedback by tracing the thermal pressure of molecular gas in the nuclear region of M83.

preprint2014arXiv

The applicability of FIR fine-structure lines as Star Formation Rate tracers over wide ranges of metallicities and galaxy types

We analyze the applicability of far-infrared fine-structure lines [CII] 158 micron, [OI] 63 micron and [OIII] 88 micron to reliably trace the star formation rate (SFR) in a sample of low-metallicity dwarf galaxies from the Herschel Dwarf Galaxy Survey and compare with a broad sample of galaxies of various types and metallicities in the literature. We study the trends and scatter in the relation between the SFR (as traced by GALEX FUV and MIPS 24 micron) and far-infrared line emission, on spatially resolved and global galaxy scales, in dwarf galaxies. We assemble far-infrared line measurements from the literature and infer whether the far-infrared lines can probe the SFR (as traced by the total-infrared luminosity) in a variety of galaxy populations. In metal-poor dwarfs, the [OI] and [OIII] lines show the strongest correlation with the SFR with an uncertainty on the SFR estimates better than a factor of 2, while the link between [CII] emission and the SFR is more dispersed (uncertainty factor of 2.6). The increased scatter in the SFR-L([CII]) relation towards low metal abundances, warm dust temperatures, large filling factors of diffuse, highly ionized gas suggests that other cooling lines start to dominate depending on the density and ionization state of the gas. For the literature sample, we evaluate the correlations for a number of different galaxy populations. The [CII] and [OI] lines are considered to be reliable SFR tracers in starburst galaxies, recovering the star formation activity within an uncertainty of factor 2. [Abridged]

preprint2014arXiv

The distribution of interstellar dust in CALIFA edge-on galaxies via oligochromatic radiative transfer fitting

We investigate the amount and spatial distribution of interstellar dust in edge-on spiral galaxies, using detailed radiative transfer modeling of a homogeneous sample of 12 galaxies selected from the CALIFA survey. Our automated fitting routine, FitSKIRT, was first validated against artificial data. This is done by simultaneously reproducing the SDSS $g$-, $r$-, $i$- and $z$-band observations of a toy model in order to combine the information present in the different bands. We show that this combined, oligochromatic fitting, has clear advantages over standard monochromatic fitting especially regarding constraints on the dust properties. We model all galaxies in our sample using a three-component model, consisting of a double exponential disc to describe the stellar and dust discs and using a Sérsic profile to describe the central bulge. The full model contains 19 free parameters, and we are able to constrain all these parameters to a satisfactory level of accuracy without human intervention or strong boundary conditions. Apart from two galaxies, the entire sample can be accurately reproduced by our model. We find that the dust disc is about 75% more extended but only half as high as the stellar disc. The average face-on optical depth in the V-band is $0.76$ and the spread of $0.60$ within our sample is quite substantial, which indicates that some spiral galaxies are relatively opaque even when seen face-on.

preprint2014arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA). VI. The distribution and properties of molecular cloud associations in M31

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.

preprint2014arXiv

The Herschel Virgo Cluster Survey. XVII. SPIRE point-source catalogs and number counts

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.

preprint2013arXiv

Herschel Exploitation of Local Galaxy Andromeda (HELGA) III: The Star Formation Law in M31

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.

preprint2013arXiv

Herschel SPIRE-FTS Observations of Excited CO and [CI] in the Antennae (NGC 4038/39): Warm and Cold Molecular Gas

We present Herschel SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205μm$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63μm$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $α_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$ in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63μm$ to be $\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength $G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.

preprint2013arXiv

The Herschel Virgo Cluster Survey - XIII. Dust in early-type galaxies

Aims. We study the dust content of a large optical input sample of 910 early-type galaxies (ETG) in the Virgo cluster, extending also to the dwarf ETG, and examine the results in relation with those on the other cold ISM components. Methods. We searched for far-infrared emission in all galaxies of the input sample using the 250 micron image of the Herschel Virgo Cluster Survey (HeViCS). This image covers a large fraction of the cluster. For the detected ETG we measured fluxes in 5 bands from 100 to 500 micron, and estimated the dust mass and temperature with modified black-body fits. Results. Dust is detected above the completeness limit of 25.4 mJy at 250 micron in 46 ETG, 43 of which are in the optically complete part of the input sample. In addition dust is present at fainter levels in another 6 ETG. We detect dust in the 4 ETG with synchrotron emission, including M 87. Dust appears to be much more concentrated than stars and more luminous ETG have higher dust temperatures. Dust detection rates down to the 25.4 mJy limit are 17% for ellipticals, about 40% for lenticulars (S0 + S0a) and around 3% for dwarf ETG. Dust mass does not correlate clearly with stellar mass and is often much more than that expected for a passive galaxy in a closed-box model. The dust-to-stars mass ratio anticorrelates with galaxy luminosity, and for some dwarf ETG reaches values as high as for dusty late-type galaxies. In the Virgo cluster slow rotators appear more likely to contain dust than fast ones. Comparing the dust results with those on HI from ALFALFA, there are only 8 ETG detected both in dust and in HI in the HeViCS area; 39 have dust but only an upper limit on HI, and 8 have HI but only an upper limit on dust. The locations of these galaxies in the cluster are different, with the dusty ETG concentrated in the densest regions, while the HI rich ETG are at the periphery.

preprint2013arXiv

The Herschel Virgo Cluster Survey XIV: transition-type dwarf galaxies in the Virgo cluster

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)

preprint2012arXiv

Can dust emission be used to map the interstellar medium in high-redshift galaxies? Results from the Herschel Reference Survey

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.

preprint2012arXiv

The dust energy balance in the edge-on spiral galaxy NGC 4565

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.

preprint2012arXiv

The Herschel Virgo Cluster Survey: X.The relationship between cold dust and molecular gas content in Virgo spirals

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.

preprint2012arXiv

The Herschel Virgo Cluster Survey. XI. Environmental effects on molecular gas and dust in spiral disks

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 μ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.

preprint2011arXiv

Efficient 3D NLTE dust radiative transfer with SKIRT

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.

preprint2011arXiv

Panchromatic radiative transfer modeling of stars and dust in the Sombrero galaxy

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.

preprint2011arXiv

The Herschel Virgo Cluster Survey. IX. Dust-to-gas mass ratio and metallicity gradients in four Virgo spiral galaxies

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.

preprint2011arXiv

The reliability of [CII] as a star formation rate indicator

The [CII] 157.74 μm line is an important coolant for the neutral interstellar gas. Since [CII] is the brightest spectral line for most galaxies, it is a potentially powerful tracer of star formation activity. In this paper we present a calibration of the star formation rate as a function of the [CII] luminosity for a sample of 24 star-forming galaxies in the nearby universe. This sample includes objects classified as HII regions or LINERs, but omits all Seyfert galaxies with a significant AGN contribution to the mid-infrared photometry. In order to calibrate the SFR against the line luminosity, we rely on both GALEX FUV data, which is an ideal tracer of the unobscured star formation, and MIPS 24 μm, to probe the dust-enshrouded fraction of star formation. In case of normal star-forming galaxies, the [CII] luminosity correlates well with the star formation rate. However, the extension of this relation to more quiescent (Halpha EW <= 10 A) or ultraluminous galaxies should be handled with caution, since these objects show a non-linearity in the L([CII])-to-L(FIR) ratio as a function of L(FIR) (and thus, their star formation activity). We provide two possible explanations for the origin of the tight correlation between the [CII] emission and the star formation activity on a global galaxy-scale. A first interpretation could be that the [CII] emission from PDRs arises from the immediate surroundings of star-forming regions. Since PDRs are neutral regions of warm dense gas at the boundaries between HII regions and molecular clouds and they provide the bulk of [CII] emission in most galaxies, we believe that a more or less constant contribution from these outer layers of photon-dominated molecular clumps to the [CII] emission provides a straightforward explanation for this close link between the [CII] luminosity and SFR. Alternatively,... (abridged)

preprint2010arXiv

The complex interplay of dust and star light in spiral galaxy discs

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.