We present an extensively updated version of the purely ray-tracing 3D dust radiation transfer code DART-Ray. The new version includes five major upgrades : 1) a series of optimizations for the ray-angular density and the scattered radiation source function; 2) the implementation of several data and task parallelizations using hybrid MPI+OpenMP schemes; 3) the inclusion of dust self-heating; 4) the ability to produce surface brightness maps for observers within the models in HEALPix format; 5) the possibility to set the expected numerical accuracy already at the start of the calculation. We tested the updated code with benchmark models where the dust self-heating is not negligible. Furthermore, we performed a study of the extent of the source influence volumes, using galaxy models, which are critical in determining the efficiency of the DART-Ray algorithm. The new code is publicly available, documented for both users and developers, and accompanied by several programmes to create input grids for different model geometries and to import the results of N-body and SPH simulations. These programmes can be easily adapted to different input geometries, and for different dust models or stellar emission libraries.

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

We describe the calculation of the stochastically heated dust emission using the 3D ray-tracing dust radiative transfer code DART-Ray, which is designed to solve the dust radiative transfer problem for galaxies with arbitrary geometries. In order to reduce the time required to derive the non-equilibrium dust emission spectra from each volume element within a model, we implemented an adaptive SED library approach, which we tested for the case of axisymmetric galaxy geometries. To show the capabilities of the code, we applied DART-Ray to a high-resolution N-body+SPH galaxy simulation to predict the appearance of the simulated galaxy at a set of wavelengths from the UV to the sub-mm. We analyse the results to determine the effect of dust on the observed radial and vertical profiles of the stellar emission as well as on the attenuation and scattering of light from the constituent stellar populations. We also quantify the proportion of dust re-radiated stellar light powered by young and old stellar populations, both bolometrically and as a function of infrared wavelength.

The probability distribution function (PDF) of the mass surface density of molecular clouds provides essential information about the structure of molecular cloud gas and condensed structures out of which stars may form. In general, the PDF shows two basic components: a broad distribution around the maximum with resemblance to a log-normal function, and a tail at high mass surface densities attributed to turbulence and self-gravity. In a previous paper, the PDF of condensed structures has been analyzed and an analytical formula presented based on a truncated radial density profile, $\rho(r) = \rho_c/(1+(r/r_0)^2)^{n/2}$ with central density $\rho_c$ and inner radius $r_0$, widely used in astrophysics as a generalization of physical density profiles. In this paper, the results are applied to analyze the PDF of self-gravitating, isothermal, pressurized, spherical (Bonnor-Ebert spheres) and cylindrical condensed structures with emphasis on the dependence of the PDF on the external pressure $p_{ext}$ and on the overpressure $q^{-1} =p_c /p_{ext}$, where $p_c$ is the central pressure. Apart from individual clouds, we also consider ensembles of spheres or cylinders, where effects caused by a variation of pressure ratio, a distribution of condensed cores within a turbulent gas, and (in case of cylinders) a distribution of inclination angles on the mean PDF are analyzed. The probability distribution of pressure ratios $q^{-1}$ is assumed to be given by $P(q^{-1}) \propto q^{-k_1}/(1+(q_0/q)^{\gamma})^{(k_1+k_2)/{\gamma}}$, where $k_1$, ${\gamma}$, $k_2$, and $q_0$ are fixed parameters.

The probability distribution function (PDF) of the mass surface density is an essential characteristic of the structure of molecular clouds or the interstellar medium in general. Observations of the PDF of molecular clouds indicate a composition of a broad distribution around the maximum and a decreasing tail at high mass surface densities. The first component is attributed to the random distribution of gas which is modeled using a log-normal function while the second component is attributed to condensed structures modeled using a simple power-law. The aim of this paper is to provide an analytical model of the PDF of condensed structures which can be used by observers to extract information about the condensations. The condensed structures are considered to be either spheres or cylinders with a truncated radial density profile at cloud radius r_cl. The assumed profile is of the form rho(r)=rho_c/(1+(r/r_0)^2)^{n/2} for arbitrary power n where rho_c and r_0 are the central density and the inner radius, respectively. An implicit function is obtained which either truncates (sphere) or has a pole (cylinder) at maximal mass surface density. The PDF of spherical condensations and the asymptotic PDF of cylinders in the limit of infinite overdensity rho_c/rho(r_cl) flattens for steeper density profiles and has a power law asymptote at low and high mass surface densities and a well defined maximum. The power index of the asymptote Sigma^(-gamma) of the logarithmic PDF (Sigma x P(Sigma)) in the limit of high mass surface densities is given by gamma = (n+1)/(n-1)-1 (spheres) or by gamma=n/(n-1)-1 (cylinders in the limit of infinite overdensity).

We investigate star formation and dust heating in the compact FIR bright sources detected in the Herschel maps of M83. We use the source extraction code GETSOURCES to detect and extract sources in the FIR, as well as their photometry in the MIR and H{\alpha}. By performing infrared SED fitting and applying an H{\alpha} based star formation rate (SFR) calibration, we derive the dust masses and temperatures, SFRs, gas masses and star formation efficiencies (SFEs). The detected sources lie exclusively on the spiral arms and represent giant molecular associations (GMAs), with gas masses and sizes of 10^6-10^8 solar masses and 200-300 pc, respectively. The inferred parameters show little to no radial dependence and there is only a weak correlation between the SFRs and gas masses, which suggests that more massive clouds are less efficient at forming stars. Dust heating is mainly due to local star formation. However, although the sources are not optically thick, the total intrinsic young stellar population luminosity can almost completely account for the dust luminosity. This suggests that other radiation sources contribute to the dust heating as well and approximately compensate for the unabsorbed fraction of UV light.

We present a detailed study of the surface brightness profiles of dense filaments in IC 5146 using recent Herschel observations done with SPIRE. We describe the profile through an equilibrium solution of a self-gravitating isothermal cylinder pressure confined by its surrounding medium. In this first analysis we applied a simple modified black body function for the emissivity, neglecting any radiative transfer effects. Overall we found a good agreement of the observed surface brightness profiles with the model. The filaments indicate strong self-gravity with mass line densities M/l\gtrsim 0.5 (M/l)_max where (M/l)_max is the maximum possible mass line density. In accordance with the model expectations we found a systematic decrease of the FWHM, a steepening of the density profile, and for filaments heated by the interstellar radiation field a decrease of the luminosity to mass ratio for higher central column density and mass line density. We illustrate and discuss the possibility of estimating the distance, external pressure, and dust opacity. For a cloud distance D\sim500 pc and a gas temperature of T_cyl=10 K the model implies an external pressure p_ext/k\sim 2x10^4 K cm^-3 and an effective dust emission coefficient at 250 microns given by delta x kappa_0^em \sim 0.0588 cm^2 g^-1 where delta is the dust-to-gas ratio. Given the largest estimate of the distance to the cloud complex, 1 kpc, the model yields an upper limit delta x kappa_0^em \sim 0.12 cm^2 g^-1.

We analyze the physical parameters of interstellar filaments that we describe by an idealized model of isothermal self-gravitating infinite cylinder in pressure equilibrium with the ambient medium. Their gravitational state is characterized by the ratio f_cyl of their mass line density to the maximum possible value for a cylinder in a vacuum. Equilibrium solutions exist only for f_cyl < 1. This ratio is used in providing analytical expressions for the central density, the radius, the profile of the column density, the column density through the cloud centre, and the fwhm. The dependence of the physical properties on external pressure and temperature is discussed and directly compared to the case of pressure-confined isothermal self-gravitating spheres. Comparison with recent observations of the fwhm and the central column density N_H(0) show good agreement and suggest a filament temperature of ~10 K and an external pressure p_ext/k in the range 1.5x10^4 K/cm^3 to 5x10^4 K/cm^3. Stability considerations indicate that interstellar filaments become increasingly gravitationally unstable with mass line ratio f_cyl approaching unity. For intermediate f_cyl>0.5 the instabilities should promote core formation through compression, with a separation of about five times the fwhm. We discuss the nature of filaments with high mass line densities and their relevance to gravitational fragmentation and star formation.

We investigated the star formation efficiency for all the dust emitting sources in Stephan's Quintet (SQ). We inferred star formation rates using Spitzer MIR/FIR and GALEX FUV data and combined them with gas column density measurements by various authors, in order to position each source in a Kennicutt-Schmidt diagram. Our results show that the bright IGM star formation regions in SQ present star formation efficiencies consistent with those observed within local galaxies. On the other hand, star formation in the intergalactic shock region seems to be rather inhibited.

We analysed the Spitzer maps of Stephan's Quintet in order to investigate the nature of the dust emission associated with the X-ray emitting regions of the large scale intergalactic shock and of the group halo. This emission can in principle be powered by dust-gas particle collisions, thus providing efficient cooling of the hot gas. However the results of our analysis suggest that the dust emission from those regions is mostly powered by photons. Nonetheless dust collisional heating could be important in determining the cooling of the IGM gas and the large scale star formation morphology observed in SQ.

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 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, {\it 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.

The aim of the paper is to give a physical explanation of the absence of the feature in the Calzetti extinction curve. We analyze the dust attenuation of a homogeneous source seen through a distant inhomogeneous distant screen. The inhomogeneities are described through an idealized isothermal turbulent medium where the probability distribution function (PDF) of the column density is log-normal. In addition it is assumed that below a certain critical column density the carriers of the extinction bump at 2175 Angstroem are being destroyed by the ambient UV radiation field. Turbulence is found to be a natural explanation not only of the flatter curvature of the Calzetti extinction curve but also of the missing bump provided the critical column density is N_H >= 10^21 cm^-2. The density contrast needed to explain both characteristics is well consistent with the Mach number of the cold neutral medium of our own Galaxy which suggests a density contrast sigma_(rho/<rho>) 6.

The dust emission spectrum and the brightness profile of passively heated condensed cores is analyzed in relation to their astrophysical environment. The cores are modeled as critically stable self-gravitating spheres embedded at the center of self-gravitating filaments that are assumed to be either spherical or cylindrical in shape. The filaments are heated by an isotropic interstellar radiation field (ISRF). The calculations are based on a physical dust model of stochastically heated grains of diffuse interstellar dust. The spectral energy distribution (SED) of the cores is calculated using a ray-tracing technique where the effects of scattered emission and re-heating by dust grains are accurately taken into account. To compare with observational studies, the dust re-emission spectrum is approximated by a modified black-body function and the brightness profile with a Gaussian source. A simplified single-zone model for cores is presented that incorporates on the basis of the derived emissivities a first order approximation of their SED. Colder dust temperatures are, independent of the core mass, related to a higher pressure both inside and around the filament. The pressure-temperature relation for given external pressure is found to be largely independent of the true shape of the filament. The calculations show that the radiative transfer leads to a lower emission coefficient at 250 micron and to a flatter emissivity law of typically beta<1.8 in the far-infrared sub-millimeter regime. These effects cause the core mass to be underestimated by more than a factor of 2 based on the typical assumptions used in observational programs. A larger uncertainty is expected for high pressure regions.

We present a self-consistent model of the spectral energy distributions (SEDs) of spiral galaxies from the ultraviolet (UV) to the mid-infrared (MIR)/far-infrared (FIR)/submillimeter (submm) based on a full radiative transfer calculation of the propagation of starlight in galaxy disks. This model predicts not only the total integrated energy absorbed in the UV/optical and re-emitted in the infrared/submm, but also the colours of the dust emission based on an explicit calculation of the strength and colour of the UV/optical radiation fields heating the dust, and incorporating a full calculation of the stochastic heating of small dust grains and PAH molecules. The geometry of the translucent components of the model is empirically constrained using the results from the radiation transfer analysis of Xilouris et al. on spirals in the middle range of the Hubble sequence, while the geometry of the optically thick components is constrained from physical considerations with a posteriori checks of the model predictions with observational data. These geometrical constraints enable the dust emission to be predicted in terms of a minimum set of free parameters: the central face-on dust opacity in the B-band tau^f_B, a clumpiness factor F for the star-forming regions, the star-formation rate SFR, the normalised luminosity of the old stellar population old and the bulge-to-disk ratio B/D. We show that these parameters are almost orthogonal in their predicted effect on the colours of the dust/PAH emission. The results of the calculations are made available in the form of a large library of simulated dust emission SEDs spanning the whole parameter space of our model, together with the corresponding library of dust attenuation calculated using the same model. (see full abstract in the paper)

We analyse a comprehensive set of MIR/FIR observations of Stephan's Quintet (SQ), taken with the Spitzer Space Observatory. Our study reveals the presence of a luminous (L_{IR}\approx 4.6x10^43 erg/s) and extended component of infrared dust emission, not connected with the main bodies of the galaxies, but roughly coincident with the X-ray halo of the group. We fitted the inferred dust emission spectral energy distribution of this extended source and the other main infrared emission components of SQ, including the intergalactic shock, to elucidate the mechanisms powering the dust and PAH emission, taking into account collisional heating by the plasma and heating through UV and optical photons. Combining the inferred direct and dust-processed UV emission to estimate the star formation rate (SFR) for each source we obtain a total SFR for SQ of 7.5 M(sun)/yr, similar to that expected for non-interacting galaxies with stellar mass comparable to the SQ galaxies. Although star formation in SQ is mainly occurring at, or external to the periphery of the galaxies, the relation of SFR per unit physical area to gas column density for the brightest sources is similar to that seen for star-formation regions in galactic disks. We also show that available sources of dust in the group halo can provide enough dust to produce up to L_{IR}\approx 10^42 erg/s powered by collisional heating. Though a minority of the total infrared emission (which we infer to trace distributed star-formation), this is several times higher than the X-ray luminosity of the halo, so could indicate an important cooling mechanism for the hot IGM and account for the overall correspondence between FIR and X-ray emission.

We derive the spectral energy distribution (SED) of dusty, isothermal, self gravitating, stable and spherical clouds externally heated by the ambient interstellar radiation field. For a given radiation field and dust properties, the radiative transfer problem is determined by the pressure of the surrounding medium and the cloud mass expressed as a fraction of the maximum stable cloud mass above which the clouds become gravitational unstable. To solve the radiative transfer problem a ray-tracing code is used to accurately derive the light distribution inside the cloud. This code considers both non isotropic scattering on dust grains and multiple scattering events. The dust properties inside the clouds are assumed to be the same as in the diffuse interstellar medium in our galaxy. We analyse the effect of the pressure, the critical mass fraction, and the ISRF on the SED and present brightness profiles in the visible, the IR/FIR and the submm/mm regime with the focus on the scattered emission and the thermal emission from PAH-molecules and dust grains.

We combine the the stellar spectral synthesis code Starburst99, the nebular modelling code MAPPINGSIII, and a 1-D dynamical evolution model of HII regions around massive clusters of young stars to generate improved models of the spectral energy distribution (SED) of starburst galaxies. We introduce a compactness parameter, C, which characterizes the specific intensity of the radiation field at ionization fronts in HII regions, and which controls the shape of the far-IR dust re-emission, often referred to loosely as the dust ``temperature''. We also investigate the effect of metallicity on the overall SED and in particular, on the strength of the PAH features. We provide templates for the mean emission produced by the young compact HII regions, the older (10 - 100 Myr) stars and for the wavelength-dependent attenuation produced by a foreground screen of the dust used in our model. We demonstrate that these components may be combined to produce a excellent fit to the observed SEDs of star formation dominated galaxies which are often used as templates (Arp 220 and NGC 6240). This fit extends from the Lyman Limit to wavelengths of about one mm. The methods presented in both this paper and in the previous papers of this series allow the extraction of the physical parameters of the starburst region (star formation rates, star formation rate history, mean cluster mass, metallicity, dust attenuation and pressure) from the analysis of the pan-spectral SED.

We build, as far as theory will permit, self consistent model HII regions around central clusters of aging stars. These produce strong emission line diagnostics applicable to either individual HII regions in galaxies, or to the integrated emission line spectra of disk or starburst galaxies. The models assume that the expansion and internal pressure of individual HII regions is driven by the net input of mechanical energy from the central cluster, be it through winds or supernova events. This eliminates the ionization parameter as a free variable, replacing it with a parameter which depends on the ratio of the cluster mass to the pressure in the surrounding interstellar medium. These models explain why HII regions with low abundances have high excitation, and demonstrate that at least part of the warm ionized medium is the result of overlapping faint, old, large, and low pressure HII regions. We present line ratios (at both optical and IR wavelengths) which provide reliable abundance diagnostics for both single HII regions or for integrated galaxy spectra, and we find a number that can be used to estimate the mean age of the cluster stars exciting individual HII regions.

We examine, from a theoretical viewpoint, how the physical parameters of HII regions are controlled in both normal galaxies and in starburst environments. These parameters are the HII region luminosity function, the time-dependent size, the covering fraction of molecular clouds, the pressure in the ionized gas and the ionization parameter. The factors which control them are the initial mass function of the exciting stars, the cluster mass function, the metallicity and the mean pressure in the surrounding interstellar medium. We investigate the sensitivity of the H$\alpha$ luminosity to the IMF, and find that this can translate to more than a factor two variation in derived star formation rates. The molecular cloud dissipation timescale is estimated from a case study of M17 to be $\sim1$ Myr for this object. Based upon HII luminosity function fitting for nearby galaxies, we suggest that the \HII region cluster mass function is fitted by a log-normal form peaking at $\sim 100 M_{\odot}$. The cluster mass function continues the stellar IMF to higher mass. The pressure in the HII regions is controlled by the mechanical luminosity flux from the central cluster. Since this is closely related to the ionizing photon flux, we show that the ionization parameter is not a free variable, and that the diffuse ionized medium may be composed of many large, faint and old HII regions. Finally, we derive theoretical probability distributions for the ionization parameter as a function of metallicity and compare these to those derived for SDSS galaxies.

The optical properties of small spherical iron grains are derived using a Kramers-Kronig-consistent model of the dielectric function including its dependence on temperature and size. Especially discussed is the effect of the size dependence, which results from the limitation of the free path of the free electrons in the metal by the size of the grain, on the absorption behaviour of small iron spheres and spheroids. The estimated absorption properties are applied to study the temperature behaviour of spherical and spheroidal grains which are heated by the interstellar radiation field.

We analyse in detail the attenuation caused by an isothermal turbulent distant foreground dust screen. The attenuation curve is well determined by two parameters, the absolute-to-relative attenuation ratio R^A_V=A_V/E(B-V) and the absolute attenuation A_V. We show quantitatively how these two observable quantities depend on the statistical properties of the local density and the mean attenuation <A_V> and how they vary with the thickness of the screen measured in units of the largest turbulent scale. The attenuation through a turbulent medium is characterised by higher transparency and a flatter attenuation curve in comparison with a homogeneous dust screen. In general, the effect of the turbulent medium on the attenuation increases with slice thickness. In the limit of thick slices, typically larger than the maximum turbulent scale, R^A_V asymptotically approaches a maximum value.

Our ability to correct the observational photometry of galaxies depends upon our knowledge of the attenuation of light produced by the dust contained in the interstellar medium. We will present a model based on the statistical properties of the local density in isothermal turbulence which might be appropriate to calculate the radiative transport through the diffuse interstellar medium. The model will be applied to study the attenuation caused by a distant foreground screen and a non scattering slab where the sources are mixed within the turbulent medium. It will be shown how the turbulent structure affects the attenuation curve and how the attenuation curve varies with the thickness of a foreground screen.

In this paper, we combine the stellar spectral synthesis code STARBURST 99, the nebular modelling code MAPPINGS IIIq, a 1-D dynamical evolution model of \HII regions around massive clusters of young stars and a simplified model of synchrotron emissivity to produce purely theoretical self-consistent synthetic spectral energy distributions (SEDs) for (solar metallicity) starbursts lasting some $10^8$ years. These SEDs extend from the Lyman Limit to beyond 21 cm. We find that two ISM parameters control the form of the SED; the pressure in the diffuse phase of the ISM (or, equivalently, its density), and the molecular cloud dissipation timescale. We present detailed SED fits to Arp 220 and NGC 6240, and we give the predicted colors for starburst galaxies derived from our models for the IRAS and the Spitzer Space Observatory MIPS and IRAC instruments. Our models reproduce the spread in observed colors of starburst galaxies. Finally, we present absolute calibrations to convert observed fluxes into star formation rates in the UV (GALEX), at optical wavelengths (H$\alpha$), and in the IR (IRAS or the Spitzer Space Observatory). (Abstract Truncated)

We present measurements of the mid-infrared (MIR) emission from SN 1987A, made using the Infrared Space Observatory (ISO) 11 years after outburst. They are the only late epoch detections of this source in the thermal IR regime. The position of the source, determined from an offset to an IR-emitting star, suggests that the emission is associated with SN 1987A or its extended supernova remnant (SNR). A predominantly circumstellar origin is however suggested by the size and orientation of the IR-emitting region, which is comparable with the extension of the inner ring seen with the Hubble Space Telescope (HST). The emission is most probably from collisionally-heated circumstellar grains embedded in shocked gas downstream of the blast wave. The MIR extent is consistent with the hypothesis that the blast wave was propagating into material of moderate density interior to the thick inner ring at the epoch of the ISOCAM observations.

Detailed models are presented for the late epoch mid infrared (MIR) emission from collisionally heated grains in the shocked circumstellar gas around SN 1987A. Thermal dust emission from a region of moderate density interior to the thick inner ring seen with the Hubble Space Telescope (HST) is found to be a natural explanation for the MIR spectral energy distribution measured by ISOCAM. The MIR-spectrum can be reproduced by a mixture of silicate-iron or silicate-graphite grains or by a composition of pure graphite grains. A composition of pure iron grains on the other hand can be excluded and a pure silicate composition does not seem to be very likely. The dust-to-gas ratio in the interaction zone is ~0.01%, an order of magnitude lower than estimates for dust abundances in the winds of red supergiant (RSG) stars in the LMC. This low dust abundance can be accounted for by a combination of evaporation through the UV-flash from the supernova outburst and subsequent sputtering in the shocked gas. For this explanation to hold, dust in the pre-supernova circumstellar medium (CSM) would have to have been predominantly composed of grains other than graphite, with a maximum size smaller than ~0.1 microns.

The empirical reddening function for starburst galaxies generated by Calzetti and her co-workers has proven very successful, and is now used widely in the observational literature. Despite its success, however, the physical basis for this extinction law, or more correctly, attenuation law remains weak. Here we provide a physical explanation for the Calzetti Law based on a turbulent interstellar medium. In essence, this provides a log-normal distribution of column densities, giving a wide range of column densities in the dusty foreground screen. Therefore, extended sources such as starburst regions or HII regions seen through it suffer a point-to-point stochastic extinction and reddening. Regions of high column densities are "black" in the UV, but translucent in the IR, which leads to a flatter extinction law, and a larger value of the total to selective extinction, R_V. We fit the Calzetti Law, and infer that the variance sigma of the log-normal distribution lies in the range 0.6<sigma<2.2. The absolute to selective extinction R_V is found to be in the range 4.3 to 5.2 consistent with R_V=4.05+/-0.80 of the Calzetti Law.

We describe a new tool for the analysis of the UV to the sub-millimeter (sub-mm) spectral energy distribution (SED) of spiral galaxies. We use a consistent treatment of grain heating and emission, solve the radiation transfer problem for a finite disk and bulge, and self-consistently calculate the stochastic heating of grains placed in the resulting radiation field. We use this tool to analyse the well-studied nearby edge-on spiral galaxy NGC891. First we investigate whether the old stellar population in NGC891, along with a reasonable assumption about the young stellar population, can account for the heating of the dust and the observed far-infrared (FIR) and sub-mm emission. The dust distribution is taken from the model of Xilouris et al. (1999), who used only optical and near-IR observations to determine it. We have found that such a simple model cannot reproduce the SED of NGC891, especially in the sub-mm range. It underestimates by a factor of 2-4 the observed sub-mm flux. A number of possible explanations exist for the missing sub-mm flux. We investigate a few of them and demonstrate that one can reproduce the observed SED in the FIR and the sub-mm quite well, as well as the observed radial profile at 850 microns. For the models calculated we give the relative proportion of the dust radiation powered by the old and young stellar populations as a function of FIR/sub-mm wavelength. In all models we find that the dust is predominantly heated by the young stellar population.

We make predictions for the diffuse far-infrared (FIR) emission from dust in the intracluster medium (ICM) of the Virgo cluster using detailed information on potential dust sources, grain heating and sputtering rates available for this cluster from recent optical and X-ray studies. We show that the infrared emission from dust injected by sources inside the core region of the Virgo cluster is a factor of 10 below the detection limit of currently available telescopes. The outer regions of dynamically young clusters like the Virgo cluster have a further potential source of intracluster grains since they are still accreting freshly infalling spiral galaxies which are presumably contained in an accreting intergalactic medium (IGM). We show that cosmic ray driven winds from the infalling spirals can inject grains into a subvirial IGM that is external to the observed X-ray-emitting ICM. Sputtering during the injection process and in the IGM is weak, so that the injected grains should accumulate in the IGM until the infall brings them into contact with the hot ICM. This effect dominates the dust injection rate from known sources embedded in the hot Virgo ICM. Thus, any detection of diffuse IR emission would probe the current dust accretion rate for the cluster, acting as an indicator of the youth and the dynamical state of the cluster. The predictions for the Virgo cluster are generalised to other clusters and the possibility of detection of dynamically young clusters at cosmological distances is discussed.