We analyze the chemical composition of the planetary nebula (PN) NGC 5315, through high-resolution (R = 40000) optical spectroscopy with UVES at the Very Large Telescope, and medium-resolution (R = 4800) near-infrared spectroscopy with FIRE at Magellan Baade Telescope, covering a wide spectral range from 0.31 to 2.50 micron. The main aim of this work is to investigate neutron (n)-capture element abundances to study the operation of the slow n-capture ("s-process") in the AGB progenitor of NGC 5315. We detect more than 700 emission lines, including ions of the n-capture elements Se, Kr, Xe, and possibly Br. We compute physical conditions from a large number of diagnostic line ratios, and derive ionic abundances for species with available atomic data. The total abundances are computed using recent ionization correction factors (ICFs) or by summing ionic abundances. Total abundances of common elements are in good agreement with previous work on this object. Based on our abundance analysis of NGC 5315, including the lack of s-process enrichment, we speculate that the most probable scenario is that the progenitor star is in a binary system as hinted at by radial velocity studies, and interactions with its companion truncated the AGB before s-process enrichment could occur. However there are other two possible scenarios for its evolution, that cannot be ruled out: i) the progenitor is a low-mass single star that did not undergo third dredge-up; ii) the progenitor star of NGC 5315 had an initial mass of 4--6 solar masses, and any s-process enhancements were heavily diluted by the massive envelope during the AGB phase.

We identify [Se III] 1.0994 micron in the planetary nebula (PN) NGC 5315 and [Kr VI] 1.2330 micron in three PNe, from spectra obtained with the FIRE spectrometer on the 6.5-m Baade Telescope. Se and Kr are the two most widely-detected neutron-capture elements in astrophysical nebulae, and can be enriched by s-process nucleosynthesis in PN progenitor stars. The detection of [Se III] 1.0994 micron is particularly valuable when paired with observations of [Se IV] 2.2858 micron, as it can be used to improve the accuracy of nebular Se abundance determinations, and allows Se ionization correction factor (ICF) schemes to be empirically tested for the first time. We present new effective collision strength calculations for Se^{2+} and Kr^{5+}, which we use to compute ionic abundances. In NGC 5315, we find that the Se abundance computed from Se^{3+}/H^+ is lower than that determined with ICFs that incorporate Se^{2+}/H^+. We compute new Kr ICFs that take Kr^{5+}/H^+ into account, by fitting correlations found in grids of Cloudy models between Kr ionic fractions and those of more abundant elements, and use these to derive Kr abundances in four PNe. Observations of [Se III] and [Kr VI] in a larger sample of PNe, with a range of excitation levels, are needed to rigorously test the ICF prescriptions for Se and our new Kr ICFs.

We have done a spectroscopical analysis of the type I planetary nebula (PN) NGC 5315, through high-resolution (R$\sim$40000) optical spectroscopy with UVES at the 8.2m Very Large Telescope, and medium-resolution (R$\sim$4800) near-IR spectroscopy with FIRE at the 6.5m Magellan Baade telescope, covering a wide spectral range from 0.31 $\mu$m to 2.50 $\mu$m. The main aim of this work is to investigate the slow neutron(n)-capture process (the s-process) in the Asymptotic Giant Branch (AGB) star progenitor of a type I PNe. We detected and identified about 700 features, including lines from the n-capture elements Kr, Se, and possibly Br and Xe. We compute physical conditions using line ratios of common ions. Ionic abundances are computed for the species with available atomic data. We calculate total abundances using recent ionization correction factors (ICFs) or by summing ionic abundances. Our results for common elements are in good agreement with previous works on the same object. We do not find a substantial s-process enrichment in NGC 5315, which is typical for type I PNe.

Based on observations of HII regions and the new computations of the recombination coefficients of the He I lines by Porter et al. (2013) we obtain a primordial helium abundance by mass of $Y_P = 0.2446\pm0.0029$. We consider thirteen sources of error for the $Y_P$ determination, some of them are mainly due to systematic effects, while the rest are mainly due to statistical effects. We compare our results with other determinations of $Y_P$ present in the literature. Combining our $Y_P$ value with computations of primordial nucleosynthesis we find a number of neutrino species $N_{eff} = 2.90\pm0.22$, and a neutron mean life $\tau_{\nu} = 872\pm14(s)$.

The chemical content of the planetary nebula NGC 3918 is investigated through deep, high-resolution (R~40000) UVES at VLT spectrophotometric data. We identify and measure more than 750 emission lines, making ours one of the deepest spectra ever taken for a planetary nebula. Among these lines we detect very faint lines of several neutron-capture elements (Se, Kr, Rb, and Xe), which enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up in the progenitor star of NGC 3918.

(Abridged) The chemical content of the planetary nebula NGC 3918 is investigated through deep, high-resolution UVES at VLT spectrophotometric data. We identify and measure more than 750 emission lines, making ours one of the deepest spectra ever taken for a planetary nebula. Among these lines we detect very faint lines of several neutron-capture elements (Se, Kr, Rb, and Xe), which enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up in the progenitor star of NGC 3918. We find that Kr is strongly enriched in NGC 3918 and that Se is less enriched than Kr, in agreement with the results of previous papers and with predicted s-process nucleosynthesis. We also find that Xe is not as enriched by the s-process in NGC 3918 as is Kr and, therefore, that neutron exposure is typical of modestly sub-solar metallicity AGB stars. A clear correlation is found when representing [Kr/O] vs. log(C/O) for NGC 3918 and other objects with detection of multiple ions of Kr in optical data, confirming that carbon is brought to the surface of AGB stars along with s-processed material during third dredge-up episodes, as predicted by nucleosynthesis models. We also detect numerous refractory element lines (Ca, K, Cr, Mn, Fe, Co, Ni, and Cu). We compute physical conditions from a large number of diagnostics. Thanks to the high ionization of NGC 3918 we detect a large number of recombination lines of multiple ionization stages of C, N, O and Ne. The abundances obtained for these elements by using recently-determined state-of-the-art ICF schemes or simply adding ionic abundances are in very good agreement, demonstrating the quality of the recent ICF scheme for high ionization planetary nebulae.

We present deep, high-resolution (R~40000) UVES at VLT spectrophotometric data of the planetary nebula NGC 3918. This is one of the deepest spectra ever taken of a planetary nebula. We have identified and measured more than 700 emission lines and, in particular, we have detected very faint lines of several neutron-capture elements (s-process elements: Kr, Xe and Rb) that enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up.

Analysis of emission lines in gaseous nebulae yields direct measures of physical conditions and chemical abundances and is the cornerstone of nebular astrophysics. Although the physical problem is conceptually simple, its practical complexity can be overwhelming since the amount of data to be analyzed steadily increases; furthermore, results depend crucially on the input atomic data, whose determination also improves each year. To address these challenges we created PyNeb, an innovative code for analyzing emission lines. PyNeb computes physical conditions and ionic and elemental abundances, and produces both theoretical and observational diagnostic plots. It is designed to be portable, modular, and largely customizable in aspects such as the atomic data used, the format of the observational data to be analyzed, and the graphical output. It gives full access to the intermediate quantities of the calculation, making it possible to write scripts tailored to the specific type of analysis one wants to carry out. In the case of collisionally excited lines, PyNeb works by solving the equilibrium equations for an n-level atom; in the case of recombination lines, it works by interpolation in emissivity tables. The code offers a choice of extinction laws and ionization correction factors, which can be complemented by user-provided recipes. It is entirely written in the python programming language and uses standard python libraries. It is fully vectorized, making it apt for analyzing huge amounts of data. The code is stable and has been benchmarked against IRAF/NEBULAR. It is public, fully documented, and has already been satisfactorily used in a number of published papers.

The propagation of a shock wave into a medium is expected to heat the material beyond the shock, producing noticeable effects in intensity line ratios such as [O III]/Halpha. To investigate the occurrence of shocks in planetary nebulae (PNe), we have used all narrowband [O III] and Halpha images of PNe available in the HST archive to build their [O III]/Halpha ratio maps and to search for regions where this ratio is enhanced. Regions with enhanced [O III]/Halpha emission ratio can be ascribed to two different types of morphological structures: bow-shock structures produced by fast collimated outflows and thin skins enveloping expanding nebular shells. Both collimated outflows and expanding shells are therefore confirmed to generate shocks in PNe. We also find regions with depressed values of the [O III]/Halpha ratio which are found mostly around density bounded PNe, where the local contribution of [N II] emission into the F656N Halpha filter cannot be neglected.

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).

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.

We present IACTalks, a free and open access seminars archive (http://iactalks.iac.es) aimed at promoting astronomy and the exchange of ideas by providing high-quality scientific seminars to the astronomical community. The archive of seminars and talks given at the Instituto de Astrofi\'isica de Canarias goes back to 2008. Over 360 talks and seminars are now freely available by streaming over the internet. We describe the user interface, which includes two video streams, one showing the speaker, the other the presentation. A search function is available, and seminars are indexed by keywords and in some cases by series, such as special training courses or the 2011 Winter School of Astrophysics, on secular evolution of galaxies. The archive is made available as an open resource, to be used by scientists and the public.

We present results from integral field optical spectroscopy with the Potsdam Multi-Aperture Spectrograph of the Herbig-Haro (HH) object HH 204, with a spatial sampling of 1 x 1 arcsec^2. We have obtained maps of different emission lines, physical conditions and ionic abundances from collisionally excited lines. The ionization structure of the object indicates that the head of the bow shock is optically thick and has developed a trapped ionization front. The density at the head is at least five times larger than in the background ionized gas. We discover a narrow arc of high T_e([N II]) values delineating the southeast edge of the head. The temperature in this zone is about 1,000 K higher than in the rest of the field and should correspond to a shock-heated zone at the leading working surface of the gas flow. This is the first time this kind of feature is observed in a photoionized HH object. We find that the O^+ and O abundance maps show anomalous values at separate areas of the bow shock probably due to: a) overestimation of the collisional de-excitation effects of the [O II] lines in the compressed gas at the head of the bow shock, and b) the use of a too high T_e([N II]) at the area of the leading working surface of the flow.

This workshop brought together scientists (including atomic physicists, theoretical astrophysicists and astronomers) concerned with the completeness and accuracy of atomic data for astrophysical applications. The topics covered in the workshop included the evaluation of uncertainties in atomic data, the propagation of such uncertainties in chemical abundances, and the feedback between observations and calculations. On a different level, we also discussed communication issues such as how to ensure that atomic data are correctly understood and used, and which forum is the best one for a fluid interaction between all communities involved in the production and use of atomic data. This paper reports on the discussions held during the workshop and introduces AstroAtom, a blog created as a platform for timely and open discussions on the needs and concerns over atomic data, and their effects on astronomical research. The complete proceedings will be published on http://astroatom.wordpress.com/.

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.

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 a review on the determination of the primordial helium abundance Yp, based on the study of hydrogen and helium recombination lines in extragalactic HII regions. We also discuss the observational determinations of the increase of helium to the increase of oxygen by mass Delta Y / Delta O, and compare them with predictions based on models of galactic chemical evolution.

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 present results of deep echelle spectrophotometry of the brightest knot of the HH202 in the Orion Nebula --HH202-S-- using the ultraviolet Visual Echelle Spectrograph (UVES). The high spectral resolution has permitted to separate the component associated with the ambient gas from that associated with the gas flow. We derive electron densities and temperatures for both components, as well as the chemical abundances of several ions and elements from collisionally excited lines, including the first determinations of Ca^{+} and Cr^{+} abundances in the Orion Nebula. We also calculate the He^{+}, C^{2+}, O^{+} and O^{2+} abundances from recombination lines. The difference between the O^{2+} abundances determined from collisionally excited and recombination lines --the so-called abundance discrepancy factor-- is 0.35 dex and 0.11 dex for the shock and nebular components, respectively. Assuming that the abundance discrepancy is produced by spatial variations in the electron temperature, we derive values of the temperature fluctuation parameter, t^2, of 0.050 and 0.016, for the shock and nebular components, respectively. Interestingly, we obtain almost coincident t^2 values for both components from the analysis of the intensity ratios of He I lines. We find significant departures from case B predictions in the Balmer and Paschen flux ratios of lines of high principal quantum number n. We analyze the ionization structure of HH202-S, finding enough evidence to conclude that the flow of HH202-S has compressed the ambient gas inside the nebula trapping the ionization front. We measure a strong increase of the total abundances of nickel and iron in the shock component, the abundance pattern and the results of photoionization models for both components are consistent with the partial destruction of dust after the passage of the shock wave in HH202-S.

We present results from integral field spectroscopy with the Potsdam multi-Aperture Spectrograph of the head of the Herbig-Haro object HH 202 with a spatial sampling of 1"x1". We have obtained maps of different emission lines, physical conditions --such as electron temperature and density-- and ionic abundances from recombination and collisionally excited lines. We present the first map of the Balmer temperature and of the temperature fluctuation parameter, t^2. We have calculated the t^2 in the plane of the sky, which is substantially smaller than that determined along the line of sight. We have mapped the abundance discrepancy factor of O^{2+}, ADF(O^{2+}), finding its maximum value at the HH 202-S position. We have explored the relations between the ADF(O^{2+}) and the electron density, the Balmer and [O III] temperatures, the ionization degree as well as the t^2 parameter. We do not find clear correlations between these properties and the results seem to support that the ADF and t^2 are independent phenomena. We have found a weak negative correlation between the O^{2+} abundance determined from recombination lines and the temperature, which is the expected behaviour in an ionized nebula, hence it seems that there is not evidence for the presence of super-metal rich droplets in H II regions.

(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 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.

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.

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.

We have derived a primordial helium abundance of Yp = 0.2477 +- 0.0029, based on new atomic physics computations of the recombination coefficients of He I and of the collisional excitation of the H I Balmer lines together with observations and photoionization models of metal-poor extragalactic H II regions. The new atomic data increase our previous determination of Yp by 0.0086, a very significant amount. By combining our Yp result with the predictions made by the standard Big Bang nucleosynthesis model, we find a baryon-to-photon ratio, \eta, in excellent agreement both with the \eta value derived by the primordial deuterium abundance value observed in damped Lyman-\alpha systems and with the one obtained from the WMAP observations.

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 presently accepted "Theory of the Universe" was pioneered 60 years ago by Gamow, Alpher and Herman. As a consequence of the, later dubbed, Hot Big-Bang, matter was neutrons, and after some decay protons, and a history of successive captures built up the elements. It wasn't until some 15 years later (with the discovery of the Cosmic Microwave Background radiation) that Gamow and colleagues theories were validated and present day Standard Big-bang Nucleosynthesis theory was developed. We will discuss the importance of state of the art observations and modelling in the quest to determine precise values of the primordial abundance of D and 4He, using observations of astrophysical objects and modern day atomic parameters. In particular, we will present the search for understanding and coping with systematic errors in such determinations.

We present a review on the determination of the primordial helium abundance, Yp, based on the study of hydrogen and helium recombination lines in extragalactic H II regions. We also discuss the observational determinations of the increase of helium to the increase of oxygen by mass Delta Y/Delta O, and compare them with predictions based on models of galactic chemical evolution.

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.

(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.

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.

We discuss the effects of collisional enhancement of Balmer lines on the determination of the primordial helium abundance. To this aim, we present a photoionization model of the metal-poor extragalactic HII region SBS 0335--052. We show that the derived helium abundance (Y) of this HII region depends on the amount of collisional excitation affecting the Balmer lines, both directly (through an underestimation of the actual He/H ratio) and indirectly (through an overestimation of the interstellar reddening). We detail how each of these effects affects the derived value of Y.

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.

Recent advances on the quest to determine the primordial or pregalactic helium abundance, Yp, are reviewed. There are seven problems affecting the He/H abundance determinations of H II regions that are briefly discussed: the underlying absorption lines in the observed spectra, the ionization, temperature, and density structures of each object, the collisional contribution to the intensity of the He I and H I lines, the optical thickness of the He I lines, and the extrapolation to derive Yp based on the Y and O/H values.

The primordial abundance of 4He, Y_P, is one of the hottest themes in present-day astronomy, mostly due to its cosmological relevance. The disagreement between different determinations has been currently reduced to the 1-2 percent level, but these differences are still large enough to have deep implications for Big-Bang nucleosynthesis. It is therefore crucial to estimate precisely the uncertainties involved in the measurement of Y_P. Here, I review the methods used in the determination of Y_P and the related uncertainties. I also discuss some recent results, and emphasize the assumptions underlying the differences among them.

We present high resolution spectroscopic data of the two giant extragalactic H II regions NGC 5461 and NGC 5471 in M101, which have been obtained with the 2.1-m telescope of the Observatorio Astron\'omico Nacional in San Pedro M\'artir, Baja California. We measured the intensities of several H and He recombination lines, and of forbidden lines of a large number of ions. We calculate the physical conditions in the two nebulae with a large number of diagnostics and determine their chemical abundances by applying ionization correction factors icf's to the observed ionic abundances. For NGC 5461, the icf's are based on a tailored photoionization model of the region (Luridiana & Peimbert 2001), while for NGC 5471 they are computed from those predicted by a photoionization model of NGC 2363 (Luridiana, Peimbert, & Leitherer 1999), a region which is similar to NGC 5471 in the ionization structure. For both regions, the icf's are compared to those computed following the prescriptions by Mathis & Rosa (1991). Such comparison shows large discrepancies for several elements, including nitrogen, neon, and chlorine.

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 study the effect that the temperature structure has on the determination of the primordial helium abundance, Yp. We provide an equation linking T(OIII), the temperature derived from the [OIII] lines, and T(HeII), the temperature of the He I lines, both for HII regions with O++ only and for HII regions where a fraction of O+ is present. By means of T(HeII), which is always smaller than T(OIII), we derive the helium abundances of 5 objects with low and very low metallicity (NGC 346, NGC 2363, Haro29, SBS 0335-052, and I Zw 18); these objects were selected from the literature because they include the 3 low metallicity objects with the best line determinations and the 2 objects with the lowest metallicity. From these abundances we obtain that Yp(nHc)=0.2356+-0.0020, a value 0.0088 lower than that derived by using T(OIII). We call this determination Yp(nHc) because the collisional contribution to the Balmer line intensities has not been taken into account. All the recent Yp determinations in the literature have not taken into account the collisional contribution to the Balmer line intensities. By considering the collisional contribution to the Balmer line intensities of these five objects we derive that Yp(+Hc)=0.2384+-0.0025.

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.

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 compute photoionization models for the giant extragalactic H II region NGC 5461, and compare their predictions to several observational constraints. Since we aim at reproducing not only the global properties of the region, but its local structure also, the models are constrained to reproduce the observed density profile, and our analysis takes into consideration the bias introduced by the shapes and sizes of the slits used by different observers. We find that an asymmetric nebula with a gaussian density distribution, powered by a young burst of 3.1 Myr, satisfactorily reproduces most of the constraints, and that the star-formation efficiency inferred from the model agrees with current estimates. Our results strongly depend on the assumed density law, since constant density models overestimate the hardness of the ionizing field, affecting the deduced properties of the central stellar cluster. We illustrate the features of our best model, and discuss the possible sources of errors and uncertainties affecting the outcome of this type of studies.

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.

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.

We quantify the energy radiated through all the collisionally excited lines in a photoionized nebula which is permeated by temperature fluctuations. We assume that these correspond to hot spots which are the results of an unknown heating process distinct from the photoelectric heating. We consider all the effects of using a higher mean temperature (as compared to the equilibrium temperature) due to the fluctuations not only on each emission line but also on the ionization state of the gas. If this yet unknown process was to radiate a fixed amount of energy, we find that the fluctuations should correlate with metallicity Z when it exceeds 0.7 solar. The excess energy radiated in the lines as a result of the fluctuations is found to scale proportionally to their amplitude t^2. When referred to the total energy absorbed through photoionization, the excess energy is comparable in magnitude to t^2.