No Arabic abstract
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 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 have acquired radio continuum data between 70,MHz and 48,GHz for a sample of 19 southern starburst galaxies at moderate redshifts ($0.067 < z < 0.227$) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework we find the radio continuum is rarely characterised well by a single power law, instead often exhibiting low frequency turnovers below 500,MHz, steepening at mid-to-high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 to 500,MHz the radio-continuum at low frequency ($ u < 200$,MHz) could be overestimated by upwards of a factor of twelve if a simple power law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be $alpha=-1.06$, which is steeper then the canonical value of $-0.8$ for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution.
The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.
High-resolution non-ideal magnetohydrodynamical simulations of the turbulent magnetized ISM, powered by supernovae types Ia and II at Galactic rate, including self-gravity and non-equilibriuim ionization (NEI), taking into account the time evolution of the ionization structure of H, He, C, N, O, Ne, Mg, Si, S and Fe, were carried out. These runs cover a wide range (from kpc to sub-parsec) of scales, providing resolution independent information on the injection scale, extended self-similarity and the fractal dmension of the most dissipative structures.
We present new calculations of the attenuation of stellar light from spiral galaxies using geometries for stars and dust which can reproduce the entire spectral energy distribution from the UV to the FIR/submm and can also account for the surface brightness distribution in both the optical/NIR and FIR/submm. The calculations are based on the model of Popescu et al. (2000), which incorporates a dustless stellar bulge, a disk of old stars with associated diffuse dust, a thin disk of young stars with associated diffuse dust, and a clumpy dust component associated with star-forming regions in the thin disk. The attenuations, which incorporate the effects of multiple anisotropic scattering, are derived separately for each stellar component, and presented in the form of easily accessible polynomial fits as a function of inclination, for a grid in optical depth and wavelength. The wavelength range considered is between 912 AA and 2.2 micron, sampled such that attenuation can be conveniently calculated both for the standard optical bands and for the bands covered by GALEX. The attenuation characteristics of the individual stellar components show marked differences between each other. A general formula is given for the calculation of composite attenuation, valid for any combination of the bulge-to-disk ratio and amount of clumpiness. As an example, we show how the optical depth derived from the variation of attenuation with inclination depends on the bulge-to-disk ratio. Finally, a recipe is given for a self-consistent determination of the optical depth from the Halpha/Hbeta line ratio.