Recent models show that TP-AGB stars should dominate the NIR spectra of populations 0.3 to 2 Gyr old, leaving unique signatures that can be used to detect young/intermediate stellar population in galaxies. However, no homogeneous database of star-forming galaxies is available in the NIR to fully explore these results. With this in mind, we study the NIR spectra of a sample of 23 HII and starburst galaxies, aimed at characterizing the most prominent spectral features and continuum shape in the 0.8-2.4 micron region of these objects. Spectral indices are derived for the relevant absorption lines/bands and a comparison with optical indices of the same sample available in the literature is made. We found no correlation between the optical and the NIR indexes. This is probably due to the differences in aperture between these two sets of data. That result is further supported by the absence or weakness of emission lines in the NIR for a subsample galaxies, while in the optical the emission lines are strong and clear, which means that the ionisation source in many of these galaxies is not nuclear, but circumnuclear or located in hot spots. We detected important signatures predicted for a stellar population dominated by the TP-AGBs, like CN 1.1 micron and CO 2.3 micron. In at least one galaxy (NGC 4102) the CN band at 1.4 micron was detected for the first time. We also detect TiO and ZrO bands that have never been reported before in extragalactic sources. The shape of the continuum emission is found to be strongly correlated to the presence/lack of emission lines. An observational template for the star-forming galaxies is derived to be used as a benchmark of stellar population(s) in starburst galaxies against which to compare near-IR spectroscopy of different types of galaxies, especially those with AGN activity and/or those at high-redshift.

Detections of molecular lines, mainly from H2$ and CO, reveal molecular material in planetary nebulae. Observations of a variety of molecules suggest that the molecular composition in these objects differs from that found in interstellar clouds or in circumstellar envelopes. The success of the models, which are mostly devoted to explain molecular densities in specific planetary nebulae, is still partial, however. The present study aims at identifying the influence of stellar and nebular properties on the molecular composition of planetary nebulae by means of chemical models. A comparison of theoretical results with those derived from the observations may provide clues to the conditions that favor the presence of a particular molecule. A self-consistent photoionization numerical code was adapted to simulate cold molecular regions beyond the ionized zone. The code was used to obtain a grid of models and the resulting column densities are compared with those inferred from observations. Our models show that the inclusion of an incident flux of X-rays is required to explain the molecular composition derived for planetary nebulae. We also obtain a more accurate relation for the N(CO)/N(H2) ratio in these objects. Molecular masses obtained by previous works in the literature were then recalculated, showing that these masses can be underestimated by up to three orders of magnitude. We conclude that the problem of the missing mass in planetary nebulae can be solved by a more accurate calculation of the molecular mass.

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H2 emission lines show that at least part of the molecular emission is produced inside the ionised region. In the best-studied case, the Helix nebula, the H2 emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionised gas. Most of the H2 emission of the CKs seems to be produced in a thin layer between the ionised diffuse gas and the neutral material of the knot, in a mini photodissociation region (PDR). However, PDR models published so far cannot fully explain all the characteristics of the H2 emission of the CKs. In this work, we use the photoionisation code \textsc{Aangaba} to study the H2 emission of the CKs, particularly that produced in the interface H^+/H^0 of the knot, where a significant fraction of the H2 1-0S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines. We find that the temperature derived from H2 observations even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H_alpha and NII peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the NII line is produced closer to the border of the CK than H_alpha.

The analysis and interpretation of the H2 line emission from planetary nebulae have been done in the literature assuming that the molecule survives only in regions where the hydrogen is neutral, as in photodissociation, neutral clumps or shocked regions. However, there is strong observational and theoretical evidence that at least part of the H2 emission is produced inside the ionized region of such objects. The aim of the present work is to calculate and analyze the infrared line emission of H2 produced inside the ionized region of planetary nebulae using a one-dimensional photoionization code. The photoionization code Aangaba was improved in order to calculate the statistical population of the H2 energy levels and the intensity of the H2 infrared emission lines in physical conditions typical of planetary nebulae. A grid of models was obtained and the results are analyzed and compared with the observational data. We show that the contribution of the ionized region to the H2 line emission can be important, particularly in the case of nebulae with high temperature central stars. This result explains why H2 emission is more frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object typically has hotter stars. Collisional excitation plays an important role on the population of the rovibrational levels of the electronic ground state of H2. Radiative mechanisms are also important, particularly for the upper vibrational levels. Formation pumping can have minor effects on the line intensities produced by de-excitation from very high rotational levels, especially in dense and dusty environments. We included the effect of the H2 on the thermal equilibrium of the gas, concluding that H2 only contributes to the thermal equilibrium in the case of a very high temperature of the central star or a high dust-to-gas ratio, mainly through collisional de-excitation.

In the present work, we use a photoionisation code to study the H2 emission of the Helix nebula (NGC 7293) cometary knots, particularly that produced in the interface H+/H0 of the knot, where a significant fraction of the H2 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines.

We performed stellar population synthesis on the nuclear and extended regions of NGC 1068 by means of near-infrared spectroscopy to disentangle their spectral energy distribution components. This is the first time that such a technique is applied to the whole 0.8 - 2.4 micron wavelength interval in this galaxy. NGC 1068 is one of the nearest and probably the most studied Seyfert 2 galaxy, becoming an excellent laboratory to study the interaction between black holes, the jets that they can produce and the medium in which they propagate. Our main result is that traces of young stellar population are found at ~ 100 south of the nucleus. The contribution of a power-law continuum in the centre is about 25%, which is expected if the light is scattered from a Seyfert 1 nucleus. We find peaks in the contribution of the featureless continuum about 100 - 150 pc from the nucleus on both sides. They might be associated with regions where the jet encounters dense clouds. Further support to this scenario is given by the peaks of hot dust distribution found around these same regions and the H2 emission line profile, leading us to propose that the peaks might be associate to regions where stars are being formed. Hot dust also has an important contribution to the nuclear region, reinforcing the idea of the presence of a dense, circumnuclear torus in this galaxy. Cold dust appears mostly in the south direction, which supports the view that the southwest emission is behind the plane of the galaxy and is extinguished very likely by dust in the plane. Intermediate age stellar population contributes significantly to the continuum, specially in the inner 200 pc.

We report the first simultaneous zJHK spectroscopy on the archetypical Seyfert 2 Galaxy NGC 1068 covering the wavelength region 0.9 to 2.4 micron. The slit, aligned in the NS direction and centred in the optical nucleus, maps a region 300 pc in radius at sub-arcsec resolution, with a spectral resolving power of 360 km s^-1. This configuration allow us to study the physical properties of the nuclear gas including that of the north side of the ionization cone, map the strong excess of continuum emission in the K-band and attributed to dust and study the variations, both in flux and profile, in the emission lines. Our results show that (1) Mid- to low-ionization emission lines are splitted into two components, whose relative strengths vary with the position along the slit and seem to be correlated with the jet. (2) The coronal lines are single-peaked and are detected only in the central few hundred of parsecs from the nucleus. (3) The absorption lines indicate the presence of intermediate age stellar population, which might be a significant contributor to the continuum in the NIR spectra. (4) Through some simple photoionization models we find photoionization as the main mechanism powering the emitting gas. (5) Calculations using stellar features point to a mass concentration inside the 100 - 200 pc of about 10^10 solar masses.

Mid-resolution spectra are used to deduce the size and kinematics of the coronal gas in a sample of Seyfert galaxies by means of observations of the [FeXI], [FeX], [FeVII], [SiVI] and [SiVII] lines. These coronal lines (CL) extend from the unresolved nucleus up to a few tens to a few hundreds of parsecs. The region of the highest ionized ions studied, [FeXI] and [FeX], is the least spatially extended, and concentrates at the center; intermediate ionization lines extend from the nucleus up to a few tens to a few hundred parsecs; lower [OIII]-like ions are known to extendin the kpc range. All together indicates a stratification in the ionized gas, usually interpreted in terms of nuclear photoionization as the driving ionization mechanism. However, CL profiles show various peculiarities: they are broader by a factor of two than lower ionization lines, the broadening being in terms of asymmetric blue wings, and their centroid position at the nucleus is blueshifted by a few hundreds of km/s. Moreover, in NGC1386 and NGC1068, a double peak [FeVII] line is detected in the nuclear and extended coronal region, this being the first report in of such type of profile in CL in active galactic nuclei. If interpreted as outflow signatures, the total broadening of the lines at zero intensity levels implies gas velocities up to 2000 km/s. Although the stratification of ions across the coronal region means that photoionization is the main power mechanism, the high velocities deduced from the profiles, the relatively large spatial extension of the emission, and the results from photoionization models indicate that an additional mechanism is at work. We suggest that shocks generated by the outflow could provide the additional required power for line formation.

When analysing HII regions, a possible source of systematic error on empirically derived physical quantities is the limited size of the slit used for the observations. A grid of photoionization models was built through the Aangaba code varying the ionizing radiation spectrum emitted by a stellar cluster, as well as the gas abundance. The calculated line surface brightness was then used to simulate slit observations and to derive empirical parameters using the usual methods described in the literature. Depending on the fraction of the object covered by the slit, the parameters can be different from those obtained from observations of the whole object, an effect that is mainly dependent on the age of the ionizing stellar cluster. The low-ionization forbidden lines are more sensitive to the size of the area covered by the slit than the high-ionization forbidden lines or recombination lines. Regarding the temperature indicator T[OIII], the slit effects are small since this temperature is derived from [OIII] lines. On the other hand, for the abundance indicator R23, which depends also on the [OII] line, the slit effect is slightly higher. Therefore, the systematic error due to slit observations on the O abundance is low, being usually less than 10%, except for HII regions powered by stellar clusters with a relative low number of ionizing photons between 13.6 and 54.4 eV, which create a smaller O++ emitting volume. In this case, the systematic error on the empirical O abundance deduced from slit observations is more than 10% when the covered area is less than 50%.

We present the results of a spatio-kinematic study of the planetary nebula Menzel 1 using spectro-photometric mapping and a 3-D photoionization code. We create several 2-D emission line images from our long-slit spectra, and use these to derive the line fluxes for 15 lines, the Halpha/Hbeta extinction map, and the [SII] line ratio density map of the nebula. We use our photoionization code constrained by these data to derive the three-dimensional nebular structure and ionizing star parameters of Menzel 1 by simultaneously fitting the integrated line intensities, the density map, and the observed morphologies in several lines, as well as the velocity structure. Using theoretical evolutionary tracks of intermediate and low mass stars, we derive a mass for the central star of 0.63+-0.05 Msolar. We also derive a distance of 1050+_150 pc to Menzel 1.

We present the results of mapping the planetary nebula NGC6369 using multiple long slit spectra taken with the CTIO 1.5m telescope. We create two dimensional emission line images from our spectra, and use these to derive fluxes for 17 lines, the Halpha/Hbeta extinction map, the [SII] line ratio density map, and the [NII] temperature map of the nebula. We use our photoionization code constrained by these data to determine the distance, the ionizing star characteristics, and show that a clumpy hour-glass shape is the most likely three-dimensional structure for NGC6369. Note that our knowledge of the nebular structure eliminates all uncertainties associated with classical distance determinations, and our method can be applied to any spatially resolved emission line nebula. We use the central star, nebular emission line, and optical+IR luminosities to show that NGC6369 is matter bound, as about 70% of the Lyman continuum flux escapes. Using evolutionary tracks from Bloecker(1995) we derive a central star mass of about 0.65 Msolar.

This paper presents an analysis of the concentration of the hydrogen molecule inside the ionized region of planetary nebulae. The equations corresponding to the ionization and chemical equilibria of H, H+, H-, H2, H2+, and H3+ are coupled with the equations of ionization and thermal balance for a photoionized atomic gas. Forty different reactions related to the formation or the destruction of these species are included. The presence of dust is taken into account, since grains act as catalysts for the production of H2, as well as shield the molecules against the stellar ionizing radiation. We analyze the effect of the stellar ionizing continuum, as well as of the gas and grain properties on the calculated H2 mass. It is shown that a significant concentration of H2 can survive inside the ionized region of planetary nebulae, mostly in the inner region of the recombination zone. The total H2 to total hydrogen mass ratio inside the ionized region increases with the central star temperature, and, depending on the PN physical conditions, it can be of the order of 10^-6 or even higher. The increase of the recombination zone with the stellar temperature can account for such correlation. This can explain why the H2 emission is more frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object has typically hotter stars. Applying our results for the planetary nebula NGC 6720, we obtain an H2 to hydrogen mass ratio similar to the value obtained from the observed H2 line emission.

The narrow line region of active galaxies is formed by gas clouds surrounded by a diluted gas. Standard one-dimensional photoionization models are usually used to model this region in order to reproduce the observed emission lines. Since the narrow line region is not homogeneous, two major types of models are used: (a) those assuming a homogeneous gas distribution and a filling factor less than unity to mimic the presence of the emitting clouds; (b) those based on a composition of single-cloud models combined in order to obtain the observed spectra. The first method is largely used but may induce to misleading conclusions as shown in this paper. The second one is more appropriate, but requires a large number of observed lines in order to limit the number of single models used. After discussing the case of an extragalactic HII region, for which the ionizing radiation spectrum is better known, we show that 1-D models for the narrow line region with a filling factor less than unit do not properly mimic the clumpiness, but just simulates an overall lower density. Multi-cloud models lead to more reliable results. Both models are tested in this paper, using the emission-line spectra of two well-known Seyfert galaxies, NGC 4151 and NGC 1068. It is shown that ionizing radiation spectra with a blue bump cannot be excluded by multi-cloud models, although excluded by Alexander et al. (1999, 2000)using homogeneous models with a filling factor less than unity.

Helium and hydrogen recombination lines observed in low-metallicity, extragalactic, HII regions provide the data used to infer the primordial helium mass fraction, Y_P. In deriving abundances from observations, the correction for unseen neutral helium or hydrogen is usually assumed to be absent; i.e., the ionization correction factor is taken to be unity (icf = 1). In a previous paper (VGS), we revisited the question of the icf, confirming a "reverse" ionization correction: icf < 1. In VGS the icf was calculated using more nearly realistic models of inhomogeneous HII regions, suggesting that the published values of Y_P needed to be reduced by an amount of order 0.003. As star clusters age, their stellar spectra evolve and so, too, will their icfs. Here the evolution of the icf is studied, along with that of two, alternate, measures of the "hardness" of the radiation spectrum. The differences between the icf for radiation-bounded and matter-bounded models are also explored, along with the effect on the icf of the He/H ratio (since He and H compete for some of the same ionizing photons). Particular attention is paid to the amount of doubly-ionized helium predicted, leading us to suggest that observations of, or bounds to, He++ may help to discriminate among models of HII regions ionized by starbursts of different ages and spectra. We apply our analysis to the Izotov & Thuan (IT) data set utilizing the radiation softness parameter, the [OIII]/[OI] ratio, and the presence or absence of He++ to find 0.95 < icf < 0.99. This suggests that the IT estimate of the primordial helium abundance should be reduced by Delta-Y = 0.006 +- 0.002, from 0.244 +- 0.002 to 0.238 +- 0.003.

The effect of the nebula optical depth on the determination of the temperature (T*) of the central stars in planetary nebulae is discussed. Based on photoionization models for planetary nebulae with different optical depths, we show, quantitatively, that the details of the distribution of the H and He II Zanstra temperatures are mainly explained by an optical depth effect; in particular, that the discrepancy is larger for low stellar temperatures. The results also show that for high stellar temperatures the He II Zanstra temperature underestimates the stellar temperature, even for high optical depths. The stellar temperature, as well as the optical depth, can be obtained from a Zanstra temperature ratio (ZR) plot ZR = Tz(He II)/Tz(H) versus Tz(He II). The effects of departures from a blackbody spectrum, as well as of the He abundance in the nebulae, are also discussed. For nebulae of very low optical depth and/or high stellar temperature the distribution ZR versus Tz(He II) only provides lower limits for T*. In order to obtain better values for the optical depth and T*, we propose the use of the line intensity ratio He II/He I versus Tz(He II) diagram.

We use a 3D photoionization modeling tool to study the morpho-kinematic properties of the Planetary Nebula (PN) NGC 3132. We show that it is possible to reproduce the low resolution observations (spectra and images) with an ellipsoidal shell. However, high resolution observations, as those showing a density variation along the nebula and the [O III]5007 velocity profiles, definitively rule out this description. We show that a bipolar Diabolo shape with a 40 degres rotation of the symmetry axis relative to the line of sight successfully reproduces the observed images, as well as the high resolution observations.

We present a new modeling tool for planetary nebulae, based on 3D photoionization calculations. Our goal is to show that all the information provided by observations, regarding kinematics and morphology, have to be consistently accounted for, in order to get a real insight of the object. Only 3D simulations offer this possibility. From models for two theoretical PNe, we show that the enhancement in the equatorial zone observed in several PNe is not necessarily due to a density gradient, as usually interpreted. It is also shown that asymmetric velocity profiles often observed (e.g., Gesicki et al. 1998) can be easily reproduced. Observations providing a better insight on the morphology of the PN are discussed.

Helium and hydrogen recombination lines observed in low-metallicity, extragalactic H II regions provide the data used to infer the primordial helium mass fraction, Y_P. The ionization corrections for unseen neutral helium (or hydrogen) are usually assumed to be absent; i.e., the ionization correction factor is taken to be unity (icf = 1). In this paper we revisit the question of the icf for H II regions ionized by clusters of young, hot, metal-poor stars. Our key result is that for the H II regions used in the determination of Y_P, there is a ``reverse'' ionization correction: icf < 1. We explore the effect on the icf of more realistic inhomogeneous H II region models and find that for those regions ionized by young stars, with ``hard'' radiation spectra, the icf is reduced further below unity. In Monte Carlos using H II region data from the literature (Izotov and Thuan 1998) we estimate a reduction in the published value of Y_P of order 0.003, which is roughly twice as large as the quoted statistical error in the Y_P determination.

We present a new modeling tool for planetary nebulae (PNe), based on 3D photoionization calculations. From models for two theoretical PNe, we show that the enhancement in the equatorial zone observed in several PNe is not necessarily due to a density gradient, as usually interpreted. We study the morpho-kinematic properties of the PN NGC 3132 and show that a bipolar Diabolo shape successfully reproduces the observed images, as well as the high resolution observations.

A hot gas halo is predicted by chemodynamical models during the early evolution of spheroidal galaxies. Cold condensations, arising from thermal instabilities in the hot gas, are expected to be embedded in the hot halo. In the early phases of the galaxy (t < 1 Gyr), a strong X-ray and EUV emission is produced by the extended hot gas distribution, ionizing the cold clouds. This self-irradiating two-phase halo model successfully explains several line ratios observed in QSO absorption-line systems, and reproduces the temperature distribution of Lyman-$\alpha$ clouds.

Type I planetary nebulae (PNe) are defined as those with high He and N abundances (Peimbert & Torres-Peimbert 1983). These objects present in general bipolar geometries and have high stellar temperatures (Corradi & Schwarz 1995, Torres-Peimbert & Peimbert 1997). In this paper we analyse the empirical methods for abundance determination in order to check if the He and N overabundances in Type I PNe are a consequence of a geometrical effect, due to the bipolarity, or the ionization stratification, due to the stellar temperature. For this, we obtain simulated spherically symmetrical as well as bipolar nebulae, using a 3D photoionization code. From the projected emission line intensities for: a) the whole nebula; b) for a slit crossing the nebula; as well as c) for different positions in the nebula, we applied the formulae used in the literature to obtain empirical abundances. These empirical abundances are then compared with the adopted ones. We show that empirical abundances depend on the particular line of sight covered by the observation and can simulate an overabundance and/or the presence of abundance gradients of He and N in planetary nebulae with high stellar temperature. The geometrical effects are also discussed. Systematic errors in abundance determinations by empirical methods are higher for the N/H ratio than for N/O. Thus, it seems better to use the N/O value when discussing N rich objects.

There is evidence for temperature fluctuations in Planetary Nebulae and in Galactic HII regions. If such fluctuations occur in the low-metallicity, extragalactic HII regions used to probe the primordial helium abundance, the derived 4He mass fraction, Y_P, could be systematically different from the true primordial value. For cooler, mainly high-metallicity HII regions the derived helium abundance may be nearly unchanged but the oxygen abundance could have been seriously underestimated. For hotter, mainly low-metallicity HII regions the oxygen abundance is likely accurate but the helium abundance could be underestimated. The net effect is to tilt the Y vs. Z relation, making it flatter and resulting in a higher inferred Y_P. Although this effect could be large, there are no data which allow us to estimate the size of the temperature fluctuations for the extragalactic HII regions. Therefore, we have explored this effect via Monte Carlos in which the abundances derived from a fiducial data set are modified by \Delta-T chosen from a distribution with 0 < \Delta-T < \Delta-T_max where \Delta-T_max is varied from 500K to 4000K. It is interesting that although this effect shifts the locations of the HII regions in Y vs. O/H plane, it does not introduce any significant additional dispersion.