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.

M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution.

Dust emission, an important diagnostic of star formation and ISM mass throughout the Universe, can be powered by sources unrelated to ongoing star formation. In the framework of the DustPedia project we have set out to disentangle the radiation of the ongoing star formation from that of the older stellar populations. This is done through detailed, 3D radiative transfer simulations of face-on spiral galaxies. In this particular study, we focus on NGC 1068, which contains an active galactic nucleus (AGN). The effect of diffuse dust heating by AGN (beyond the torus) was so far only investigated for quasars. This additional dust heating source further contaminates the broadband fluxes on which classic galaxy modelling tools rely to derive physical properties. We aim to fit a realistic model to the observations of NGC 1068 and quantify the contribution of the several dust heating sources. Our model is able to reproduce the global spectral energy distribution of the galaxy. It matches the resolved optical and infrared images fairly well, but deviates in the UV and the submm. We find a strong wavelength dependency of AGN contamination to the broadband fluxes. It peaks in the MIR, drops in the FIR, but rises again at submm wavelengths. We quantify the contribution of the dust heating sources in each 3D dust cell and find a median value of 83% for the star formation component. The AGN contribution is measurable at the percentage level in the disc, but quickly increases in the inner few 100 pc, peaking above 90%. This is the first time the phenomenon of an AGN heating the diffuse dust beyond its torus is quantified in a nearby star-forming galaxy. NGC 1068 only contains a weak AGN, meaning this effect can be stronger in galaxies with a more luminous AGN. This could significantly impact the derived star formation rates and ISM masses for such systems.

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)

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.

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.

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.

This paper contains the results from the first 30 months of the NELIOTA project for impacts of NEOs/meteoroids on the lunar surface. Our analysis on the statistics concerning the efficiency of the campaign and the parameters of the projectiles and the impacts is presented. The parameters of the lunar impact flashes based on simultaneous observations in two wavelength bands are used to estimate the distributions of the masses, sizes and frequency of the impactors. These statistics can be used both in space engineering and science. The photometric fluxes of the flashes are measured using aperture photometry and their apparent magnitudes are calculated using standard stars. Assuming that the flashes follow a black body law of irradiation, the temperatures can be derived analytically, while the parameters of the projectiles are estimated using fair assumptions on their velocity and luminous efficiency of the impacts. 79 lunar impact flashes have been observed with the 1.2 m Kryoneri telescope in Greece. The masses of the meteoroids range between 0.7 g and 8 kg and their respective sizes between 1-20 cm depending on their assumed density, impact velocity, and luminous efficiency. We find a strong correlation between the observed magnitudes of the flashes and the masses of the meteoroids. Moreover, an empirical relation between the emitted energies of each band has been derived allowing the estimation of the physical parameters of the meteoroids that produce low energy impact flashes. The NELIOTA project has so far the highest detection rate and the faintest limiting magnitude for lunar impacts compared to other ongoing programs. Based on the impact frequency distribution on Moon, we estimate that sporadic meteoroids with typical masses less than 100 g and sizes less than 5 cm enter the mesosphere of the Earth with a rate ~108 meteoroids/hr and also impact Moon with a rate of ~8 meteoroids/hr.

Observations of evolution in the dust-to-metal ratio allow us to constrain the dominant dust processing mechanisms. In this work, we present a study of the dust-to-metal and dust-to-gas ratios in a sub-sample of ~500 DustPedia galaxies. Using literature and MUSE emission line fluxes, we derived gas-phase metallicities (oxygen abundances) for over 10000 individual regions and determine characteristic metallicities for each galaxy. We study how the relative dust, gas, and metal contents of galaxies evolve by using metallicity and gas fraction as proxies for evolutionary state. The global oxygen abundance and nitrogen-to-oxygen ratio are found to increase monotonically as galaxies evolve. Additionally, unevolved galaxies (gas fraction > 60%, metallicity 12 + log(O/H) < 8.2) have dust-to-metal ratios that are about a factor of 2.1 lower (a factor of six lower for galaxies with gas fraction > 80%) than the typical dust-to-metal ratio (Md/MZ ~ 0.214) for more evolved sources. However, for high gas fractions, the scatter is larger due to larger observational uncertainties as well as a potential dependence of the dust grain growth timescale and supernova dust yield on local conditions and star formation histories. We find chemical evolution models with a strong contribution from dust grain growth describe these observations reasonably well. The dust-to-metal ratio is also found to be lower for low stellar masses and high specific star formation rates (with the exception of some sources undergoing a starburst). Finally, the metallicity gradient correlates weakly with the HI-to-stellar mass ratio, the effective radius and the dust-to-stellar mass ratio, but not with stellar mass.

The dust mass absorption coefficient, $\kappa_{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 $\kappa_{d}$ in galaxies. We determine $\kappa_{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 $\kappa_{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 $\kappa_{d}$. We find values of $\kappa_{d}$ at 500 $\mu$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 $\kappa_{d}$ shows a distinct inverse correlation with the local density of the interstellar medium. This inverse correlation is the opposite of what is predicted by standard dust models. However, we find this relationship to be robust against a large range of changes to our method - only the adoption of unphysical or highly unusual assumptions would be able to suppress it.

Aims: We compare the far-infrared to sub-millimetre dust emission properties measured in high Galactic latitude cirrus with those determined in a sample of 204 late-type DustPedia galaxies. The aim is to verify if it is appropriate to use Milky Way dust properties to derive dust masses in external galaxies. Methods: We used Herschel observations and atomic and molecular gas masses to estimate the disc-averaged dust emissivity at 250 micrometres, and from this, the absorption cross section per H atom and per dust mass. The emissivity requires one assumption, which is the CO-to-H_2 conversion factor, and the dust temperature is additionally required for the absorption cross section per H atom; yet another constraint on the dust-to-hydrogen ratio D/H, depending on metallicity, is required for the absorption cross section dust mass. Results: We find epsilon(250) = 0.82 +/- 0.07 MJy sr^-1 (1E20 H cm^-2)^-1 for galaxies with 4 < F(250)/F(500) < 5. This depends only weakly on the adopted CO-to-H_2 conversion factor. The value is almost the same as that for the Milky Way at the same colour ratio. Instead, for F(250)/F(500) > 6, epsilon(250) is lower than predicted by its dependence on the heating conditions. The reduction suggests a variation in dust emission properties for spirals of earlier type, higher metallicity, and with a higher fraction of molecular gas. When the standard emission properties of Galactic cirrus are used for these galaxies, their dust masses might be underestimated by up to a factor of two. Values for the absorption cross sections at the Milky Way metallicity are also close to those of the cirrus. Mild trends of the absorption cross sections with metallicity are found, although the results depend on the assumptions made.

We use a sub-set of the DustPedia galaxy sample (461 galaxies) to investigate the effect the environment has had on galaxies. We consider Virgo cluster and field samples and also assign a density contrast parameter to each galaxy, as defined by the local density of SDSS galaxies. We consider their chemical evolution (using M_{Dust}/M_{Baryon} and M_{Gas}/M_{Baryon}), their specific star formation rate (SFR/M_{Stars}), star formation efficiency (SFR/M_{Gas}), stars-to-dust mass ratio (M_{Stars}/M_{Dust}), gas-to-dust mass ratio (M_{Gas}/M_{Dust}) and the relationship between star formation rate per unit mass of dust and dust temperature (SFR/M_{Dust} and T_{Dust}). Late type galaxies (later than Sc) in all of the environments can be modelled using simple closed box chemical evolution and a simple star formation history (SFR(t) \propto t\exp{-t/\tau}). For earlier type galaxies the physical mechanisms that give rise to their properties are clearly much more varied and require a more complicated model (mergers, gas in or outflow). However, we find little or no difference in the properties of galaxies of the same morphological type within the cluster, field or with different density contrasts. It appears that it is morphology, how and whenever this is laid down, and consistent internal physical processes that primarily determine the derived properties of galaxies in the DustPedia sample and not processes related to differences in the local environment.

Within the framework of the DustPedia project we investigate the properties of cosmic dust and its interaction with the stellar radiation (originating from different stellar populations) for 814 galaxies in the nearby Universe, all observed by the Herschel Space Observatory. We take advantage of the widely used galaxy SED fitting code CIGALE, properly adapted to include the state-of-the-art dust model THEMIS. Using the DustPedia photometry we determine the physical properties of the galaxies, such as, the dust and stellar mass, the star-formation rate, the bolometric luminosity as well as the unattenuated and the absorbed by dust stellar light, for both the old (> 200 Myr) and young (<= 200 Myr) stellar populations. We show how the mass of stars, dust, and atomic gas, as well as the star-formation rate and the dust temperature vary between galaxies of different morphologies and provide recipes to estimate these parameters given their Hubble stage (T). We find a mild correlation between the mass fraction of the small a-C(:H) grains with the specific star-formation rate. On average, young stars are very efficient in heating the dust, with absorption fractions reaching as high as ~77% of the total, unattenuated luminosity of this population. On the other hand, the maximum absorption fraction of old stars is ~24%. Dust heating in early-type galaxies is mainly due to old stars, up to a level of ~90%. Young stars progressively contribute more for `typical' spiral galaxies and they become the dominant source of dust heating for Sm type and irregular galaxies, donating up to ~60% of their luminosity to this purpose. Finally, we find a strong correlation of the dust heating fraction by young stars with morphology and the specific star-formation rate.

NELIOTA (NEO Lunar Impacts and Optical TrAnsients) is an ESA-funded lunar monitoring project, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. A prime focus, high-speed, twin-camera system providing simultaneous observations in two photometric bands at a rate of 30 frames-per-second on the 1.2 m Kryoneri telescope of the National Observatory of Athens was commissioned for this purpose. A dedicated software processes the images and automatically detects candidate lunar impact flashes, which are then validated by an expert user. The four year observing campaign began in February 2017 and has so far detected more than 40 lunar impact events. The software routinely detects satellites, which typically appear as streaks or dots crossing the lunar disk. To avoid confusing these events with real flashes, we check different available catalogs with spacecraft orbital information and exclude spacecraft identifications.

Most radiative transfer models assume that dust in spiral galaxies is distributed exponentially. In this paper our goal is to verify this assumption by analysing the two-dimensional large-scale distribution of dust in galaxies from the DustPedia sample. For this purpose, we make use of Herschel imaging in five bands, from 100 to 500{\mu}m, in which the cold dust constituent is primarily traced and makes up the bulk of the dust mass in spiral galaxies. For a subsample of 320 disc galaxies, we successfully perform a simultaneous fitting with a single S\'ersic model of the Herschel images in all five bands using the multiband modelling code GALFITM. We report that the S\'ersic index $n$, which characterises the shape of the S\'ersic profile, lies systematically below 1 in all Herschel bands and is almost constant with wavelength. The average value at 250{\mu}m is $0.67\pm0.37$ (187 galaxies are fitted with $n_{250}\leq0.75$, 87 galaxies have $0.75<n_{250}\leq1.25$, and 46 - with $n_{250}>1.25$). Most observed profiles exhibit a depletion in the inner region (at $r<0.3-0.4$ of the optical radius $r_{25}$ ) and are more or less exponential in the outer part. We also find breaks in the dust emission profiles at longer distances $(0.5-0.6)r_{25}$ which are associated with the breaks in the optical and near-infrared. We assume that the observed deficit of dust emission in the inner galaxy region is related to the depression in the radial profile of the HI surface density in the same region because the atomic gas reaches high enough surface densities there to be transformed into molecular gas. If a galaxy has a triggered star formation in the inner region (for example, because of a strong bar instability, which transfers the gas inwards to the centre, or a pseudobulge formation), no depletion or even an excess of dust emission in the centre is observed.

We study the fraction of stellar radiation absorbed by dust, f_abs, in 814 galaxies of different morphological types. The targets constitute the vast majority (93%) of the DustPedia sample, including almost all large (optical diameter larger than 1'), nearby (v <= 3000 km/s) galaxies observed with the Herschel Space Observatory. For each object, we model the spectral energy distribution from the ultraviolet to the sub-millimetre using the dedicated, aperture-matched DustPedia photometry and the fitting code CIGALE. The value of f_abs is obtained from the total luminosity emitted by dust and from the bolometric luminosity, which are estimated by the fit. On average, 19% of the stellar radiation is absorbed by dust in DustPedia galaxies. The fraction rises to 25% if only late-type galaxies are considered. The dependence of f_abs on morphology, showing a peak for Sb-Sc galaxies, is weak; it reflects a stronger, yet broad, positive correlation with the bolometric luminosity, which is identified for late-type, disk-dominated, high-specific-star-formation rate, gas-rich objects. We find no variation of f_abs with inclination, at odds with radiative transfer models of edge-on galaxies. These results call for a self-consistent modelling of the evolution of the dust mass and geometry along the build-up of the stellar content. We also provide template spectral energy distributions in bins of morphology and luminosity and study the variation of f_abs with stellar mass and specific star formation rate. We confirm that the local Universe is missing the high f_abs}, luminous and actively star-forming objects necessary to explain the energy budget in observations of the extragalactic background light.

We present the technical specifications and first results of the ESA-funded, lunar monitoring project "NELIOTA" (NEO Lunar Impacts and Optical TrAnsients) at the National Observatory of Athens, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. For the purposes of this project a twin camera instrument was specially designed and installed at the 1.2 m Kryoneri telescope utilizing the fast-frame capabilities of scientific Complementary Metal-Oxide Semiconductor detectors (sCMOS). The system provides a wide field-of-view (17.0' $\times$ 14.4') and simultaneous observations in two photometric bands (R and I), reaching limiting magnitudes of 18.7 mag in 10 sec in both bands at a 2.5 signal-to-noise level. This makes it a unique instrument that can be used for the detection of NEO impacts on the Moon, as well as for any astronomy projects that demand high-cadence multicolor observations. The wide field-of-view ensures that a large portion of the Moon is observed, while the simultaneous, high-cadence, monitoring in two photometric bands makes possible, for the first time, the determination of the temperatures of the impacts on the Moon's surface and the validation of the impact flashes from a single site. Considering the varying background level on the Moon's surface we demonstrate that the NELIOTA system can detect NEO impact flashes at a 2.5 signal-to-noise level of ~12.4 mag in the I-band and R-band for observations made at low lunar phases ~0.1. We report 31 NEO impact flashes detected during the first year of the NELIOTA campaign. The faintest flash was at 11.24 mag in the R-band (about two magnitudes fainter than ever observed before) at lunar phase 0.32. Our observations suggest a detection rate of $1.96 \times 10^{-7}$ events $km^{-2} h^{-1}$.

We report the first scientific results from the NELIOTA (NEO Lunar Impacts and Optical TrAnsients) project, which has recently begun lunar monitoring observations with the 1.2-m Kryoneri telescope. NELIOTA aims to detect faint impact flashes produced by near-Earth meteoroids and asteroids and thereby help constrain the size-frequency distribution of near-Earth objects in the decimeter to meter range. The NELIOTA setup, consisting of two fast-frame cameras observing simultaneously in the $R$ and $I-$bands, enables - for the first time - direct analytical calculation of the flash temperatures. We present the first 10 flashes detected, for which we find temperatures in the range ~1,600-3,100 K, in agreement with theoretical values. Two of these flashes were detected on multiple frames in both filters and therefore yield the first measurements of the temperature drop for lunar flashes. In addition, we compute the impactor masses, which range between ~100 g and ~50 kg.

The purpose of this work is the characterization of the radial distribution of dust, stars, gas, and star-formation rate (SFR) in a sub-sample of 18 face-on spiral galaxies extracted from the DustPedia sample. This study is performed by exploiting the multi-wavelength, from UV to sub-mm bands, DustPedia database, in addition to molecular (12CO) and atomic (HI) gas maps and metallicity abundance information available in the literature. We fitted the surface brightness profiles of the tracers of dust and stars, the mass surface density profiles of dust, stars, molecular gas, and total gas, and the SFR surface density profiles with an exponential curve and derived their scale-lengths. We also developed a method to solve for the CO-to-H2 conversion factor (alpha_CO) per galaxy by using dust and gas mass profiles. Although each galaxy has its own peculiar behaviour, we identified a common trend of the exponential scale-lengths vs. wavelength. On average, the scale-lengths normalized to the B-band 25 mag/arcsec^2 radius decrease from UV to 70 micron, from 0.4 to 0.2, and then increase back up to 0.3 at 500 microns. The main result is that, on average, the dust mass surface density scale-length is about 1.8 times the stellar one derived from IRAC data and the 3.6 micron surface brightness, and close to that in the UV. We found a mild dependence of the scale-lengths on the Hubble stage T: the scale-lengths of the Herschel bands and the 3.6 micron scale-length tend to increase from earlier to later types, the scale-length at 70 micron tends to be smaller than that at longer sub-mm wavelength with ratios between longer sub-mm wavelengths and 70 micron that decrease with increasing T. The scale-length ratio of SFR and stars shows a weak increasing trend towards later types.

We provide the basic integrated physical properties of all the galaxies contained in the full CASSIS with available broad-band photometry from UV to 22 $\mu$m. We have collected photometric measurements in 14 wavelengths from available public surveys in order to study the SED of each galaxy in CASSIS, thus constructing a final sample of 1,146 galaxies in the redshift range 0<z<2.5. The SEDs are modelled with the CIGALE code which relies on the energy balance between the absorbed stellar and the dust emission while taking into account the possible contribution due to the presence of an AGN. We split the galaxies in three groups, a low-redshift, a mid-redshift and a high-redshift and find that the vast majority of the Spitzer/IRS galaxies are star-forming and lie on or above the star-forming MS of the corresponding redshift. Moreover, the emission of Spitzer/IRS galaxies with z<0.1 is mostly dominated by star-formation, galaxies in the mid-redshift bin are a mixture of star forming and AGN galaxies, while half of the galaxies with z<0.5 show moderate or high AGN activity. Additionally, using rest-frame NUV-r colour, S\'ersic indices, optical OIII and NII emission lines we explore the nature of these galaxies by investigating their structure as well as their star-formation and AGN activity. Using a colour magnitude diagram we confirm that 97% of the galaxies with redshift smaller than 0.5 have experienced a recent star-formation episode. For a sub-sample of galaxies with available structural information and redshift smaller than 0.3 we find that early-type galaxies are placed below the MS, while late-type galaxies are found on the MS as expected. Finally, for all the galaxies with z<0.5 and available optical spectral line measurements we find that galaxies with high AGN luminosity, as calculated by CIGALE, are most likely to be classified as composite or AGNs by optical spectral lines.

Several dedicated surveys focusing on early-type galaxies (ETGs) reveal that significant fractions of them are detectable in all interstellar medium phases studied to date. We select ETGs from the Herschel Reference Survey that have both far-infrared Herschel and either HI or CO detection (or both). We derive their star formation rates (SFR), stellar masses and dust masses via modelling their spectral energy distributions. We combine these with literature information on their atomic and molecular gas properties, in order to relate their star formation, total gas mass and dust mass on global scales. The ETGs deviate from the dust mass-SFR relation and the Schmidt-Kennicutt relation that SDSS star forming galaxies define: compared to SDSS galaxies, ETGs have more dust at the same SFR, or less SFR at the same dust mass. When placing them in the M*-SFR plane, ETGs show a much lower specific SFR as compared to normal star-forming galaxies. ETGs show a large scatter compared to the Schmidt-Kennicutt relation found locally within our Galaxy, extending to lower SFRs and gas mass surface densities. Using an ETG's SFR and the Schmidt-Kennicutt law to predict its gas mass leads to an underestimate. ETGs have similar observed-gas-to-modelled-dust mass ratios to star forming-galaxies of the same stellar mass, as well as they exhibit a similar scatter.

Images of dust continuum and carbon monoxide (CO) line emission are powerful tools for deducing structural characteristics of galaxies, such as disk sizes, H$_2$ gas velocity fields and enclosed H$_2$ and dynamical masses. We report on a fundamental constraint set by the cosmic microwave background (CMB) on the observed structural and dynamical characteristics of galaxies, as deduced from dust continuum and CO-line imaging at high redshifts. As the CMB temperature rises in the distant Universe, the ensuing thermal equilibrium between the CMB and the cold dust and H$_2$ gas progressively erases all spatial and spectral contrasts between their brightness distributions and the CMB. For high-redshift galaxies, this strongly biases the recoverable H$_2$ gas and dust mass distributions, scale lengths, gas velocity fields and dynamical mass estimates. This limitation is unique to mm/submm wavelengths and unlike its known effect on the global dust continuum and molecular line emission of galaxies, it cannot be addressed simply. We nevertheless identify a unique signature of CMB-affected continuum brightness distributions, namely an increasing rather than diminishing contrast between such brightness distributions and the CMB when the cold dust in distant galaxies is imaged at frequencies beyond the Raleigh-Jeans limit. For the molecular gas tracers, the same effect makes the atomic carbon (CI) lines maintain a larger contrast than the CO lines against the CMB.

Previous studies have shown the existence of an excess of emission at submillimeter (submm) and millimeter (mm) wavelengths in the spectral energy distribution (SED) of many low-metallicity galaxies. The goal of the present study is to model separately the emission from the star forming (SF) component and the emission from the diffuse interstellar medium (ISM) in the nearby spiral galaxy M33. We decomposed the observed SED of M33 into its SF and diffuse components. Mid-infrared (MIR) and far-infrared (FIR) fluxes were extracted from Spitzer and Herschel data. At submm and mm wavelengths, we used ground-based observations from APEX to measure the emission from the SF component and data from the Planck space telescope to estimate the diffuse emission. Both components were separately fitted using radiation transfer models based on standard dust properties and a realistic geometry. Both modeled SEDs were combined to build the global SED of M33. In addition, the radiation field necessary to power the dust emission in our modeling was compared with observations from GALEX, Sloan, and Spitzer. Our modeling is able to reproduce the observations at MIR and FIR wavelengths, but we found a strong excess of emission at submm and mm wavelengths, where the model expectations severely underestimate the LABOCA and Planck fluxes. We also found that the ultraviolet (UV) radiation escaping the galaxy is 70% higher than the model predictions. We determined a gas-to-dust mass ratio Gdust~100, significantly lower than the value expected from the sub-solar metallicity of M33. We discussed different hypotheses to explain the discrepancies found in our study (i.e., excess of emission at submm and mm wavelengths, deficit of UV attenuation, and abnormally low value of Gdust), concluding that different dust properties in M33 is the most plausible explanation.

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $\mu$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $\mu$m, [OI] 63 $\mu$m, and [OI] 145 $\mu$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $\mu$m, and [CI] 609 $\mu$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $\mu$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $\mu$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = \alpha \log L'_{\rm CO} + \beta$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($\alpha \sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($\beta \sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

Understanding the heating and cooling mechanisms in nearby (Ultra) luminous infrared galaxies can give us insight into the driving mechanisms in their more distant counterparts. Molecular emission lines play a crucial role in cooling excited gas, and recently, with Herschel Space Observatory we have been able to observe the rich molecular spectrum. CO is the most abundant and one of the brightest molecules in the Herschel wavelength range. CO transitions are observed with Herschel, and together, these lines trace the excitation of CO. We study Arp 299, a colliding galaxy group, with one component harboring an AGN and two more undergoing intense star formation. For Arp 299 A, we present PACS spectrometer observations of high-J CO lines up to J=20-19 and JCMT observations of $^{13}$CO and HCN to discern between UV heating and alternative heating mechanisms. There is an immediately noticeable difference in the spectra of Arp 299 A and Arp 299 B+C, with source A having brighter high-J CO transitions. This is reflected in their respective spectral energy line distributions. We find that photon-dominated regions (PDRs) are unlikely to heat all the gas since a very extreme PDR is necessary to fit the high-J CO lines. In addition, this extreme PDR does not fit the HCN observations, and the dust spectral energy distribution shows that there is not enough hot dust to match the amount expected from such an extreme PDR. Therefore, we determine that the high-J CO and HCN transitions are heated by an additional mechanism, namely cosmic ray heating, mechanical heating, or X-ray heating. We find that mechanical heating, in combination with UV heating, is the only mechanism that fits all molecular transitions. We also constrain the molecular gas mass of Arp 299 A to 3e9 Msun and find that we need 4% of the total heating to be mechanical heating, with the rest UV heating.

We used the SPIRE/FTS instrument aboard the Herschel Space Observatory (HSO) to obtain the Spectral Line Energy Distributions (SLEDs) of CO from J=4-3 to J=13-12 of Arp 193 and NGC 6240, two classical merger/starbursts selected from our molecular line survey of local Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11} L_{sol}). The high-J CO SLEDs are then combined with ground-based low-J CO, {13}CO, HCN, HCO+, CS line data and used to probe the thermal and dynamical states of their large molecular gas reservoirs. We find the two CO SLEDs strongly diverging from J=4-3 onwards, with NGC6240 having a much higher CO line excitation than Arp193, despite their similar low-J CO SLEDs and L_{FIR}/L_{CO,1-0}, L_{HCN}/L_{CO} (J=1-0) ratios (proxies of star formation efficiency and dense gas mass fraction). In Arp193, one of the three most extreme starbursts in the local Universe, the molecular SLEDs indicate a small amount ~(5-15)% of dense gas (n>=10^{4}cm^{-3}) unlike NGC6240 where most of the molecular gas (~(60-70)%) is dense n~(10^4-10^5)cm^{-3}. Strong star-formation feedback can drive this disparity in their dense gas mass fractions, and also induce extreme thermal and dynamical states for the molecular gas.In NGC6240, and to a lesser degree in Arp193, we find large molecular gas masses whose thermal states cannot be maintained by FUV photons from Photon Dominated Regions (PDRs). We argue that this may happen often in metal-rich merger/starbursts, strongly altering the initial conditions of star formation. ALMA can now directly probe these conditions across cosmic epoch, and even probe their deeply dust-enshrouded outcome, the stellar IMF averaged over galactic evolution.

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

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 aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%.

Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher.

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 study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).

We report results from a large molecular line survey of Luminous Infrared Galaxies (L_{IR} >= 10^{11} L_sol) in the local Universe (z<=0.1), conducted during the last decade with the James Clerk Maxwell Telescope (JCMT) and the IRAM 30-m telescope. This work presents the CO and {13}CO line data for 36 galaxies, further augmented by multi-J total CO luminosities available for other IR-bright galaxies from the literature. This yields a sample of N=70 galaxies with the star-formation (SF) powered fraction of their IR luminosities spanning L_{IR} (10^{10}-2x10^{12}) L_sol and a wide range of morphologies. Simple comparisons of their available CO Spectral Line Energy Distributions (SLEDs) with local ones, as well as radiative transfer models discern a surprisingly wide range of average ISM conditions, with most of the surprises found in the high-excitation regime. These take the form of global CO SLEDs dominated by a very warm (T_{kin}>=100 K) and dense (n>=10^4 cm^{-3}) gas phase, involving galaxy-sized (~(few)x10^9 M_sol) gas mass reservoirs under conditions that would otherwise amount only ~1% of mass per typical SF molecular cloud in the Galaxy. Some of the highest excitation CO SLEDs are found in the so-called Ultra Luminous Infrared Galaxies and seem irreducible to ensembles of ordinary SF-powered regions. Highly supersonic turbulence and high cosmic ray (CR) energy densities rather than far-UV/optical photons or SNR-induced shocks from individual SF sites can globally warm the large amounts of dense gas found in these merger-driven starbursts and easily power their extraordinary CO line excitation.....

In this work we conclude the analysis of our CO line survey of Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11}L_{sol}) in the local Universe (Paper\,I), by focusing on the influence of their average ISM properties on the total molecular gas mass estimates via the so-called X_{co}=M(H_2)/L_{co,1-0} factor. One-phase radiative transfer models of the global CO Spectral Line Energy Distributions (SLEDs) yield an X_{co} distribution with: <X_{co}>\sim(0.6+/-0.2) M_{sol}(K km s^{-1} pc^2)^{-1} over a significant range of average gas densities, temperatures and dynamical states. The latter emerges as the most important parameter in determining X_{co}, with unbound states yielding low values and self-gravitating states the highest ones. Nevertheless in many (U)LIRGs where available higher-J CO lines (J=3--2, 4--3, and/or J=6--5) or HCN line data from the literature allow a separate assessment of the gas mass at high densities (>=10^{4} cm^{-3}) rather than a simple one-phase analysis we find that {\it near-Galactic X_{co} (3-6)\, M_sol\,(K\,km^{-1}\,pc^2)^{-1} values become possible.} We further show that in the highly turbulent molecular gas in ULIRGs a high-density component will be common and can be massive enough for its high X_{co} to dominate the average value for the entire galaxy. ......... ...this may have thus resulted to systematic underestimates of molecular gas mass in ULIRGs.

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.

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

Within the framework of the HerM33es Key Project for Herschel and in combination with multi-wavelength data, we study the Spectral Energy Distribution (SED) of a set of HII regions in the Local Group Galaxy M33. Using the Halpha emission, we perform a classification of a selected HII region sample in terms of morphology, separating the objects in filled, mixed, shell and clear shell objects. We obtain the SED for each HII region as well as a representative SED for each class of objects. We also study the emission distribution of each band within the regions. We find different trends in the SEDs for each morphological type that are related to properties of the dust and their associated stellar cluster. The emission distribution of each band within the region is different for each morphological type of object.

Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 \mu m. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 \mu m to 500 \mu m, along with HI, H\alpha\ maps, and GALEX ultraviolet data we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations.

We analyse Herschel/SPIRE images of the edge-on spiral galaxy NGC 891 at 250, 350 and 500 micron. Using a 3D radiative transfer model we confirm that the dust has a radial fall-off similar to the stellar disk. The dust disk shows a break at about 12 kpc from the center, where the profile becomes steeper. Beyond this break, emission can be traced up to 90% of the optical disk in the NE side. On the SW, we confirm dust emission associated with the extended, asymmetric HI disk, previously detected by the Infrared Space Observatory (ISO). This emission is marginally consistent with the large diffuse dust disk inferred from radiative transfer fits to optical images. No excess emission is found above the plane beyond that of the thin, unresolved, disk.

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

We report the results of a pilot study with the EVLA of 12CO J=1-0 emission from four SMGs at z=2.2-2.5, each with an existing detection of CO J=3-2. Using the EVLA's most compact configuration we detect strong, broad J=1-0 line emission from all of our targets. The median line width ratio, sigma(1-0)/sigma(3-2) = 1.15 +/- 0.06, suggests that the J=1-0 is more spatially extended than the J=3-2 emission, a situation confirmed by our maps which reveal velocity structure in several cases and typical sizes of ~16 kpc FWHM. The median Tb ratio is r(3-2/1-0) = 0.55 +/- 0.05, noting that our value may be biased high because of the J=3-2-based sample selection. Naively, this suggests gas masses ~2x higher than estimates made using higher-J transitions of CO, with the discrepency due to the difference in assumed Tb ratio. We also estimate masses using the 12CO J=1-0 line and the observed global Tb ratios, assuming standard underlying Tb ratios as well as a limiting SFE, i.e. without calling upon X(CO). Using this new method, we find a median molecular gas mass of (2.5 +/- 0.8) x 10^10 Msun, with a plausible range stretching 3x higher. Even larger masses cannot be ruled out, but are not favoured by dynamical constraints: the median dynamical mass for our sample is (2.3 +/- 1.4) x 10^11 Msun. We examine the Schmidt-Kennicutt relation for all the distant galaxy populations for which CO J=1-0 or J=2-1 data are available, finding small systematic differences. These have previously been interpreted as evidence for different modes of star formation, but we argue that these differences are to be expected, given the still considerable uncertainties. Finally, we discuss the morass of degeneracies surrounding molecular gas mass estimates, the possibilities for breaking them, and the future prospects for imaging and studying cold, quiescent molecular gas at high redshifts [abridged].

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.

Recently, Papadopoulos et al., 2010 using sub-mm CO molecular line observations of nearby ultra-luminous IRAS galaxies, (U)LIRGs, have found that exceptionally large gas column densities (N_H > 10^25 cm-2) can be present across some of the very dense gaseous disks that are typically found in these objects. They also proposed a diagnostic for finding such sources using CO and HCN molecular lines. Given that such high column densities are expected to absorb any X-ray luminous AGN, yielding Compton-thick sources, we set out toexplore whether this can be discerned using X-ray observations. More specifically we examine X-ray spectral observations of 14 sources in their sample, using public Chandra observations (0.5-10 keV) for eleven sources as well as BeppoSAX results (2-100 keV) from the literature for another three sources. Our goal is to find candidate Compton-thick AGN and to check whether the molecular line selection criterion is successful in selecting such systems. X-ray spectroscopy reveals four candidate Compton-thick AGN of which half fall within the high obscuration region in the molecular line ratio diagnostics. Of the remaining five sources falling into the `high dust obscuration' box, one (Mrk273) is highly obscured (N_H ~4x10^23 cm-2) while in the other four the X-ray emission is most probably associated with star-forming processes rather than an AGN on the basis of their X-ray and mid-infrared properties. Overall, we argue that although this method as expected cannot recover all Compton-thick AGN, there are no examples of X-ray luminous AGN inside that region that have low obscuration, suggesting that this method is efficient in finding heavily obscured AGN in dust-enshrouded star-forming galaxies. The above results bear important implications for future joint ALMA and X-ray observations for the detection of Compton-thick AGN.

Hickson Compact Groups (HCGs) are among the densest galaxy environments of the local universe. To examine the effects of the environment on the infrared properties of these systems, we present an analysis of Spitzer and ISO mid-infrared imaging as well as deep ground based near-infrared imaging of 14 HCGs containing a total of 69 galaxies. Based on mid-infrared color diagnostics we identify the galaxies which appear to host an active nucleus, while using a suite of templates, we fit the complete infrared spectral energy distribution for each group member. We compare our estimates of galaxy mass, star formation rate, total infrared luminosities, and specific star formation rates (sSFR) for our HCG sample, to samples of isolated galaxies and interacting pairs and find that overall there is no discernible difference among them. However, HCGs which can be considered as dynamically "old", host late-type galaxies with a slightly lower sSFR than the one found in dynamically "young" groups. This could be attributed to multiple past interactions among the galaxies in old groups, that have led to the build up of their stellar mass. It is also consistent with our prediction for the presence of diffuse cold dust in the intergalactic medium of 9 of the dynamically "old" groups.

We examine the properties of the X-ray detected, Infrared Excess AGN or Dust Obscured Galaxies (DOGs) in the Chandra Deep Fields (CDF). We find 26 X-ray selected sources which obey the 24 micron to R-band flux ratio criterion f_24/f_R>1000. These are at a median redshift of 2.3 while their IR luminosities are above 10^12 solar. Their X-ray luminosities are all above a few times 10^42 erg s-1 in the 2-10 keV band unambiguously arguing that these host AGN. Nevertheless, their IR Spectral Energy Distributions are split between AGN (Mrk231) and star-forming templates (Arp220). Our primary goal is to examine their individual X-ray spectra in order to assess whether this X-ray detected DOG population contains heavily obscured or even Compton-thick sources. The X-ray spectroscopy reveals a mixed bag of objects. We find that four out of the 12 sources with adequate photon statistics and hence reliable X-ray spectra, show evidence for a hard X-ray spectral index (~1) or harder,consistent with a Compton-thick spectrum. In total 12 out of the 26 DOGs show evidence for flat spectral indices. However, owing to the limited photon statistics we cannot differentiate whether these are flat because they are reflection-dominated or because they show moderate amounts of absorption. Seven DOGs show relatively steep spectra (>1.4) indicative of small column densities. All the above suggest a fraction of Compton-thick sources that does not exceed 5%. The average X-ray spectrum of all 26 DOGs is hard (~1.1) or even harder (~0.6) when we exclude the brightest sources. These spectral indices are well in agreement with the stacked spectrum of X-ray undetected sources (~0.8 in the CDFN). This could suggest (but not necessarily prove) that X-ray undetected DOGs, in a similar fashion to the X-ray detected ones presented here, are hosting a moderate fraction of Compton-thick sources.

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

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

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust.

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

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.

We report on new sensitive CO J=6-5 line observations of several luminous infrared Galaxies (LIRGs: L$_{\rm IR}$(8-1000$\mu $m$)\ga 10^{11}$ L$_{\odot}$), 36% (8/22) of them ULIRGs (L$_{\rm IR}$$>10^{12}$ L$_{\odot}$), and two powerful local AGN: the optically luminous QSO PG 1119+120, and the powerful radio galaxy 3C 293 using the James Clerk Maxwell Telescope (JCMT) on Mauna Kea in Hawaii. We combine these observations with existing low-J CO data and dust emission Spectral Energy Distributions (SEDs) in the far-infrared - submillimetre from the literature to constrain the properties of the star-forming ISM in these systems. We then build the first {\it local} CO Spectral Line Energy Distributions (SLEDs) for the {\it global} molecular gas reservoirs that reach up to high J-levels. These CO SLEDs are neither biased by strong lensing (which affects many of those constructed for high-redshift galaxies), nor suffer from undersampling of CO-bright regions (as most current high-J CO observations of nearby extended systems do). We find: ...

We present a new observational study of the CO(1-0) line emission as an H2 gas mass tracer under extreme conditions in extragalactic environments. Our approach is to study the full neutral interstellar medium (H2, HI and dust) of two galaxies whose bulk interstellar medium (ISM) resides in environments that mark (and bracket) the excitation extremes of the ISM conditions found in infrared luminous galaxies, the starburst NGC3310 and the quiescent spiral NGC157. Our study maintains a robust statistical notion of the so-called X factor (i.e. a large ensemble of clouds is involved) while exploring its dependency on the very different average ISM conditions prevailing within these two systems. These are constrained by fully-sampled CO(3-2) and CO(1-0) observations, at a matched beam resolution of Half Power Beam Width 15", obtained with the JCMT the Nobeyama 45-m telescope, combined with sensitive 850 and 450 micron dust emission and HI interferometric images which allow a complete view of all the neutral ISM components. We found an X factor varying by a factor of 5 within the spiral galaxy NGC157 and about 2 times lower than the Galactic value in NGC3310. In addition, the dust emission spectrum in NGC3310 shows a pronounced submm "excess". We explain and fit this by an enhanced abundance of Very Small Grains (VSGs). The alternative approach to interpret this excess as very cold dust would yield a cold dust temperature of 5-11 K with a correspondingly very low gas-to-dust mass ratio of 5-43, very unlikely for this starburst galaxy with low metallicity.

Deep optical CCD images of the supernova remnant G 32.8-0.1 were obtained where filamentary and diffuse emission was discovered. The images were acquired in the emission lines of Halpha+[N II] and [S II]. Filamentary and diffuse structures are detected in most areas of the remnant, while no significant [O III] emission is present. The flux-calibrated images suggest that the optical emission originates from shock-heated gas since the [S II]/Halpha ratio is greater than 1.2. The Spitzer images at 8 micron and 24 micron show a few filamentary structures to be correlated with the optical filaments, while the radio emission at 1.4 GHz in the same area is found to be very well correlated with the brightest optical filaments. Furthermore, the results from deep long-slit spectra also support the origin of the emission to be from shock-heated gas ([S II]/Halpha > 1.5). The absence of [O III] emission indicates slow shocks velocities into the interstellar "clouds" (< 100 km/s), while the [S II] 6716/6731 ratio indicates electron densities up to ~200 cm^{-3}. Finally, the Halpha emission is measured to lie between 1.8 to 4.6 x 10^{-17} erg/s/cm^2/arcsec^2, while from VGPS HI images a distance to the SNR is estimated to be between 6 to 8.5 kpc.

We present the results from a study of the X-ray variability and the near-IR to X-ray spectral energy distribution of four low-luminosity, Seyfert 1 galaxies. We compared their variability amplitude and broad band spectrum with those of more luminous AGN in order to investigate whether accretion in low-luminosity AGN operates as in their luminous counterparts. We used archival XMM-Newton and, in two cases, ASCA data to estimate their X-ray variability amplitude and determine their X-ray spectral shape and luminosity. We also used archival HST data to measure their optical nuclear luminosity, and near-IR measurements from the literature, in order to construct their near-IR to X-ray spectra. The X-ray variability amplitude of the four Seyferts is what one would expect, given their black hole masses. Their near-IR to X-ray spectrum has the same shape as the spectrum of quasars which are 10^2-10^5 times more luminous. The objects in our sample are optically classified as Seyfert 1-1.5. This implies that they host a relatively unobscured AGN-like nucleus. They are also of low luminosity and accrete at a low rate. They are therefore good candidates to detect radiation from an inefficient accretion process. However, our results suggest that they are similar to AGN which are 10^2-10^5 times more luminous. The combination of a "radiative efficient accretion disc plus an X-ray producing hot corona" may persist at low accretion rates as well.

Deep optical CCD images of the supernova remnant G 15.1-1.6 were obtained and filamentary and diffuse emission has been discovered. The images, taken in the emission lines of Halpha+[N II], [S II] and [O III], reveal filamentary and diffuse structures all around the remnant. The radio emission at 4850 MHz in the same area is found to be well correlated with the brightest optical filaments. The IRAS 60 micron emission may also be correlated with the optical emission but to a lesser extent. The flux calibrated images suggest that the optical emission originates from shock-heated gas ([S II]/Halpha > 0.4), while there is a possible HII region ([S II]/Halpha ~0.3) contaminating the supernova remnant's emission to the east. Furthermore, deep long-slit spectra were taken at two bright filaments and also show that the emission originates from shock heated gas. An [O III] filamentary structure has also been detected further to the west but it lies outside the remnant's boundaries and possibly is not associated to it. The [O III] flux suggests shock velocities into the interstellar "clouds" ~100 km/s, while the [S II] 6716/6731 ratio indicates electron densities up to ~250 cm^{-3}. Finally, the Halpha emission has been measured to be between 2 to 7 x 10^{-16} erg/s/cm^2/arcsec^2, while the lower limit to the distance is estimated at 2.2 kpc.

Galactic dust constitutes approximately half of the elements more massive than helium produced in stellar nucleosynthesis. Notwithstanding the formation of dust grains in the dense, cool atmospheres of late-type stars, there still remain huge uncertainties concerning the origin and fate of galactic stardust. In this paper, we identify the intergalactic medium (i.e. the region between gravitationally-bound galaxies) as a major sink for galactic dust. We discover a systematic shift in the colour of background galaxies viewed through the intergalactic medium of the nearby M81 group. This reddening coincides with atomic, neutral gas previously detected between the group members. The dust-to-HI mass ratio is high (1/20) compared to that of the solar neighborhood (1/120) suggesting that the dust originates from the centre of one or more of the galaxies in the group. Indeed, M82, which is known to be ejecting dust and gas in a starburst-driven superwind, is cited as the probable main source.

We use the COBE/DIRBE (1.2, 2.2, 60, 100, 140, and 240 $\mu$m) maps and the COBE/FIRAS spectra (for the wavelength range 100 - 1000 $\mu$m) to constrain a model for the spatial distribution of the dust, the stars, and the gas in the Milky Way. By assuming exponential axisymmetric distributions for the dust and the stars and by performing the corresponding radiative transfer calculations we closely (given the simple geometry of the model) reproduce the FIR and NIR maps of the Milky Way. Similar distributions for the atomic and molecular hydrogen in the disk are used (with an inner cut-off radius for the atomic hydrogen) to fit the gas data. The star formation rate as a function of the Galactic radius is derived from the FIR emission and is well in agreement with existing estimates from various star formation tracers. The gas surface density is plotted against the star formation rate density and an ``intrinsic'' Galactic Schmidt law is derived with excellent agreement with the ``external'' Schmidt law found for spiral galaxies. The Milky Way is found to consume $\sim 1%$ and $\sim 10%$ of its gas in the outer and inner regions respectively (for a period of 0.1 Gyr) to make stars. The dust-induced B-V color excess observed in various directions and distances (up to $\sim 6.5$ kpc) with well-studied Cepheid stars is compared with the model predictions showing a good agreement. The simple assumption of exponential distributions of stars and dust in the Galaxy is found to be quite instructive and adequate in modeling all the available data sets from 0.45 $\mu$m (B-band) to 1000 $\mu$m.

The presentation of new results from an [O III] 5007 A survey in a search for planetary nebulae (PNe) in the galactic bulge is continued. A total of 60 objects, including 19 new PNe, have been detected in the remaining 34 per cent of the survey area, while 41 objects are already known. Deep Halpha+[N II] CCD images as well as low resolution spectra have been acquired for these objects. Their spectral signatures suggest that the detected emission originates from photoionized nebulae. In addition, absolute line fluxes have been measured and the electron densities are given. Accurate optical positions and optical diameters are also determined.

The vertical profiles of disc galaxies are built by the material trapped around stable periodic orbits, which form their "skeletons". According to this, the knowledge of the stability of the main families of periodic orbits in appropriate 3D models, can predict possible morphologies for edge-on disc galaxies. In a pilot survey we compare the orbital structures which lead to the appearance of "peanuts" and "X"-like features with the edge-on profiles of three disc galaxies (IC 2531, NGC 4013 and UGC 2048). The subtraction from the images of a model representing the axisymmetric component of the galaxies reveals the contribution of the non-axisymmetric terms. We find a direct correspondence between the orbital profiles of 3D bars in models and the observed main morphological features of the residuals. We also apply a simple unsharp masking technique in order to study the sharpest features of the images. Our basic conclusion is that the morphology of the boxy "bulges" of these galaxies can be explained by considering disc material trapped around stable 3D periodic orbits. In most models these building-blocks periodic orbits are bifurcated from the planar central family of a non-axisymmetric component, usually a bar, at low order vertical resonances. In such a case the boxy "bulges" are parts of bars seen edge-on. For the three galaxies we study the families associated with the "peanut" or "X"-shape morphology are most probably bifurcations at the vertical 2/1 or 4/1 resonance.

We have observed S5 2007+777 and 3C371 in the B and I bands for 13 and 8 nights, respectively, during various observing runs in 2001, 2002 and 2004. The observations resulted in almost evenly sampled light curves, 6-9 hours long. We do not detect any flares within the observed light curves, but we do observe small amplitude, significant variations, in both bands, on time scales of hours and days. The average variability amplitude on time scales of minutes/hours is 2.5% and 1-1.5% in the case of S5 2007+777 and 3C371, respectively. The average amplitudes increase to 5-12% and 4-6%, respectively, on time scales of days. We find that the B and I band variations are highly correlated, on both short and long time scales. During the 2004 observations, which resulted in the longest light curves, we observe two well defined flux-decay and rising trends in the light curves of both objects. When the flux decays, we observe significant delays, with the B band flux decaying faster than the flux in the I band. As a result, we also observe significant, flux related spectral variations as well. The flux-spectral relation is rather complicated, with loop-like structures forming during the flux evolution. The presence of spectral variations imply that the observed variability is not caused by geometric effects. On the other hand, our results are fully consistent with the hypothesis that the observed variations are caused by perturbations which affect different regions in the jet of the sources.

Optical CCD imaging and spectroscopic observations of three supernova remnants are presented. Optical emission from G 54.4-0.3 and G 59.8+1.2 is detected for the first time, while the first flux calibrated CCD images of the supernova remnant G 126.2+1.6 were performed in the optical emission lines of Halpha+[N II], [O III] and [S II]. A mixture of filamentary and diffuse structures is observed in G 54.4-0.3 and G 59.8+1.2, mainly in Halpha+[N II], while the deep optical images of G 126.2+1.6 reveal several new filamentary and diffuse structures inside the extent of the remnant as defined by its known radio emission. In all cases, the radio emission is found to be well correlated with the optical filaments. [O III] emission was not detected at G 54.4-0.3 and G 59.8+1.2 while in G 126.2+1.6, significant morphological differences between the low and medium ionization images are present suggesting incomplete shock structures. Deep long-slit spectra were taken at different positions of the remnants. Both the flux calibrated images and the long-slit spectra clearly show that the emission originates from shock-heated gas, while some spectra of G 126.2+1.6 are characterized by large [O III]/Hbeta ratios. This remnant's [O III] flux suggests shock velocities into the interstellar "clouds" between 100 and 120 km/s, while the [O III] absence in the other two remnants indicates slower shock velocities. For all remnants, the [S II]6716/6731 ratio indicates electron densities below 600 cm^{-3} with particularly low densities for G 54.4-0.3 (below 50 cm^{-3}). Finally, the Halpha emission has been measured to be between 3.0 to 15.2x10^{-17} erg/s/cm^2/ arcsec^2, 3.2x10^{-17} erg/s/cm^2/ arcsec^2 and between 6.5 to 16.8x10^{-17} erg/s/cm^2/ arcsec^2 for G 54.4-0.3, G 59.8+1.2 and G 126.2+1.6, respectively.

We present detailed modeling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared and submm wavelengths. The model takes into account the emission produced by the diffuse dust and the star forming HII complexes. The dust mass is constrained by radiative transfer simulations in the optical (Xilouris et al. 1999). We find that the submm emission predicted by our model cannot account for the observed fluxes. Two scenarios may account for the "missing" submm flux. In the first scenario (Popescu et al. 2000), additional dust (to that derived from the optical, and associated with young stars) is embedded in the galaxy in the form of a thin disk and gives rise to additional submm emission. The other scenario investigates whether the average submm emissivity of the dust grains is higher than the values widely used in Galactic environments. In this case, the dust mass is equal to that derived from the optical observations, and the submm emissivity is treated as a free parameter calculated by fitting our model to the observed SED. We find the submm emissivity value to be ~3 times that often used for our Galaxy. While both scenarios reproduce the observed 850 micron surface brightness, the extra embedded dust model is not supported by the near infrared observations. We, thus, find that the enhanced dust submm emissivity scenario is the most plausible. [abridged]

The extended supernova remnant G 116.5+1.1 was observed in the optical emission lines of Halpha+[N II], [S II] and [O III]; deep long-slit spectra were also obtained. The morphology of the remnant's observed emission is mainly diffuse and patchy in contrast to the known filamentary emission seen along the western limb. The bulk of the detected emission in the region appears unrelated to the remnant but there is one area of emission in the south-east which is characterized by a [S II]/Halpha ratio of ~0.5, implying a possible relation to G 116.5+1.1. If this is actually the case, it would imply a more extended remnant than previously realized. Emission in the [O III] 5007 A line image is not detected, excluding moderate or fast velocity shocks running into ionized interstellar clouds. Our current estimate of the distance to G 116.5+1.1 of ~3 kpc is in agreement with earlier estimates and implies a very extended remnant (69 pc x 45 pc). Observations further to the north-east of G 116.5+1.1 revealed a network of filamentary structures prominent in Halpha+[N II] and [S II] but failed to detect [O III] line emission. Long-slit spectra in a number of positions provide strong evidence that this newly detected emission arises from shock heated gas. Typical Halpha fluxes lie in the range of 9 to 17 x10^{-17} erg/s/cm^2/ arcsec^2, while low electron densities are implied by the intensities of the sulfur lines. Weak emission from the medium ionization line at 5007 A is detected in only one spectrum. Cool dust emission at 60 and 100 microns may be correlated with the optical emission in a limited number of positions. Surpisingly, radio emission is not detected in published surveys suggesting that the new candidate remnant may belong to the class of "radio quiet" supernova remnants.

The dust in galaxies makes radiation transport calculations in them absolutely necessary. It is not only common practice in Astrophysics, but also wisdom, to try to make as simple models as possible to simulate physical systems. For spiral galaxies, however, this turned out to be catastrophic. For years, the major question of the opacity of spiral galaxies kept the community divided, because the models were too simple. A spiral galaxy appears, to first order, to have exponential distributions of stars and dust, which cannot be approximated with uniform distributions. We will review the radiative transfer methods used in galaxies and we will comment on their pluses and minuses. We will also present some of the main results of the application of one of the methods to the observations.

We use the radiation transfer simulation of Xilouris et al. (1999) to constrain the quantity of dust in three nearby spiral galaxies (NGC 4013, NGC 5907 and NGC 891). The predicted visual optical depth from the model is compared with the thermal continuum radiation detected from these galaxies at 850 microns. This comparison yields the emissivity of dust grains in the submillimeter waveband which is a factor 4 higher than the benchmark, semi-empirical model of Draine & Lee (1984). Our estimates are more closely aligned with recent measurements carried out in the laboratory on amorphous carbon and silicate particulates. A comparison between the distribution of 850 microns surface brightness and the intensity levels in the ^{12}CO(1-0) and 21 cm lines underlines the spatial association between dust detected in the submillimeter waveband and molecular gas clouds. We suggest that the relatively high emissivity values that we derive may be attributable to amorphous, fluffy grains situated in denser gas environments.

We investigate the cold and warm dust properties during galaxy interactions using a merging galaxy sample ordered into a chronological sequence from pre- to post-mergers. Our sample comprises a total of 29 merging systems selected to have far-infrared and sub-millimeter observations. We use the 100-to-850 micron flux density ratio, f100/f850, as a proxy to the mass fraction of the warm and the cold dust in these systems. We find evidence for an increase in f100/f850 along the merging sequence from early to advanced mergers and interpret this trend as an increase of the warm relative to the cold dust mass. We argue that the two key parameters affecting the f100/f850 flux ratio is the star-formation rate and the dust content of individual systems relative to the stars. Using a sophisticated model for the absorption and re-emission of the stellar UV radiation by dust we show that these parameters can indeed explain both the increase and the observed scatter in the f100/f850 along the merging galaxy sequence.

We present detailed modelling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared (FIR) and sub-millimeter (submm) wavelengths. The model takes into account the emission of the diffuse dust component, which is heated by the UV and optical radiation fields produced by the stars, as well as the emission produced locally in star forming HII complexes. The radiative transfer simulations of Xilouris et al. (1999) in the optical bands are used to constrain the stellar and dust geometrical parameters, as well as the total amount of dust. We find that the submm emission predicted by our model can not account for the observed fluxes at these wavelengths. We examine two cases, one having more dust embedded in a second thin disk and another allowing for an enhanced submillimeter emissivity of the dust grains. We argue that both cases can equally well reproduce the observed SED. The case of having more dust embedded in a second disk though, is not supported by the near-infrared observations and thus more realistic distributions of the dust (i.e., in spiral arms and clumps) have to be examined in order to better fit the surface brightness of each galaxy.

We investigate the cold and warm dust properties during galaxy interactions using a merging galaxy sample ordered into a chronological sequence from pre- to post-mergers. Our sample comprises a total of 29 merging systems selected to have far-infrared and sub-millimeter observations. The sub-millimeter data are mainly culled from the literature while for 5 galaxies (NGC 3597, NGC 3690, NGC 6090, NGC 6670 and NGC 7252) the sub-millimeter observations are presented here for the first time. We use the 100-to-850 micron flux density ratio, f_{100}/f_{850}, as a proxy to the mass fraction of the warm and the cold dust in these systems. We find evidence for an increase in f_{100}/f_{850} along the merging sequence from early to advanced mergers and interpret this trend as an increase of the warm relative to the cold dust mass. We argue that the two key parameters affecting the f_{100}/f_{850} flux ratio is the star-formation rate and the dust content of individual systems relative to the stars. Using a sophisticated model for the absorption and re-emission of the stellar UV radiation by dust we show that these parameters can indeed explain both the increase and the observed scatter in the f_{100}/f_{850} along the merging galaxy sequence. We also discuss our results under the hypothesis that elliptical galaxies are formed via disc galaxy mergers.

We use the radiation transfer simulation of Xilouris et al. (1999) to constrain the quantity of dust in three nearby spiral galaxies (NGC 4013, NGC 5907 and NGC 4565). The predicted visual optical depth from the model is compared with the thermal continuum radiation detected from NGC 4013 and NGC 5907 at 850 microns and from NGC 4565 at 1.2 mm. The former is based on SCUBA images of NGC 4013 and NGC 5907, reduced and presented for the first time in this work. The comparison of visual optical depth and 850 microns (1.2 mm) emission yields the emissivity of dust grains in the submillimeter (millimeter) waveband. We infer a value of 1.2x10^{4} for the emissivity at 850 microns which is a factor 4 higher than the benchmark, semi-empirical model of Draine & Lee (1984). At 1.2 mm our values are a factor 1.5 higher than this model. Our estimates are more closely aligned with recent measurements carried out in the laboratory on amorphous carbon and silicate particulates. A comparison between the distribution of 850 microns (1.2 mm) surface brightness and the intensity levels in the ^{12}CO(1-0) and 21 cm lines underlines the spatial association between dust detected in the submillimeter/millimeter waveband and molecular gas clouds. We suggest that the relatively high emissivity values that we derive may be attributable to amorphous, fluffy grains situated in denser gas environments.

We present mid-infrared (MIR) maps for a sample of 18 early-type galaxies observed at 4.5, 6.7 and 15 microns with the ISOCAM instrument on board the ISO satellite with a 6'' spatial resolution. We model the Spectral Energy Distribution (SED) of these galaxies using the stellar evolutionary synthesis model PEGASE and we derive the MIR excess over the stellar component. We then explore the nature of this excess in terms of dust and Polycyclic Aromatic Hydrocarbon molecules (PAHs). We find that out of 18 galaxies, 10 show excess at 6.7 microns (due to the presence of PAH features) and 14 show excess at 15 microns (due to the presence of warm dust). In two galaxies, where a more complete wavelength coverage exists, an excess around 9.7 microns is seen (presumably due to silicate dust emission), while two other galaxies are totally devoid of dust. We also examine the morphology of the galaxies in these wavelengths by plotting the azimuthally averaged radial profiles as well as the MIR color profiles. We find that for the majority of the galaxies the 4.5 micron emission is well described by a de Vaucouleurs profile. The 6.7 and 15 micron emission is smoothly distributed throughout the galaxy while only a few galaxies show MIR emission which is more concentrated close to the center. Two dwarf galaxies in our sample show patchy distributions of the MIR emission while two other galaxies show edge-on disks. With color-color diagrams we delineate the regions occupied by late-type and early-type galaxies. Finally we show that the MIR excess found in strong radio galaxies like NGC 4486 (M87) can be explained by synchrotron emission.

We study the morphology of the central region of a sample of Active Galactic Nuclei (AGN) and a ``control'' sample of normal galaxies using archival observations of the WFPC2 instrument onboard the Hubble Space Telescope (HST). We use the ellipse fitting technique in order to get a good description of the inner ``smooth'' light distribution of the galaxy. We then divide the observed galaxy image by the artificial image from the fitted ellipses in order to detect morphological signatures in the central region around the nucleus of the galaxy. We perform quantitative comparisons of different subgroups of our sample of galaxies (according to the Hubble type and the nuclear activity of the galaxies) by calculating the average amplitude of the structures that are revealed with the ellipse fitting technique. Our main conclusions are as follows: 1) All AGNs show significant structure in their inner 100 pc and 1 kpc regions whose amplitude is similar in all of them, independent of the Hubble type of the host galaxy. 2) When considering early-type galaxies, non-AGN galaxies show no structure at all, contrary to what we find for AGN. 3) When considering late-type galaxies, both AGN and non-AGN galaxies show significant structure in their central region. Our results are consistent with the hypothesis that all early-type galaxies host a supermassive black hole, but only those that have enough material in the central regions to fuel it show an active nucleus. The situation is more complicated in late-type galaxies. Either not all of them host a central black hole, or, in some of them, the material inside the innermost 100 pc region is not transported to the scales of the central engine for some reason, or the large amount of gas and dust hides the active nucleus from our sight.

We compare 850um SCUBA images of NGC 891 with the corresponding V-band optical depth predicted from radiation transfer simulations. These two tracers of dust show a very similar distribution along the minor axis and a reasonable agreement along the major axis. Assuming that the grains responsible for optical extinction are also the source of 850um emission we derive a submillimeter emissivity (emission efficiency) for dust in the NGC 891 disk. This quantity is found to be a factor of 2-3 higher than the generally-accepted (but highly uncertain) values adopted for the Milky Way. It should be stated, however, that if a substantial fraction of dust in NGC 891 is clumped, the emissivity in the two galaxies may be quite similar. We use our newly-acquired emissivity to convert our 850um images into detailed maps of dust mass and, utilizing 21cm and CO-emission data for NGC 891, derive the gas-to-dust ratio along the disk. We compute an average ratio of 260 -- a value consistent with the Milky Way and external spirals within the uncertainties in deriving both the dust mass and the quantity of molecular gas. The bulk of dust in NGC 891 appears to be closely associated with the molecular gas phase although it may start to follow the distribution of atomic hydrogen at radii >9 kpc (i.e. >0.5 R_25). Using the optical depth of the NGC 891 disk, we quantify how light emitted at high redshift is attenuated by dust residing in foreground spirals. For B-band observations of galaxies typically found in the Hubble Deep Field, the amount of light lost is expected to be small (~ 5%). This value depends critically on the maximum radial extent of cold dust in spiral disks (which is poorly known). It may also represent a lower limit if galaxies expel dust over time into the intergalactic medium.

We investigate how significant the spiral structure is on calculations concerning radiative transfer in dusty spiral galaxies seen edge-on. The widely adopted exponential disk model (i.e. both the stars and the dust are distributed exponentially in the radial direction and also perpendicular to the plane of the disk) is now subject to a detailed comparison with a realistic model that includes spiral structure for the stars and the dust in the disk. In particular, model images of galaxies with logarithmic spiral arms are constructed, such that the azimuthally averaged disk is exponential in radius and in height, as the observations suggest. Then, pure exponential disk models (i.e. with no spiral structure) are used to fit the edge-on appearance of the model images. As a result, the parameters derived after the fit are compared to the real values used to create the spiral-structured images. It turns out that the plain exponential disk model is able to give a very good description of the galactic disk with its parameters varying only by a few percent from their true values.

The opacity of spiral galaxies is examined by modelling the dust and stellar content of individual galaxies. The model is applied to five late-type spiral galaxies (NGC 4013, IC 2531, UGC 1082, NGC 5529 and NGC 5907). Having analyzed a total of seven galaxies thus far, the five galaxies mentioned above plus UGC 2048 and NGC 891 presented in (Xilouris et al. 1997, 1998), we are able to draw some general conclusions, the most significant of which are: 1) The face-on central optical depth is less than one in all optical bands indicating that typical spiral galaxies like the ones that we have modelled would be completely transparent if they were to be seen face-on. 2) The dust scaleheight is about half that of the stars, which means that the dust is more concentrated near the plane of the disk. 3) The dust scalelength is about 1.4 times larger than that of the stars and the dust is more radially extended than the stars. 4) The dust mass is found to be about an order of a magnitude more than previously measured using the IRAS fluxes, indicating the existence of a cold dust component. The gas-to-dust mass ratio calculated is close to the value derived for our Galaxy. 5) The derived extinction law matches quite well the Galactic extinction law, indicating a universal dust behaviour.