Optical/UV photons and even harder radiation components in galaxies are absorbed and scattered by dust and re-emitted at infrared wavelengths. For a better understanding of the obscured regions of the galaxies detailed models of the interaction of photons with dust grains and the propagation of light are required. A problem which can only be solved by means of numerical solution of the radiative transfer equation. As a prologue we present high angular mid IR observations of galactic nuclei in the spirit of future ELT instrumentation. Dust models are discussed, which are suited to fit the extinction curves and relevant to compute the emission of external galaxies. Self-consistent radiative transfer models have been presented in spherical symmetry for starburst nuclei, in two dimensions for disk galaxies (spirals) and, more recently, in three dimensional configuration of the dust density distribution. For the latter, a highlighting example is the clumpy dust tori around AGN. Modern advances in the field are reviewed which are either based on a more detailed physical picture or progress in computational sciences.

The Sombrero galaxy (M104) is an interesting object for a dust energy balance study due to its very symmetric dust lane, its proximity and its (nearly edge-on) inclination of 84°. From a panchromatic radiative transfer analysis, including scattering, absorption and thermal dust re-emission, we construct a standard model for M104 accounting for observations in the optical wave bands (stellar SED, images and extinction profiles in the V and RC band). This standard model underestimates the observed dust emission at infrared wavelengths by a factor of ~ 3, similar to the discrepancy found in other energy balance studies of edge-on spirals. Supplementing this standard model with a young stellar component of low star formation activity in both the inner disk (SFR ~ 0.21 M⊙ yr−1) and dust ring (SFR ~ 0.05 M⊙ yr−1), we are capable of solving the discrepancy in the dust energy budget of the Sombrero galaxy at wavelengths shortwards of 100 μm. To account for the remaining discrepancy at longer wavelengths, we propose a secondary dust component distributed in quiescent clumps. This model with a clumpy dust structure predicts three-quarters of the total dust content to reside in compact dust clouds with no associated embedded sources.

We use Herschel Space Observatory and Spitzer Space Telescope 70-500 μm data along with ground-based optical and near-infrared data to understand how dust heating in the nearby face-on spiral galaxies M81, M83, and NGC 2403 is affected by the starlight from all stars and by the radiation from star-forming regions. We find that 70/160 μm flux density ratios tend to be more strongly influenced by star-forming regions. However, the 250/350 and 350/500 μm micron flux density ratios are more strongly affected by the light from the total stellar populations, suggesting that the dust emission at > 250 μm originates predominantly from a component that is colder than the dust seen at <160 μm and that is relatively unaffected by star formation activity. We conclude by discussing the implications of this for modelling the spectral energy distributions of both nearby and more distant galaxies and for using far-infrared dust emission to trace star formation.

Late-type low surface brightness (LSB) disk galaxies are common in the local universe and appear dynamically and chemically under evolved compared to their high surface brightness (HSB) counterparts. We have utilized multi-wavelength imaging and photometry of three edge-on, low-mass LSB disk galaxies to investigate the dust distribution in such systems. Through the use of Monte Carlo radiation transfer models to interpret the data, we find that the dust disk appears to have a vertical scale height similar to the stellar disk. This is in contrast to previous findings from HSB galaxies, where the dust is believed to be more concentrated in the galactic mid-plane. We believe the change in the relative scale heights of the dust and stellar disks is likely associated with the increased stability of the ISM against vertical collapse and the thin nature of the stellar disks.

We investigate the three-dimensional structure of the nearby edge-on spiral galaxy NGC 891 using 3D Monte Carlo radiative transfer models, with realistic spiral structure and fractally clumped dust. Using the spiral and clumpiness parameters found from recently completed scattered light models we produce lower resolution SED models which reproduce the spatially-integrated UV-to-FIR SED of NGC 891. Our models contain a color gradient across the major axis of the galaxy - similar to what is seen in images of the NGC 891. With minor adjustment our SED models are able to match the SED of M51, a similar galaxy at a near face-on inclination.

Some 300-500 Myr ago, the Whirlpool galaxy (NGC 5194/M51a) and its nearby post-starburst galaxy neighbour, NGC 5195/M51b closely interacted, resulting in significant changes to their star formation activity. Both galaxies display colors indicative of enhanced star formation during closest passage, but since then, star formation has ceased in NGC 5195 yet remained ongoing in the spiral NGC 5194. With a wealth of multi-wavelength (0.2–500 μm for this study) observations available, this nearby (10 Mpc) system, whose star formation history is well constrained through optical colors of individual stars and its dynamical history, provides the optimal laboratory to test the relation between dust emission and stellar emission within the fundamental framework of today's stellar population synthesis and dust emission models.

Due to its proximity, the Andromeda galaxy (M31, NGC 224) offers a unique insight into how the spectra of stars, dust, and gas combine to form the integrated Spectral Energy Distribution (SED) of galaxies. We introduce here Herschel Space Observatory PACS and SPIRE photometric observations of M31 which cover the far-infrared to sub-mm wavelengths (70-500 μm). These new observations reveal that the total IR luminosity of M31 is relatively weak, with LIR=109.65L⊙, only 10% of the total luminosity of M31. However, as seen in the previous studies of M31, the IR luminosity is dominated by a 10 kpc ring in all Herschel bands. This is distinct from the optical, where the bulge in the central 2kpc, dominates the luminosity, clearly demonstrating how different components at distinct positions in a galaxy contribute to make the integrated SED.

A large body of evidence has demonstrated that the global rest-frame optical and IR colours of galaxies correlate well with each other, which can be readily interpreted as a sign of typically smooth star formation histories. However the processes that lead to the observed correlations are contrary: the stellar light that contributes to the optical is readily absorbed by dust which emits in the IR. Thus on small scales we expect these correlations to break down. In this contribution we present our recent results (Zibetti & Groves 2011) from a pixel-by-pixel multi-wavelength (u-band to 8μm) analysis of seven nearby galaxies ranging from early- to late-types. We show that such a break-down occurs already on scales on few 100 pc, as a result of the different physical conditions in spatially distinct regions inside the galaxy, as we demonstrate by means of a Principal Component Analysis. Despite the lack of internal correlation between optical and IR within individual galaxies, when the pixels of all galaxies are compared the well known optical-IR colour correlations return, demonstrating that the variance observed within galaxies is limited around a mean which follows the well-known trends. We also examine the extremely strong correlations between the mid IR (Spitzer-IRAC)-NIR colours which extend continuously across all galaxies. These correlations arise from the differing contribution of stellar light and dust to the IRAC bands, enabling us to determine pure stellar colours for these bands, but still demonstrating the need for dust (or stellar) corrections in these bands when being used as stellar (dust) tracers.

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 H ii regions in the Local Group Galaxy M33. Using the Hα emission, we perform a classification of a selected H ii region sample in terms of morphology, separating the objects in filled, mixed, shell and clear shell objects. We obtain the SED for each H ii 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.

Using the Herschel PACS and SPIRE FIR/submm data, we investigate variations in the dust spectral index β in the nearby spiral galaxy M33 at a linear resolution of 160 pc. We use an iteration method in two different approaches, single and two-component modified black body models. In both approaches, β is higher in the central disk than in the outer disk similar to the dust temperature. There is a positive correlation between β and Hα as well as with the molecular gas traced by CO(2-1). A Monte-Carlo simulation shows that the physical parameters are better constrained when using the two-component model.

Edge-on spiral galaxies offer a unique perspective on the vertical structure of spiral disks, both stars and the iconic dark dustlanes. The thickness of these dustlanes can now be resolved for the first time with Herschel in far-infrared and sub-mm emission. We present NHEMESES, an ongoing project that targets 12 edge-on spiral galaxies with the PACS and SPIRE instruments on Herschel. These vertically resolved observations of edge-on spirals will impact on several current topics.

First and foremost, these Herschel observations will settle whether or not there is a phase change in the vertical structure of the ISM with disk mass. Previously, a dramatic change in dustlane morphology was observed as in massive disks the dust collapses into a thin lane. If this is the case, the vertical balance between turbulence and gravity dictates the ISM structure and consequently star-formation and related phenomena (spiral arms, bars etc.). We specifically target lower mass nearby edge-ons to complement existing Herschel observations of high-mass edge-on spirals (the HEROES project).

Secondly, the combined data-set, together with existing Spitzer observations, will drive a new generation of spiral disk Spectral Energy Distribution models. These model how dust reprocesses starlight to thermal emission but the dust geometry remains the critical unknown.

And thirdly, the observations will provide an accurate and unbiased census of the cold dusty structures occasionally seen extending out of the plane of the disk, when backlit by the stellar disk. To illustrate the NHEMESES project, we present early results on NGC 4244 and NGC 891, two well studies examples of a low and high-mass edge-on spiral.

We revisit the question of the ionization of the diffuse medium in late type galaxies, by studying NGC 891. The most important challenge for the models considered so far was the observed increase of [O iii]/Hβ, [O iii]/Hβ, and [N ii]\Hα with increasing distance to the galactic plane. We propose a scenario based on the expected population of massive OB stars and hot low-mass evolved stars (HOLMES) in this galaxy to explain this observational fact. In the framework of this scenario we construct a finely meshed grid of photoionization models. For each value of the galactic latitude z we look for the models which simultaneously fit the observed values of the [O iii]/Hβ, [O ii]/Hβ, and [N ii]\Hα ratios. For each value of z we find a range of solutions which depends on the value of the oxygen abundance. The models which fit the observations indicate a systematic decrease of the electron density with increasing z. They become dominated by the HOLMES with increasing z only when restricting to solar oxygen abundance models, which argues that the metallicity above the galactic plane should be close to solar. They also indicate that N/O increases with increasing z.

Significant discrepancies have been found between the dust masses derived from various tracers (optical/near-IR, far-IR/sub-millimeter observations, and the variation of dust attenuation with viewing angle). Here we report the first detection of the extended far-UV (FUV) and near-UV (NUV) haloes perpendicular to the galactic plane of NGC 891, which can be interpreted as scattered stellar light from the galactic plane. An additional “geometrically thick” dust disk, which contains about the same mass as the standard thin dust disk, is needed to reproduce the vertically extended UV profiles

UGC 12281 has been classified as having a pure disk and being a low surface brightness galaxy (LSBG), thus being an obvious member of the so-called superthin galaxies. At the same time it represents an extremely untypical type of LSBG due to its remarkable amount of current star formation and evidence for extraplanar ionized gas. This makes it become a perfect tool to investigate the triggering of star formation in LSB galaxies, located in an alleged isolated area. By means of deep photometry and long-slit spectroscopy we analyse the Hα halo and verify the existence of a potential dwarf companion which we found on processed SDSS images.

The low metallicity interstellar medium of dwarf galaxies gives a different picture in the far infrared(FIR)/submillimetre(submm)wavelengths than the more metal-rich galaxies. Excess emission is often found in the submm beginning at or beyond 500 μm. Even without taking this excess emission into account as a possible dust component, higher dust-to-gas mass ratios (DGR) are often observed compared to that expected from their metallicity for moderately metal-poor galaxies. The Spectral Energy Distributions (SEDs) of the lowest metallicity galaxies, however, give very low dust masses and excessively low values of DGR, inconsistent with the amount of metals expected to be captured into dust if we presume the usual linear relationship holding for all metallicities, including the more metal-rich galaxies. This transition seems to appear near metalllicities of 12 + log(O/H) 8.0 - 8.2. These results rely on accurately quantifying the total molecular gas reservoir, which is uncertain in low metallicity galaxies due to the difficulty in detecting CO(1-0) emission. Dwarf galaxies show an exceptionally high [CII] 158 μm/CO (1-0) ratio which may be indicative of a significant reservoir of ‘CO-free’ molecular gas residing in the photodissociated envelope, and not traced by the small CO cores.

The Herschel Space Observatory is revolutionizing our view of dust in galaxies with high sensitivity observations in the far infrared(FIR)/submillimeter (submm) regime from 70 to 500 μm. Herschel is confirming the submm excess that has been noted previously in low metallicity dwarfs. We present here the Dwarf Galaxies Survey sample through a Herschel colour-colour diagram. We will then focus on two galaxies, Haro11 and NGC4449 presenting different interesting behaviours in the FIR as revealed by Herschel.

We have carried out a detailed modelling of the dust heating and emission in the nearby, starbursting dwarf galaxy NGC 4214. Due to its proximity and the great wealth of data from the UV to the millimeter range (from GALEX, HST, Spitzer, Herschel, Planck and IRAM) it is possible to separately model the emission from HII regions and their associated photodissociation regions (PDRs) and the emission from diffuse dust. Furthermore, most model parameters can be directly determined from the data leaving very few free parameters. We can fit both the emission from HII+PDR regions and the diffuse emission in NGC 4214 with these models with “normal” dust properties and realistic parameters.

We have carried out a detailed modeling of the dust Spectral Energy Distribution (SED) of the nearby, starbursting dwarf galaxy NGC 4214. A key point of our modeling is that we distinguish the emission from (i) HII regions and their associated photodissociation regions (PDRs) and (ii) diffuse dust. For both components we apply templates from the literature calculated with a realistic geometry and including radiation transfer. The large amount of existing data from the ultraviolet (UV) to the radio allows the direct measurement of most of the input parameters of the models. We achieve a good fit for the total dust SED of NGC 4214. In the present contribution we describe the available data, the data reduction and the determination of the model parameters, whereas a description of the general outline of our work is presented in the contribution of Lisenfeld et al. in this volume.

We use the photoionisation and dust radiative transfer code MOCASSIN to create a model of the dwarf irregular galaxy NGC 4449. The best-matching model reproduces the global optical emission line fluxes and the observed spectral energy distribution (SED) spanning wavelengths from the UV to sub-mm, and requires the bolometric luminosity of 6.25 × 109L⊙ for the underlying stellar component, Mdust/Mgas of 1/680 and Mdust of 2.2 × 106M⊙.

The SAGE-LMC, SAGE-SMC and HERITAGE surveys have mapped the Magellanic Clouds in the infrared using the Spitzer and Herschel Space Telescopes. Over 8.5 million point sources were detected and catalogued in the LMC alone. Staring mode observations using the InfraRed Spectrograph (IRS) on board Spitzer have been obtained for 1,000 positions in the LMC and ~250 in the SMC. From the infrared spectroscopy we have identified the nature of the sources for which spectroscopy is available. These IRS staring mode targets represent an important contribution to the SED of these dwarf galaxies. Here we report on our latest results.