Search

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., Comment: A&A accepted. DART-Ray can be downloaded at https://github.com/gnatale/DART-Ray/; Code documentation: http://www.star.uclan.ac.uk/~gn/dartray_doc/; N-body/SPH galaxy simulation examples: http://www.star.uclan.ac.uk/mgsv/

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., Comment: 22 pages, 9 figures, accepted for publication in A&A

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., Comment: Accepted for publication in MNRAS. Appendix available only in the journal online version

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., Comment: 24 pages, 14 figures

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)., Comment: 12 pages, 6 figures, appendices will be published only online

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., Comment: paper accepted for publication, 14 pages, 7 figures

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., Comment: 18 pages, 12 figures accepted for publication (13/4/2012)

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/) 6., Comment: 6 pages, 6 figures accepted for publication in A&A, section 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., Comment: 24 pages, 20 figures

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), Comment: 39 pages; accepted for publication in Astronomy & Astrophysics; For a higher resolution version of Fig.1 and Fig.20 see http://www.star.uclan.ac.uk/~ccp/index.shtml

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., Comment: accepted for publication in ApJS, May 2008, v176n1 issue

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., Comment: 22 pages. 18 figures. Accepted for publication in Astrophysical journal Supplements. Electronic tabular material is available on request to Michael.Dopita@anu.edu.au

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., Comment: 36 pages, 8 figures. To be published in Astrophysical Journal, August 2006

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 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., Comment: 38 pages, 48 figures, accepted for publication in ApJ, Jan. 2005

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., Comment: 10 pages, 11 figures, conference proceedings paper, "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu & R.J. Tuffs, AIP Conf. Ser., in press

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., Comment: 12 pages, 16 figures

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), Comment: 56 pages, 16 figures, accepted by The Apstrophysical Journal For version with full, colour figures go to http://www.mso.anu.edu.au/~bgroves/starburst/

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.

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

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., Comment: 15 pages, accepted for publication in Astronomy and Astrophysics

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., Comment: 19 pages, accepted for publication in Astronomy and Astrophysics