Using a representative sample of 14 star-forming dwarf galaxies in the local Universe, we show the existence of a spaxel-to-spaxel anti-correlation between the index N2 (log([NII]6583/Halpha)) and the Halpha flux. These two quantities are commonly employed as proxies for gas-phase metallicity and star formation rate (SFR), respectively. Thus, the observed N2 to Halpha relation may reflect the existence of an anti-correlation between the metallicity of the gas forming stars and the SFR it induces. Such an anti-correlation is to be expected if variable external metal-poor gas fuels the star-formation process. Alternatively, it can result from the contamination of the star-forming gas by stellar winds and SNe, provided that intense outflows drive most of the metals out of the star-forming regions. We also explore the possibility that the observed anti-correlation is due to variations in the physical conditions of the emitting gas, other than metallicity. Using alternative methods to compute metallicity, as well as previous observations of HII regions and photoionization models, we conclude that this possibility is unlikely. The radial gradient of metallicity characterizing disk galaxies does not produce the correlation either.

The optical integrated spectra of three LMC young stellar clusters (NGC 1984, NGC 1994 and NGC 2011) exhibit concave continua and prominent molecular bands which deviate significantly from the predictions of single stellar population (SSP) models. In order to understand the appearance of these spectra, we create a set of young stellar population (MILES) models, which we make available to the community. We use archival International Ultraviolet Explorer integrated UV spectra to independently constrain the cluster masses and extinction, and rule out strong stochastic effects in the optical spectra. In addition, we also analyze deep colour-magnitude diagrams of the clusters to provide independent age determinations based on isochrone fitting. We explore hypotheses including age-spreads in the clusters, a top-heavy initial mass function, different SSP models and the role of red supergiant stars (RSG). We find that the strong molecular features in the optical spectra can only be reproduced by modeling an increased fraction of about 20 per cent by luminosity of RSG above what is predicted by canonical stellar evolution models. Given the uncertainties in stellar evolution at Myr ages, we cannot presently rule-out the presence of Myr age-spreads in these clusters. Our work combines different wavelengths as well as different approaches (resolved data as well as integrated spectra for the same sample) in order to reveal the complete picture. We show that each approach provides important information but in combination can we better understand the cluster stellar populations.

Star Formation Rate (SFR) inferences are based in the so-called constant SFR approximation, where synthesis models are require to provide a calibration; we aims to study the key points of such approximation to produce accurate SFR inferences. We use the intrinsic algebra used in synthesis models, and we explore how SFR can be inferred from the integrated light without any assumption about the underling Star Formation history (SFH). We show that the constant SFR approximation is actually a simplified expression of more deeper characteristics of synthesis models: It is a characterization of the evolution of single stellar populations (SSPs), acting the SSPs as sensitivity curve over different measures of the SFH can be obtained. As results, we find that (1) the best age to calibrate SFR indices is the age of the observed system (i.e. about 13Gyr for z=0 systems); (2) constant SFR and steady-state luminosities are not requirements to calibrate the SFR; (3) it is not possible to define a SFR single time scale over which the recent SFH is averaged, and we suggest to use typical SFR indices (ionizing flux, UV fluxes) together with no typical ones (optical/IR fluxes) to correct the SFR from the contribution of the old component of the SFH, we show how to use galaxy colors to quote age ranges where the recent component of the SFH is stronger/softer than the older component. Particular values of SFR calibrations are (almost) not affect by this work, but the meaning of what is obtained by SFR inferences does. In our framework, results as the correlation of SFR time scales with galaxy colors, or the sensitivity of different SFR indices to sort and long scale variations in the SFH, fit naturally. In addition, the present framework provides a theoretical guide-line to optimize the available information from data/numerical experiments to improve the accuracy of SFR inferences.

In this document and the accompanying documents we describe a data model (Simulation Data Model) describing numerical computer simulations of astrophysical systems. The primary goal of this standard is to support discovery of simulations by describing those aspects of them that scientists might wish to query on, i.e. it is a model for meta-data describing simulations. This document does not propose a protocol for using this model. IVOA protocols are being developed and are supposed to use the model, either in its original form or in a form derived from the model proposed here, but more suited to the particular protocol. The SimDM has been developed in the IVOA Theory Interest Group with assistance of representatives of relevant working groups, in particular DM and Semantics.

Our understanding of stellar systems depends on the adopted interpretation of the IMF, phi(m). Unfortunately, there is not a common interpretation of the IMF, which leads to different methodologies and diverging analysis of observational data.We study the correlation between the most massive star that a cluster would host, mmax, and its total mass into stars, M, as an example where different views of the IMF lead to different results. We assume that the IMF is a probability distribution function and analyze the mmax-M correlation within this context. We also examine the meaning of the equation used to derive a theoretical M-char_mmax relationship, N x int[Char_mmax-mup] phi(m) dm = 1 with N the total number of stars in the system, according to different interpretations of the IMF. We find that only a probabilistic interpretation of the IMF, where stellar masses are identically independent distributed random variables, provides a self-consistent result. Neither M nor N, can be used as IMF scaling factors. In addition, Char_mmax is a characteristic maximum stellar mass in the cluster, but not the actual maximum stellar mass. A <M>-Char_mmax correlation is a natural result of a probabilistic interpretation of the IMF; however, the distribution of observational data in the N (or M)-cmmax plane includes a dependence on the distribution of the total number of stars, N (and M), in the system, Phi(N), which is not usually taken into consideration. We conclude that a random sampling IMF is not in contradiction to a possible mmax-M physical law. However, such a law cannot be obtained from IMF algebraic manipulation or included analytically in the IMF functional form. The possible physical information that would be obtained from the N (or M)-mmax correlation is closely linked with the Phi(M) and Phi(N) distributions; hence it depends on the star formation process and the assumed.

In a probabilistic framework of the interpretation of the initial mass function (IMF), the IMF cannot be arbitrarily normalized to the total mass, M, or number of stars, N, of the system. Hence, the inference of M and N when partial information about the studied system is available must be revised. (i.e., the contribution to the total quantity cannot be obtained by simple algebraic manipulations of the IMF). We study how to include constraints in the IMF to make inferences about different quantities characterizing stellar systems. It is expected that including any particular piece of information about a system would constrain the range of possible solutions. However, different pieces of information might be irrelevant depending on the quantity to be inferred. In this work we want to characterize the relevance of the priors in the possible inferences. Assuming that the IMF is a probability distribution function, we derive the sampling distributions of M and N of the system constrained to different types of information available. We show that the value of M that would be inferred must be described as a probability distribution Phi[M; m_a, N_a, Phi(N)] that depends on the completeness limit of the data, m_a, the number of stars observed down to this limit, N_a, and the prior hypothesis made on the distribution of the total number of stars in clusters, Phi(N).

Aims. In this paper we study the influence of the ionizing cluster mass on the emission line spectrum of Hii regions in order to determine the influence of low mass clusters on the integrated emission line spectra of galaxies. Methods. For this purpose, we present a grid of photoionization models that covers metallicities from Z = 0.001 to Z = 0.040, ages from 0.1 to 10 Ma (with a time step of 0.1 Ma), and cluster initial masses from 1 to 107 Mo. The stellar masses follow a Salpeter initial mass function (IMF) in an instantaneous burst mode of star formation. We obtain power-law scale-relations between emission-line luminosities and ionizing cluster masses from the grids and we evaluate the dependences on the ionizing cluster mass for some line luminosities, equivalent widths and line ratios. Results. Power-law scale-relations are shown to be useful tools to obtain robust diagnostics, as examples: (a) H?/H? ratio varies from the usually assumed value of 2.86, these variations imply the existence of a lower limit to the attainable precision in extinction estimations of ?E(B - V) ~ 0.1.; (b) EW(H?) is a good age indicator with a small dependence on cluster mass, while EW([O iii] 5007) shows a noteworthy mass dependence; (c) abundance estimations from R23 are practically unaffected by variations of the cluster mass; (d) estimations from S 23 and ?' would improve if the cluster mass dependences were considered and (e) [Oii] 3727/H? is a good star formation rate indicator for ages older than -4.5 Ma. We also show that the ionizing cluster mass dependence explains why empirical calibrations produce more reliable diagnostics of some emission lines than photoionization models grids. Finally, we show preliminary results about the contribution of low mass clusters (M < 104 Mo) to the integrated emission line spectra of galaxies, which can be as high as 80% for some relevant lines.

We present the underlying relations between colour-magnitude diagrams (CMDs) and synthesis models through the use of stellar luminosity distribution func- tions. CMDs studies make a direct use of the stellar luminosity distribution function while, in general, synthesis models only use its mean value, even though high-order moments can also be obtained. We show that the mean, high-order moments and in- tegrated luminosity distribution functions of stellar ensembles are related to the stellar luminosity distribution function, within the formalism of probabilistic synthesis mod- els. More details have been yet presented in Cervin ~ o & Luridiana (2006) and references therein. As a direct application of this formalism, we discuss two key issues. First, in- ferences on the upper mass limit of the initial mass function as a function of the total mass of clusters. Second, we apply extreme value theory to show that that the cluster mass obtained from normalising the IMF between mmax and mup does not provide the cluster mass in the case where only one star in this mass range is present, as assumed in the IGIMF theory. It provides instead the cluster mass with a 60% probability to have a star with mass larger than mmax, and we argue that in light of this result the basic formulation ofthe IGIMF theory must be revised.

Aims. We study how IMF sampling affects the ionizing flux and emission line spectra of low mass stellar clusters. Methods. We performed 2 x 10^6 Monte Carlo simulations of zero-age solar-metallicity stellar clusters covering the 20 - 10^6 Mo mass range. We study the distribution of cluster stellar masses, Mclus, ionizing fluxes, Q(H0), and effective temperatures, Tclus. We compute photoionization models that broadly describe the results of the simulations and compare them with photoionization grids. Results. Our main results are: (a) A large number of low mass clusters (80% for Mclus = 100 Mo) are unable to form an H ii region. (b) There are a few overluminous stellar clusters that form H ii regions. These overluminous clusters preserve statistically the mean value of <Q(H0)> obtained by synthesis models, but the mean value cannot be used as a description of particular clusters. (c) The ionizing continuum of clusters with Mclus < 10^4 Mo is more accurately described by an individual star with self-consistent effective temperature(T*) and Q(H0) than by the ensemble of stars (or a cluster Tclus) produced by synthesis models. (d)Photoionization grids of stellar clusters can not be used to derive the global properties of low mass clusters. Conclusions. Although variations in the upper mass limit, mup, of the IMF would reproduce the effects of IMF sampling, we find that an ad hoc law that relates mup to Mclus in the modelling of stellar clusters is useless, since: (a) it does not cover the whole range of possible cases, and (b) the modelling of stellar clusters with an IMF is motivated by the need to derive the global properties of the cluster: however, in clusters affected by sampling effects we have no access to global information of the cluster but only particular information about a few individual stars.

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

This paper presents the characterization of the optical range of the ALHAMBRA photometric system, a 20 contiguous, equal-width, medium-band CCD system with wavelength coverage from 3500A to 9700A. The photometric description of the system is done by presenting the full response curve as a product of the filters, CCD and atmospheric transmission curves, and using some first and second order moments of this response function. We also introduce the set of standard stars that defines the system, formed by 31 classic spectrophotometric standard stars which have been used in the calibration of other known photometric systems, and 288 stars, flux calibrated homogeneously, from the Next Generation Spectral Library (NGSL). Based on the NGSL, we determine the transformation equations between Sloan Digital Sky Survey (SDSS) ugriz photometry and the ALHAMBRA photometric system, in order to establish some relations between both systems. Finally we develop and discuss a strategy to calculate the photometric zero points of the different pointings in the ALHAMBRA project.

We analyzed the resolved stellar population of the C component of the extremely metal-poor dwarf galaxy Izw18 in order to evaluate its distance and star formation history as accurately as possible. In particular, we aimed at answering the question of whether this stellar population is young. We developed a probabilistic approach to analyzing high-quality photometric data obtained with the Advanced Camera for Surveys of the Hubble Space Telescope. This approach gives a detailed account of the various stochastic aspects of star formation. We carried out two successive models of the stellar population of interest, paying attention to how our assumptions could affect the results. We found a distance to the C component of I Zw 18 as high as 27 Mpc, a significantly higher value than those cited in previous works. The star formation history we inferred from the observational data shows various interesting features: a strong starburst that lasted for about 15 Myr, a more moderate one that occurred approx 100 Myr ago, a continuous process of star formation between both starbursts, and a possible episode of low level star formation at ages over 100 Myr. The stellar population studied is likely approx 125 Myr old, although ages of a few Gyr cannot be ruled out. Furthermore, nearly all the stars were formed in the last few hundreds of Myr.

We developed a new population synthesis code for groups of massive stars, where we model the emission of different forms of energy and matter from the stars of the association. In particular, the ejection of the two radioactive isotopes 26Al and 60Fe is followed, as well as the emission of hydrogen ionizing photons, and the kinetic energy of the stellar winds and supernova explosions. We investigate various alternative astrophysical inputs and the resulting output sensitivities, especially effects due to the inclusion of rotation in stellar models. As the aim of the code is the application to relatively small populations of massive stars, special care is taken to address their statistical properties. Our code incorporates both analytical statistical methods applicable to small populations, as well as extensive Monte Carlo simulations. We find that the inclusion of rotation in the stellar models has a large impact on the interactions between OB associations and their surrounding interstellar medium. The emission of 26Al in the stellar winds is strongly enhanced, compared to non-rotating models with the same mass-loss prescription. This compensates the recent reductions in the estimates of mass-loss rates of massive stars due to the effects of clumping. Despite the lower mass-loss rates, the power of the winds is actually enhanced for rotating stellar models. The supernova power (kinetic energy of their ejecta) is decreased due to longer lifetimes of rotating stars, and therefore the wind power dominates over supernova power for the first 6 Myr after a burst of star-formation. For populations typical of nearby star-forming regions, the statistical uncertainties are large and clearly non-Gaussian.

The ALHAMBRA survey aims to cover 4 square degrees using a system of 20 contiguous, equal width, medium-band filters spanning the range 3500 A to 9700 A plus the standard JHKs filters. Here we analyze deep near-IR number counts of one of our fields (ALH08) for which we have a relatively large area (0.5 square degrees) and faint photometry (J=22.4, H=21.3 and K=20.0 at the 50% of recovery efficiency for point-like sources). We find that the logarithmic gradient of the galaxy counts undergoes a distinct change to a flatter slope in each band: from 0.44 at [17.0, 18.5] to 0.34 at [19.5, 22.0] for the J band; for the H band 0.46 at [15.5, 18.0] to 0.36 at [19.0, 21.0], and in Ks the change is from 0.53 in the range [15.0, 17.0] to 0.33 in the interval [18.0, 20.0]. These observations together with faint optical counts are used to constrain models that include density and luminosity evolution of the local type-dependent luminosity functions. Our models imply a decline in the space density of evolved early-type galaxies with increasing redshift, such that only 30% - 50% of the bulk of the present day red-ellipticals was already in place at z~1.

In the next years, several cosmological surveys will rely on imaging data to estimate the redshift of galaxies, using traditional filter systems with 4-5 optical broad bands; narrower filters improve the spectral resolution, but strongly reduce the total system throughput. We explore how photometric redshift performance depends on the number of filters n_f, characterizing the survey depth through the fraction of galaxies with unambiguous redshift estimates. For a combination of total exposure time and telescope imaging area of 270 hrs m^2, 4-5 filter systems perform significantly worse, both in completeness depth and precision, than systems with n_f >= 8 filters. Our results suggest that for low n_f, the color-redshift degeneracies overwhelm the improvements in photometric depth, and that even at higher n_f, the effective photometric redshift depth decreases much more slowly with filter width than naively expected from the reduction in S/N. Adding near-IR observations improves the performance of low n_f systems, but still the system which maximizes the photometric redshift completeness is formed by 9 filters with logarithmically increasing bandwidth (constant resolution) and half-band overlap, reaching ~0.7 mag deeper, with 10% better redshift precision, than 4-5 filter systems. A system with 20 constant-width, non-overlapping filters reaches only ~0.1 mag shallower than 4-5 filter systems, but has a precision almost 3 times better, dz = 0.014(1+z) vs. dz = 0.042(1+z). We briefly discuss a practical implementation of such a photometric system: the ALHAMBRA survey.

(abridged) Non-ionizing stellar continua are a source of photons for continuum pumping in the hydrogen Lyman transitions. In the environments where these transitions are optically thick, deexcitation occurs through higher series lines, so that the flux in these lines has a fluorescent contribution in addition to recombination; in particular, Balmer emissivities are systematically enhanced above case B. The effectiveness of such mechanism in HII regions and the adequacy of photoionization models as a tool to study it are the two main focuses of this work. We find that photoionization models of H II regions illuminated by low-resolution population synthesis models significantly overpredict the fluorescent contribution to the Balmer lines. Conversely, photoionization models in which the non-ionizing part of the continuum is omitted or is not transferred underpredict the fluorescent contribution to the Balmer lines, producing a bias of similar amplitude in the opposite direction. In this paper, we carry out realistic estimations of the fluorescent Balmer intensity and discuss the variations to be expected as the simulated observational setup and the stellar population's parameters are varied. In all the cases explored, we find that fluorescent excitation provides a significant contribution. We also show that differential fluorescent enhancement may produce line-of-sight differences in the Balmer decrement, mimicking interstellar extinction. Fluorescent excitation emerges from our study as a small but important mechanism for the enhancement of Balmer lines, which should be taken into account in the abundance analysis of photoionized regions, particularly in the case of high-precision applications such as the determination of primordial helium.

This work aims to provide a theoretical formulation of Surface Brightness Fluctuations (SBF) in the framework of probabilistic synthesis models, and to distinguish between the different distributions involved in the SBF definition. RESULTS: We propose three definitions of SBF: (i) stellar population SBF, which can be computed from synthesis models and provide an intrinsic metric of fit for stellar population studies; (ii) theoretical SBF, which include the stellar population SBF plus an additional term that takes into account the distribution of the number of stars per resolution element psi(N); theoretical SBF coincide with Tonry & Schneider (1998) definition in the very particular case that psi(N) is assumed to be a Poisson distribution. However, the Poisson contribution to theoretical SBF is around 0.1% of the contribution due to the stellar population SBF, so there is no justification to include any reference to Poisson statistics in the SBF definition; (iii) observational SBF, which are those obtained in observations that are distributed around the theoretical SBF. Finally, we show alternative ways to compute SBF and extend the application of stellar population SBF to defining a metric of fitting for standard stellar population studies. CONCLUSIONS: We demostrate that SBF are observational evidence of a probabilistic paradigm in population synthesis, where integrated luminosities have an intrinsic distributed nature, and they rule out the commonly assumed deterministic paradigm of stellar population modeling.

We analyse high-resolution Hubble Space Telescope/Advanced Camera for Surveys imaging of the nuclear starburst region of M82, obtained as part of the Hubble Heritage mosaic made of this galaxy, in four filters (Johnson-Cousins equivalent B, V, and I broad bands, and an Halpha narrow-band filter), as well as subsequently acquired U-band images. We find a complex system of ~150 star clusters in the inner few 100 pc of the galaxy. We do not find any conclusive evidence of a cluster-formation epoch associated with the most recent starburst event, believed to have occurred about 4-6 Myr ago. This apparent evidence of decoupling between cluster and field-star formation is consistent with the view that star cluster formation requires special conditions. However, we strongly caution, and provide compelling evidence, that the `standard' simple stellar population analysis method we have used significantly underestimates the true uncertainties in the derived ages due to stochasticity in the stellar initial mass function and the corresponding sampling effects.

The fundamental properties of stellar clusters, such as the age or the total initial mass in stars, are often inferred from population synthesis models. The predicted properties are then used to constrain the physical mechanisms involved in the formation of such clusters in a variety of environments. Population synthesis models cannot, however, be applied blindy to such systems. We show that synthesis models cannot be used in the usual straightforward way to small-mass clusters (say, M < few times 10**4 Mo). The reason is that the basic hypothesis underlying population synthesis (a fixed proportionality between the number of stars in the different evolutionary phases) is not fulfilled in these clusters due to their small number of stars. This incomplete sampling of the stellar mass function results in a non-gaussian distribution of the mass-luminosity ratio for clusters that share the same evolutionary conditions (age, metallicity and initial stellar mass distribution function). We review some tests that can be carried out a priori to check whether a given cluster can be analysed with the fully-sampled standard population synthesis models, or, on the contrary, a probabilistic framework must be used. This leads to a re-assessment in the estimation of the low-mass tail in the distribution function of initial masses of stellar clusters.

We present here the VO access to the results of an analysis of the spectra of Sloan Digital Sky Survey (SDSS) galaxies performed with the STARLIGHT code by Cid Fernandes et al. (2005). The results include for each galaxy the original SDSS spectrum, the best-fit synthetic spectrum, the star formation history, the pure emission line spectrum corrected from underlying stellar population (in SDSS emission line galaxies) and the intensity of several emission/absorption lines. The database will be accessible from the PGos3 at the end of summer 2007.

Although full interoperativity between theoretical and observational data in the framework of the Virtual Observatory would be a very desirable achievement, the current status of VO offers few approaches to handle theoretical models. TSAP (Theoretical Spectra Access Protocol) has been proposed as a tool to fill this void, providing a simple scheme to easily operate with th is kind of data. TSAP is useful not only for synthetic spectra but also for other types of theoretical data. As an example we show an Isochrone and Evolutionary Tracks server using TSAP. Finally, we pay special attention to the correct treatment of the credits an important issue in the field of theoretical models.

In this contribution we investigate how to describe the results and usage of evolutionary synthesis models. In particular, we look for an explicit and quantitative description of the parameter space of synthesis models and the evaluation of their associated uncertainties and dispersion. First, we need to understand what synthesis models actually compute: we show that a synthetic stellar population with fixed physical parameters (age, metallicity, star formation history, initial mass function and size of the system) can only be described in terms of probability distributions (i.e. there is an intrinsic dispersion in any model). Second, we need to identify and characterize the coverage in the parameter space of the models (i.e. the combinations of input parameters that yield meaningful models) and the different sources of systematic errors. Third, we need a way to describe quantitatively the intrinsic dispersion, the systematic error and the parameter space coverage of the models. Up to now, the parameter space coverage and uncertainties have been described qualitatively in the models' reference papers, with potential misinterpretations by models' users. We show how Virtual Observatory developments enable a correct use of synthesis models to obtain accurate (and not simply precise) results.

In general, synthesis models provide the mean value of the distribution of possible integrated luminosities, this distribution (and not only its mean value) being the actual description of the integrated luminosity. Therefore, to obtain the closest model to an observation only provides confi- dence about the precision of such a fit, but not information about the accuracy of the result. In this contribution we show how to overcome this drawback and we propose the use of the theoretical mean-averaged dispersion that can be produced by synthesis models as a metric of fitting to infer accurate physical parameters of observed systems.

In this poster we present the analysis of the CMD of M67 (proposed in the Stellar Population Challenge) performed with VO applications. We found that, although the VO environment is still not ready to perform a complete analysis, its use provides highly useful additional information for the analysis. Thanks to the current VO framework, we are able to identify stars in the provided CMD that are not suitable for isochrone fitting. Additionally, we can complete our knowledge of this cluster extending the analysis to IR colors, which were not provided in the original data but that are available thanks to the VO. On the negative side, we find it difficult to access theoretical data from VO applications, so, currently, it is not possible to perform completely the analysis of the cluster inside the VO framework. However it is expected that the situation will improve in a near future.

The theory interest group in the International Virtual Observatory Alliance (IVOA) has the goal of ensuring that theoretical data and services are taken into account in the IVOA standards process. In this poster we present some of the efforts carried out by this group to include evolutionary synthesis models in the VO framework. In particular we present the VO tool PGos3, developed by the INAOE (Mexico) and the Spanish Virtual Observatory which includes most of public SSP models in the VO framework (e.g. VOSpec). We also describe the problems related with the inclusion of synthesis models in the VO framework and we try to encourage people to define the way in which synthesis models should be described. This issue has implications not only for the inclusion of synthesis models in the the VO framework but also for a proper usage of synthesis models.

Evolutionary synthesis models (ESM) have been extensively used to obtain the star formation history in galaxies by means of SED fitting. Implicit in this use of ESM is that (a) for given evolutionary parameters, the shape of the SED is fixed whatever the size of the observed cluster (b) all regions of the observed SED have the same weight in the fit. However, Nature does not follow these two assumptions, as is implied by the existence of Surface Brightness Fluctuations in galaxies and as can be shown by simple logical arguments.

The combination of stellar population synthesis models and photoionization models allows a better understanding of the spectral features of HII regions and HII galaxies. In this work we show that sampling effects in the initial mass function (IMF) are very important in the low cluster mass case. To this aim, we compute photoionization models ionized by realistic clusters made up of various combinations of individual stars and clusters made up with a synthesis model. We discuss the differences in the position on diagnostic diagrams and their implications.

The surface brightness fluctuations (SBF) method is a statistical method applied on image pixels in different bands. This contribution aims to distinguish between the observational (statistical) method and the theoretical (probabilistic) method based on stellar population synthesis and needed for the calibration of observational SBF. We find that the commonly used SBF theoretical definition as the "mean luminosity-weighted luminosity of the stellar population" is only compatible with the observational method under quite strong hypotheses, and that it is not compatible with stellar population theory results.

For a number of reasons, the properties of integrated stellar populations are distributed. Traditional synthesis models usually return the mean value of such distribution, and a perfect fitting to observational data is sought for to infer the age and metallicity of observed stellar populations. We show here that, while this is correct on average, it is not in individual cases because the mean may not be representative of actual values. We present a simple mathematical formalism to derive the shape of the population's luminosity distribution function (pLDF), and an abridged way to estimate it without computing it explicitly. This abridged treatment can be used to establish whether, for a specific case, the pLDF is Gaussian and the application of Gaussian tools, such as the chi^2 test, is correct. More in general, our formalism permits to compute weights to be attributed to different properties (spectral features or band luminosities) in the fitting process. We emphasize that our formalism does not supersede the results of traditionaly synthesis models, but permits to reinterpret and extend them into more powerful tools. The reader is referred to the original paper for further details.

Synthesis models predict the integrated properties of stellar populations. Several problems exist in this field, mostly related to the fact that integrated properties are distributed. To date, this aspect has been either ignored (as in standard synthesis models, which are inherently deterministic) or interpreted phenomenologically (as in Monte Carlo simulations, which describe distributed properties rather than explain them). We approach population synthesis as a problem in probability theory, in which stellar luminosities are random variables extracted from the stellar luminosity distribution function (sLDF). We derive the population LDF (pLDF) for clusters of any size from the sLDF, obtaining the scale relations that link the sLDF to the pLDF. We recover the predictions of standard synthesis models, which are shown to compute the mean of the sLDF. We provide diagnostic diagrams and a simplified recipe for testing the statistical richness of observed clusters, thereby assessing whether standard synthesis models can be safely used or a statistical treatment is mandatory. We also recover the predictions of Monte Carlo simulations, with the additional bonus of being able to interpret them in mathematical and physical terms. We give examples of problems that can be addressed through our probabilistic formalism. Though still under development, ours is a powerful approach to population synthesis. In an era of resolved observations and pipelined analyses of large surveys, this paper is offered as a signpost in the field of stellar populations.

In this review we address the uncertainties implicit in evolutionary synthesis model computations. After describing the general structure of synthesis codes, we discuss several source of uncertainties that may affect their results. In particular, we discuss the uncertainties arising in the computation of isochrones from evolutionary tracks; those related to atmosphere models; those that are a consequence of the incompleteness of the input ingredients; and those associated with the computational aspect used in synthesis codes. We also discuss the issue of the inclusion of distributed properties in synthesis models; as a paradigm of this case, we illustrate the difficulties implied by the inclusion of tracks with rotation in synthesis models. Finally, we describe several examples of the statistical approach to population synthesis. We report on the failure of the fuel consumption theorem (FCT) and the isochrone synthesis code to produce mutually consistent results. However, we argue that FCT and isochrone synthesis results are reliable for application to real systems in the wavelength range where they coincide. On the constructive side, we derive several useful survival strategies to bypass uncertainties. We show that single stellar populations at the turn-off ages of the tabulated tracks can be safely compared, as they are scarcely affected by the interpolation scheme used to compute isochrones. Finally, we suggest to use derivative quantities, such as the SN-rate, as bug detectors. On the recommendation side, we advocate for greater transparency and more documentation in synthesis modeling. We also ask stellar model makers to think of us and include more mass values in the tracks.

We present a library of 1654 high-resolution stellar spectra, with a sampling of 0.3 A and covering the wavelength range from 3000 A ~ to 7000 A The library was computed with the latest improvements in stellar atmospheres, incorporating non-LTE line-blanketed models for hot, massive (Teff >= 27500 K) and line-blanketed models for cool (3000 K =< Teff =< 4500 K) stars. The total coverage of the grid is 3000 K =< Teff =< 55000 K and -0.5 =< log g =< 5.5, for four chemical abundance values: twice solar, solar, half solar and 1/10 solar. Evolutionary synthesis models using this library are presented in a companion paper. We tested the general behavior of the library by calculating and comparing equivalent widths of numerous H and HeI lines, and some of the commonly used metallic indices. We also compared the library with the empirical libraries STELIB and Indo-US. The full set of the synthetic stellar spectra is available from our websites (http://www.iaa.csic.es/~rosa and http://www.astro.iag.usp.br/~lucimara/library.html).

We present the single stellar population (SSP) synthesis results of our new synthetic stellar atmosphere models library with a spectral sampling of 0.3 A, covering the wavelength range from 3000 A to 7000 A for a wide range of metallicities (twice solar, solar, half and 1/10 solar). The stellar library is composed of 1650 spectra computed with the latest improvements in stellar atmospheres. In particular it incorporates non-LTE line-blanketed models for hot (Teff >= 27500 K) and LTE line-blanketed models (Phoenix) for cool (3000 K =< Teff =< 4500 K) stars. Because of the high spectral resolution of this library, evolutionary synthesis models can be used to predict the strength of numerous weak absorption lines, and the evolution of the profiles of the strongest lines over a wide range of ages. The SSP results have been calculated for ages 1 Myr to 17 Gyr using the stellar evolutionary provided by the Geneva and Padova-tracks groups. For young stellar populations, our results have a very detailed coverage of high-temperature stars with similar results for Padova and Geneva isochrones. For intermediate and old stellar populations, our results, once degraded to a lower resolution, are similar to the ones obtained by other groups apart from limitations imposed by the stellar evolutionary physics. The limitations and advantages of our library for the analysis of integrated populations are described. The full set of the stellar library and the evolutionary models are available for retrieval at the websites http://www.iaa.csic.es/~rosa and http://www.iaa.csic.es/~mcs/Sed@, or on request from the first two authors.

We present new yields of 26Al from Wolf--Rayet stellar winds based on rotating stellar models which account well for numerous observed properties of massive stars. We study the impacts on the yields of a change of initial mass, metallicity and initial rotation velocity. We also consider the effects of a change of mass loss rates during the Wolf--Rayet phase. We show that for surface rotation velocities during the core H--burning phase matching to the observed ones, the quantity of 26Al ejected by a star of a given initial mass and metallicity is roughly doubled when the effects of rotation are taken into account. The metallicity dependence of the yield is, on the other hand, very similar to that obtained from non--rotating models. We estimate that at least about 20% to 50% (e.g. 0.6 --1.4 Msun) of the live 26Al detected in the Milky--Way originates from Wolf--Rayet stellar winds. We show the importance of a good knowledge of the present metallicity gradient and star formation rate in our galaxy for modeling both the variation of the 26Al surface density with the galactocentric distance and the global contribution of the Wolf--Rayet stellar winds to the present galactic mass of 26Al.

We have characterized the physical properties (electron temperature, density, metallicity) of the ionized gas and the ionizing population (age, metallicity, presence of WR stars) in the Lynx arc at z=3.36. This is a low metallicity (~10% solar) HII galaxy which is undergoing a burst of star formation of <5 Myr age. One possible scenario that explains the emission line spectrum of the Lynx arc, the large strength of the nitrogen lines and the HeII emission is that the object has experienced a merger event that has triggered a burst of star formation. WR stars have formed that contribute to a fast enrichment of the ISM and its overabundance of nitrogen.

We present the results of a double analysis of the ionizing cluster in NGC 588, a giant HII region (GHR) in the outskirts of the nearby galaxy M33. For this purpose, we obtained ground based long-slit spectroscopy and combined it with archival ground based and space borne imaging and spectroscopy, in the wavelength range 1100-9800 A. A first modeling of the cluster was performed using integrated properties, such as the spectral energy distribution (SED), broad band colors, nebular emission Hbeta equivalent width, the main ultraviolet resonance lines, and the presence of Wolf-Rayet star features. By applying standard assumptions about the initial mass function (IMF), we were unable to fit satisfactorily these observational data. This contradictory result led us to carry out a second modeling, based on a resolved photometric analysis of individual stars in Hubble Space Telescope (HST) images, by means of finding the best fit isochrone in color-magnitude diagrams (CMD), and assigning a theoretical SED to each individual star. The overall SED of the cluster, obtained by integrating the individual stellar SEDs, is found to fit better the observed SED than the best solution found through the integrated first analysis, but at a significantly later stage of evolution of the cluster of 4.2 Myr, as obtained from the best fit to the CMD. A comparative analysis of both methods traces the different results to the effects of statistical fluctuations in the upper end of the IMF, which are significant in NGC 588, with a computed cluster mass of 5600 Msun, as predicted by Cervino et al (2002). We discuss the results in terms of the strong influence of the few most massive stars, six in the case of NGC 588, that dominate the overall SED and, in particular, the ionizing far ultraviolet range beyond the Lyman limit.

We present new $^{26}$Al stellar yields from rotating Wolf--Rayet stellar models which, at solar metallicity, well reproduce the observed properties of the Wolf-Rayet populations. These new yields are enhanced with respect to non--rotating models, even with respect to non--rotating models computed with enhanced mass loss rates. We briefly discuss some implications of the use of these new yields for estimating the global contribution of Wolf-Rayet stars to the quantity of $^{26}$Al now present in the Milky Way.

(abriged) In this paper we establish a Lowest Luminosity Limit (LLL) for the use of synthesis codes. The limit is defined by the following statement: The total luminosity of the cluster modeled must be larger than the individual contribution of any of the stars included in the assumed isochrones. This limit is independent of the assumptions on the IMF and almost independent on the star formation history. We have obtained the LLL for a wide range of ages (5 Myr to 20 Gyr) and metallicities (Z=0 to Z=0.019) from Girardi et al. (2002) isochrones. Using the LLL and the results of evolutionary synthesis models we have also obtained the minimal cluster mass, Mmin, for which the results of synthesis models may suffer from a severe bias in the computation of colors. We show that the results of synthesis models for clusters with mass equal to Mmin have a relative dispersion about or larger than 32% (i.e. a dispersion of 0.35 mag) in the corresponding photometric bands. This effect is more important for near infrared bands (except in the case where Asymptotic Giant Branch stars dominate the luminosity in optical bands). In particular, Mmin takes values between 10^4 and 10^5 Mo for the K band. From the observational point of view, we show that in surveys that reach luminosities near the LLL, the color distributions will be skewed toward the color with lower effective sources (toward the red in general), and that the skewness is a signature of the distribution of initial cluster masses in the survey. We also apply the LLL to Globular Clusters showing that sampling effects are relevant in these clusters and that they can explain (at least partially) the bias of the observed colors with respect to the predictions of synthesis models.

(abriged) Stellar clusters with the same general physical properties (e.g., total mass, age, and star-formation mode) may have very different stellar mass spectra due to the incomplete sampling of the underlying mass function; such differences are especially relevant in the high-mass tail due to the smaller absolute number of massive stars. The dispersion in the number of massive stars also produces a dispersion in the properties of the corresponding ionising spectra. In this paper, we lay the bases for the future analysis of this effect by evaluating the dispersion in the ionising fluxes of synthetic spectra. As an important consequence, we found that the intensities of synthetic fluxes at different ionisation edges are strongly correlated, a fact suggesting that no additional dispersion will result from the inclusion of sampling effects in the analysis of diagnostic diagrams; this is true for HII regions on all scales. Additionally, we find convincing suggestions that the He II lines are strongly affected by sampling, and so cannot be used to constrain the evolutionary status of stellar clusters. We also establish the range of applicability of synthesis models set by the Lowest Luminosity Limit for the ionising flux, that is the lowest limit in cluster mass for which synthesis models can be applied to predict ionising spectra. This limit marks the boundary between the situations in which the ionising flux is better modeled with a single star as opposed to a star cluster; this boundary depends on the metallicity and age, ranging from 10^3 to more than 10^6 Mo. As a consequence, synthesis models should not be used to try to account for the properties of clusters with smaller masses.

We present evidence of the co-existence of either an AGN or an ultraluminous X-ray source (ULX), together with a young super stellar cluster in the 3 central parsecs of NGC4303. The galaxy contains a low luminosity AGN and hosts a number of starburst regions in a circumnuclear spiral, as well as in the nucleus itself. A high spatial resolution Chandra image of this source reveals that the soft X-ray emission traces the ultraviolet nuclear spiral down to a core, which is unresolved both in soft and hard X-rays. The astrometry of the X-ray core coincides with the UV core within the Chandra positioning accuracy. The total X-ray luminosity of the core, 1.5*10^{39} erg/s, is similar to that from some LINERs or from the weakest Seyferts detected so far. The soft X-rays in both the core and the extended structure surrounding it can be well reproduced by evolutionary synthesis models (which include the emission expected from single stars, the hot diffuse gas, supernova remnants and binary systems), consistent with the properties of the young stellar clusters identified in the UV. The hard X-ray tail detected in the core spectrum, however, most likely requires the presence of an additional source. This additional source could either be a weak active nucleus black hole or an ultraluminous X-ray object. The implications of these results are discussed.

This paper describes a new determination of the primordial helium abundance (Y_P), based on the abundance analysis of five metal-poor extragalactic HII regions. For three regions of the sample (SBS 0335-052, I Zw 18, and H29) we present tailored photoionization models based on improved calculations with respect to previous models. In particular, we use the photoionization models to study quantitatively the effect of collisional excitation of Balmer lines on the determination of the helium abundance (Y) in the individual regions. This effect is twofold: first, the intensities of the Balmer lines are enhanced with respect to the pure recombination value, mimicking a higher hydrogen abundance; second, the observed reddening is larger than the true extinction, due to the differential effect of collisions on different Balmer lines. In addition to these effects, our analysis takes into account the following features of HII regions: (i) the temperature structure, (ii) the density structure, (iii) the presence of neutral helium, (iv) the collisional excitation of the HeI lines, (v) the underlying absorption of the HeI lines, and (vi) the optical thickness of the HeI lines. The object that shows the highest increase in Y after the inclusion of collisional effects in the analysis is SBS 0335-052, whose helium abundance has been revised by Delta Y = +0.0107. The revised Y values for the five objects in our sample yield an increase of +0.0035 in Y_P, giving Y_P = 0.2391 +/- 0.0020.

We determine the physical conditions in NGC 6543, obtaining two different estimates of the electron temperature (T_e) and one estimate of the electron density (N_e). The electron temperature is computed by means of the [O III] ratio lambda 5007/lambda 4363, and of the diagnostic diagram combining lambda 5007 to the [O III] infrared lines 52m and 88m. The continuum intensity measured on slit spectra is much higher than theoretically predicted under the simplest assumptions. After considering several possibilities, we suggest enhanced 2-photon emission as the most probable source of the additional continuum. While T_e and N_e derived from the diagnostic diagram agree with the most recent determination, T_e derived from lambda 5007/lambda 4363 is smaller than previously published values, probably due to the bias in lambda 4363 introduced by the uncertainty in the continuum. Our main conclusion is that it is not possible, with present-day data, to derive accurate temperature maps of photoionized nebulae.

In this contribution we evaluate the minimal cluster mass for which the effects of an incomplete sampling of the Initial Mass Function (IMF) cannot be neglected. This minimal cluster mass corresponds to the situation in which the integrated luminosity of the cluster modeled equals the luminosity of the most luminous individual star included in the model, and it takes values between 3x10^2 and 6x10^5 Mo depending on the age and the observed band. We show different examples for young (t < 10 Myr) and old (t up to 10 Gyr) stellar populations. We also make here a first release of the spectral energy distribution (SED) from 160 mic. (1.8x10^12 Hz) to 25 keV (6.0x10^18 Hz) for star-forming regions with ages between 0.1 and 10 Myr, and metallicities between Z=0.001 and Z=0.040. The SEDs are available at our WWW server and include the corresponding quantities for the evaluation of sampling effects. These SEDs can be directly used to obtain colors (including the X-rays ones), as input of photoionization codes, or for chi^2 fitting with observed data taking into account the intrinsic uncertainty of the models due to the IMF sampling.

We revise the physical conditions in the O++ zone of the planetary nebula NGC 6543, obtaining two different estimates of the electron temperature (Te) and one estimate of the electron density (Ne). The electron temperature is computed by means of two independent methods, the nebular-to-auroral ratio [O III] 5007/4363, and the diagnostic diagram that combines 5007 to the [O III] infrared lines 52m and 88m. The optical and infrared fluxes have been obtained from archive HST/WFPC2 images and ISO LWS spectra respectively, and the continuum intensity in the optical has been measured from narrow-slit spectra obtained with the Isaac Newton Telescope at La Palma. The measured continuum intensity is higher than predicted by recombination theory under the hypothesis that all the Ly-alpha photons either escape or are destroyed. This fact can be explained in terms of an enhancement of the 2-photon continuum due to Ly-alpha conversion, a process that depends strongly on the local structure of the nebula. Alternative possibilities, outside the framework of recombination theory, have also been considered: e.g., the optical tail in the X-ray emission of a very hot plasma, and dust scattering of stellar radiation, but these hypotheses are not supported by quantitative estimates. While the electron temperature and density derived from the diagnostic diagram agree with the most recent determination, the temperature derived from 5007/4363,, Te(O++)_opt, is somewhat smaller than previously published values. We discuss several technical issues that contribute to the overall uncertainty in our results, focusing on the instrumental effects that might bias the [O III] 4363 intensity. We also discuss the effects of the collisional de-excitation of the O++ ground terms on the relation between Te(O++)_opt and Te(O++)_IR.

In this contribution I show the physical limitations to the application of standard synthesis models from 0 to 10 Gyr. I also present the multiwavelength spectrum (from gamma-rays to radio) of young star forming regions based on evolutionary population synthesis models that takes into account the stochasticity of the star formation proces. I show how the correlation of different observables can allow to establish a better understanding of the star formation processes. Finally I compare the spectral energy distribution of AGNs and star forming regions and I show the possible connection of Star Forming Galaxies with Seyfert 2 AGNs.

Evolutionary synthesis models have been used to study the physical properties of unresolved populations in a wide range of scenarios. Unfortunately, their self-consistency are difficult to test and there are some theoretical open questions without an answer: (1) The change of the homology relations assumed in the computation of isochrones due to the effect of stellar winds (or rotation) and the discontinuities in the stellar evolution are not considered. (2) There is no a consensus about how the isochrones must be integrated. (3) The discreteness of the stellar populations (that produce an intrinsic statistical dispersion) usually are not taken into account, and model results are interpreted in a deterministic way instead a statistical one... The objective of this contribution is to present some inconsistencies in the computation and some cautions in the application of the results of such codes.

Sampling fluctuations in stellar populations give rise to dispersions in observables when a small number of sources contribute effectively to the observables. This is the case for a variety of linear functions of the spectral energy distribution (SED) in small stellar systems, such as galactic and extragalactic HII regions, dwarf galaxies or stellar clusters. In this paper we show that sampling fluctuations also introduce systematic biases and multi-modality in non-linear functions of the SED, such as luminosity ratios, magnitudes and colours. In some cases, the distribution functions of rational and logarithmic quantities are bimodal, hence complicating considerably the interpretation of these quantities in terms of age or evolutionary stages. These biases can be only assessed by Monte Carlo simulations. We find that biases are usually negligible when the effective number of stars, Neff, which contribute to a given observable is larger than 10. Bimodal distributions may appear when Neff is between 10 and 0.1. Predictions from any model of stellar population synthesis become extremely unreliable for small Neff values, providing an operational limit to the applicability of such models for the interpretation of integrated properties of stellar systems. In terms of stellar masses, assuming a Salpeter Initial Mass Function in the range 0.08 -- 120 Mo, Neff=10 corresponds to about 10**5 Mo (although the exact value depends on the age and the observable). This bias may account, at least in part, for claimed variations in the properties of the stellar initial mass function in small systems, and arises from the discrete nature of small stellar populations.

We examine the impact of discrete numbers of stars in stellar populations on the results of Chemical Evolution Models. We explore the resulting dispersion in the true yields and their possible relation with the dispersion in observational data based on a Simple Closed-Box model. In this framework we find that the dispersion is larger for the less evolved or low abundance regions. Thus, the age-metallicity relation may be a tracer of the Star Formation History of our Galaxy. This theoretical dispersion is especially high for the relative abundance log(N/O) in regions where the total number of stars created is still low. This may explain part of the scatter in the N/O ratio observed in star forming galaxies. We have also found a first order theoretical estimation for the goodness of a linear fit of the Helium abundance vs. 12 + log (O/H) with values of the regression coefficient between 0.9 and 0.7 (independent of sampling effects). We conclude that it is necessary to include these sampling effects in a more realistic Chemical Evolution Model in order that such a model reproduces, at the same time, the mean value and the dispersion of observed abundances.

In this work we investigate the evolution of the X-ray emission of a cluster of single young massive stars with different metallicities. We have considered the X-ray contribution originated by the diffuse gas heated by the mechanical energy released by stellar winds and SN explosions as well as the X-ray contribution from SN remnants. The resulting ionizing spectrum (i.e. lambda < 912 A) has been used to compute the expected intensity of the nebular He II 4686 A. The observational ratio He II/Hbeta could be reproduced by the models assuming that a fraction of the mechanical energy produced by the star-formation episode is reprocessed by interaction with the ISM as soft X-ray radiation, contributing to the He ionization. However, the discreteness of the stellar populations affects the ionizing flux and may be responsible for the observed dispersion of the ratio. We have finally used the synthesis models to estimate the contribution of circumnuclear star-forming regions to the multiwavelength energy distribution in Active Galactic Nuclei, finding that the UV to soft X-ray continuum in many Seyfert 2 galaxies seems to be dominated by star-formation processes.

In this contribution, I present the 90% confidence limits on the diagnostic diagrams of EW(WR bump) and L(WR bump)/L(Hb) ratio vs. EW(Hb) resulting from evolutionary synthesis models that include the statistical dispersion due to finite stellar populations in real star forming regions.

Gamma-ray and microwave observations of the Cygnus region reveal an intense signal of 1.809 MeV line emission, attributed to radioactive decay of 26Al, that is closely correlated with 53 GHz free-free emission, originating from the ionised interstellar medium. We modelled both emissions using a multi-wavelength evolutionary synthesis code for massive star associations that we applied to the known massive star populations in Cygnus. For all OB associations and young open clusters in the field, we determined the population age, distance, and richness as well as the uncertainties in all these quantities from published photometric and spectroscopic data. We propagate the population uncertainties in model uncertainties by means of a Bayesian method. The young globular cluster Cyg OB2 turns out to be the dominant 26Al nucleosynthesis and ionisation source in Cygnus. Our model reproduces the ionising luminosity of the Cygnus region very well, yet it underestimates 26Al production by about a factor of 2. We attribute this underestimation to shortcomings of current nucleosynthesis models, and suggest the inclusion of stellar rotation as possible mechanism to enhance 26Al production. We also modelled 60Fe nucleosynthesis in the Cygnus region, yet the small number of recent supernova events suggests only little 60Fe production. Consequently, a detection of the 1.137 MeV and 1.332 MeV decay lines of 60Fe from Cygnus by the upcoming INTEGRAL observatory is not expected.

In this contribution I present my current work in a new generation of evolutionary synthesis models that compute the multiwavelength energy distribution (from gamma-rays to radio) as well as the associated dispersion for young stellar systems. I will also show some statistical effects that may appear in the analysis of surveys, like bimodal or multi-modal distributions and bias when color indices computed by the codes are compared with observations. Such new generation of synthesis models may be useful for the analysis of the data expected from GTC.

The probabilistic nature of the IMF in stellar systems implies that clusters of the same mass and age do not present the same unique values of their observed parameters. Instead they follow a distribution. We address the study of such distributions in terms of their confidence limits that can be obtained by evolutionary synthesis models. These confidence limits can be understood as the inherent uncertainties of synthesis models. We will compare such confidence limits arising from the discreteness of the number of stars obtained with Monte Carlo simulations with the dispersion resulting from an analytical formalism. We give some examples of the effects on the kinetic energy, V--K, EW(H$\beta$) and multiwavelength continuum.

We compare the rate of kinetic energy injected by stellar winds into the extragalactic H II region NGC 2363 to the luminosity needed to feed the observed temperature fluctuations. The kinetic luminosity associated to the winds is estimated by means of two different evolutionary synthesis codes, one of which takes into account the statistical fluctuations expected in the Initial Mass Function. We find that, even in the most favorable conditions considered by our model, such luminosity is much smaller than the luminosity needed to account for the observed temperature fluctuations. The assumptions underlying our study are emphasized as possible sources of uncertainty affecting our results.

In terms of statistical fluctuations, stellar population synthesis models are only asymptotically correct in the limit of a large number of stars, where sampling errors become asymptotically small. When dealing with stellar clusters, starbursts, dwarf galaxies or stellar populations within pixels, sampling errors introduce a large dispersion in the predicted integrated properties of these populations. We present here an approximate but generic statistical formalism which allows a very good estimation of the uncertainties and confidence levels in any integrated property, bypassing extensive Monte Carlo simulations, and including the effects of partial correlations between different observables. Tests of the formalism are presented and compared with proper estimates. We derive the minimum mass of stellar populations which is required to reach a given confidence limit for a given integrated property. As an example of this general formalism, which can be included in any synthesis code, we apply it to the case of young (t < 20 Myr) starburst populations. We show that, in general, the UV continuum is more reliable than other continuum bands for the comparison of models with observed data. We also show that clusters where more than 10^5 Mo have been transformed into stars have a relative dispersion of about 10% in Q(He+) for ages smaller than 3 Myr. During the WR phase the dispersion increases to about 25% for such massive clusters. We further find that the most reliable observable for the determination of the WR population is the ratio of the luminosity of the WR bump over the Hbeta luminosity. A fraction of the observed scatter in the integrated properties of clusters and starbursts can be accounted for by sampling fluctuations.

Ranging from track interpolation techniques through model atmospheres to the stochastic nature of the IMF, there are many uncertainties which need to be taken into account when modelling HR diagrams or performing population synthesis, particularly if comparison with actual data is sought. In this paper, we highlight (1) the problem of discontinuities along evolutionary tracks of massive stars (M> 8 Mo), showing that inconsistencies appear in the computation of the corresponding isochrones, and (2) the sampling fluctuations produced by the stochastic nature of the IMF, presenting a statistical formalism to estimate the dispersion in any given observable of a stellar population due to sampling effects which bypasses the need of performing Monte Carlo simulations.

The COMPTEL observations of the diffuse galactic 1.809 MeV emission attributed to the radioactive decay of 26Al have confirmed the diffuse nature of this interstellar emission line. One of the most significant features of the reconstructed intensity pattern is a flux enhancement in the direction of the Cygnus region. This region is fairly young and contains a wealth of massive stars, most of them grouped in the Cygnus OB associations. Multi-frequency model fitting strongly supports the hypothesis of massive stars and their descendent supernovae being the dominant sources of interstellar 26Al as observed by COMPTEL. Massive stars and supernovae are known to impart a large amount of kinetic energy into the surrounding ISM which lead to shockregions and large cavities. In addition, a large fraction of the electro-magnetic radiation of massive stars lies in the extreme ultra-violet causing photoionisation of the surrounding interstellar medium. We applied a population synthesis model in combination with an analytic model of the expansion of superbubbles to the Cygnus OB associations. The model predicts the expected 1.809 MeV flux and the gamma-ray line intensity due to interstellar 60Fe. We compute the sizes and expansion parameters of the expected HI-structures and the free-free emission intensities due to the photoionizing radiation from massive stars within this region of the sky. We discuss our present understanding of the Cygnus region with respect to the massive star census and history. Our model assigns about 70% of the 1.809 MeV intensity to six known OB associations, about 20% to known isolated sources and roughly 10% to an unkown diffuse component.

Evolutionary synthesis models are a fundamental tool to interpret the properties of observed stellar systems. In order to achieve a meaningful comparison between models and real data, it is necessary to calibrate the models themselves, i.e. to evaluate the dispersion due to the discreteness of star formation as well as the possible model errors. In this paper we show that linear interpolations in the log M - log t_k plane, that are customary in the evaluation of isochrones in evolutionary synthesis codes, produce unphysical results. We also show that some of the methods used in the calculation of time-integrated quantities (kinetic energy, and total ejected masses of different elements) may produce unrealistic results. We propose alternative solutions to solve both problems. Moreover, we have quantified the expected dispersion of these quantities due to stochastic effects in stellar populations. As a particular result, we show that the dispersion in the 14N/12C ratio increases with time.

We construct three extreme different scenarios of the star formation histories applicable to a sample of dwarf galaxies, based either on their present metallicity or their luminosity. The three possible scenarios imply different mechanical energy input rates and these we compare with the theoretical lower limits established for the ejection of processed matter out of dwarf galaxies. The comparison strongly points at the existence of extended gaseous haloes in these galaxies, acting as the barrier that allows galaxies to retain their metals and enhance their abundance. At the same time our findings strongly point at a continuous star-forming process, rather than to coeval bursts, as the main contributors to the overall metallicity in our galaxy sample.

The detection of 1.809 MeV gamma-ray line emission from the Cygnus X complex by the COMPTEL telescope is one of the most convincing proves of massive star nucleosynthesis in our Galaxy. The Cygnus X complex is an extremely active nearby region, containing several star forming regions, OB associations and young open star clusters. It houses some of the most massive stars known in our Galaxy and concentrates a large number of Wolf-Rayet stars. Thermal radio continuum emission and intense [C II] line emission reveals widespread ionisation, and at least 60 individual H II regions have been identified. In order to understand the 1.809 MeV gamma-ray line emission from the Cygnus X complex, and to compare the observations to theoretical nucleosynthesis calculations, we modelled the multi-wavelength spectrum of the region by means of an evolutionary synthesis model. Our investigation leads us to the following conclusions: 1. Stellar wind ejection is the dominant mechanism for the observed 26Al enrichment in the Cygnus X region. 2. Cyg OB2 is by far the dominating massive star association in Cygnus X and 1.809 MeV line emission from 26Al produced in this association should be detectable by the spectrometer SPI on INTEGRAL. 3. There is only low supernova activity in the Cygnus X complex and the 60Fe lines should be below INTEGRAL's detection sensitivity. The detectability of an individual massive star cluster (Cyg OB2) by INTEGRAL provides a fantastic opportunity for future nucleosynthesis studies using gamma-ray line spectroscopy. We will explore the scientific potential of such observations, and provide estimates for gamma-ray line intensity distributions based on the massive star census of the Cygnus X region.

When evolutionary synthesis models take into account the stochastic nature of the IMF together with the discrete number of stars in real stellar clusters, typical output turns to dispersion band (where real data can be placed) instead of narrow lines. We present here a qualitative analysis of such dispersion in the SN rate, the kinetic energy and the 14N/12C ratio for different amounts of mass transformed into stars.

We have developed a new diagnostic tool for the study of gamma-ray emission lines from radioactive isotopes (such as 26Al and 60Fe in conjunction with other multi-wavelength observables of Galactic clusters, associations, and alike objects. Our evolutionary synthesis models are based on the code of Cervino & Mas-Hesse (1994), which has been updated to include recent stellar evolution tracks, new stellar atmospheres for OB and WR stars, and nucleosynthetic yields from massive stars during hydrostatic burning phases and explosive SN II and SN Ib events. The temporal evolution of 26Al and 60Fe production, the equivalent yield of 26Al per ionising O7 V star (yal), and other observables are predicted for a coeval population. We also employed our model to estimate the steady state nucleosynthesis of a population of solar metallicity. To allow for a fully quantitative analysis of existing and future multi-wavelength observations, we propose a Bayesian approach that allows the inclusion of IMF richness effects and observational uncertainties in the analysis. In particular, a Monte Carlo technique is adopted to estimate probability distributions for all observables of interest. We outline the procedure of exploiting these distributions by applying our model to a fictive massive star association. Applications to existing observations of the Cygnus and Vela regions will be discussed in companion papers.

The stochastic nature of the IMF in young stellar clusters implies that clusters of the same mass and age do not present the same unique values of their observed parameters. Instead they follow a distribution. We address the study of such distributions, parameterised in terms of their confidence limits, in evolutionary synthesis models. These confidence limits can be understood as the inherent uncertainties of the synthesis models. Here we concentrate on some parameters such as EW(Hb) in emission. For instance, we show that for a cluster where 10^5 Mo have been transformed into stars, the dispersion of EW(Hb) is about 18% within the 90% confidence levels at ages between 3.5 and 5 Myrs.

Wolf-Rayet stars have been detected in a large number of galaxies experiencing intense bursts of star formation. All stars initially more massive than a certain, metallicity-dependent, value are believed to experience the Wolf-Rayet phase at the end of their evolution, just before collapsing in supernova explosion. The detection of Wolf-Rayet stars puts therefore important constraints on the evolutionary status of starbursts, the properties of their Initial Mass Functions and their star formation regime. In this contribution we review the properties of galaxies hosting Wolf-Rayet stars, with special emphasis on the factors that determine their presence and evolution, as well as their impact on the surrounding medium.

OB associations and young open clusters constitute the most prolific nucleosynthesis sites in our Galaxy. The combined activity of stellar winds and core-collapse supernovae ejects significant amounts of freshly synthesised nuclei into the interstellar medium. Radioactive isotopes, such as 26Al or 60Fe, that have been co-produced in such events may eventually be observed by gamma-ray instruments through their characteristic decay-line signatures. In order to study such signatures and to derive constraints on the involved nucleosynthesis processes, we developed a multi-wavelength evolutionary synthesis model for massive star associations. We apply this model to associations and clusters in the well-studied Cygnus region for which we re-determined the stellar census based on photometric and spectroscopic data. In particular we study the relation between 1.809 MeV gamma-ray line emission and ionising flux, since the latter has turned out to provide an excellent tracer of the global galactic 1.809 MeV emission. We compare our model to COMPTEL 1.8 MeV gamma-ray line observations from which we derive limits on the relative contributions from massive stars and core-collapse supernovae to the actual 26Al content in this region. Based on our model we make predictions about the expected 26Al and 60Fe line signatures in the Cygnus region. These predictions make the Cygnus region a prime target for the future INTEGRAL mission.

We present in this contribution our set of multiwavelength synthesis models including the evolution of single and binary stars. The main results we have obtained can be summarized as follows: (a) massive close-binary systems will start to experience mass transfer episodes after the first 4 Myr of the starburst evolution; (b) as a result of these mass transfer processes, stars of relatively low initial mass can loose completely their envelope and become a WR. In this way, the formation of WR stars is extended over longer than 15 Myr, and does not stop at 6 Myr as predicted by models including only single stars; (c) WR stars can thus be coeval with red supergiants, which peak at around 10 Myr for solar metallicities; (d) the accretion of mass will originate relatively massive stars at ages for which they should have already disappeared; these stars, together with the WR stars formed in rather evolved clusters, increase the production of ionizing photons, so that the H$\beta$ equivalent width will not drop as rapidly as predicted by models considering only individual stars; and (e) the mass transfer to compact companions will produce an additional source of high-energy radiation in the form of high-mass X-ray binaries, not predicted either by standard synthesis models.

We present in this contribution the predictions on the multiwavelength spectral energy distribution of our evolutionary population synthesis models including single and binary stellar systems. The high energy computations include the emission associated to X-ray binaries and supernovae remnants, as well as the mechanical energy released into the interstellar medium, which can be partially reprocessed into thermal X-rays. With these components we compute the spectral energy distribution of starburst galaxies from X-ray to radio ranges, and analyze finally the effects of the high energy emission on the H and He ionizing continuum.