Modelling Intermediate Age and Old Stellar Populations in the Infrared

We have investigated the spectro-photometric properties of the Asymptotic Giant Branch (AGB) stars and their contribution to the integrated infrared emission in simple stellar populations (SSP). Adopting analytical relations describing the evolution of these stars in the HR diagram and empirical relations for the mass-loss rate and the wind terminal velocity, we were able to model the effects of the dusty envelope around these stars, with a minimal number of parameters. We computed isochrones at different age and initial metal content. We compare our models with existing infrared colors of M giants and Mira stars and with IRAS PSC data. Contrary to previous models, in the new isochrones the mass-loss rate, which establishes the duration of the AGB phase, also determines the spectral properties of the stars. The contribution of these stars to the integrated light of the population is thus obtained in a consistent way. We find that the emission in the mid infrared is about one order of magnitude larger when dust is taken into account in an intermediate age population, irrespective of the particular mixture adopted. The dependence of the integrated colors on the metallicity and age is discussed, with particular emphasis on the problem of age-metallicity degeneracy. We show that, contrary to the case of optical or near infrared colors, the adoption of a suitable pass-band in the mid infrared allows a fair separation of the two effects. We suggest intermediate redshift elliptical galaxies as possible targets of this method of solving the age-metallicity dilemma. The new SSP models constitute a first step in a more extended study aimed at modelling the spectral properties of the galaxies from the ultraviolet to the far infrared.Comment: 16 pages, 10 figures, to appear in A&

SPITZER IRS spectra of Virgo early type galaxies: detection of stellar silicate emission

We present high signal to noise ratio Spitzer Infrared Spectrograph observations of 17 Virgo early-type galaxies. The galaxies were selected from those that define the colour-magnitude relation of the cluster, with the aim of detecting the silicate emission of their dusty, mass-losing evolved stars. To flux calibrate these extended sources we have devised a new procedure that allows us to obtain the intrinsic spectral energy distribution and to disentangle resolved and unresolved emission within the same object. We have found that thirteen objects of the sample (76%) are passively evolving galaxies with a pronounced broad silicate feature which is spatially extended and likely of stellar origin, in agreement with model predictions. The other 4 objects (24%) are characterized by different levels of activity. In NGC 4486 (M 87) the line emission and the broad silicate emission are evidently unresolved and, given also the typical shape of the continuum, they likely originate in the nuclear torus. NGC 4636 shows emission lines superimposed on extended (i.e. stellar) silicate emission, thus pushing the percentage of galaxies with silicate emission to 82%. Finally, NGC 4550 and NGC 4435 are characterized by polycyclic aromatic hydrocarbon (PAH) and line emission, arising from a central unresolved region. A more detailed analysis of our sample, with updated models, will be presented in a forthcoming paper.Comment: 6 pages; ApJ Letters, accepte

Probing the Galaxy I. The galactic structure towards the galactic pole

Observations of (B-V) colour distributions towards the galactic poles are compared with those obtained from synthetic colour-magnitude diagrams to determine the major constituents in the disc and spheroid. The disc is described with four stellar sub-populations: the young, intermediate, old, and thick disc populations, which have respectively scale heights of 100 pc, 250 pc, 0.5 kpc, and 1.0 kpc. The spheroid is described with stellar contributions from the bulge and halo. The bulge is not well constrained with the data analyzed in this study. A non-flattened power-law describes the observed distributions at fainter magnitudes better than a deprojected R^{1/4}-law. Details about the age, metallicity, and normalizations are listed in Table 1. The star counts and the colour distributions from the stars in the intermediate fields towards the galactic anti-centre are well described with the stellar populations mentioned above. Arguments are given that the actual solar offset is about 15 pc north from the galactic plane.Comment: 11 pages TeX, 4 separate pages with additional figures, accepted for publication in A&

Star Formation History and Extinction in the central kpc of M82-like Starbursts

We report on the star formation histories and extinction in the central kpc region of a sample of starburst galaxies that have similar far infrared (FIR), 10 micron and K-band luminosities as those of the archetype starburst M82. Our study is based on new optical spectra and previously published K-band photometric data, both sampling the same area around the nucleus. Model starburst spectra were synthesized as a combination of stellar populations of distinct ages formed over the Hubble time, and were fitted to the observed optical spectra and K-band flux. The model is able to reproduce simultaneously the equivalent widths of emission and absorption lines, the continuum fluxes between 3500-7000 Ang, the K-band and the FIR flux. We require a minimum of 3 populations -- (1) a young population of age < 8 Myr, with its corresponding nebular emission, (2) an intermediate-age population (age < 500 Myr), and (3) an old population that forms part of the underlying disk or/and bulge population. The contribution of the old population to the K-band luminosity depends on the birthrate parameter and remains above 60% in the majority of the sample galaxies. Even in the blue band, the intermediate age and old populations contribute more than 40% of the total flux in all the cases. A relatively high contribution from the old stars to the K-band nuclear flux is also apparent from the strength of the 4000 Ang break and the CaII K line. The extinction of the old population is found to be around half of that of the young population. The contribution to the continuum from the relatively old stars has the effect of diluting the emission equivalent widths below the values expected for young bursts. The mean dilution factors are found to be 5 and 3 for the Halpha and Hbeta lines respectively.Comment: 20 pages, uses emulateapj.cls. Scheduled to appear in ApJ Jan 1, 200

The role of the synchrotron component in the mid infrared spectrum of M 87

We study in detail the mid-infrared Spitzer-IRS spectrum of M 87 in the range 5 to 20 micron. Thanks to the high sensitivity of our Spitzer-IRS spectra we can disentangle the stellar and nuclear components of this active galaxy. To this end we have properly subtracted from the M 87 spectrum, the contribution of the underlying stellar continuum, derived from passive Virgo galaxies in our sample. The residual is a clear power-law, without any additional thermal component, with a zero point consistent with that obtained by high spatial resolution, ground based observations. The residual is independent of the adopted passive template. This indicates that the 10 micron silicate emission shown in spectra of M 87 can be entirely accounted for by the underlying old stellar population, leaving little room for a possible torus contribution. The MIR power-law has a slope alpha ~ 0.77-0.82 (S$_\nu\propto\nu^{-\alpha}$), consistent with optically thin synchrotron emission.Comment: 5 pages, 4 figures, accepted for publication in ApJ main journa

Panchromatic models of galaxies: GRASIL

We present here a model for simulating the panchromatic spectral energy distribution of galaxies, which aims to be a complete tool to study the complex multi-wavelength picture of the universe. The model take into account all important components that concur to the SED of galaxies at wavelengths from X-rays to the radio. We review the modeling of each component and provide several applications, interpreting observations of galaxy of different types at all the wavelengths.Comment: 10 pages, 4 figures, invited talk, to appear in the proceedings of: "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu and R.J. Tuffs, AIP Conf. Ser., in pres

Source parameters and stress release of seismic sequences occurred in the Friuli Venezia Giulia region (Northeastern Italy) and in Western Slovenia

The source parameters of the major events of a swarm and of two seismic sequences, occurred in the Friuli area (Northeastern Italy) and in Western Slovenia, were estimated. The Claut swarm (C96) occurred since the end of January to June 1996, with a MD 4.3 major shock and it appears composed of three sub-sequences. The two sequences are the Kobarid sequence (K98) started on April 12, 1998 with a MD 5.6 mainshock and the M.te Sernio (S02) sequence caused by the February 14, 2002 earthquake (MD = 4.9). Acceleration and velocity data recorded by the local seismic network of the Istituto Nazionale di Oceanografia e di Geofisica Sperimenale (OGS) and corrected for attenuation, were employed to estimate seismic moments and radiated energies. Source dimensions were inferred from the computed corner frequencies and the stress release was estimated from the Brune stress drop, the apparent stress and the RMS stress drop. On the whole, seismic moments range from 1.7×1012 to 1.1×1017 Nm, and radiated energies are in the range 106–1013 J. Brune stress drops are scattered and do not show any evidence of a self-similarity breakdown for sources down to 130m radius. The radiated seismic energy scales as a function of seismic moment, with a slope of the scaling relation that decreases for increasing seismic moments. The mechanism of stress release was analyzed by computing the ε parameter of Zuniga [Zuniga, R., 1993. Frictional overshoot and partial stress drop. Which one? Bull. Seismol. Soc. Am. 83, 939–944]. The K98 and S02 sequences are characterized by a wide range of the ε parameter with stress drop mechanism varying from partial locking to overshoot cases. The ε values of the C96 swarm are more homogeneous and close to the Orowan’s condition. The radiated seismic energy and the ratio of stress drop between mainshock and aftershocks appear different among the analyzed cases.We therefore investigated the relationship between the stress parameters of the main shock and the energy radiated by the aftershock sequences. For this purpose, we also estimated the source parameters of two other sequences occurred in the area, with mainshocks of MD 4.1 and 5.1, respectively. We found a positive correlation between the Brune stress drop of the mainshock and the ratio between the radiated energy of the mainshock and the summation of the energies radiated by the aftershocks

Cepheid Mass-loss and the Pulsation -- Evolutionary Mass Discrepancy

I investigate the discrepancy between the evolution and pulsation masses for Cepheid variables. A number of recent works have proposed that non-canonical mass-loss can account for the mass discrepancy. This mass-loss would be such that a 5Mo star loses approximately 20% of its mass by arriving at the Cepheid instability strip; a 14Mo star, none. Such findings would pose a serious challenge to our understanding of mass-loss. I revisit these results in light of the Padova stellar evolutionary models and find evolutionary masses are ($17\pm5$)% greater than pulsation masses for Cepheids between 5<M/Mo<14. I find that mild internal mixing in the main-sequence progenitor of the Cepheid are able to account for this mass discrepancy.Comment: 15 pages, 3 figures, ApJ accepte