Modelling galaxy spectra in presence of interstellar dust. I. The model of ISM and the library of dusty SSPs

In this paper, the first of a series of two devoted to modelling the spectra of galaxies of different morphological type in presence of dust, we present our description of the dust both in the diffuse ISM and the molecular clouds. Our model for the dust takes into account three components, i.e. graphite, silicates and polycyclic aromatic hydrocarbons. We consider and adapt to our aims two prescriptions for the size distribution of the dust grains and two models for the emission of the dusty ISM. We cross-check the emission and extinction models of the ISM by calculating the extinction curves and the emission for the typical environments of the Milky Way and the Large and Small Magellanic Clouds and by comparing the results with the observational data. The final model we have adopted is an hybrid one which stems from combining the analysis of Guhathakurta et al (1989) for the emission of graphite and silicates and Puget et al. (1985) for the PAH emission, and using the distribution law of Weingartner et al. (2001a) and the ionization model for PAHs of Weingartner et al. (2001b). We apply the model to calculate the SED of SSPs of different age and chemical composition, which may be severely affected by dust in young, massive stars while they are still embedded in their parental MCs. We use the "Ray Tracing" method to solve the problem of radiative transfer and to calculate libraries of SSP SEDs. Particular care is paid to model the contribution from PAHs, introducing different abundances of C in the population of very small carbonaceous grains (VSGs) and different ionization states in PAHs. The SEDs of young SSPs are then compared with observational data of star forming regions of four local galaxies successfully reproducing their SEDs from the UV-optical regions to the mid and far infrared.Comment: 24 pages, 29 figure

New response functions for absorption-line indices from high-resolution spectra

Basing on the huge library of 1-A resolution spectra calculated by Munari et al. over a large range of logT, log g, [Fe/H] and both for solar and a-enhanced abundance ratios [a/Fe], we present theoretical absorption-line indices on the Lick system. First we derive the so-called response functions (RFs) of Tripicco & Bell for a wide range of logT, log g, [Fe/H] and [a/Fe]=+0.4 dex. The RFs are commonly used to correct indices with solar [a/Fe] ratios to indices with [a/Fe]>0. Not only the RFs vary with the type of star but also with the metallicity. Secondly, with the aid of this and the fitting functions (FFs) of Worthey et al., we derive the indices for single stellar populations and compare them with those obtained by previous authors, e.g. Tantalo & Chiosi. The new RFs not only supersede the old ones by Tripicco & Bell, but also show that Hb increases with the degree of enhancement in agreement with the results by Tantalo & Chiosi. The new indices for single stellar populations are used to derive with aid of the recursive Minimum Distance method the age, metallicity and degree of enhancement of a sample of Galactic Globular Clusters for which these key parameters have been independently derived from the Colour-Magnitude Diagram and/or spectroscopic studies. The agreement is remarkably good.Comment: 18 pages, to appear in Astronomy & Astrophysic

Detailed AGB evolutionary models and near infrared colours of intermediate-age stellar populations: Tests on star clusters

We investigate the influence of Asymptotic Giant Branch stars on integrated colours of star clusters of ages between ~100 Myr and a few gigayears, and composition typical for the Magellanic Clouds. We use state-of-the-art stellar evolution models that cover the full thermal pulse phase, and take into account the influence of dusty envelopes on the emerging spectra. We present an alternative approach to the usual isochrone method, and compute integrated fluxes and colours using a Monte Carlo technique that enables us to take into account statistical fluctuations due to the typical small number of cluster stars. We demonstrate how the statistical variations in the number of Asymptotic Giant Branch stars and the temperature and luminosity variations during thermal pulses fundamentally limit the accuracy of the comparison (and calibration, for population synthesis models that require a calibration of the Asymptotic Giant Branch contribution to the total luminosity) with star cluster integrated photometries. When compared to observed integrated colours of individual and stacked clusters in the Magellanic Clouds, our predictions match well most of the observations, when statistical fluctuations are taken into account, although there are discrepancies in narrow age ranges with some (but not all) set of observations.Comment: 12 pages, 14 figures, accepted for publication in A&

The Role of Dust in Models of Population Synthesis

We have employed state-of-the-art evolutionary models of low and intermediate-mass AGB stars, and included the effect of circumstellar dust shells on the spectral energy distribution (SED) of AGB stars, to revise the Padua library of isochrones (Bertelli et al. 1994). The major revision involves the thermally pulsing AGB phase, that is now taken from fully evolutionary calculations by Weiss & Ferguson (2009). Two libraries of about 600 AGB dust-enshrouded SEDs each have also been calculated, one for oxygen-rich M-stars and one for carbon-rich C-stars. Each library accounts for different values of input parameters like the optical depth {\tau}, dust composition, and temperature of the inner boundary of the dust shell. These libraries of dusty AGB spectra have been implemented into a large composite library of theoretical stellar spectra, to cover all regions of the Hertzsprung-Russell Diagram (HRD) crossed by the isochrones. With the aid of the above isochrones and libraries of stellar SEDs, we have calculated the spectro-photometric properties (SEDs, magnitudes, and colours) of single-generation stellar populations (SSPs) for six metallicities, more than fifty ages (from 3 Myr to 15 Gyr), and nine choices of the Initial Mass Function. The new isochrones and SSPs have been compared to the colour-magnitude diagrams (CMDs) of field populations in the LMC and SMC, with particular emphasis on AGB stars, and the integrated colours of star clusters in the same galaxies, using data from the SAGE (Surveying the Agents of Galaxy Evolution) catalogues. We have also examined the integrated colours of a small sample of star clusters located in the outskirts of M31. The agreement between theory and observations is generally good. In particular, the new SSPs reproduce the red tails of the AGB star distribution in the CMDs of field stars in the Magellanic Clouds.Comment: Accepted for publication in MNRA

Study of ISM tracers in galaxies

We collected data for two samples of normal and interacting galaxies for a total of 2953 galaxies having fluxes in one or more of the following wavebands: FIR, 21 cm line, CO(1-0) lines and soft X-ray. The large set of data obtained allowed us to revisit some of the already known relations between the different tracers of the interstellar medium (ISM), such as the link between the FIR flux and the CO line emission, the relation between X-ray emission and the blue or FIR luminosity. The relation lacking from observations for early-type galaxies has been discussed and explained in detail in the frame of a suitable theoretical model, obtained by coupling chemo-dynamical N-body simulations with a dusty spectrophotometric code of population synthesis.Comment: 2 pages, o appear in the Proceedings of the Conf. "From Stars to Galaxies: Building the Pieces to Build Up the Universe", Vallenari et al. eds., ASP Conf. Serie