We present the first radiative transfer (RT) model of a non-edge-on disk galaxy in which the large-scale geometry of stars and dust is self-consistently derived through fitting of multiwavelength imaging observations from the UV to the submm. To this end we used the axi-symmetric RT model of Popescu et al. and a new methodology for deriving geometrical parameters, and applied this to decode the{spectral energy distribution (SED) of M33. We successfully account for both the spatial and spectral energy distribution, with residuals typically within $7\%$ in the profiles of surface brightness and within $8\%$ in the spatially-integrated SED. We predict well the energy balance between absorption and re-emission by dust, with no need to invoke modified grain properties, and we find no submm emission that is in excess of our model predictions. We calculate that $80\pm8\%$ of the dust heating is powered by the young stellar populations. We identify several morphological components in M33, a nuclear, an inner, a main and an outer disc, showing a monotonic trend in decreasing star-formation surface-density ($\Sigma_{\rm SFR}$) from the nuclear to the outer disc. In relation to surface density of stellar mass, the $\Sigma_{\rm SFR}$ of these components define a steeper relation than the "main sequence" of star-forming galaxies, which we call a "structurally resolved main sequence". Either environmental or stellar feedback mechanisms could explain the slope of the newly defined sequence. We find the star-formation rate to be ${\rm SFR}=0.28^{+0.02}_{-0.01}{\rm M}_{\odot}{\rm yr}^{-1}$.

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 present a solution for the ultraviolet (UV) - submillimeter (submm) interstellar radiation fields (ISRFs) of the Milky Way, derived from modelling COBE, IRAS and Planck maps of the all-sky emission in the near-, mid-, far-infrared and submm.The analysis uses the axisymmetric radiative transfer (RT) model that we have previously implemented to model the panchromatic spectral energy distributions (SEDs) of star forming galaxies in the nearby universe, but with a new methodology allowing for optimisation of the radial and vertical geometry of stellar emissivity and dust opacity, as deduced from the highly resolved emission seen from the vantage point of the Sun. As such, this is the first self-consistent model of the broad-band continuum emission from the Milky Way. In this paper, we present model predictions for the spatially integrated SED of the Milky Way as seen from the Sun, showing good agreement with the data, and give a detailed description of the solutions for the distribution of ISRFs, as well as their physical origin, throughout the volume of the galaxy. We explore how the spatial and spectral distribution of our new predictions for the ISRF in the Milky Way affects the amplitude and spectral distribution of the gamma-rays produced via Inverse Compton scattering for cosmic ray electrons situated at different positions in the galaxy, as well as the attenuation of the gamma-rays due to interactions of the gamma-ray photons with photons of the ISRF. We also compare and contrast our solutions for the ISRF with those incorporated in the GALPROP package used for modelling the high energy emission from cosmic rays in the Milky Way.

The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.

A peculiar infrared ring-like structure was discovered by {\em Spitzer} around the strongly magnetised neutron star SGR 1900$+$14. This infrared structure was suggested to be due to a dust-free cavity, produced by the SGR Giant Flare occurred in 1998, and kept illuminated by surrounding stars. Using a 3D dust radiative transfer code, we aimed at reproducing the emission morphology and the integrated emission flux of this structure assuming different spatial distributions and densities for the dust, and different positions for the illuminating stars. We found that a dust-free ellipsoidal cavity can reproduce the shape, flux, and spectrum of the ring-like infrared emission, provided that the illuminating stars are inside the cavity and that the interstellar medium has high gas density ($n_H\sim$1000 cm$^{-3}$). We further constrain the emitting region to have a sharp inner boundary and to be significantly extended in the radial direction, possibly even just a cavity in a smooth molecular cloud. We discuss possible scenarios for the formation of the dustless cavity and the particular geometry that allows it to be IR-bright.

We present a new theoretical population synthesis model (the Galaxy Model) to examine and deal with large amounts of data from surveys of the Milky Way and to decipher the present and past structure and history of our own Galaxy. We assume the Galaxy to consist of a superposition of many composite stellar populations belonging to the thin and thick disks, the stellar halo and the bulge, and to be surrounded by a single dark matter halo component. A global model for the Milky Way's gravitational potential is built up self-consistently with the density profiles from the Poisson equation. In turn, these density profiles are used to generate synthetic probability distribution functions (PDFs) for the distribution of stars in colour-magnitude diagrams (CMDs). Finally, the gravitational potential is used to constrain the stellar kinematics by means of the moment method on a (perturbed)-distribution function. Spiral arms perturb the axisymmetric disk distribution functions in the linear response framework of density-wave theory where we present an analytical formula of the so-called `reduction factor' using Hypergeometric functions. Finally, we consider an analytical non-axisymmetric model of extinction and an algorithm based on the concept of probability distribution function to handle colour magnitude diagrams with a large number of stars. A genetic algorithm is presented to investigate both the photometric and kinematic parameter space. This galaxy model represents the natural framework to reconstruct the structure of the Milky Way from the heterogeneous data set of surveys such as Gaia-ESO, SEGUE, APOGEE2, RAVE and the Gaia mission.

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.

A large fraction of Galactic very-high energy (VHE; E$\gtrsim$100 GeV) $\gamma$-ray sources is cataloged as unidentified. In this work we explore the possibility that these unidentified sources are located in ambients particularly rich in material content unaccounted by traditional tracers. In a scenario where the VHE emission is due to the interaction of the accelerated particles with a target mass, a large mass of untraced material could be substantially contributing to the VHE emission from these regions. Here, we use three tracers for the commonly explored components: intensity of the $^\textrm{12}$CO(1$\rightarrow$0) line to trace the molecular material, HI hyperfine transition at 21cm to trace atomic hydrogen, and dust emission to trace the total hydrogen content. We show that the estimates of material content from these three tracers are compatible if the uncertainty on the respective conversion factors is taken into account. No additional gas component is found in these regions. However, a simple mass estimation from the $^\textrm{12}$CO(1$\rightarrow$0) line intensity might underestimate the total mass component in some locations.

We present a multiwavelength study of a sample of far-infrared (FIR) sources detected on the Herschel broad--band maps of the nearby galaxy M33. We perform source photometry on the FIR maps as well as mid-infrared (MIR), H$\alpha$, far-ultraviolet and integrated HI and CO line emission maps. By fitting MIR/FIR dust emission spectra, the source dust masses, temperatures and luminosities are inferred. The sources are classified based on their H$\alpha$ morphology (substructured versus not-substructured) and on whether they have a significant CO detection ($S/N>$3$\sigma$). We find that the sources have dust masses in the range 10$^2$-10$^4$~M$_\odot$ and that they present significant differences in their inferred dust/star formation/gas parameters depending on their H$\alpha$ morphology and CO detection classification. The results suggests differences in the evolutionary states or in the number of embedded HII regions between the subsamples. The source background--subtracted dust emission seems to be predominantly powered by local star formation, as indicated by a strong correlation between the dust luminosity and the dust-corrected H$\alpha$ luminosity and the fact that the extrapolated young stellar luminosity is high enough to account for the observed dust emission. Finally, we do not find a strong correlation between the dust-corrected H$\alpha$ luminosity and the dust mass of the sources, consistent with previous results on the breakdown of simple scaling relations at sub-kpc scales. However, the scatter in the relation is significantly reduced by correcting the H$\alpha$ luminosity for the age of the young stellar populations in the star--forming regions.

We present DART-Ray, a new ray-tracing 3D dust radiative transfer (RT) code designed specifically to calculate radiation field energy density (RFED) distributions within dusty galaxy models with arbitrary geometries. In this paper we introduce the basic algorithm implemented in DART-Ray which is based on a pre-calculation of a lower limit for the RFED distribution. This pre-calculation allows us to estimate the extent of regions around the radiation sources within which these sources contribute significantly to the RFED. In this way, ray-tracing calculations can be restricted to take place only within these regions, thus substantially reducing the computational time compared to a complete ray-tracing RT calculation. Anisotropic scattering is included in the code and handled in a similar fashion. Furthermore, the code utilizes a Cartesian adaptive spatial grid and an iterative method has been implemented to optimize the angular densities of the rays originated from each emitting cell. In order to verify the accuracy of the RT calculations performed by DART-Ray, we present results of comparisons with solutions obtained using the DUSTY 1D RT code for a dust shell illuminated by a central point source and existing 2D RT calculations of disc galaxies with diffusely distributed stellar emission and dust opacity. Finally, we show the application of the code on a spiral galaxy model with logarithmic spiral arms in order to measure the effect of the spiral pattern on the attenuation and RFED.

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 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 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 present a theoretical investigation of the light, radius and radial velocity variations of the prototype $\delta$ Cephei. We find that the best fit model accounts for luminosity and velocity amplitudes with an accuracy better than $0.8\sigma$, and for the radius amplitude with an accuracy of $1.7\sigma$. The chemical composition of this model suggests a decrease in both helium (0.26 vs 0.28) and metal (0.01 vs 0.02) content in the solar neighborhood. Moreover, distance determinations based on the fit of light curves agree at the $0.8\sigma$ level with the trigonometric parallax measured by the Hubble Space Telescope (HST). On the other hand, distance determinations based on angular diameter variations, that are independent of interstellar extinction and of the $p$-factor value, indicate an increase of the order of 5% in the HST parallax.